• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1581-1593
• /
• 2001

An efficient variable-reordering method for finite element meshes is used and the effect of variable-reordering is investigated. For the element renumbering of unstructured meshes, Cuthill-McKee ordering is adopted. The newsy reordered global matrix has a much narrower bandwidth than the original one, making the ILU preconditioner perform bolter. The effect of variable reordering on the convergence behaviour of saddle point type matrix it studied, which results from P2/P1 element discretization of the Navier-Stokes equations. We also propose and test 'level(0) preconditioner'and 'level(2) ILU preconditioner', which are another versions of the existing 'level(1) ILU preconditioner', for the global matrix generated by P2/P1 finite element method of incompressible Navier-Stokes equations. We show that 'level(2) ILU preconditioner'performs much better than the others only with a little extra computations.

• Journal of the Korean Mathematical Society
• /
• v.48 no.3
• /
• pp.529-550
• /
• 2011

In this work, a numerical solution of the incompressible Navier-Stokes equations is proposed. The method suggested is based on an algorithm of discretization by mixed finite elements with a posteriori error estimation of the computed solutions. In order to evaluate the performance of the method, the numerical results are compared with some previously published works or with others coming from commercial code like Adina system.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.3
• /
• pp.1015-1027
• /
• 1996

In this paper, a numerical procedure for the calculation of the incompressible Navier-Stokes equations in a generalized nonorthogonal body fitted coordinate system is proposed and is validated through three test problems. Present numerical procedure derives the pressure equation by using the pressure substitution method on the regular grid system, and discretized momentum equations are based on the covariant velocity components. Cavity flow, backward facing step flow, and two dimensional channel flow with a sinusoidal wavy wall are chosen as three test problems. Numerical solutions obtained by present procedure shows a good agreement with previous numerical and/or experimental results. Convergence rate is also satisfactory.

• Korean Journal of Mathematics
• /
• v.27 no.1
• /
• pp.193-219
• /
• 2019

The Navier-Stokes equations with variable density are challenging problems in numerical analysis community. We recently built the 2nd order stabilized Gauge-Uzawa method [SGUM] to solve the Navier-Stokes equations with constant density and have estimated theoretically optimal accuracy. Also we proved that SGUM is unconditionally stable. In this paper, we apply SGUM to the Navier-Stokes equations with nonconstant variable density and find out the stability condition of the algorithms. Because the condition is rather strong to apply to real problems, we consider Allen-Cahn scheme to construct unconditionally stable scheme.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.8
• /
• pp.1950-1963
• /
• 1995

A modification of the approximate-factorization method is made to accelerate the convergency rate and to take sufficiently large Courant number without loss of accuracy. And a stable implicit finite-difference scheme for solving the incompressible Navier-Stokes equations employed above modified method is developed. In the present implicit scheme, the volume fluxes with contravariant velocity components and the pressure formulation in curvilinear coordinates is adopted. In order to satisfy the continuity condition completely and to remove spurious errors for the pressure, the Navier-Stokes equations are solved by a modified SMAC scheme using a staggered gird. The upstream-difference scheme such as the QUICK scheme is also employed to the right hand side. The implicit scheme is unconditionally stable and satisfies a diagonally dominant condition for scalar diagonal linear systems of implicit operator on the left hand side. Numerical results for some test calculations of the two-dimensional flow in a square cavity and over a backward-facing step are obtained using both usual approximate-factorization method and the modified one, and compared with each other. It is shown that the present scheme allows a sufficiently large Courant number of O(10$^{2}$) and reduces the computing time.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.397-403
• /
• 2010

The Immersed boundary method(IBM) is one of CFD techniques which can simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In these cases, pressure buildup near IB was found to occur when linear interpolation and stadard mass conservation is used and the interpolation scheme became complicated when higher order of interpolation is adopted. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. Bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies.

• Korean Journal of Mathematics
• /
• v.25 no.3
• /
• pp.405-435
• /
• 2017

We deal with a sensitivity analysis of an optimal shape control problem for the stationary Navier-Stokes system. A two-dimensional channel flow of an incompressible, viscous fluid is examined to determine the shape of a bump on a part of the boundary that minimizes the viscous drag. By using the material derivative method and adjoint variables for a shape sensitivity analysis, we derive the shape gradient of the design functional for the model problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.392-404
• /
• 1997

The flow induced by the oscillatory motion of a solid body is important in a number of practical problems. As the solid boundary oscillates harmonically, there is steady streaming motion invoked by the Reynolds stresses, which could cause extensive migration of the fluid during a period of fluid motion. We here analyzed the flow in a square cavity with an oscillating top wall for the parameters which make the time derivatives and the convective terms equally important in the entire cavity flow. The full Navier-Stokes equations are solved by the second-order time accurate Momentum Coupling Method which is devised by the authors. The particular numerical scheme does not need subiteration at each time step which is usually a required process to calculate the incompressible Navier-Stokes equations. The effect of two parameters, the Reynolds number and the frequency parameter, on the oscillatory flow has been investigated.

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

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