• 한국전산유체공학회지
• /
• 제17권2호
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• 한국전산유체공학회지
• /
• 제1권1호
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• 대한기계학회논문집
• /
• 제12권3호
• /
• pp.575-581
• /
• 1988

본 연구에서는 앞서의 '압력의 단가성' 개념을 재도입하여 유동함수의 경계치 설정에 수치적으로 유용하게 쓸 수 있음을 보였고, 이를 다중연결된 물체에 관한 2차 원 응용문제들에서 수렴성이 나쁘고 많은 계산시간을 요구하는 종래의 원시변수들을 이용한 수식화 과정을 벗어나, 이제는 간편하고 명확한 유동함수-와도의 방법으로 비 교적 적은 컴퓨터 시간으로 유동장의 계산이 가능하게 되었음을 보였다.

• 대한기계학회논문집B
• /
• 제30권3호
• /
• pp.262-269
• /
• 2006

Segregation algorithms of the incompressible Wavier-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 the midpoint nodes. For a comparison of both the elapsed time and the accuracy between 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-segregation algorithm performs better than the conventional P1P1 segregation algorithm in terms of both accuracy and computation time.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권2호
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• 한국유체기계학회 논문집
• /
• 제14권3호
• /
• pp.39-44
• /
• 2011

In this study, performance analyses have been conducted for a 5MW class wind turbine blade model. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate detailed dynamic responsed of wind turbine blade. Reynolds-averaged Navier-Stokes (RANS) equations with K-${\epsilon}$ turbulence model are solved for unsteady flow problems of the rotating turbine blade model. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3D turbine blade for fluid-structure interaction (FSI) problems. Predicted aerodynamic performance considering structural deformation effect of the blade show different results compared to the case of rigid blade model.

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

To investigate the flow inside the centrifugal impeller, computer program which can solve Three-dimensional compressible turbulent flow has been developed. The Navier-Stokes equations were chosen as the governing equation for viscous flow while Euler equations for inviscid case. Time marching method was incorporated with the Flux Difference Splitting method suggested by Roe to capture the steep gradients such as a shock. For high order of accuracy, MUSCL approach was adopted while differentiable limiter to ensure TVD property. For turbulence closure, Baldwin- Lomax model was applied due to its simplicity. To demonstrate the capabilities of present program, several validation problems have been solved and compared with experiments and other available data. From the above calculations generally good agreements were obtained. Finally, the developed code was applied to Eckardt's impeller and the performance prediction was carried out. Some important aspects on boundary condition for successful simulation were discussed and the remedy was also introduced.

• 대한수학회지
• /
• 제55권2호
• /
• pp.471-506
• /
• 2018

In this article, some projection methods (or fractional-step methods) are proposed and analyzed for the micropolar Navier-Stokes equations (MNSE). These methods allow us to decouple the MNSE system into two sub-problems at each timestep, one is the linear and angular velocities system, the other is the pressure system. Both first-order and second-order projection methods are considered. For the classical first-order projection scheme, the stability and error estimates for the linear and angular velocities and the pressure are established rigorously. In addition, a modified first-order projection scheme which leads to some improved error estimates is also proposed and analyzed. We also present the second-order projection method which is unconditionally stable. Ample numerical experiments are performed to confirm the theoretical predictions and demonstrate the efficiency of the methods.

• 대한토목학회논문집
• /
• 제13권5호
• /
• pp.53-65
• /
• 1993

본 연구에서는 강성기초의 동적응답을 얻기 위해서 비완화 경계조건(non-relaxed boundary condition)을 적용한 3차원 경계요소를 사용하였다. 경계요소는 장래의 비선형문제의 확장을 위해서 시간영역에서 형식화되었으며 기본해는 무한영역의 Stokes 해를 사용하였다. 본 연구는 검증되었으며 지반기초 및 임의 형상의 지하구조물외 동적응답을 얻는데 이용할 수 있다.

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

A two-dimensional Navier-Stokes solver using hybrid meshes is parallelized with a domain decompostion method. The focus of this paper is placed on minimizing the amount of effort in parallelizing the serial version of the solver, and this is achieved by adding an additional layer of cells to each decomposed domain. Most subroutines of the serial solver are used without modification, and the information exchange between neighboring domains is achieved using MPI(Message Passing Interface) library. Load balancing among the processors and scheduling of the message passing are implemented to reduce the overhead of parallelization, and the speed-up achieved by parallelization is measured on the transonic invisicd and turbulent flow problems. The parallelization efficiencies of the explicit Runge-Kutta scheme and the implicit point-SGS scheme are compared and the effects of various factors on the results are also studied.