气相爆轰高阶中心差分-WENO组合格式自适应网格方法

研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点.     

A class of finite difference schemes with low dispersion and controllable dissipation for DNS of compressible turbulence

In this paper, a class of finite difference schemes which achieves low dispersion and controllable dissipation in smooth region and robust shock-capturing capabilities in the vicinity of discontinuities is presented. Firstly, a sufficient condition for semi-discrete finite difference schemes to have independent dispersion and dissipation is derived. This condition enables a novel approach to separately optimize the dissipation and dispersion properties of finite difference schemes and a class of schemes with minimized dispersion and controllable dissipation is thus obtained. Secondly, for the purpose of shock-capturing, one of these schemes is used as the linear part of the WENO scheme with symmetrical stencils to constructed an improved WENO scheme. At last, the improved WENO scheme is blended with its linear counterpart to form a new hybrid scheme for practical applications. The proposed scheme is accurate, flexible and robust. The accuracy and resolution of the proposed scheme are tested by the solutions of several benchmark test cases. The performance of this scheme is further demonstrated by its application in the direct numerical simulation of compressible turbulent channel flow between isothermal walls.     

Increasing the accuracy in locally divergence-preserving finite volume schemes for MHD

It is of utmost interest to control the divergence of the magnetic flux in simulations of the ideal magnetohydrodynamic equations since, in general, divergence errors tend to accumulate and render the schemes unstable. This paper presents a higher-order extension of the locally divergence-preserving procedure developed in Torrilhon [M. Torrilhon, Locally divergence-preserving upwind finite volume schemes for magnetohydrodynamic equations, SIAM J. Sci. Comput. 26 (2005) 1166–1191]; a fourth-order accurate local redistribution of the numerical magnetic field fluxes of a finite volume base scheme is introduced. The redistribution ensures that a fourth-order accurate discrete divergence operator is preserved to round off errors when applied to the cell averages of the magnetic flux density. The developed procedure is applicable to generic semi-discrete finite volume schemes and its purpose is to stabilize the schemes using a local procedure that respects the accuracy of the base scheme to a greater extent than the previous second-order achievements. Numerical experiments that demonstrate the properties of the new procedure are also presented.     

三维边界层内定常横流涡的感受性研究

层流向湍流转捩的预测与控制一直是研究的前沿热点问题之一,其中感受性阶段是转捩过程中的初始阶段,它决定着湍流产生或形成的物理过程.但是有关三维边界层内感受性问题的数值和理论研究都比较少;实际工程问题中大部分转捩过程都是发生在三维边界层流中,所以研究三维边界层中的感受性问题显得尤为重要.本文以典型的后掠角45?无限长平板为例,数值研究了在三维壁面局部粗糙作用下的三维边界层感受性问题,探讨了三维边界层感受性问题与三维壁面局部粗糙长、宽和高之间的关系;然后,考虑在后掠平板上设计不同的三维壁面局部粗糙的分布状态、几何形状、距离后掠平板前缘的位置以及流向和展向设计多个三维壁面局部粗糙对三维边界层感受性问题有何影响;最后,讨论两两三维壁面局部粗糙中心点之间的距离以及后掠角的改变对三维边界层感受性的物理过程将会发生何种影响等.这一问题的深入研究将为三维边界层流中层流向湍流转捩过程的认识和理解提供理论依据.     

有限差分格式的一个通用色散-耗散条件

通过分析显式有限差分格式的数值色散和数值耗散,导出一个适于有限差分格式的通用色散-耗散条件.根据群速度和耗散率之间的物理关系,确定了用以抑制数值解中伪高波数波所需要的适度耗散.在以往发展的低耗散加权基本无振荡格式WENO-CU6-M2上的应用表明,该条件可用作优化线性或非线性有限差分格式的色散和耗散的通用指导准则.此外,满足色散-耗散条件的改进WENO-CU6-M2格式还可选作低分辨率数值模拟,以三维Taylor-Green涡向湍流转捩和自相似能量衰减问题展现了它的这种能力.与经典的动态Smagorinsky亚网格尺度模型相比,在Reynolds数Re=400~3000条件下,无黏和黏性Taylor-Green涡的数值模拟结果均得到明显改善.在保持激波捕捉特性同时,与最新的隐式大涡模拟模型的计算效果相当.      

