对流扩散方程的迎风变换及相应有限差分方法

本文提出所谓迎风变换,将对流扩散方程分解为对流迎风函数和扩散方程,并构造相应的有限差分格式。对流迎风函数以简明的指数解析形式反映对流扩散现象的迎风效应,原则上消除了源于不对称对流算子的困难,能够便利对流扩散方程的数值求解。有限差分格式具有二阶精度和无条件稳定性,算例表明其准确性、收敛速度及对边界层效应的适应能力均明显优于中心差分格式和迎风差分格式。     

二维对流扩散方程的欧拉—拉格朗日分裂格式

本文在[1]基础上发展了一种有效的处理大P_e(R_e)数、非定常二维对流扩散方程的欧拉-拉格朗日(E-L)分裂格式,由于方法本质上与区域形状无关,且不需再分网格,因此是一种无网格的E-L方法,特别对于定常流动,E.-L.分裂格式可以导致比一阶迎风格式更精确的单调、无振荡格式,文中对于常系数、变系数和非线性的二维非定常和定常对流扩散方程的(初)边值问题进行了数值计算,数值结果与精确解的比较表明,本方法具有很好的精度,解是单调无振荡的,比通常一阶迎风格式具有较少的数值扩散,最大计算网格P-e(R-e)数可达100—500。     

在蒸发条件下土壤水盐运动的数值模拟

本研究对土壤中水盐运动对流扩散方程进行了数值计算。利用省时的一阶精度的特征-分步法及使对流项作保留迎风特性的二阶精度的差分法,对在蒸发条件下土壤盐分向上运动积盐过程进行模拟,两个方法的计算结果与实验结果总体上都重合较好。     

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

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

求解温盐双扩散系统的一种高精度方法

引入边界拟合坐标系来研究温盐双扩散系统。为了提高求解的精确性 ,对流项采用四阶精度的迎风格式 ,扩散项和涡量方程的浮力项则采用四阶精度的中心差 ,因此本文的方法是高精度的方法。首先针对温度占优( Rρ=0 .32 )和盐度占优 ( Rρ=1 .6 8)的情形进行了验证性计算 ,得到了与前人一致的结果。进一步 ,本文系统给出了不同的盐度通量强度下的流动形态 ,包括对称结构 ,不对称结构 ,反转结构等 ,结果与前人的吻合。     

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

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

高阶多维半离散中心迎风格式及其应用

提出了求解多维对流-扩散方程的四阶半离散中心迎风格式。该格式以中心加权基本无振荡(CWENO)重构为基础,同时考虑到在Riemann扇内波传播的局部速度,从而更加准确地估计出了局部Riemann扇的宽度,最终既回避了网格的交错,又降低了格式的数值粘性,建立了介于迎风格式和中心格式之间的半离散中心迎风格式。本文还将该四阶半离散中心迎风格式与涡度-流函数方法相结合,有效地求解了二维不可压Euler方程组和Navier-Stokes方程组。     

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

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

一种新的求解对流占优问题的自适应网格细化方法

在均匀网格上求解对流占优问题时,往往会产生数值震荡现象,因此需要局部加密网格来提高解的精度。针对对流占优问题,设计了一种新的自适应网格细化算法。该方法采用流线迎风SUPG(Petrov-Galerkin)格式求解对流占优问题,定义了网格尺寸并通过后验误差估计子修正来指导自适应网格细化,以泡泡型局部网格生成算法BLMG为网格生成器,通过模拟泡泡在区域中的运动得到了高质量的点集。与其他自适应网格细化方法相比,该方法可在同一框架内实现网格的细化和粗化,同时在所有细化层得到了高质量的网格。数值算例结果表明,该方法在求解对流占优问题时具有更高的数值精度和更好的收敛性。     

对流扩散方程的绝对稳定高阶中心差分格式

将作者提出的数值摄动算法改进为区分离散单元内上游和下游并分别对通量进行高精度重构的双重数值摄动算法,与原(单重)摄动算法相比,双重摄动算法既提高了格式精度又明显扩大了格式的稳定域范围.利用双重摄动算法,即分别利用上游和下游基点变量的摄动重构将高阶流体力学关系及迎风机制耦合进二阶中心格式之中,由此构建了对流扩散方程的对网格Reynolds数的任意值均稳定(绝对稳定)高精度(四阶和八阶精度)三基点中心TVD差分格式,通过解析分析以及3个算例计算证实了构建格式的优良性能;3个算例包括一维线性、非线性(Burgers方程)和二维变系数对流扩散方程.数值计算表明:构建的格式在粗网格下不振荡,构建格式在粗网格时的最大误差L_∞和均方误差L_2与二阶中心格式在细网格时的相应误差一致,对线性方程,构建格式在细网格下可达到L_2精度阶.     

线性定常对流占优对流扩散问题的无网格解法

应用无网格Galerkin方法求解对流占优对流扩散问题时会出现非物理现象的数值伪振荡,本文将SUPG方法、GLS方法、SGS方法与无网格Galerkin方法相耦合,成功解决了对流扩散方程中对流项占优时的数值伪振荡问题。运用本文构造的方法,采用线性基和具有C2连续的权函数,应用移动最小二乘法可容易地构造高阶导数连续的形函数,从而避免了有限元方法中当采用线性元插值时,因忽略稳定项中二阶导数项而降低计算精度和稳定性的问题。数值实验表明:本文构造的方法具有计算精度高、稳定性好、计算算法实施简单、前后处理方便的优点,这些方法不仅能适用于对流项占优问题,而且也能很好地消除反应项占优时的数值伪振荡问题。     

