Boundary conformed co-ordinate systems for selected two-dimensional fluid flow problems. Part I: Generation of BFGs

In this paper the generation of general curvilinear co-ordinate systems for use in selected two-dimensional fluid flow problems is presented. The curvilinear co-ordinate systems are obtained from the numerical solution of a system of Poisson equations. The computational grids obtained by this technique allow for curved grid lines such that the boundary of the solution domain coincides with a grid line. Hence, these meshes are called boundary fitted grids (BFG). The physical solution area is mapped onto a set of connected rectangles in the transformed (computational) plane which form a composite mesh. All numerical calculations are performed in the transformed plane. Since the computational domain is a rectangle and a uniform grid with mesh spacings Δξ = Δη = 1 (in two-dimensions) is used, the computer programming is substantially facilitated. By means of control functions, which form the r.h.s. of the Poisson equations, the clustering of grid lines or grid points is governed. This allows a very fine resolution at certain specified locations and includes adaptive grid generation. The first two sections outline the general features of BFGs, and in section 3 the general transformation rules along with the necessary concepts of differential geometry are given. In section 4 the transformed grid generation equations are derived and control functions are specified. Expressions for grid adaptation arc also presented. Section 5 briefly discusses the numerical solution of the transformed grid generation equations using sucessive overrelaxation and shows a sample calculation where the FAS (full approximation scheme) multigrid technique was employed. In the companion paper (Part II), the application of the BFG method to selected fluid flow problems is addressed.     

Prediction of incompressible flow separation with the approximate factorization technique

This paper describes one application of the approximate factorization technique to the solution of incompressible steady viscous flow problems in two dimensions. The velocity-pressure formulation of the Navier-Stokes equations written in curvilinear non-orthogonal co-ordinates is adopted. The continuity equation is replaced with one equation for the pressure by means of the artificial compressibility concept to obtain a system parabolic in time. The resulting equations are discretized in space with centred finite differences, and the steady state solution obtained by a time-marching ADI method requiring to solve 3 x 3 block tridiagonal linear systems. An optimized fourth-order artificial dissipation is introduced to damp the numerical instabilities of the artificial compressibility equation and ensure convergence. The resulting solver is applied to the prediction of a wide variety of internal flows, including both streamlined boundaries and sharp corners, and fast convergence and good results obtained for all the configurations investigated.     

NUMERICAL SIMULATION OF UNSTEADY VISCOUS FLOWS USING AN IMPLICIT PROJECTION METHOD

In this paper an implicit projection method for the solution of the two-dimensional, time-dependent, incompressible Navier– Stokes equations is presented. The basic principle of this method is that the evaluation of the time evolution is split into intermediate steps. The computational method is based on the approximate factorization technique. The coupled approach is used to link the equations of motion and the turbulence model equations. The standard k-ϵ turbulence model is used. The current methodology, which has been tested extensively for steady problems, is now applied for the numerical simulation of unsteady flows. Several cases were tested, such as plane or axisymmetric channels, a backward-facing step, a square cavity and an axisymmetric stenosis.     

Numerical solution of coupled burgers equations in inhomogeneous form

A finite difference scheme based on the operator-splitting technique with cubic spline functions is derived for solving the two-dimensional Burgers equations in ‘inhomogeneous’ form. The scheme is of first-order accuracy in time and second-order accuracy in space direction and is unconditionally stable. The numerical results are obtained with severe/moderate gradients in the initial and boundary conditions and the steady state solutions are plotted for different values of the parameters. It is concluded that the resulting scheme works very well even in the case of very severe gradient in the solution. Also, the general nature of the scheme provides a wider application in the solution of non-linear problems.     

Computational technique for flow in blood vessels with porous effects

IntroductionThehemodynamicsofflowsthroughbloodvesselsisofgreatinterest,becausethesevesselspresentasubstantialhealthriskandareamajorcauseofmortalityandmorbidityintheindustrializedworld .Researchpapersonthebloodflowhaveappearedbutmostofthemhaveneglectedtheporosityeffectsduetovesselwalls.Inthisstudyweareinterestedintheflowthroughabloodvesseltakingintoaccounttheporosityeffectsofthevessels.Fluidflowthroughaporousmediumisoffundamentalimportancetowiderangeofdisciplinesinthevariousbranchesofnaturalsci…     

