BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Kaibo Hu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
The report is based on an extended abstract for the MFO workshop "Hilbert Complexes: Analysis, Applications, and Discretizations", held at Oberwolfach during 19-25 June 2022. The aim is to provide an overview of some aspects of discretization of Hilbert complexes with an emphasis on conforming finite elements.
We investigate discretization of $H(\mathrm{curl})$ and $H(\mathrm{div})$ in two and three space dimensions by partially discontinuous nodal finite elements, i.e., vector-valued Lagrange finite elements with discontinuity in certain directions. These spaces can be implemented as a combination of continuous and discontinuous Lagrange elements and fit in de~Rham complexes. We construct well-conditioned nodal bases.
We develop structure-preserving finite element methods for the incompressible, resistive Hall magnetohydrodynamics (MHD) equations. These equations incorporate the Hall current term in Ohm's law and provide a more appropriate description of fully ionized plasmas than the standard MHD equations on length scales close to or smaller than the ion skin depth. We introduce a stationary discrete variational formulation of Hall MHD that enforces the magnetic Gauss's law exactly (up to solver tolerances) and prove the well-posedness and convergence of a Picard linearization. For the transient problem, we present time discretizations that preserve the energy and magnetic and hybrid helicity precisely in the ideal limit for two types of boundary conditions. Additionally, we present an augmented Lagrangian preconditioning technique for both the stationary and transient cases. We confirm our findings with several numerical experiments.
We consider spline functions over simplicial meshes in $\RR^n$. We assume that the spline pieces join together with some finite order of smoothness but the pieces themselves are infinitely smooth. Such splines can have extra orders of smoothness at a vertex, a property known as \emph{supersmoothness}, which plays a role in the construction of multivariate splines and in the finite element method. In this paper we characterize supersmoothness in terms of the degeneracy of spaces of polynomial splines over the cell of simplices sharing the vertex, and use it to determine the maximal order of supersmoothness of various cell configurations.
We construct conforming finite element elasticity complexes on the Alfeld splits of tetrahedra. The complex consists of vector fields and symmetric tensor fields, interlinked via the linearized deformation operator, the linearized curvature operator, and the divergence operator, respectively. The construction is based on an algebraic machinery that derives the elasticity complex from de~Rham complexes, and smoother finite element differential forms.
We develop a theory of Finite Element Systems, for the purpose of discretizing sections of vector bundles, in particular those arizing in the theory of elasticity. In the presence of curvature we prove a discrete Bianchi identity. In the flat case we prove a de Rham theorem on cohomology groups. We check that some known mixed finite elements for the stress-displacement formulation of elasticity fit our framework. We also define, in dimension two, the first conforming finite element spaces of metrics with good linearized curvature, corresponding to strain tensors with Saint-Venant compatibility conditions. Cochains with coefficients in rigid motions are given a key role in relating continuous and discrete elasticity complexes.
The purpose of this paper is to discuss representations of high order $C^0$ finite element spaces on simplicial meshes in any dimension. When computing with high order piecewise polynomials the conditioning of the basis is likely to be important. The main result of this paper is a construction of representations by frames such that the associated $L^2$ condition number is bounded independently of the polynomial degree. To our knowledge, such a representation has not been presented earlier. The main tools we will use for the construction is the bubble transform, introduced previously in [Falk and Winther, Found Comput Math (2016) 16: 297], and properties of Jacobi polynomials on simplexes in higher dimensions. We also include a brief discussion of preconditioned iterative methods for the finite element systems in the setting of representations by frames.
In this paper, we develop two classes of robust preconditioners for the structure-preserving discretization of the incompressible magnetohydrodynamics (MHD) system. By studying the well-posedness of the discrete system, we design block preconditioners for them and carry out rigorous analysis on their performance. We prove that such preconditioners are robust with respect to most physical and discretization parameters. In our proof, we improve the existing estimates of the block triangular preconditioners for saddle point problems by removing the scaling parameters, which are usually difficult to choose in practice. This new technique is not only applicable to the MHD system, but also to other problems. Moreover, we prove that Krylov iterative methods with our preconditioners preserve the divergence-free condition exactly, which complements the structure-preserving discretization. Another feature is that we can directly generalize this technique to other discretizations of the MHD system. We also present preliminary numerical results to support the theoretical results and demonstrate the robustness of the proposed preconditioners.
This paper is devoted to the design and analysis of some structure-preserving finite element schemes for the magnetohydrodynamics (MHD) system. The main feature of the method is that it naturally preserves the important Gauss law, namely $\nabla\cdot\boldsymbol{B}=0$. In contrast to most existing approaches that eliminate the electrical field variable $\boldsymbol{E}$ and give a direct discretization of the magnetic field, our new approach discretizes the electric field $\boldsymbol{E}$ by Nédélec type edge elements for $H(\mathrm{curl})$, while the magnetic field $\boldsymbol{B}$ by Raviart-Thomas type face elements for $H(\mathrm{div})$. As a result, the divergence-free condition on the magnetic field holds exactly on the discrete level. For this new finite element method, an energy stability estimate can be naturally established in an analogous way as in the continuous case. Furthermore, well-posedness is rigorously established in the paper for both the Picard and Newton linearization of the fully nonlinear systems by using the Brezzi theory for both the continuous and discrete cases. This well-posedness naturally leads to robust (and optimal) preconditioners for the linearized systems.