• Journal of Ocean Engineering and Technology
• /
• v.15 no.4
• /
• pp.20-27
• /
• 2001

Recently numerical approaches for wave loads acting on the vertical caisson of breakwater, and resulting wave reflection and transmission coefficients have been performed. Although the numerical studies by Sulisz's(1997) and Kim et al.(2000) are suggested representatively, theoretical formulation for nonlinear wave pressure is not developed yet. And experimental results of Sulisz(1997) revealed that nonlinear uplift pressure on the caisson may be produced largely on the case of caisson founded on the high rubble mound. From the results of this study, the nonlinear theory for the uplift wave pressure acting on the caisson by applying boundary integral method of Green theorem is formulated, and also the characteristics of nonlinear uplift pressure and run-up height on the caisson are evaluated numerically, according to the variations of hydraulic properties of the rubble mound.

• International Journal of Safety
• /
• v.1 no.1
• /
• pp.18-23
• /
• 2002

A mixture fraction formulation is used to numerically simulate the structure of diluted axisymmetric methane-air nonpremixed counterflow flames. The effects of global strain rate and gravity wert! investigated and results were compared. Fuel of a mixture of 20% methane and 80% nitrogen by volume and oxidizer of pure air at low and moderate global strain rates $a_g= 20, 40, 80 s^{-1}$ in normal and zero gravity were computed. It is shown that the numerical method is capable of predicting the structure of counterflow flames in normal and microgravity environments at low and moderate global strain rates.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2000.11a
• /
• pp.63-71
• /
• 2000

In this study a drift prediction model based on fluid dynamics theory is introduced. The essential effects of environmental loads and target characteristics are taken into account from a fluid dynamics point of view. The governing equations of motion are derived from Newton's law of dynamics. In the mathematical formulation only three degrees of freedom(surge, sway, yaw) of the drifting object are assumed and the environmental loads considered are the forces and moments by wind and current. A computer algorithm for this model is implemented to obtain the numerical result in the time domain. The preliminary tests for model verification are conducted and the results are compared with the field experiment data as well as leeway formula suggested from the field test data.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.8-14
• /
• 2001

The phenomena of two-phase suspension flows appear widely in nature and industrial processes. Hence, it is of great importance to understand the mechanism of the gas-solid two-phase flows. In the present study, the numerical simulation has been approached by utilizing the Eulerian-Lagrangian methodology for describing the characteristics of the fluid and particulate phases in a vertical pipe and a 90°square-sectioned bend. The continuous phase(gas phase) is described by the Eulerian formulation and a κ-ε turbulence model is employed to find mean and turbulent properties of the gas phase. The particle properties(velocity and trajectory) are then described by a Lagrangian approach and computed using the mean velocity and turbulent fluctuating velocity of the gas phase. The predictions are compared with measurements by laser-Doppler velocimeter for the validation. As a result, the calculated results show good agreements.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.8-17
• /
• 2004

A Semi-Lagrangian method based on CIP(Cubic Interpolated Pseudoparticle)method is proposed and it is applied to solve the two dimensional advection equation. Especially the attentions are given to settle the pole problem and to enhance the accuracy in solving the advection equation on the spherical coordinate system. Tn this algorithm, the CU method is employed as the Semi-Lagrangian method and extended to the spherical coordinate system. To enhance the accuracy of the solution, the spatial discretization is made by CIP method. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional rotating cosine bell problem and the frontogenesis problem are simulated by the present scheme. As results, it is confirmed that the present scheme gives an accurate solution and settles the pole problem in the advection equation on the sphere.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.578-583
• /
• 2007

Many papers have studied computer-aided simulations of elastic bodies undergoing large deflections and large deformations. But there have been few attempts to validate their numerical formulations used in these studies. The main aim of this paper is to validate the absolute nodal coordinate formulation (ANCF) by comparing the results to experimental measurements on beams. Physical experiments with a high-speed camera were carried out to capture the large displacement of the beam and to verify the results of computer simulations. To consider the damping forces, the Rayleigh's damping and quadratic damping are employed and compared to the experimental results, respectively. Numerical results obtained from computer simulations were compared with the results from the physical experiments according to the $1^{st}$ mode and the $2^{nd}$ mode of the beam, respectively.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.06a
• /
• pp.775-778
• /
• 2007

The present study numerically investigate the detailed structure of the syngas diffusion flames. In order to realistically represent the turbulence-chemistry interaction, the transient flamelet model has been applied to simulate the combustion processes and $NO_X$ formation in the syngas turbulent nonpremixed flames. The single mixture fraction formulation is extended to account for the effects of the secondary inlet mixture. Computations are the wide range of syngas compositions and oxidizer dilutions. Based on numerical results, the detailed discussion has been made for the effects of syngas composition and oxidizer dilution on the structure of the syngas-air and syngas-oxygen turbulent nonpremixed flames.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• pp.200-205
• /
• 1998

Although discretization methods such as the transfer matrix method (TMM) and the finite element method (FEM) have played an important role in the design or analysis of rotor-bearing systems, continuous system modeling and analysis are often desirable especially for sensitivity analysis or design. The present paper proposes a comprehensive modeling procedure to obtain exact solution of general rotor-bearing systems. The proposed method considers a Timoshenko beam model and makes use of complex coordinate in the formulation. The proposed method provides exact eigensolutions and frequency response functions (FRFS) of general multi-stepped rotor-bearing systems. The first numerical example compares the proposed method with FEM. The numerical study proves that the proposed method is very efficient and useful for the analysis of rotor-bearing systems.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.1023-1026
• /
• 2007

Optimal damping layout of the constrained viscoelastic damping layer on beam is identified with temperatures by using a gradient-based numerical search algorithm. An optimal design problem is defined in order to determine the constrained damping layer configuration. A finite element formulation is introduced to model the constrained damping layer beam. The four-parameter fractional derivative model and the Arrhenius shift factor are used to describe dynamic characteristics of viscoelastic material with respect to frequency and temperature. Frequency-dependent complex-valued eigenvalue problems are solved by using a simple resubstitution algorithm in order to obtain the loss factor of each mode and responses of the structure. The results of the numerical example show that the proposed method can reduce frequency responses of beam at peaks only by reconfiguring the layout of constrained damping layer within a limited weight constraint.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2001.05a
• /
• pp.245-249
• /
• 2001

This paper is concerned with numerical analyses for bulk forming of a single crystal milli-product, whose typical size ranges from a few hundreds ${\mu}m$ to a few mm. The numerical formulation invoked in this paper combines the crystal plasticity theory considering texture development and the ductile damage mechanics for growth of micro voids, since orientation development and growth of micro voids become the primary factors for bulk forming of milli-size products. As applications, milli-extrusion of a single crystal round bar and milli-rolling of a single crystal plate are simulated and the results are discussed in detail.