Numerical simulation and homogenization of two-phase flow in heterogeneous porous media

A mathematically rigorous method of homogenization is presented and used to analyze the equivalent behavior of transient flow of two incompressible fluids through heterogeneous media. Asymptotic expansions and H-convergence lead to the definition of a global or effective model of an equivalent homogeneous reservoir. Numerical computations to obtain the homogenized coefficients of the entire reservoir have been carried out via a finite element method. Numerical experiments involving the simulation of incompressible two-phase flow have been performed for each heterogeneous medium and for the homogenized medium as well as for other averaging methods. The results of the simulations are compared in terms of the transient saturation contours, production curves, and pressure distributions. Results obtained from the simulations with the homogenization method presented show good agreement with the heterogeneous simulations.     

Static Model Describing the Inhomogeneous Deformation of Unidirectional Fibrous Composites under Plane Strain

A mathematical model is proposed to describe the static deformation of a unidirectional fiber-reinforced composite under plane-strain conditions in a plane perpendicular to the fiber direction. The model is based on the stochastic static equations for an elastic homogeneous two-component medium with nonzero body forces. Applying the method of conditional moments and double Fourier transform yields a nonlocal model in the form of integro-differential equations. Expanding the integral kernels into series in expansion coefficients yields differential equations whose order depends on the number of terms in the series. The zero-order approximation leads to the theory of effective moduli, the first-order approximation to the refined theory of effective moduli, and the second-order approximation to fourth-order equilibrium equations for mean displacements (mathematical expectations) and formulas for displacements, strains, and stresses averaged over the composite and its components. All the coefficients in theses formulas are expressed in terms of the elastic constants of the components and the geometric parameters of the structure. The model is valid for heavy stress gradients. In a particular case, the static theory of two-component elastic mixtures follows from the model__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 2, pp. 41–51, February 2005.     

An Asymptotic Model of Seismic Reflection from a Permeable Layer

Analysis of compression wave propagation in a poroelastic medium predicts a peak of reflection from a high-permeability layer in the low-frequency end of the spectrum. An explicit formula expresses the resonant frequency through the elastic moduli of the solid skeleton, the permeability of the reservoir rock, the fluid viscosity and compressibility, and the reservoir thickness. This result is obtained through a low-frequency asymptotic analysis of Biot’s model of poroelasticity. A review of the derivation of the main equations from the Hooke’s law, momentum and mass balance equations, and Darcy’s law suggests an alternative new physical interpretation of some coefficients of the classical poroelasticity. The velocity of wave propagation, the attenuation factor, and the wave number are expressed in the form of power series with respect to a small dimensionless parameter. The absolute value of this parameter is equal to the product of the kinematic reservoir fluid mobility and the wave frequency. Retaining only the leading terms of the series leads to explicit and relatively simple expressions for the reflection and transmission coefficients for a planar wave crossing an interface between two permeable media, as well as wave reflection from a thin highly permeable layer (a lens). Practical applications of the obtained asymptotic formulae are seismic modeling, inversion, and attribute analysis.     

A semi-analysis method of differential equations with variable coefficients under complicated boundary conditions

IntroductionInthemostcases,asfarasthesolutiontothedifferentialequationwithvariablecoefficientsisconcerned ,therehasbeennoanysatisfactoryanswer.Forthisreason ,thescholarsbothathomeandabroadhavedonesomeresearchjustasindicatedinRefs.[1 ] ,[2 ]and [3 ] ,etc.Theauthorsinthispaperhaveusedthefiniteelementmodelandthemetricfunctiontheoryincombinationwithadjustableparametermodelwithvariablecoefficientstosimplifythesolutiontothedifferentialequationwithvariablecoefficientsintothefeaturevalueaswellasthefea…     

一类偏微分方程的Hamilton正则表示

