铝熔体在多孔介质中的渗流过程

在流体流动相似原理的基础上，提出了铝熔体在多孔介质间隙（简称多孔介质）中渗流模拟的原理和方法．设计建造了用于铝液和模拟液在外加压力下进行单向渗流的两套试验装置模拟试验揭示渗流时铝液在多孔介质中渗流的若干规律实践表明理论模型与试验吻合良好．     

离心渗铸工艺中铝熔体在SiC多孔介质内的渗流传热过程特征及最弱环损伤模型

针对离心力场中铝熔体在SiC多孔介质内的渗流传热现象,考虑离心力对渗流传热过程的影响,根据局部非热平衡假设建立了多孔介质渗流传热模型。采用全隐格式TDMA算法和第一类迎风差分方法对渗流过程的温度场进行了数值计算。研究分析了不同复合层厚度下离心渗透过程中的流场和温度场瞬态变化规律。计算结果表明,在渗透区域,熔体与SiC颗粒存在着一定温差,而在渗透前沿,这种温差相对较大。渗流速度变化存在两个十分明显的阶段,渗流速度较高且急剧下降的初始渗透阶段以及渗流速度相当平稳的后续阶段。渗流速度的这种瞬态变化规律主要是多孔介质内流体流动与离心压力相互作用的结果。渗透初期形成的紊流状态,是导致熔体卷吸空气、使复合材料内部形成气孔的主要原因之一。选择合适的工艺参数对于确保铸件质量是十分关键的。     

铝熔体在多孔介财中的渗流过程

在流体流动相似原理的基础上，提出了铝熔体在多孔介南间隙中渗流模拟的原理和方法。设计建造了用于铝液和模拟液在外加压国下进行单向渗流的两套试验装置。     

多孔介质的一种流-固耦合动态边界理论

基于Biot理论,推导了考虑渗流作用的可变形多孔介质流-固耦合问题的基本方程,建立了本问题的渗流模型,给出了所考虑问题的流体动力弥散分布概率及系数表达式,进一步建立了多孔介质中微压液体位移场模型,讨论了流体动力弥散因素对多孔介质边界的影响,建立了描述多孔介质的动态边界方程,并分别对所建立的四种边界模拟了动态结果及算例。     

高温熟料非稳态非热平衡渗流换热模型

回转窑卸入篦冷机的高温水泥熟料为红热半透明的多孔介质,其复杂的气固换热机理给篦冷机的工艺改进带来较大难度。针对这一问题,将高温红热颗粒等效为光学厚介质,推导了一种高温熟料颗粒间传导与辐射综合换热系数,基于渗流力学与传热学理论,建立了考虑高温熟料颗粒间热辐射效应的水泥熟料非稳态非热平衡渗流换热模型。通过对所建模型进行求解,得到了料层内熟料温度与气体温度的分布规律,比较了不同区域冷却速率的差异,获得了辐射传热因素对料层温度分布的影响。     

饱和多孔介质一维瞬态波动问题的解析分析

采用基于混合物理论的多孔介质模型,提出了饱和多孔介质一维动力响应的初边值问题。利用拉氏变换和卷积定理,分别得到了边界自由排水时在任意应力边界条件和任意位移边界条件下瞬态波动过程的解析表达。几种典型的数值算例同时给出了两类边界条件下瞬态波动过程中多孔固体的位移场、应力场和孔隙流体的速度场、压力场。结果表明,饱和多孔介质的波动过程是多孔固体和孔隙流体中以同一速度传播的两种波动的耦合过程,时效特性分析也揭示了饱和多孔介质固有的表观粘弹性性质。     

饱和多孔介质中动力渗流耦合分析的Biot-Cosserat连续体模型与应变局部化有限元模拟

提出了适用于饱和多孔介质中应变局部化分析及动力渗流耦合分析的Biot-Cosserat连续体模型。基于饱和多孔介质动力渗流耦合分析的Biot理论,将固体骨架看作Cosserat连续体,并考虑旋转惯性,建立了饱和多孔介质动力渗流耦合分析的Biot-Cosserat连续体模型。基于Galerkin加权余量法,对所发展的模型推导了以固体骨架广义位移(包含旋转)及孔隙水压力为基本未知量的有限元公式。利用所发展的数值模型,对包含压力相关弹塑性固体骨架材料的饱和多孔介质进行了动力渗流耦合分析与应变局部化有限元模拟,结果表明,所发展的两相饱和多孔介质动力渗流耦合分析的Biot-Cosserat连续体模型能保持饱和两相介质应变局部化问题的适定性及模拟饱和多孔介质中由应变软化引起的应变局部化现象的有效性。     

