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.  相似文献   

西安地区单桩桩土相互作用数值模拟分析   总被引：1，自引：0，他引：1  

魏静  王建华  李永林 《地球科学与环境学报》2003,25(3):63-66

随着计算机技术和数值方法的不断发展和完善,使有限元法得到越来越多的应用。根据西安地区原型摩擦桩的地层特点和土性条件,运用有限元基础理论,利用相应的边界条件与假设条件,建立单桩数学模型。通过模型研究了单桩的沉降特性与土共同作用的影响因素,对西安地区单桩桩土相互作用机理进行了探讨。  相似文献   

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.  相似文献   

Prediction of local scour around circular piles under waves using a novel artificial intelligence approach     

Iman Ahmadianfar  Xuefeng Chu 《Marine Georesources & Geotechnology》2021,39(1):44-55

Abstract

The scour phenomena around vertical piles in oceans and under waves may influence the structure stability. Therefore, accurately predicting the scour depth is an important task in the design of piles. Empirical approaches often do not provide the required accuracy compared with data mining methods for modeling such complex processes. The main objective of this study is to develop three data-driven methods, locally weighted linear regression (LWLR), support vector machine (SVR), and multivariate linear regression (MLR) to predict the scour depth around vertical piles due to waves in a sand bed. It is the first effort to develop the LWLR to predict scour depth around vertical piles. The models simulate the scour depth mainly based on Shields parameter, pile Reynolds number, grain Reynolds number, Keulegan–Carpenter number, and sediment number. 111 laboratory datasets, derived from several experimental studies, were used for the modeling. The results indicated that the LWLR provided highly accurate predictions of the scour depths around piles (R?=?0.939 and RMSE = 0.075). Overall, this study demonstrated that the LWLR can be used as a valuable tool to predict the wave-induced scour around piles.  相似文献   

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

彭临慧  赵燕鹏  郁高坤 《中国海洋大学学报(自然科学版)》2007,37(4):671-675

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

港口靠船桩工作性状计算的双参数法   总被引：8，自引：0，他引：8  

谢耀峰 《海洋工程》2002,20(2):38-42,48

港口靠船桩是承受重复性荷载或多循环荷载的一种横向承载桩。双参数法能较好地描述靠船桩 (单桩、群桩 )的桩土工作性状。用双参数法进行计算时 ,所给的地面条件必须由重复荷载或多循环荷载作用下实测出来。根据现场试桩资料 ,标定桩土参数 ,即k =mx1/n的指数l/n、桩土相对柔度系数α、桩土综合刚度EI等值 ,可以用来设计试桩附近局部地区同类条件下的其它长桩。对于没有进行试桩的某些工程 ,可以联合运用双参数法和p -y曲线法进行桩的分析  相似文献   

A General Linear Wave Theory for Water Waves Propagating over Uneven Porous Bottoms   总被引：1，自引：0，他引：1  

锁要红  黄虎 《中国海洋工程》2004,18(1):163-171

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.  相似文献   

Experimental study of strong reflection of regular water waves over submerged breakwaters in tandem   总被引：3，自引：0，他引：3  

Yong-Sik Cho  Jong-In Lee  Young-Taek Kim   《Ocean Engineering》2004,31(10):88-1335

A series of laboratory experiments was carried out to investigate the strong reflection of regular water waves over a train of submerged breakwaters. Rectangular and trapezoidal shapes of submerged breakwaters are employed and compared for reflecting capability of incident waves. Measured reflection coefficients of regular waves over impermeable submerged breakwaters are verified by comparing with those of the eigenfunction expansion method. A very good agreement is observed. Reflection coefficients of permeable submerged breakwaters are less than those of impermeable breakwaters. The trapezoidal shape is recommended for a submerged breakwater in terms of reflecting capability and practical application.  相似文献   

Wave interaction with a perforated wall breakwater with a submerged horizontal porous plate   总被引：1，自引：0，他引：1  

Yong Liu  Yu-cheng Li  Bin Teng 《Ocean Engineering》2007,34(17-18):2364-2373

This study examines the hydrodynamic performance of a new perforated-wall breakwater. The breakwater consists of a perforated front wall, a solid back wall and a submerged horizontal porous plate installed between them. The horizontal porous plate enhances the stability and wave-absorbing capacity of the structure. An analytical solution based on linear potential theory is developed for the interaction of water waves with the new proposed breakwater. According to the division of the structure, the whole fluid domain is divided into three sub-domains, and the velocity potential in each domain is obtained using the matched eigenfunction method. Then the reflection coefficient and the wave forces and moments on the perforated front wall and the submerged horizontal porous plate are calculated. The numerical results obtained for limiting cases are exactly the same as previous predictions for a perforated-wall breakwater with a submerged horizontal solid plate [Yip, T.L., Chwang, A.T., 2000. Perforated wall breakwater with internal horiontal plate. Journal of Engineering Mechanics ASCE 126 (5), 533–538] and a vertical wall with a submerged horizontal porous plate [Wu, J.H., Wan, Z.P., Fang, Y., 1998. Wave reflection by a vertical wall with a horizontal submerged porous plate. Ocean Engineering 25 (9), 767–779]. Numerical results show that with suitable geometric porosity of the front wall and horizontal plate, the reflection coefficient will be always rather small if the relative wave absorbing chamber width (distance between the front and back walls versus incident wavelength) exceeds a certain small value. In addition, the wave force and moment on the horizontal plate decrease significantly with the increase of the plate porosity.  相似文献