主要给出一系列关于力学中的偏微分方程的无穷维Ｈａｍｉｌｔｏｎ正则表示．其中包括变系数线性偏微分方程,ＫｄＶ方程,ＭＫｄＶ方程,ＫＰ方程,Ｂｏｕｓｉｎｅｓｑ方程等的无穷维Ｈａｍｉｌｔｏｎ正则表示．     

Influence of interfaces on effective properties of nanomaterials with stochastically distributed spherical inclusions

The method of conditional moments is generalized to include evaluation of the effective elastic properties of particulate nanomaterials and to investigate the size effect in those materials. Determining the effective constants necessitates finding a stochastically averaged solution to the fundamental equations of linear elasticity coupled with surface/interface conditions (Gurtin–Murdoch model). To obtain such a solution the system of governing stochastic differential equations is first transformed to an equivalent system of stochastic integral equations. Using statistical averaging, the boundary-value problem is then converted to an infinite system of linear algebraic equations. A two-point approximation is considered and the stress fluctuations within the inclusions are neglected in order to obtain a finite system of algebraic equations in terms of component-average strains. Closed-form expressions are derived for the effective moduli of a composite consisting of a matrix and randomly distributed spherical inhomogeneities. As a numerical example a nanoporous material is investigated assuming a model in which the interface effects influence only the bulk modulus of the material. In that model the resulting shear modulus is the same as for the material without surface effects. Dependence of the bulk moduli on the radius of nanopores and on the pore volume fraction is analyzed. The results are compared to, and discussed in the context of other theoretical predictions.     

裂缝性油藏流固耦合渗流

本文给出了考虑介质变形的双重孔隙介质流固耦合渗流模型,并考虑渗流参数随有效应力而变化的非线性双重孔隙介质流固耦合渗流,在此基础上,本文还推导了双重孔隙介质非线性系数非线性等流固耦合流流计算,并给出了算例。     

Mathematical simulation of two-phase fluid flow through a water-flooded porous reservoir with sediment formation

The results of numerical investigations of the problem of flow through a porous reservoir in the process of its water flooding with addition of a sediment-forming component are given. The mathematical model is based on the mass conservation laws for each of the considered phases and components supplemented with the equations of motion and constitutive relations necessary to close the system of equations. In solving the problem, an empirical dependence of the sediment formation intensity on the content of the sediment-forming component in the aqueous solution with allowance for variation in the effective porosity of the medium is used. The main features in solving the sediment-formation problem are distinguished using the empirical dependence and a contrastive analysis of the effect of choosing this dependence on the solution results is carried out. It is shown that the neglect of the experimental results in the mathematical formulation can lead to not only unjustified overestimated results in realizing the method but also give a distorted pattern of the entire process of sediment formation in fluid flow through a water-flooded porous reservoir.     

裂缝性低渗透油藏流-固耦合理论与数值模拟

根据裂缝性低渗油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型。将渗流力学与弹塑性力学相结合,建立裂缝性低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.通过数值模拟对一实际井网开发过程中孔隙度、渗透率的变化以及开发指标进行计算,并和刚性模型以及双重介质模型的计算结果进行了分析比较.     

Effective characteristics of corrugated plates

