Analysis of coupled thermo-hydro-mechanical behavior of unsaturated soils based on theory of mixtures I

This paper analyzes a coupled thermo-hydro-mechanical behavior of unsaturated soils based on the theory of mixtures. Unsaturated soil is considered as a mixture composed of soil skeleton, liquid water, vapor, dry air, and dissolved air. In addition to the mass and momentum conservation equations of each component and the energy conservation equation of the mixture, the system is closed using other 37 constitutive （or restriction） equations. As the change in water chemical potential is identical to the change in vapor chemical potential, a thermodynamic restriction relationship for the phase transition between pore water and pore vapor is formulated, in which the impact of the change in gas pressure on the phase transition is taken into account. Six final govern- ing equations are given in incremental form in terms of six primary variables, i.e., three displacement components of soil skeleton, water pressure, gas pressure, and temperature. The processes involved in the coupled model include thermal expansions of soil skeleton and soil particle, Soret effect, phase transition between water and vapor, air dissolution in pore water, and deformation of soil skeleton.     

AN ALGEBRAIC MULTIGRID METHOD FOR COUPLED THERMO-HYDRO-MECHANICAL PROBLEMS

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ａｇｒｅａｔｄｅａｌｏｆａｔｔｅｎｔｉｏｎｈａｓｂｅｅｎｄｉｒｅｃｔｅｄｔｏｗａｒｄｓｆｕｌｌｃｏｕｐｌｅｄｔｈｅｒｍｏ＿ｈｙｄｒｏ＿ｍｅｃｈａｎｉｃａｌａｎａｌｙｓｉｓｉｎｄｅｆｏｒｍｉｎｇｐｏｒｏｕｓｍｅｄｉａｄｕｅｔｏｅｘｔｒａｃｔｉｏｎｏｆｕｎｄｅｒｇｒｏｕｎｄｆｌｕｉｄｓ(ｗａｔｅｒ,ｐｅｔｒｏｌｅｕｍ ,ｎａｔｕｒａｌｇａｓ)ｉｎｒｅｓｅｒｖｏｉｒ,ｅｎｖｉｒｏｎｍｅｎｔａｎｄｃｏｎｓｔｒｕｃｔｉｏｎｅｎｇｉｎｅｅｒｉｎｇ (Ｌｅｗｉｓａｎ…     

固体炸药爆轰与惰性介质相互作用的一种扩散界面模型

提出一种保持热力学一致性的扩散界面模型，用来数值模拟固体炸药爆轰与惰性介质的相互作用问题。基于混合网格内各组分物质间可以达到力学平衡状态而不能达到热学平衡状态的假设，由混合网格能量守恒以及压力相等条件，推导出每种组分物质的体积分数演化方程。由此获得的扩散界面模型包括组分物质的质量守恒方程、混合物质的动量及总能量守恒方程，同时包括组分物质的体积分数演化方程和混合物质的压力演化方程。该扩散界面模型的主要特点是考虑了化学反应以及热学非平衡的影响。提出的扩散界面模型在物质界面附近不会出现物理量的非物理振荡现象、适用于任意表达形式的物质状态方程以及任意数目的惰性介质。     

Boiling in Porous Media: Model and Simulations

We present a modelization of the heat and mass transfers within a porous medium, which takes into account phase transitions. Classical equations are derived for the mass conservation equation, whereas the equation of energy relies on an entropy balance adapted to the case of a rigid porous medium. The approximation of the solution is obtained using a finite volume scheme coupled with the management of phase transitions. This model is shown to apply in the case of an experiment of heat generation in a porous medium. The vapor phase appearance is well reproduced by the simulations, and the size of the two-phase region is correctly predicted. A result of this study is the evidence of the discrepancy between the air – water capillary and relative permeability curves and water – water vapor ones.     

Diffuse interface model for high speed cavitating underwater systems

High speed underwater systems involve many modelling and simulation difficulties related to shocks, expansion waves and evaporation fronts. Modern propulsion systems like underwater missiles also involve extra difficulties related to non-condensable high speed gas flows. Such flows involve many continuous and discontinuous waves or fronts and the difficulty is to model and compute correctly jump conditions across them, particularly in unsteady regime and in multi-dimensions. To this end a new theory has been built that considers the various transformation fronts as ‘diffuse interfaces’. Inside these diffuse interfaces relaxation effects are solved in order to reproduce the correct jump conditions. For example, an interface separating a compressible non-condensable gas and compressible water is solved as a multiphase mixture where stiff mechanical relaxation effects are solved in order to match the jump conditions of equal pressure and equal normal velocities. When an interface separates a metastable liquid and its vapor, the situation becomes more complex as jump conditions involve pressure, velocity, temperature and entropy jumps. However, the same type of multiphase mixture can be considered in the diffuse interface and stiff velocity, pressure, temperature and Gibbs free energy relaxation are used to reproduce the dynamics of such fronts and corresponding jump conditions. A general model, based on multiphase flow theory is thus built. It involves mixture energy and mixture momentum equations together with mass and volume fraction equations for each phase or constituent. For example, in high velocity flows around underwater missiles, three phases (or constituents) have to be considered: liquid, vapor and propulsion gas products. It results in a flow model with 8 partial differential equations. The model is strictly hyperbolic and involves waves speeds that vary under the degree of metastability. When none of the phase is metastable, the non-monotonic sound speed is recovered. When phase transition occurs, the sound speed decreases and phase transition fronts become expansion waves of the equilibrium system. The model is built on the basis of asymptotic analysis of a hyperbolic total non-equilibrium multiphase flow model, in the limit of stiff mechanical relaxation. Closure relations regarding heat and mass transfer are built under the examination of entropy production. The mixture equation of state (EOS) is based on energy conservation and mechanical equilibrium of the mixture. Pure phases EOS are used in the mixture EOS instead of cubic one in order to prevent loss of hyperbolicity in the spinodal zone of the phase diagram. The corresponding model is able to deal with metastable states without using Van der Waals representation.     

混凝土中化学-热-湿-力耦合过程的数值方法

提出了一个火灾下混凝土中化学-热-湿-力耦合过程分析的两级数学模型. 混凝 土模型化为充满两种非混溶孔隙流体的非饱和变形多孔多相介质. 数学模型基于控制干空 气、湿份及基质溶解物的质量守恒、混凝土介质混合体的动量守恒和焓(能量)守恒的耦合 偏微分方程组. 模型中特别考虑到了高温下的脱盐过程. 构造了一个用于数值模拟 化学-热-湿-力耦合行为的有限元求解过程的混合弱形式. 并且针对其中具有非自伴随算子特性的 双曲线控制方程的空间离散进行了特殊考虑. 数值结果例题显示所发展的数学模型和数值方 法在重现火灾条件下的混凝土中化学-热-湿-力耦合行为的有效性.     

Nucleation and flashing in nozzles—2. Comparison with experiments using a five-equation model for vapor void development

A quasi-one-dimensional, five-equation, homogeneous, nonequilibrium model has been developed and utilized on a microcomputer to calculate the behavior of flowing, initially subcooled, flashing water systems. Equations for mixture and vapor mass conservation, mixture momentum conservation, liquid energy conservation and bubble transport were discretized and linearized semi-implicitly, and solved using a successive iteration Newton method. Closure was obtained through simple constitutive equations for friction and spherical bubble growth, and a new nucleation model for wall nucleation in small nozzles combined with an existing model for bulk nucleation in large geometries to obtain the thermal nonequilibrium between phases. The model described was applied to choked nozzle flow with subcooled water inlets based on specified inlet conditions of pressure and temperature, and vanishing inlet void fraction and bubble number density. Good qualitative and quantitative agreement with the experiment confirms the adequacy of the nucleation models in determining both the initial size and number density of nuclei, and indicates that mechanical nonequilibrium between phases is not an important factor in these flows. It is shown that bulk nucleation becomes important as the volume-to-surface ratio of the geometry is increased.     

The prediction of two-phase mixture level and hydrodynamically-controlled dryout under low flow conditions

A theory is given for the thermal-hydraulic phenomena during uncovery of a flow channel. This is relevant to a reactor core under typical small break or operational transient conditions. A distinct equivalent collapsed liquid level and a two-phase-mixture level are defined in the model. The former represents the liquid inventory in the channel, while the latter characterizes the heat transfer regimes. The definition of these levels are coupled through the mass and energy conservation equations, and the constitutive relations for void fraction and net vapor generation location.Analytical solutions are obtained for the transient variation of both the collapsed liquid and the two-phase mixture levels.The analyses have been compared with existing single-tube data with uniform heat flux, and rod bundle experiments with an axial power profile and inlet feedwater flow. The results demonstrate the potential of the present model for application to reactor conditions.     

一种全耦合多相流分析的并行计算方法1)

