Search

This paper proposes and analyzes a novel efficient high-order finite volume method for the ideal magnetohydrodynamics (MHD). As a distinctive feature, the method simultaneously preserves two critical physical constraints: a discretely divergence-free (DDF) constraint on the magnetic field and the positivity-preserving (PP) property, which ensures the positivity of density, pressure, and internal energy. To enforce the DDF condition in each cell, we design a new discrete projection approach that projects the reconstructed point values at the cell interface into a DDF space, without using any approximation polynomials. This projection method is highly efficient, easy to implement, and particularly suitable for the high-order finite volume methods that return only the point values (no explicit approximation polynomials) in the reconstruction. Moreover, we also develop a new finite volume framework for constructing provably PP schemes for the ideal MHD system. The framework comprises the discrete projection technique, a suitable approximation to the Godunov--Powell source terms, and a simple PP limiter. We provide rigorous analysis of the PP property of the proposed finite volume method, demonstrating that the DDF condition and the proper approximation to the source terms eliminate the impact of magnetic divergence terms on the PP property. The analysis is challenging due to the internal energy function's nonlinearity and the intricate relationship between the DDF and PP properties. To address these challenges, we adopt the recently developed geometric quasilinearization approach [K. Wu and C.-W. Shu, SIAM Rev., 65 (2023), pp. 1031--1073], which transforms a nonlinear constraint into a family of linear constraints. Finally, we validate the effectiveness of the proposed method through several benchmark and demanding numerical examples. The results demonstrate that the proposed method is robust, accurate, and highly effective, confirming the significance of the proposed DDF projection and PP techniques. [ABSTRACT FROM AUTHOR]

In the numerical simulation of ideal magnetohydrodynamics (MHD), keeping the pressure and density always positive is essential for both physical considerations and numerical stability. This is however a challenging task, due to the underlying relation between such a positivity-preserving (PP) property and the magnetic divergence-free (DF) constraint as well as the strong nonlinearity of the MHD equations. In this paper, we present the first rigorous PP analysis of the central discontinuous Galerkin (CDG) methods and construct arbitrarily high-order provably PP CDG schemes for ideal MHD. By the recently developed geometric quasilinearization (GQL) approach, our analysis reveals that the PP property of standard CDG methods is closely related to a discrete magnetic DF condition, whose form was unknown prior to our analysis and differs from that for the noncentral discontinuous Galerkin (DG) and finite volume methods in [K. Wu, SIAM J. Numer. Anal., 56 (2018), pp. 2124-2147]. The discovery of this relation lays the foundation for the design of our PP CDG schemes. In the one-dimensional (1D) case, the discrete DF condition is naturally satisfied, and we rigorously prove that the standard CDG method is PP under a condition that can be enforced easily with an existing PP limiter. However, in the multidimensional cases, the corresponding discrete DF condition is highly nontrivial, yet critical, and we analytically prove that the standard CDG method, even with the PP limiter, is not PP in general, as it generally fails to meet the discrete DF condition. We address this issue by carefully analyzing the structure of the discrete divergence terms and then constructing new locally DF CDG schemes for Godunov's modified MHD equations with an additional source term. The key point is to find out the suitable discretization of the source term such that it exactly cancels out all the terms in the discovered discrete DF condition. Based on the GQL approach, we prove in theory the PP property of the new multidimensional CDG schemes under a CFL condition. The robustness and accuracy of the proposed PP CDG schemes are further validated by several demanding 1D, two-dimensional, and three-dimensional numerical MHD examples, including the high-speed jets and blast problems with very low plasma beta. [ABSTRACT FROM AUTHOR]

This paper is the first in a series of two articles, aiming to prove the convergence of the HX preconditioner originally proposed by Hiptmair and Xu [SIAM J. Numer. Anal., 45 (2007), pp. 2483-2509] for Maxwell's equations with jump coefficients. In this paper, we establish several extensions of the discrete regular decomposition for edge finite element functions defined in threedimensional domains. The functions defined by the new discrete regular decompositions can inherit zero degrees of freedom of the considered edge finite element function on some faces and edges of polyhedral domains as well as of some non-Lipschitz domains and possess nearly optimal stability with only a logarithmic factor. [ABSTRACT FROM AUTHOR]

We correct assumption (4.6) in Theorem 4.1 in [F. Gozzi and F. Masiero, SIAM J. Control Optim., 55 (2017), pp. 2981-3012], and we consequently change the proof of Theorem 4.1. We correct the proof of Lemma 5.4 in [F. Gozzi and F. Masiero, SIAM J. Control Optim., 55 (2017), pp. 2981-3012]. Moreover, throughout the paper, starting from formulas (3.13) and (3.14), we notice that the current cost of the control problem cannot depend on the state for the results of the paper to work. This change in the current cost must be made in the paper [F. Gozzi and F. Masiero, SIAM J. Control Optim., 55 (2017), pp. 3013-3038], too. [ABSTRACT FROM AUTHOR]

