Multiscale finite-volume method for parabolic problems arising from compressible multiphase flow in porous media

The multiscale finite-volume (MSFV) method was originally developed for the solution of heterogeneous elliptic problems with reduced computational cost. Recently, some extensions of this method for parabolic problems have been proposed. These extensions proved effective for many cases, however, they are neither general nor completely satisfactory. For instance, they are not suitable for correctly capturing the transient behavior described by the parabolic pressure equation. In this paper, we present a general multiscale finite-volume method for parabolic problems arising, for example, from compressible multiphase flow in porous media. Opposed to previous methods, here, the basis and correction functions are solutions of full parabolic governing equations in localized domains. At the same time, to enhance the computational efficiency of the scheme, the basis functions are kept pressure independent and do not have to be recalculated as pressure evolves. This general approach requires no additional assumptions and its good efficiency and high accuracy is demonstrated for various challenging test cases. Finally, to improve the quality of the results and also to extend the scheme for highly anisotropic heterogeneous problems, it is combined with the iterative MSFV (i-MSFV) method for parabolic problems. As one iterates, the i-MSFV solutions of compressible multiphase problems (parabolic problems) converge to the corresponding fine-scale reference solutions in the same way as demonstrated recently for incompressible cases (elliptic problems). Therefore, the proposed MSFV method can also be regarded as an efficient linear solver for parabolic problems and studies of its efficiency are presented for many test cases.     

裂缝性介质多尺度深度学习模型

结合人工神经网络建立裂缝介质多尺度深度学习流动模型.基于一套粗网格和一套细网格,通过在粗网格上训练数据,多尺度神经网络能够以较少的自由度训练出准确的神经网络.并在粗网格上通过求解局部流动问题获得多尺度基函数,结合神经网络进一步得到精细网格的解.基于离散裂缝的流动方程可视为多层网络,网络层数依赖于求解时间步数.阐述裂缝介质多尺度机器学习数值计算格式的建立,介绍如何使用多尺度算法构建离散裂缝模型的多尺度基函数,并采用超样本技术进一步提高计算准确性.数值结果表明,多尺度有限元算法与机器学习结合是一种有效的流体流动模拟算法.     

多孔复合介质周期结构热传导和质扩散问题的多尺度数值方法

本文针对一类复杂的多孔复合介质的热传导和质扩散问题,给出具体的多尺度渐近展开公式,并在此基础上设计了有限元算法格式,它是宏观和细观相结合的数值方法。理论分析和数值实验均表明：多尺度数值方法对求解多孔复合介质周期结构的热传导和质扩散问题是可行的和有效的。     

Local–global multiscale model reduction for flows in high-contrast heterogeneous media

In this paper, we study model reduction for multiscale problems in heterogeneous high-contrast media. Our objective is to combine local model reduction techniques that are based on recently introduced spectral multiscale finite element methods (see [19]) with global model reduction methods such as balanced truncation approaches implemented on a coarse grid. Local multiscale methods considered in this paper use special eigenvalue problems in a local domain to systematically identify important features of the solution. In particular, our local approaches are capable of homogenizing localized features and representing them with one basis function per coarse node that are used in constructing a weight function for the local eigenvalue problem. Global model reduction based on balanced truncation methods is used to identify important global coarse-scale modes. This provides a substantial CPU savings as Lyapunov equations are solved for the coarse system. Typical local multiscale methods are designed to find an approximation of the solution for any given coarse-level inputs. In many practical applications, a goal is to find a reduced basis when the input space belongs to a smaller dimensional subspace of coarse-level inputs. The proposed approaches provide efficient model reduction tools in this direction. Our numerical results show that, only with a careful choice of the number of degrees of freedom for local multiscale spaces and global modes, one can achieve a balanced and optimal result.     

Numerical methods for multiscale transport equations and application to two-phase porous media flow

We discuss numerical methods for linear and nonlinear transport equations with multiscale velocity fields. These methods are themselves multiscaled in nature in the sense that they use macro and micro grids, multiscale test functions. We demonstrate the efficiency of these methods and apply them to two-phase flow in heterogeneous porous media.     

Adaptive iterative multiscale finite volume method

