解弹性力学第二类边界积分方程的求积法与分裂外推

利用Sidi奇异求积公式,提出了解曲边多角形域上线性弹性力学第二类边界积分方程的求积法,即离散矩阵的每个元素的生成只需赋值不需计算任何奇异积分.通过估计离散矩阵的特征值和利用Anselone聚紧收敛理论,证明了近似解的收敛性;同时得到了误差的多参数渐近展开式;通过并行地解粗网格上的离散方程,利用分裂外推获得了高精度近似解和后验误差.     

稳态问题混合边界积分方程的高精度求积法与分裂外推

提出了求积法解稳态问题的混合边界积分方程,它拥有高精度,低复杂度.通过并行地解粗网格上的离散方程,根据误差的多参数渐近展开,应用分裂外推算法得到高精度的近似解,同时获得后验误差估计.     

Comparison of Several Spatial Discretizations for the Navier–Stokes Equations

Grid convergence studies for subsonic and transonic flows over airfoils are presented in order to compare the accuracy of several spatial discretizations for the compressible Navier–Stokes equations. The discretizations include the following schemes for the inviscid fluxes: (1) second-order-accurate centered differences with third-order matrix numerical dissipation, (2) the second-order convective upstream split pressure scheme (CUSP), (3) third-order upwind-biased differencing with Roe's flux-difference splitting, and (4) fourth-order centered differences with third-order matrix numerical dissipation. The first three are combined with second-order differencing for the grid metrics and viscous terms. The fourth discretization uses fourth-order differencing for the grid metrics and viscous terms, as well as higher-order approximations near boundaries and for the numerical integration used to calculate forces and moments. The results indicate that the discretization using higher-order approximations for all terms is substantially more accurate than the others, producing less than two percent numerical error in lift and drag components on grids with less than 13,000 nodes for subsonic cases and less than 18,000 nodes for transonic cases. Since the cost per grid node of all of the discretizations studied is comparable, the higher-order discretization produces solutions of a given accuracy much more efficiently than the others.     

Difference equation solutions for non-steady-state diffusion of charged particles in solids

The hopping transport of charged particles through a solid is described by means of difference equations based on the concept of classical thermal motion over discrete energy barriers. Homogeneous electric fields and concentration gradients are considered to be the driving forces for transport. Transient and steady-state currents are derived, and the concentration profiles are obtained for the mobile charged defect species. For the case of a slab geometry the discreteness of the potential barriers leads to a nonlinear dependence of current on voltage in the high electric field limit, with a more rapid increase of current with voltage than would be expected from an extrapolation of the low field linear dependence. The field-dependent relaxation of a non-steady-state defect concentration profile to the corresponding steady-state profile can be nearly exponential in the limit of large fields. Tracer distributions for the cases of semi-infinite and unbounded diffusion mediums are likewise affected appreciably in the high field limit. The velocity of the peak is increased over that obtained by linear extrapolation from the low field limit. It is concluded that the combination of high electric fields with the natural microscopic discreteness of a solid-state diffusion medium can result in readily observable nonlinear electric field effects which increase approximately exponentially with the atomic separation distances of the discrete barriers in hopping transport. Some of this nonlinear behavior can be retained in differential equations derived from the difference equations by means of Taylor series expansions of the carrier concentration with respect to position.     

Development of arbitrary-order multioperators-based schemes for parallel calculations.: Part 2: Families of compact approximations with two-diagonal inversions and related multioperators

One-parameter families of compact approximations to grid functionals with inverses of two-point operators and their properties are described. As particular examples, interpolation/extrapolations operators, quadratures formulas and approximations to derivatives are presented. Using operators from the families with fixed parameters values as basis operators, their linear combinations providing formally arbitrary-order approximations (multioperators) are constructed. Numerical illustrations are presented. Special emphasis is placed on first derivatives discretizations in the context of conservation laws. As an example, a highly accurate tenth-order scheme is outlined and tested against the Burgers’ equation. It is shown how extrapolation multioperators can be used to create boundary closures.     

高密度和高粘度比率下气液两相流动的数值模拟

