Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConvergence of a decoupled mixed FEM for miscible displacement in interfacial porous media
137
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper, we study the stability and convergence of a decoupled and linearized mixed finite element method (FEM) for incompressible miscible displacement in a porous media whose permeability and porosity are discontinuous across some interfaces. We show that the proposed scheme has optimal-order convergence rate unconditionally, without restriction on the grid ratio (between the time-step size and spatial mesh size). Previous works all required certain restrictions on the grid ratio except for the problem with globally smooth permeability and porosity. Our idea is to introduce an intermediate system of elliptic interface problems, whose solution is uniformly regular in each subdomain separated by the interfaces and its finite element solution coincides with the fully discrete solution of the original problem. In order to prove the boundedness of the fully discrete solution, we study the finite element discretization of the intermediate system of elliptic interface problems.
rate research
Read More
Fully discrete Galerkin finite element methods are studied for the equations of miscible displacement in porous media with the commonly-used Bear--Scheidegger diffusion-dispersion tensor: $$ D({bf u}) = gamma d_m I + |{bf u}|bigg( alpha_T I + (alpha_L - alpha_T) frac{{bf u} otimes {bf u}}{|{bf u}|^2}bigg) , . $$ Previous works on optimal-order $L^infty(0,T;L^2)$-norm error estimate required the regularity assumption $ abla_xpartial_tD({bf u}(x,t)) in L^infty(0,T;L^infty(Omega))$, while the Bear--Scheidegger diffusion-dispersion tensor is only Lipschitz continuous even for a smooth velocity field ${bf u}$. In terms of the maximal $L^p$-regularity of fully discrete finite element solutions of parabolic equations, optimal error estimate in $L^p(0,T;L^q)$-norm and almost optimal error estimate in $L^infty(0,T;L^q)$-norm are established under the assumption of $D({bf u})$ being Lipschitz continuous with respect to ${bf u}$.
We develop a theory for the problem of high pressure air injection into deep reservoirs containing light oil. Under these conditions, the injected fluid (oxygen + inert components) is completely miscible with the oil in the reservoir. Moreover, exothermic reactions between dissolved oxygen and oil are possible. We use Kovals model to account for the miscibility of the components, such that the fractional-flow functions resemble the ones from Buckley-Leverett flow. This allows to decompose the solution of this problem into a series of waves. We then proceed to obtain full analytical solutions in each wave. Of particular interest {is the case where} the combustion wave presents a singularity in its internal wave profile. Evaluation of the variables of the problem at the singular point determines the macroscopic parameters of the wave, i.e., combustion temperature, wave speed and downstream oil fraction. The waves structure was observed previously for reactive immiscible displacement and we describe it here for the first time for reactive miscible displacement of oil. We validate the developed theory using numerical simulations.
Accurate numerical simulations of interaction between fluid and solid play an important role in applications. The task is challenging in practical scenarios as the media are usually highly heterogeneous with very large contrast. To overcome this computational challenge, various multiscale methods are developed. In this paper, we consider a class of linear poroelasticity problems in high contrast heterogeneous porous media, and develop a mixed generalized multiscale finite element method (GMsFEM) to obtain a fast computational method. Our aim is to develop a multiscale method that is robust with respect to the heterogeneities and contrast of the media, and gives a mass conservative fluid velocity field. We will construct decoupled multiscale basis functions for the elastic displacement as well as fluid velocity. Our multiscale basis functions are local. The construction is based on some suitable choices of local snapshot spaces and local spectral decomposition, with the goal of extracting dominant modes of the solutions. For the pressure, we will use piecewise constant approximation. We will present several numerical examples to illustrate the performance of our method. Our results indicate that the proposed method is able to give accurate numerical solutions with a small degree of freedoms.
We propose and analyze novel adaptive algorithms for the numerical solution of elliptic partial differential equations with parametric uncertainty. Four different marking strategies are employed for refinement of stochastic Galerkin finite element approximations. The algorithms are driven by the energy error reduction estimates derived from two-level a posteriori error indicators for spatial approximations and hierarchical a posteriori error indicators for parametric approximations. The focus of this work is on the mathematical foundation of the adaptive algorithms in the sense of rigorous convergence analysis. In particular, we prove that the proposed algorithms drive the underlying energy error estimates to zero.
We present a novel approach to the simulation of miscible displacement by employing adaptive enriched Galerkin finite element methods (EG) coupled with entropy residual stabilization for transport. In particular, numerical simulations of viscous fingering instabilities in heterogeneous porous media and Hele-Shaw cells are illustrated. EG is formulated by enriching the conforming continuous Galerkin finite element method (CG) with piecewise constant functions. The method provides locally and globally conservative fluxes, which is crucial for coupled flow and transport problems. Moreover, EG has fewer degrees of freedom in comparison with discontinuous Galerkin (DG) and an efficient flow solver has been derived which allows for higher order schemes. Dynamic adaptive mesh refinement is applied in order to save computational cost for large-scale three dimensional applications. In addition, entropy residual based stabilization for high order EG transport systems prevents any spurious oscillations. Numerical tests are presented to show the capabilities of EG applied to flow and transport.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
