基于蒙特卡罗随机有限元方法的随机多孔介质内流体自然对流不确定性研究

为分析孔隙率不确定性对多孔介质方腔内自然对流换热的影响,发展了一种基于KL(Karhunen-Loeve展开)-蒙特卡罗随机有限元算法的随机多孔介质内自然对流不确定性分析数理模型及有限元数值模拟程序框架。通过K-L展开及基于拉丁抽样法生成多孔介质孔隙率随机实现,并耦合多孔介质自然对流有限元程序,进行随机多孔介质内自然对流传热数值模拟,得出了多孔介质内流场与温度场平均值与标准偏差,并分析了孔隙率不确定性条件下Da数对Nu数的影响。结果表明,孔隙率不确定性对多孔介质方腔内自然对流有重要影响。随机多孔介质内流场及温度场与确定性条件下的流场及温度场存在一定偏差,Nu数标准偏差随着Da的增大先增大后减小。     

防风网PM-IBM模型及其在堆场扬尘庇护区湍流模拟中的应用

防风网是煤炭和矿石等物料堆场的重要抑尘设施，能够利用多孔网板形成风速庇护效应，减少扬尘带来的大气污染。基于计算流体力学的数值风洞是预测防风网减风抑尘效率的先进手段，然而网板边界模型及其复杂阻风作用机制一直是风场湍流模拟的难点。本文将浸入边界法与多孔介质模型相结合，通过在NS方程中引入与渗流压降对应的力源项，建立PM-IBM（Porous Medium-Immersed Boundary Method）模型，实现了结构化网格下的防风网数值边界，并应用于堆场扬尘庇护区湍流风场模拟。结果表明，防风网网高大于1.4倍堆场高时，堆顶扬尘风速随网高呈非线性增长；PM-IBM模型结合通用计算流体力学求解器FLUENT，能够快速、高精度地模拟防风网的阻风作用，为堆场扬尘抑制效率评估提供数据支持。     

基于蒙特卡罗随机有限元方法的随机多孔介质内流体自然对流不确定性研究

为分析孔隙率不确定性对多孔介质方腔内自然对流换热的影响,发展了一种基于KL（Karhunen-Loeve展开）-蒙特卡罗随机有限元算法的随机多孔介质内自然对流不确定性分析数理模型及有限元数值模拟程序框架。通过K-L展开及基于拉丁抽样法生成多孔介质孔隙率随机实现,并耦合多孔介质自然对流有限元程序,进行随机多孔介质内自然对流传热数值模拟,得出了多孔介质内流场与温度场平均值与标准偏差,并分析了孔隙率不确定性条件下Da数对Nu数的影响。结果表明,孔隙率不确定性对多孔介质方腔内自然对流有重要影响。随机多孔介质内流场及温度场与确定性条件下的流场及温度场存在一定偏差,Nu数标准偏差随着Da的增大先增大后减小。     

多孔介质壁面剪切湍流速度时空关联的研究

作为一个基础统计量,时空关联函数在湍流问题的研究中有着广泛的应用,是研究湍流噪声、湍流中物质扩散和大涡模拟亚格子模型等问题的重要参考.本文通过建立三维多孔结构壁面剪切湍流模型,采用含Darcy-Brinkman-Forchheimer作用力项的格子Boltzmann方程对无穷大多孔介质平行板之间壁湍流进行了数值模拟,进而研究其速度脉动时空关联函数的统计特性.一方面,根据计算得到的流场数据,对比分析了常规槽道湍流与多孔介质壁面槽道湍流的时间关联函数.另一方面,计算并讨论了不同孔隙率和渗透率的多孔介质壁面对速度脉动时空关联性的影响.通过研究表明:多孔结构壁面剪切湍流的时空关联函数等值线与椭圆理论相符;在研究参数范围内,多孔介质壁面的速度时空关联系数随着孔隙率增大而增大,随着渗透率增大而减小.同时发现在槽道壁面的近壁区、过渡区、对数律区和中心区等不同位置处,速度时空关联呈现较大差异性:越远离壁面位置(对数律区和中心区),其时空关联函数所呈现的关联等值线椭圆越细长,高值相关等值线越集中.多孔介质主要改变速度时空关联椭圆图像的椭圆率,说明多孔介质壁面主要影响湍流横扫速度.     

