A representation of curved boundaries for the solution of the Navier–Stokes equations on a staggered three-dimensional Cartesian grid

A method is presented for representing curved boundaries for the solution of the Navier–Stokes equations on a non-uniform, staggered, three-dimensional Cartesian grid. The approach involves truncating the Cartesian cells at the boundary surface to create new cells which conform to the shape of the surface. We discuss in some detail the problems unique to the development of a cut cell method on a staggered grid. Methods for calculating the fluxes through the boundary cell faces, for representing pressure forces and for calculating the wall shear stress are derived and it is verified that the new scheme retains second-order accuracy in space. In addition, a novel “cell-linking” method is developed which overcomes problems associated with the creation of small cells while avoiding the complexities involved with other cell-merging approaches. Techniques are presented for generating the geometric information required for the scheme based on the representation of the boundaries as quadric surfaces. The new method is tested for flow through a channel placed oblique to the grid and flow past a cylinder at Re=40 and is shown to give significant improvement over a staircase boundary formulation. Finally, it is used to calculate unsteady flow past a hemispheric protuberance on a plate at a Reynolds number of 800. Good agreement is obtained with experimental results for this flow.     

Compact fourth-order finite volume method for numerical solutions of Navier–Stokes equations on staggered grids

The development of a compact fourth-order finite volume method for solutions of the Navier–Stokes equations on staggered grids is presented. A special attention is given to the conservation laws on momentum control volumes. A higher-order divergence-free interpolation for convective velocities is developed which ensures a perfect conservation of mass and momentum on momentum control volumes. Three forms of the nonlinear correction for staggered grids are proposed and studied. The accuracy of each approximation is assessed comparatively in Fourier space. The importance of higher-order approximations of pressure is discussed and numerically demonstrated. Fourth-order accuracy of the complete scheme is illustrated by the doubly-periodic shear layer and the instability of plane-channel flow. The efficiency of the scheme is demonstrated by a grid dependency study of turbulent channel flows by means of direct numerical simulations. The proposed scheme is highly accurate and efficient. At the same level of accuracy, the fourth-order scheme can be ten times faster than the second-order counterpart. This gain in efficiency can be spent on a higher resolution for more accurate solutions at a lower cost.     

A consistent and conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part III: On a staggered mesh

The consistent and conservative scheme developed on a rectangular collocated mesh [M.-J. Ni, R. Munipalli, N.B. Morley, P. Huang, M.A. Abdou, A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part I: on a rectangular collocated grid system, Journal of Computational Physics 227 (2007) 174–204] and on an arbitrary collocated mesh [M.-J. Ni, R. Munipalli, P. Huang, N.B. Morley, M.A. Abdou, A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part II: on an arbitrary collocated mesh, Journal of Computational Physics 227 (2007) 205–228] has been extended and specially designed for calculation of the Lorentz force on a staggered grid system (Part III) by solving the electrical potential equation for magnetohydrodynamics (MHD) at a low magnetic Reynolds number. In a staggered mesh, pressure (p) and electrical potential (φ) are located in the cell center, while velocities and current fluxes are located on the cell faces of a main control volume. The scheme numerically meets the physical conservation laws, charge conservation law and momentum conservation law. Physically, the Lorentz force conserves the momentum when the magnetic field is constant or spatial coordinate independent. The calculation of current density fluxes on cell faces is conducted using a scheme consistent with the discretization for solution of the electrical potential Poisson equation, which can ensure the calculated current density conserves the charge. A divergence formula of the Lorentz force is used to calculate the Lorentz force at the cell center of a main control volume, which can numerically conserve the momentum at constant or spatial coordinate independent magnetic field. The calculated cell-center Lorentz forces are then interpolated to the cell faces, which are used to obtain the corresponding velocity fluxes by solving the momentum equations. The “conservative” is an important property of the scheme, which can guarantee computational accuracy of MHD flows at high Hartmann number with a strongly non-uniform mesh employed to resolve the Hartmann layers and side layers. 2D fully developed MHD flows with analytical solutions available have been conducted to validate the scheme at a staggered mesh. 3D MHD flows, with the experimental data available, at a constant magnetic field in a rectangular duct with sudden expansion and at a varying magnetic field in a rectangular duct are conducted on a staggered mesh to verify the computational accuracy of the scheme. It is expected that the scheme for the Lorentz force can be employed together with a fully conservative scheme for the convective term and the pressure term [Y. Morinishi, T.S. Lund, O.V. Vasilyev, P. Moin, Fully conservative higher order finite difference schemes for incompressible flow, Journal of Computational Physics 143 (1998) 90–124] for direct simulation of MHD turbulence and MHD instability with good accuracy at a staggered mesh.     

