成型充填过程的ALE有限元模拟

在ALE框架中提出了一个用于成型充填过程有限元数值模拟的模型。应用ALE参考构形及ALE参考粒子速度描写充填过程中的熔体质量运动。摒弃了Hele-Shaw近似假定,因而所提出的模型能用于非薄壁型腔中高分子材料充填过程的数值模拟。应用基于时域分步算法的Taylor-Galerkin方法,对控制成型充填过程的守恒方程建立了弱形式。对移动自由面附近的充填材料区构造了网格生成算法与网格重划分方案。给出了在几种不同形状的典型腔体中充填过程的数值模拟结果,表明了所提出的ALE有限元模型模拟充填过程的有效性。     

基于计算流体力学的“虚拟飞行”技术及初步应用

以动态混合网格技术为基础, 通过耦合求解刚体动力学方程、流体力学控制方程以及飞行控制律, 建立了适用于飞行器"虚拟飞行" 过程研究的一体化数值模拟技术. 通过典型的外挂物投放算例对流体力学控制方程/动力学方程的耦合算法进行了测试, 并对某导弹的姿态角控制过程、过载控制过程以及变马赫数条件下的控制过程进行了数值模拟, 得到了与实验非常一致的结果. 这些算例证明该一体化算法已经初步具备了针对复杂飞行器"数值虚拟飞行" 的应用能力.     

饱和颗粒土-维冰分凝模型及数值模拟

基于相平衡和力甲衡提出了一个关于饱和颗粒土的冰分凝模型，引入了冰分凝时未冻水迁移的驱动势，基于Clapeyron方程和吸附膜的观点简要叙述了相变区内未冻水含量与温度及孔隙水压的关系，数值模拟了饱和颗粒士的一维冻结过程，得到了冰分层分凝的计算结果，与已有实验所观测到的现象较为吻合，计算所得的冻胀量的变化趋势与已有实验所测结果是相符的．     

饱和土动力学有限元分析的改进稳定分步算法

基于Blot理论，控制饱和变形多孔介质中固相位移u和孔隙压力pw演变的场方程的空间半离散化导致u-pw型混合有限元方程。在固体颗粒和孔隙液体不可压缩以及零渗透性情况下，基本未知量u，pw的近似插值函数必须满足Babuska—Brezzi条件或者与之等价的Zienkiewicz和Taylor分片试验。采用相同低阶u-pw插值的有限元(如线性三角形单元和双线性四边形单元)不能满足B—B条件。分步算法作为一种稳定技术的引入可以绕开B—B条件，但现有分步算法在瞬态问题中仍存在虚假数值振荡和不稳定现象。本文在现有分步算法的基础上引入迭代过程，有效地缓解和克服了数值振荡现象，使低阶u-pw单元得以正常应用。应用双线性四节点u-pw单元的数值结果表明了所提出的包含迭代过程的改进分步算法的有效性。     

植物固沙区土壤水热运移耦合模型研究

根据Ｐｈｉｌｉｐ与Ｖｒｉｅｓ提出的土壤中水热交换的耦合理论，建立了植物固沙区土壤水热运移的耦合模型，考虑了液态水和汽态水运移对温度的变化，分析了植物根系吸水对土壤水热交换的影响，给出了植物蒸腾量、土壤蒸发量、根系汲水率及土壤表面热通量等有关变量的计算公式．利用本模式并采用具有二阶精度的Ｃｒａｎｋ－Ｎｉｃｏｌｓｏｎ格式对非线性扩散方程进行离散，对沙坡头植物固沙区土壤水热交换过程进行了数值模拟，模拟结果和实测值进行了比较，总体上符合较好，证明本模型具有实用价值，可为改造沙漠提供科学依据     

植物固沙区土壤水热运移耦合模型研究

