• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.64-69
• /
• 1995

A computer code for solving the Reynolds averaged full Navier-Stokes equations has bent developed for analysis of gas and steam turbine cascade flows with the option of using one of two types of turbulence model. One is the Baldwin-Lomax model and the other is standard $k-{\varepsilon}$ model. The numerical integration is based on the explicit four stage Runge-Kutta scheme and finite volume method. To be verified, the resulting code is applied to VKI turbine cascade and compared with the previous experimental results. Finally, the flowfield around a steam turbine cascade is analyzed. Comparisons with experimental data show that present numerical scheme is an accurate Navier-Stokes solver and can give very good predictions for both gas and steam turbine cascade flow.

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

The screech tone of an underexpanded jet is numerically calculated without any specific modeling for the screech tone itself. A fourth-order optimized compact scheme and fourth-order Runge-Kutta method are used to solve the 2D axisymmetric Euler equation. The Fourier transform of pressure signal at upstream shows the directivity pattern of the screech tone very clearly. Pressure signal is shown to observe the generation of the screech tone. Most importantly, we can simulate the axisymmetric mode change of the screech tone very precisely with the proposed method. It can be concluded that the basic phenomenon of the screech tone including its frequency can be calculated and its mode change can be simulated with inviscid Euler equations.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.783-786
• /
• 2002

The unsteady turbulent flow characteristics of NACA4406 airfoil at close proximity to the pound are numerically investigated, especially focused on the noise generation mechanism near the blunt trailing edge. The unsteady two-dimensional compressible Wavier-Stokes equations with a Spalart-Allmaras turbulence closure model are solved by the 6th-order compact scheme and the 4th-order Runge-Kutta scheme. The computation shows a noise generation by a feedback mechanism at the blunt tailing edge, where the acoustic-fluidic coupling occurs between the wall-reflected sound waves and the periodically disturbed turbulent shear layer.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.115-126
• /
• 1999

A two-dimensional Navier-Stokes solver using hybrid meshes is parallelized with a domain decompostion method. The focus of this paper is placed on minimizing the amount of effort in parallelizing the serial version of the solver, and this is achieved by adding an additional layer of cells to each decomposed domain. Most subroutines of the serial solver are used without modification, and the information exchange between neighboring domains is achieved using MPI(Message Passing Interface) library. Load balancing among the processors and scheduling of the message passing are implemented to reduce the overhead of parallelization, and the speed-up achieved by parallelization is measured on the transonic invisicd and turbulent flow problems. The parallelization efficiencies of the explicit Runge-Kutta scheme and the implicit point-SGS scheme are compared and the effects of various factors on the results are also studied.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• Journal of Mechanical Science and Technology
• /
• v.17 no.8
• /
• pp.1219-1225
• /
• 2003

Acoustic sounds generated by uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method. A third-order-accurate up-wind scheme is used for the spatial derivatives. A second-order-accurate Runge-Kutta scheme is also used for time marching. Very small acoustic pressure fluctuation, with same frequency as that of Karman vortex street, is compared with pressure fluctuation around a circular cylinder. The propagation velocity of acoustic sound shows that acoustic approaching the upstream, due to the Doppler effect in uniform flow, slowly propagates. For the downstream, on the other hand, it quickly propagates. It is also apparent that the size of sound pressure is proportional to the central distance ${\gamma}$$\^$-1/2/ of the circular cylinder.

• Journal of Ship and Ocean Technology
• /
• v.5 no.1
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• Proceedings of the SAREK Conference
• /
• 2006.06a
• /
• pp.116-121
• /
• 2006

Numerical simulations for dendritic growth of crystals are conducted in this study by the level set method. The effect of order of difference is tested for reinitialization error in simple problems and authors founded in case of 1st order of difference that very fine grids have to be used to minimize the error and higher order of difference is desirable to minimize the reinitialization error The 2nd and 4th order Runge-Kutta scheme in time and 3rd and 5th order of WENO schemes with Godunov scheme are applied for space discretization. Numerical results are compared with the analytical theory, phase-field method and other researcher's level set method.

• The Bulletin of The Korean Astronomical Society
• /
• v.45 no.1
• /
• pp.55.1-55.1
• /
• 2020

In an attempt to investigate the nonlinear dynamics such as shock, shear, and turbulence associated with ultra-relativistic jets, we develop a new relativistic hydrodynamics (RHD) code based on the weighted essentially non-oscillatory (WENO) scheme. It is a 5th-order accurate, finite-difference scheme, which has been widely used for solving hyperbolic systems of conservation equations. The code is parallelized with MPI and OpenMP. Through an extensive set of tests, the accuracy and efficiency of different WENO reconstructions, and different time discretizations are assessed. Different implementations of the equation of state (EOS) for relativistic fluid are incorporated, As the fiducial setup for simulations of ultra-relativistic jets, we adopt the EOS in Ryu et al. (2006) to treat arbitrary adiabatic index of relativistic fluid, the WENO-Z reconstructions to minimize numerical dissipation without loss of stability, and the strong stability preserving Runge-Kutta (SSPRK) method to achieve stable time stepping with large CFL numbers. In addition, the code includes a high-order flux averaging along the transverse directions for multi-dimensional problems, and the modified eigenvalues for the acoustic modes to effectively control the carbuncle instability. We find that the new code performs satisfactorily simulations of ultra-relativistic jets.