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.     

Biostratigraphic and aminostratigraphic constraints on the age of the Middle Pleistocene glacial succession in north Norfolk,UK

Considerable debate surrounds the age of the Middle Pleistocene glacial succession in East Anglia following some recent stratigraphical reinterpretations. Resolution of the stratigraphy here is important since it not only concerns the glacial history of the region but also has a bearing on our understanding of the earliest human occupation of north‐western Europe. The orthodox consensus that all the tills were emplaced during the Anglian (Marine Isotope Stage (MIS) 12) has recently been challenged by a view assigning each major till to a different glacial stage, before, during and after MIS 12. Between Trimingham and Sidestrand on the north Norfolk coast, datable organic sediments occur immediately below and above the glacial succession. The oldest glacial deposit (Happisburgh Till) directly overlies the ‘Sidestrand Unio‐bed’, here defined as the Sidestrand Hall Member of the Cromer Forest‐bed Formation. Dating of these sediments therefore has a bearing on the maximum age of the glacial sequence. This paper reviews the palaeobotany and describes the faunal assemblages recovered from the Sidestrand Unio‐bed, which accumulated in a fluvial environment in a fully temperate climate with regional deciduous woodland. There are indications from the ostracods for weakly brackish conditions. Significant differences are apparent between the Sidestrand assemblages and those from West Runton, the type site of the Cromerian Stage. These differences do not result from contrasting facies or taphonomy but reflect warmer palaeotemperatures at Sidestrand and a much younger age. This conclusion is suggested by the higher proportion of thermophiles at Sidestrand and the occurrence of a water vole with unrooted molars (Arvicola) rather than its ancestor Mimomys savini with rooted molars. Amino acid racemisation data also indicate that Sidestrand is significantly younger than West Runton. These data further highlight the stratigraphical complexity of the ‘Cromerian Complex’ and support the conventional view that the Happisburgh Till was emplaced during the Anglian rather than the recently advanced view that it dates from MIS 16. Moreover, new evidence from the Trimingham lake bed (Sidestrand Cliff Formation) above the youngest glacial outwash sediments (Briton's Lane Formation) indicates that they also accumulated during a Middle Pleistocene interglacial – probably MIS 11. All of this evidence is consistent with a short chronology placing the glacial deposits within MIS 12, rather than invoking multiple episodes of glaciation envisaged in the ‘new glacial stratigraphy’ during MIS 16, 12, 10 and 6. Copyright © 2009 John Wiley & Sons, Ltd.     

Statistical analysis of sand grain/bed collision process recorded by high-speed digital camera

A high‐speed digital camera was employed to record the sand grain/bed collision process. With image processing and a statistical method, a series of parameters of the collision process were obtained. The results show that the collision process of a grain with rebounding can be represented by two parameters: the kinetic energy restitution coefficient and the collision angle. Both parameters satisfy a normal distribution, and they are dependent on one another. With an increase of the collision angle, the distribution of the kinetic energy restitution gradually reduces from a broad to a narrow range with low values. The percentage of vertical velocity restitution coefficients greater than 1 can reach 70% or more, which ensures that the settling time of the sand grains in the air increases and that they receive more energy from the air to progress the saltation movement.     

河口区治江围涂对杭州湾水动力及海床影响分析

自二十世纪六十年代起，钱塘江河口进行了大规模治江围涂，人类活动对杭州湾水域的影响越来越为人们关注。基于实测资料的分析和潮汐水流数学模型计算，探讨了钱塘江河口区围涂对杭州湾水域潮汐水流和海床的影响。研究成果可为钱塘江河口和杭州湾的进一步开发治理和环境保护提供科学依据。     

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

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

纳米碳酸钙材料的工业合成与应用

该文对目前工业上生产纳米碳酸钙的方法作了较为详细的综述。对一步碳化法、两步碳化法、多段喷雾碳化法、旋转填充床碳化反应器碳化技术制备链形、纺锤形、球形、立方形等不同形状的纳米碳酸钙材料的方法进行了比较和总结     

A General Linear Wave Theory for Water Waves Propagating over Uneven Porous Bottoms

Starting from the widespread phenomena of porous bottoms in the near shore region, considering fully the diversity of bottom topography and wave number variation, and including the effect of evanescent modes, a general linear wave theory for water waves propagating over uneven porous bottoms in the near shore region is established by use of Green‘s scond identity. This theory can be reduced to a number of the most typical mild-slope equations curreutly in use and provide a reliable research basis for follow-up development of nonlinear water wave theory involving porous bottoms.     

半潜式超大浮体连接器动力特性的一种时间序列分析方法

根据刚性模块弹性连接器 (RMFC)假设研究由 3个模块连接而成的半潜式超大浮体MobileOffshoreBase(MOB)的波浪载荷动力响应。模块和连接器动力响应根据时间序列法计算 ,这种方法包括海况模拟 ,波浪载荷计算和结构响应计算三个主要部分。通过对结构的适当简化 ,根据Airy波理论 ,Morrison公式以及流体中的结构运动方程、弹性连接器变形方程 ,计算求得MOB在不规则海况下的动力响应 ,并将其与相关文献中的结果进行比较 ,表明该方法有良好的精确性