流体力学问题的三次样条配置法

本文给出了三次样条配置法在流体力学问题数值解中的应用以及在这一领域的新进展。给出了流体力学方程中主要的样条函数关系和解算步骤。所有情形都是便于反演的三对角形矩阵。简要评述了SADI方法和样条方法在每一坐标方向的分步计算方法、截断误差和稳定性。给出了处理混合边界条件的一般公式。最后简要讨论了样条近似引起的数值弥散和耗散。      

Determining the limiting solutions of nonstationary axisymmetric Hele-Shaw problems

Numerical solution of the Hele-Shaw problem reduces to solution of three boundary-value problems of determining analytic functions of a complex variable in each time step: conformal mapping of the range of the parametric variable to the physical plane, the Dirichlet problems for determining the electric-field strength, and the Riemann-Hilbert problem for calculating partial time derivatives of the coordinates of points of the interelectrode space (the images of the points on the boundary of the parametric plane are fixed). Unlike in the two-dimensional problem, the electric-field strength is determined using integral transformations of an analytic function. Approximation by spline function is performed, and more accurate and steady (than the well-known ones) general solution algorithms for the nonstationary axisymmetric problems are described. Results of a numerical study of the formation of stationary and self-similar configurations are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 87–99, July–August, 2009.     

Shock capturing for steady,supersonic, two-dimensional isentropic flow

A finite difference scheme is presented for the solution of the two-dimensional equations of steady, supersonic, isentropic flow. The scheme incorporates numerical characteristic decomposition, is shock-capturing by design and incorporates space marching as a result of the assumption that the flow is wholly supersonic in at least one space dimension. Results are shown for problems involving oblique hydraulic jumps and reflection from a wall.     

Calculation of laminar recirculating flows using a local non-staggered grid refinement system

A grid-embedding technique for the solution of two-dimensional incompressible flows governed by the Navier-Stokes equations is presented. A finite volume method with collocated primitive variables is employed to ensure conservation at the interfaces of embedding grids as well as global conservation. The discretized equations are solved simultaneously for the whole domain, providing a strong coupling between regions of different refinement. The formulation presented herein is applicable to uniform or non-uniform Cartesian meshes. The method was applied to the solution of two scalar transport equations, to cavity flows driven by body and shear forces and to a sudden plane contraction flow. The numerical predictions are compared with the exact solutions when available and with experimental data. The results show that neither the convergence rate nor the stability of the method is affected by the presence of embedded grids. Embedded grids provide a better distribution of grid nodes over the computational domain and consequently the solution accuracy was improved. The grid-embedding technique proved also that significant savings in computing time could be achieved.     

Meshless analysis of plasticity with application to crack growth problems

The governing equations for classical rate-independent plasticity are formulated in the framework of meshless method. The special J₂ flow theory for three-dimensional, two-dimensional plane strain and plane stress problems are presented. The numerical procedures, including return mapping algorithm, to obtain the solutions of boundary-value problems in computational plasticity are outlined. For meshless analysis the special treatment of the presence of barriers and mirror symmetries is formulated. The crack growth process in elastic–plastic solid under plane strain and plane stress conditions is analyzed. Numerical results are presented and discussed.     

Two‐dimensional implicit time‐dependent calculations on adaptive unstructured meshes with time evolving boundaries

An implicit multigrid‐driven algorithm for two‐dimensional incompressible laminar viscous flows has been coupled with a solution adaptation method and a mesh movement method for boundary movement. Time‐dependent calculations are performed implicitly by regarding each time step as a steady‐state problem in pseudo‐time. The method of artificial compressibility is used to solve the flow equations. The solution mesh adaptation method performs local mesh refinement using an incremental Delaunay algorithm and mesh coarsening by means of edge collapse. Mesh movement is achieved by modeling the computational domain as an elastic solid and solving the equilibrium equations for the stress field. The solution adaptation method has been validated by comparison with experimental results and other computational results for low Reynolds number flow over a shedding circular cylinder. Preliminary validation of the mesh movement method has been demonstrated by a comparison with experimental results of an oscillating airfoil and with computational results for an oscillating cylinder. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite element computations for unsteady fluid and elastic membrane interaction problems