根据Ｐｈｉｌｉｐ与Ｖｒｉｅｓ提出的土壤中水热交换的耦合理论，建立了植物固沙区土壤水热运移的耦合模型，考虑了液态水和汽态水运移对温度的变化，分析了植物根系吸水对土壤水热交换的影响，给出了植物蒸腾量、土壤蒸发量、根系汲水率及土壤表面热通量等有关变量的计算公式．利用本模式并采用具有二阶精度的Ｃｒａｎｋ－Ｎｉｃｏｌｓｏｎ格式对非线性扩散方程进行离散，对沙坡头植物固沙区土壤水热交换过程进行了数值模拟，模拟结果和实测值进行了比较，总体上符合较好，证明本模型具有实用价值，可为改造沙漠提供科学依据     

高速冲击过程数值分析的再生核质点法

利用新型的无网格方法-再生核质点方法对高速冲击过程进行数值模拟，给出了控制 方程，分析中引入Bordner-Partom本构模型来实现材料高速冲击条件下的大应变和高应变 率的特性，提出了新的接触面处理方法，介绍了速度配点法来实现接触面及其它本质边界的 速度条件. 数值模拟表明，再生核质点法能够准确模拟高速冲击过程，具有较高的计算精度.     

除尘管道内铝粉运移及沉积的数值研究

基于DPM(Discrete Phase Model)模型，研究了长直通风管道内粒径服从Rosin-Rammler 分布的铝粉的运移与沉积规律．基于颗粒与壁面的碰撞过程中的能量分析，建立了粉尘沉积-回弹模型，得出了粉尘沉积的判定准则及脱离壁面时的回弹速度．利用UDF 将沉积-回弹模型嵌入Fluent，完成了对管道内粉尘运移和沉积的数值模拟．粉尘沉积的数值结果与实验结果符合得较好，验证了所提模型的有效性．数值结果表明风速的增大使管道内粉尘浓度明显降低，管壁粉尘沉积率也降低；粒径的增大对粉尘浓度的大小影响不明显，主要影响粉尘浓度在管道内的分布情况，同时会增大粉尘在管壁的沉积率．     

可展桁架结构展开过程分析

提出了一种分析构架式结构展开过程的有效算法。基于含多余广义坐标的动力学普遍方程 ,利用约束雅可比矩阵的零空间基引入一组准速率 ,得到独立的展开过程分析的动力学微分方程。为提高展开模拟的数值精度 ,文中提出了一种控制展开过程几何违约、速度违约和能量违约的数值稳定算法。该算法求解效率高 ,能和任意数值积分方法结合使用 ,能分析大型的构架式可展结构的展开过程     

高温下混凝土化学塑性-损伤耦合本构模拟及破坏分析

提出了一个用于模拟高温下混凝土中化学-热-湿-力学耦合行为的化学塑性-损伤耦合本构模型.发展了用于积分率形式的耦合本构方程的三步算子分裂算法.导出了化学-热-湿-力学(CTHM)耦合本构模型的一致性切线模量矩阵,以保证对于全局耦合控制方程Newton迭代过程的二阶收敛率.算例数值结果显示了本文发展的化学塑性-损伤耦合本构模型在重现火灾和热辐射条件下的混凝土中化学-热-湿-力学耦合行为的有效性.     

An Approach to Transport in Heterogeneous Porous Media Using the Truncated Temporal Moment Equations: Theory and Numerical Validation

In the last decade, the characterization of transport in porous media has benefited largely from numerical advances in applied mathematics and from the increasing power of computers. However, the resolution of a transport problem often remains cumbersome, mostly because of the time-dependence of the equations and the numerical stability constraints imposed by their discretization. To avoid these difficulties, another approach is proposed based on the calculation of the temporal moments of a curve of concentration versus time. The transformation into the Laplace domain of the transport equations makes it possible to develop partial derivative equations for the calculation of complete moments or truncated moments between two finite times, and for any point of a bounded domain. The temporal moment equations are stationary equations, independent of time, and with weaker constraints on their stability and diffusion errors compared to the classical advection–dispersion equation, even with simple discrete numerical schemes. Following the complete theoretical development of these equations, they are compared firstly with analytical solutions for simple cases of transport and secondly with a well-performing transport model for advective–dispersive transport in a heterogeneous medium with rate-limited mass transfer between the free water and an immobile phase. Temporal moment equations have a common parametrization with transport equations in terms of their parameters and their spatial distribution on a grid of discretization. Therefore, they can be used to replace the transport equations and thus accelerate the achievement of studies in which a large number of simulations must be carried out, such as the inverse problem conditioned with transport data or for forecasting pollution hazards.     

