Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAnalysis of fully discrete FEM for miscible displacement in porous media with Bear--Scheidegger diffusion tensor
67
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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}$.
rate research
Read More
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.
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.
We analyse a PDE system modelling poromechanical processes (formulated in mixed form using the solid deformation, fluid pressure, and total pressure) interacting with diffusing and reacting solutes in the medium. We investigate the well-posedness of the nonlinear set of equations using fixed-point theory, Fredholms alternative, a priori estimates, and compactness arguments. We also propose a mixed finite element method and rigorously demonstrate the stability of the scheme. Error estimates are derived in suitable norms, and numerical experiments are conducted to illustrate the mechano-chemical coupling and to verify the theoretical rates of convergence.
Homogenization in terms of multiscale limits transforms a multiscale problem with $n+1$ asymptotically separated microscales posed on a physical domain $D subset mathbb{R}^d$ into a one-scale problem posed on a product domain of dimension $(n+1)d$ by introducing $n$ so-called fast variables. This procedure allows to convert $n+1$ scales in $d$ physical dimensions into a single-scale structure in $(n+1)d$ dimensions. We prove here that both the original, physical multiscale problem and the corresponding high-dimensional, one-scale limiting problem can be efficiently treated numerically with the recently developed quantized tensor-train finite-element method (QTT-FEM). The method is based on restricting computation to sequences of nested subspaces of low dimensions (which are called tensor ranks) within a vast but generic virtual (background) discretization space. In the course of computation, these subspaces are computed iteratively and data-adaptively at runtime, bypassing any offline precomputation. For the purpose of theoretical analysis, such low-dimensional subspaces are constructed analytically to bound the tensor ranks vs. error $tau>0$. We consider a model linear elliptic multiscale problem in several physical dimensions and show, theoretically and experimentally, that both (i) the solution of the associated high-dimensional one-scale problem and (ii) the corresponding approximation to the solution of the multiscale problem admit efficient approximation by the QTT-FEM. These problems can therefore be numerically solved in a scale-robust fashion by standard (low-order) PDE discretizations combined with state-of-the-art general-purpose solvers for tensor-structured linear systems. We prove scale-robust exponential convergence, i.e., that QTT-FEM achieves accuracy $tau$ with the number of effective degrees of freedom scaling polynomially in $log tau$.
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
