Three dimensional simulation of fluid flow in X‐ray CT images of porous media     

A. Al‐Omari  E. Masad 《国际地质力学数值与分析法杂志》2004,28(13):1327-1360

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

On The Concept Of Periapsis In Hill’s Problem     

B.F. Villac  D.J. Scheeres 《Celestial Mechanics and Dynamical Astronomy》2004,90(1-2):165-178

The concept of closest approach is analyzed in Hill’s problem, resulting in a partitioning of the position space. The different behavior between the direct and retrograde motion is explained analytically, resulting in a simple estimate of the variation of Hill’s periodic and quasi-circular orbits as a function of the Jacobi constant. The local behavior of the orbits on the zero velocity surfaces and an analytical definition of local escape and capture in Hill’s problem are also given.  相似文献   

About Darcy's law in non-Galilean frame     

C. Geindreau  E. Sawicki  J.-L. Auriault  P. Royer 《国际地质力学数值与分析法杂志》2004,28(3):229-249

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

Structures in the phase space of a four dimensional symplectic map     

Ch. Skokos  G. Contopoulos  C. Polymilis 《Celestial Mechanics and Dynamical Astronomy》1996,65(3):223-251

We numerically investigate the projections of non periodic orbits in a 4-dimensional (4-D) symplectic map composed of two coupled 2-dimensional (2-D) maps. We describe in detail the structures that are produced in different planes of projection and we find how the morphology of the 4-D orbits is influenced by the features of the 2-D maps as the coupling parameter increases. We give an empirical law that describes this influence.  相似文献   

Minor planet short periodic perturbations: The indirect part of the disturbing function     

Zoran Knežević 《Celestial Mechanics and Dynamical Astronomy》1993,55(4):387-404

Leverrier's development of the indirect part of the disturbing function has been extended to include terms up to degree 4 in eccentricity and inclination; the resulting series has been expressed with respect to a fixed plane, and in a computer readable form (a list of integers). Tests have been performed for the relative significance of the terms of degrees 2, 3 and 4, and estimates have been obtained for the accuracy of the short periodic perturbations of a minor planet, and of the corresponding mean orbital elements. It was found that: (i) even in extreme cases, the indirect part of the disturbing function gives rise to very small short periodic perturbations; (ii) bodies of very high eccentricity/inclination and those close to mean motion resonances are most significantly affected; (iii) indirect perturbations for minor planets can be computed up to the degree 2 terms only, without any significant loss of accuracy; and (iv) higher degree indirect perturbations appear to be important only for their contribution to the long periodic effects of higher order (with respect to the perturbing mass).  相似文献   

The Homoclinic Tangle of A 1:2 Resonance in a 2-D Hamiltonian System     

C. Polymilis  G. Contopoulos  A. Dokoumetzidis 《Celestial Mechanics and Dynamical Astronomy》2003,85(2):105-144

We study in great detail the geometry of the homoclinic tangle, with respect to the energy, corresponding to an unstable periodic orbit of type 1:2, on a surface of section representing a 2-D Hamiltonian system. The tangle consists of two resonance areas, in contrast with the tangles of type-l or -{l, m, k, x = 0} considered in previous studies, that consist of only one resonance area. We study the intersections of the inner and outer lobes of the same resonance area and of the two resonance areas. The intersections of the lobes follow certain rules. The detailed study of these rules allows us to derive quantitative relations about the number of intersections and to understand the complex behavior of the higher order lobes by studying the lower order lobes. We find 1st, 2nd, 3rd, etc. order intersections formed by lobes making 1, 2, 3, etc. turns around an island. After a sufficiently high order of iterations a lobe may intersect its image and thus produce a Poincaré recurrence. Numerical results for a wide interval of energies are presented. The number of intersections changes through tangencies. In any finite interval of the energy between two tangencies of 1st order, an infinite number of higher order tangencies occur and thus, according to the Newhouse theorem, there exist nearby islands of stability.  相似文献   

Resonant Periodic Orbits of Trans-Neptunian Objects     

Kotoulas  Thomas A.  Hadjidemetriou  John D. 《Earth, Moon, and Planets》2002,91(2):63-93

We study two and three-dimensional resonant periodic orbits, usingthe model of the restricted three-body problem with the Sun andNeptune as primaries. The position and the stability character ofthe periodic orbits determine the structure of the phase space andthis will provide useful information on the stability and longterm evolution of trans-Neptunian objects. The circular planarmodel is used as the starting point. Families of periodic orbitsare computed at the exterior resonances 1/2, 2/3 and 3/4 withNeptune and these are used as a guide to select the energy levelsfor the computation of the Poincaré maps, so that all basicresonances are included in the study. Using the circular planarmodel as the basic model, we extend our study to more realisticmodels by considering an elliptic orbit of Neptune and introducingthe inclination of the orbit. Families of symmetric periodicorbits of the planar elliptic restricted three-body problem andthe three-dimensional problem are found. All these orbitsbifurcate from the families of periodic orbits of the planarcircular problem. The stability of all orbits is studied. Althoughthe resonant structure in the circular problem is similar for allresonances, the situation changes if the eccentricity of Neptuneor the inclination of the orbit is taken into account. All theseresults are combined to explain why in some resonances there aremany bodies and other resonances are empty.  相似文献   

用随机介质模型方法描述孔洞型油气储层   总被引：5，自引：0，他引：5  

朱生旺  魏修成  曲寿利  刘春园  吴开龙 《地质学报》2008,82(3):420-427

孔洞型储集层具有强烈的非均质性,在空间分布上变化剧烈,规律性差,其地球物理参数(如速度、密度、弹性参数等)在空间上的变化很难用传统的层状介质模型进行刻画.本文借助随机模型方法,研究孔洞型储层中不同孔洞尺度、分布密度与模型参数的关系,建立反映实际孔洞油气储层介质地球物理参数空间分布统计特征的随机介质模型,以更方便地对孔洞型油气储层进行描述,并为碳酸盐岩孔洞型储层地震反射特征的研究提供有效的方法途径.  相似文献   

珠江三角洲网河区水位变化趋势研究   总被引：4，自引：0，他引：4  

沈焕庭  杨清书  罗宪林  罗章仁  杨干然  欧素英 《海洋学报》2002,24(2):30-38

根据珠江三角洲网河区29个验潮站的实测验潮记录,应用傅氏变换与最平滤波器串联的方法来消除月均序列的周期波动对确定水位变化趋势的影响,由低通序列一元线性回归分析确定各站水位的变化趋势;结果表明,周期波动对确定海平面变化趋势的影响是显著的.应用经验正交函数(EOF)对网河区的水位变化场进行分解,由相互独立的时间函数和空间特征函数表征网河区区域的水位变化特征;应用时间特征函数计算区域水位的平均变化率为0.02mm/a.根据验潮站的水位变化趋势,探讨网河区水位变化与河床冲淤的关系.  相似文献