本文在VOF方法的基础上,采用粗细两套网格对高密度和高粘度比率下的气液两相流动模拟进行了研究分析.在细网格中求解流体体积函数方程,在粗网格中采用交错网格求解动量方程和压力修正方程,通过粗细网格间的数据传递获得求解动量方程时需要的准确的界面密度和粘度及控制体密度,克服了高密度和高粘度比率下通过插值方法计算界面密度和粘度及控制体密度带来较大误差的困难,保证了质量和动量同时守恒.高密度和高粘度比率下气液两相流动中气液交界面处密度、速度和压力急剧变化,为了保证格式的有界性和稳定性,采用稳定的有界高阶组合格式STOIC.最后模拟了不同工况下气泡在液体中的运动,并通过实验和模拟结果验证了方法的可行性及准确性.     

Reply to comment on “On the controversy around Daganzo’s requiem for and Aw-Rascle’s resurrection of second-order traffic flow models” by H.M. Zhang

This contribution compares several different approaches allowing one to derive macroscopic traffic equation directly from microscopic car-following models. While it is shown that some conventional approaches lead to theoretical problems, it is proposed to use an approach reminding of smoothed particle hydrodynamics to avoid gradient expansions. The derivation circumvents approximations and, therefore, demonstrates the large range of validity of macroscopic traffic equations, without the need of averaging over many vehicles. It also gives an expression for the “traffic pressure”, which generalizes previously used formulas. Furthermore, the method avoids theoretical inconsistencies of macroscopic traffic models, which have been criticized in the past by Daganzo and others.     

Derivation of non-local macroscopic traffic equations and consistent traffic pressures from microscopic car-following models

This contribution compares several different approaches allowing one to derive macroscopic traffic equation directly from microscopic car-following models. While it is shown that some conventional approaches lead to theoretical problems, it is proposed to use an approach reminding of smoothed particle hydrodynamics to avoid gradient expansions. The derivation circumvents approximations and, therefore, demonstrates the large range of validity of macroscopic traffic equations, without the need of averaging over many vehicles. It also gives an expression for the “traffic pressure”, which generalizes previously used formulas. Furthermore, the method avoids theoretical inconsistencies of macroscopic traffic models, which have been criticized in the past by Daganzo and others.     

Generalized conservative approximations of split convective derivative operators

We propose a strategy to design locally conservative finite-difference approximations of convective derivatives for shock-free compressible flows with arbitrary order of accuracy, that generalizes the approach of Ducros et al. (2000) [1], and that can be applied as a building block of low-dissipative, hybrid shock-capturing methods. The approximations stem from application of standard central difference formulas to split forms of the convective terms in the compressible Euler equations, which guarantee strong numerical stability and (near) energy preservation in the inviscid limit. A convenient implementation of the high-order fluxes is suggested, which guarantees improved computational efficiency over existing methods. Numerical tests performed for isotropic turbulence at zero viscosity show stability of schemes with order of accuracy up to ten, and effectiveness of convective splitting of Kennedy and Gruber (2008) [2] in providing extra stability in the presence of strong density variations. Numerical simulations of compressible turbulent boundary layer flow indicate suitability of the method for non-uniform grids, and overall support superior computational efficiency of high-order schemes.     

A high order moving boundary treatment for compressible inviscid flows

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax–Wendroff procedure proposed in [17] for conservation laws in static geometries. This procedure helps us obtain normal spatial derivatives at inflow boundaries from Lagrangian time derivatives and tangential derivatives by repeated use of the Euler equations. Together with high order extrapolation at outflow boundaries, we can impose accurate values of ghost points near the boundaries by a Taylor expansion. To maintain high order accuracy in time, we need some special time matching technique at the two intermediate Runge–Kutta stages. Numerical examples in one and two dimensions show that our boundary treatment is high order accurate for problems with smooth solutions. Our method also performs well for problems involving interactions between shocks and moving rigid bodies.     

求解Euler方程的高精度的空间--时间守恒格式

空间-时间守恒(STC)格式是近年来发展出的一种计算格式,在现有的STC格式构造过程中,流动变量在解元中的分布都用其一阶Taylor展开式来表示.STC格式的精度与所采用的Taylor展开式的阶数有关.该文采用流动变量的二阶Taylor展开式来表示其在解元上的分布、构造出了求解一维Euler方程的STC格式.用该格式对几个问题进行了计算,将计算结果与精确解进行了比较,比较表明该格式有较高的精度.     