The numerical prediction of planar viscoelastic contraction flows using the pom–pom model and higher-order finite volume schemes

This study investigates the numerical solution of viscoelastic flows using two contrasting high-order finite volume schemes. We extend our earlier work for Poiseuille flow in a planar channel and the single equation form of the extended pom–pom (SXPP) model [M. Aboubacar, J.P. Aguayo, P.M. Phillips, T.N. Phillips, H.R. Tamaddon-Jahromi, B.A. Snigerev, M.F. Webster, Modelling pom–pom type models with high-order finite volume schemes, J. Non-Newtonian Fluid Mech. 126 (2005) 207–220], to determine steady-state solutions for planar 4:1 sharp contraction flows. The numerical techniques employed are time-stepping algorithms: one of hybrid finite element/volume type, the other of pure finite volume form. The pure finite volume scheme is a staggered-grid cell-centred scheme based on area-weighting and a semi-Lagrangian formulation. This may be implemented on structured or unstructured rectangular grids, utilising backtracking along the solution characteristics in time. For the hybrid scheme, we solve the momentum-continuity equations by a fractional-staged Taylor–Galerkin pressure-correction procedure and invoke a cell-vertex finite volume scheme for the constitutive law. A comparison of the two finite volume approaches is presented, concentrating upon the new features posed by the pom–pom class of models in this context of non-smooth flows. Here, the dominant feature of larger shear and extension in the entry zone influences both stress and stretch, so that larger stretch develops around the re-entrant corner zone as Weissenberg number increases, whilst correspondingly stress levels decline.     

A systematic methodology for constructing high-order energy stable WENO schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by the authors of this paper [N.K. Yamaleev, M.H. Carpenter, Third-order energy stable WENO scheme, J. Comput. Phys. 228 (2009) 3025–3047] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables “energy stable” modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes including one sixth-order central scheme; ESWENO schemes up to eighth order are presented in the Appendix. We also develop new weight functions and derive constraints on their parameters, which provide consistency, much faster convergence of the high-order ESWENO schemes to their underlying linear schemes for smooth solutions with arbitrary number of vanishing derivatives, and better resolution near strong discontinuities than the conventional counterparts.     

A Linear Energy-Preserving Finite Volume Element Method for the Improved Korteweg−de Vries Equation

In this paper, we introduce a linearized energy-preserving scheme which preserves the discrete global energy of solutions to the improved Korteweg?deVries equation. The method presented is based on the finite volume element method, by resorting to the variational derivative to transform the improved Korteweg?deVries equation into a new form, and then designing energy-preserving schemes for the transformed equation. The proposed scheme is much more efficient than the standard nonlinear scheme and has good stability. To illustrate its efficiency and conservative properties, we also compare it with other nonlinear schemes. Finally, we verify the efficiency and conservative properties through numerical simulations.     

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.     

Energy-Conserving Lattice Boltzmann Thermal Model in Two Dimensions

A discrete velocity model is presented for lattice Boltzmann thermal fluid dynamics. This model is implemented and tested in two dimensions with a finite difference scheme. Comparison with analytical solutions shows an excellent agreement even for wide temperature differences. An alternative approximate approach is then presented for traditional lattice transport schemes     

DGM-FD: A finite difference scheme based on the discontinuous Galerkin method applied to wave propagation