数值模拟水环境污染的一种L∞稳定的分步杂交方法

本文提出了一种数值求解对流扩散方程的分步杂交方法。在不规则的三角形网格上,采用迎风离散格式或改型特征线方法处理对流算子;采用集中质量的有限元方法处理扩散算子。详细分析了这种算法的稳定性同题,在数学上严格证明了在满足①Δt≤min((2d)/v,(d²)/(3K)),其中d是三角形网格中最短垂线的长度,V和K分别为流场中的最大速度和扩散系数。②所有三角形的内角θ≤π/2的条件下,整个计算格式是L_∞稳定的,从而保证了在海洋环境和水质的数值模拟中海水的盐度、污染物的浓度和核电站冷却水系统中的超温不会出现负值。应用非线性的对流扩散方程对此方法的精度和收敛性进行了检验。通过数值解与精确解的比较,表明本方法的数值耗散很小,用改型特征线方法处理对流算子较迎风离散格式有更高的精度;两种处理对流算子的方法都没有伪振荡现象发生。本方法由于具有算法简单、L_∞稳定、计算网格灵活等优点,可推广使用于实际的海洋环境(潮波、海流、海洋污染)、港口和海湾的数值模拟以及不可压粘性流和对流传热同题的数值计算。     

A new combined stable and dispersion relation preserving compact scheme for non-periodic problems

A new compact scheme is presented for computing wave propagation problems and Navier–Stokes equation. A combined compact difference scheme is developed for non-periodic problems (called NCCD henceforth) that simultaneously evaluates first and second derivatives, improving an existing combined compact difference (CCD) scheme. Following the methodologies in Sengupta et al. [T.K. Sengupta, S.K. Sircar, A. Dipankar, High accuracy schemes for DNS and acoustics, J. Sci. Comput. 26 (2) (2006) 151–193], stability and dispersion relation preservation (DRP) property analysis is performed here for general CCD schemes for the first time, emphasizing their utility in uni- and bi-directional wave propagation problems – that is relevant to acoustic wave propagation problems. We highlight: (a) specific points in parameter space those give rise to least phase and dispersion errors for non-periodic wave problems; (b) the solution error of CCD/NCCD schemes in solving Stommel Ocean model (an elliptic p.d.e.) and (c) the effectiveness of the NCCD scheme in solving Navier–Stokes equation for the benchmark lid-driven cavity problem at high Reynolds numbers, showing that the present method is capable of providing very accurate solution using far fewer points as compared to existing solutions in the literature.     

高阶辛算法的稳定性与数值色散性分析

利用Maxwell方程的哈密尔顿函数,导出对应的欧拉-哈密尔顿方程.利用辛积分技术与高阶交错差分技术,建立求解三维时域Maxwell方程的高阶辛算法;结合电磁场中的物理概念,借助矩阵分析和张量分析理论,获得高阶时域方法及高阶辛算法的稳定性和数值色散性的统一处理新方法.用数值结果证实方法的正确性,与FDTD算法和其它时域高阶方法相比,高阶辛算法具有较大的计算优势,为电磁计算提供了新的途径.     

Simulation of natural convection under high magnetic field by means of the thermal lattice Boltzmann method

The thermal lattice Boltzmann method (TLBM), which was proposed by J. G. M. Eggels and J. A. Somers previously, has been improved in this paper. The improved method has introduced a new equilibrium solution for the temperature distribution function on the assumption that flow is incompressible, and it can correct the effect of compressibility on the macroscopic temperature computed. Compared to the previous method, where the half-way bounce back boundary condition was used for non-slip velocity and temperature, a non-equilibrium extrapolation scheme has been adopted for both velocity and temperature boundary conditions in this paper. Its second-order accuracy coincides with the ensemble accuracy of lattice Boltzmann method. In order to validate the improved thermal scheme, the natural convection of air in a square cavity is simulated by using this method. The results obtained in the simulation agree very well with the data of other numerical methods and benchmark data. It is indicated that the improved TLBM is also successful for the simulations of non-isothermal flows. Moreover, this thermal scheme can be applied to simulate the natural convection in a non-uniform high magnetic field. The simulation has been completed in a square cavity filled with the aqueous solutions of KCl (11wt%), which is considered as a diamagnetic fluid with electrically low-conducting, with Grashof number Gr=4.64× 10^4 and Prandtl number Pr=7.0. And three cases, with different cavity locations in the magnetic field, have been studied. In the presence of a high magnetic field, the natural convection is quenched by the body forces exerted on the electrically low-conducting fluids, such as the magnetization force and the Lorentz force. From the results obtained, it can be seen that the quenching efficiencies decrease with the variation of location from left, symmetrical line, to the right. These phenomena originate from the different distributions of the magnetic field strengths in the zones of the symmetrical central line of the magnetic fields. The results are also compared with those without a magnetic field. Finally, we can conclude that the improved TLBM will enable effective simulation of the natural convection under a high magnetic field.     

