一种新型的高分辨率通量差分裂格式

传统的Roe格式不满足熵条件并且在计算激波问题时会遭遇不同形式的不稳定现象,如慢行激波的波后振荡和红玉(carbuncle)现象.基于Zha-Bilgen对流-压力通量分裂方法,构造一种新型的通量差分裂格式.利用约旦标准型理论,通过添加广义特征向量构造通量差分裂方法来计算对流子系统.压力子系统具有一组完备的线性无关特征向量,因此可以构造传统的通量差分裂格式进行计算.为了提高接触间断的分辨率,利用界面变差下降(BVD)算法来重构对流通量耗散项中的密度差.激波稳定性分析表明,新格式可以有效地衰减数值误差,从而抑制不稳定现象的发生.一系列数值实验证明了本文构造的新型通量差分裂格式比Roe格式具有更高的分辨率和更好的鲁棒性.     

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.     

On the flow in underexpanded supersonic water jets

The problem of a supersonic water jet exhausting from a cylindrical nozzle in the underexpansion regime is solved. The calculations are made with allowance lor splitting. The distributions of the parameters over the jet are given. It is shown that the splitting region depends on the degree of underexpansion of the regime. The solution is found by the method of stabilization based on the difference scheme of Godunov, Zabrodin, and Prokopov [3] generalized to the case of flow with splitting.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 155–157, March–April, 1982.I thank G. A. Atanov for discussing the work.     

Development of an artificial compressibility methodology using flux vector splitting

An implicit, upwind arithmetic scheme that is efficient for the solution of laminar, steady, incompressible, two-dimensional flow fields in a generalised co-ordinate system is presented in this paper. The developed algorithm is based on the extended flux-vector-splitting (FVS) method for solving incompressible flow fields. As in the case of compressible flows, the FVS method consists of the decomposition of the convective fluxes into positive and negative parts that transmit information from the upstream and downstream flow field respectively. The extension of this method to the solution of incompressible flows is achieved by the method of artificial compressibility, whereby an artificial time derivative of the pressure is added to the continuity equation. In this way the incompressible equations take on a hyperbolic character with pseudopressure waves propagating with finite speed. In such problems the ‘information’ inside the field is transmitted along its characteristic curves. In this sense, we can use upwind schemes to represent the finite volume scheme of the problem's governing equations. For the representation of the problem variables at the cell faces, upwind schemes up to third order of accuracy are used, while for the development of a time-iterative procedure a first-order-accurate Euler backward-time difference scheme is used and a second-order central differencing for the shear stresses is presented. The discretized Navier–Stokes equations are solved by an implicit unfactored method using Newton iterations and Gauss–Siedel relaxation. To validate the derived arithmetical results against experimental data and other numerical solutions, various laminar flows with known behaviour from the literature are examined. © 1997 John Wiley & Sons, Ltd.     

Composite high resolution localized relaxation scheme based on upwinding for hyperbolic conservation laws

In this work we present an upwind‐based high resolution scheme using flux limiters. Based on the direction of flow we choose the smoothness parameter in such a way that it leads to a truly upwind scheme without losing total variation diminishing (TVD) property for hyperbolic linear systems where characteristic values can be of either sign. Here we present and justify the choice of smoothness parameters. The numerical flux function of a high resolution scheme is constructed using wave speed splitting so that it results into a scheme that truly respects the physical hyperbolicity property. Bounds are given for limiter functions to satisfy TVD property. The proposed scheme is extended for non‐linear problems by using the framework of relaxation system that converts a non‐linear conservation law into a system of linear convection equations with a non‐linear source term. The characteristic speed of relaxation system is chosen locally on three point stencil of grid. This obtained relaxation system is solved using composite scheme technique, i.e. using a combination of proposed scheme with the conservative non‐standard finite difference scheme. Presented numerical results show higher resolution near discontinuity without introducing spurious oscillations. Copyright © 2008 John Wiley & Sons, Ltd.     

Transonic flow of a gas jet from a vessel with plane walls

The problem of transonic flow of a gas jet with equilibrium excitation of the vibrational degrees of freedom of the molecules from an infinite symmetrical vessel with plane walls is reduced to a generalized Tricomi boundary-value problem for an equation of Chaplygin type. It is solved using a difference scheme [13] based on the decomposition of the difference operator in accordance with the type of differential operator. Calculation results are presented for a mixture of oxygen and nitrogen simulating air. The effect of the angle between the walls, the stagnation temperature and the back-pressure coefficient on the flow coefficient is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 128–135, November–December, 1987.The authors are grateful to F. Yu. Stepanov for his useful comments.     

Numerical solution of 2D and 3D turbulent internal flow problems

The paper describes a method for solving numerically two-dimensional or axisymmetric, and three-dimensional turbulent internal flow problems. The method is based on an implicit upwinding relaxation scheme with an arbitrarily shaped conservative control volume. The compressible Reynolds-averaged Navier-Stokes equations are solved with a two-equation turbulence model. All these equations are expressed by using a non-orthogonal curvilinear co-ordinate system. The method is applied to study the compressible internal flow in modern power installations. It has been observed that predictions for two-dimensional and three-dimensional channels show very good agreement with experimental results.     

