流体饱和标准线性粘弹性多孔介质中的平面波

研究了流体饱和不可压标准线性粘弹性多孔介质中平面波的传播和反射问题.在固相骨架小变形的假定下,得到了粘弹性多孔介质中波动方程的一般解,讨论了弥散关系和波的衰减特性.结果表明:在流体饱和不可压粘弹性多孔介质中,仅存在一个耦合纵波和一个耦合横波,纵波和横波的波速、衰减率等取决于孔隙流体与固相骨架间的相互作用以及固相骨架本身的粘性.同时,研究了半空间自由边界上入射波(纵波、横波)的反射问题,得到了非均匀反射波的波速、反射系数、衰减率等的表达式及其相关的数值结果.     

一维流体饱和粘弹性多孔介质层的动力响应

本文研究了不可压流体饱和粘弹性多孔介质层的一维动力响应问题。基于粘弹性理论和多孔介质理论,在流相和固相微观不可压、固相骨架服从粘弹性积分型本构关系和小变形的假定下,建立了不可压流体饱和粘弹性多孔介质层一维动力响应的数学模型,利用Laplace变换,求得了原初边值问题在变换空间中的解析解,并利用Laplace逆变换的Crump数值反演方法,得到原动力响应问题的数值解。数值研究了饱和标准线性粘弹性多孔介质层的动力响应,分析了固相位移、渗流速度、孔隙压力及固相有效应力等的响应特征。结果表明,与不可压流体饱和弹性多孔介质相同,不可压流体饱和粘弹性多孔介质中亦只存在一个纵波,并且固相骨架的粘性对动力行为有显著的影响。     

不可压饱和粘弹性多孔介质固结问题的Gurtin型变分原理

基于多孔介质理论,在固相骨架和孔隙流体微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,利用卷积积分的性质,本文首先建立了以固相骨架位移、孔隙流体相对速度和孔隙流体压力为宗量的流体饱和粘弹性多孔介质固结问题的一个Gurtin型变分原理.其次,利用Lagrange乘子法解除相关的变分约束条件,建立了流体饱和粘弹性多孔介质固结问题的若干广义Gurtin型变分原理,包括第三类的Hu-Washizu型变分原理.最后,简单讨论了等价初边值问题的相应变分原理.这些Gurtin型变分原理的建立不仅丰富了饱和粘弹性多孔介质的相关理论,而且为相关数值模拟方法,如有限元法、无网格法等的建立奠定了理论基础.     

热局部非平衡条件下含柱形空洞横观各向同性饱和多孔介质的热应力分析

提出了一个热局部非平衡条件下横观各向同性多孔介质的弹性力学模型,研究了含柱形空洞的无限大多孔介质在冷对流边界条件下的孔隙压力和固相热应力.应用Laplace变换反演法数值分析了热局部非平衡及材料的各向异性对孔隙压力和固相热应力的影响.结果表明:孔隙压力和固相热应力随着横向热膨胀系数的增大而显著增大,而弹性模量的各向异性对孔隙压力和固相热应力的影响则是微弱的.同时,结果还证实了在Biot数为中等值范围内,热局部非平衡效应是非常明显的.     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

不可压饱和多孔弹性简支梁的动力响应

在杆件弯曲小变形的假定下,考虑杆件的侧向变形因素,根据多孔介质理论,本文首先建立了不可压饱和多孔弹性梁弯曲变形时动力响应的控制方程。其次,基于所建立的控制微分方程,利用变量分离法,研究了两端可渗透的饱和多孔弹性简支梁在梁中间集中载荷作用下的动力响应,得到了不同物性参数下简支梁动态弯曲时挠度和孔隙流体压力等效力偶等随时间的响应曲线。研究发现由于孔隙流体和固相骨架的相互作用,不可压饱和多孔弹性梁挠度的动力响应具有粘性特征,同时,随着时间的增加,饱和多孔弹性梁的挠度、弯矩等最终趋于经典弹性梁的静挠度、弯矩,此时,孔隙流体压力为零,梁的固相骨架承担所有的外载荷。     

饱和黏弹性多孔介质中的平面波及能量耗散