A local directional ghost cell approach for incompressible viscous flow problems with irregular boundaries

An immersed boundary method for the incompressible Navier–Stokes equations in irregular domains is developed using a local ghost cell approach. This method extends the solution smoothly across the boundary in the same direction as the discretization it will be used for. The ghost cell value is determined locally for each irregular grid cell, making it possible to treat both sharp corners and thin plates accurately. The time stepping is done explicitly using a second order Runge–Kutta method. The spatial derivatives are approximated by finite difference methods on a staggered, Cartesian grid with local grid refinements near the immersed boundary. The WENO scheme is used to treat the convective terms, while all other terms are discretized with central schemes. It is demonstrated that the spatial accuracy of the present numerical method is second order. Further, the method is tested and validated for a number of problems including uniform flow past a circular cylinder, impulsively started flow past a circular cylinder and a flat plate, and planar oscillatory flow past a circular cylinder and objects with sharp corners, such as a facing square and a chamfered plate.     

自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

一种基于AUSM分裂的真正多维HLL格式

文章给出了一种真正多维的HLL Riemann解算器.采用AUSM分裂将通量分解成为对流通量和压力通量, 其中对流通量的计算采用迎风格式, 压力通量的计算采用HLL格式, 且将HLL格式的耗散项中的密度差用压力差代替, 从而使得格式能够分辨接触间断.为了实现数值格式真正多维的特性, 分别计算了网格界面中点和角点上的数值通量, 并且采用Simpson公式加权组合中点和角点上的数值通量得到网格界面的数值通量.为了减少重构角点处状态时的模板宽度, 计算中采用基于SDWLS梯度的线性重构获得2阶空间精度, 而时间离散采用2阶保强稳Runge-Kutta方法.数值实验表明, 相比于传统的一维HLL格式, 文章的真正多维HLL格式具有能够分辨接触间断, 以及更大的时间步长等优点.与其他能够分辨接触间断的格式(例如HLLC格式)不同, 真正多维的HLL格式在计算二维问题时不会出现激波不稳定现象.      

High order finite difference methods with subcell resolution for advection equations with stiff source terms

A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers.     

基于SIMPLE的紧致方法

通过泰勒展式系数匹配的方法发展了基于非等距网格的有限容积紧致格式,采用延迟修正的方法建立了基于SIMPLE的紧致方法,,该方法能够得到高精度的数值解,增加迭代求解代数方程组的稳定性。对底部加热的方腔内自然对流换热问题进行数值模拟,结果表明,紧致方法比二阶中心差分方法具有更高的精度。     

一种基于波动方程的新高阶FDTD方法

 提出了一种基于二阶波动方程的（2M,4）高阶时域有限差分(FDTD)方法,通过使用辛积分传播子(SIP)在时域上获得4阶精度,使用离散奇异卷积(DSC)方法在空域上达到2M阶精度。与已有的(2M,4) 阶FDTD方法相比,虽然两者都采用SIP和DSC方法,但是此二者的不同点在于：第一,新方法基于二阶波动方程;第二,在离散计算空间时使用单一网格而不是传统的Yee网格;第三,单独计算某一场分量从而节约内存并减少计算量。数值计算结果表明,与传统高阶算法相比,基于波动方程的高阶FDTD方法耗费的机时只有它的50%,内存消耗下降10%, 而两者的计算结果之间相对误差小于5‰。     

求解Navier-Stokes方程组的组合紧致迎风格式

