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

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

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.     

About Darcy's law in non-Galilean frame

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.     

声波在水-多孔介质海底界面上的反射与透射

声波在海底界面的反射和透射是海底散射、海底混响、海底目标探测的重要问题。利用Biot多孔介质声传播理论对声波在水-多孔海底界面上的反射和透射进行了分析,具体给出反射声波的反射系数,3种透射声波的透射系数以及声能透射系数随入射波入射角和频率(10~40 kHz)的变化关系,分析了各种透射波对透射声能的贡献。多孔海底介质参数分别采用Stoll和Chotiros给出的2组参数进行理论计算。     

现代黄河三角洲土地利用变化分析

人类对海岸带的开发利用影响着海岸的变迁。土地利用和土地覆被变化可以说是自然与人文过程交叉最为密切的问题，是人地关系的具体体现。通过对黄河三角洲土地利用变化的研究，可了解人为因素对三角洲海岸带的开发利用状况及其对海岸侵蚀的影响。根据1987年和2001年两期卫星遥感图像和2001年野外调查数据，运用地理信息系统方法及数理统计方法，分析了现代黄河三角洲地区1987年至2001年土地利用／覆被的数量变化和空间变化特征；揭示了该地区各类土地利用数量变化的幅度、速度和土地利用空间变化的主要类型；阐明了该区各类土地利用变化的区域特点。该地区土地利用变化速度比整个环渤海地区速度快，草地、滩地和不明地类这3种类型大量减少，其他类型增多。从变化率来看，养殖区、库塘、油田用地增长最快。库塘、海滩养殖区等类型的增加，滨海草地的减少，使海岸侵蚀加剧。     

波浪作用下海床动力反应的数值分析

近海和离岸建筑物的基通常处于连续不断的小风浪作用之下，可将其变形视为弹性。基于二维广义Biot理论，提出了线性或非线性波浪作用下饱和弹性海床动力应应的时域有限元数值解法。静力平衡条件和Biot方程组成的边值方程可视为其特例，在比较算例中，数值计算得到的孔压和有效应力幅值沿海床深度的分布与解析解十分吻合。土骨架和孔隙流体的加速度对海床动力反应的影响很小。具体算例表明，线性波沿缓坡海床传播时，土层中超静孔压和有效应力幅值随之增大，有可能发生滑动坡坏。所提出的数值解法能够灵活地处理非线性波浪荷载，海底复杂地质条件和波浪沿缓坡传播等复杂情况。     

Modelling linear wave transformation induced by dissipative structures—Random waves

The relevant theory is presented and numerical results are compared with the analytical solution for the interaction of non-breaking waves with an array of vertical porous circular cylinders on a horizontal bed. The extension to the cases of unidirectional and multidirectional waves is obtained by means of a transfer function. The influence of the mechanical properties of porous structures and wave irregularity on wave transformation is analysed. Results for unidirectional and multidirectional wave spectra are compared to those obtained for regular waves. The model presented reproduces well the analytical results and provides a tool for analysing several engineering problems.     

Laboratory Culture and Preliminary Characterization of the Nitrogen-Fixing Rhizosolenia-Richelia Symbiosis

Abstract. The oceanic diatom Rhizosolenia clevei and its cyanobactcrial symbiont Richelia intracel-luluris were isolated into laboratory culture and their biology and growth dynamics examined. Nitrogen-fixation by the symbiont could sustain the symbiosis. Growth and nitrogenase activity demonstrated light saturation kinetics, with no photoinhibition up to 315 μE m^-1 s^-1 for growth and 780μE m^-1 s^-1 for nitrogen fixation. The symbiosis is not obligate for Rhizosolenia , which is capable of growth independent of the symbiont if a nitrogen source is available. The symbiont is contained in the pcriplasmic space between the Rhizosolenia plasmalcmma and frustulc, and preliminary evidence suggests excretion of fixed nitrogen into the medium may be occurring.     

黄河尾闾河道1996年改道的意义及黄河三角洲演化趋势

现代黄河三角洲由几个亚三角洲组成。 1 85 5 -1 93 4年 ,形成了以宁海为顶点的第一亚三角洲 ;1 93 4-1 996年 ,三角洲冲积扇顶点进一步下移 ,形成了以渔洼为顶点的第二亚三角洲。最近的一次改道是在 1 996年 6月 ,改道后的河流沿清水沟流路的清 8剖面附近东流入海。这次改道可能是第三亚三角洲形成的开始。以后数十年内的黄河尾闾摆动将在清 8剖面附近进行。虽然人类活动将对三角洲的发展演化产生越来越明显的影响 ,但这些亚三角洲的发育演化仍有明显的规律性