No abstract uploaded!

We present a new stabilization technique for multiscale convection diffusion problems. Stabilization for these problems has been a challenging task, especially for the case with high Peclet numbers. Our method is based on a constraint energy minimization idea and the discontinuous Petrov-Galerkin formulation. In particular, the test functions are constructed by minimizing an appropriate energy subject to certain orthogonality conditions, and are related to the trial space. The resulting test functions have a localization property, and can therefore be computed locally. We will prove the stability, and present several numerical results. Our numerical results confirm that our test space gives a good stability, in the sense that the solution error is close to the best approximation error.

To easily generalize the maximum-principle-satisfying schemes for scalar conservation laws to convection diffusion equations, we propose a non-conventional high order finite volume weighted essentially non-oscillatory (WENO) scheme which can be proven maximum-principle-satisfying. Two-dimensional extensions are straightforward. We also show that the same idea can be used to construct high order schemes preserving the maximum principle for two-dimensional incompressible Navier-Stokes equations in the vorticity stream-function formulation. Numerical tests for the fifth order WENO schemes are reported.

A class of upwind flux splitting methods in the Euler equations of compressible flow is considered in this paper. Using the property that Euler flux F(U) is a homogeneous function of degree one in U, we reformulate the splitting fluxes with F^{+}=A^{+} U, F^{-}=A^{-} U, and the corresponding matrices are either symmetric or symmetrizable and keep only non-negative and non-positive eigenvalues. That leads to the conclusion that the first order schemes are positive in the sense of Lax-Liu [18], which implies that it is L^2-stable in some suitable sense. Moreover, the second order scheme is a stable perturbation of the first order scheme, so that the positivity of the second order schemes is also established, under a CFL-like condition. In addition, these splitting methods preserve the positivity of density and energy.

No abstract uploaded!