Sixth order compact scheme combined with multigrid method and extrapolation technique for 2D poisson equation

We develop a sixth order finite difference discretization strategy to solve the two dimensional Poisson equation, which is based on the fourth order compact discretization, multigrid method, Richardson extrapolation technique, and an operator based interpolation scheme. We use multigrid V-Cycle procedure to build our multiscale multigrid algorithm, which is similar to the full multigrid method (FMG). The multigrid computation yields fourth order accurate solution on both the fine grid and the coarse grid. A sixth order accurate coarse grid solution is computed by using the Richardson extrapolation technique. Then we apply our operator based interpolation scheme to compute sixth order accurate solution on the fine grid. Numerical experiments are conducted to show the solution accuracy and the computational efficiency of our new method, compared to Sun–Zhang’s sixth order Richardson extrapolation compact (REC) discretization strategy using Alternating Direction Implicit (ADI) method and the standard fourth order compact difference (FOC) scheme using a multigrid method.     

Numerical experiments on the efficiency of local grid refinement based on truncation error estimates

Local grid refinement aims to optimise the relationship between accuracy of the results and number of grid nodes. In the context of the finite volume method no single local refinement criterion has been globally established as optimum for the selection of the control volumes to subdivide, since it is not easy to associate the discretisation error with an easily computable quantity in each control volume. Often the grid refinement criterion is based on an estimate of the truncation error in each control volume, because the truncation error is a natural measure of the discrepancy between the algebraic finite-volume equations and the original differential equations. However, it is not a straightforward task to associate the truncation error with the optimum grid density because of the complexity of the relationship between truncation and discretisation errors. In the present work several criteria based on a truncation error estimate are tested and compared on a regularised lid-driven cavity case at various Reynolds numbers. It is shown that criteria where the truncation error is weighted by the volume of the grid cells perform better than using just the truncation error as the criterion. Also it is observed that the efficiency of local refinement increases with the Reynolds number. The truncation error is estimated by restricting the solution to a coarser grid and applying the coarse grid discrete operator. The complication that high truncation error develops at grid level interfaces is also investigated and several treatments are tested.     

Compact local integrated-RBF approximations for second-order elliptic differential problems

This paper presents a new compact approximation method for the discretisation of second-order elliptic equations in one and two dimensions. The problem domain, which can be rectangular or non-rectangular, is represented by a Cartesian grid. On stencils, which are three nodal points for one-dimensional problems and nine nodal points for two-dimensional problems, the approximations for the field variable and its derivatives are constructed using integrated radial basis functions (IRBFs). Several pieces of information about the governing differential equation on the stencil are incorporated into the IRBF approximations by means of the constants of integration. Numerical examples indicate that the proposed technique yields a very high rate of convergence with grid refinement.     

Finite difference methods for second order in space,first order in time hyperbolic systems and the linear shifted wave equation as a model problem in numerical relativity

Motivated by the problem of solving the Einstein equations, we discuss high order finite difference discretizations of first order in time, second order in space hyperbolic systems. Particular attention is paid to the case when first order derivatives that can be identified with advection terms are approximated with non-centered finite difference operators. We first derive general properties of these discrete operators, then we extend a known result on numerical stability for such systems to general order of accuracy. As an application we analyze the shifted wave equation, including the behavior of the numerical phase and group speeds at different orders of approximations. Special attention is paid to when the use of off-centered schemes improves the accuracy over the centered schemes.     

Three-dimensional finite-difference lattice Boltzmann model and its application to inviscid compressible flows with shock waves

In this paper, a three-dimensional (3D) finite-difference lattice Boltzmann model for simulating compressible flows with shock waves is developed in the framework of the double-distribution-function approach. In the model, a density distribution function is adopted to model the flow field, while a total energy distribution function is adopted to model the temperature field. The discrete equilibrium density and total energy distribution functions are derived from the Hermite expansions of the continuous equilibrium distribution functions. The discrete velocity set is obtained by choosing the abscissae of a suitable Gauss–Hermite quadrature with sufficient accuracy. In order to capture the shock waves in compressible flows and improve the numerical accuracy and stability, an implicit–explicit finite-difference numerical technique based on the total variation diminishing flux limitation is introduced to solve the discrete kinetic equations. The model is tested by numerical simulations of some typical compressible flows with shock waves ranging from 1D to 3D. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.     