The interaction between the hydrodynamic forces of a flow field and the elastic forces of adjacent deformable boundaries is described by elastohydrodynamics, a coupled fluid–elastic membrane problem. Direct numerical solution of the unsteady, highly non-linear equations requires that the dynamic evolution of both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. This paper describes a numerical algorithm based on the deformable spatial domain space–time (DSD/ST) finite element method for the unsteady motion of an incompressible, viscous fluid with elastic membrane interaction. The unsteady Navier–Stoke and elastic membrane equations are solved separately using an iterative procedure by the GMRES technique with an incomplete lower-upper (ILU) decomposition at every time instant. One-dimensional, two-dimensional and deformable domain model problems are used to demonstrate the capabilities and accuracy of the present algorithm. Both steady state and transient problems are studied. © 1997 John Wiley & Sons, Ltd.     

Variational iteration method for two-dimensional steady slip flow in micro-channels

In this paper, the two-dimensional steady slip flow in microchannels is investigated. Research on micro flow, especially on micro slip flow, is very important for designing and optimizing the micro electromechanical system (MEMS). The Navier-Stokes equations for two-dimensional steady slip flow in microchannels are reduced to a nonlinear third-order differential equation by using similarity solution. The variational iteration method (VIM) is used to solve this nonlinear equation analytically. Comparison of the result obtained by the present method with numerical solution reveals that the accuracy and fast convergence of the new method.     

Heat transfer analysis of a particle-containing channel flow

This paper describes an analytical model of heat transfer in a two-dimensional, steady, nonreacting particle-containing channel flow. An idealized gas flow of specified uniform velocity between insulated parallel plates is assumed and the nonvaporizing particles are conceptualized as contained within an thin sheet injected at the symmetry plane. Two dimensionless parameters that affect the solution are described. These are the effective gas diffusivityK and the dimensionless particle number densityP. The linear, coupled differential equations governing the energy exchange between the gas and liquid phases are solved by means of the Green's function technique. This procedure yields a Volterra integral-series equation as the solution of the gas-phase energy equation. A series solution of this integral equation is obtained by the method of successive substitutions and terms up to second order are calculated.     

Prediction of thermal stratification in a curved duct with 3D body-fitted co-ordinates

This paper concerns a numerical prediction method for buoyancy-influenced flows using three-dimensional non-orthogonal curvilinear co-ordinates. The numerical analysis of the transformed governing equations for thermal hydraulics is based on a Lagrangian method, in which advected physical values are evaluated by local cubic spline interpolations with third-order accuracy in the three-dimensional computational domain. In addition, the buoyancy and diffusion terms are discretized in the Lagrangian scheme so as to have second-order accuracy with respect to time and space. The Neumann boundary conditions, which have been rather difficult for non-orthogonal co-ordinates to deal with, can be implemented by making use of normal vectors on the physical boundary surfaces and cubic spline interpolations. The developed numerical method is applied to the steady isothermal flow in a curved pipe and the unsteady stratified flow in a curved duct. Both of the predicted values are in good agreement with the experimental results and the validity of the prediction method is confirmed.     

Numerical solution of three‐dimensional velocity–vorticity Navier–Stokes equations by finite difference method

This paper describes the finite difference numerical procedure for solving velocity–vorticity form of the Navier–Stokes equations in three dimensions. The velocity Poisson equations are made parabolic using the false‐transient technique and are solved along with the vorticity transport equations. The parabolic velocity Poisson equations are advanced in time using the alternating direction implicit (ADI) procedure and are solved along with the continuity equation for velocities, thus ensuring a divergence‐free velocity field. The vorticity transport equations in conservative form are solved using the second‐order accurate Adams–Bashforth central difference scheme in order to assure divergence‐free vorticity field in three dimensions. The velocity and vorticity Cartesian components are discretized using a central difference scheme on a staggered grid for accuracy reasons. The application of the ADI procedure for the parabolic velocity Poisson equations along with the continuity equation results in diagonally dominant tri‐diagonal matrix equations. Thus the explicit method for the vorticity equations and the tri‐diagonal matrix algorithm for the Poisson equations combine to give a simplified numerical scheme for solving three‐dimensional problems, which otherwise requires enormous computational effort. For three‐dimensional‐driven cavity flow predictions, the present method is found to be efficient and accurate for the Reynolds number range 100?Re?2000. Copyright © 2004 John Wiley & Sons, Ltd.     

