Related papers: An algorithm for one-dimensional Generalized Porou…
Mixed-dimensional elliptic equations exhibiting a hierarchical structure are commonly used to model problems with high aspect ratio inclusions, such as flow in fractured porous media. We derive general abstract estimates based on the theory…
The present paper deals with the numerical solution of time-fractional advection-diffusion equations involving the Caputo derivative with source term by means of an unconditionally stable implicit finite difference method on quasi-uniform…
We consider fluid flow across a permeable interface within a deformable porous medium. We use mixture theory. The mixture's constituents are assumed to be incompressible in their pure form. We use Hamilton's principle to obtain the…
In this paper, we propose a unified algorithmic framework for solving many known variants of \mds. Our algorithm is a simple iterative scheme with guaranteed convergence, and is \emph{modular}; by changing the internals of a single…
We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems…
We analyze the continuum limit of a thresholding algorithm for motion by mean curvature of one dimensional interfaces in various space-time discrete regimes. The algorithm can be viewed as a time-splitting scheme for the Allen-Cahn equation…
This work outlines a diffuse interface method for the study of fracture and fragmentation in ductile metals at high strain-rates in Eulerian finite volume simulations. The work is based on an existing diffuse interface method capable of…
One of the major issues in the computational mechanics is to take into account the geometrical complexity. To overcome this difficulty and to avoid the expensive mesh generation, geometrically unfitted methods, i.e. the numerical methods…
This paper presents a high-order method for solving an interface problem for the Poisson equation on embedded meshes through a coupled finite element and integral equation approach. The method is capable of handling homogeneous or…
The $hp$-adaptive finite element method (FEM) - where one independently chooses the mesh size ($h$) and polynomial degree ($p$) to be used on each cell - has long been known to have better theoretical convergence properties than either $h$-…
We propose a method that morphs high-orger meshes such that their boundaries and interfaces coincide/align with implicitly defined geometries. Our focus is particularly on the case when the target surface is prescribed as the zero…
Developing high-order numerical schemes for two-phase flow in porous media that preserve key physical properties remains a significant challenge in numerical analysis. In this article, we propose a general framework to construct fully…
We fully generalize a previously-developed computational geometry tool [1] to perform large-scale simulations of arbitrary two-dimensional faceted surfaces $z = h(x,y)$. Our method uses a three-component facet/edge/junction storage model,…
We consider the wave equation $\varepsilon^2(-\partial_t^2 + \Delta)u + f(u) = 0$ for $0<\varepsilon\ll 1$, where $f$ is the derivative of a balanced, double-well potential, the model case being $f(u) = u-u^3$. For equations of this form,…
We design an adaptive unfitted finite element method on the Cartesian mesh with hanging nodes. We derive an hp-reliable and efficient residual type a posteriori error estimate on K-meshes. A key ingredient is a novel hp-domain inverse…
In this paper, we study the numerical approximation of a system of partial dif-ferential equations describing the corrosion of an iron based alloy in a nuclear waste repository. In particular, we are interested in the convergence of a…
Imbibition is a commonly encountered multiphase problem in various fields, and exact prediction of imbibition processes is a key issue for better understanding capillary flow in heterogeneous porous media. In this work, a numerical…
The finite volume correction for a mean-field monomer-dimer system with an attractive interaction are computed for the pressure density, the monomer density and the susceptibility. The results are obtained by introducing a two-dimensional…
In this paper we propose the use of a continuous data assimilation algorithm for miscible flow models in a porous medium. In the absence of initial conditions for the model, observed sparse measurements are used to generate an approximation…
This paper makes the first attempt to apply newly developed upwind GFDM for the meshless solution of two-phase porous flow equations. In the presented method, node cloud is used to flexibly discretize the computational domain, instead of…