Simulation of axisymmetric jets with a finite element Navier–Stokes solver and a multilevel VOF approach

A multilevel VOF approach has been coupled to an accurate finite element Navier–Stokes solver in axisymmetric geometry for the simulation of incompressible liquid jets with high density ratios. The representation of the color function over a fine grid has been introduced to reduce the discontinuity of the interface at the cell boundary. In the refined grid the automatic breakup and coalescence occur at a spatial scale much smaller than the coarse grid spacing. To reduce memory requirements, we have implemented on the fine grid a compact storage scheme which memorizes the color function data only in the mixed cells. The capillary force is computed by using the Laplace–Beltrami operator and a volumetric approach for the two principal curvatures. Several simulations of axisymmetric jets have been performed to show the accuracy and robustness of the proposed scheme.     

Upwind schemes for the wave equation in second-order form

We develop new high-order accurate upwind schemes for the wave equation in second-order form. These schemes are developed directly for the equations in second-order form, as opposed to transforming the equations to a first-order hyperbolic system. The schemes are based on the solution to a local Riemann-type problem that uses d’Alembert’s exact solution. We construct conservative finite difference approximations, although finite volume approximations are also possible. High-order accuracy is obtained using a space-time procedure which requires only two discrete time levels. The advantages of our approach include efficiency in both memory and speed together with accuracy and robustness. The stability and accuracy of the approximations in one and two space dimensions are studied through normal-mode analysis. The form of the dissipation and dispersion introduced by the schemes is elucidated from the modified equations. Upwind schemes are implemented and verified in one dimension for approximations up to sixth-order accuracy, and in two dimensions for approximations up to fourth-order accuracy. Numerical computations demonstrate the attractive properties of the approach for solutions with varying degrees of smoothness.     

Classical and semirelativistic magnetohydrodynamics with anisotropic ion pressure

We study the magnetohydrodynamics (MHD) equations with anisotropic ion pressure and isotropic electron pressure under both the classical and semirelativistic approximations in order to develop a numerical model. The dispersion relation as well as the characteristic wave speeds are derived. In addition to the exact wave speed solutions, we also provide efficient approximate formulas for the semirelativistic magnetosonic speeds. The equations are discretized with the Rusanov and Harten-Lax-van Leer numerical schemes and implemented into the BATS-R-US MHD code. We perform a set of verification tests.     

Stochastic diagrams for critical point spectra

A new technique for calculating the time-evolution, correlations and steady state spectra for nonlinear stochastic differential equations is presented. To illustrate the method, we consider examples involving cubic nonlinearities in an N-dimensional phase-space. These serve as a useful paradigm for describing critical point phase transitions in numerous equilibrium and non-equilibrium systems, ranging from chemistry, physics and biology, to engineering, sociology and economics. The technique consists in developing the stochastic variable as a power series in time, and using this to compute the short time expansion for the correlation functions. This is then extrapolated to large times, and Fourier transformed to obtain the spectrum. Stochastic diagrams are developed to facilitate computation of the coefficients of the relevant power series expansion. Two different types of long-time extrapolation technique, involving either simple exponentials or logarithmic rational approximations, are evaluated for third-order diagrams. The analytical results thus obtained are compared with numerical simulations, together with exact results available in special cases. The agreement is found to be excellent up to and including the neighborhood of the critical point. Exponential extrapolation works especially well even above the critical point at large N values, where the dynamics is one of phase-diffusion in the presence of a spontaneously broken symmetry. This method also enables the calculation of the steady state spectra of polynomial functions of the stochastic variables. In these cases, the final correlations can be non-bistable even above threshold. Here logarithmic rational extrapolation has the greater accuracy of the two extrapolation methods. Stochastic diagrams are also applicable to more general problems involving spatial variation, in addition to temporal variation. Received: 12 January 1998