通孔金属泡沫渗透率解析模型

针对通孔金属泡沫中的渗透率预测及现有理论模型的局限性,发展了一种新的全解析渗透率模型. 该模型以立方体结构作为代表单元,采用基于追踪流体微团轨迹的分支算法解析求解代表单元内的流动迂曲度. 渗透率的表达形式简单且不含任何拟合或经验参数,仅是孔隙率与平均孔径的函数. 采用实验测量和文献数据对模型预测进行了验证. 结果表明:提出的模型能够在较为宽广的孔隙率(0.55~0.98) 和孔密度(5~100 PPI) 范围内预测孔通孔金属泡沫的渗透率;采用分支算法得到的流动迂曲度能够较好地描述流体在通孔金属泡沫中的流动特征;采用开孔率修正的解析模型亦能对半开孔泡沫材料的渗透率提供良好预测.      

Permeability of open-cell foams under compressive strain

A model for the behavior of low-density, open-cell foam under compressive strain is proposed. Using this model, a tractable relationship between the normalized permeability and the applied strain is developed. An experimental study of the effect of strain on the permeability of open-cell polyurethane foams is presented. The experiments are performed using a Newtonian fluid in the fully laminar regime, where viscous forces are assumed to dominate. The model is found to describe the experimental data well and be independent of the foam cell size, the direction of flow with respect to the foam rise direction, and the properties of the saturating fluid. In a companion paper, the model for the permeability of open-cell foam is combined with Darcy’s law to give the contribution of viscous fluid flow to the stress–strain response of a reticulated foam under dynamic loading.     

不同滑移边界下的页岩渗透率修正模型

页岩中的孔隙直径通常为纳米量级,基于连续流的达西定律已不能描述纳米级孔隙内的气体流动规律,一般采用附加滑移边界条件的Navier-Stokes方程对其进行描述. 由此可导出与压力相关的渗透率公式(称为"视渗透率"),并用来修正达西定律.因而,渗透率修正方法研究成为页岩气流动研究的热点之一.首先,基于Hagen-Poiseuille 流推导出一般形式二阶滑移模型下的速度分布和流量公式,并推导出相应的渗透率修正公式.该渗透率修正公式基本能将现有的滑移速度模型统一表达为对渗透率的修正. 基于一般形式的渗透率修正公式,重点研究了Maxwell, Hsia, Beskok与Ng 滑移模型速度分布渗透率修正系数、及其对井底压力的影响;提出了基于Ng 滑移速度模型的渗透率修正公式. 基于页岩实际储层温压系统及孔隙分布,计算了Kn 范围及储层条件下页岩气的流动形态,表明页岩气流动存在滑移流、过渡流与分子自由流. 而Ng 模型能描述Kn<88 的滑移流、过渡流、自由分子流的流量规律,因此可以用于描述页岩实际储层中页岩气的流动特征. 计算表明,随着Kn 的增加,不同滑移模型下的渗透率修正系数差异增大.Maxwell与Hsia模型适用于滑移流与过渡流早期,Beskok与Ng 模型可描述自由分子流下的流动规律,但二者在虚拟的孔径均为10nm页岩中,井底压力的差别开始显现;在虚拟的孔径均为1nm页岩中,井底压力的差别开始明显.      

A Theoretical Model with Experimental Verification to Predict Hydrodynamics of Foams

An experimentally validated theoretical model, based on hydraulic resistance network and scale analysis at the pore level, is developed to predict the pressure drop for flow through foams. The complex microstructure of the foams is modeled as a matrix of interconnected solid ligaments forming simple cubic arrays of cylinders. New correlations for permeability and form drag (inertia) coefficient are presented as functions of the mean pore and ligament diameter as well as the foam porosity. The present model makes it possible to conduct parametric studies. Results obtained from the proposed model are successfully compared with our experimental data as well those found in the literature to observe good agreement.     

通孔泡沫铝的动态压缩行为

在SHPB装置上对渗流法制备的通孔泡沫铝进行了动态压缩实验,研究了相对密度为0.341~0.419的通孔泡沫铝在10-3~2000 s-1应变率范围内的压缩响应特征和应变率相关性,并用扫描电镜（scanning electron microscope,SEM)分析了泡沫铝的压缩变形特征。实验结果表明,通孔泡沫铝有明显应变率效应,随应变率上升,泡沫铝流动应力提高。SEM观察结果揭示,在动态压缩下,通孔泡沫铝宏观上均匀变形,微观变形机制以泡孔横向伸展坍塌为主。     

