• Journal of computational fluids engineering
• /
• v.4 no.3
• /
• pp.21-27
• /
• 1999

In this paper, tile characteristics of flow resulting from the configurations of piston head and intake-port of the cylinder in a gasoline-direct-injection engine are investigated numerically. Calculations are carried out from intake process to the end of compression. GTT code which includes the third order upwind Chakravarthy-Osher TVD scheme and κ-ε turbulence model with the law of wall as a boundary condition. As a result, a piston head with a smaller radius of curvature and larger radius gives stronger reverse tumble. It is also shown that as the maximum tumble ratio increases by the configuration of the intake-port the tumble ratio at the end of compression stroke increases. It is concluded that flows at the end of compression stroke can be controlled by the optimum design of intake-port and piston head.

• 유체기계공업학회:학술대회논문집
• /
• 1998.12a
• /
• pp.223-232
• /
• 1998

To investigate the flow inside the centrifugal impeller, computer program which can solve Three-dimensional compressible turbulent flow has been developed. The Navier-Stokes equations were chosen as the governing equation for viscous flow while Euler equations for inviscid case. Time marching method was incorporated with the Flux Difference Splitting method suggested by Roe to capture the steep gradients such as a shock. For high order of accuracy, MUSCL approach was adopted while differentiable limiter to ensure TVD property. For turbulence closure, Baldwin- Lomax model was applied due to its simplicity. To demonstrate the capabilities of present program, several validation problems have been solved and compared with experiments and other available data. From the above calculations generally good agreements were obtained. Finally, the developed code was applied to Eckardt's impeller and the performance prediction was carried out. Some important aspects on boundary condition for successful simulation were discussed and the remedy was also introduced.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.597-601
• /
• 2008

Three-dimensional structures of detonation wave propagating in circular tube were investigated. Inviscid fluid dynamics equations coupled with a conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Variable-$\gamma$ formulation was used to account for the variable properties between unburned and burned states and the chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The computational code was parallelized based on domain decomposition technique using MPI-II message passing library. The computations were carried out using a home made Windows based PC cluster having 160 AMD AthloxXP and Athlon64 processor. The computational domain consisted of through a roundshaped tube with wall conditions. As an initial condition, analytical ZND solution was distributed over the computational domain with disturbances. The disturbances has circumferential large gradient. The unsteady computational results in three-dimension show the detailed mechanisms of multi-cell mode of detonation wave instabilities resulting diamond shape in smoked-foil record.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• The Bulletin of The Korean Astronomical Society
• /
• v.40 no.1
• /
• pp.65.2-65.2
• /
• 2015

We developed a three-dimensional magnetohydrodynamic (MHD) code to reproduce the structure of a solar wind, the properties of a coronal mass ejection (CME) and the interaction between them. This MHD code is based on the finite volume method incorporating total variation diminishing (TVD) scheme with an unstructured grid system. In particular, this grid system can avoid the singularity at the north and south poles and relax tight CFL conditions around the poles, both of which would arise in a spherical coordinate system (Tanaka 1994). In this model, we first apply an MHD tomographic method (Hayashi et al. 2003) to interplanetary scintillation (IPS) observational data and derive a solar wind from the physical values obtained at 50 solar radii away from the Sun. By comparing the properties of this solar wind to observational data obtained near the Earth orbit, we confirmed that our model captures the velocity, temperature and density profiles of a solar wind near the Earth orbit. We then insert a spheromak-type CME (Kataoka et al. 2009) into the solar wind to reproduce an actual CME event. This has been done by introducing a time-dependent boundary condition to the inner boundary of our simulation domain. On the basis of a comparison between a simulated CME and observations near the Earth, we discuss the physics involved in an ICME interacting with a solar wind.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.6
• /
• pp.461-468
• /
• 2012

In this study, a Cartesian-grid method based on finite volume approach is applied to simulate the ship motions in large amplitude waves. Fractional step method is applied for pressure-velocity coupling and TVD limiter is used to interpolate the cell face value for the discretization of convective term. Water, air, and solid phases are identified by using the concept of volume-fraction function for each phase. In order to capture the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with weighed line interface calculation (WLIC) method which considers multidimensional information. The volume fraction of solid body embedded in the Cartesian grid system is calculated using a level-set based algorithm, and the body boundary condition is imposed by a volume weighted formula. Numerical simulations for the two-dimensional barge type model and Wigley hull in linear waves have been carried out to validate the newly developed code. To demonstrate the applicability for highly nonlinear wave-body interactions such as green water on the deck, numerical analysis on the large-amplitude motion of S175 containership is conducted and all computational results are compared with experimental data.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.2
• /
• pp.1-7
• /
• 2002

In the modern propulsion systems, requited thrust is obtained using a nozzle. Sometimes shock and induced boundary layer separation is generated in an over-expanded convergent-divergent supersonic nozzle. It occurs because the nozzle expansion ratio is too large for a given nozzle pressure ratio (NPR). This phenomenon can be explained that it redefines effective nozzle geometry, shorer nozzle geometry and lower pressure ratio, in a given pressure ratio. Numerical studies were conducted about a fixed geometry 2D nozzle in overexpanded condition and compared with Hunter's experimental result. For the numerical simulation of the supersonic nozzle, Navier-Stokes equations are considered and as a turbulent model, $\kappa$-$\varepsilon$ /$\kappa$-$\omega$ blended SST two equation turbulent model is used. The characteristics of $\lambda$-shape shock systems due to the interaction of shock and boundary layer was investigated in a low NPR. And the result of comparison of thrust value shows that a fixed geometry nozzle can cover required flight mission.