• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.33-39
• /
• 1996

An incompressible flow stover based on the unstructured finite volume method has been developed. The flow domain is discretized by triangles in 2D or tetrahedra in 3D. The convective and viscous fluxes are obtained using edge connectivities of the unstructured 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. Laminar test flow problems are computed and presented with a comparison against other numerical solutions or experimental results.

• Journal of computational fluids engineering
• /
• v.17 no.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.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• 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.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.71-80
• /
• 1996

An iterative time marching procedure for solving incompressible turbulent flow has been applied to the flows around a high speed train including cross-wind effects. This procedure solves three-dimensional unsteady incompressible Reynolds-averaged Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using first-order accurate schemes for the time derivatives and third/second-order accurate schemes for the spatial derivatives. Turbulent flows have been modeled by Baldwin-Lomax turbulent model. To validate present procedure, the flow around a high speed train at zero yaw angle was simulated and compared with experimental data. Generally good agreement with experiments was achieved. The flow fields around the high speed train at 9.2°, 16.7°, and 45° of yaw angle were also simulated.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.55-63
• /
• 1995

The flow field around a high-speed train including cross-wind effects has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations in the inertial frame using the iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives, 3th-order accurate QUICK scheme for the convective terms and 2nd-order accurate central difference scheme for the viscous terms. The Marker-and-Cell concept was applied to efficiently solve continuity equation, which is differenced with 2nd-order accurate central difference scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A C-H type of elliptic grid system is generated around a high-speed train including ground. The Baldwin-Lomax turbulent model was implemented to simulate the turbulent flows. To validate the present procedure, the flow around a high speed train at constant yaw angle of $45^{\circ}\;and\;90^{\circ}$ has been simulated. The simulation shows 3-D vortex generation in the lee corner. The flow separation is also observed around the rear of the train. It has concluded that the results of present study properly agree with physical flow phenomena.

• The KSFM Journal of Fluid Machinery
• /
• v.2 no.2 s.3
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.414-426
• /
• 1997

Numerical computations were conducted to characterize the three-dimensional laminar flow through an injector orifice having an inclined angle of 30 .deg.. For this study, the incompressible Navier-Stokes equations in generalized curvilinear coordinates, using a pseudocompressibility approach for continuity equation, were solved. The computations were performed using the finite difference implicit, approximately factored scheme of Beam and Warming and multi-block grids of complete continuity at block interfaces. The multi-block computations were validated for the steady state using direct comparison of multi-block solutions with equivalent single-block ones, including 2-D 180.deg. TAD and 3-D 90.deg. pipe bend. The comparisons between the numerical solutions and the flow field measurements for a tube with sudden contraction were presented in this work for solution validation. Computational results showed the nature of complex flow fields within the inclined injector orifice, including strong pressure-driven secondary flows in the cross stream induced by the effect of streamline curvature. In addition, asymmetric secondary flows were induced in the Reynolds number range above assumed laminar flow regime considered. However, turbulence calculations and grid dependency studies are needed for more accurate computations.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.153-158
• /
• 2005

Numerical simulations of two-dimensional steady incompressible lid-driven flow in a square cavity are presented to verify the validity of a new solution code(PowerCFD) with unstructured grids. The code uses the non-staggered(collocated) grid approach which is very popular for incompressible flow analysis because of its numerical efficiency on the curvilinear or unstructured grids. Solutions are obtained for configurations with a Reynolds number as high as 10,000 with both rectangular and hybrid types of unstructured grid mesh. Interesting features of the flow are presented in detail and comparisons are made with benchmark solutions found in the literature. It is found that the code is capable of producing accurately the nature of the lid-driven cavity flow at high Reynolds numbers.

• Journal of applied mathematics & informatics
• /
• v.41 no.2
• /
• pp.401-411
• /
• 2023

The effects of longitudinal velocity dusty fluid flow in a weak magnetic field are investigated in this paper. An external uniform magnetic field parallel to the flow of dusty fluid influences the flow of dusty fluid. Besides that, the problem under investigation is completely defined in terms of identifying parameters such as longitudinal velocity (u), Hartmann number (M), dust particle interactions β, stock resistance γ, Reynolds number (Re) and magnetic Reynolds number (Rm). While using suitable transformations of resemblance, The governing partial differential equations are transformed into a system of ordinary differential equations. The Hankel Transformation is used to solve these equations numerically. The effects of representing parameters on the fluid phase and particle phase velocity flow are investigated in this analysis. The magnitude of the fluid particle is reduced significantly. The result indicates the magnitude of the particle reduced significantly. Although some of our numerical solutions agree with some of the available results in the literature review, other results differs because of the effect of the introduced magnetic field.