A novel slightly compressible model for low Mach number perfect gas flow calculation

By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved very efficient, for it is kept within thep-v frame but does not have to satisfy the time consuming divergence-free condition in order to get the incompressible Navier-Stokes equation solution. Writing the equations in the form of conservation laws, we have derived the characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used for the equation solutions and a new Exact Newton Relaxation (ENR) implicit method is developed. Various computed results are presented to validate the present model. Laminar flow solutions over a circular cylinder with wake developing and vortex shedding are presented. Results for inviscid flow over a sphere are compared in excellent agreement with the exact analytic incompressible solution. Three-dimensional viscous flow solutions over sphere and prolate spheroid are also calculated and compared well with experiments and other incompressible solutions. Finally, good convergent performances are shown for sphere viscous flows. The project supported by the Basic Research on Frontier Problems in Fluid and Aerodynamics in China and the National Natural Science Foundation of China (19772069)     

An efficient shock-capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.     

A spatial laminar boundary layer on a streamline in conical external flow of a uniform gas in the absence of sources and sinks

Theoretical study of a three-dimensional laminar boundary layer is a complex problem, but it can be substantially simplified in certain particular cases and even reduced to the solution of ordinary differential equations.One such particular case is the flow of a compressible gas on a streamline in conical external flow. The case is of considerable practical importance because the local heat fluxes may take extremal values on such lines.Such flow, except for the conical case, has been examined [1–4], and an approximate method has been given [1] on the basis of integral relationships and a special form for the approximating functions. A numerical solution has been given [2, 3] for such flow around an infinite cylinder. It was assumed in [1–3] that the Prandtl number and the specific heats were constant, and that the dynamic viscosity was proportional to temperature. Heat transfer has been examined [4] near a cylinder exposed to a flow of dissociated air.Here we give results from numerical solution of a system of ordinary differential equations for the flow of a compressible gas in a laminar boundary layer on streamlines in conical external flow, with or without influx or withdrawal of a homogeneous gas. It is assumed that the gas is perfect and that the dynamic viscosity has a power-law temperature dependence.     

Extension of a high‐resolution scheme to 1D liquid–gas flow

This paper presents the extension of a high‐resolution conservative scheme to the one‐dimensional one‐pressure six‐equation two‐fluid flow model. Only mixtures of water and air have been considered in this study, both fluids have been characterized using simple equations of state, namely stiffened gas for the liquid phase and perfect gas for the gas phase. The resulting scheme is explicit and first‐order accurate in space and time. A second‐order version of the scheme has also been derived using the MUSCL strategy and slope limiters. Some numerical results show the good capabilities of this type of schemes in the solution of discontinuities in two‐fluid flow problems, all of them are based on water/air numerical benchmarks widely used in the two‐phase flow literature. Copyright © 2005 John Wiley & Sons, Ltd.     

A shock-capturing scheme for body-fitted meshes

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional Euler equations of gas dynamics in a generalized co-ordinate system. The scheme is based on numerical characteristic decomposition and solves locally linearized Riemann problems using upwind differencing. The decomposition is for a generalized co-ordinate system and a convex equation of state. This ensures good shock-capturing properties when incorporated with operator splitting and the advantage of using body-fitted co-ordinates. The resulting scheme is applied to supersonic flow of real air' past a circular cylinder.     

Numerical investigation of hypersonic step-flows

Hypersonic aerospace vehicles are exposed to extreme flight conditions with heavy contour loads during their mission. Especially at ridges and sharp corners, the wall heat flux and pressure may cause serious damage to the body. Sometimes, the surface material cannot resist the high loading and fails completely. In this work the laminar hypersonic flow over forward and backward facing steps is investigated by CFD techniques and the results are compared with experimental data. The selected flow conditions correspond to cold hypersonic flow according to the availability of experimental data. The Navier-Stokes equations in the high temperature gas approximation of a thermally perfect gas in local equilibrium serve as the model for the physical problem. A multiblock finite-volume method is used to discretize consistently all spatial derivatives appearing in the balance equations. A second order in space Godunov-type method is utilized for the non-diffusive part of the governing equations whereas centered differences are used for the diffusive part. Time integration is performed by a second order implicit scheme. In each time step, the resulting nonlinear system of equations is solved by Newton's method employing a relaxation scheme based on conjugate gradients for the linear equation system. The results obtained permit a close insight into the physics of the flow problems under consideration and by this provide valuable information for construction concepts of hypersonic vehicles. Besides a careful comparison of the numerical results with experimental data, numerical aspects like the grid influence are addressed. Received 9 November 1998 / Accepted 2 December 1999     

The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

非定常不可压N-S方程的最小二乘算子分裂有限元法数值求解

