Related papers: A monolithic space-time temporal multirate finite …
For most finite element simulations, boundary-conforming meshes have significant advantages in terms of accuracy or efficiency. This is particularly true for complex domains. However, with increased complexity of the domain, generating a…
In ecological studies of pattern formation, models of the competitive-diffusion type are generally singularly perturbed, and the numerical approximation of such models is challenging. In this paper, we present finite element discretization…
The Stokes equations play an important role in the incompressible flow simulation. In this paper, a novel divergence-free parametric mixed finite element method is proposed for solving three-dimensional Stokes equations on domains with…
This paper proposes a parallel numerical algorithm to simulate the flow and the transport in a discrete fracture network taking into account the mass exchanges with the surrounding matrix. The discretization of the Darcy fluxes is based on…
Multiscale methods for second order elliptic equations based on non-overlapping domain decomposition schemes have great potential to take advantage of multi-core, state-of-the-art parallel computers. These methods typically involve solving…
We implement a stabilized finite element method for steady Darcy-Brinkman-Forchheimer model within the continuous Galerkin framework. The nonlinear fluid model is first linearized using a standard \textit{Newton's method. The sequence of…
We study the discretisation of a uniaxial (rank-one) reduction of the Oldroyd-B model for dilute polymer solutions, in which the conformation tensor is represented as $\sig = \vec b \otimes \vec b$. Building on structural analogies with…
Robust and accurate schemes are designed to simulate the coupling between subsurface and overland flows. The coupling conditions at the interface enforce the continuity of both the normal flux and the pressure. Richards' equation governing…
In many applications, one wants to model physical systems consisting of two different physical processes in two different domains that are coupled across a common interface. A crucial challenge is then that the solutions of the two…
We propose a boundary-corrected weak Galerkin mixed finite element method for solving elliptic interface problems in 2D domains with curved interfaces. The method is formulated on body-fitted polygonal meshes, where interface edges are…
We present a high-order hybridizable discontinuous Galerkin method for the numerical solution of time-dependent three-phase flow in heterogeneous porous media. The underlying algorithm is a semi-implicit operator splitting approach that…
In this paper we present a framework for solving two phase flow problems in porous media. The discretization is based on a Discontinuous Galerkin method and includes local grid adaptivity and local choice of polynomial degree. The method is…
We leverage the proximal Galerkin algorithm (Keith and Surowiec, Foundations of Computational Mathematics, 2024, DOI: 10.1007/s10208-024-09681-8), a recently introduced mesh-independent algorithm, to obtain a high-order finite element…
A diffusion interface two-phase magnetohydrodynamic model has been used for matched densities in our previous work [1,2], which may limit the applications of the model. In this work, we derive a thermodynamically consistent diffuse…
A unified fluid-structure interaction (FSI) formulation is presented for solid, liquid and mixed membranes. Nonlinear finite elements (FE) and the generalized-alpha scheme are used for the spatial and temporal discretization. The membrane…
In this work we investigate the advantages of multiscale methods in Petrov-Galerkin (PG) formulation in a general framework. The framework is based on a localized orthogonal decomposition of a high dimensional solution space into a low…
A Galerkin finite element method for the membrane elasticity problem on a meshed surface is constructed by using two-dimensional elements extended into three dimensions. The membrane finite element model is established using the intrinsic…
In this paper, we consider the numerical approximation for a diffuse interface model of the two-phase incompressible inductionless magnetohydrodynamics problem. This model consists of Cahn-Hilliard equations, Navier-Stokes equations and…
As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a…
The multimesh finite element method enables the solution of partial differential equations on a computational mesh composed by multiple arbitrarily overlapping meshes. The discretization is based on a continuous--discontinuous function…