研究了孔隙介质中热、水和汽流全耦合分析的并行计算方法.模型中采用了考虑毛细压力关系的修正有效应力概念,并考虑了相变和潜热传递.基本变量为位移、毛细压力、汽压和温度.并行程序是在国家高性能计算中心(北京)的曙光1000A上借助PVM(Parallel Virtual Machine)软件系统实现的,考题显示出较高的并行加速比和效率.     

一种全耦合多相流分析的并行计算方法

研究了孔隙介质中热、水和汽流全耦合分析的并行计算方法．模型中采用了考虑毛细压力关系的修正有效应力概念,并考虑了相变和潜热传递．基本变量为位移、毛细压力、汽压和温度．并行程序是在国家高性能计算中心（北京）的曙光１０００Ａ上借助ＰＶＭ（ＰａｒａｌｅｌＶｉｒｔｕａｌＭａｃｈｉｎｅ）软件系统实现的,考题显示出较高的并行加速比和效率     

Multicomponent, Multiphase Thermodynamics of Swelling Porous Media with Electroquasistatics: I. Macroscale Field Equations

A systematic development of the macroscopic field equations (conservation of mass, linear and angular momentum, energy, and Maxwell's equations) for a multiphase, multicomponent medium is presented. It is assumed that speeds involved are much slower than the speed of light and that the magnitude of the electric field significantly dominates over the magnetic field so that the electroquasistatic form of Maxwell's equations applies. A mixture formulation for each phase is averaged to obtain the macroscopic formulation. Species electric fields are considered, however it is assumed that it is the total electric field which contributes to the electrically induced forces and energy. The relationships between species and bulk phase variables and the macroscopic and microscopic variables are given explicitly. The resulting field equations are of relevance to many practical applications including, but not limited to, swelling clays (smectites), biopolymers, biological membranes, pulsed electrophoresis, and chromatography.     

SOLUTIONS FOR CYLINDRICAL CAVITY IN SATURATED THERMOPOROELASTIC MEDIUM

Based on the thermodynamics of irreversible processes, the mass conservation equation and heat energy balance equation are established. The governing equations of thermal consolidation for homogeneous isotropic materials are presented, accounting for the coupling effects of the temperature, stress and displacement fields. The case of a saturated medium with a long cylindrical cavity subjected to a variable thermal loading and a variable hydrostatic pressure （or a variable radial water flux） with time is considered. The analytical solutions are derived in the Laplace transform space. Then, the time domain solutions are obtained by a numerical inversion scheme. The results of a typical example indicate that thermodynamically coupled effects have considerable influences on thermal responses.     

中心舱体与薄壁梁刚-柔耦合系统的热弹性-结构动力学分析

空间柔性结构受太阳热流冲击而诱发的振动是导致航天器失效的典型模式之一,准确预测结构热致振动的响应及稳定性是卫星设计的基础。针对常见的中心舱体与附属薄壁杆件组成的空间结构,提出了考虑刚-柔耦合、耦合热弹性和耦合热-结构三重耦合效应的热致振动分析理论模型。其中,刚-柔耦合是指舱体姿态角、顶端集中质量转动与柔性附件运动的耦合;耦合热弹性是指应变率与温度场的耦合;耦合热-结构是指舱体转动及结构变形与薄壁杆件吸收太阳热流的耦合。基于热弹性理论和Lagrange方程,推导了传热和运动的耦合方程;采用Laplace变换方法并使用Routh-Hurwitz稳定判据推导了稳定性边界方程。结果表明,该模型能够更为准确的给出热致振动响应及稳定性预测。     

Mixture modeling of jointed fluidsaturated porous media

The quasi-static equations of motion are studied for bi-laminated fluid-saturated porous media within the framework of non-phenomenological mixture theories. The flow-deformation coupled behavior of the media is governed by Biot's theory for which all constituents are considered compressible. The asymptotic analysis for a periodic microstructure with multiple scales, developed by Hegemier and Murakami, is adopted to obtain the equations of equilibrium and mass conservation in a binary saturated porous medium. The multiscale analysis appears to be advantageous for dealing with consolidation phenomena because it is capable of transforming a coupled, transient problem into two decoupled, steady-state ones. Various models with different degrees of approximation are generated, and among them a theory for saturated rocks with a single joint system is described. Mixture properties are expressed explicitly in terms of characteristics of intact and joint material. The most distinctive feature of this model comes from the fact that some cross terms, that have not been included in previous models, appear in the constitutive equations for fluid mass change and fluid flux. These cross terms are physically understood because they simply take into account effects occurring on the local level: the deformation-flow coupled phenomenon, the stress continuity and displacement compatibility conditions. These novel results may have far-reaching consequences for future theoretical modeling and experimental programs in two-phase fluid-filled porous media.     

Two-fluid model for transient analysis of slug flow in oil wells

In this work it is presented a transient, one-dimensional, adiabatic model for slug flow simulation, which appears when liquid (mixture of oil and water) and gas flow simultaneously through pipes. The model is formed by space and time averaged conservation equations for mass, momentum and energy for each phase, the numerical solution is based on the finite difference technique in the implicit scheme. Velocity, pressure, volumetric fraction and temperature profiles for both phases were predicted for inclination angles from the horizontal to the vertical position (unified model) and ascendant flow. Predictions from the model were validated using field data and ten correlations commonly used in the oil industry. The effects of gas heating or cooling, due to compression and expansion processes, on the predictions and numerical stability, were studied. It was found that when these effects are taken into account, a good behavior of temperature predictions and numerical stability are obtained. The model presents deviations lower than 14% regarding field data and it presents better predictions than most of the correlations.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅵ)—CONSERVATION LAWS OF MASS AND INERTIA

