An official website of the United States government
We present the lowest-order hybridizable discontinuous Galerkin schemes with numerical integration (quadrature), denoted as HDG-P0 for the reaction-diffusion equation and the generalized Stokes equations on conforming simplicial meshes in two- and three-dimensions. Here by lowest order, we mean that the (hybrid) finite element space for the global HDG facet degrees of freedom (DOFs) is the space of piecewise constants on the mesh skeleton. A discontinuous piecewise linear space is used for the approximation of the local primal unknowns. We give the optimal a priori error analysis of the proposed HDG-P0 schemes, which hasn’t appeared in the literature yet for HDG discretizations as far as numerical integration is concerned. Moreover, we propose optimal geometric multigrid preconditioners for the statically condensed HDG-P0 linear systems on conforming simplicial meshes. In both cases, we first establish the equivalence of the statically condensed HDG system with a (slightly modified) nonconforming Crouzeix–Raviart (CR) discretization, where the global (piecewise-constant) HDG finite element space on the mesh skeleton has a natural one-to-one correspondence to the nonconforming CR (piecewise-linear) finite element space that live on the whole mesh. This equivalence then allows us to use the well-established nonconforming geometry multigrid theory to precondition the condensed HDG system. Numerical results in two- and three-dimensions are presented to verify our theoretical findings.
more »
« less
Optimal error estimates of ultra-weak discontinuous Galerkin methods with generalized numerical fluxes for multi-dimensional convection-diffusion and biharmonic equations
In this paper, we study ultra-weak discontinuous Galerkin methods with generalized numerical fluxes for multi-dimensional high order partial differential equations on both unstructured simplex and Cartesian meshes. The equations we consider as examples are the nonlinear convection-diffusion equation and the biharmonic equation. Optimal error estimates are obtained for both equations under certain conditions, and the key step is to carefully design global projections to eliminate numerical errors on the cell interface terms of ultra-weak schemes on general dimensions. The well-posedness and approximation capability of these global projections are obtained for arbitrary order polynomial space based on a wide class of generalized numerical fluxes on regular meshes. These projections can serve as general analytical tools to be naturally applied to a wide class of high order equations. Numerical experiments are conducted to demonstrate these theoretical results.
more »
« less
- PAR ID:
- 10537895
- Publisher / Repository:
- American Mathematical Society
- Date Published:
- Journal Name:
- Mathematics of Computation
- Volume:
- 93
- Issue:
- 349
- ISSN:
- 0025-5718
- Page Range / eLocation ID:
- 2135 to 2183
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Simulation of flow and transport in petroleum reservoirs involves solving coupled systems of advection-diffusion-reaction equations with nonlinear flux functions, diffusion coefficients, and reactions/wells. It is important to develop numerical schemes that can approximate all three processes at once, and to high order, so that the physics can be well resolved. In this paper, we propose an approach based on high order, finite volume, implicit, Weighted Essentially NonOscillatory (iWENO) schemes. The resulting schemes are locally mass conservative and, being implicit, suited to systems of advection-diffusion-reaction equations. Moreover, our approach gives unconditionally L-stable schemes for smooth solutions to the linear advection-diffusion-reaction equation in the sense of a von Neumann stability analysis. To illustrate our approach, we develop a third order iWENO scheme for the saturation equation of two-phase flow in porous media in two space dimensions. The keys to high order accuracy are to use WENO reconstruction in space (which handles shocks and steep fronts) combined with a two-stage Radau-IIA Runge-Kutta time integrator. The saturation is approximated by its averages over the mesh elements at the current time level and at two future time levels; therefore, the scheme uses two unknowns per grid block per variable, independent of the spatial dimension. This makes the scheme fairly computationally efficient, both because reconstructions make use of local information that can fit in cache memory, and because the global system has about as small a number of degrees of freedom as possible. The scheme is relatively simple to implement, high order accurate, maintains local mass conservation, applies to general computational meshes, and appears to be robust. Preliminary computational tests show the potential of the scheme to handle advection-diffusion-reaction processes on meshes of quadrilateral gridblocks, and to do so to high order accuracy using relatively long time steps. The new scheme can be viewed as a generalization of standard cell-centered finite volume (or finite difference) methods. It achieves high order in both space and time, and it incorporates WENO slope limiting.more » « less
-
Abstract In this paper we study the biharmonic equation with Navier boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the fourth-order problem as a system of Poisson equations. Our method differs from the naive mixed method that leads to two Poisson problems but only applies to convex domains; our decomposition involves a third Poisson equation to confine the solution in the correct function space, and therefore can be used in both convex and nonconvex domains. A $C^0$ finite element algorithm is in turn proposed to solve the resulting system. In addition, we derive optimal error estimates for the numerical solution on both quasi-uniform meshes and graded meshes. Numerical test results are presented to justify the theoretical findings.more » « less
-
The generalized nonlinear Schr\"odinger equation with full dispersion (FDNLS) is considered in the semiclassical regime. The Whitham modulation equations are obtained for the FDNLS equation with general linear dispersion and a generalized, local nonlinearity. Assuming the existence of a four-parameter family of two-phase solutions, a multiple-scales approach yields a system of four independent, first-order, quasi-linear conservation laws of hydrodynamic type that correspond to the slow evolution of the two wavenumbers, mass, and momentum of modulated periodic traveling waves. The modulation equations are further analyzed in the dispersionless and weakly nonlinear regimes. The ill-posedness of the dispersionless equations corresponds to the classical criterion for modulational instability (MI). For modulations of linear waves, ill-posedness coincides with the generalized MI criterion, recently identified by Amiranashvili and Tobisch [New J. Phys., 21 (2019), 033029]. A new instability index is identified by the transition from real to complex characteristics for the weakly nonlinear modulation equations. This instability is associated with long wavelength modulations of nonlinear two-phase wavetrains and can exist even when the corresponding one-phase wavetrain is stable according to the generalized MI criterion. Another interpretation is that while infinitesimal perturbations of a periodic wave may not grow, small but finite amplitude perturbations may grow, hence this index identifies a nonlinear instability mechanism for one-phase waves. Classifications of instability indices for multiple FDNLS equations with higher-order dispersion, including applications to finite-depth water waves and the discrete NLS equation, are presented and compared with direct numerical simulations.more » « less
-
null (Ed.)Pressure Poisson equation (PPE) reformulations of the incompressible Navier–Stokes equations (NSE) replace the incompressibility constraint by a Poisson equation for the pressure and a suitable choice of boundary conditions. This yields a time-evolution equation for the velocity field only, with the pressure gradient acting as a nonlocal operator. Thus, numerical methods based on PPE reformulations are representatives of a class of methods that have no principal limitations in achieving high order. In this paper, it is studied to what extent high-order methods for the NSE can be obtained from a specific PPE reformulation with electric boundary conditions (EBC). To that end, implicit–explicit (IMEX) time-stepping is used to decouple the pressure solve from the velocity update, while avoiding a parabolic time-step restriction; and mixed finite elements are used in space, to capture the structure imposed by the EBC. Via numerical examples, it is demonstrated that the methodology can yield at least third order accuracy in space and time.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1090/mcom/3927