Experimental study on the two-phase pressure drop in copper foams

Experiments were performed to study pressure drops in copper foams embedded in a rectangular copper channel. De-ionized water was used as the working fluid with mass fluxes of 30–200 kg/m² s, and inlet temperature of 40–80°C. The copper foam has the porosity of 0.88 and the pore densities of 30, 60 and 90 ppi (pores per inch). Both single-phase liquid flow and boiling two-phase flow are studied. Effects of mass fluxes, vapor mass qualities, and average pore diameters of metallic foams are investigated. It is found that friction factors for the single-phase liquid flow are mainly dependent on the Reynolds number and the average pore diameter of metallic foams. The friction factors are decreased with increases in the Reynolds numbers, and will approach 0.22 at high Reynolds numbers. For the boiling two-phase flow, two-phase pressure drops are increased with increases in the outlet vapor mass qualities, mass fluxes, and ppi values. The two-phase multiplier is increased with increases in the outlet vapor mass qualities and mass fluxes, and it is decreased with increases in the Martinelli parameter and will attain a constant value depending on the mass fluxes. The larger the mass fluxes, the larger the constant value is. An experimental correlation considering the effects of vapor mass qualities, mass fluxes, and average pore diameters of metallic foams is recommended, showing good accuracy to predict the two-phase pressure drops in metallic foams.     

Effect of Compressible Foam Properties on Pressure Amplification During Shock Wave Impact

A comprehensive study is made of the influence of the physical properties of compressible open-cell foam blocks exposed to shock-wave loading, and particularly on the pressure distribution on the shock tube walls. Seven different foams are used, with three different shock Mach numbers, and three different slab lengths. Foam properties examined include permeability, density, stiffness, tortuosity and cell characteristics. The investigations concentrate on both side-wall and back-wall pressures, and the peak pressures achieved, as well as collapse velocities of the front face and the strength and nature of the reflected shock wave. The consequences of deviations from one-dimensionality are identified; primarily those due to wall friction and side-wall leakage. The results presented are the most comprehensive and wide ranging series conducted in a single facility and are thus a significant resource for comparison with theoretical and numerical studies. The different foams show significant differences in behavior, both in terms of peak pressure and duration, depending primarily on their density and permeability.This paper was based on work presented at the 2nd International Symposium on Interdisciplinary Shock Wave Research, Sendai, Japan on March 1–3, 2005.     

Dependency of Tortuosity and Permeability of Porous Media on Directional Distribution of Pore Voids

Single-phase fluid flow in porous media is usually direction dependent owing to the tortuosity associated with the internal structures of materials that exhibit inherent anisotropy. This article presents an approach to determine the tortuosity and permeability of porous materials using a structural measure quantifying the anisotropic distribution of pore voids. The approach uses a volume averaging method through which the macroscopic tortuosity tensor is related to both the average porosity and the directional distribution of pore spaces. The permeability tensor is derived from the macroscopic momentum balance equation of fluid in a porous medium and expressed as a function of the tortuosity tensor and the internal structure of the material. The analytical results generally agree with experimental data in the literature.     

Analytical models for the mechanical behavior of closed and open-cell Al foams

Two 3D analytical models are proposed for the determination of mechanical properties of Al closed and open-cell foams under compression load. The first model, referring to closed cell foams, is symmetrical, considering ellipsoid cells equally arranged in a rectangular plate, whereas the second one, related to open-cell foams, consists of a simple unit parallelepiped cell.The closed cell model produces much higher values of the plateau stress than comparable experimental results, mainly due to the associated conditions of symmetry, contrary to the open-cell model which yields values close to experimental and theoretical results of other investigators. Additionally, in the latter case, a unit cubic cell is also considered for comparison reasons. Both models are solved by the finite element method using a commercial program. The open-cell model is simple, time sparing and easy to use. Finally, a fracture analysis of the model is conducted based on the energy density concept and results are given for distortion and dilatational effects.     

Effects of microstructure on uniaxial strength asymmetry of open-cell foams

