对流项占优问题的MLPG/SUPG方法数值模拟

在计算对流项占优问题时易产生假扩散,本文把流线型迎风格式应用于MLPG方法中可以减少对流项的影响,通过两个典型例子（旋转流场问题和Brezzi问题）验证该格式的精度与有效性,并与文献中的迎风格式的计算结果进行比较,计算结果表明,该方法能有效地克服假扩散现象,有较好的稳定性和较高的计算精度。     

用物理黏性构建高阶不振荡对流扩散差分格式

利用数值摄动算法, 通过扩散格式数值摄动重构把对流扩散方程的2阶中心差分格式(2-CDS)重构为高精度高分辨率格式, 解析分析和模型方程计算证实了新格式的高精度不振荡性质. 新格式是把物理黏性使流动光滑化的扩散运动规律引入2-CDS 中的结果. 该法显然与构建高级离散格式的常见方法不同. 证实: 数值摄动重构中引入扩散运动规律的结果格式与引入对流运动规律(下游不影响上游的规律)的结果格式一致, 说明对离散方程的数值摄动运算, 在维持原格式结构形式不动的条件下, 不仅能提高格式精度和稳健性, 且可揭示对流离散运动规律与扩散离散运动规律之间的内在关联;同时证实, 文中提出和使用的上、下游分裂方法是构建高精度不振荡离散格式的一个有效方法.     

饱和非饱和土壤中溶质运移的数值模拟

本文结合拉格朗日观点的运动坐标与欧拉观点的固定坐标,分别求解对流—弥敌方程的对流效应和弥散效应,并在从固定坐标向运动坐标投影时利用了三次样条插值方法,使模型更为精确而简单。通过与特定情况下的理论解比较,说明模型是准确可靠的。通过与实验室室内试验结果的比较,模型的可靠性得到了进一步的验证。利用模型对淡水洗盐结合浅群井抽排咸水治理盐碱地过程的模拟分析,认为结合浅群井抽排可较大地提高淡水洗盐的效果。     

无量纲混相驱对流-弥散方程的数值求解

给出了混相驱对流-弥散方程的无量纲形式；采用高阶精度差分格式对其进行数值求解并应用于几个具体算例中。本文数值结果和算例给出的结果对比，充分验证了本文数值模型的可靠性和准确性。该数值模型可以用来研究诸如土壤或地下水污染物输运以及空气污染等领域里的传质问题．     

入水问题的Euler-Lagrange流固耦合数值模以技术研究

针对结构入水过程的特点,建立了一种把Eulerian计算过程和Lagrangian计算过程耦合起来的数值技术,编写了计算程序,把结构对水作用的过程和结构动力学响应过程的计算结合起来,实现了结构与水之间流固耦合作用过程的模拟.其中水、空气的运动利用Eulerian过程的流体弹塑性有限差分计算程序进行模拟,结构的运动变形利用Lagrangian过程的结构动力学有限元计算程序模拟结构响应.通过板入水问题的计算结果与商业软件计算结果进行比较,验证了计算方法对流固间相互作用处理的正确性.本文计算程序可用于一般入水问题的结构响应-流场运动分析中.     

非定常流函数涡量方程的一种数值解法的研究

对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解．以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数（Ra）条件下的数值研究．结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟．     

高精度边界格式的研究

利用有精确解的Ｒｉｎｇｌｅｂ流动，构造了对流场数值解精度进行检验的“Ｒｉｎｇｌｅｂ机器”．重点讨论了边界格式对流场数值解的影响及高精度边界格式的建立．计算表明，在场内应用二阶精度格式情况下，采用二阶精度的边界格式所得到的流场解精度将大大高于采用一阶边界格式所得到的精度．为提高流场解的精度，不仅需要高精度的边界格式，还必须注意边界格式与场内格式的匹配．计算还表明，采用特征线修正的方法能有效地提高边界处理的精度     

一种基于Associated Hermite正交函数求解对流扩散方程的算法