不可压渗流驱动问题混合有限元/间断有限元耦合格式的误差分析

多孔介质中渗流驱动问题与环境污染和油藏开采等问题密切相关,是当今的研究热点.对具有分子扩散和弥散效应的不可压渗流驱动问题,本文用混合有限元/间断有限元耦合格式来求解,即用混合有限元方法求解压力方程,用对称内罚间断有限元方法逼近浓度方程.运用比剪切算子更为便捷的归纳假设和插值投影,导出了先验hp误差估计.     

多孔介质模型镁合金流变铸轧LBM模拟

根据REV尺度多孔介质格子-Boltzmann方法基本理论,建立求解高固相率半固态浆料渗流过程的热流耦合数值模型.应用模型对变孔隙率Poiseuille流动及纯扩散凝固进行模拟.LBM模拟结果与有限差分模拟吻合.分析了压差作用下浆料在平板间渗流冷却与热流耦合过程.模拟结果表明浆料渗流速度随温度降低而减慢,相同时刻下凝固潜热较大的浆料渗流Re数高.利用提出的模型初步探索了AZ91D镁合金半固态流变铸轧过程的热流场分布.     

多孔介质材料在低温下的传热特性实验研究

对多孔介质材料在低温下的传热特性进行了实验研究,在填充液氮以后其低温维持时间明显增加,主要原因是由于多孔材料的参与改变了传热特性;采用连续介质管束模型,用有限元分析软件对其整体温度场分布进行了数值模拟计算,计算结果和实验数据吻合。     

An efficient finite element procedure for analyzing three‐phase porous media based on the relaxed Picard method

Effective simulation of the solid‐liquid‐gas coupling effect in unsaturated porous media is of great significance in many diverse areas. Because of the strongly nonlinear characteristics of the fully coupled formulations for the three‐phase porous media, an effective numerical solution scheme, such as the finite element method with an efficient iterative algorithm, has to be employed. In this paper, an efficient finite element procedure based on the adaptive relaxed Picard method is developed for analyzing the coupled solid‐liquid‐gas interactions in porous media. The coupled model and the finite element analysis procedure are implemented into a computer code PorousH2M, and the proposed procedure is validated through comparing the numerical simulations with the experimental benchmarks. It is shown that the adaptive relaxed Picard method has salient advantage over the traditional one with respect to both the efficiency and the robustness, especially for the case of relatively large time step sizes. Compared with the Newton‐Raphson scheme, the Picard method successfully avoids the unphysical ‘spurious unloading’ phenomenon under the plastic deformation condition, although the latter shows a better convergence rate. The proposed procedure provides an important reference for analyzing the fully coupled problems related to the multi‐phase, multi‐field coupling in porous media. Copyright © 2014 John Wiley & Sons, Ltd.     

恒压下树脂在双尺度多孔介质中的流动特性

基于复合材料液态模塑(LCM)工艺过程中存在半饱和区域的实验现象以及对预制体双尺度效应的逐步认识, 一些学者提出用沉浸模型来研究双尺度多孔介质的不饱和流动。通过体积均匀化方法描述了双尺度多孔介质复合材料液态模塑工艺模型的特征, 得到含有沉浸项的双尺度多孔介质的质量守恒方程, 并采用有限元法对方程进行数值求解, 通过具体算例计算了考虑双尺度效应时恒压树脂注射下不同时段的压力分布状态, 得到树脂在填充过程中流动前沿半饱和区域从出现到消失的过程, 采用不同注射压力进行模拟并比较。结果表明, 与单尺度多孔介质模型不同, 双尺度多孔介质模型更能反映实际树脂填充过程中出现的半饱和区域现象。     

A consistent u‐p formulation for porous media with hysteresis

