An important property for finite difference schemes designed on curvilinear meshes is the exact preservation of free-stream solutions. This property is difficult to fulfill for high order conservative essentially non-oscillatory (WENO) finite difference schemes. In this paper we explore an alternative flux formulation for such finite difference schemes which can preserve free-stream solutions, based on the numerical technique for the metric terms, which can be applied to this alternative flux formulation but is difficult to be applied to the standard finite difference formulation. Free-stream and vortex preservation properties are investigated, and comparison with standard finite difference WENO schemes is made. Theoretical derivation and numerical results show that the finite difference WENO schemes based on the alternative flux formulation can preserve free-stream and vortex solutions on both stationary and dynamically generalized coordinate systems, hence giving much better performance than the standard finite difference WENO schemes for such problems.

No abstract uploaded!

How to properly specify boundary conditions for pressure is a longstanding problem for the incompressible Navier–Stokes equations with no-slip boundary conditions. An analytical resolution of this issue stems from a recently developed formula for the pressure in terms of the commutator of the Laplacian and Leray projection operators. Here we make use of this formula to (a) improve the accuracy of computing pressure in two kinds of existing time-discrete projection methods implicit in viscosity only, and (b) devise new higher-order accurate time-discrete projection methods that extend a slip-correction idea behind the well-known finite-difference scheme of Kim and Moin. We test these schemes for stability and accuracy using various combinations of C0 finite elements. For all three kinds of time discretization, one can obtain third-order accuracy for both pressure and velocity without a time-step stability restriction of diffusive type. Furthermore, two kinds of projection methods are found stable using piecewise-linear elements for both velocity and pressure.

This paper gives a classification of first order polynomial differential operators of form $\mathscr{X} = X_1(x_1,x_2)\delta_1 + X_2(x_1,x_2)\delta_2$, $(\delta_i = \partial/\partial x_i)$. The classification is given through the order of an operator that is defined in this paper. Let $X=\mathscr{X}y$ to be the differential polynomial associated with $\mathscr{X}$, the order of $\mathscr{X}$, $\mathrm{ord}(\mathscr{X})$, is defined as the order of a differential ideal $\Lambda$ of differential polynomials that is a nontrivial expansion of the ideal $\{X\}$ and with the lowest order. In this paper, we prove that there are only four possible values for the order of a differential operator, $0$, $1$, $2$, $3$, or $\infty$. Furthermore, when the order is finite, the expansion $\Lambda$ is generated by $X$ and a differential polynomial $A$, which can be obtained through a rational solution of a partial differential equation that is given explicitly in this paper. When the order is infinite, the expansion $\Lambda$ is just the unit ideal. In additional, if, and only if, the order of $\mathscr{X}$ is $0$, $1$, or $2$, the polynomial differential equation associating with $\mathscr{X}$ has Liouvillian first integrals. Examples and connections with Godbilon-Vey sequences are also discussed.

No abstract uploaded!