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

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.110-115
• /
• 1998

An axisymmetric Navier-Stokes procedure has been developed to analyze the pulse jet flow in a ceramic filter unit for the dust dislodging process. Using Baldwin-Lomax turbulence model as a closure relationship, the SIAF(Scalar Implicit Approximate Factorization) algorithm together with the ${\delta}^k-Correction$ iterative time marching scheme is adopted to solve the unsteady compressible Navier-Stokes equations. After some validation tests, the code has been applied to solve the pulse jet flow and examine the effects of geometry and reservoir pressure condition on the pressure level inside the filter unit. To avoid dealing with the uncertainty of such factors as the cohesion of the collected dust and the adhesion of the dust to the medium and also to simplify the analysis, the filter wall is assumed to be impermeable. The results for various test cases are presented.

• Journal of Energy Engineering
• /
• v.10 no.2
• /
• pp.166-175
• /
• 2001

핀-관 열교환기의 효율을 증대시키기 위하여는 열저항을 결정하는데 있어서 중요한 역할을 하는 공기측 열전달 특성의 향상이 필요하다. 본 연구에서는 슬릿핀과 같은 핀 형상의 개선으로 인한 높은 압력손실을 감소시키기 위하여 타원관형 열교환기를 고려하였으며 이를 자체 개발한 3차원 비압축성 Navier-Stokes 코드를 이용하여 공기측 성능을 해석하였다. 이 코드는 시간항에 스칼라 내재적 근사분해법(scalar implicit approximate factorization)절차, 공간항에 유한체적법과 2차의 풍상차분법을 사용한다. 핀 형상이 평판핀인 타원관에 대하여 세장비가 0.3, 0.5, 0.75인 경우와 핀 형상이 슬릿핀이고 세장비가 0.5인 타원관에 대하여 공기측 유동을 해석하였으며 이를 원관의 성능과 비교하였다. 평판핀의 경우, 압력강하량은 원관과 비교하여 세장비가 작아질수록 감소하는 것을 볼 수 있었으며, 열전달 계수는 고려한 유입속도 범위에서 원관을 기준으로 약 5%이하의 변화로 거의 동일하게 유지되는 것을 볼 수 있었다. 슬릿핀의 경우 역시 문제되는 높은 압력강하량은 타원관을 도입함으로써 효과적으로 낮출 수 있음을 볼 수 있었다.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.179-186
• /
• 1999

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67), and to evaluate the performances of k-$\epsilon$ and Baldwin-Lomax turbulence models. A finite volume method is used for spatial discretization. And, the equations are solved implicitly in time with the use of approximate factorization. Upwind difference scheme is used for inviscid terms, but viscous terms are centrally differenced. The flux-difference-splitting of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. And, the results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, we conclude that k-$\epsilon$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost same.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.459-471
• /
• 1994

A numerical scheme solving the quasi-3D wave-induced circulation is presented. The governing equations have been solve implicitly using a fractional step method in conjunction with the approximate factroization techniques. The equation of each step was discretized by the finite volume scheme which yields more accurate and conservative approximations than the schemes based on finite differences. Examples of computed nearshore current patterns are presented to demonstrate the validity of the model for typical situations through comparison with laboratory experimental data (Gourlay. 1974).

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

The low-speed preconditioning technique is applied to solve the compressible Reynolds averaged Navier-Stokes equations for low-speed flows. The space discretization is based on Roe's flux-difference splitting with third-order-accurate MUSCL extrapolation. Time integration is performed employing a diagonal approximate factorization algorithm. The dual-time stepping has been incorporated to solve the unsteady flows. Full multigrid method is implemented to accelerate the convergence rate. To verify the algorithms several cases have been tested. Demonstrated the improvement on convergence and quality of the solution.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. The results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, it is concluded that Abid'k-$\omega$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost the same.

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.4
• /
• pp.181-189
• /
• 2000

A weakly nonlinear mild-slope equation has been derived directly from the continuity equation with the aid of the Galerkin's method. The equation is combined with the momentum equations defined at the mean water level. A single component model has also been obtained in terms of the surface displacement. The linearized form is completely identical with the time-dependent mild-slope equation proposed by Smith and Sprinks(1975). For the verification purposes of the present nonlinear model, the degenerate forms were compared with Airy(1845)'s non-dispersive nonlinear wave equation, classical Boussinesq equation, andsecond¬order permanent Stokes waves. In this study, the present nonlinear wave equations are discretized by the approximate factorization techniques so that a tridiagonal matrix solver is used for each direction. Through the comparison with physical experiments, nonlinear wave model capacity was examined and the overall agreement was obtained.