Robust Multi-D Transport Schemes with Reduced Grid Orientation Effects

In this paper we investigate truly multi-D upwind schemes for simulating adverse mobility ratio displacements in porous media. Due to an underlying physical instability at the simulation scale, numerical results are highly sensitive to discretization errors and hence the orientation of the underlying computational grid. We use modified equations analysis to predict preferred flow angles on structured grids for several popular methods and present a conservative, multi-D framework for designing positive upwind schemes for general velocity fields. After placing the common schemes in this framework, we go on to develop a novel scheme with “minimal” constant transverse (cross-wind) diffusion. Results for miscible gas injection into homogeneous and heterogeneous media demonstrate that truly multi-D schemes, and in particular our new scheme, greatly reduce grid orientation effects and numerical biasing as compared to dimensional upwinding.     

A dynamic coupling numerical model for pollutant transport under wave condition

In order to study the diffusion, migration, and distribution of pollutants among overlying water-body and porous seabed under wave conditions, a dynamic coupling numerical model is proposed. In this model, the coupling between wave field of overlying water-body and seepage of porous bed, the capture and release of pollutants in porous media, and the transport process between the two different regions are taken into account. We use the unified equations for pressure correction and pollutant concentration to solve the numerical model, which avoids repeated iteration on the interface boundary. The model is verified by several case studies. Afterwards, the processes involving release of pollutant from porous seabed and transportation to overlying water-body under different wave conditions are investigated. The results show that the water depth, wave height,and wave period have great influences on the release, capture, and transport processes for phosphorus pollutant.     

Coupling between shallow water and solute flow equations: analysis and management of source terms in 2D

A two‐dimensional model for the simulation of solute transport by convection and diffusion into shallow water flow over variable bottom is presented. It is based on a finite volume method over triangular unstructured grids. A first order upwind technique is applied to solve the flux terms in both the flow and solute equations and the bed slope source terms and a centred discretization is applied to the diffusion and friction terms. The convenience of considering the fully coupled system of equations is indicated and the methodology is well explained. Three options are suggested and compared in order to deal with the diffusion terms. Some comparisons are carried out in order to show the performance in terms of accuracy and computational effort of the different options. Copyright © 2005 John Wiley & Sons, Ltd.     

Solution of the 2‐D shallow water equations with source terms in surface elevation splitting form

A vertex‐centred finite‐volume/finite‐element method (FV/FEM) is developed for solving 2‐D shallow water equations (SWEs) with source terms written in a surface elevation splitting form, which balances the flux gradients and source terms. The method is implemented on unstructured grids and the numerical scheme is based on a second‐order MUSCL‐like upwind Godunov FV discretization for inviscid fluxes and a classical Galerkin FE discretization for the viscous gradients and source terms. The main advantages are: (1) the discretization of SWE written in surface elevation splitting form satisfies the exact conservation property (??‐Property) naturally; (2) the simple centred‐type discretization can be used for the source terms; (3) the method is suitable for both steady and unsteady shallow water problems; and (4) complex topography can be handled based on unstructured grids. The accuracy of the method was verified for both steady and unsteady problems, including discontinuous cases. The results indicate that the new method is accurate, simple, and robust. Copyright © 2007 John Wiley & Sons, Ltd.     

Analysis of a second‐order upwind method for the simulation of solute transport in 2D shallow water flow

