More Like this

1. Discrete elasticity exact sequences on Worsey–Farin splits

   https://doi.org/10.1051/m2an/2023084  

   Gong, Sining; Gopalakrishnan, Jay; Guzmán, Johnny; Neilan, Michael (November 2023, ESAIM: Mathematical Modelling and Numerical Analysis)

   We construct conforming finite element elasticity complexes on Worsey–Farin splits in three dimensions. Spaces for displacement, strain, stress, and the load are connected in the elasticity complex through the differential operators representing deformation, incompatibility, and divergence. For each of these component spaces, a corresponding finite element space on Worsey–Farin meshes is exhibited. Unisolvent degrees of freedom are developed for these finite elements, which also yields commuting (cochain) projections on smooth functions. A distinctive feature of the spaces in these complexes is the lack of extrinsic supersmoothness at subsimplices of the mesh. Notably, the complex yields the first (strongly) symmetric stress finite element with no vertex or edge degrees of freedom in three dimensions. Moreover, the lowest order stress space uses only piecewise linear functions which is the lowest feasible polynomial degree for the stress space. 

   more » « less
2. A divergence-free and $H\left(div\right)$ -conforming embedded-hybridized DG method for the incompressible resistive MHD equations

   https://doi.org/10.1016/j.cma.2024.117415  

   Chen, Jau-Uei; Horváth, Tamás L; Bui-Thanh, Tan (December 2024, Computer Methods in Applied Mechanics and Engineering)

   We present a divergence-free and $$\Hsp\LRp{div}$$-conforming hybridized discontinuous Galerkin (HDG) method and a computationally efficient variant called embedded-HDG (E-HDG) for solving stationary incompressible viso-resistive magnetohydrodynamic (MHD) equations. The proposed E-HDG approach uses continuous facet unknowns for the vector-valued solutions (velocity and magnetic fields) while it uses discontinuous facet unknowns for the scalar variable (pressure and magnetic pressure). This choice of function spaces makes E-HDG computationally far more advantageous, due to the much smaller number of degrees of freedom, compared to the HDG counterpart. The benefit is even more significant for three-dimensional/high-order/fine mesh scenarios. On simplicial meshes, the proposed methods with a specific choice of approximation spaces are well-posed for linear(ized) MHD equations. For nonlinear MHD problems, we present a simple approach exploiting the proposed linear discretizations by using a Picard iteration. The beauty of this approach is that the divergence-free and $$\Hsp\LRp{div}$$-conforming properties of the velocity and magnetic fields are automatically carried over for nonlinear MHD equations. We study the accuracy and convergence of our E-HDG method for both linear and nonlinear MHD cases through various numerical experiments, including two- and three-dimensional problems with smooth and singular solutions. The numerical examples show that the proposed methods are pressure robust, and the divergence of the resulting velocity and magnetic fields is machine zero for both smooth and singular problems. 

   more » « less
3. Uniform block-diagonal preconditioners for divergence-conforming HDG Methods for the generalized Stokes equations and the linear elasticity equations

   https://doi.org/10.1093/imanum/drac021  

   Fu, Guosheng; Kuang, Wenzheng (June 2022, IMA Journal of Numerical Analysis)

   Abstract We propose a uniform block-diagonal preconditioner for condensed $$H$$(div)-conforming hybridizable discontinuous Galerkin schemes for parameter-dependent saddle point problems, including the generalized Stokes equations and the linear elasticity equations. An optimal preconditioner is obtained for the stiffness matrix on the global velocity/displacement space via the auxiliary space preconditioning technique (Xu (1994) The Auxiliary Space Method and Optimal Multigrid Preconditioning Techniques for Unstructured Grids, vol. 56. International GAMM-Workshop on Multi-level Methods (Meisdorf), pp. 215–235). A spectrally equivalent approximation to the Schur complement on the element-wise constant pressure space is also constructed, and an explicit computable exact inverse is obtained via the Woodbury matrix identity. Finally, the numerical results verify the robustness of our proposed preconditioner with respect to model parameters and mesh size. 

   more » « less
4. Construction of Scalar and Vector Finite Element Families on Polygonal and Polyhedral Meshes

   https://doi.org/10.1515/cmam-2016-0019  

   Gillette, Andrew; Rand, Alexander; Bajaj, Chandrajit (January 2016, Computational Methods in Applied Mathematics)

   Abstract We combine theoretical results from polytope domain meshing, generalized barycentric coordinates, and finite element exterior calculus to construct scalar- and vector-valued basis functions for conforming finite element methods on generic convex polytope meshes in dimensions 2 and 3. Our construction recovers well-known bases for the lowest order Nédélec, Raviart–Thomas, and Brezzi–Douglas–Marini elements on simplicial meshes and generalizes the notion of Whitney forms to non-simplicial convex polygons and polyhedra. We show that our basis functions lie in the correct function space with regards to global continuity and that they reproduce the requisite polynomial differential forms described by finite element exterior calculus. We present a method to count the number of basis functions required to ensure these two key properties. 

   more » « less
5. The VEM for time-harmonic Maxwell equations in inhomogeneous media with Lipschitz interface

   https://doi.org/10.1142/S0218202525500289  

   Chen, Chunyu; Guo, Ruchi; Wei, Huayi (May 2025, Mathematical Models and Methods in Applied Sciences)

   It has been extensively studied in the literature that solving Maxwell equations is very sensitive to mesh structures, space conformity and solution regularity. Roughly speaking, for almost all the methods in the literature, optimal convergence for low-regularity solutions heavily relies on conforming spaces and highly regular simplicial meshes. This can be a significant limitation for many popular methods based on broken spaces and non-conforming or polytopal meshes often used for inhomogeneous media, as the discontinuity of electromagnetic parameters can lead to quite low regularity of solutions near media interfaces. This very issue can be potentially worsened by geometric singularities, making those methods particularly challenging to apply. In this paper, we present a lowest-order virtual element method for solving an indefinite time-harmonic Maxwell equation in 2D inhomogeneous media with quite arbitrary polytopal meshes, and the media interface is allowed to have geometric singularity to cause low regularity. We employ the “virtual mesh” technique originally invented in [S. Cao, L. Chen and R. Guo, A virtual finite element method for two-dimensional Maxwell interface problems with a background unfitted mesh, Math. Models Methods Appl. Sci. 31 (2021) 2907–2936] for error analysis. This work admits three key novelties: (i) the proposed method is theoretically guaranteed to achieve robust optimal convergence for solutions with merely [Formula: see text] regularity, [Formula: see text]; (ii) the polytopal element shape can be highly anisotropic and shrinking, and an explicit formula is established to describe the relationship between the shape regularity and solution regularity; (iii) we show that the stabilization term is needed to produce optimal convergent solutions for indefinite problems. Extensive numerical experiments will be given to demonstrate the effectiveness of the proposed method. 

   more » « less