The equations of immiscible viscous displacement in a porous medium

1981 ◽  
Vol 59 (1) ◽  
pp. 45-56 ◽  
Author(s):  
T. J. T. Spanos

A statistical theory for the construction of the equations of viscous displacement in a porous medium is considered. This yields a continuum theory for immiscible displacement which can be applied to either a homogeneous or inhomogeneous porous medium. The relative motions of the fluid are considered in terms of the motion of surfaces of constant saturation which are smoothed surfaces at the macroscopic scale considered. The boundary conditions and initial conditions at the injection boundary are considered as well as the boundary conditions and breakthrough conditions at the recovery boundary and the side boundary conditions. The inertial terms are included in the equations and shown to be of importance in describing these initial conditions and the breakthrough conditions.

1981 ◽  
Vol 59 (5) ◽  
pp. 678-687 ◽  
Author(s):  
T. J. T. Spanos

A perturbation of an immiscible displacement process causes relative motion of the two fluids involved. At the macroscopic scale such relative motions are considered to propagate throughout the porous medium in the form of fluid waves. A description of these waves is given on surfaces of constant saturation in a similar fashion to the description of a surface wave propagating on the interface between two fluids. In the porous medium, however, the wave propagation is not restricted to the surface of constant saturation and as a result one obtains a wave equation that is both dissipative and diffusive.A stability analysis is also considered for the immiscible displacement process. Here, a characteristic time for instability to occur can be calculated when the inertial terms are included in the equations of motion. Also a generalization of the wave equations and stability criteria are considered for an inhomogeneous porous medium.


2021 ◽  
Vol 2119 (1) ◽  
pp. 012048
Author(s):  
V V Kuznetsov ◽  
S A Safonov

Abstract This paper presents the results of numerical study of the relationship between micro-and macroscale flows during immiscible displacement in a two-layer porous medium. A feature of the proposed approach is the allowance for large-scale capillarity induced flow due to curvature of the displacement front in macro-inhomogeneous porous medium. The physical mechanisms determining the development of viscous instability in a layer-inhomogeneous porous medium are considered, the methods for suppressing viscous fingers formation based on the stabilization of the displacement front due the action of capillary forces are proposed.


1979 ◽  
Vol 57 (10) ◽  
pp. 1738-1745 ◽  
Author(s):  
T. J. T. Spanos

The viscous surface conditions between two fluids are considered for fluid displacement in a homogeneous porous medium. The term viscous surface is defined here as a mathematical abstraction used to approximate the macroscopic shape of the boundary layer between two fixed saturations of displacing fluid by continuum theory. In the limit as the saturations approach each other, one then obtains a convergence of the viscous surface to a constant saturation contour. This yields a mathematical description of immiscible displacement in porous media which contains a theory of viscous fingering.


2020 ◽  
Vol 4 (41) ◽  
pp. 57-62
Author(s):  
SHAVKAT KLYCHEV ◽  
◽  
BAKHRAMOV SAGDULLA ◽  
VALERIY KHARCHENKO ◽  
VLADIMIR PANCHENKO ◽  
...  

There are needed energy (heat) accumulators to increase the efficiency of solar installations, including solar collectors (water heaters, air heaters, dryers). One of the tasks of designing heat accumulators is to ensure its minimal heat loss. The article considers the problem of determining the distribution of temperatures and heat losses by convection and radiation of the heat insulation-accumulating body (water) system for a ball heat accumulator under symmetric boundary conditions. The problem is solved numerically according to the program developed on the basis of the proposed «gap method». (Research purpose) The research purpose is in determining heat losses by convection and radiation of a two-layer ball heat accumulator with symmetric boundary conditions. (Materials and methods) Authors used the Fourier heat equation for spherical bodies. The article presents the determined boundary and initial conditions for bodies and their surfaces. (Results and discussion) The thickness of the insulation and the volume of the heat accumulator affect the dynamics and values of heat loss. The effect of increasing the thickness of the thermal insulation decreases with increasing its thickness, starting with a certain volume of the heat accumulator or with R > 0.3 meters, the heat losses change almost linearly over time. The dynamics of heat loss decreases with increasing shelf life, but the losses remain large. (Conclusions) Authors have developed a method and program for numerical calculation of heat loss and temperature over time in a spherical two-layer heat accumulator with symmetric boundary conditions, taking into account both falling and intrinsic radiation. The proposed method allows to unify the boundary conditions between contacting bodies.


2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Ji Lin ◽  
Yuhui Zhang ◽  
Chein-Shan Liu

AbstractFor nonlinear third-order three-point boundary value problems (BVPs), we develop two algorithms to find solutions, which automatically satisfy the specified three-point boundary conditions. We construct a boundary shape function (BSF), which is designed to automatically satisfy the boundary conditions and can be employed to develop new algorithms by assigning two different roles of free function in the BSF. In the first algorithm, we let the free functions be complete functions and the BSFs be the new bases of the solution, which not only satisfy the boundary conditions automatically, but also can be used to find solution by a collocation technique. In the second algorithm, we let the BSF be the solution of the BVP and the free function be another new variable, such that we can transform the BVP to a corresponding initial value problem for the new variable, whose initial conditions are given arbitrarily and terminal values are determined by iterations; hence, we can quickly find very accurate solution of nonlinear third-order three-point BVP through a few iterations. Numerical examples confirm the performance of the new algorithms.


2008 ◽  
Vol 2008 ◽  
pp. 1-7 ◽  
Author(s):  
Mantas Povilaitis ◽  
Egidijus Urbonavičius

An issue of the stratified atmospheres in the containments of nuclear power plants is still unresolved; different experiments are performed in the test facilities like TOSQAN and MISTRA. MASPn experiments belong to the spray benchmark, initiated in the containment atmosphere mixing work package of the SARNET network. The benchmark consisted of MASP0, MASP1 and MASP2 experiments. Only the measured depressurisation rates during MASPn were available for the comparison with calculations. When the analysis was performed, the boundary conditions were not clearly defined therefore most of the attention was concentrated on MASP0 simulation in order to develop the nodalisation scheme and define the initial and boundary conditions. After achieving acceptable agreement with measured depressurisation rate, simulations of MASP1 and MASP2 experiments were performed to check the influence of sprays. The paper presents developed nodalisation scheme of MISTRA for the COCOSYS code and the results of analyses. In the performed analyses, several parameters were considered: initial conditions, loss coefficient of the junctions, initial gradients of temperature and steam volume fraction, and characteristic length of structures. Parametric analysis shows that in the simulation the heat losses through the external walls behind the lower condenser installed in the MISTRA facility determine the long-term depressurisation rate.


2016 ◽  
Vol 792 ◽  
pp. 5-35 ◽  
Author(s):  
Giuseppe A. Zampogna ◽  
Alessandro Bottaro

The interaction between a fluid flow and a transversely isotropic porous medium is described. A homogenized model is used to treat the flow field in the porous region, and different interface conditions, needed to match solutions at the boundary between the pure fluid and the porous regions, are evaluated. Two problems in different flow regimes (laminar and turbulent) are considered to validate the system, which includes inertia in the leading-order equations for the permeability tensor through a Oseen approximation. The components of the permeability, which characterize microscopically the porous medium and determine the flow field at the macroscopic scale, are reasonably well estimated by the theory, both in the laminar and the turbulent case. This is demonstrated by comparing the model’s results to both experimental measurements and direct numerical simulations of the Navier–Stokes equations which resolve the flow also through the pores of the medium.


Sign in / Sign up

Export Citation Format

Share Document