Corrugated plates are widely used in modern constructions and structures, because they, in contrast to plane plates, possess greater rigidity. In many cases, such a plate can be modeled by a homogeneous anisotropic plate with certain effective flexural and tensional rigidities. Depending on the geometry of corrugations and their location, the equivalent homogeneous plate can also have rigidities of mutual influence. These rigidities allow one to take into account the influence of bending moments on the strain in the midplane and, conversely, the influence of longitudinal strains on the plate bending [1]. The behavior of the corrugated plate under the action of a load normal to the midsurface is described by equations of the theory of flexible plates with initial deflection. These equations form a coupled system of nonlinear partial differential equations with variable coefficients [2]. The dependence of the coefficients on the coordinates is determined by the corrugation geometry. In the case of a plate with periodic corrugation, the coefficients significantly vary within one typical element and depend on the values of local variables determined in each of the typical elements. There is a connection between the local and global variables, and therefore, the functions of local coordinates are simultaneously functions of global coordinates, which are sometimes called rapidly oscillating functions [3].One of the methods for solving the equations with rapidly oscillating coefficients is the asymptotic method of small geometric parameter. The standard procedure of this method usually includes preparatory stages. At the first stage, as a rule, a rectangular periodicity cell is distinguished to be a typical element. At the second stage, the scale of global coordinates is changed so that the rectangular structure periodicity cells became square cells of size l × l. The third stage consists in passing to the dimensionless global coordinates relative to the plate characteristic dimension L. As a result, the dependence between the new local variables and the new scaled dimensionless variables is such that the factor 1/α, where α=l/L ? 1 is a small geometric parameter, appears in differentiating any function of the local coordinate with respect to the global coordinate. After this, the solution of the problem in new coordinates is sought as an asymptotic expansion in the small geometric parameter [1], [4–10].We note that, in the small geometric parameter method, the asymptotic series simultaneously have the form of expansions in the gradients of functions depending only on the global coordinates. This averaging procedure can be applied to linear and nonlinear boundary value problems for differential equations with variable periodic coefficients for which the periodicity cell can be affinely transformed into the periodicity cube. In the case of an arbitrary dependence of the coefficients on the coordinates (including periodic dependence), another averaging technique can be used in linear problems. This technique is based on the possibility of the integral representation of the solution of the original problem for the linear equation with variables coefficients in terms of the solution of the same problem for an equation of the same type but with constant coefficients [11–13]. The integral representation implies that the solution of the original problem can be represented in the form of the series in the gradients of the solution of the problem for the equation with constant coefficients [13].The aim of the present paper is to develop methods for calculating effective characteristics of corrugated plates. To this end, we first write out the equilibrium equations for a flexible anisotropic plate, which is inhomogeneous in the thickness direction and in the horizontal projection, with an initial deflection. We write these equations in matrix form, which allows one to significantly reduce the length of the expressions and to simplify further calculations. After this, we average the initial matrix equations with variable coefficients. The averaging procedure implies the statement of problems such that, after solving them, we can calculate the desired effective characteristics. By way of example, we consider the case of a corrugated plate made of a homogeneous isotropic material whose corrugations are hexagonal in the horizontal projection. In this case, we obtain approximate expressions for the components of the effective tensors of flexural rigidity and longitudinal compliance and expressions for the effective plate thickness.     

多相材料传热微结构的多目标优化设计

提出基于散热弱度的材料微结构热传导性能的预测方法,分别从理论和数值上验证该方法与均匀化方法的等效性;推导出微结构等效热传导系数的灵敏度计算格式,建立传热微结构拓扑优化的数学模型.以二维、三维多相材料等效热传导系数的加权组合为目标,采用凸规划对偶优化算法和二次型周长约束进行材料微结构的设计和材料分布的棋盘格控制.数值算例表明基于散热弱度的传热材料微结构设计是可行、有效的,可以为实际的材料设计提供依据.     

Construction of Dynamically Equivalent Time-Invariant Forms for Time-Periodic Systems

In this study dynamically equivalent time-invariant forms are obtained for linear and non-linear systems with periodically varying coefficients via Lyapunov–Floquet (L–F) transformation. These forms are equivalent in the sense that the local stability and bifurcation characteristics are identical for both systems in the entire parameter space. It is well known that the L–F transformation converts a linear periodic first order system into a time-invariant one. In the first part of this study a set of linear second order periodic equations is converted into an equivalent set of time-independent second order equations through a sequence of linear transformations. Then the transformations are applied to a time-periodic quadratic Hamiltonian to obtain its equivalent time-invariant form. In the second part, time-invariant forms of nonlinear equations are studied. The application of L–F transformation to a quasi-linear periodic equation converts the linear part to a time-invariant form and leaves the non-linear part with time-periodic coefficients. Dynamically equivalent time-invariant forms are obtained via time-periodic center manifold reduction and time-dependent normal form theory. Such forms are constructed for general hyperbolic systems and for some simple critical cases, including that of one zero eigenvalue and a purely imaginary pair. As a physical example of these techniques, a single and a double inverted pendulum subjected to periodic parametric excitation are considered. The results thus obtained are verified by numerical simulation.     