In this paper we consider a numerical homogenization technique for curl-curl problems that is based on the framework of the localized orthogonal decomposition and which was proposed in [D. Gallistl, P. Henning, and B. Verfürth, SIAM J. Numer. Anal., 56 (2018), pp. 1570-1596] for problems with essential boundary conditions. The findings of the aforementioned work establish quantitative homogenization results for the time-harmonic Maxwell's equations that hold beyond assumptions of periodicity; however, a practical realization of the approach was left open. In this paper, we transfer the findings from essential boundary conditions to natural boundary conditions, and we demonstrate that the approach yields a computable numerical method. We also investigate how boundary values of the source term can affect the computational complexity and accuracy. Our findings will be supported by various numerical experiments, both in 2D and 3D. [ABSTRACT FROM AUTHOR]

In this paper, we present a new regularization scheme for the linearized Poisson{ Boltzmann equation (PBE) which models the electrostatic potential of biomolecules in a solvent. This scheme is based on the splitting of the target potential into the short- and long-range components localized in the molecular region by using the range-separated (RS) tensor format [P. Benner, V. Khoromskaia, and B. N. Khoromskij, SIAM J. Comput., 2 (2018), pp. A1034{A1062] for representation of the discretized multiparticle Dirac delta [B. N. Khoromskij, J. Comput. Phys., 401 (2020), 108998] constituting the highly singular right-hand side in the PBE. From the computational point of view our regularization approach requires only the modification of the right-hand side in the PBE so that it can be implemented within any open-source grid-based software package for solving PBE that already includes some FEM/FDM disretization scheme for elliptic PDE and solver for the arising linear system of equations. The main computational benefits are twofold. First, one applies the chosen PBE solver only for the smooth long-range (regularized) part of the collective potential with the regular right-hand side represented by a low-rank RS tensor with a controllable precision. Thus, we eliminate the numerical treatment of the singularities in the right-hand side and do not change the interface and boundary conditions. And second, the elliptic PDE need not be solved for the singular part in the right-hand side at all, since the short-range part of the target potential of the biomolecule is precomputed independently on a computational grid by simple one-dimensional tensor operations. The total potential is then obtained by adding the numerical solution of the PBE for the smooth long-range part to the directly precomputed tensor representation for the short-range contribution. Numerical tests illustrate that the new regularization scheme, implemented by a simple modification of the right-hand side in the chosen PBE solver, improves the accuracy of the approximate solution on rather coarse grids. The scheme also demonstrates good convergence behavior on a sequence of refined grids. [ABSTRACT FROM AUTHOR]

The aim of this paper is to develop and analyze high-order time stepping schemes for approximately solving semilinear subdiffusion equations. We apply the convolution quadrature generated by k-step backward differentiation formula (BDFk) to discretize the time-fractional derivative with order α ∈ (0, 1) and modify the starting steps in order to achieve optimal convergence rate. This method has already been well-studied for the linear fractional evolution equations in Jin, Li, and Zhou [SIAM J. Sci. Comput., 39 (2017), pp. A3129-A3152], while the numerical analysis for the nonlinear problem is still missing in the literature. By splitting the nonlinear potential term into an irregular linear part and a smoother nonlinear part and using the generating function technique, we prove that the convergence order of the corrected BDFk scheme is O(τmin(k,1+2α-ε)) without imposing further assumption on the regularity of the solution. Numerical examples are provided to support our theoretical results. [ABSTRACT FROM AUTHOR]

This erratum clarifies that equations (4.21) and (4.23) in Böttcher and Porter [SIAM J. Appl. Math., 81 (2021), pp. 2704-2724] assume that the underlying evolution operator has an orthonormal basis. The results and conclusions of the paper are unaffected by this clarification. [ABSTRACT FROM AUTHOR]

In [L. Li, M. Ni, and W. Zheng, SIAM J. Sci. Comput., 41 (2019), pp. B796-B815] a charge-conservative finite element method is proposed for solving inductionless and incompressible magnetohydrodynamic (MHD) equations. The purpose of this paper is to propose a robust solver for the discrete problem. Using the framework of field-of-values-equivalence, we first study the preconditioned Krylov space method for the continuous problem in the setting of Hilbert spaces. The algebraic preconditioner for the discrete problem is then obtained by representing the preconditioner for the continuous problem in finite element spaces. By three numerical examples, the optimality of the solver to the number of unknowns is demonstrated for both stationary and time-dependent MHD problems. [ABSTRACT FROM AUTHOR]

In this paper we study the existence of sufficiently regular representations of Hamilton{Jacobi equations in optimal control theory with the compact control set. We introduce a new method to construct representations for a wide class of Hamiltonians, wider than was achieved before. Our result is proved by means of these conditions on the Hamiltonian that are necessary for the existence of a representation. In particular, we solve an open problem of Rampazzo [SIAM J. Control Optim., 44 (2005), pp. 867{884]. We apply the obtained results to reduce a variational problem to an optimal control problem. [ABSTRACT FROM AUTHOR]