Face-Centered and Volume-Centered Discrete Analogs of the Exterior Differential Equations Governing Porous Medium Flow I: Theory

This paper presents a physics-oriented approach to approximate the continuum equations governing porous media flow by discrete analogs. To that end, the continuity equation and Darcys law are reformulated using exterior differential forms. This way the derivation of a system of algebraic equations (the discrete analog) on a finite-volume mesh can be accomplished by simple and elegant translation rules. In the discrete analog the information about the conductivities of the porous medium and the metric of the mesh are represented in one matrix: the discrete dual. The discrete dual of the block-centered finite difference method is presented first. Since this method has limited applicability with respect to anisotropy and non-rectangular grid blocks, the finite element dual is introduced as an alternative. Application of a domain decomposition technique yields the face-centered finite element method. Since calculations based on pressures in volume centers are sometimes preferable, a volume-centered approximation of the face-centered approximation is presented too.     

非饱和多孔介质有限元分析的基本控制方程与变分原理

本文在对问题研究现状进行阐述的基础上较系统地给出了骨架可变形非饱和多孔介质的全耦合分析模型，模型中考虑了孔隙气体，水（油）流动对介质力学性能的影响，多孔介质的饱和度，渗透系数与毛吸压力的关系，由实验给出，所导出的控制方程以固体骨架的位移与孔隙流体压力为基本未知量，由于问题的非自共轭特征，文中构造了非饱和介质动力问题的参数变分形式，并在此基础上给出有限元离散方程。     

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.     

Inertial and Second Grade Effects on Macroscopic Equations for the Flow of a Fluid in a Porous Medium

Homogenization techniques are used to upscale from pore to laboratory or field scale viscous and second grade nonNewtonian flow in a porous medium. Nonlinear forms of Darcy's law are obtained and analysed under a series of symmetry properties. The general case of displacement of one of these fluids by another with different properties is considered and a linear stability analysis is performed.     

On the Trajectory Method for the Reconstruction of Differential Equations from Time Series

This work investigates the trajectory method [1] for thereconstruction of ordinary differential equations (ODEs) from timeseries. The potentials of the method are analyzed for dynamical systemsdescribed by second- and third-order ODEs, focusing in particular on therole of the parameters of the method and on the influence of the qualityof the time series in terms of noise, length and sampling frequency.Typical models are investigated, such as the van der Pol, the linearmechanical, the Duffing and the Rössler equations, resulting in arobust and versatile method which is capable of allowing interestingapplications to experimental cases. The method is then applied to themeasured time series of a nonlinear mechanical oscillator, a typicalvelocity oscillation of the bursting phenomenon in near-wall turbulenceand the averaged annual evolution of rainfall, temperature andstreamflow over a hydrological basin.     

Chaos in the Convective Flow of a Fluid Mixture in a Porous Medium

The results of the study of the global behaviour of the convective flow of a binary mixture in a porous medium are presented. Bifurcation diagram, fixed points, periodic, chaotic solutions, stable and unstable manifolds, and basins of attraction have been calculated. Different behaviours (chaos, undecidable behaviour, etc.) have been found.     

Wave Dynamics of Saturated Porous Media and Evolutionary Equations

Nonlinear wave dynamics of an elastically deformed saturated porous media is investigated following the Biot approach. Mathematical models under research are the Biot model and its generalization by consideration of viscous stresses inside liquids. Using two-scales and linear WKB methods, the classical Biot system is transformed to a first-order wave equation. To construct the solution of the other system, an asymptotic modified two-scales method is developed. Initial system of equations is transformed to a nonlinear generalized Korteweg–de Vries–Burgers equation for quick elastic wave. Distinctions of wave propagation in the context of the Biot model and its generalization are shown.     

Remarks on the Perturbation Methods in Solving the Second-Order Delay Differential Equations

