No Arabic abstract
On the idea of mapped WENO-JS scheme, properties of mapping methods are analyzed, uncertainties in mapping development are investigated, and new rational mappings are proposed. Based on our former understandings, i.e. mapping at endpoints {0, 1} tending to identity mapping, an integrated Cm,n condition is summarized for function development. Uncertainties, i.e., whether the mapping at endpoints would make mapped scheme behave like WENO or ENO, whether piecewise implementation would entail numerical instability, and whether WENO3-JS could preserve the third-order at first-order critical points by mapping, are analyzed and clarified. A new piecewise rational mapping with sufficient regulation capability is developed afterwards, where the flatness of mapping around the linear weights and its endpoint convergence toward identity mapping can be coordinated explicitly and simultaneously. Hence, the increase of resolution and preservation of stability can be balanced. Especially, concrete mappings are determined for WENO3,5,7-JS. Numerical cases are tested for the new mapped WENO-JS, which regards numerical stability including that in long time computation, resolution and robustness. In purpose of comparison, some recent mappings such as IM by [App. Math. Comput. 232, 2014:453-468], RM by [J. Sci. Comput. 67, 2016:540-580] and AIM by [J. Comput. Phys. 381, 2019:162-188] are chosen; in addition, some recent WENO-Z type scheme are selected also. Proposed new schemes can preserve optimal orders at corresponding critical points, achieve numerical stability and indicate overall comparative advantages regarding accuracy, resolution and robustness.
In this paper, high order semi-implicit well-balanced and asymptotic preserving finite difference WENO schemes are proposed for the shallow water equations with a non-flat bottom topography. We consider the Froude number ranging from O(1) to 0, which in the zero Froude limit becomes the lake equations for balanced flow without gravity waves. We apply a well-balanced finite difference WENO reconstruction, coupled with a stiffly accurate implicit-explicit (IMEX) Runge-Kutta time discretization. The resulting semi-implicit scheme can be shown to be well-balanced, asymptotic preserving (AP) and asymptotically accurate (AA) at the same time. Both one- and two-dimensional numerical results are provided to demonstrate the high order accuracy, AP property and good performance of the proposed methods in capturing small perturbations of steady state solutions.
In this paper, a fifth-order Hermite weighted essentially non-oscillatory (HWENO) scheme with artificial linear weights is proposed for one and two dimensional hyperbolic conservation laws, where the zeroth-order and the first-order moments are used in the spatial reconstruction. We construct the HWENO methodology using a nonlinear convex combination of a high degree polynomial with several low degree polynomials, and the associated linear weights can be any artificial positive numbers with only requirement that their summation equals one. The one advantage of the HWENO scheme is its simplicity and easy extension to multi-dimension in engineering applications for we can use any artificial linear weights which are independent on geometry of mesh. The another advantage is its higher order numerical accuracy using less candidate stencils for two dimensional problems. In addition, the HWENO scheme still keeps the compactness as only immediate neighbor information is needed in the reconstruction and has high efficiency for directly using linear approximation in the smooth regions. In order to avoid nonphysical oscillations nearby strong shocks or contact discontinuities, we adopt the thought of limiter for discontinuous Galerkin method to control the spurious oscillations. Some benchmark numerical tests are performed to demonstrate the capability of the proposed scheme.
Mode-based model-reduction is used to reduce the degrees of freedom of high dimensional systems, often by describing the system state by a linear combination of spatial modes. Transport dominated phenomena, ubiquitous in technical and scientific applications, often require a large number of linear modes to obtain a small representation error. This difficulty, even for the most simple transports, originates from the inappropriateness of the decomposition structure in time dependent amplitudes of purely spatial modes. In this article an approach is discussed, which decomposes a flow field into several fields of co-moving frames, where each one can be approximated by a few modes. The method of decomposition is formulated as an optimization problem. Different singular-value-based objective functions are discussed and connected to former formulations. A boundary treatment is provided. The decomposition is applied to generic cases and to a technically relevant flow configuration of combustion physics.
We solve the Boltzmann equation whose collision term is modeled by the hybridization of the binary collision and the BGK approximation. The parameter controlling the ratio of these two collision mechanisms is selected adaptively on every grid cell at every time step. This self-adaptation is based on a heuristic error indicator describing the difference between the model collision term and the original binary collision term. The indicator is derived by controlling the quadratic terms in the modeling error with linear operators. Our numerical experiments show that such error indicator is effective and computationally affordable.
The blood flow model maintains the steady state solutions, in which the flux gradients are non-zero but exactly balanced by the source term. In this paper, we design high order finite difference weighted non-oscillatory (WENO) schemes to this model with such well-balanced property and at the same time keeping genuine high order accuracy. Rigorous theoretical analysis as well as extensive numerical results all indicate that the resulting schemes verify high order accuracy, maintain the well-balanced property, and keep good resolution for smooth and discontinuous solutions.