采用最小二乘算子分裂有限元法求解非定常不可压N-S(Navier-Stokes)方程,即在每个时间层上采用算子分裂法将N-S方程分裂成扩散项和对流项,这样既能考虑对流占优特点又能顾及方程的扩散性质。扩散项是一个抛物型方程,时间离散采用向后差分格式,空间离散采用标准Galerkin有限元法。对流项的时间项采用后向差分格式,非线性部分用牛顿法进行线性化处理,再用最小二乘有限元法进行空间离散,得到对称正定的代数方程组系数矩阵。采用Re=1000的方腔流对该算法的有效性进行检验,表明其具有较高的精度,能够很好地捕捉流场中的涡结构。同时,对圆柱层流绕流进行了数值研究,通过流线图、压力场、阻力系数、升力系数及斯特劳哈数等结果的分析与对比,表明本文算法对于模拟圆柱层流绕流是准确和可靠的。     

Finite analytic boundary condition for a higher‐order finite analytic scheme

A higher‐order finite analytic scheme based on one‐dimensional finite analytic solutions is used to discretize three‐dimensional equations governing turbulent incompressible free surface flow. In order to preserve the accuracy of the numerical scheme, a new, finite analytic boundary condition is proposed for an accurate numerical solution of the partial differential equation. This condition has higher‐order accuracy. Thus, the same order of accuracy is used for the boundary. Boundary conditions were formulated and derived for fluid inflow, outflow, impermeable surfaces and symmetry planes. The derived boundary conditions are treated implicitly and updated with the solution of the problem. The basic idea for the derivation of boundary conditions was to use the discretized form of the governing equations for the fluid flow simplified on the boundaries and flow information. To illustrate the influence of the higher‐order effects at the boundaries, another, lower‐order finite analytic boundary condition, is suggested. The simulations are performed to demonstrate the validity of the present scheme and boundary conditions for a Wigley hull advancing in calm water. Copyright © 2005 John Wiley & Sons, Ltd.     

绕Apollo飞船的高超声速化学非平衡流动的数值模拟

利用混合通量分裂方法，建立了很方便求解的隐式ＮＮＤ格式，求解了完全气体和化学非平衡空气绕Ａｐｏｌｌｏ飞船的流动，计算结果和实验值作了比较，应用拓扑分析方法，研究了背风区和尾迹内的流动结构。     

Flux difference splitting for open-channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Supersonic low-concentration dusty-gas flow past a plane cylinder in the presence of an oblique shock wave interacting with the bow shock

The motion of an inertial dispersed admixture near a plane cylinder immersed in a steady-state hypersonic dusty flow in the presence of an oblique shock wave interacting with the bow shock is considered. It is assumed that the free-stream particle mass concentration is small and the particles do not affect the carrier flow. The III and IV shock wave interaction regimes are considered. The gas flow parameters in the shock layer are calculated from the numerical solution of the full Navier-Stokes equations for the perfect gas. A TVD second-order finite-difference scheme constructed on the basis of a finite volume method is used. For calculating the dispersed-phase parameters, including the concentration, the full Lagrangian method is used. On a wide range of variation of the particle inertia parameters, the patterns of the particle trajectories, velocity, concentration, and temperature in the shock layer are studied. The possibility of aerodynamic focusing of the particles behind the shock wave intersection point and the formation of narrow beams with a high particle concentration is revealed. These beams impinge on the cylinder surface and result in a sharp increase in the local heat fluxes. The maximal possible increase in the heat fluxes caused by the particles colliding with the cylinder surface is estimated for the flows with and without the incident oblique shock wave.     

Simulation of supersonic flows on the basis of splitting algorithms

For the numerical simulation of aerodynamics problems, the Euler and Navier — Stokes equations written in integral form are used to construct an implicit finite-volume predictor-corrector scheme. At the predictor stage, the splitting of equations into physical processes and spatial directions is introduced, which makes it possible to reduce the solution of the original system to the solution of individual equations on fractional steps by the scalar sweep method and ensure the stability of the algorithm as a whole. The paper also describes the supersonic gas flows in a narrowing channel with regular and non-regular reflection of the compression shock from the symmetry plane and the numerical substantiation of the existence of pulsating flow with a supersonic flow past a cylinder with a needle.     

A numerical study of projection algorithms in the finite element simulation of three‐dimensional viscous incompressible flow

In this work a comparative study of two versions of the projection algorithm used either for time integration or as an iterative method to solve the three‐dimensional incompressible Navier–Stokes equations is presented. It is also shown that these projection algorithms combined with the finite element method are particularly suited for the treatment of outflow boundary conditions in the context of external flows. This assertion is illustrated by means of some numerical examples in which five types of boundary conditions are compared. The scheme is applied to simulate the flow past a cylinder clamped on two fixed parallel solid walls. Comparison with experimental data available for this problem shows good agreement of the velocity and pressure fields, both computed with continuous piecewise linear elements. Copyright © 1999 John Wiley & Sons, Ltd.     

