• Journal of Korea Water Resources Association
• /
• v.45 no.6
• /
• pp.545-555
• /
• 2012

A numerical modeling has become increasingly popular and more important to the study of water waves with a rapid advancement of computer technology. However, different types of problems are induced during simulating wave motion. One of the key problems is re-reflection to a computation domain at the incident boundary. The internal wave generating-absorbing boundary conditions have been commonly used in numerical wave models to prevent re-reflection. For the Navier-Stokes equations model, the internal wave maker using a mass source function of the continuity equation has been used to generate various types of waves. Nonetheless, almost every numerical experiment is performed in two dimensions and only a few tests have been expanded to three dimensions. More recently, a momentum source function of the Boussinesq equations is applied to generate essentially directional waves in the three dimensional Navier-Stokes equations model. In this study, the internal wave maker using a momentum source function is employed to generate targeted linear waves in the three-dimensional LES model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.9
• /
• pp.1325-1334
• /
• 1998

The objectives of the present study are to develop a systematic method rather than a conventional trial-and-error method for an optimal shape design using a mathematical theory, and to apply it to engineering problems. In the present study, an optimal condition for a minimum pressure loss in a two-dimensional curved duct flow is derived and then an optimal shape of the curved duct is designed from the optimal condition. In the design procedure, one needs to solve the adjoint Navier-Stokes equations which are derived from the Navier-Stokes equations and the cost function. Therefore, a computer code of solving both the Navier-Stokes and adjoint Navier-Stokes equations together with an automatic grid generation is developed. In a curved duct flow, flow separation occurs due to an adverse pressure gradient, resulting in an additional pressure loss. Optimal shapes of a curved duct are obtained at three different Reynolds numbers of 100, 300 and 800, respectively. In the optimally shaped curved ducts, the separation region does not exist or is significantly reduced, and thus the pressure loss along the curved duct is significantly reduced.

• Journal of the Korean Mathematical Society
• /
• v.57 no.5
• /
• pp.1239-1266
• /
• 2020

In this paper, some novel discrete formulations for stabilizing the mixed finite element method Q1-Q0 (bilinear velocity and constant pressure approximations) are introduced and discussed for the generalized Stokes problem. These are based on stabilizing discontinuous pressure approximations via local jump operators. The developing idea consists in a reduction of terms in the local jump formulation, introduced earlier, in such a way that stability and convergence properties are preserved. The computer implementation aspects and numerical evaluation of these stabilized discrete formulations are also considered. For illustrating the numerical performance of the proposed approaches and comparing the three versions of the local jump methods alongside with the global jump setting, some obtained results for two test generalized Stokes problems are presented. Numerical tests confirm the stability and accuracy characteristics of the resulting approximations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.4
• /
• pp.458-465
• /
• 2003

The main objective of the research is to develop a research code solving transient incompressible Navier-Stokes equation. In this research code, Adams-Bashforth method was applied to the convective terms of the navier stokes equation and the splitted equations were discretized spatially by finite element methods to solve the complex geometry problems easily. To reduce the divergence on the boundaries of pressure poisson equation due to the unsuitable pressure boundary conditions, multi step approximation pressure boundary conditions derived from the boundary linear momentum equations were used. Simulations of Lid Driven Flow and Flow over Cylinder were conducted to prove the accuracy by means of the comparison with results of the previous workers.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.103-121
• /
• 2015

In this paper, we briefly review and describe a projection algorithm for numerically computing the two-dimensional time-dependent incompressible Navier-Stokes equation. The projection method, which was originally introduced by Alexandre Chorin [A.J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comput., 22 (1968), pp. 745-762], is an effective numerical method for solving time-dependent incompressible fluid flow problems. The key advantage of the projection method is that we do not compute the momentum and the continuity equations at the same time, which is computationally difficult and costly. In the projection method, we compute an intermediate velocity vector field that is then projected onto divergence-free fields to recover the divergence-free velocity. Numerical solutions for flows inside a driven cavity are presented. We also provide the source code for the programs so that interested readers can modify the programs and adapt them for their own purposes.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.115-138
• /
• 2019