基于数值模拟方法的自然通风器风载体型系数研究

基于Fluent软件平台,采用混合网格划分,选用Realizable k-ε湍流模型,对某自然通风器的风场进行三维数值模拟研究.选择0°、45°、90°、135°和180°等风向,每个风向考虑39m/s风速,此外对0°、180°两风向还考虑了57m/s风速.通过模拟,得到了主通风器导流片的风载体型系数.经计算,0°,180°两个风向的风场对导流片抗风影响较大,基于此对这两个风向的计算结果进行了详细分析.在数值模拟过程中引人多孔介质阶跃模型模拟次通风器导流片,在降低模型复杂程度和减少计算量的同时,保证了数值模拟结果的可靠性.分析表明,基于数值模拟方法得到的风载体型系数在不同风向和风速的情况下呈现一定的变化规律性,为确定自然通风器类结构的风载体型系数提供了合理依据.     

基于多孔介质模型的养殖网箱周围流场特性研究

为探讨鱼类养殖网箱对周围海域流场的影响,基于雷诺时均方程和Realizable k-ε湍流模型,耦合多孔介质模型,建立了研究鱼类养殖网箱周围流场的三维水动力模型.根据典型养殖场的实际数据,分别对仅网衣作用、仅鱼类作用及网衣和鱼类共同作用下网箱周围流场进行数值模拟.结果表明综合考虑网衣和鱼类的作用时,网箱下游流速衰减最大...     

风吹雪廓线的风洞实验研究

利用颗粒图像测速仪对新雪形成的风吹雪进行风洞实验研究, 给出了不同高度处雪粒粒径分布函数以及平均粒径廓线、雪粒数通量廓线的分布规律. 发现当摩阻风速大于0.5 m/s 时单宽输雪率与摩阻风速满足指数函数的关系, 小于0.5 m/s 时两者满足幂函数的关系, 总体而言, 单宽输雪率与摩阻风速呈线性关系.     

无限冰介质中爆炸裂纹扩展机理与数值模拟

根据冰体力学性能,提出了无限淡水冰介质在爆炸冲击载荷下的压碎区及裂隙区的理论计算方法.利用大型有限元软件ANSYS/LS-DYNA模拟了爆炸载荷下无限区域冰体内裂纹扩展过程,模拟结果与实验结果基本一致,表明采用的损伤断裂模型能够较准确的反映无限冰介质动态力学性能及裂纹成形过程.数值模型可近似计算冰体中裂隙区的半径.数值计算值与理论计算值误差不超过5%,数值模拟结果与实验结果吻合,从而确定了无限区域冰体在-15℃时的基本破坏特征及范围.     

基于孔隙弹性耦合的水力压裂数值模拟研究

系统地阐述基于多孔介质渗流-损伤耦合原理,进行水力压裂的FEM数值实现方法.基本架构为:(1)引入孔隙流体压力膨胀系数将孔隙流体压力与应力场进行耦合;(2)基于损伤局部化模型,提出裂缝张开度表达式;(3)提出水力压裂引起的多孔介质水-力学属性的各向异性表达式;(4)提出全流量加载的耦合分析方案.作为实例,模拟三维地层水压裂缝扩展形态,通过比较模型的数值解和经典理论解,验证该方法的正确性.     

生物脏器多孔介质的物理特征

生物渗流是生物力学的组成部分,也是渗流力学的组成部分。本文简述作者在肾、肺、肝三种脏器的血管体系多孔介质物理特征——孔隙率研究的部分结果。1.实验研究方法本文所说的孔隙率是指某种生物脏器内的某种多孔介质系统的孔道空间总体积与该脏器的外表总体积之比值,以小数或百分数表之。为测定生物脏器多孔介质的孔隙率,只须测定该多孔介质的孔道空间总体积和该脏器的外表总体积。     

Tortuosity based on Anisotropic Diffusion Process in Structured Plate-like Obstacles by Monte Carlo Simulation