In this paper we formulate a numerical method that is high order with strong accuracy for numerical wave numbers, and is adaptive to non-uniform grids. Such a method is developed based on the discontinuous Galerkin method (DGM) applied to the hyperbolic equation, resulting in finite difference type schemes applicable to non-uniform grids. The schemes will be referred to as DGM-FD schemes. These schemes inherit naturally some features of the DGM, such as high-order approximations, applicability to non-uniform grids and super-accuracy for wave propagations. Stability of the schemes with boundary closures is investigated and validated. Proposed scheme is demonstrated by numerical examples including the linearized acoustic waves and solutions of non-linear Burger’s equation and the flat-plate boundary layer problem. For non-linear equations, proposed flux finite difference formula requires no explicit upwind and downwind split of the flux. This is in contrast to existing upwind finite difference schemes in the literature.     

A splitting approach for the fully nonlinear and weakly dispersive Green–Naghdi model

The fully nonlinear and weakly dispersive Green–Naghdi model for shallow water waves of large amplitude is studied. The original model is first recast under a new formulation more suitable for numerical resolution. An hybrid finite volume and finite difference splitting approach is then proposed, which could be adapted to many physical models that are dispersive corrections of hyperbolic systems. The hyperbolic part of the equations is handled with a high-order finite volume scheme allowing for breaking waves and dry areas. The dispersive part is treated with a classical finite difference approach. Extensive numerical validations are then performed in one horizontal dimension, relying both on analytical solutions and experimental data. The results show that our approach gives a good account of all the processes of wave transformation in coastal areas: shoaling, wave breaking and run-up.     

Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations

In this paper, we propose a finite difference AMR-WENO method for hyperbolic conservation laws. The proposed method combines the adaptive mesh refinement (AMR) framework  and  with the high order finite difference weighted essentially non-oscillatory (WENO) method in space and the total variation diminishing (TVD) Runge–Kutta (RK) method in time (WENO-RK)  and  by a high order coupling. Our goal is to realize mesh adaptivity in the AMR framework, while maintaining very high (higher than second) order accuracy of the WENO-RK method in the finite difference setting. The high order coupling of AMR and WENO-RK is accomplished by high order prolongation in both space (WENO interpolation) and time (Hermite interpolation) from coarse to fine grid solutions, and at ghost points. The resulting AMR-WENO method is accurate, robust and efficient, due to the mesh adaptivity and very high order spatial and temporal accuracy. We have experimented with both the third and the fifth order AMR-WENO schemes. We demonstrate the accuracy of the proposed scheme using smooth test problems, and their quality and efficiency using several 1D and 2D nonlinear hyperbolic problems with very challenging initial conditions. The AMR solutions are observed to perform as well as, and in some cases even better than, the corresponding uniform fine grid solutions. We conclude that there is significant improvement of the fifth order AMR-WENO over the third order one, not only in accuracy for smooth problems, but also in its ability in resolving complicated solution structures, due to the very low numerical diffusion of high order schemes. In our work, we found that it is difficult to design a robust AMR-WENO scheme that is both conservative and high order (higher than second order), due to the mass inconsistency of coarse and fine grid solutions at the initial stage in a finite difference scheme. Resolving these issues as well as conducting comprehensive evaluation of computational efficiency constitute our future work.     

非结构网格上的迎风有限元格式

在非结构网格上提出一种基于修正积分区域的迎风有限元格式,它与一阶迎风差分格式相当,可应用于构造各种不同的数值格式。     

Skew-symmetric form of convective terms and fully conservative finite difference schemes for variable density low-Mach number flows

The form of convective terms for compressible flow equations is discussed in the same way as for an incompressible one by Morinishi et al. [Y. Morinishi, T.S. Lund, O.V. Vasilyev, P. Moin, Fully conservative higher order finite difference schemes for incompressible flow, J. Comput. Phys. 124 (1998) 90], and fully conservative finite difference schemes suitable for shock-free unsteady compressible flow simulations are proposed. Commutable divergence, advective, and skew-symmetric forms of convective terms are defined by including the temporal derivative term for compressible flow. These forms are analytically equivalent if the continuity is satisfied, and the skew-symmetric form is secondary conservative without the aid of the continuity, while the divergence form is primary conservative. The relations between the present and existing quasi-skew-symmetric forms are also revealed. Commutable fully discrete finite difference schemes of convection are then derived in a staggered grid system, and they are fully conservative provided that the corresponding discrete continuity is satisfied. In addition, a semi-discrete convection scheme suitable for compact finite difference is presented based on the skew-symmetric form. The conservation properties of the present schemes are demonstrated numerically in a three-dimensional periodic inviscid flow. The proposed fully discrete fully conservative second-order accurate scheme is also used to perform the DNS of compressible isotropic turbulence and the simulation of open cavity flow.     

