• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.3
• /
• pp.109-124
• /
• 2017

The numerical solution of the Stokes equation with discontinuous viscosity and singular force term is challenging, due to the discontinuity of pressure, non-smoothness of velocity, and coupled discontinuities along interface.In this paper, we give an efficient algorithm to solve this problem by employing spectral Legendre and Chebyshev approximations.First, we present the algorithm for a problem defined in rectangular domain with straight line interface. Then it is generalized to a domain with smooth curve boundary and interface by employing spectral element method. Numerical experiments demonstrate the accuracy and efficiency of our algorithm and its spectral convergence.

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

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler , thin layer Navier- Stokes and full Navier-Stokes ones, are solved using implicit LU-ADI decomposition scheme. The feasible direction algorithm with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• Journal of the Korean Mathematical Society
• /
• v.46 no.5
• /
• pp.1007-1025
• /
• 2009

First-order least-squares method of a distributed optimal control problem for the incompressible Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing the vorticity and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $H^{-1}$ and $L^2$ norms of the residual equations of the system. Finite element approximations are studied and optimal error estimates are obtained. Resulting linear system of the optimality system is symmetric and positive definite. The V-cycle multigrid method is applied to the system to test computational efficiency.

• 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.20 no.5
• /
• pp.671-680
• /
• 2007

The objective of this work is to reduce drag on a bluff body within a viscous flow by applying suction or injection of fluid along the surface of the body. In addition to minimizing drag, the optimal solution tends to reduce boundary layer separation and flow recirculation. When discretized by finite elements, the optimal control problem can be posed as a large-scale nonlinearly-constrained optimization problem. The constraints correspond to the discretized form of the Navier-Stokes equations. Unfortunately, solving such large-scale problems directly is essentially intractable. We developed several Sequential Quadratic Programming methods that are tailored to the structure of the control problem. Example problems of laminar flow around an infinite cylinder in two dimensions are solved to demonstrate the methodology. We use these optimal control techniques to study the influence of number of suction/injection holes and location of holes on the resulting optimized flow. We compare the proposed SQP methods against one another, as well as against available methods from the literature, from the point of view of efficiency and robustness. The most efficient of the proposed methods is two orders of magnitude faster than existing methods.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.193-198
• /
• 2006

It takes a lot of time and effort to generate grids for numerical analysis of problems with ground effect because the relative attitude and height of airfoil should be maintained to the ground as well as the inflow. A low Mach number preconditioned turbulent flow solver using the overlap grid technique has been developed and applied to the ground effect simulation. It has been validated that the present method using the multi-block grid gives us highly accurate solutions comparing with the experimental data of the RAE 101 airfoil in an unbounded condition. Present numerical method has been extended to simulate ground effect problems by using the overlapped grid system to avoid tedious work in generating multi-block grid system. An extended method using the overlapped grid has been verified and validated by comparing with results of multi-block method and experimental data as well. Consequently, the overlapped grid method can provide not only sufficiently accurate solutions but also the efficiency to simulate ground effect problems. It is shown that the pressure and aerodynamic centers move backward by the ground effect as the airfoil approaches to the ground.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.4
• /
• pp.37-47
• /
• 1999

A Finite Volume Method for the discretization of the two-dimensional incompressible Navier-Stokes equation is used to analyse water entry & exit problems in a generalized coordinate system. The free-surface deformations generated by the water entry or exit of a rigid body are simulated by the Level-Set scheme[11]. In the water entry problems for a wedged section and a flared-ship section, the calculation results of water impact force are compared with the experimental results[5] and the time varying free-surface deformations and flow characteristics of the water exit of a cylinder are investigated.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.10
• /
• pp.1139-1146
• /
• 2012

An alternative approach for topology optimization of steady incompressible Navier-Stokes flow problems is presented by using P1 nonconforming finite elements. This study is the extended research of the earlier application of P1 nonconforming elements to topology optimization of Stokes problems. The advantages of the P1 nonconforming elements for topology optimization of incompressible materials based on locking-free property and linear shape functions are investigated if they are also valid in fluid equations with the inertia term. Compared with a mixed finite element formulation, the number of degrees of freedom of P1 nonconforming elements is reduced by using the discrete divergence-free property; the continuity equation of incompressible flow can be imposed by using the penalty method into the momentum equation. The effect of penalty parameters on the solution accuracy and proper bounds will be investigated. While nodes of most quadrilateral nonconforming elements are located at the midpoints of element edges and higher order shape functions are used, the present P1 nonconforming elements have P1, {1, x, y}, shape functions and vertex-wisely defined degrees of freedom. So its implentation is as simple as in the standard bilinear conforming elements. The effectiveness of the proposed formulation is verified by showing examples with various Reynolds numbers.

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

Splitting algorithms of the incompressible Navier-Stokes equations using P1P1/P2P1 finite element formulation are newly proposed. P1P1 formulation allocates velocity and pressure at the same nodes, while P2P1 formulation allocates pressure only at the vertex nodes and velocity at both the vertex and mid nodes. For comparison of the elapsed time and accuracy of the two methods, they have been applied to the well-known benchmark problems. The three cases chosen are the two-dimensional steady and unsteady flows around a fixed cylinder, decaying vortex, and impinging slot jet. It is shown that the proposed P2P1 semi-splitting method performs better than the conventional P1P1 splitting method in terms of both accuracy and computation time.