给出一种新的至少有四阶精度的组合紧致迎风(CCU)格式,该格式有较高的逼近解率,利用该组合迎风格式,提出一种新的适合于在交错网格系统下求解Navier-Stokes方程组的高精度紧致差分投影算法.用组合紧致迎风格式离散对流项,粘性项、压力梯度项以及压力Poisson方程均采用四阶对称型紧致差分格式逼近,算法的整体精度不低于四阶.通过对Taylor涡列、对流占优扩散问题和双周期双剪切层流动问题的计算表明,该算法适合于对复杂流体流动问题的数值模拟.     

A method of solving the stiffness problem in Biot＇s poroelastic equations using a staggered high-order finite-difference

In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.     

用低精度CCD获得高精度测量方法的研究

为了大幅提高线阵CCD的测量精度,提出了一种全新的CCD使用方法。该方法是将N个像元间距为H的线阵CCD器件许排组合住一起,并沿像元线性分布方向以距离为H/N依次均匀错开排列。多个线阵CCD的感光电信号经多通道模一数同步采集,保存到存储器中指定位置。然后,通过对所有CCD测量数据的分析计算来获得精确的测量值。分别采用单CCD和双CCD错排对长为30mm,直径为5．000mm、8．000mm、12．000mm的三个标准杆件的直径进行了测量。结果表明,蚁CCD错排可获得两倍于单CCD的测量精度。该方法可从理论上彻底打破CCD像元问距的限制,并使线阵CCD的测量精度大幅度地提高。     

A new high-order discontinuous Galerkin spectral finite element method for Lagrangian gas dynamics in two-dimensions

This paper presents a new high-order cell-centered Lagrangian scheme for two-dimensional compressible flow. The scheme uses a fully Lagrangian form of the gas dynamics equations, which is a weakly hyperbolic system of conservation laws. The system of equations is discretized in the Lagrangian space by discontinuous Galerkin method using a spectral basis. The vertex velocities and the numerical fluxes through the cell interfaces are computed consistently in the Eulerian space by virtue of an improved nodal solver. The nodal solver uses the HLLC approximate Riemann solver to compute the velocities of the vertex. The time marching is implemented by a class of TVD Runge–Kutta type methods. A new HWENO (Hermite WENO) reconstruction algorithm is developed and used as limiters for RKDG methods to maintain compactness of RKDG methods. The scheme is conservative for the mass, momentum and total energy. It can maintain high-order accuracy both in space and time, obey the geometrical conservation law, and achieve at least second order accuracy on quadrilateral meshes. Results of some numerical tests are presented to demonstrate the accuracy and the robustness of the scheme.     

三维定常问题的高阶紧致差分格式向非定常问题的推广

将已经建立的求解三维定常对流扩散方程的高阶紧致差分格式直接推广到三维非定常对流扩散方程的数值求解,时间导数项利用二阶向后欧拉差分公式,所得到的高阶隐式紧致差分格式时间为二阶精度,空间为四阶精度,并且是无条件稳定的.数值实验结果验证了本文方法的精确性和稳健性.     

紧致方法对流动换热及静态分岔的模拟

发展了基于投影法的紧致方法求解流动换热问题,对顶盖驱动流和侧壁加热的方腔内自然对流换热问题进行了数值模拟。与其它传统方法相比,紧致方法能在较少的网格结点下获得精度较高的计算结果。进一步,采用所发展的紧致方法对不同工况下的Rayleigh-Benard对流及其静态分岔现象进行了数值模拟。数值计算结果表明当长宽比变大时,底部努塞尔数会有小幅度增加。当长宽比为8时,用所发展的紧致方法不同的初场可以得出三种不同的流场和温度场。与基于QUICK格式的SIMPLE算法相比,所发展的紧致方法可以多预测一种静态分岔现象。     

铸件在凝固过程中的传热和自然对流计算

论文中以流函数、涡度和温度为基本变量,建立了铸件凝固时的传热和自然对流的数学模型。控制方程采用半隐式有限差分在固定网格上求解。在模型建立过程中,着重讨论了潜热的处理和固液相区(界面)速度条件。对于潜热处理,采用了基于合金平衡相图的比热焓法;对于固液相区(界面)速度条件,提出了一种新的、基于固相率进行修正的处理方法。另外还探讨了将本文提出的解法应用于复杂几何条件的可能性。最后,通过对一个典型例子进行计算,对自然对流对铸件凝固的影响作了深入的讨论。