• 한국항해항만학회지
• /
• 제27권5호
• /
• pp.473-478
• /
• 2003

평행한 벽 사이에 원형배관(circular pipeline)을 놓고 그 주위의 유동특성에 대한 수치연구를 수행하였다. 비압축성 유체를 가지고, Navier-Stokes 방정식을 풀었고 3차 풍상(upwind) 차분의 수치해법을 이용하였다. 한쪽 벽과의 거리가 매우 작아질 때, 볼텍스 떨어짐이 상당히 억압되는데 이것은 벽 경계와의 상호 박리 작용 때문으로 간수된다. 본 연구는 레이놀드 수의 변함과 물체가 벽에 접근함에 따른 볼텍스 떨어짐의 특성을 규명하는데 있다. 원형배관 후류와 평행벽내 유기된 박리의 상호작용을 집중적으로 다루며 서로 다른 조건에서 박리와 와역한(vorticity dynamics)의 특성을 해석하였다.

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

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• 한국전산유체공학회지
• /
• 제4권3호
• /
• pp.1-11
• /
• 1999

This paper describes the flow field analysis of an automobile with crosswind effects of 15°, 30° 45° and 60° of yaw angles. The governing equations of the 3-D incompressible Navier-Stokes equations are solved by the iterative time marching scheme. The Chimera grid technique has been applied to efficiently simulate the flow around the side-view mirror. The computated surface pressure coefficients have been compared with experimental results and a good agreement has been achieved. The A- and C-pillar vortex and other flow phenomena around the ground vehicle are evidently shown. The variation of aerodynamic coefficients of drag, lift, side force and moments with respect to yaw angle is systematically studied.

• 대한기계학회논문집B
• /
• 제24권10호
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

• 대한조선학회논문집
• /
• 제41권4호
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.

• 대한기계학회논문집
• /
• 제19권11호
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• Journal of Mechanical Science and Technology
• /
• 제20권10호
• /
• pp.1730-1740
• /
• 2006

A nodeless variables finite element method for analysis of two-dimensional, steady-state viscous incompressible flow is presented. The finite element equations are derived from the governing Navier-Stokes differential equations and a corresponding computer program is developed. The proposed method is evaluated by solving the examples of the lubricant flow in journal bearing and the flow in the lid-driven cavity. An adaptive meshing technique is incorporated to improve the solution accuracy and, at the same time, to reduce the analysis computational time. The efficiency of the combined adaptive meshing technique and the nodeless variables finite element method is illustrated by using the example of the flow past two fences in a channel.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2006년 제4회 한국유체공학학술대회 논문집
• /
• pp.349-352
• /
• 2006

The conventional segregated finite element formulation produces a small and simple matrix at each step than in an integrated formulation. And the memory and cost requirements of computations are significantly reduced because the pressure equation for the mass conservation of the Navier-Stokes equations is constructed only once if the mesh is fixed. However, segregated finite element formulation solves Poisson equation of elliptic type so that it always needs a pressure boundary condition along a boundary even when physical information on pressure is not provided. On the other hand, the conventional integrated finite element formulation in which the governing equations are simultaneously treated has an advantage over a segregated formulation in the sense that it can give a more robust convergence behavior because all variables are implicitly combined. Further it needs a very small number of iterations to achieve convergence. However, the saddle-paint-type matrix (SPTM) in the integrated formulation is assembled and preconditioned every time step, so that it needs a large memory and computing time. Therefore, we newly proposed the P2PI semi-segregation formulation. In order to utilize the fact that the pressure equation is assembled and preconditioned only once in the segregated finite element formulation, a fixed symmetric SPTM has been obtained for the continuity constraint of the present semi-segregation finite element formulation. The momentum equation in the semi-segregation finite element formulation will be separated from the continuity equation so that the saddle-point-type matrix is assembled and preconditioned only once during the whole computation as long as the mesh does not change. For a comparison of the CPU time, accuracy and condition number between the two methods, they have been applied to the well-known benchmark problem. It is shown that the newly proposed semi-segregation finite element formulation performs better than the conventional integrated finite element formulation in terms of the computation time.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.404-409
• /
• 2001

An account of second-order fractional-step methods and boundary conditions for the incompressible Navier-Stokes equations is presented. The present work has aimed at (i) identification and analysis of all possible splitting methods of second-order splitting accuracy; and (ii) determination of consistent boundary conditions that yield second-order accurate solutions. It has been found that only three types (D, P and M) of splitting methods called the canonical methods are non-degenerate so that all other second-order splitting schemes are either degenerate or equivalent to them. Investigation of the properties of the canonical methods indicates that a method of type D is recommended for computations in which the zero divergence is preferred, while a method of type P is better suited to the cases when highly-accurate pressure is more desirable. The consistent boundary conditions on the tentative velocity and pressure have been determined by a procedure that consists of approximation of the split equations and the boundary limit of the result. The pressure boundary condition is independent of the type of fractional-step methods. The consistent boundary conditions on the tentative velocity were determined in terms of the natural boundary condition and derivatives of quantities available at the current timestep (to be evaluated by extrapolation). Second-order fractional-step methods that admit the zero pressure-gradient boundary condition have been derived. The boundary condition on the new tentative velocity becomes greatly simplified due to improved accuracy built in the transformation.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.55-58
• /
• 2002

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.