研究了流体饱和不可压黏弹性多孔介质中的非均匀平面波及其能量流和能量耗散规律. 在流 相和固相物质微观不可压、固相骨架宏观服从积分型本构关系和小变形的假定下，利用 Helmholtz分解，得到了饱和黏弹性多孔介质中非均匀平面波的一般解以及纵波、横波相速 度和衰减率等的解析表达式，分析了平面波传播矢量和衰减矢量之间的关系. 数值结果表明 孔隙流体与固相骨架间的相互作用以及固相骨架的黏性对波的相速度、衰减率等有着显著的 影响. 同时，得到了饱和黏弹性多孔介质的能量方程，给出了能量流矢量和能量耗散率. 对 非均匀平面纵波和横波，推导了平均能量流矢量和平均能量耗散率的解析表达式.     

热局部非平衡下一维半无限长饱和多孔柱体的温度场分析

对于端部受温度载荷的一维半无限长多孔介质柱体,给出了热局部非平衡下固相和流相温度场在Laplace变换域中的解析表达式.对于冲击温度载荷的情况,获得了温度场在短时间内的Laplace逆变换渐近解析解.数值分析了流、固两相热扩散系数之比以及热交换系数对固相和流相温度场的影响,比较了热局部非平衡下加权温度与热局部平衡下温度之间的差别.     

饱和多孔弹性Timoshenko悬臂梁的拟静力弯曲

在经典单相Timoshenko梁变形和孔隙流体仅沿多孔梁轴向运动的假定下,基于不可压饱和多孔介质的三维理论,论文首先建立了横观各向同性饱和多孔弹性Timoshenko悬臂梁拟静力弯曲的一维数学模型,并给出了相应的边界条件.其次,利用Laplace变换及其数值逆变换,分析了端部不同渗透条件下,饱和多孔弹性Timoshenko悬臂梁在端部梯载荷作用下的拟静力响应,给出了饱和多孔Timoshenko悬臂梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等随时间的响应曲线,并与饱和多孔Euler-Bernoulli悬臂梁的响应进行了比较,考察了梁长细比对弯曲响应的影响.数值结果表明:固相骨架与孔隙流体的相互作用具有粘性效应,梁弯曲的拟静态挠度具有蠕变行为,端部渗透条件对梁的弯曲变形有显著的影响,并且,饱和多孔弹性Timoshenko悬臂梁的拟静态响应亦存在Mandel-Cryer现象.     

饱和多孔介质化学-热-弹性模型及应用

本文基于热局部非平衡(LTNE)条件和加权平均温度概念,并假设孔隙流体由溶质和溶剂两组元组成,对页岩（饱和多孔介质）,推导给出了一种LTNE条件下的化学-热-弹性模型,同时讨论了耦合方程组的解耦求解问题．作为模型的应用,考虑无限大平面含一圆形孔的情况,研究了冷/热对流以及溶质摩尔分数突变边界条件下圆孔附近的孔隙压力和化-热应力问题,用Laplace变换得到了平面轴对称情况下有关力学变量的表达式．数值分析了圆孔边界上冷/热对流的Biot数和溶质摩尔分数改变量对圆孔附近孔隙压力和化-热应力的影响．结果表明：在Biot数为中等值(1~5)范围内,LTNE效应是非常明显的;化学作用对孔隙压力和固相应力的影响不可忽视．     

GURTIN-TYPE VARIATIONAL PRINCIPLES FOR DYNAMICS OF A NON-LOCAL THERMAL EQUILIBRIUM SATURATED POROUS MEDIUM

Based on the porous media theory and by taking into account the effects of the pore fluid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fluid phases, a mathematical model for the dynamic-thermo-hydro-mechanical coupling of a non-local thermal equilibrium fluid-saturated porous medium, in which the two constituents are assumed to be incompressible and immiscible, is established under the assumption of small deformation of the solid phase, small velocity of the fluid phase and small temperature changes of the two constituents. The mathematical model of a local thermal equilibrium fluid-saturated porous medium can be obtained directly from the above one. Several Gurtin-type variational principles,especially Hu-Washizu type variational principles, for the initial boundary value problems of dynamic and quasi-static responses are presented. It should be pointed out that these variational principles can be degenerated easily into the case of isothermal incompressible fluid-saturated elastic porous media, which have been discussed previously.     