卷积型GD半解析法及矩形薄板瞬态响应解

卷积型的Gurtin变分原理是目前在数学上唯一能和动力学初值问题完全等价的变分原理,它完全反映了有关初值问题的全部特征。GD法（General Differential Method）是从泰勒展开式出发,推出的一种求解偏微分方程的数值方法,本文系统地介绍了GD法的基本原理,以及权系数的推导。本文通过卷积将矩形薄板原始控...     

考虑跃迁的指数型炸药空爆荷载等效静载动力系数

为对比抗爆设计规范采用的线性荷载计算模式,建立了考虑跃迁的指数型衰减荷载表达式,通过爆炸荷载等效单自由度微分方程,求解了关于跃迁时长、超压峰值、指数型形状调整参数、结构自振频率与荷载作用时长的等效静载抗力动力系数表达式。根据跃迁时长与形状调整参数,分析了四种典型计算工况,结果表明:现行结构抗爆设计规范等冲量线性衰减荷载可设计范围明显偏小,动力系数在延性比β<3.0下偏保守,而β≥3.0,wt_+>1.4δ时偏不安全,最大偏低17.4%;跃迁时长比值越大,动力系数越大,跃迁时长比为1%～2%时,对动力系数影响程度为0.4～0.7%,指数型荷载形状调整参数对柔度特别大的结构动力系数无影响,对其它结构动力系数增大或减少影响程度不一。     

An analytical model to predict residual thermal stress in 2D orthogonal plain weave fabric composites

This paper seeks to address a practical rectangular truss model to predict residual thermal stress in a 2 D plain weave fabric (PWF) composite. The two orthogonal yarns in a micromechanical unit cell are idealized as straight rods subjected to tensile or compression loading resulting in extension or shortening deformation. The residual thermal stresses and equivalent thermal expansion coefficients in a PWF layer are derived from the thermal constitutive equations and the deformation compatibility condition. Based on the deformation compatibility equations, the thermal constitutive relationships for PWF composites are obtained to derive the residual thermal stresses between PWF plies and pure resin. In order to validate the model, experiments have been performed to investigate the mechanical properties of two-dimensional (2D) orthogonal EW220/5284 PWF composites fabricated by resin transfer moulding (RTM). It is shown that the experimental results correlate well with predictions from the new model.     

Comparison of Two Mathematical Models for 3D Groundwater Flow: Block-Centered Heads and Edge-Based Stream Functions

Traditionally, groundwater flow models, as well as oil reservoir models, are based on the block-centered finite difference method. Well-known models based on this approach are MODFLOW (groundwater) and ECLIPSE (oil and gas). Such models are well proven and robust; their underlying principles are well understood by hydrologists and petroleum reservoir engineers. Nevertheless, the desire to improve the block-centered finite difference paradigm has always been alive, for instance, to be able to apply deformed grid blocks, or to model anisotropy that is not aligned along the coordinate axes. This article introduces the edge-based stream function as a potential alternative to the paradigmatic model, not only to mitigate the above mentioned limitations, but especially for its promise to inverse modeling. Computer programs have been developed for the discrete analog equations of the stream function method and the conventional method. The two methods are tested by using synthetic forward modeling problems of uniform and radial flow. The theoretical formulation and the numerical results show that the two methods are algebraically equivalent and yield the same flux output. However, for rectangular grid blocks and anisotropy aligned along the coordinate axes, the block-centered method is shown to be computationally more efficient than the edge-based stream function method. The major advantage of the stream function method is that it is linear in the resistivities, proving it an ideal candidate for direct inverse modeling. Moreover, any arbitrary specification of stream functions yields a solution that satisfies the mass balance.     