This paper presents a continuum formulation based on the theory of porous media for the mechanics of liquid unsaturated porous media. The hysteresis of the liquid retention model is carefully modelled, including the derivation of the corresponding consistent tangent moduli. The quadratic convergence of Newton's method for solving the highly nonlinear system with an implicit finite element code is demonstrated. A u‐p formulation is proposed where the time discretisation is carried out prior to the space discretisation. In this way, the derivation of all consistent moduli is fairly straightforward. Time integration is approximated with the Theta and Newmark's methods, and hence the fully coupled nonlinear dynamics of porous media is considered. It is shown that the liquid retention model requires also the consistent second‐order derivative for quadratic convergence. Some predictive simulations are presented illustrating the capabilities of the formulation, in particular to the modelling of complex porous media behaviour. Copyright © 2014 John Wiley & Sons, Ltd.     

率相关饱和多孔介质动力响应的数值分析

在基于混合物理论的多孔介质模型的基础上,将固体相视为弹粘塑性体,建立了饱和多孔介质的弹粘塑性模型。模型的基本思想是在无粘弹塑性本构关系中引入-时间参数,使固体骨架具备了粘性效应。利用Galerkin加权残值法推导得到了罚有限元格式,并采用Newmark预估校正法求解率相关饱和多孔介质的非线性有限元动力方程,此算法可以很...     

DYNAMIC ANALYSIS OF A FULLY SATURATED POROUS MEDIUM ACCOUNTING FOR GEOMETRICAL AND MATERIAL NON-LINEARITIES

Based on the theory of porous media (mixture theories extended by the concept of volume fractions), a model describing the dynamical behaviour of a saturated binary porous medium is presented including both geometrical and material non-linearities. Transformed toward a weak formulation, the model equations are solved by use of the finite element method. Applications of the model range from one-dimensional linear problems to two-dimensional problems including the full dynamics and non-linearities.     

A stabilized volume‐averaging finite element method for flow in porous media and binary alloy solidification processes

A stabilized equal‐order velocity–pressure finite element algorithm is presented for the analysis of flow in porous media and in the solidification of binary alloys. The adopted governing macroscopic conservation equations of momentum, energy and species transport are derived from their microscopic counterparts using the volume‐averaging method. The analysis is performed in a single domain with a fixed numerical grid. The fluid flow scheme developed includes SUPG (streamline‐upwind/Petrov–Galerkin), PSPG (pressure stabilizing/Petrov–Galerkin) and DSPG (Darcy stabilizing/Petrov–Galerkin) stabilization terms in a variable porosity medium. For the energy and species equations a classical SUPG‐based finite element method is employed. The developed algorithms were tested extensively with bilinear elements and were shown to perform stably and with nearly quadratic convergence in high Rayleigh number flows in varying porosity media. Examples are shown in natural and double diffusive convection in porous media and in the directional solidification of a binary‐alloy. Copyright © 2004 John Wiley & Sons, Ltd.     

Fully coupled hydromechanical multiscale model with microdynamic effects

In this paper, a multiscale finite element framework is developed based on the first‐order homogenization method for fully coupled saturated porous media using an extension of the Hill‐Mandel theory in the presence of microdynamic effects. The multiscale method is employed for the consolidation problem of a 2‐dimensional saturated soil medium generated from the periodic arrangement of circular particles embedded in a square matrix, which is compared with the direct numerical simulation method. The effects of various issues, including the boundary conditions, size effects, particle arrangements, and the integral domain constraints for the microscale boundary value problem, are numerically investigated to illustrate the performance of a representative volume element in the proposed computational homogenization method of fully coupled saturated porous media. This study is aimed to clarify the effect of scale separation and size dependence, and to introduce characteristics of a proper representative volume element in multiscale modeling of saturated porous media.     

A computational model for fracturing porous media

This paper presents a new computational model for simulating a fracturing process in a porous medium using the finite element method. Two independent numerical techniques are used for describing this process. The partition of unity method is used for describing the fracturing process, and the double porosity model is used for describing the resulting fluid flow. A key feature of the model is the coupling of these two independent numerical techniques, which provide the means for a better simulation of the involved physical and mechanical processes. The paper focuses on the numerical formulation of the model. The capability of the model is illustrated by means of numerical examples, which examine the behaviour of a 1D porous medium under different boundary conditions. The numerical results show that the very complicated physical and mechanical processes of the fracturing porous media can be simulated properly and efficiently. Copyright © 2006 John Wiley & Sons, Ltd.