Examination of the Thermal Equilibrium Assumption in Transient Natural Convection Flow in Porous Channel

The local thermal equilibrium assumption in the transient natural convection channel flow is investigated numerically. The Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous domain. The effect of different parameters on the validity of the local thermal equilibrium assumption is examined. It is found that the volumetric Nusselt number has the most significant effect on the local thermal equilibrium assumption.     

Validation of the Local Thermal Equilibrium Assumption in Forced Convection of Non-Newtonian Fluids through Porous Channels

In this paper, we assess the validity of the local thermal equilibrium assumption in the non-Newtonian forced convection flow through channels filled with porous media. For this purpose, the problem is solved numerically using local thermal non-equilibrium and non-Darcian models. Numerical solutions obtained over broad ranges of representative dimensionless parameters are utilized to map conditions at which the local thermal equilibrium assumption can or cannot be employed. The circumstances of a higher modified Peclet number, a lower modified Biot number, a lower fluid-to-solid thermal conductivity ratio, a lower power-law fluid index, and a lower microscopic and macroscopic frictional flow resistance coefficients, are identified as unfavorable circumstances for the local thermal equilibrium (LTE) condition to hold. Quantitative LTE validity maps that reflect the proportional effect of each parameter as related to others are presented.     

Thermo-Osmosis Effect in Saturated Porous Medium

In this paper, the thermo-poroelasticity theory is used to investigate the quasi-static response of temperatures, pore pressure, stress, displacement, and fluid flux around a cylindrical borehole subjected to impact thermal and mechanical loadings in an infinite saturated poroelastic medium. It has been reported in literatures that coupled flow known as thermo-osmosis by which flux is driven by temperature gradient, can significantly change the fluid flux in clay, argillaceous and many other porous materials whose permeability coefficients are very small. This study presents a mathematical model to investigate the coupled effect of thermo-osmosis in saturated porous medium. The energy balance equations presented here fulfill local thermal non-equilibrium condition (LTNE) which is different from the local thermal equilibrium transfer theory, accounting for that temperatures of solid and fluid phases are not the same and governed by different heat transfer equations. Analytical solutions of temperatures, pore pressure, stress, displacement, and fluid flux are obtained in Laplace transform space. Numerical results for a typical clay are used to investigate the effect of thermo-osmosis. The effects of LTNE on temperatures, pore pressure, and stress are also studied in this paper.     

A (Dual) Network Model for Heat Transfer in Porous Media

Transport in Porous Media - We present a dual network model to simulate coupled single-phase flow and energy transport in porous media including conditions under which local thermal equilibrium...     

Numerical investigation of transpiration and ablation cooling

To predict the integral performance of transpiration and ablation cooling during the reentry of hypersonic vehicles, an unsteady numerical model based on the assumption of thermal equilibrium is presented. The non-thermal equilibrium model and the thermal equilibrium model are coupled by the effective thermal properties of the porous matrix and the coolant. The calculation using the thermal equilibrium model shows the influence of the variation of the effective thermal properties on the numerical results by a comparison between constant and variable thermal properties. The comparison indicates that near the melting temperature of the porous matrix, the position of the moving boundary due to ablation is sensitive to the temperature, therefore, the variation of the thermal properties are considered in this paper. The process of ablation and transpiration cooling is simulated under different numerical conditions. The simulations demonstrate that the injection rate of coolant mass flow and initial temperature of cooling are important parameters for the control of the ablation process.     

Steady and Unsteady Heat Transfer in a Channel Partially Filled with Porous Media Under Thermal Non-Equilibrium Condition