The diffusion of fluids in porous media, composed of regularly aligned plate-like obstacles, was studied by Monte Carlo simulation. The diffusion coefficients and all diagonal components of the diffusion tensor were estimated for these media. The calculated tortuosities were modeled as a function of porosity by using the Koponen’s equation related to percolation threshold. These results indicated that a media with a homogeneous porosity has a heterogeneous tortuosity, is affected by the alignments of the plate-like obstacles. Furthermore, the calculation results were compared with the experimental results for fixed perpendicular plates of Comiti and Renaud as a function of porosity. The results for tortuosity compared well for porosity larger than 0.86.     

Multiphase transfer model for intermittent microwave-convective drying of food: Considering shrinkage and pore evolution

Intermittent microwave convective (IMCD) drying is an advanced drying technology that improves both energy efficiency and food quality during drying. Although many experimental studies on IMCD have been conducted, there is no complete multiphase porous media model describing the physics of the process available in the literature. A multiphase porous media model considering liquid water, gases and the solid matrix of food during drying can provide in-depth understanding of IMCD process. Currently there is no IMCD model that have taken shrinkage and pore evolution during drying into consideration. In this study, first a multiphase porous media model with shrinkage (IMCD2) has been developed for IMCD. Then the model has been compared with IMCD model without shrinkage (IMCD1). Simulated temperature, moisture content, density, porosity from IMCD2 are then validated against experimental data. The profile of vapour pressures and evaporation during IMCD are also presented and discussed.     

Calculating the Anisotropic Permeability of Porous Media Using the Lattice Boltzmann Method and X-ray Computed Tomography

A lattice Boltzmann (LB) method is developed in this article in a combination with X-ray computed tomography to simulate fluid flow at pore scale in order to calculate the anisotropic permeability of porous media. The binary 3D structures of porous materials were acquired by X-ray computed tomography at a resolution of a few microns, and the reconstructed 3D porous structures were then combined with the LB model to calculate their permeability tensor based on the simulated velocity field at pore scale. The flow is driven by pressure gradients imposed in different directions. Two porous media, one gas diffusion porous layer used in fuel cells industry and glass beads, were simulated. For both media, we investigated the relationship between their anisotropic permeability and porosity. The results indicate that the LB model is efficient to simulate pore-scale flow in porous media, and capable of giving a good estimate of the anisotropic permeability for both media. The calculated permeability is in good agreement with the measured date; the relationship between the permeability and porosity for the two media is well described by the Kozeny–Carman equation. For the gas diffusion layer, the simulated results showed that its permeability in one direction could be one order of magnitude higher than those in other two directions. The simulation was based on the single-relaxation time LB model, and we showed that by properly choosing the relaxation time, it could give similar results to those obtained using the multiple-relaxation time (MRT) LB method, but with only one third of the computational costs of MRTLB model.     

海洋地震工程流固耦合问题统一计算框架

海底地震动的模拟以及海洋工程结构的地震反应分析中,涉及到海水、饱和海床、弹性基岩、结构之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩和结构,分别考虑相互之间的耦合,十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后,以饱和多孔介质方程为基础,经集中质量有限元离散,考虑不同孔隙率的饱和多孔介质之间耦合的一般情形,建立了该耦合情形的求解方法.进一步论证了该一般情形的耦合计算方法可分别退化到流体与固体、流体与饱和多孔介质、固体与饱和多孔介质之间的耦合计算,从而将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.以P-SV波垂直入射时,半无限层状海水-饱和海床、海水-弹性基岩、海水-饱和海床-弹性基岩三种情形的动力分析为例,采用统一计算框架结合透射边界条件进行求解,并与传递矩阵方法得到的解进行对比,验证了该统一计算框架的有效性以及并行计算的可行性.      

Double porous screen element for gas–liquid phase separation

We consider an assembly of two parallel porous screens suspended in a tube at a distance L. The screens are connected by wicking aids. If one screen is brought into contact with a wetting liquid, the other screen will be wetted as well enclosing gas in between. Due to surface tension in the screen pores, the gas can only be removed from the chamber when the pressure difference across one screen exceeds the bubble point. With such a double porous screen element it is therefore possible to block liquid flow using trapped gas as plug. We present a model approach, experiments and numerical calculations on the performance of such a screen element. The model is based on capillary transport in vertical and radial capillaries and allows to predict how fast the element will trap the gas to become operational. For the experiments, we have built such an element using Dutch Twilled weaves made of stainless steel. Placed in a vertical tube and initially dry, it is wetted from below or above and submitted to an increasing pressure difference until breakthrough occurs where the element fails. Corresponding numerical calculations elucidate what happens within the element when it fails. Our results confirm the concept of the double porous screen element and encourage its application as liquid management device.     

