Cavitating/non-cavitating flows simulation by third-order finite volume scheme and power-law preconditioning method

Equations of steady inviscid and laminar flows are solved by means of a third-order finite volume （FV） scheme. For this purpose, a cell-centered discretization technique is employed. In this technique, the flow parameters at the cell faces are computed using a third-order weighted averages procedure. A fourth-order artificial dissipation is used for stability of the solution. In order to achieve the steady-state situation, four-step Runge-Kutta explicit time integration method is applied. An advanced progressive preconditioning method, named the power-law preconditioning method, is used for faster convergence. In this method, the preconditioning matrix is adjusted automatically from the velocity and/or pressure flow-field by a power-law relation. Attention is directed towards accuracy and convergence of the schemes. The results presented in the paper focus on steady inviscid and laminar flows around sheet-cavitating and fully-wetted bodies including hydrofoils and circular/elliptical cylinder. Excellent agreements are obtained when numerical predictions are compared with other available experimental and numerical results. In addition, it is found that using the power-law preconditioner significantly increases the numerical convergence speed.

Existence and iteration of monotone positive solutions for a third-order two-point boundary value problem

The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″（t） ＋ q（t）f（t,u（t）,u′（t）） = 0, 0 〈 t 〈 1, u（0） = u″（0） = u′（1） = 0 is studied. The iterative schemes for approximating the solutions are obtained by applying a monotone iterative method.

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.

一种三阶WENO-Z格式改进方法

高分辨率激波捕捉格式对含激波流场的数值模拟具有重要意义。在三阶WENO-Z格式(WENO-Z3)基础上,通过构造不同形式的全局光滑因子得到WENO-Z3N1、WENO-Z3N2、WENO-Z3N3格式。选取Sod激波管、双爆轰波碰撞、激波与熵波相互作用等经典算例,考察了3种格式(WENO-Z3N1、WENO-Z3N2、WENO-Z3N3)的计算性能。根据泰勒级数展开,理论推导给出3种格式的精度分析。通过探讨各格式理论精度与实际计算精度之间的关系得到如下结论:3种格式在连续解非极值点处的理论精度对实际计算性能起决定性的作用,并通过双马赫反射问题进一步验证该结论的可靠性。本文的研究给出一种三阶WENO-Z格式的改进方法,合理构造全局光滑因子使得格式在连续解非极值点处满足设计精度的要求。

Third-order unconditional positivity-preserving schemes for reactive flows keeping both mass and mole balance

In this paper,the previously proposed second-order process-based modified Patankar Runge-Kutta schemes are extended to the third order of accuracy.Owing to the process-based implicit handling of reactive source terms,the mass conservation,mole balance and energy conservation are kept simultaneously while the positivity for the density and pressure is preserved unconditionally even with stiff reaction networks.It is proved that the first-order truncation terms for the Patankar coefficients must be zero to achieve a prior third order of accuracy for most cases.A twostage Patankar procedure for each Runge-Kutta step is designed to eliminate the first-order truncation terms,accomplish the prior third order of accuracy and maximize the Courant number which the total variational diminishing property requires.With the same approach as the second-order schemes,the third-order ones are applied to Euler equations with chemical reactive source terms.Numerical studies including both 1D and 2D ordinary and partial differential equations are conducted to affirm both the prior order of accuracy and the positivity-preserving property for the density and pressure.

DISCRETE ENERGY ANALYSIS OF THE THIRD-ORDER VARIABLE-STEP BDF TIME-STEPPING FOR DIFFUSION EQUATIONS

This is one of our series works on discrete energy analysis of the variable-step BDF schemes.In this part,we present stability and convergence analysis of the third-order BDF(BDF3)schemes with variable steps for linear diffusion equations,see,e.g.,[SIAM J.Numer.Anal.,58:2294-2314]and[Math.Comp.,90:1207-1226]for our previous works on the BDF2 scheme.To this aim,we first build up a discrete gradient structure of the variable-step BDF3 formula under the condition that the adjacent step ratios are less than 1.4877,by which we can establish a discrete energy dissipation law.Mesh-robust stability and convergence analysis in the L^(2) norm are then obtained.Here the mesh robustness means that the solution errors are well controlled by the maximum time-step size but independent of the adjacent time-step ratios.We also present numerical tests to support our theoretical results.

An improved third-order HWCNS for compressible flow simulation on curvilinear grids

Due to the very high requirements on the quality of computational grids,stability property and computational efficiency,the application of high-order schemes to complex flow simulation is greatly constrained.In order to solve these problems,the third-order hybrid cell-edge and cell-node weighted compact nonlinear scheme(HWCNS3)is improved by introducing a new nonlinear weighting mechanism.The new scheme uses only the central stencil to reconstruct the cell boundary value,which makes the convergence of the scheme more stable.The application of the scheme to Euler equations on curvilinear grids is also discussed.Numerical results show that the new HWCNS3 achieves the expected order in smooth regions,captures discontinuities sharply without obvious oscillation,has higher resolution than the original one and preserves freestream and vortex on curvilinear grids.