An official website of the United States government
We propose an operator-splitting scheme to approximate scalar conservation equations with stiff source terms having multiple (at least two) stable equilibrium points. The scheme com- bines a (reaction-free) transport substep followed by a (transport-free) reaction substep. The transport substep is approximated using the forward Euler method with continuous finite elements and graph viscosity. The reaction substep is approximated using an exponential integrator. The crucial idea of the paper is to use a mesh-dependent cutoff of the reaction time-scale in the reaction substep. We establish a bound on the entropy residual motivating the design of the scheme. We show that the proposed scheme is invariant-domain preserv- ing under the same CFL restriction on the time step as in the nonreactive case. Numerical experiments in one and two space dimensions using linear, convex, and nonconvex fluxes with smooth and nonsmooth initial data in various regimes show that the proposed scheme is asymptotic preserving.
more »
« less
A Coupled Method Combining Crouzeix-Raviart Nonconforming and Node Conforming Finite Element Spaces for Boit Consolidation Model
A mixed finite element method is presented for the Biot consolidation problem in poroelasticity. More precisely, the displacement is approximated by using the Crouzeix-Raviart nonconforming finite elements, while the fluid pressure is approximated by using the node conforming finite elements. The well-posedness of the fully discrete scheme is established, and a corresponding priori error estimate with optimal order in the energy norm is also derived. Numerical experiments are provided to validate the theoretical results.
more »
« less
- PAR ID:
- 10536777
- Publisher / Repository:
- Global Science Press
- Date Published:
- Journal Name:
- Journal of Computational Mathematics
- Volume:
- 42
- Issue:
- 4
- ISSN:
- 0254-9409
- Page Range / eLocation ID:
- 911 to 931
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
An H(div)-conforming finite element method for the Biot’s consolidation mo- del is developed, with displacements and fluid velocity approximated by elements from BDM_k space. The use of H(div)-conforming elements for flow variables ensures the local mass conservation. In the H(div)-conforming approximation of displacement, the tan- gential components are discretised in the interior penalty discontinuous Galerkin frame- work,and the normal components across the element interfaces are continuous. Having introduced a spatial discretisation, we develop a semi-discrete scheme and a fully dis- crete scheme,prove their unique solvability and establish optimal error estimates for each variable.more » « less
-
Several well-known open questions (such as: are all groups sofic/hyperlinear?) have a common form: can all groups be approximated by asymptotic homomorphisms into the symmetric groups $$\text{Sym}(n)$$ (in the sofic case) or the finite-dimensional unitary groups $$\text{U}(n)$$ (in the hyperlinear case)? In the case of $$\text{U}(n)$$ , the question can be asked with respect to different metrics and norms. This paper answers, for the first time, one of these versions, showing that there exist finitely presented groups which are not approximated by $$\text{U}(n)$$ with respect to the Frobenius norm $$\Vert T\Vert _{\text{Frob}}=\sqrt{\sum _{i,j=1}^{n}|T_{ij}|^{2}},T=[T_{ij}]_{i,j=1}^{n}\in \text{M}_{n}(\mathbb{C})$$ . Our strategy is to show that some higher dimensional cohomology vanishing phenomena implies stability , that is, every Frobenius-approximate homomorphism into finite-dimensional unitary groups is close to an actual homomorphism. This is combined with existence results of certain nonresidually finite central extensions of lattices in some simple $$p$$ -adic Lie groups. These groups act on high-rank Bruhat–Tits buildings and satisfy the needed vanishing cohomology phenomenon and are thus stable and not Frobenius-approximated.more » « less
-
Finite element de Rham complexes and finite element Stokes complexes with varying degrees of smoothness in three dimensions are systematically constructed in this paper. Smooth scalar finite elements in three dimensions are derived through a non-overlapping decomposition of the simplicial lattice. H(div)-conforming finite elements and H(curl)-conforming finite elements with varying degrees of smoothness are devised based on these smooth scalar finite elements. The finite element de Rham complexes with corresponding smoothness and commutative diagrams are induced by these elements. The div stability of the H(div)-conforming finite elements is established, and the exactness of these finite element complexes is proven.more » « less
-
The convergence property of a stochastic algorithm for the self-consistent field (SCF) calculations of electron structures is studied. The algorithm is formulated by rewriting the electronic charges as a trace/diagonal of a matrix function, which is subsequently expressed as a statistical average. The function is further approximated by using a Krylov subspace approximation. As a result, each SCF iteration only samples one random vector without having to compute all the orbitals. We consider the common practice of SCF iterations with damping and mixing. We prove that the iterates from a general linear mixing scheme converge in a probabilistic sense when the stochastic error has a second finite moment.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.4208/jcm.2212-m2021-0231