Based on the elongated Kelvin model, the effect of microstructure on the uniaxial strength asymmetry of open-cell foams is investigated. The results indicate that this asymmetry depends on the relative density, the solid material, the cell morphology, and the strut geometry of open-cell foams. Even though the solid material has the same tensile and compressive strength, the tensile and compressive strength of open-cell foams with asymmetrical sectional struts are still different. In addition, with the increasing degree of anisotropy, the uniaxial strength as well as the strength asymmetry increases in the rise direction but reduces in the transverse direction. Moreover, the plastic collapse ratio between two directions is verified to depend mainly on the cell morphology. The predicted results are compared with Gibson and Ashby''s theoretical results as well as the experimental data reported in the literature, which validates that the elongated Kelvin model is accurate in explaining the strength asymmetry presented in realistic open-cell foams.     

Modeling of the role of defects in sintered FeCrAlY foams

The metal sintering approach offers a cost- effective means for the mass-production of open-cell foams from a range of materials, including high-temperature steel alloys, which offer novel mechanical and acoustic properties. In a separate experimental study, the mechanical properties of open-celled steel alloy (FeCrAlY) foams have been characterized under uniaxial compression and shear loading. Compared to predictions from established models, a significant knockdown in material properties was observed. This knockdown was attributed to the presence of defects throughout the microstructure that result from the unique fabrication process. In the present paper, the microstructure of sintered FeCrAlY foams was modeled by using a finite element (FE) model. In particular, microstructural variations were introduced to a base lattice, and the effects on the strength and stiffness calculated. A range of defects identified under scanning electronic microscope (SEM) imaging were considered including broken ligaments, thickness variations, and pore blockages, which are the three primary imperfections observed in sintered foams. The corresponding levels of defect present in the material were subsequently input into the FE model, with the resulting predictions correlating well with experimental data. The project supported by the National Basic Research Program of China (2006CB601202), the National Natural Science Foundation of China (10328203, 10572111, 10632060), the National 111 Project of China (B06024), and the US Office of Naval Research (N000140210117).     

Geometrical modelling of foam structures using implicit functions

Mechanical, thermo-mechanical, and fluid dynamic simulations of open-cell foams require an accurate geometry model. Usually, models are derived from computer- tomography (CT) data which do not allow analysing systematically variation and optimisation of the geometry. On the other hand, entirely computer generated models are mostly assembled of primitive objects like cylinders. This disregards strut thickness variations and node rounding which are observed in real open-cell foams. This paper presents an approach to generate models of ceramic open-cell foams using simple objects with variable thickness generated by implicit functions. This approach can also account for cavities within struts and nodes, which are observed in many real foam structures. The specific rounding at the foam nodes can be modelled by applying the transformation of Blinn. The quality of the generated foam models is verified using CT data of real foams.     

Hydrodynamic Characterization of Nickel Metal Foam,Part 1: Single-Phase Permeability

This article presents results of the investigation of the fluid dynamic behavior in CVD processed nickel metal foams. An experimental facility was developed to measure the single-phase permeability in nickel metal foams in Darcian flow regime. Data on permeability values of seven different nickel foam samples was obtained. The pore sizes of the foam were obtained with scanning electron microscope. By defining friction factor and Reynolds number based on the permeability length scale a correlation was obtained for the foam permeability in Darcian flow regime. The result from this study was compared with the correlations reported by other researchers, and was found to be in good agreement.     

Pressure drop and friction factor of steady and oscillating flows in open-cell porous media

Open-cell metal foams are often used in heat exchangers and absorption equipment because they exhibit large specific surface area and present tortuous coolant flow paths. However, published research works on the characteristics of fluid flow in metal foams are relatively scarce, especially for the flow oscillation condition. The present experimental investigation attempts to uncover the behavior of steady and oscillating flows through metal foams with a tetrakaidecahedron structure. In the experiments, steady flow was supplied by an auto-balance compressor and flow oscillation was provided by an oscillating flow generator. The pressure drop and velocity were measured by the differential pressure transducer and hot-wire sensor, respectively. The friction factor of steady flow in metal foam channel was analyzed through the permeability and inertia coefficient of the porous medium. The results show that flow resistance in the metal foams increases with increasing form coefficient and decreasing permeability. The empirical equation obtained by the present study indicates that the maximum friction factor of oscillating flow through the tested aluminum foams with specific structure is governed by the hydraulic ligament diameter-based kinetic Reynolds number and the dimensionless flow amplitude.     

Micro-structural Impact of Different Strut Shapes and Porosity on Hydraulic Properties of Kelvin-Like Metal Foams