The paper presents a study on the validity of perturbation methods, suchas the method of multiple scales, the Lindstedt–Poincaré method and soon, in seeking for the periodic motions of the delayed dynamic systemsthrough an example of a Duffing oscillator with delayed velocityfeedback. An important observation in the paper is that the method ofmultiple scales, which has been widely used in nonlinear dynamics, worksonly for the approximate solutions of the first two orders, and givesrise to a paradox for the third-order approximate solutions of delaydifferential equations. The same problem appears when theLindstedt–Poincaré method is implemented to find the third-orderapproximation of periodic solutions for delay differential equations,though it is effective in seeking for any order approximation ofperiodic solutions for nonlinear ordinary differential equations. Apossible explanation to the paradox is given by the results obtained byusing the method of harmonic balance. The paper also indicates thatthese perturbation methods, despite of some shortcomings, are stilleffective in analyzing the dynamics of a delayed dynamic system sincethe approximate solutions of the first two orders already enable one togain an insight into the primary dynamics of the system.     

Investigation of the Relative Motion of a Viscous Fluid and a Porous Medium Using the Brinkman Equations

Exact solutions are obtained for the following three problems in which the Brinkman filtration equations are used: laminar fluid flow between parallel plane walls, one of which is rigid while the other is a plane layer of saturated porous medium, motion of a plane porous layer between parallel layers of viscous fluid, and laminar fluid flow in a cylindrical channel bounded by an annular porous layer.     

Exponential Dichotomies,the Fredholm Alternative,and Transverse Homoclinic Orbits in Partial Functional Differential Equations

In this paper we investigate the transversality of homoclinic orbits in partial functional differential equations. We first discuss the exponential dichotomies for linear operator equations. Then we show that the Fredholm Alternative holds if the homogeneous equation has exponential dichotomies on R. Transversality of homoclinic orbits for periodically perturbed partial functional differential equations is studied using the Liapunov-Schmidt method and the Melnikov integral. Ams Subject Classifications: 35R10; 58F14.     

The Conjugacy of Stochastic and Random Differential Equations and the Existence of Global Attractors

We consider stochastic differential equations in d-dimensional Euclidean space driven by an m-dimensional Wiener process, determined by the drift vector field f₀ and the diffusion vector fields f₁,...,f_m, and investigate the existence of global random attractors for the associated flows . For this purpose is decomposed into a stationary diffeomorphism given by the stochastic differential equation on the space of smooth flows on R^d driven by m independent stationary Ornstein Uhlenbeck processes z¹,...,z^m and the vector fields f₁,...,f_m, and a flow generated by the nonautonomous ordinary differential equation given by the vector field (_t/x)^–1[f₀(_t)+ _i=1 ¹ f_i(_t)z _t ⁱ ]. In this setting, attractors of are canonically related with attractors of . For , the problem of existence of attractors is then considered as a perturbation problem. Conditions on the vector fields are derived under which a Lyapunov function for the deterministic differential equation determined by the vector field f₀ is still a Lyapunov function for , yielding an attractor this way. The criterion is finally tested in various prominent examples.     

砌体结构破坏倒塌过程的离散单元模拟

基于传统颗粒离散单元模型（DEM）,开发了一种绑定式离散单元模型（BTDEM）用于砌体结构崩塌机制的研究,并对一个典型砌体结构的破坏倒塌过程进行数值模拟。结果证实BTDEM方法能够详细地模拟出砌体建筑结构的崩塌过程,同时还可以模拟出砌块的破碎,是一种普适性较强的模拟砌体结构建筑的高效数值方法。BTDEM方法具有计算速度...     

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.     

Derivation and Consideration of Governing Equations of the Envelope Surface Created by Directional,Nearly Monochromatic Waves

In this paper we expand frequency-based nearly monochromatic waves to propagation-direction-based nearly monochromatic waves, and givesthe new concept of directional, nearly monochromaticwaves, whose energy is almost concentrated in a singlepropagation direction. We show directional,nearly monochromatic waves also create the envelopesurface, which is featured by time-invariant in alinear system. The fact that a non-linear systemmakes the envelope surface time-variant is shown.     

Lyapunov Stability and Adaptive Regulation of a Class of Nonlinear Nonautonomous Second-Order Differential Equations

