• Journal of The Korean Astronomical Society
• /
• v.28 no.2
• /
• pp.223-243
• /
• 1995

We describe the implementation of a multi-dimensional numerical code to solve the equations for idea! magnetohydrodynamics (MHD) in cylindrical geometry. It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. Curvature and source terms are included in a way to insure the formal accuracy of the code to be second order. The constraint of a divergence-free magnetic field is enforced exactly by adding a correction, which involves solving a Poisson equation. The Fourier Analysis and Cyclic Reduction (FACR) method is employed to solve it. Results from a set of tests show that the code handles flows in cylindrical geometry successfully and resolves strong shocks within two to four computational cells. The advantages and limitations of the code are discussed.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.60-64
• /
• 2009

Laminar separation bubble and transitional flow over an airfoil are investigated at a moderate range of Reynolds numbers. In this research, a Reynolds-Averaged Navier-Stokes code is coupled with an empirical transition model that can predict transition onset points and the length of transition region. Without solving the boundary layer equations, approximated e-N method is directly applied to the RANS code and iteratively solved together. The computational results are compared with the experimental data for NACA0012 airfoil. Results of transition onset point and length are compared well with experimental and XFOIL prediction. In high angle of attack the present RANS results show better agreement than XFOIL results using the boundary layer equations.

• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.62-69
• /
• 2004

The present work solved 3D incompressible RANS equations on a rotating, multi-blocked grid system to efficiently analyze ducted marine propulsor with the interaction of propeller guidevane and annular duct. To handle the interface boundary between the guidevane and the propeller, a sliding multiblock technique based on the cubic spline interpolation was applied. To validate the present code, a turbine flow was simulated and the time-averaged pressure coefficients were compared with experiment. After the code validation, the flowfield around a ducted marine propeller with pre-swirl guidevane was simulated.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.93-98
• /
• 2001

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.1-8
• /
• 2010

Laminar separation bubble and transitional flow over the NACA0012 are investigated at a moderate range of Reynolds numbers. A Reynolds-Averaged Navier-Stokes code is coupled with an empirical transition model that can predict transition onset points and the length of transition region. Without solving the boundary layer equations, approximated e-N method is directly applied to the RANS code and iteratively solved together. The computational results are compared with the experimental data for the NACA0012 airfoil. Results of transition onset point and the length are compared well with experimental data and Xfoil prediction. The present RANS results show at high angles of attack better agreement with experimental data than Xfoil results using the boundary layer equations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.4
• /
• pp.458-465
• /
• 2003

The main objective of the research is to develop a research code solving transient incompressible Navier-Stokes equation. In this research code, Adams-Bashforth method was applied to the convective terms of the navier stokes equation and the splitted equations were discretized spatially by finite element methods to solve the complex geometry problems easily. To reduce the divergence on the boundaries of pressure poisson equation due to the unsuitable pressure boundary conditions, multi step approximation pressure boundary conditions derived from the boundary linear momentum equations were used. Simulations of Lid Driven Flow and Flow over Cylinder were conducted to prove the accuracy by means of the comparison with results of the previous workers.

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.47-54
• /
• 2000

A robust Navier-Stokes code has been developed to efficiently predict hypersonic flows in chemical nonequilibrium. The HLLE+ flux discretization scheme is used to improve accuracy and robustness of hypersonic flow analysis. An efficient LU approximate factorization method is also used to solve the flow equations and species continuity equations in fully coupled fashion to implicitly treat stiff source terms of chemical reactions. The HLLE+ scheme shows lower grid dependency for the wall heating rates than other schemes. The developed code has been used to compute chemical nonequilibrium air flow through expanding hypersonic nozzle and past two and three dimensional blunt-nosed bodies. The results are in good agreement with existing numerical and experimental results.

• Nuclear Engineering and Technology
• /
• v.29 no.4
• /
• pp.310-319
• /
• 1997

A finite element analysis of a trunnion pipe anchor is presented. The structure is analyzed for the case of internal pressure and moment loadings. The stress results are categorized into the average (membrane) stress, the linearly varying (bending) stress and the peak stress through the thickness. The resulting stresses are interpreted per Section III of the ASME Boiler and Pressure Vessel Code from which the Primary(B$_1$), Secondary(C$_1$) and Peak(K$_1$) stress indices for pressure, the Primary (B$_2$), Secondary(C$_2$) and Peak(K$_2$) stress indices for moment are developed. Based on the comparison between stress value by stress indices derived in this paper and stress value represented by the ASME Code Case N-391-1, the empirical equations for stress indices are effectively used in the piping stress analysis. Therefore, the use of empirical equations can simplify the procedure of evaluating the local stress in the piping design stage.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2401-2405
• /
• 2008

In the present study, a three-dimensional least square/level set based two-phase flow code was developed for the simulation of three-dimensional sloshing problems using finite element discretization. The present method can be utilized for the analysis of a free surface flow problem in a complex geometry due to the feature of FEM. Since the finite element method is employed for the spatial discretization of governing equations, an unstructured mesh can be naturally adopted for the level set simulation of a free surface flow without an additional load for the code development except that solution methods of the hyperbolic type redistancing and advection equations of the level set function should be devised in order to give a bounded solution on the unstructured mesh. From the numerical experiments of the present study, it is shown that the proposed method is both robust and accurate for the simulation of three-dimensional sloshing problems.