混合通量差分格式的比较

系统研究了几种混合通量差分格式的构造方法和耗散模型,分别对低速平板绕流、二维跨音速喷管流动和高超音速钝头体无粘绕流进行了数值模拟,结合先进的EASM湍流模型对格式的粘性分辨率和激波稳定性进行了细致的比较分析.结果表明混合通量差分格式兼顾了FDS和FVS格式的优点,具有较高的间断分辨率和数值稳定性.     

Analysis of numerical dispersion properties of SFDTD and MRTD schemes

In this paper, the Maxwell's equations are written as Hamilton canonical equations by using Hamilton functional variation method. Maxwell's equations can be discretized with symplectic propagation technique combined with high-order difference schemes approximations to construct symplectic finite difference time domain (SFDTD) method. The high-order dispersion equations of the scheme for space is deduced. The numerical dispersion analysis is included, and it is compared with the multiresolution time-domain (MRTD) method based on the Daubechies scaling functions. Numerical results show high efficiency and accuracy of the SFDTD method.     

A New Finite Volume Element Formulation for the Non-Stationary Navier-Stokes Equations

A semi-discrete scheme about time for the non-stationary Navier-Stokes equations is presented firstly, then a new fully discrete finite volume element (FVE) formulation based on macroelement is directly established from the semi-discrete scheme about time. And the error estimates for the fully discrete FVE solutions are derived by means of the technique of the standard finite element method. It is shown by numerical experiments that the numerical results are consistent with theoretical conclusions. Moreover, it is shown that the FVE method is feasible and efficient for finding the numerical solutions of the non-stationary Navier-Stokes equations and it is one of the most effective numerical methods among the FVE formulation, the finite element formulation, and the finite difference scheme.     

基于局部基本解法的静电场仿真分析

将局部基本解方法应用于静电场问题的模拟与分析。局部基本解方法是利用控制方程的基本解,基于局部理论和移动最小二乘原理提出的一种无网格算法。相比于有限元和有限差分等传统网格类方法,该方法仅需离散节点,避免了复杂的网格剖分难题。作为一种半解析数值技术,物理问题的基本解被作为插值基函数建立数值离散模型,从而保证了算法的较高精度。此外,与具有全局离散格式的无网格方法相比,局部基本解法更适用于高维复杂几何和大尺度模拟。二维和三维数值试验表明,该方法具有实施方便灵活,计算精度高和计算速度快等优势。为静电场仿真研究开辟新的途径,拓展了局部基本解方法的应用领域。     

Full time-domain nonlinear coupled dynamic analysis of a truss spar and its mooring/riser system in irregular wave

A new full time-domain nonlinear coupled method has been established and then applied to predict the responses of a Truss Spar in irregular wave.For the coupled analysis,a second-order time-domain approach is developed to calculate the wave forces,and a finite element model based on rod theory is established in three dimensions in a global coordinate system.In numerical implementation,the higher-order boundary element method(HOBEM)is employed to solve the velocity potential,and the 4th-order Adams-Bashforth-Moultn scheme is used to update the second-order wave surface.In deriving convergent solutions,the hull displacements and mooring tensions are kept consistent at the fairlead and the motion equations of platform and mooring-lines/risers are solved simultaneously using Newmark-integration scheme including Newton-Raphson iteration.Both the coupled quasi-static analysis and the coupled dynamic analysis are performed.The numerical simulation results are also compared with the model test results,and they coincide very well as a whole.The slow-drift responses can be clearly observed in the time histories of displacements and mooring tensions.Some important characteristics of the coupled responses are concluded.