IntroductionThisworkisadirectcontinuationandasupplementofRefs .[1～8] .InRefs.[1～8]thecoupledbalancelawsandequationsofmomentum ,angularmomentumandenergyaswellasthenewHamiltonprinciple,principleofvirtualpowerandNoethertheoremhavebeenpresented .However,thecoupledconservationlawsofmassandinertiahavenotbeenreestablishedyet.Thepurposeofthispaperistoreestablishtheconservationlawsandequationsofmassandinertiaandtocombinethemwiththecoupledbalancelawsandequationsofmomentum ,angularmomentum ,energyand…     

Characteristics of flow field and water concentration in a horizontal deadleg

Deadlegs are defined as the inactive portion of the pipe where the flow is stagnant. Corrosion in deadlegs occurs as a result of water separation due to the very low flow velocity. The present work provides an investigation of the effect of deadleg geometry and average flow velocity on flow field and oil/water separation in deadlegs. The investigation is based on the solution of the mass and momentum conservation equations of an oil/water mixture together with the volume fraction equation for the secondary phase. A fluid flow model based on the time-averaged governing equations of 3-D turbulent flow has been developed. An algebraic slip mixture model is utilized for the calculation of the two immiscible fluids (water and crude oil). The model solves the continuity and momentum equations for the mixture, and the volume fraction equation for the secondary phase utilizing an algebraic expression for the relative velocity. Flow visualization experiments were conducted in order to validate the numerical procedure. Good agreement was obtained between the calculated and measured flow patterns. Results are obtained for different lengths of the deadleg. The considered fluid mixture contains 90% oil and 10% water (by volume). The inlet flow velocity ranges from 0.2 to 10 m/s and the deadleg length to diameter ratio (L/D) ranges from 1 to 10. The results showed that the size of the stagnant fluid region increases with the increase of L/D and decreases with the increase of inlet velocity. The results also indicated that the water volumetric concentration increases with the increase of L/D and influenced by the deadleg geometry.     