Steady and pulsatile flow and heat transfer in a channel lined with two porous layers subject to constant wall heat flux under local thermal non-equilibrium (LTNE) condition is numerically investigated. To do this, a physical boundary condition in the interface of porous media and clear region of the channel is derived. The objective of this work is, first, to assess the effects of local solid-to-fluid heat transfer (a criterion indicating on departure from local thermal equilibrium (LTE) condition), solid-to-fluid thermal conductivity ratio and porous layer thickness on convective heat transfer in steady condition inside a channel partially filled with porous media; second, to examine the impact of pulsatile flow on heat transfer in the same channel. The effects of LTNE condition and thermal conductivity ratio in pulsatile flow are also briefly discussed. It is observed that Nusselt number inside the channel increases when the problem is tending to LTE condition. Therefore, careless consideration of LTE may lead to overestimation of heat transfer. Solid-to-fluid thermal conductivity ratio is also shown to enhance heat transfer in constant porous media thickness. It is also revealed that an increase in the amplitude of pulsation may result in enhancement of Nusselt number, while Nusselt number has a minimum in a certain frequency for each value of amplitude.     

Lattice Boltzmann Pore Scale Simulation of Natural Convection in a Differentially Heated Enclosure Filled with a Detached or Attached Bidisperse Porous Medium

In this research, pore scale simulation of natural convection in a differentially heated enclosure filled with a conducting bidisperse porous medium is investigated using the thermal lattice Boltzmann method. For the first time, the effect of connection of the bidisperse porous medium to the enclosure walls is studied by considering the attached geometry in addition to the detached one. Effect of most relevant parameters on the streamlines and isotherms as well as hot wall average Nusselt number is studied for two of the bidisperse porous medium configurations. It is observed that effect of geometrical and thermo-physical parameters of the bidisperse porous medium on the heat transfer characteristics is more complicated for the attached configuration. To assess the validity of the local thermal equilibrium condition in the micro-porous media, the pore scale results are used to compute the percentage of the local thermal non-equilibrium for two of the bidisperse porous medium configurations. It is concluded that for the detached configuration, the local thermal equilibrium condition is confirmed in the entire micro-porous media for the ranges of the parameters studied here. However, for the attached geometry, it is shown that departure from the local thermal equilibrium condition is observed for the higher values of the Rayleigh number, micro-porous porosity, solid–fluid thermal conductivity ratio, and the smaller values of the macro-pores volume fraction.     

Modeling and Simulation of Local Thermal Non-equilibrium Effects in Porous Media with Small Thermal Conductivity

The two-equation model in porous media can describe the local thermal non-equilibrium (LTNE) effects between fluid and solid at REV scale, with the temperature differences in a solid particle neglected. A multi-scale model has been proposed in this study. In the model, the temperature differences in a solid particle are considered by the coupling of the fluid energy equation at REV scale with the heat conduction equation of a solid particle at pore scale. The experiments were conducted to verify the model and numerical strategy. The multi-scale model is more suitable than the two-equation model to predict the LTNE effects in porous media with small thermal conductivity. The effects of particle diameter, mass flow rate, and solid material on the LTNE effects have been investigated numerically when cryogenic nitrogen flows through the porous bed with small thermal conductivity. The results indicate that the temperature difference between solid center and fluid has the same trend at different particle diameters and mass flow rates, while the time to reach the local thermal equilibrium is affected by solid diameter dramatically. The results also show that the temperature difference between solid center and surface is much greater than that between solid surface and fluid. The values of \( \rho {\text{c}} \) for different materials have important influence on the time to reach the local thermal equilibrium between solid and fluid.     

Plane stress state problems for notched bodies whose plastic properties depend on the form of the stress state

We consider one possible approach to the problem of describing the dependence of material plastic strain characteristics on the stress hydrostatic component arising in many porous, fractured, and other inhomogeneous materials. The plastic strain of the media under study is investigated under the plasticity assumption in the corresponding generalized form with the use of the form parameter of the stress state. The plasticity constitutive relations are stated on the basis of the plastic flow law associated with the accepted plasticity condition. For the conditions of plane stress state in the framework of the material rigid-plastic model, a system of partial differential equations is obtained and conditions for its hyperbolicity are determined. The relations for determining the stress fields and velocity fields in plastic domains are obtained, and their properties are investigated. The problem of tension of a strip with symmetric angular notches is solved, where the stress fields are determined and the continuous displacement rate field is constructed. The problem of uniform symmetric tension of a plane with a circular hole is considered. The stress fields in a strip with symmetric circular notches are examined. A comparison with solutions for plastically incompressible media whose properties are invariant with respect to the form of the stress state is performed.