The multiscale finite volume (MSFV) method is a computationally efficient numerical method for the solution of elliptic and parabolic problems with heterogeneous coefficients. It has been shown for a wide range of test cases that the MSFV results are in close agreement with those obtained with a classical (computationally expensive) technique. The method, however, fails to give accurate results for highly anisotropic heterogeneous problems due to weak localization assumptions. Recently, a convergent iterative MSFV (i-MSFV) method was developed to enhance the quality of the multiscale results by improving the localization conditions. Although the i-MSFV method proved to be efficient for most practical problems, it is still favorable to improve the localization condition adaptively, i.e. only for a sub-domain where the original MSFV localization conditions are not acceptable, e.g. near shale layers and long coherent structures with high permeability contrasts. In this paper, a space–time adaptive i-MSFV (ai-MSFV) method is introduced. It is shown how to improve the MSFV results adaptively in space and simulation time. The fine-scale smoother, which is necessary for convergence of the i-MSFV method, is also applied locally. Finally, for multiphase flow problems, two criteria are investigated for adaptively updating the MSFV interpolation functions: (1) a criterion based on the total mobility change for the transient coefficients and (2) a criterion based on the pressure equation residual for the accuracy of the results. For various challenging test cases it is demonstrated that iterations in order to obtain accurate results even for highly anisotropic heterogeneous problems are required only in small sub-domains and not everywhere. The findings show that the error introduced in the MSFV framework can be controlled and improved very efficiently with very little additional computational cost compared to the original, non-iterative MSFV method.     

XFEM modeling of hydraulic fracture in porous rocks with natural fractures

Hydraulic fracture (HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method (XFEM) coupling with Biot theory is developed to study the HF in permeable rocks with natural fractures (NFs). In the recent XFEM based computational HF models, the fluid flow in fractures and interstitials of the porous media are mostly solved separately, which brings difficulties in dealing with complex fracture morphology. In our new model the fluid flow is solved in a unified framework by considering the fractures as a kind of special porous media and introducing Poiseuille-type flow inside them instead of Darcy-type flow. The most advantage is that it is very convenient to deal with fluid flow inside the complex fracture network, which is important in shale gas extraction. The weak formulation for the new coupled model is derived based on virtual work principle, which includes the XFEM formulation for multiple fractures and fractures intersection in porous media and finite element formulation for the unified fluid flow. Then the plane strain Kristianovic-Geertsma-de Klerk (KGD) model and the fluid flow inside the fracture network are simulated to validate the accuracy and applicability of this method. The numerical results show that large injection rate, low rock permeability and isotropic in-situ stresses tend to lead to a more uniform and productive fracture network.     

A multiscale coupling method for the modeling of dynamics of solids with application to brittle cracks

We present a multiscale model for numerical simulations of dynamics of crystalline solids. The method combines the continuum nonlinear elasto-dynamics model, which models the stress waves and physical loading conditions, and molecular dynamics model, which provides the nonlinear constitutive relation and resolves the atomic structures near local defects. The coupling of the two models is achieved based on a general framework for multiscale modeling – the heterogeneous multiscale method (HMM). We derive an explicit coupling condition at the atomistic/continuum interface. Application to the dynamics of brittle cracks under various loading conditions is presented as test examples.     

A stochastic mixed finite element heterogeneous multiscale method for flow in porous media

A computational methodology is developed to efficiently perform uncertainty quantification for fluid transport in porous media in the presence of both stochastic permeability and multiple scales. In order to capture the small scale heterogeneity, a new mixed multiscale finite element method is developed within the framework of the heterogeneous multiscale method (HMM) in the spatial domain. This new method ensures both local and global mass conservation. Starting from a specified covariance function, the stochastic log-permeability is discretized in the stochastic space using a truncated Karhunen–Loève expansion with several random variables. Due to the small correlation length of the covariance function, this often results in a high stochastic dimensionality. Therefore, a newly developed adaptive high dimensional stochastic model representation technique (HDMR) is used in the stochastic space. This results in a set of low stochastic dimensional subproblems which are efficiently solved using the adaptive sparse grid collocation method (ASGC). Numerical examples are presented for both deterministic and stochastic permeability to show the accuracy and efficiency of the developed stochastic multiscale method.     

Fractional dynamics of tracer transport in fractured media from local to regional scales

Tracer transport through fractured media exhibits concurrent direction-dependent super-diffusive spreading along high-permeability fractures and sub-diffusion caused by mass transfer between fractures and the rock matrix. The resultant complex dynamics challenge the applicability of conventional physical models based on Fick’s law. This study proposes a multi-scaling tempered fractional-derivative (TFD) model to explore fractional dynamics for tracer transport in fractured media. Applications show that the TFD model can capture anomalous transport observed in small-scale single fractures, intermediate-scale fractured aquifers, and two-dimensional large-scale discrete fracture networks. Tracer transport in fractured media from local (0.255-meter long) to regional (400-meter long) scales therefore can be quantified by a general fractional-derivative model. Fractional dynamics in fractured media can be scale dependent, owning to 1) the finite length of fractures that constrains the large displacement of tracers, and 2) the increasing mass exchange capacity along the travel path that enhances sub-diffusion.     

