• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.41-52
• /
• 1999

이 논문에서는, QUICK해법의 불안정성을 개량하므로써, 수치계산에 있어서 수렴이 빠르고, 수치적으로 안정한 계산을 할 수 있는 새로운 MQUICK 상류해법을 제안하고, 이를 비압축성 층류유동의 계산에 적용하였다. 또한, 해법의 정확성, 안정성, 수렴속도에 대한 검토를 통하여 본 MQUICK 상류해법의 유효성과 타당성이 평가되었다. 이 해법에서는 인공산일의 가감을 조절하기 위하여 가중계수 α를 써서 정식화 하였고, 위의 검토를 통하여 α의 최적값을 조사하였다. 이 해법을 SMAC 음해법에 적용하여 2 차원 공동유동, 3 차원 덕트유동과 같은 몇몇 표준문제를 계산하고, 계산된 결과를 실험값 또는, 3 차 정확도의 상류해법 및 QUICK해법에 의한 결과 들과 비교 하므로써, 본 MQUICK 상류해법이 위의 다른 해법에 비하여 안정하고, 유효성이 높은 해법임을 확인 하였다.

• Journal of Power System Engineering
• /
• v.3 no.4
• /
• pp.16-21
• /
• 1999

Numerical simulation on two-dimensional turbine cascade flow has been performed using compressible Navier-Stokes equations. The flow equations are written in a cartesian coordinate system, then mapped into a generalized body-fitted ones. All direction of viscous terms are incoporated and turbulent effects are modeled using the extended ${\kappa}-{\epsilon}$ model. Equations are discretized using control volume SIMPLE algorithm on the nonstaggered grid sysetem. Applications are made at a VKI turbine cascade flow in atransonic wind-tunnel and compared to experimental data. Present numerical results are shown to be in good agreement with the experimental results and simulate the compressible viscous flow characteristics inside the turbine blade passage.

• Korean Journal of Mathematics
• /
• v.21 no.2
• /
• pp.125-150
• /
• 2013

The stabilized Gauge-Uzawa method (SGUM), which is a second order projection type algorithm to solve the time-dependent Navier-Stokes equations, has been newly constructed in 2013 Pyo's paper. The accuracy of SGUM has been proved only for time discrete scheme in the same paper, but it is crucial to study for fully discrete scheme, because the numerical errors depend on discretizations for both space and time, and because discrete spaces between velocity and pressure can not be chosen arbitrary. In this paper, we find out properties of the fully discrete SGUM and estimate its errors and stability to solve the evolution Navier-Stokes equations. The main difficulty in this estimation arises from losing some cancellation laws due to failing divergence free condition of the discrete velocity function. This result will be extended to Boussinesq equations in the continuous research (part II) and is essential in the study of part II.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.177-185
• /
• 1999

The airfoil design optimization procedures based on the Navier-Stokes equations were developed, This procedure enables more realistic and practical transonic airfoil designs. The modified Hicks-Henne functions were used to generate the shape of airfoils. Five Hick-Henne functions were used to design upper surface of airfoil only. To enhance the ability of Hick-Henne function to generate various airfoil shape with limited number of functions, the positions of control points were adjusted through optimization procedure. The design procedure was applied to the single-point design for the drag minimization problem with lift and area constraints. The result shows the capability of the procedure to generate much realistic airfoils with very small drag-creep in the low transonic regime. This is mainly due to the viscosity effect of Navier-Stokes flow analysis. However, in the higher transonic range tile drag-creep appears. The multi-point design is shown to be an effective way to avoid the drag-creep and improve off-design performance which is very similar in the Euler design.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.2 s.152
• /
• pp.130-138
• /
• 2007

A coupled variational equation for fluid-structure interaction (FSI) problems is derived from a steady state Navier-Stokes equation for incompressible Newtonian fluid and an equilibrium equation for geometrically nonlinear structures. For a fully coupled FSI formulation, between fluid and structures, a traction continuity condition is considered at interfaces where a no-slip condition is imposed. Under total Lagrange formulation in the structural domain, finite rotations are well described by using the second Piola-Kirchhoff stress and Green-Lagrange strain tensors. An adjoint shape design sensitivity analysis (DSA) method based on material derivative approach is applied to the FSI problem to develop a shape design optimization method. Demonstrating some numerical examples, the accuracy and efficiency of the developed DSA method is verified in comparison with finite difference sensitivity. Also, for the FSI problems, a shape design optimization is performed to obtain a maximal stiffness structure satisfying an allowable volume constraint.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.10
• /
• pp.1457-1463
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a multi-blade centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k - c turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.54-61
• /
• 2015

In the present study, a flow solver using a kinetic BGK scheme was developed for the compressible Navier-Stokes equation. The kinetic BGK scheme was used to simulate flow field from the continuum up to the transitional regime, because the kinetic BGK scheme can take into account the statistical properties of the gas particles in a non-equilibrium state. Various numerical simulations were conducted by the present flow solver. The laminar flow around flat plate and the hypersonic flow around hollow cylinder of flare shape in the continuum regime were numerically simulated. The numerical results showed that the flow solver using the kinetic BGK scheme can obtain accurate and robust numerical solutions. Also, the present flow solver was applied to the hypersonic flow problems around circular cylinder in the transitional regime and the results were validated against available numerical results of other researchers. It was found that the kinetic BGK scheme can similarly predict a tendency of the flow variables in the transitional regime.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2157-2161
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k-e turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.474-484
• /
• 2018

In this study, a new numerical modeling system was proposed to predict oil spills, which increasingly occur at sea as a result of abnormal weather conditions such as global warming. The hydrodynamic conditions such as the flow velocity needed to calculate oil dispersion were estimated using a three dimensional hydrodynamic model based on the Navier-Stokes equation, which considered all of the physical variations in the vertical direction. This improved the accuracy compared to those estimated by the conventional shallow water equation. The advection-diffusion model for the spilled oil was combined with the hydrodynamic model to predict the movement and fate of the oil. The effects of absorption, weathering, and wind were also considered in the calculation process. The combined model developed in this study was then applied to various test cases to identify the characteristics of oil dispersion over time. It is expected that the developed model will help to establish initial response and disaster prevention plans in the event of a nearshore oil spill.