The dual singular function method [DSFM] is a solver for corner sigulaity problem. We already construct DSFM in previous reserch to solve the Stokes equations including one singulairity at each reentrant corner, but we find out a crucial incorrection in the proof of well-posedness and regularity of dual singular function. The goal of this paper is to prove accuracy and well-posdness of DSFM for Stokes equations including two singulairities at each corner. We also introduce new applicable algorithms to slove multi-singulrarity problems in a complicated domain.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.54-61
• /
• 2015

In the present study, a flow solver using a kinetic BGK scheme was developed for the compressible Navier-Stokes equation. The kinetic BGK scheme was used to simulate flow field from the continuum up to the transitional regime, because the kinetic BGK scheme can take into account the statistical properties of the gas particles in a non-equilibrium state. Various numerical simulations were conducted by the present flow solver. The laminar flow around flat plate and the hypersonic flow around hollow cylinder of flare shape in the continuum regime were numerically simulated. The numerical results showed that the flow solver using the kinetic BGK scheme can obtain accurate and robust numerical solutions. Also, the present flow solver was applied to the hypersonic flow problems around circular cylinder in the transitional regime and the results were validated against available numerical results of other researchers. It was found that the kinetic BGK scheme can similarly predict a tendency of the flow variables in the transitional regime.

• Journal of the Korean Institute of Navigation
• /
• v.9 no.2
• /
• pp.43-63
• /
• 1985

As a circular cylinder has a comparatively simple shape and becomes a basic problem for flows around other various shapes of bodies, the problem of two-dimensional viscous flow around the circular cylinder has been investigated, both theoretically and experimentally. But not a few problems are left unsolved. It is well known that the calculations are successfully made with the approximations of Stokes or Oseen for very low Reynolds numbers, but as Reynolds number is increased, Oseen's approximations as well as Stokes's ones become more and more remote from the exact solution of the Navier-Stokes equations. Therefore, in this paper, the authors transform the Navier-Stokes equations into the finite difference equations in the steady two-dimensional viscous flow at Reynolds number up to 45, and then solve the solution of the Navier-Stokes equations numerically. Also, the authors examine the accuracy of the solution by means of flow visualization with aluminum powder. The main results are as follows; (1) The critical Reynolds number at which twin vortices begin to form in the rear of the circular cylinder is found to be 6 in the experiment and 4 in the numerical solution. (2) As Reynolds number is increased, it is proved that the ratio of the length of the twin vortices to the diameter is grown almost linearly, both experimentally and numerically. (3) Separation angle is also increased according to reynolds number. But it is found that it would converge into 101.3 degrees, both experimentally and numerically.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.67-72
• /
• 2005

This research is based on the iterative space-marching method for incompressible and compressible Navier-Stokes equations[1-4]. A principle of parallel computational schemes construction for steady and unsteady problems is suggested. It is analytically proven that convergence of these schemes is unconditional for incompressible case. When the parallel scheme is used the total volume of computations is the sum of a large number of independent and equal parts. Estimation of the speed-up K shows that K > 1000 in ideal case. First results of using the parallel schemes are presented.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.49-55
• /
• 2009

This paper evaluates performances of a recently developed divergence-free finite element method based on Hermite interpolated stream functions. Velocity bases are derived from Hermite interpolated stream functions to form divergence-free basis functions. These velocity basis functions constitute a solenoidal function space, and the simple gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into a solenoidal and an irrotational parts, and the decoupled Navier-Stokes equations are projected onto their corresponding spaces to form proper variational formulations. To access accuracy and convergence of the present algorithm, three test problems are selected. They are lid-driven cavity flow, flow over a backward-facing step and buoyancy-driven flow within a square enclosure. Hermite interpolation functions from cubic to quintic are chosen to run the test problems. Numerical results are shown. In all cases it has shown that the present method has performed well in accuracies and convergences. Moreover, the present method does not require an upwinding or a stabilized term.