Various ideal periodic isotropic structures of foams (tetrakaidecahedron) with constant ligament cross section are studied. Different strut shapes namely circular, square, diamond, hexagon, star, and their various orientations are modeled using CAD. We performed direct numerical simulations at pore scale, solving Navier–Stokes equation in the fluid space to obtain various flow properties namely permeability and inertia coefficient for all shapes in the porosity range, \(0.60<\varepsilon <0.95\) for wide range of Reynolds numbers, \(10^{-6} . We proposed an analytical model to obtain pressure drop in metallic foams in order to correlate the resulting macroscopic pressure and velocity gradients with the Ergun-like approach. The analytical results are fully compared with the available numerical data, and an excellent agreement is observed.     

储层含水条件下致密砂岩/页岩无机质纳米孔隙气相渗透率模型

页岩及致密砂岩储层富含纳米级孔隙,且储层条件下页岩孔隙(尤其无机质孔隙)及致密砂岩孔隙普遍含水,因此含水条件下纳米孔隙气体的流动能力的评价对这两类气藏的产能分析及生产预测具有重要意义.本文首先基于纳米孔隙内液态水及汽态水热力学平衡理论,量化了储层孔隙含水饱和度分布特征;进一步在纳米孔隙单相气体传质理论的基础上,考虑了孔隙含水饱和度对气体流动的影响;最终建立了含水饱和度与气相渗透率的关系曲线. 基于本文岩心孔隙分布特征,计算结果表明:储层含水饱和度对气体流动能力的影响不容忽视,在储层含水饱和度20%的情况下,气相流动能力与干燥情况相比将降低约10%;在含水饱和度40% 的情况下,气相流动能力将降低约20%.      

Hydrodynamic Characterization of Nickel Metal Foam,Part 2: Effects of Pore Structure and Permeability

In this article, the structural characterization of chemical vapor deposition (CVD) nickel metal foam is presented. Scanning electron microscope and post image processing were used to carefully analyze the surface of the nickel metal foams. Data on foam unit cell, ligament thickness, projected pore diameter, and averaged porosity was obtained. Unit cell and projected pore diameters of CVD nickel metal foam possess Gaussian-like distribution. Characteristics of pore structure and its effect on permeability in Darcian flow regime were analyzed. The relations between the permeability, pore size, and porosity are presented. Present and previous data are compared with these relations. Measurement results indicate that the permeability or the viscous conductivity of the CVD processed metal foam is affected not only by the pore size, and porosity but also by the ligament structure.     

A constitutive approach to 3-d nonlinear fluid flow through finite deformable porous continua

The present paper proposes a thermodynamically consistent Forchheimer-type filter law for application in macroscopic porous media theories. The constitutive flow equation is thereby capable of describing the essential nonlinearities during 3-d fluid percolation through deformable porous solids. In particular, tortuosity effects, anisotropic properties, and the indispensable influence of finite distortions of the interconnected pore space are accounted for. However, the common shape of a Darcy-type relation is retained by assigning all nonlinearities to a general permeability tensor. Finally, to show the validity and applicability of the proposed formulation, the filter law is correlated with the data of permeability experiments on a high-porosity polyurethane foam and is used in a 3-d finite element analysis to simulate the pneumatic damping properties of the material.     

BIOT’S CONSOLIDATION ANALYSIS FOR CYLINDRICAL SOIL-SAMPLE BASED ON HANSBO’S FLOW1)

诸多黏性土渗流试验表明,在低水力梯度下渗流会出现明显偏离Darcy定律的现象. 为了分析渗流的非Darcy特性对固结过程 的影响,引入Hansbo渗流方程描述圆柱土样内的渗流,重新推导轴对称条件下的Biot固结方程,并给出方程的Crank--Nicolson有限 差分格式. 通过与Darcy渗流条件下轴对称Biot固结方程解析解的对比,验证计算方法的有效性. 然后分析Hansbo模型参数对圆柱 土样固结过程的影响. 计算结果表明：与Darcy渗流相比,Hansbo渗流会延缓圆柱土样的固结过程. 随着Hansbo渗流参数m或I¹的增大,在固结前期,Mandel--Cryer效应会更加显著,即孔隙水压力峰值将提高,且达到该峰值的时间 会延迟;在固结中后期,孔压消散滞后的现象也更加明显. 不过,Hansbo渗流对位移的影响很小.