提出了一种基于AH(Associated Hermite)正交基函数求解对流扩散方程的无条件稳定算法。该算法将方程的时间项通过Hermite多项式作为正交基函数进行展开,利用Galerkin方法消除时间变量项,从而导出有限维AH域隐式差分方程,突破了传统显式差分格式稳定性条件的限制,最后通过对AH域展开系数的求解得到该对流扩散方程的数值解。在数值算例中,将该算法与传统显示差分法和交替方向隐式差分法进行对比分析,数值计算结果表明,算法无条件稳定且其计算精度与时间步长无关,对于具有精细结构的对流换热问题,该算法具有明显的效率优势,且保持了较高的精度。     

膛口反应流并行数值模拟

采用轴对称多组分N-S方程对含有高速运动弹丸的膛口反应流进行了数值模拟.控制方程采用时间分裂方法并在大型计算机上采用MPI方法进行多核并行求解,其中对流项采用二阶AUSM+格式和MUSCL插值方法进行处理,燃气采用氢气-空气混合气,反应机理为9组分19步基元反应.对于弹丸引起的网格运动,采用嵌套网格法处理.并行验证算例与串行计算结果一致,采用20个CPU计算时效率为64％.根据数值结果详细讨论了发射过程中的气体动力学和化学动力学过程,并且通过对两种条件下的计算结果比较分析了化学反应对膛口流场发展的影响.结果表明,上述算法能够较为正确地模拟弹丸和化学反应对膛口流场的影响,并大大提高了计算速度.     

竖直流道宽度对气泡运动行为影响的数值模拟

用数值方法模拟了竖直通道宽度对气泡在液体中的非定常运动、变形以及传热特性的影响。在这个模拟中,界面跟踪采用了VOF方法,并采用PL IC进行界面重构。主流场计算采用有限容积方法将控制方程离散,其中扩散项采用中心差分格式,对流项采用一阶迎风格式。用成熟的S IM PLE算法求解N-S方程的速度与压力的耦合问题。引入CSF模型处理运动界面的表面张力。利用所编制的程序计算了竖直流道中的单个气泡的形状、运动特性以及气泡内外流场与传热特性,并对竖直通道宽度在不同情况下,对气泡的形状、运动特性以及传热特性进行了进一步的研究。得到了一系列有价值的结果,并与实验结果比较。表明数值模拟结果与实验结果吻合的较好。     

含离子类溶质黏土的 C-H-M 本构模型简

在水力-力学本构模型基础上,采用化学软化公式引入孔隙水中含盐、碱类等离子类溶质对土体前期固结压力的影响,得到土体的化学-水力-力学耦合本构模型（C-H-M 模型）. 根据试验结果对化学软化公式的参数进行数值拟合,通过算例验证了C-H-M 本构模型的有效性     

含离子类溶质黏土的C-H-M本构模型

在水力-力学本构模型基础上,采用化学软化公式引入孔隙水中含盐、碱类等离子类溶质对土体前期固结压力的影响,得到土体的化学-水力-力学耦合本构模型（C-H-M 模型）. 根据试验结果对化学软化公式的参数进行数值拟合,通过算例验证了C-H-M 本构模型的有效性     

Regimes of haline convection during the evaporation of groundwater containing a dissolved admixture

The problem of formation of salt concentration profile in high-permeability soil duringwater evaporation and solution upflow is considered. The numerical experiments performed showed that the salt concentration profile may be either stable or unstable. As instability develops, there arises natural haline convection whose different regimes are described and analyzed. If the evaporation intensity is moderate, in soil the curvilinear upward or circulatory flow that fills the entire layer is established. The intense evaporation leads to the formation of a small-scale structure of salt “fingers”. Boundaries between regimes are determined.     

盐渍土工程与力学性质研究进展

介绍我国西部地区大开发高等级公路、铁路网建设中盐渍土病害的力学相关问题. 综述了对盐胀、溶陷和腐蚀3大工程危害的机理研究以及对盐渍土的工程及力学特性的研究进展, 回顾了盐渍土水、盐、热的数学运移模型的发展, 盐渍土水、盐、热、力4场耦合模型的数值与实验探索, 对盐渍土研究领域中亟待进一步研究的问题进行了展望.     