Cylinder in viscous incompressible fluid flow

We present a technique for calculating the temperature field in the vicinity of a cylinder in a viscous incompressible fluid flow under given conditions for the heat flux or the cylinder surface temperature. The Navier-Stokes equations and the energy equation for the steady heat transfer regime form the basis of the calculations. The numerical calculations are made for three flow regimes about the cylinder, corresponding to Reynolds numbers of 20, 40, and 80. The pressure distribution, voracity, and temperature distributions along the cylinder surface are found.It is known that for a Reynolds number R>1 the calculation of cylinder drag within the framework of the solution of the Oseen and Stokes equations yields a significant deviation from the experimental data. In 1933 Thom first solved this problem [1] on the basis of the Navier-Stokes equations. Subsequently several investigators [2, 3] studied the problem of viscous incompressible fluid flow past a cylinder.It has been established that a stable solution of the Navier-Stokes equations exists for R40 and that in this case the calculation results are in good agreement with the experimental data. According to [2], a stable solution also exists for R=44. The possibility of obtaining a steady solution for R>44 is suggested.Analysis of the results of [2] permits suggesting that the questions of constructing a difference scheme with a given order of approximation of the basic differential relations which will permit obtaining the sought solution over the entire range of variation of the problem parameters of interest are still worthy of attention.Calculation of the velocity field in the vicinity of a cylinder also makes possible the calculation of the cylinder temperature regime for given conditions for the heat flux or the temperature on its surface. However, we are familiar only with experience in the analytic solution of several questions of cylinder heat transfer with the surrounding fluid for large R within the framework of boundary layer theory [4].     

Calculation of supersonic viscid flow near stagnation line of a blunt body

A numerical study is made of supersonic flow of a viscous gas in the vicinity of the stagnation line of plane and axisymmetric blunt bodies (cylinder, sphere). As in [1–5], which consider the compressed layer of a viscous gas in the vicinity of the stagnation point, use is made of the locally self-similar approximation, which is used to transform the Navier-Stokes equations into a system of ordinary differential equations. In the present paper the solution is sought with the simplifications of [5] and with more general conditions, which makes it possible to study a broad class of flows. The proposed numerical algorithm permits obtaining the structure of the compressed layer near the stagnation line, including the shock wave and the boundary layer. The calculations made on a computer for different flow conditions are illustrated by graphs.The author wishes to thank G. I. Petrov, G. F. Telenin, and L. A., Chudov for their interest in the study and for their helpful discussions. discussions.     

Detached flows past a cylinder with a flat end

In the framework of the Navier-Stokes equations, a numerical solution is found to the problem of longitudinal (axisymmetric) flow of a viscous perfect heat conducting gas past a cylinder with a flat front end at subsonic and transonic velocities of the oncoming flow. The flow in the neighborhood of the corner of the cylinder is investigated and a study made of the occurrence and development of detached flow along the side of the cylinder. Tie heat fluxes to the surface of the cylinder are determined, and the influence of the temperature of the wall of the cylinder on the development of the separation is found. The investigations yield the critical Reynolds number Re₀ for the occurrence of a detached flow as a function of the Mach number of the oncoming flow for the investigated bodies.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 3–10, March–April, 1979.I am grateful to Yu. A. Dem'yanov for discussing the work and for a number of valuable comments.     

Modelling of two‐dimensional laminar flow using finite element method

In this paper, a Galerkin weighted residual finite element numerical solution method, with velocity material time derivative discretisation, is applied to solve for a classical fluid mechanics system of partial differential equations modelling two‐dimensional stationary incompressible Newtonian fluid flow. Classical examples of driven cavity laminar flow and laminar flow past a cylinder are presented. Numerical results are compared with data found in the literature. Copyright © 1999 John Wiley & Sons, Ltd.     

THE APPLICATION OF A FINITE DIFFERENCE/FINITE VOLUME ALGORITHM TO SOLVE MAXWELL'S EQUATIONS IN THE TIME DOMAIN

Abstract

A finite volume/finite difference method based on Ni's multigrid formulation is introduced for the solution of Maxwell's equations. The scheme is presented for the cases of transverse magnetic scattering from two-dimensional circular and square cylinders, as well as from NACA 0012 airfoil. The codes are validated against the traditional Method of Moments, which is analogous to a panel method in CFD. The circular cylinder scattering is compared to the analytical series solution for better understanding how the roles of numerical dispersion and dissipation errors affect the solution. The reflecting boundary conditions are modeled by the idea of inducing fields inside the conductor and a method of modeling the singularities that arise at a sharp corner is presented. Absorbing boundary conditions are modeled by integrating along the characteristic compatibility equations in the direction of the outgoing wave.