Multigrid flow solutions in complex two-dimensional geometries

A finite difference solution algorithm is described for use on two-dimensional curvilinear meshes generated by the solution of the transformed Laplace equation. The efficiency of the algorithm is improved through the use of a full approximation scheme (FAS) multigrid algorithm using an extended pressure correction scheme as smoother. The multigrid algorithm is implemented as a fixed V-cycle through the grid levels with a constant number of sweeps being performed at each grid level. The accuracy and efficiency of the numerical code are validated using comparisons of the flow over two backward step configurations. Results show close agreement with previous numerical predictions and experimental data. Using a standard Cartesian co-ordinate flow solver, the multigrid efficiency obtainable in a rectangular system is shown to be reproducible in two-dimensional body-fitted curvilinear co-ordinates. Comparisons with a standard one-grid method show the multigrid method, on curvilinear meshes, to give reductions in CPU time of up to 93%.     

Flux difference splitting for the Euler equations in one spatial co-ordinate with area variation

An approximate (linearized) Riemann solver is presented for the solution of the Euler equations of gas dynamics in one spatial co-ordinate. This includes cylindrically and spherically symmetric geometries and also applies to narrow ducts with area variation. The method is Roe's flux difference splitting with a technique for dealing with source terms. The results of two problems, a spherically divergent infinite shock and a converging cylindrical shock, are presented. The numerical results compare favourably with those of Noh's recent survey and also with those of Ben-Artzi and Falcovitz using a more complicated Riemann solver.     

Development of the generalized MacCormack scheme and its extension to low Mach number flows

The low Mach number performance of the MacCormack scheme is examined. The inherent dissipation in the scheme is found to suffer from the degradation in accuracy observed with traditional, density‐based methods for compressible flows. Two specific modifications are proposed, leading to the formation of the generalized MacCormack scheme within a dual‐time framework (called GMC‐PC). The first modification involves reformulating the flux by splitting it into particle convection and acoustic parts, with the former terms treated using the traditional MacCormack discretization and the latter terms augmented by the addition of a pressure‐based artificial dissipation. The second modification involves a reformulation of the traditional nonlinear fix introduced by MacCormack in 1971, which is found to be necessary to suppress pressure oscillations at low Mach numbers. The new scheme is demonstrated to have superior performance, independent of Mach number, compared with standard MacCormack implementations using several canonical test problems. Copyright © 2017 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid

This paper presents a finite element method for incompressible multiphase flows with capillary interfaces based on a (formally) second‐order projection scheme. The discretization is on a fixed Eulerian grid. The fluid phases are identified and advected using a level set function. The grid is temporarily adapted around the interfaces in order to maintain optimal interpolations accounting for the pressure jump and the discontinuity of the normal velocity derivatives. The least‐squares method for computing the curvature is used, combined with piecewise linear approximation to the interface. The time integration is based on a formally second order splitting scheme. The convection substep is integrated over an Eulerian grid using an explicit scheme. The remaining generalized Stokes problem is solved by means of a formally second order pressure‐stabilized projection scheme. The pressure boundary condition on the free interface is imposed in a strong form (pointwise) at the pressure‐computation substep. This allows capturing significant pressure jumps across the interface without creating spurious instabilities. This method is simple and efficient, as demonstrated by the numerical experiments on a wide range of free‐surface problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Analyzing dynamic response of non-homogeneous string fixed at both ends

To survey the local conservation properties of complex dynamic problems, a structure-preserving numerical method, named as generalized multi-symplectic method, is proposed to analyze the dynamic response of a non-homogeneous string fixed at both ends. Firstly, based on the multi-symplectic idea, a generalized multi-symplectic form derived from the vibration equation of a non-homogeneous string is presented. Secondly, several conservation laws are deduced from the generalized multi-symplectic form to illustrate the local properties of the system. Thirdly, a centered box difference scheme satisfying the discrete local momentum conservation law exactly, named as a generalized multi-symplectic scheme, is constructed to analyze the dynamic response of the non-homogeneous string fixed at both ends. Finally, numerical experiments on the generalized multi-symplectic scheme are reported. The results illustrate the high accuracy, the good local conservation properties as well as the excellent long-time numerical behavior of the generalized multi-symplectic scheme well.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid: 3D validation

A three‐dimensional finite element method for incompressible multiphase flows with capillary interfaces is developed based on a (formally) second‐order projection scheme. The discretization is on a fixed (Eulerian) reference grid with an edge‐based local h‐refinement in the neighbourhood of the interfaces. The fluid phases are identified and advected using the level‐set function. The reference grid is then temporarily reconnected around the interface to maintain optimal interpolations accounting for the singularities of the primary variables. Using a time splitting procedure, the convection substep is integrated with an explicit scheme. The remaining generalized Stokes problem is solved by means of a pressure‐stabilized projection. This method is simple and efficient, as demonstrated by a wide range of difficult free‐surface validation problems, considered in the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

High‐resolution,monotone solution of the adjoint shallow‐water equations