含盐量对水泥土强度影响的室内试验研究

通过含盐量对非有机质土加固强度影响的试验研究,得到了含盐量对水泥土强度的提高或减小的阈值为3.5%。当盐渍土的含盐量低于这个阈值时,盐渍土的加固强度会因可溶性盐的结晶膨胀作用,提高水泥土的强度;相反当盐渍土的含盐量高于该阈值时,盐渍土的强度会因可溶性盐的过多的结晶膨胀作用,使水泥土的结构遭到破坏,从而使水泥土的强度大大降低。同时分析了可溶性硫酸盐、镁盐和氯盐对水泥土的浸蚀性作用,并从盐类对水泥土强度的影响从机理上进行了阐释,提出了高含盐量对水泥土破坏作用的对策。     

Heterogeneity Effects on Evaporation-Induced Halite and Gypsum Co-precipitation in Porous Media

Understanding evaporation from porous media in the presence of soluble salts plays a key role in describing many environmental processes. Several studies done in this field led to a wide acceptance that the prediction of soil salinization driven by evaporation from the unsaturated zone and mainly encountered in arid and semiarid regions is a challenging task because of the temporal and spatial variabilities of soil rocks combined with different interactions between the porous medium and the atmosphere. In this work, we present a reactive transport model developed with the aim of describing the processes of evaporation, salt accumulation and precipitation. We took the model presented in our previous paper (Jambhekar et al. in Transp Porous Media 114:341–369, 2016) as the basis and developed it with the required geochemical model to account for evaporative salt co-precipitation. The salts considered in this work are halite (NaCl) and gypsum (CaSO\(_{4}\cdot 2\mathrm{H}_{2}\)O). We focus particularly on the influence of spatial heterogeneities in the porous medium on the dynamics of the physical processes. In the numerical simulations performed in this work, we distinguished different heterogeneity configurations. The results show that the drying from heterogeneous porous media initially affects the coarser pores, where the salt crystals first appear.     

Water and salt movements in simultaneous flood-irrigation and well-drainage operations

This paper describes a new technology for solonchak soil reclamation in which surface flood irrigation of fresh water and pumped wells drainage of salty groundwater are combined. The comprehensive investigation of water and salt movement has been conducted through field test, laboratory simulation and numerical calculation. The dependence of desalination on irrigation water quantity, drainage quantity, leaching time and other parameters is obtained based on the field tests. The entire desalination process under the flood-irrigation and well-drainage operations was experimentally simulated in a vertical soil column. The water and salt movement has been numerically analyzed for both the field and laboratory conditions. The present work indicates that this new technology can greatly improve the effects of desalination. The project supported by the National Natural Science Foundation of China     

Approximate models of salt transfer with washings of soils with water-impermeable base

On the basis of an approximate model of filtration with drainage washings of a soil layer with a water-impermeable base [1], shown schematically in Fig. 1, an investigation is made of the desalinization process of soils for several forms of the equations of salt transfer between the moving solution and the fixed phase. Neglect of the external diffusion makes it possible to reduce the investigation to the solution of one-dimensional boundary-value problems, characterizing salt transfer along fixed streamlines. Analogous problems were discussed in [2] with application to the one-dimensional filtration of wash waters in a dry soil.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 44–51, January–February, 1978.The author thanks V. N. Émikh and V. I. Pen'kovskii for their repeated and useful evaluations of the work.     

Numerical simulation of convective flows in a soil during the evaporation of water containing a dissolved admixture

The concentration profile is investigated for a solution that saturates a low-permeability soil. The simulation results showed the presence of three flow regimes. The salt accumulated near the phase transition boundary increases the solution density and may lead to the development of natural concentration-induced convection which interacts with the rising flow (forced convection). The stability threshold of the forced flow and the effect on it of natural convection that arises are determined. It is shown that at intense flow to the evaporation surface the admixture concentration increases at this boundary rapidly and reaches the saturation concentration. In this case, the admixture precipitates. In the slow evaporation regime the admixture diffuses from the high-concentration region, which prevents the development of convective flow.