海洋地震工程流固耦合问题统一计算框架——————不规则界面情形

海底地震动场及海洋声场的模拟中,需要考虑复杂海床介质及海底地形的影响,涉及到海水、饱和海床、弹性基岩之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩,分别进行空间离散和界面耦合, 十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后, 以饱和多孔介质方程为基础, 经集中质量有限元离散,严格考虑不同孔隙率的饱和多孔介质在不规则界面的耦合条件,通过求解法向和切向界面力的途径,建立了不同孔隙率的饱和多孔介质耦合情形的求解方法,将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.考虑海水--弹性基岩、海水--饱和海床--弹性基岩体系中凹陷地形情形,采用本文提出的统一计算框架, 结合透射边界条件,分析了P波入射时的动力反应, 并通过结果是否满足界面条件,验证了该统一计算框架的有效性以及并行计算的可行性.      

Air shock wave interaction with an obstacle covered by porous material

The interaction between an air shock wave and a rigid wall covered by a porous screen is investigated numerically and experimentally. A mathematical two velocity with two stress tensors model is used for studying the wave processes in saturated porous media. The process of reflection of a step-type wave from a rigid wall covered with a porous layer is considered, the effect of the porous medium and wave parameters on the reflection is analyzed, and the numerical results are compared with the experimental data.Received: 30 July 2002, Accepted: 24 December 2002, Published online: 27 May 2003     

Moisture Removal from a Two-Layer Porous Media: A Conceptual Model and Experimental Results

Moisture removal from a two-layer porous media in which air is circulated through one layer and moisture is removed from the second has not been well studied due to the emphasis given to single-layer systems. This two-layer configuration is common in natural and engineered systems and can be used as a means to create a barrier to downward migrating fluids and to remove liquids and gases that may be present in the finer layer. However, there is little data on moisture removal from a two-layer porous media in which air is circulated through one layer parallel to the interface and moisture is removed from the finer second layer by evaporation. A conceptual model of the moisture removal from a two-layer porous media system was developed and compared to experimental moisture removal rates from laboratory scale dry barriers. The limited experimental data agrees well with the results predicted by the conceptual model, providing an initial validation.     

Forced convective heat transfer in porous medium of wire screen meshes

The hydrodynamic and heat transfer characteristics of a porous medium consisting of 20 wire screen meshes are examined theoretically and experimentally. The hydrodynamic experiments are conducted for the range of Reynolds number based on mean velocity and wire diameter from 1.5 to 12. The Ergun's constants and thermal dispersion coefficients are calculated in this range. Nusselt number variation is determined in both thermally developing and fully developed flows by the help of forced convection heat transfer experiments conducted for the uniform heat flux boundary condition. Correlation functions of Nusselt number in the range of fully developed and thermally developing, and of thermal entrance length are obtained from experimental data. Solutions of momentum and energy equations simulating the experimental model are obtained numerically with variable porosity and the anticipated thermal dispersion coefficients. The thermal dispersion coefficients well-adjusted to the experimental data are determined by numerical solution of the energy equation. Received on 22 November 1996     

A Water Retention Curve Model for the Simulation of Coupled Thermo-Hydro-Mechanical Processes in Geological Porous Media

This paper presents a new water retention curve (WRC) model for the simulation of coupled thermo-hydro-mechanical processes in geological porous media. The model simultaneously considers the impact of porosity and temperature on suction, for both wetting processes and drying processes. The model is based on an idealization of porous geological media as having an isotropic and homogeneous microscopic pore structure. Suction is expressed as a function of degree of saturation, porosity, surface tension of the water–air interface, and the length of air bubble perimeter of the pores per unit area on a random 2D cross-section of the medium. The tension of water–air interface is written as a function of temperature, and the length of perimeter of the water–air interface of the pores becomes a function of porosity and degree of saturation. The final equation of the new WRC is a function of suction, effective degree of saturation, temperature, porosity, pore-gas pressure, and the rate of degree of saturation change with time for both wetting and drying processes. The model was used to fit experimental data of the FEBEX bentonite, with good agreements between measured and calculated results.