A monotone, second‐order accurate numerical scheme is presented for solving the differential form of the adjoint shallow‐water equations in generalized two‐dimensional coordinates. Fluctuation‐splitting is utilized to achieve a high‐resolution solution of the equations in primitive form. One‐step and two‐step schemes are presented and shown to achieve solutions of similarly high accuracy in one dimension. However, the two‐step method is shown to yield more accurate solutions to problems in which unsteady wave speeds are present. In two dimensions, the two‐step scheme is tested in the context of two parameter identification problems, and it is shown to accurately transmit the information needed to identify unknown forcing parameters based on measurements of the system response. The first problem involves the identification of an upstream flood hydrograph based on downstream depth measurements. The second problem involves the identification of a long wave state in the far‐field based on near‐field depth measurements. Copyright © 2002 John Wiley & Sons, Ltd.     

On the numerical solution of the turbulence energy equations for wave and tidal flows

This paper deals with the numerical solution, using finite difference methods, of the hydrodynamic and turbulence energy equations which describe wind wave and tidally induced flow. Calculations are performed using staggered and non-staggered finite difference grids in the vertical, with various time discretizations of the production and dissipation terms in the turbulence energy equations. It is shown that the time discretization of these terms can significantly influence the stability of the solution. The effect of time filtering on the numerical stability of the solution is also considered. The form of the mixing length is shown to significantly influence the bed stress in wind wave problems. A no-slip condition is applied at the sea bed, and the associated high-shear bottom boundary layer is resolved by transforming the equations onto a logarithmic or log-linear co-ordinate system before applying the finite difference scheme. A computationally economic method is developed which remains stable even when a very fine vertical grid (over 200 points) is used with a time step of up to 30 min.     

A numerical method for 3D viscous incompressible flows using non-orthogonal grids

This paper presents a numerical method for fluid flow in complex three-dimensional geometries using a body-fitted co-ordinate system. A new second-order-accurate scheme for the cross-derivative terms is proposed to describe the non-orthogonal components, allowing parts of these terms to be treated implicitly without increasing the number of computational molecules. The physical tangential velocity components resulting from the velocity expansion in the unit tangent vector basis are used as dependent variables in the momentum equations. A coupled equation solver is used in place of the complicated pressure correction equation associated with grid non-orthogonality. The co-ordinate-invariant conservation equations and the physical geometric quantities of control cells are used directly to formulate the numerical scheme, without reference to the co-ordinate derivatives of transformation. Several two- and three-dimensional laminar flows are computed and compared with other numerical, experimental and analytical results to validate the solution method. Good agreement is obtained in all cases.     

高精度迎风紧致差分格式与热流计算研究的注记

利用高精度差分格式求解了可压缩 N-S方程球头热流问题。分析了不同差分格式在对球头粘性绕流热流计算中存在的问题 ,并分析了相应的网格雷诺数。在利用高精度迎风紧致 [1 ] 格式求解粘性绕流热流问题时 ,采用 Steger-Warming[2 ]的通量分裂技术将守恒型方程中的流通向量分裂成两部分 ,在此基础上据风向构造逼近于无粘项的高精度迎风格式。对方程中的粘性部分采用中心差分格式。数值结果表明 :高精度差分格式能在较大的网格雷诺数下较好地计算球头驻点热流     

Solution to Euler equations by high‐resolution upwind compact scheme based on flux splitting

A high‐resolution upwind compact method based on flux splitting is developed for solving the compressive Euler equations. The convective flux terms are discretized by using the modified advection upstream splitting method (AUSM). The developed scheme is used to compute the one‐dimensional Burgers equation and four different example problems of supersonic compressible flows, respectively. The results show that the high‐resolution upwind compact scheme based on modified AUSM+ flux splitting can capture shock wave and other discontinuities, obtain higher resolution and restrain numerical oscillation. Copyright © 2007 John Wiley & Sons, Ltd.     

Two-dimensional interactions due to moving load in generalized thermoelastic solid with diffusion

The present paper is concerned with the investigation of disturbances in＇a homogeneous, isotropic elastic medium with generalized thermoelastic diffusion, when a moving source is acting along one of the co-ordinate axis on the boundary of the medium. Eigen value approach is applied to study the disturbance in Laplace-Fourier transform domain for a two dimensional problem. The analytical expressions for displacement components, stresses, temperature field, concentration and chemical potential are obtained in the physical domain by using a numerical technique for the inversion of Laplace transform based on Fourier expansion techniques. These expressions are calculated numerically for a copper like material and depicted graphically. As special cases, the results in generalized thermoelastic and elastic media are obtained. Effect of presence of diffusion is analyzed theoretically and numerically.     

How to define the storage and loss moduli for a rheologically nonlinear material?

A large amplitude oscillatory shear (LAOS) is considered in the strain-controlled regime, and the interrelation between the Fourier transform and the stress decomposition approaches is established. Several definitions of the generalized storage and loss moduli are examined in a unified conceptual scheme based on the Lissajous–Bowditch plots. An illustrative example of evaluating the generalized moduli from a LAOS flow is given.