• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.6
• /
• pp.1-8
• /
• 2004

The objective of this study is convergence acceleration of multigrid Navier-Stokes solvers. This study has been performed to enhance the performance of preconditioned multi-stage time stepping method which is a popular smoother for the multigrid Navier-Stokes solvers. Comparative study on the convergence characteristics of the ADI and DDADI preconditioners has been conducted. It is shown that the DDADI preconditioner has better performance than the ADI by numerical tests on the 2-D compressible turbulent flows past airfoils. The Spalart-Allmaras turbulent model and the Baldwin-Lomax turbulent model have been compared with the multigrid calculations.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.2
• /
• pp.113-120
• /
• 2009

To numerically simulate aerodynamics of rotor-airframe interaction in a rigorous manner, we need to solve the Navier-Stokes system for a rotor-airframe combination in a single computational domain. This imposes a computational burden since rotating blades and a stationary body have to be simultaneously dealt with. An efficient alternative is a momentum source method in which the action of rotor is approximated as momentum source in a stationary mesh system built around the airframe. This makes the simulation much easier. The magnitude of the momentum source is usually evaluated by the blade element theory, which often results in a poor accuracy. In the present work, we evaluate the momentum source from the simulation data by using the Navier-Stokes equations only for a rotor system. Using this data, we simulated the time-averaged steady rotor-airfame interaction and developed the unsteady rotor-airframe interaction. Computations were carried out for the simplified rotor-airframe model (the Georgia Tech configuration) and the results were compared with experimental data. The results were in good agreement with experimental data, suggesting that the present approach is a usefull method for rotor-airframe interaction analysis.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2004.11a
• /
• pp.221-231
• /
• 2004

A kinetic theory analysis is made of low-speed gas flows in a microfluidic system consisted of three microchannels in series. The Boitzmann equation simplified by a collision model is solved by means of a finite difference approximation with the discrete ordinate method. For the evaluation of the present method results are compared with those from the DSMC method and an analytical solution of the Navier-Stokes equations with slip boundary conditions. Calculations are made for flows at various Knudsen numbers and pressure ratios across the channel. The results compared well with those from the DSMC method. It is shown that the analytical solution of the Navier-Stokes equations with slip boundary conditions which is suited fur fully developed flows can give relatively good results. In predicting the geometrically complex flows up to a Knudsen number of about 0.06. It is also shown that the present method can be used to analyze extremely low-speed flow fields for which the DSMC method is Impractical.

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

In this article, combination of the FAS-FMG multi-grid method and the Krylov subspace method was presented in solving two dimensional driven-cavity flows. Three algorithms of the Krylov subspace method, CG, CGSTAB(Bi-CG Stabilized) and GMRES method were tested with MILU preconditioner. As a smoother of the pressure correction equation, the MILU-CG is recommended rather than MILU-GMRES(k) or MILU-CGSTAB, since the MILU-GMRES(k) preconditioner has too much computation on the coarse grid compared to the MILU-CG one. As for the momentum equation, relatively cheap smoother like SIP solver may be sufficient.

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

Three-dimensional flow analysis and numerical optimization methods are presented for the design of an axial-flow fan. Steady, Incompressible, three-dimensional Reynolds-averaged Wavier-Stokes equations are used as governing equations, and standard k-$\epsilon$ turbulence model is chosen as a turbulence model. Governing equations are discretized using finite volume method. Steepest descent method, conjugate gradient method and BFGS method are compared to determine the searching directions. Golden section method and quadratic fit-sectioning method are tested for one dimensional search. Objective function is defined as a ratio of generation rate of the turbulent kinetic energy to pressure head. Sweep angle distributions are used as design variables.

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

• The KSFM Journal of Fluid Machinery
• /
• v.7 no.3 s.24
• /
• pp.32-38
• /
• 2004

Shape of a multi-blades centrifugal fan is optimized by response surface method based on three-dimensional Navier-Stokes analysis. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard $k-{epsilon}$ 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. Optimizations with and without constraints are carried out. 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. The correlation of efficiency with relative size of inactive zone at the exit of impeller is discussed as well as with average momentum fluxes in the scroll.