Managing false diffusion during second‐order upwind simulations of liquid micromixing

Liquid mixing is an important component of many microfluidic concepts and devices, and computational fluid dynamics (CFD) is playing a key role in their development and optimization. Because liquid mass diffusivities can be quite small, CFD simulation of liquid micromixing can over predict the degree of mixing unless numerical (or false) diffusion is properly controlled. Unfortunately, the false diffusion behavior of higher‐order finite volume schemes, which are often used for such simulations, is not well understood, especially on unstructured meshes. To examine and quantify the amount of false diffusion associated with the often recommended and versatile second‐order upwind method, a series of numerical simulations was conducted using a standardized two‐dimensional test problem on both structured and unstructured meshes. This enabled quantification of an ‘effective’ false diffusion coefficient (D_false) for the method as a function of mesh spacing. Based on the results of these simulations, expressions were developed for estimating the spacing required to reduce D_false to some desired (low) level. These expressions, together with additional insights from the standardized test problem and findings from other researchers, were then incorporated into a procedure for managing false diffusion when simulating steady, liquid micromixing. To demonstrate its utility, the procedure was applied to simulate flow and mixing within a representative micromixer geometry using both unstructured (triangular) and structured meshes. Copyright © 2016 John Wiley & Sons, Ltd.     

饱和-非饱和土壤中污染物运移过程的数值模拟

本文提出了一个模拟饱和 非饱和土壤中溶和污染物运移过程的数值模型．模拟的控制污染物运移的物理 化学现象包括：对流，机械逸散，分子弥散，吸附，蜕变，不动水效应．发展了一个修正的特征线Ｇａｌｅｒｋｉｎ方法以离散污染物运移过程的控制方程并导出了一个用于有限元方程求解的显式算法．数值例题结果表明所提出模型和算法的功能     

超音速粘性流动的SUPG有限元数值解法

本文构造了准简化 N-S 方程组的 SUPG(Streamline Upwind/Petrov-Galerk-in)加权剩余式,并利用该方法对 Burgers 方程、无粘性激波反射问题、以及超音速平板和压缩拐角的层流流动作了数值求解。计算结果表明,本文方法是精确、收敛和稳定的。     

A MESH-FREE ALGORITHM FOR DYNAMIC IMPACT ANALYSIS OF HYPERELASTICITY

A mesh-free method based on local Petrov-Galerkin formulation is presented to solve dynamic impact problems of hyperelastic material.In the present method,a simple Heaviside test function is chosen for simplifying domain integrals.Trial function is constructed by using a radial basis function(RBF)coupled with a polynomial basis function,in which the shape function possesses the kronecker delta function property.So,additional treatment is not required for imposing essential boundary conditions.Governing equations of impact problems are established and solved node by node by using an explicit time integration algorithm in a local domain,which is very similar to that of the collocation method except that numerical integration can be implemented over local domain in the present method.Numerical results for several examples show that the present method performs well in dealing with the dynamic impact problem of hyperelastic material.     

A comparative study of central and upwind difference schemes using the primitive variables

The use of the velocity-pressure formulation of the Navier-Stokes equations for the numerical solution of fluid flow problems is favoured for free-surface problems, more involved flow configurations, and three-dimensional flows. Many engineering problems involve such features in addition to strong inertial effects. The computational instabilities arising from central-difference schemes for the convective terms of the governing equations impose serious limitations on the range of Reynolds numbers that can be investigated by the numerical method. Solutions for higher Reynolds numbers Re > 1000 could be reached using upwind-difference schemes. A comparative study of both schemes using a method based on the primitive variables is presented. The comparison is made for the model problem of the driven flow in a square cavity. Using a central scheme stable solutions of the pressure and velocity fields were obtained for Reynolds numbers up to 5000. The streamfunction and vorticity fields were calculated from the resulting velocity field and compared with previous solutions. It is concluded that total upwind differencing results in a considerable change in the flow pattern due to the false diffusion. For practical calculations, by a proper choice of a small amount of partial upwind differencing the vorticity diffusion near the walls and the global features of the solutions are not sigificantly altered.     

Stabilization meshless method for convection-dominated problems

It is weN-known that the standard Galerkin is not ideally suited to deal with the spatial discretization of convection-dominated problems. In this paper, several techniques are proposed to overcome the instabilitY issues in convection-dominated problems in the simulation with a meshless method. These stable techniques included nodal refinement, enlargement of the nodal influence domain, full upwind meshless technique and adaptive upwind meshless technique. Numerical results for sample problems show that these techniques are effective in solving convection-dominated problems, and the adaptive upwind meshless technique is the most effective method of all.     

Application of streamwise diffusion to time-dependent free convection of liquid metals

A numerical analysis is given for the application of streamwise diffusion to high-intensity flows with marginal spatial resolution. Terms are added to the momentum equation which are similar to those used in the Petrov-Galerkin, Taylor-Galerkin and balancing tensor diffusivity methods. Values for the streamwise viscosity are obtained from mesh refinement studies. An illustration is given for the time-dependent free convection of a liquid metal in a cavity with differentially heated sided walls. The spatial problem is solved with the Galerkin finite element method and the time integration is performed with the backward Euler method. Solution quality and computation time are compared for results with and without added streamwise diffusion. For the cases presented, streamwise diffusion eliminates spurious oscillations and saves computation time without compromising the solution.