A two‐dimensional model for the simulation of solute transport by convection and diffusion into shallow water flow over variable bottom is presented. It is based on a finite volume method over triangular unstructured grids. A first‐order upwind technique, a second order in space and time and an extended first‐order method are applied to solve the non‐diffusive terms in both the flow and solute equations and a centred implicit discretization is applied to the diffusion terms. The stability constraints are studied and the form to avoid oscillatory results in the solute concentration in the presence of complex flow situations is detailed. Some comparisons are carried out in order to show the performance in terms of accuracy of the different options. Copyright © 2007 John Wiley & Sons, Ltd.     

A coupled lattice Boltzmann model for the shallow water‐contamination system

A coupled numerical method for the direct simulation of shallow water dynamics and pollutant transport is formulated and implemented. The conservation equations of shallow water dynamics equations and the convection–diffusion equations are solved using the lattice Boltzmann (LB) method. The local equilibrium distribution of the pollutant has no terms of second order in flow velocity. And the relaxation time of the pollutant deviates from a constant for the flows with variable free surface water depth. The numerical tests show that this scheme strictly obeys the conservation law of mass and momentum. Excellent agreement is obtained between numerical predictions and analytical solutions in the pure diffusion problem and convection–diffusion problem. Furthermore, the influences on the accuracy of the lattice size and the diffusivity are also studied. The results indicate that the variation in the free surface water depth cannot affect the conservation of the model, and the model has the ability to simulate the complex topography problem. The comparison shows that the LB scheme has the capacity to solve the complex convection–diffusion problem in shallow water. Copyright © 2008 John Wiley & Sons, Ltd.     

Co-located equal-order control-volume finite element method for two-dimensional axisymmetric incompressible fluid flow

The formulation of a control-volume-based finite element method (CVFEM) for axisymmetric, two-dimensional, incompressible fluid flow and heat transfer in irregular-shaped domains is presented. The calculation domain is discretized into torus-shaped elements and control volumes. In a longitudinal cross-sectional plane, these elements are three-node triangles, and the control volumes are polygons obtained by joining the centroids of the three-node triangles to the mid-points of the sides. Two different interpolation schemes are proposed for the scalar-dependent variables in the advection terms: a flow-oriented upwind function, and a mass-weighted upwind function that guarantees that the discretized advection terms contribute positively to the coefficients in the discretized equations. In the discretization of diffusion transport terms, the dependent variables are interpolated linearly. An iterative sequential variable adjustment algorithm is used to solve the discretized equations for the velocity components, pressure and other scalar-dependent variables of interest. The capabilities of the proposed CVFEM are demonstrated by its application to four different example problems. The numerical solutions are compared with the results of independent numerical and experimental investigations. These comparisons are quite encouraging.     

极坐标系下Fourier-Legendre谱元方法与有限差分法数值扩散的比较

提出一种Fourier-Legendre谱元方法用于求解极坐标系下的Navier-Stokes方程,其中极点所在单元的径向采用Gauss-Radau积分点,避免了r=0处的1/r坐标奇异性。时间离散采用时间分裂法,引入数值同位素模型跟踪同位素的输运过程验证数值模拟的精度,分别利用谱元法和有限差分法的迎风差分格式求解匀速和加速坩埚旋转流动中的同位素方程。计算结果表明,有限差分法中的一阶迎风差分格式存在严重的数值假扩散,二阶迎风差分格式的数值结果较精确,增加节点可以有效地缓解数值扩散。然而,谱元法具有以较少节点得到高精度解的优势。     

A well‐balanced weighted essentially non‐oscillatory scheme for pollutant transport in shallow water

In this paper, a well‐balanced finite difference weighted essentially non‐oscillatory scheme is presented for modeling transport and diffusion of pollutant in shallow water flows. The scheme balances exactly the flux gradients and the source terms. Extensive one‐dimensional and two‐dimensional numerical experiments on uniform and curvilinear meshes strongly suggest that high resolution results are achieved for both water depth and pollutant concentration. The scheme is efficient and robust and can be applied to practical numerical simulation of pollutant transport phenomena in shallow water flows. Copyright © 2012 John Wiley & Sons, Ltd.