On steady,inviscid shock waves at continuously curved,convex surfaces

An accurate and efficient numerical method for steady, two-dimensional Euler equations is applied to study steady shock waves perpendicular to smooth, convex surfaces. The main subject of study is the flow near both ends of the shock wave: the shock-foot and shock-tip flow. A known analytical model of the inviscid shock-foot flow is critically investigated, analytically and numerically. The results obtained agree with those of the existing analytical model. For the inviscid shock-tip flow, two existing analytical solutions are reviewed. Numerical results are presented which agree with one of these two solutions. Good numerical accuracy is achieved through a monotone, second-order accurate, finite-volume discretization. Good computational efficiency is obtained through iterative defect correction iteration and a multigrid acceleration technique which employs local grid refinement.This work was performed as part of a BRITE/EURAM Area 5 project, under Contract No. AERO-0003/1094.     

A high‐order discontinuous Galerkin method for all‐speed flows

In this article, we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of steady solutions of the compressible fully coupled Reynolds‐averaged Navier–Stokes and k ? ω turbulence model equations for solving all‐speed flows. The system of equations is iterated to steady state by means of an implicit scheme. The DG solution is extended to the incompressible limit by implementing a low Mach number preconditioning technique. A full preconditioning approach is adopted, which modifies both the unsteady terms of the governing equations and the dissipative term of the numerical flux function by means of a new preconditioner, on the basis of a modified version of Turkel's preconditioning matrix. At sonic speed the preconditioner reduces to the identity matrix thus recovering the non‐preconditioned DG discretization. An artificial viscosity term is added to the DG discretized equations to stabilize the solution in the presence of shocks when piecewise approximations of order of accuracy higher than 1 are used. Moreover, several rescaling techniques are implemented in order to overcome ill‐conditioning problems that, in addition to the low Mach number stiffness, can limit the performance of the flow solver. These approaches, through a proper manipulation of the governing equations, reduce unbalances between residuals as a result of the dependence on the size of elements in the computational mesh and because of the inherent differences between turbulent and mean‐flow variables, influencing both the evolution of the Courant Friedrichs Lewy (CFL) number and the inexact solution of the linear systems. The performance of the method is demonstrated by solving three turbulent aerodynamic test cases: the flat plate, the L1T2 high‐lift configuration and the RAE2822 airfoil (Case 9). The computations are performed at different Mach numbers using various degrees of polynomial approximations to analyze the influence of the proposed numerical strategies on the accuracy, efficiency and robustness of a high‐order DG solver at different flow regimes. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical simulation of vortical ideal fluid flow through curved channel

A numerical algorithm to study the boundary‐value problem in which the governing equations are the steady Euler equations and the vorticity is given on the inflow parts of the domain boundary is developed. The Euler equations are implemented in terms of the stream function and vorticity. An irregular physical domain is transformed into a rectangle in the computational domain and the Euler equations are rewritten with respect to a curvilinear co‐ordinate system. The convergence of the finite‐difference equations to the exact solution is shown experimentally for the test problems by comparing the computational results with the exact solutions on the sequence of grids. To find the pressure from the known vorticity and stream function, the Euler equations are utilized in the Gromeka–Lamb form. The numerical algorithm is illustrated with several examples of steady flow through a two‐dimensional channel with curved walls. The analysis of calculations shows strong dependence of the pressure field on the vorticity given at the inflow parts of the boundary. Plots of the flow structure and isobars, for different geometries of channel and for different values of vorticity on entrance, are also presented. Copyright © 2003 John Wiley & Sons, Ltd.     

Analysis of high-order-accurate fluctuation splitting schemes for steady advection

This article provides an analysis of high-order-accurate two-dimensional fluctuation splitting schemes for steady advection. Using Lagrangian elements, a residual distribution scheme is formulated and its properties are assessed. Distributing the residuals over the sub-triangles of each element allows one to obtain a well-posed scheme. It is also shown that the standard elemental approach is ill-posed, insofar as it produces an undetermined linear system. A steady scalar-advection problem is used to verify that the numerical schemes do obey the analysis.