• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.48-54
• /
• 1998

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• 한국가시화정보학회:학술대회논문집
• /
• 2004.11a
• /
• pp.72-75
• /
• 2004

The present study describes a computational work to investigate detailed behaviors of the twin shock waves discharged from the exits of two-parallel ducts. In computations, the Yee-Roe-Davis's TVD scheme was used to solve the unsteady, three-dimensional, inviscid, compressible, Euler equations. The distance between two ducts is varied and the Mach number of the incident shock wave is changed below 2.0. The results obtained show that on the symmetric axis between two-parallel ducts, the maximum pressure achieved by the merge of twin shock waves and its location strongly depend upon the distance between two-parallel ducts and the Mach number of the incident shock wave. It is also found that the twin shock waves discharged from the exits of two-parallel ducts leads to the complicated flow fields, such as Mach stem, spherical waves, and vertical structures.

• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.

• Journal of Mechanical Science and Technology
• /
• v.14 no.8
• /
• pp.867-878
• /
• 2000

In order to study unsteady aerodynamic loads on high speed trains passing by each other 350km/h, three-dimensional flow fields around trains during the crossing event are numerically simulated using three-dimensional Euler equations. Roe's FDS with MUSCL interpolation is employed to simulate wave phenomena. An efficient moving grid system based on domain decomposition techniques is developed to analyze the unsteady flow field induced by the restricted motion of a train on a rail. Numerical simulations of the trains passing by on the double-track are carried out to study the effect of the train nose-shape, length and the existence of a tunnel on the crossing event. Unsteady aerodynamic loads-a side force and a drag force-acting on the train during the crossing are numerically predicted and analyzed. The side force mainly depends on the nose-shape, and the drag force depends on tunnel existence. Also. a push-pull (i.e.impluse force) force successively acts on each car and acts in different directions between the neighborhood cars. The maximum change of the impulsive force reaches about 3 tons. These aerodynamic force data are absolutely necessary to evaluate the stability of high speed multi-car trains. The results also indicate the effectiveness of the present numerical method for simulating the unsteady flow fields induced by bodies in relative motion.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.249-257
• /
• 2007

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• Journal of the Korean Society for Railway
• /
• v.5 no.4
• /
• pp.260-266
• /
• 2002

A numerical simulation has been performed to estimate the transient pressure variation in the tunnel when G7 test train passes through the test tunnel in the Kyoeng-Bu high-speed railway. A modified patched grid scheme is developed to handle the relative motion between a train and a tunnel. Also, a hybrid dimensional approach is proposed to calculate the train-tunnel interaction problem efficiently. An axi-symmetric unsteady Euler solve using the Roe's FDS is used for analyzing a complicated pressure field in tunnel during the test train is passing through the tunnel. Usually, this complex phenomenon depends ell the train speed, train length, tunnel length, blockage ratio between train and tunnel cross-sectional area, relative position between train and tunnel, etc. Therefore, numerical simulation should be done carefully in consideration of these factors. Numerical results in this study would be good guidance to make test plans, test equipments selection and to decide their measuring locations. They will also supply important information to the pressurization equipment for high-speed train.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.599-604
• /
• 2000

In order to study unsteady aerodynamic loads on high speed trains passing by each other at the speed of 350km/h, three-dimensional flow fields around trains during the crossing event are numerically simulated using the three-dimensional Euler equations. The Roe's FDS with MUSCL interpolation is employed to simulate wave phenomena properly. An efficient moving grid system based on domain decomposition techniques is developed to analyze the unsteady flow field induced by the restricted motion of a train on a rail. The numerical simulations of the trains passing by on the double-track are carried out to study the effect of the train nose-shape, the train length and the existence of tunnel when the crossing event occur. Unsteady aerodynamic loads side force and drag force-acting on the train during the crossing are numerically predicted and anlayzed. It is found that the strength of the side force mainly depends on the nose-shape, and that of drag force on tunnel existence. And it is observed that the push-pull like impulsive force successively acts on each car and acts in different directions between the neighborhood cars. The maximum change of the impulsive force reaches about 3 tons. These aerodynamic force data are absolutely necessary for the evaluation of the stability of the high speed multi-car train. The results also indicate the effectiveness of the present numerical method for the simulation of unsteady flow field induced by the bodies in the relative motion.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1595-1600
• /
• 2004

When a shock wave arrives at a duct, an impulsive wave is discharged from the duct exit and causes serious noise and vibration problems. In the current study, the characteristics of the impulsive wave discharged from a partial closed duct exit is numerically investigated using a CFD method. The Yee-Roe- Davis's total variation diminishing(TVD) scheme is used to solve the axisymmetric, unsteady, compressible Euler equations. With several partial closed duct exits, the Mach number of the incident shock wave $M_s$ and the distance L/D between the duct exit and a flat plate are varied in the range of $M_s$ = 1.01 ${\sim}$ 1.50 and L/D = 1.0 ${\sim}$ 4.0, respectively. The results obtained show that the magnitude of the impulsive wave impinging upon the flat plate strongly depends upon $M_s$, L/D and the partial closure of duct exit. The impulsive wave on the flat plate can be considerably alleviated by the partial closure of duct exit and, thus, the present method can be a passive control for the impulsive wave.