Search

We present a divergence-free and $H(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 $H(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.

In (J. Comput. Phys., 417, 109577, 2020) we introduced a space-time embedded-hybridizable discontinuous Galerkin method for the solution of the incompressible Navier-Stokes equations on time-dependent domains of which the motion of the domain is prescribed. This discretization is exactly mass conserving, locally momentum conserving, and energy-stable. In this manuscript we extend this discretization to fluid-rigid body interaction problems in which the motion of the fluid domain is not known a priori. To account for large rotational motion of the rigid body, we present a novel conforming space-time sliding mesh technique. We demonstrate the performance of the discretization on various numerical examples.

We introduce and analyze a space-time hybridized discontinuous Galerkin method for the evolutionary Navier--Stokes equations. Key features of the numerical scheme include point-wise mass conservation, energy stability, and pressure robustness. We prove that there exists a solution to the resulting nonlinear algebraic system in two and three spatial dimensions, and that this solution is unique in two spatial dimensions under a small data assumption. A priori error estimates are derived for the velocity in a mesh-dependent energy norm.

This paper presents a space-time embedded-hybridized discontinuous Galerkin (EHDG) method for the Navier--Stokes equations on moving domains. This method uses a different hybridization compared to the space-time hybridized discontinuous Galerkin (HDG) method we presented previously in (Int. J. Numer. Meth. Fluids 89: 519--532, 2019). In the space-time EHDG method the velocity trace unknown is continuous while the pressure trace unknown is discontinuous across facets. In the space-time HDG method, all trace unknowns are discontinuous across facets. Alternatively, we present also a space-time embedded discontinuous Galerkin (EDG) method in which all trace unknowns are continuous across facets. The advantage of continuous trace unknowns is that the formulation has fewer global degrees-of-freedom for a given mesh than when using discontinuous trace unknowns. Nevertheless, the discrete velocity field obtained by the space-time EHDG and EDG methods, like the space-time HDG method, is exactly divergence-free, even on moving domains. However, only the space-time EHDG and HDG methods result in divergence-conforming velocity fields. An immediate consequence of this is that the space-time EHDG and HDG discretizations of the conservative form of the Navier--Stokes equations are energy stable. The space-time EDG method, on the other hand, requires a skew-symmetric formulation of the momentum advection term to be energy-stable. Numerical examples will demonstrate the differences in solution obtained by the space-time EHDG, EDG, and HDG methods.

This paper presents the first analysis of a space--time hybridizable discontinuous Galerkin method for the advection--diffusion problem on time-dependent domains. The analysis is based on non-standard local trace and inverse inequalities that are anisotropic in the spatial and time steps. We prove well-posedness of the discrete problem and provide a priori error estimates in a mesh-dependent norm. Convergence theory is validated by a numerical example solving the advection--diffusion problem on a time-dependent domain for approximations of various polynomial degree.

We present a finite element method for the incompressible Navier--Stokes problem that is locally conservative, energy-stable and pressure-robust on time-dependent domains. To achieve this, the space--time formulation of the Navier--Stokes problem is considered. The space--time domain is partitioned into space--time slabs which in turn are partitioned into space--time simplices. A combined discontinuous Galerkin method across space--time slabs, and space--time hybridized discontinuous Galerkin method within a space--time slab, results in an approximate velocity field that is $H({\rm div})$-conforming and exactly divergence-free, even on time-dependent domains. Numerical examples demonstrate the convergence properties and performance of the method.