随机多孔复相介质稳态温度场的多尺度数值模拟

本文就随机多孔复相介质稳态温度场,提出了多尺度分析的基本思想和具体的计算公式,它们是建立在确定性温度场的多尺度算法和Ｋａｒｈｕｎｅｎ－Ｌｏｅｖｅ展开式基础上．文中算例表明这一算法的有效性。     

Characterization of fractured permeable porous media using relaxation-weighted imaging techniques

The inversion time (TI) weighted inversion-recovery spin-echo (IRSE) imaging technique was used for laboratory characterization of fractures and porous matrix in permeable media. The method is based on the fact that relaxation rates in fractures and those in surrounding porous matrix are substantially different. Thus, by selectively suppressing imaging signals from either fractures or porous matrix, one can highlight fluid distributions in either of these regions so that fractures can be unambiguously identified and surrounding porous structures can be characterized. The advantages over conventional NMR spin-echo imaging are demonstrated with various fractured porous rock types. A technique to speed acquisitions of images of fluid distributions in porous matrix is also demonstrated.     

非均一多孔介质中的水热迁移研究

孔隙裂隙非均一多孔介质中水热迁移的研究，国际上只是近年来才开始取得明显的进展。在发育裂隙的孔隙岩层中同时存在着两种渗流系统：孔隙总体积较大、渗透性相对弱的多孔岩块和总体积较小、渗透性却相对较强的分割多孔块体的裂隙。从而提出了“孔隙一裂隙二重性”假定，即地下水主要贮存在孔隙中，而水的运动主要在裂隙中进行，用一个一阶量描述孔隙裂隙间水流（热流）的传递耦合项。本文导出描述孔隙─裂隙岩层中的水流和热迁移的基本微分方程，建立起相应的数学模型，并成功地用于我国西藏羊八井热田分布参数模型的水热迁移研究。     

An efficient multi-point flux approximation method for Discrete Fracture–Matrix simulations

We consider a control volume discretization with a multi-point flux approximation to model Discrete Fracture–Matrix systems for anisotropic and fractured porous media in two and three spatial dimensions. Inspired by a recently introduced approach based on a two-point flux approximation, we explicitly account for the fractures by representing them as hybrid cells between the matrix cells. As well as simplifying the grid generation, our hybrid approach excludes small cells in the intersection of the fractures and hence avoids severe time-step restrictions associated with small cells. Excluding the small cells also reduces the condition number of the discretization matrix. For examples involving realistic anisotropy ratios in the permeability, numerical results show significant improvement compared to existing methods based on two-point flux approximations. We also investigate the hybrid method by studying the convergence rates for different apertures and fracture/matrix permeability ratios. Finally, the effect of removing the cells in the intersections of the fractures are studied. Together, these examples demonstrate the efficiency, flexibility and robustness of our new approach.     

生物大分子多尺度理论和计算方法

分子模拟是研究生物大分子的重要手段. 过去二十年来, 人们将分子模拟与实验研究相结合, 揭示出生物大分子结构和动力学方面的诸多重要性质. 传统分子模拟主要采用全原子分子模型或各种粗粒化的分子模型. 在实际应用中, 传统分子模拟方法通常存在精度或效率瓶颈, 一定程度上限制了其应用范围. 近年来, 多尺度分子模型越来越受到人们的关注. 多尺度分子模型基于统计力学原理, 将全原子模型和粗粒化模型相耦合, 有望克服传统分子模拟方法中的精度/效率瓶颈, 进而拓展分子模拟在生物大分子研究中的应用范围. 根据模型之间的耦合方式, 近年来发展起来的多尺度分子模拟方法可归纳为如下四种类型: 混合分辨多尺度模型、并行耦合多尺度模型、单向耦合多尺度模型、以及自学习多尺度模型. 本文将对上述四类多尺度模型做简要介绍, 并讨论其主要优缺点、应用范围以及进一步发展方向.     