The use of Hugoniot analysis for the propagation of vapor explosion waves

Experimental studies have shown that a mixture of molten metal and water can support the propagation of a quasi-steady vapor explosion wave. Analysis of steadily propagating vapor explosion waves has been carried out by applying the one-dimensional conservation laws of mass, momentum and energy and appropriate equations of state to a homogeneous mixture of molten tin, water and steam. The effects of void fraction, melt/water volume fraction and melt temperature on the Hugoniot curves have been considered. For low temperature melts, the Hugoniot curve lies partially inside the saturation dome and a Chapman-Jouguet (CJ) detonation point occurs only for low void fractions. For high melt temperatures, the downstream states lie entirely outside the saturation region. Increasing the volatility of the coolant or the addition of chemical reactions increases the predicted CJ detonation pressure and velocity. CJ deflagration solutions were obtained in all cases. The existence of a CJ detonation or CJ deflagration for a multiphase fuel-coolant mixture has yet to be substantiated experimentally and nonequilibrium effects may play a role in the divergence between theory and experiments.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

CONCENTRATION OF COUPLED CUBIC NONLINEAR SCHRODINGER EQUATIONS

A coupled nonlinear Schrodinger equations is considered in 2-D space. Based upon the conservation of mass and energy, local identities is established by the study of the limit behavior of the solutions, and L^2-concentration for the blow-up solutions with radially symmetry is obtained.     

Heat and mass transfer in unsaturated porous media with solid–liquid change

 A model on heat and mass transfer in unsaturated porous media with solid/liquid phase change was developed with extending the three-variable model previously proposed. The movement of air phase and its effect on the motion of water is considered. The model has been checked with comparison of the experimental results of the temperature distribution for two dimensional freezing process. The evolution of air pressure, water and ice saturation were predicted by solving the governing equations. The ice segregation and moisture movement toward the front of freezing were numerically simulated. Received on 8 December 1999