The Edge-Based Face Element Method for 3D-Stream Function and Flux Calculations in Porous Media Flow

We present a velocity-oriented discrete analog of the partial differential equations governing porous media flow: the edge-based face element method. Conventional finite element techniques calculate pressures in the nodes of the grid. However, such methods do not satisfy the requirement of flux continuity at the faces. In contrast, the edge-based method calculates vector potentials along the edges, leading to continuity of fluxes. The method is algebraically equivalent with the popular block-centered finite difference method and with the mixed-hybrid finite element method, but is algorithmically different and has the same robustness as the more conventional node-based velocity-oriented method. The numerical examples are computed analytically and may, therefore, be considered as an 'heuristic proof' of the theory and its practical applicability for reservoir engineering and geohydrology.     

含间隙裂缝的RC结构滞回系统混沌振动分析

研究了含间隙裂缝的钢筋混凝土结构对称滞回非线性问题. 建立了 一种分段线性的对称滞回模型, 利用一次谐波线性方法求解结构系统的等效阻尼和 等效刚度系数,得到了对称滞回非线性系统的等价线性方程. 通过数值分析比较了考虑和不考虑间隙与碰撞 影响的两种情况下系统的混沌动力特性,研究表明: 不考虑间隙与碰撞影响的系统出现周期运 动, 考虑间隙与碰撞影响系统更容易出现混沌运动; 在特定的参数范围内系统一定会出现无 序的混沌运动.     

AN APPROXIMATE MODEL OF UNSTEADY AERODYNAMICS FOR HYPERSONIC PROBLEMS AT HIGH ALTITUDE

针对高空高马赫数飞行环境和强黏性干扰的物理特性, 在当地流活塞理论的基础上引入有效外形修正, 发展了黏性修正当地流活塞理论, 结合定常N-S方程解给出了高空高马赫数下针对该方法的有效外形的判据, 并通过数值算例对该判据进行了验证.通过对典型尖头薄翼和典型钝头翼的一系列二维非定常算例, 将该方法与一阶活塞理论、基于欧拉(Euler)方程的当地流活塞理论和非定常N-S方程数值解进行了对比. 结果显示在高度为40~70 km、马赫数为10~20范围内, 通过该方法计算得到的非定常气动力与非定常N-S方程数值解吻合较好, 明显优于活塞理论和基于Euler方程的当地流活塞理论.该方法克服了传统的活塞理论和当地流活塞理论不能用于高空高马赫数这类强黏性效应情况的弊端, 在较宽的马赫数、攻角、飞行高度范围内都有良好的适用性, 同时其计算效率远高于非定常N-S方程.     

环形桁架结构径向振动的等效圆环模型

在大型环形网架式可展天线中,环形桁架结构的动力学性能对于整个天线的工作状态至关重要.针对大型空间桁架结构,基于连续体等效的思想,将其动力学模型简化为简单的弹性连续体模型一直是动力学研究的热点.将环形桁架结构看作由重复的平面桁架单元构成的环形周期结构,在周期桁架单元等效梁模型的基础上,提出采用不计剪切变形和转动惯量的等效圆环模型分析环形桁架结构的径向振动,并对等效圆环模型的偏微分运动方程进行了解析求解.首先通过变量代换将描述圆环径向振动的四阶偏微分方程组降阶为一阶偏微分方程组,然后通过对降阶后的偏微分方程组进行Laplace变换将其转化为常微分方程组,并采用微分方程组的Green函数解法,获得了等效圆环模型在复频域下动力响应的解析表达式,进而得到等效圆环模型固有振动的特征方程及传递函数的表达式.最后通过数值算例对环形桁架有限元模型与等效圆环模型的固有频率和振型以及传递函数进行了对比分析,验证了等效圆环模型用于环形桁架结构径向振动分析的可行性.