This paper presents an adaptive regulation scheme for a class ofordinary nonlinear nonautonomous second-order differential equationswhich includes as particular cases a number of particular differentialequations which occur in applications. The unforced reference model isproposed to be a stable differential parametrization within the generalclass dealt with. Therefore, some sufficient Lyapunov's stabilityconditions for such a class are previously investigated which can beused, in particular, to set an appropriate reference model. Theresulting closed-loop adaptive scheme is proved to be stable and itinvolves a parameter estimation scheme of least-squares type which isproved to possess all suitable properties in terms of estimatesboundedness and asymptotic convergence of the estimates to finite limitsas well as time-integrability of the squared adaptation error.     

Sensitivity Analysis and Quantification of the Role of Governing Transport Mechanisms and Parameters in a Gas Flow Model for Low-Permeability Porous Media

Recent models represent gas (methane) migration in low-permeability media as a weighted sum of various contributions, each associated with a given flow regime. These models typically embed numerous chemical/physical parameters that cannot be easily and unambiguously evaluated via experimental investigations. In this context, modern sensitivity analysis techniques enable us to diagnose the behavior of a given model through the quantification of the importance and role of model input uncertainties with respect to a target model output. Here, we rely on two global sensitivity analysis approaches and metrics (i.e., variance-based Sobol’ indices and moment-based AMA indices) to assess the behavior of a recent interpretive model that conceptualizes gas migration as the sum of a surface diffusion mechanism and two weighted bulk flow components. We quantitatively investigate the impact of (i) each uncertain model parameter and (ii) the type of their associated probability distribution on the evaluation of methane flow. We then derive the structure of an effective diffusion coefficient embedding all complex mechanisms of the model considered and allowing quantification of the relative contribution of each flow mechanism to the overall gas flow.

     

Generalized Macroscopic Models for Fluid Flow in Deformable Porous Media II: Numerical Results and Applications

In Part I of this study, generalized mathematical models were developed to describe the motion of fluids in porous media. The second part of this study solved the problem of fluid flow in small channels of a periodic elastic solid matrix at the pore scale numerically, and applied the volume-averaging technique to predict the macroscopic behavior of reservoirs. The numerical results demonstrated different macroscopic behavior of a porous medium due to cyclic excitation at various frequencies corresponding to the five separate characteristic macroscopic models identified in Part I. The results emphasize the need to use an appropriate model to interpret the corresponding responses of a saturated porous medium.     

考虑颗粒形状和破碎的胶结钙质砂力学行为离散元模拟研究

质砂作为一种建筑材料,近年来广泛应用于我国南海岛礁工程建设中。本文通过建立考虑钙质砂真实颗粒形状和颗粒破碎的胶结钙质砂离散元模型,研究了二维剪切条件下试样的宏微观力学行为,包括应力-应变行为、颗粒破碎、胶结破坏、位移场和裂纹随剪应变的演化规律,讨论了颗粒形状、颗粒粒径范围、颗粒强度和水泥胶结强度对胶结钙质砂力学行为的影响规律。结果表明,钙质砂颗粒粒径区间越宽,胶结钙质砂的强度越高。同一级配条件下,考虑真实颗粒形状的胶结钙质砂试样比圆颗粒试样的强度更高,试样总体颗粒破碎率也更高。钙质砂颗粒的强度越高,胶结钙质砂的性能越好。但是提高水泥的强度对胶结钙质砂力学性能的影响并不显著。本文的研究结果可为实际工程中钙质砂的加固提供理论依据。     

Effect of the Source Term on Steady Free Convection Boundary Layer Flows over an Vertical Plate in a Porous Medium. Part I

The Darcy free convection boundary layer flow over a vertical flat plate is considered in the presence of volumetric heat generation/absorption. In the present first part of the paper it is assumed that the heat generation/absorption takes place in a self-consistent way, the source term q ^′′′≡ S of the energy equation being an analytical function of the local temperature difference T − T _∞. In a forthcoming second part, the case of the externally controlled source terms S = S(x,y ) will be considered. It is shown that due to the presence of S, the physical equivalence of the up- and downflows gets in general broken, in the sense that the free convection flow over the upward projecting hot plate (“upflow”) and over its downward projecting cold counterpart (“downflow”) in general become physically distinct. The consequences of this circumstance are examined for different forms of S. Several analytical solutions are given. Some of them describe algebraically decaying boundary layers which can also be recovered as limiting cases of exponentially decayingones. This asymptotic phenomenon is discussed in some detail.