• Journal of the Korean Mathematical Society
• /
• v.36 no.1
• /
• pp.173-192
• /
• 1999

We study a boundary optimal control problem of the fluid flow governed by the Navier-Stokes equations. the control problem is formulated with the flow through a channel with steps. The first-order optimality condition of the optimal control is derived. Finite element approximations of the solutions of the optimality system are defined and optimal error estimates are derived. finally, we present some numerical results.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.10-18
• /
• 2016

Numerical boundary conditions are as important as the governing equations when analyzing the fluid flows numerically. An explicit boundary condition method updates the solutions at the boundaries with extrapolation from the interior of the computational domain, while the implicit boundary condition method in conjunction with an implicit time integration method solves the solutions of the entire computational domain including the boundaries simultaneously. The implicit boundary condition method, therefore, is more robust than the explicit boundary condition method. In this paper, steady compressible 2-Dimensional Navier-Stokes solver is developed. We present the implicit boundary condition method coupled with LU-SGS(Lower Upper Symmetric Gauss Seidel) method. Also, the explicit boundary condition method is implemented for comparison. The preconditioning Navier-Stokes equations are solved on unstructured meshes. The numerical computations for a number of flows show that the implicit boundary condition method can give accurate solutions.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.111-116
• /
• 2005

Two coupling methods for the Navier-Stokes equations and a two-equation turbulence model equations are compared. They are the strongly coupled method and the loosely coupled method. The strongly coupled method solves the Navier-Stokes equations and the two-equation turbulence model equations simultaneously, while the loosely coupled method solves the Navier-Stokes equation with the turbulence viscosity fixed and subsequently solves the turbulence model equations with all the flow quantities fixed. In this paper, performances of two coupling methods are compared for two and three-dimensional problems.

• Journal of the Korean Chemical Society
• /
• v.40 no.6
• /
• pp.409-414
• /
• 1996

The limit of applicability of Navier-Stokes equation to stationary plane shock-waves is examined by using the principle of minimum entropy production of linear irreversible thermodynamics. In order to obtain analytic results, the equation is linearized near the equilibrium of downstream. Results show that the solution of Navier-Stokes equation which fits the boundary condition of far downstream flow is consistent with the thermodynamic requirement within the first order when the solution is expanded around the M=1, where M is the Mach number of upstream speed.

• Journal of the Korean Mathematical Society
• /
• v.50 no.4
• /
• pp.771-795
• /
• 2013

In this paper we consider the Navier-Stokes equations in the half space. Our aim is to construct a mild solution for initial data in $B^{-\alpha}_{{\infty},{\infty}}(\mathbb{R}^n_+)$, 0 < ${\alpha}$ < 1. To do this, we derive the estimate of the Stokes flow with singular initial data in $B^{-\alpha}_{{\infty},q}(\mathbb{R}^n_+)$, 0 < ${\alpha}$ < 1, 1 < $q{\leq}{\infty}$.

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.675-691
• /
• 2002

The focus of this work is on the development of large-scale numerical optimization methods for optimal control of steady incompressible Navier-Stokes flows. The control is affected by the suction or injection of fluid on portions of the boundary, and the objective function of fluid on portions of the boundary, and the objective function represents the rate at which energy is dissipated in the fluid. We develop reduced Hessian sequential quadratic programming. Both quasi-Newton and Newton variants are developed and compared to the approach of eliminating the flow equations and variables, which is effectively the generalized reduced gradient method. Optimal control problems we solved for two-dimensional flow around a cylinder. The examples demonstrate at least an order-of-magnitude reduction in time taken, allowing the optimal solution of flow control problems in as little as half an hour on a desktop workstation.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.90-95
• /
• 2011

The Navier-Stokes equations are solved to study the flow characteristics of a micro viscous pump. The viscous micropump is consisted of a stationary disk with a spiral shaped channel and a rotating disk. A simple geometrical model for the tip clearance is proposed and validated by comparing computed flow rate with corresponding experimental data. Present numerical solutions show satisfactory agreement with the corresponding experimental data. The tip clearance effect is found to become significant as the rotational speed increases. As the pressure load increases, a reversed flow region is seen to form near the stationary disk. The height of the channel is shown to be optimized in terms of the flow rate for a given rotational speed and pressure load. The optimal height of the channel becomes small as the rotational speed decreases or the pressure load increases. The flow rate of the pump is found to be in proportion to the width of channel.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.11 no.1
• /
• pp.31-39
• /
• 2007

A pedagogic model for blood flow is introduced to help medicine majors understand a simplified version of Navier-Stokes equations which is known to be a good tool for interpreting the phenomena in blood flow. The pressure gradient consists of a time-independent part known as Hagen-Poiseuille's gradient and a time-dependent part known as Sexl's, and the model formula for the volume rate of blood flow is reduced to a very simple form. For demonstration, the blood rate in human aorta system is analyzed in connection with the time-dependence of pressure gradient. It is shown for Sexl's part that the flow rate lags the pressure gradient by ${\pi}/2$, which is thought to be due to the relaxation process involved.