• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1999년도 춘계 학술대회논문집
• /
• pp.98-105
• /
• 1999

The Pressure-based Unstructured-grid Navier-Stokes Solver PUNS-2/3D for incompressible steady and unsteady viscous flows has been developed. It is based on nonstaggered cell-centered finite volume method. Second-order upwind scheme with least-square reconstruction is used for convective fluxes. The SIMPLE method is implemented to couple the pressure and velocity fields. And the time derivatives in the momentum equations are discretised using a second-order Euler backward-differencing scheme. The discretised linear equations are solved by the preconditioned Biconjugate Gradient Stabilized method(Bi-CGSTAB). The developed solver is applied to validation problems using hybrid meshes.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1999년도 유체기계 연구개발 발표회 논문집
• /
• pp.163-169
• /
• 1999

The numerical procedure has been developed for simulating incompressible viscous flow around a turbine stage with rotor-stator interaction. This study solves 2-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system. The Marker-and-Cell concept is applied to efficiently solve continuity equation. To impose an accurate boundary condition, O-H multiblocked grid system is generated. O-type grid and H-type grid is generated near and outer rotor-stator The cubic-spline interpolation is applied to handle a relative motion of a rotor to the stator. Turbulent flows have been modeled by the Baldwin- Lomax turbulent model. To validate present procedure, the time averaged pressure coefficients around the rotor and stator are compared with experiment and a good agreement obtained.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2012년도 제38회 춘계학술대회논문집
• /
• pp.187-195
• /
• 2012

본 연구에서는 터보펌프 인듀서 주변의 극저온 유동장에 대한 수치해석을 주로 다룬다. 이를 위해 수치해석에 사용된 수치기법에 대한 내용을 우선 다룬다. 수치기법에 대한 검증 문제로 2차원/축대칭 극저온 유동문제를 해석하였고 이 과정에서 다양한 공동모델에 대한 비교 연구를 수행하였다. 마지막으로 3차원 터보펌프 인듀서에 대한 수치해석을 수행하고, 실험 결과와 타 연구결과와 비교함과 동시에 극저온 유체에서 나타나는 공동현상의 특징을 살펴보았다.

• Korean Journal of Mathematics
• /
• 제27권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.

• 한국음향학회:학술대회논문집
• /
• 한국음향학회 1996년도 제11회 수중음향학 학술발표회 논문집 11th Underwater Acoustics Symposium Proceedings
• /
• pp.59-62
• /
• 1996

A numerical study of a two-dimensional acoustic field is carride ort for a spinning vortex pair located neat a wall to investigate the effect of the wall from the spinning quadrupole source in unsteady vortical flows. Based on the known incompressible flow field, the perturbed compressible acoustic terms derived from the Euler equations are calculated. Non-reflecting boundary conditions on the free field and the solid boundary conditions are developed for a generalized curvilinear coordinates system to investigate the effect of a curced wall. It is concluded that the sound generated by the quadrupole sources of unsteady vortical flows in the presence of a flat wall or a circular cylinder can be calculated by using the source terms of hydrodynamic flow fluctuations in both near and far acoustic fields simultaneously.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.15-16
• /
• 2003

An implicit pressure correction method on unstructured Cartesian grid is developed for the incompressible Navier-Stokes equations. An immersed boundary method is also incorporated to treat the body geometry. Tests show that with an appropriate amount of dissipation, the method is second order accurate both in time and space. The driven cavity flows with and without immersed bodies are computed to demonstrate the capability of the present scheme.

• 한국해양공학회지
• /
• 제17권5호
• /
• pp.11-18
• /
• 2003

In two-dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure (HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint hence, achieving an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared (CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper and efficient iterative procedure for higher order finite element formulation has not been available, thus far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient, and useful in flow analyses, in conjunction with hierarchical elements.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제19권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.

• 한국해양공학회지
• /
• 제25권5호
• /
• pp.9-14
• /
• 2011

The divergence-free finite elements introduced in this paper are derived from Hermite functions, which interpolate stream functions. Velocity bases are derived from the curl of the Hermite functions. These velocity basis functions constitute a solenoidal function space, and the gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into its solenoidal and irrotational parts, and the decoupled Navier-Stokes equations are then projected onto their corresponding spaces to form appropriate variational formulations. The degrees of the Hermite functions we introduce in this paper are bi-cubis, quartic, and quintic. To verify the accuracy and convergence of the present method, three well-known benchmark problems are chosen. These are lid-driven cavity flow, flow over a backward facing step, and buoyancy-driven flow within a square enclosure. The numerical results show good agreement with the previously published results in all cases.