• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1205-1219
• /
• 2011

In this work, a stabilized characteristic finite volume method for the time-dependent Navier-Stokes equations is investigated based on the lowest equal-order finite element pair. The temporal differentiation and advection term are dealt with by characteristic scheme. Stability of the numerical solution is derived under some regularity assumptions. Optimal error estimates of the velocity and pressure are obtained by using the relationship between the finite volume and finite element methods.

• 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.17 no.1
• /
• pp.44-53
• /
• 2012

Immersed boundary method(IBM) is a numerical scheme proposed to 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 this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. The weight coefficients of the 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. For an analysis of moving boundary, we smulated an oscillating circular cylinder with Re=100 and KC(Keulegan-Carpenter) number of 5. The predicted flow fields were compared with experimental data and they also showed good agreements.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.51-67
• /
• 1999

The numerical solution of the Navier-Stokes equations in a constricted stepped channel problem has been obtained using a fourth order numerical method. Transformations are made to have a fine grid near the sharp corner and a long channel downstream. The derivatives in the Navier-Stokes equations are replaced by fourth order central differences which result a 29-point computational stencil. A procedure is used to avoid extra numerical boundary conditions near the solid walls. Results have been obtained for Reynolds numbers up to 1000.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1331-1345
• /
• 2010

This paper describes a numerical solution of Navier-Stokes equations. It includes algorithms for discretization by finite element methods and 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.21 no.2
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.

• Bulletin of the Korean Mathematical Society
• /
• v.48 no.5
• /
• pp.1079-1092
• /
• 2011

We consider the Newton's method for an direct solver of the optimal control problems of the Navier-Stokes equations. We show that the finite element solutions of the optimal control problem for Stoke equations may be chosen as the initial guess for the quadratic convergence of Newton's algorithm applied to the optimal control problem for the Navier-Stokes equations provided there are sufficiently small mesh size h and the moderate Reynold's number.

• Bulletin of the Korean Mathematical Society
• /
• v.50 no.5
• /
• pp.1737-1751
• /
• 2013

An efficient and stable point collocation scheme based on a meshfree method is studied for the stationary incompressible Navier-Stokes equations. We describe the diffuse derivatives associated with the moving least square method. Using these diffuse derivatives, we propose a point collocation method to fit in solving the Navier-Stokes equations which improves the stability of the direct point collocation scheme. The convergence of the numerical solution is investigated from numerical examples. The driven cavity ow and the backward facing step ow are implemented for the reliability of the scheme. Also, the viscous ow on complicated geometry is successfully calculated such as the ow past a circular cylinder in duct.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

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