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

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

• 대한기계학회논문집B
• /
• 제31권9호
• /
• pp.814-821
• /
• 2007

A hybrid method is presented for efficient computation of turbulent flow noise at low Mach numbers. In this method, the turbulent flow field is computed by incompressible large eddy simulation (LES), while the acoustic field is computed with the linearized perturbed compressible equations (LPCE) derived in this study. Since LPCE is computed on the rather coarse acoustic grid with the flow variables and source term obtained by the incompressible LES, the computational efficiency of calculation is greatly enhanced. Furthermore, LPCE suppress the instability of perturbed vortical mode and therefore secure consistent and stable acoustic solutions. The proposed LES/LPCE hybrid method is applied to three low Mach number turbulent flow noise problems: i) circular cylinder, ii) isolated flat plate, and iii) interaction between cylinder wake and airfoil. The computed results are closely compared with the experimental measurements.

• 한국전산유체공학회지
• /
• 제12권4호
• /
• pp.20-27
• /
• 2007

Two-dimensional stagnation flow toward a plane wall coated with magnetic fluid of uniform thickness is investigated. The flow field is represented as a similarity solution of the Navier-Stokes equation for this incompressible laminar flow. The resulting third order ordinary differential equation is solved numerically by using the shooting method and by determining two shooting parameters so as to satisfy the boundary and interface conditions. Features of the flow including streamline patterns are investigated for the varying values of density ratio, viscosity ratio, and Reynolds number. An adverse flow with double eddy pair in magnetic fluid region is found to emerge as the Reynolds number becomes higher than a threshold value. The results for the interface velocity, interface and wall shear stress, and boundary layer and displacement thickness are also presented.

• 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.

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

Substructuring methods are usually used in finite element structural analyses. In this study a multi-level substructuring algorithm is developed and proposed as a possible candidate for incompressible fluid solves. Finite element formulation for incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et.al.[5]. The present algorithm consists of four stages such as a gathering stage, a condensing stage, a solving stage and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At highest level, each subdomain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each subdomain has been replaced by a sequential static condensation. The global algebraic system arising feom the assembly of each subdomains is solved using Conjugate Gradient Squared(CGS) method. In this case, pre-conditioning techniques usually accompanied by iterative solvers are not needed.

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

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

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

In this paper, the new higher order wall boundary conditions are proposed for solving the incompressible flows. The square driven cavity flows are simulated by using the variable order method of lines with the present wall boundary conditions. The variable order method of lines is constructed by the spatial discretization, i.e., the variable order proper convective scheme for convective terms and the modified differential quadrature method for diffusive terms, and time integration. The 2nd, 4th, and 6th order solutions are presented and these results show this higher order boundary conditions are very promising for the incompressible flow simulations.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.289-294
• /
• 2010

Numerical stud of an oscillating body in incompressible fluid is performed. Stabilized finite element method comprising of Streamline-Upwind/Petrov-Galerkin (SUPG) and Pressure-Stabilizing/Petrov-Galerkin (PSPG) formulations of linear triangular elements was employed to solve 2D incompressible Navier-Stokes equations whereas the motion of the body was considered by incorporating the arbitrary Langrangian-Eulerian(ALE) formulation. An algebraic moving mesh strategy is utilized for obtaining body conforming mesh deformation at each time step. Two tests cases, namely motion of a circular cylinder and of an airfoil in incompressible flow were analyzed. The model is first validated against the stationary cases and then the capability to handle moving boundaries is demonstrated.

• 한국전산유체공학회지
• /
• 제14권4호
• /
• pp.93-100
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved successively in the moving least square sense. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as a FVM.

• 한국유체기계학회 논문집
• /
• 제2권2호
• /
• pp.57-63
• /
• 1999

The flow through multistage turbomachinery is affected by the interaction between a rotor and a stator. The interaction is due to the inviscid potential effect and viscous effect between closely spaced rotor and stator airfoils. Three-dimensional, unsteady, incompressible Navier-Stokes equations with a standard $k-{\epsilon}$ model are solved using a non-staggered grid system. This method is applied to the flow through a multistage compressor measured by Stauter et al. The results have shown strong interaction between the rotating and stationary flow field. The decay of rotor wake and the pressure profiles agree very well with experimental data. The wake produced by rotor causes unsteady pressure on the surface of a stator. The rotor/stator interaction produces the unsteady pressure force on the rotor and stator blades.