اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدA pressure-correction and bound-preserving discretization of the phase-field method for variable density two-phase flows
82
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
In this paper, we present an efficient numerical algorithm for solving the time-dependent Cahn--Hilliard--Navier--Stokes equations that model the flow of two phases with different densities. The pressure-correction step in the projection method consists of a Poisson problem with a modified right-hand side. Spatial discretization is based on discontinuous Galerkin methods with piecewise linear or piecewise quadratic polynomials. Flux and slope limiting techniques successfully eliminate the bulk shift, overshoot and undershoot in the order parameter, which is shown to be bound-preserving. Several numerical results demonstrate that the proposed numerical algorithm is effective and robust for modeling two-component immiscible flows in porous structures and digital rocks.
قيم البحث
اقرأ أيضاً
A two-phase, low-Mach-number flow solver is proposed for variable-density liquid and gas with phase change. The interface is captured using a split Volume-of-Fluid method, which solves the advection of the reference phase, generalized for the case wh
ere the liquid velocity is not divergence-free and both phases exchange mass. A sharp interface is identified by using PLIC. Mass conservation is achieved in the limit of incompressible liquid, but not with the liquid compressibility and mass exchange. This is a relevant modeling choice for two-phase mixtures at near-critical and supercritical pressure conditions for the liquid but away from the mixture critical temperature. Under this thermodynamic environment, the dissolution of lighter gas species into the liquid phase is enhanced and vaporization or condensation can occur simultaneously at different interface locations. The numerical challenge of solving two-phase, supercritical-pressure flows is greater than simpler two-phase solvers because: a) local phase equilibrium is imposed at each interface cell to determine temperature, composition, or surface tension coefficient; b) a real-fluid thermodynamic model is used to obtain fluid properties; and c) necessary phase-wise values for certain variables are obtained via extrapolation techniques. To alleviate the increased numerical cost, the pressure Poisson equation (PPE) used to solve the low-Mach-number flow is split into a constant-coefficient implicit part and a variable-coefficient explicit part. Thus, a Fast Fourier Transform method can be used for the PPE. Various verification tests are performed to show the accuracy and viability of the present approach. The growth of surface instabilities in a binary system composed of liquid n-decane and gaseous oxygen at supercritical pressures for n-decane is analyzed. Other features of supercritical liquid injection are also shown.
In the presence of strong heterogeneities, it is well known that the use of explicit schemes for the transport of species in a porous medium suffers from severe restrictions on the time step. This has led to the development of implicit schemes that a
re increasingly favoured by practitioners for their computational efficiency. The transport equation requires knowledge of the velocity field, which results from an elliptic problem (Darcy problem) that is the most expensive part of the computation. When considering large reservoirs, a cost-effective way of approximating the Darcy problems is using multiscale domain decomposition (MDD) methods. They allow for the pressure and velocity fields to be computed on coarse meshes (large scale), while detailed basis functions are defined locally, usually in parallel, in a much finer grid (small scale). In this work we adopt the Multiscale Robin Coupled Method (MRCM, [Guiraldello, et al., J. Comput. Phys., 355 (2018) pp. 1-21], [Rocha, et al., J. Comput. Phys., (2020) 109316]), which is a generalization of previous MDD methods that allows for great flexibility in the choice of interface spaces. In this article we investigate the combination of the MRCM with implicit transport schemes. A sequentially implicit strategy is proposed, with different trust-region algorithms ensuring the convergence of the transport solver. The method is assessed on several very stringent 2D two-phase problems, demonstrating its stability even for large time steps. It is also shown that the best accuracy is achieved by considering recently introduced non-polynomial interface spaces, since polynomial spaces are not optimal for high-contrast channelized permeability fields.
We present a finite element based variational interface-preserving and conservative phase-field formulation for the modeling of incompressible two-phase flows with surface tension dynamics. The preservation of the hyperbolic tangent interface profile
of the convective Allen-Cahn phase-field formulation relies on a novel time-dependent mobility model. The mobility coefficient is adjusted adaptively as a function of gradients of the velocity and the order parameter in the diffuse interface region in such a way that the free energy minimization properly opposes the convective distortion. The ratio of the convective distortion to the free energy minimization is termed as the convective distortion parameter, which characterizes the deviation from the hyperbolic tangent shape due to the convection effect. The mass conservation is achieved by enforcing a Lagrange multiplier with both temporal and spatial dependence on the phase-field function. We integrate the interface-preserving and conservative phase-field formulation with the incompressible Navier-Stokes equations and the continuum surface tension force model for the simulation of incompressible two-phase flows. A positivity preserving scheme designed for the boundedness and stability of the solution is employed for the variational discretization using unstructured finite elements. We examine the convergence and accuracy of the Allen-Cahn phase-field solver through a generic one-dimensional bistable diffusion-reaction system in a stretching flow. We quantify and systematically assess the relative interface thickness error and the relative surface tension force error with respect to the convective distortion parameter. Two- and three-dimensional rising bubble cases are further simulated to examine the effectiveness of the proposed model on the volume-preserving mean curvature flow and the interface-preserving capability.
We introduce a phase-field method for continuous modeling of cracks with frictional contacts. Compared with standard discrete methods for frictional contacts, the phase-field method has two attractive features: (1) it can represent arbitrary crack ge
ometry without an explicit function or basis enrichment, and (2) it does not require an algorithm for imposing contact constraints. The first feature, which is common in phase-field models of fracture, is attained by regularizing a sharp interface geometry using a surface density functional. The second feature, which is a unique advantage for contact problems, is achieved by a new approach that calculates the stress tensor in the regularized interface region depending on the contact condition of the interface. Particularly, under a slip condition, this approach updates stress components in the slip direction using a standard contact constitutive law, while making other stress components compatible with stress in the bulk region to ensure non-penetrating deformation in other directions. We verify the proposed phase-field method using stationary interface problems simulated by discrete methods in the literature. Subsequently, by allowing the phase field to evolve according to brittle fracture theory, we demonstrate the proposed methods capability for modeling crack growth with frictional contact.
In this work we formulate and test a new procedure, the Multiscale Perturbation Method for Two-Phase Flows (MPM-2P), for the fast, accurate and naturally parallelizable numerical solution of two-phase, incompressible, immiscible displacement in porou
s media approximated by an operator splitting method. The proposed procedure is based on domain decomposition and combines the Multiscale Perturbation Method (MPM) with the Multiscale Robin Coupled Method (MRCM). When an update of the velocity field is called by the operator splitting algorithm, the MPM-2P may provide, depending on the magnitude of a dimensionless algorithmic parameter, an accurate and computationally inexpensive approximation for the velocity field by reusing previously computed multiscale basis functions. Thus, a full update of all multiscale basis functions required by the MRCM for the construction of a new velocity field is avoided. There are two main steps in the formulation of the MPM-2P. Initially, for each subdomain one local boundary value problem with trivial Robin boundary conditions is solved (instead of a full set of multiscale basis functions, that would be required by the MRCM). Then, the solution of an inexpensive interface problem provides the velocity field on the skeleton of the decomposition of the domain. The resulting approximation for the velocity field is obtained by downscaling. We consider challenging two-phase flow problems, with high-contrast permeability fields and water-oil finger growth in homogeneous media. Our numerical experiments show that the use of the MPM-2P gives exceptional speed-up - almost 90% of reduction in computational cost - of two-phase flow simulations. Hundreds of MRCM solutions can be replaced by inexpensive MPM-2P solutions, and water breakthrough can be simulated with very few updates of the MRCM set of multiscale basis functions.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