表征体元尺度渗流的离散统一动理学格式

结合表征体元尺度的通用渗流模型,提出离散统一动理学格式（DUGKS）渗流方法,分别用均匀网格和非均匀网格计算二维Poiseuille、Couette、方腔流等经典渗流问题,检验DUGKS渗流方法的有效性和非均匀网格应用的优势,将DUGKS渗流方法应用到裂缝系统中.     

Multiscale combustion and turbulence

Multiscale physics is the interaction of different physical processes occurring at largely separated scales. In combustion, many elementary reactions combine to only a few, but still have separated time scales. In flames, owing to the presence of diffusion, time scales manifest themselves as length scales, i.e. thicknesses of reaction layers embedded within each other. For premixed flames there results a single velocity scale, the laminar burning velocity, which in turn defines a flame thickness and a flame time as global length and time scales, respectively. The laminar burning velocity represents the simplest microscale model to be used at a premixed combustion interface.While combustion is a multiscale process, this is not so evident for turbulence. Based on the picture of a cascade process traditional turbulent closure approximations treat turbulence as a single-scale problem. Attempts to model turbulent combustion in the same way by using methods developed for non-reacting turbulent flows therefore must fail, because they ignore the multiscale nature of combustion.There is, however, a long tradition and much progress in multiscale modeling of combustion, both on the macroscale as well as on the microscale level. Unfortunately much of that work is conceived only in its particular context, not as part of a multiscale approach. For instance, papers in the TURBULENT FLAMES Colloquium and the FIRE RESEARCH Colloquium at this and at previous Combustion Symposia often take the viewpoint of macroscale modeling only, while REACTION KINETICS and LAMINAR FLAMES concentrate on microscale aspects. What seems to be needed is a more explicit reference to the needs of models developed in the other parts of the community. Furthermore, research is needed to develop suitable definitions of the interface between macroscale and microscale models.     

Dimension reduction method for ODE fluid models

We develop a new dimension reduction method for large size systems of ordinary differential equations (ODEs) obtained from a discretization of partial differential equations of viscous single and multiphase fluid flow. The method is also applicable to other large-size classical particle systems with negligibly small variations of particle concentration. We propose a new computational closure for mesoscale balance equations based on numerical iterative deconvolution. To illustrate the computational advantages of the proposed reduction method, we use it to solve a system of smoothed particle hydrodynamic ODEs describing single-phase and two-phase layered Poiseuille flows driven by uniform and periodic (in space) body forces. For the single-phase Poiseuille flow driven by the uniform force, the coarse solution was obtained with the zero-order deconvolution. For the single-phase flow driven by the periodic body force and for the two-phase flows, the higher-order (the first- and second-order) deconvolutions were necessary to obtain a sufficiently accurate solution.     

缝洞型介质流动模拟的多尺度分解法

缝洞型介质通常具有非均质性强、结构多尺度的特征.传统数值方法在解决此类多尺度流动问题时,难以兼顾计算精度与计算效率,无法实际应用.对此,本文提出了多孔介质流体流动的多尺度分解法,并应用于缝洞介质流动模拟,能够大幅减小计算的复杂度,同时,可以通过控制均化程度控制计算精度.该方法将求解空间分为若干个子空间的正交直和,从而获得一个近线性的计算复杂度;以分层计算的方式实现了快速计算,另外这种方法是一种无网格方法,具有较好的地层适应性.同时,采用离散缝洞模型简化缝洞结构,进一步提高了计算效率.详细阐述了基于多尺度分解法的多孔介质流体流动数值计算格式的建立,重点介绍了如何在不同的层次上计算基函数.数值结果表明,本文提出的计算方法不仅能够准确捕捉多孔介质中的精细流动特征,而且具有很高的计算效率,是一种有效的流动模拟方法.     

Variational multiscale turbulence modelling in a high order spectral element method

In the variational multiscale (VMS) approach to large eddy simulation (LES), the governing equations are projected onto an a priori scale partitioning of the solution space. This gives an alternative framework for designing and analyzing turbulence models. We describe the implementation of the VMS LES methodology in a high order spectral element method with a nodal basis, and discuss the properties of the proposed scale partitioning. The spectral element code is first validated by doing a direct numerical simulation of fully developed plane channel flow. The performance of the turbulence model is then assessed by several coarse grid simulations of channel flow at different Reynolds numbers.