• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.7 s.250
• /
• pp.671-681
• /
• 2006

This study focuses on the development of numerical procedure to analyze the nonlinear combustion instabilities in liquid rocket engine. Nonlinear behaviors of acoustic instabilities are characterized by the existence of limit cycle in linearly unstable engines and nonlinear or triggering instability in linearly stable engines. To discretize convective fluxes with high accuracy and robustness, approximated Riemann solver based on characteristics and Euler-characteristic boundary conditions are employed. The present procedure predicts well the transition processes from initial harmonic pressure disturbance to N-like steep-fronted shock wave in a resonant pipe. Longitudinal pressure oscillations within the SSME(Space Shuttle Main Engine) engine have been analyzed using the pressure-sensitive time lag model to account for unsteady combustion response. It is observed that the pressure oscillations reach a limit cycle which is independent of the characteristics of the initial disturbances and depends only on combustion parameters and operating conditions.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.330-340
• /
• 2005

Compressible two-phase flow is analyzed based on the arbitrary Lagrangian-Eulerian (ALE) formulation. For water, Tamman type stiffened equation of state is used. Numerical fluxes are calculated using the ALE two-phase Godunov scheme which assumes only that the speed of sound and pressure can be provided whenever density and internal energy are given. Effects of the approximations of a material interface speed are Investigated h method Is suggested to assign a rigid body boundary condition effectively To validate the developed code, several well-known problems are calculated and the results are compared with analytic or other numerical solutions including a single material Sod shock tube problem and a gas/water shock tube problem The code is applied to analyze the refraction and transmission of shock waves which are impacting on a water-gas interface from gas or water medium.

• International Journal of Aeronautical and Space Sciences
• /
• v.12 no.4
• /
• pp.318-330
• /
• 2011

This paper investigates the convergence characteristics of the modified artificial compressibility method proposed by Turkel. In particular, a focus is mode on the convergence characteristics due to variation of the preconditioning factor (${\alpha}_u$) and the artificial compressibility (${\beta}$) in conjunction with an upwind method. For the investigations, a code using the modified artificial compressibility is developed. The code solves the axisymmetric incompressible Reynolds averaged Navier-Stokes equations. The cell-centered finite volume method is used in conjunction with Roe's approximate Riemann solver for the inviscid flux, and the central difference discretization is used for the viscous flux. Time marching is accomplished by the approximated factorization-alternate direction implicit method. In addition, Menter's k-${\omega}$ shear stress transport turbulence model is adopted for analysis of turbulent flows. Inviscid, laminar, and turbulent flows are solved to investigate the accuracy of solutions and convergence behavior in the modified artificial compressibility method. The possible reason for loss of robustness of the modified artificial compressibility method with ${\alpha}_u$ >1.0 is given.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.3
• /
• pp.307-316
• /
• 2012

Numerical simulations of 3D aircraft configurations are performed in order to understand the effects of turbulence models on the prediction of aircraft's aerodynamic characteristics. An in-house CFD code that solves 3D RANS equations and two-equation turbulence model equations are used. The code applies Roe's approximated Riemann solver and an AF-ADI scheme. Van Leer's MUSCL extrapolation with van Albada's limiter is also adopted. Various versions of Menter's $k-{\omega}$ SST turbulence models as well as Coakley's $q-{\omega}$ model are incorporated into the CFD code. Menter's $k-{\omega}$ SST models include the standard model, the 2003 model, the model incorporating the vorticity source term, and the model containing controlled decay. Turbulent flows over a wing are simulated in order to validate the turbulence models contained in the CFD code. The results from these simulations are then compared with computational results from the $3^{rd}$ AIAA CFD Drag Prediction Workshop. Numerical simulations of the DLR-F6 wing-body and wing-body-nacelle-pylon configurations are conducted and compared with computational results of the $2^{nd}$ AIAA CFD Drag Prediction Workshop. Aerodynamic characteristics as well as flow features are scrutinized with respect to the turbulence models. The results obtained from each simulation incorporating Menter's $k-{\omega}$ SST turbulence model variations are compared with one another.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.