基于无结构网格单元中心有限体积法的二维对流扩散方程离散

在无结构网格单元中心有限体积法二维水流模型基础上,建立物质输运对流扩散方程离散模式.通过通量重构法和SOM(Support Operators Method),分别对输运方程的对流项和扩散项进行离散.该离散模式具有空间二阶精度,并适用于任意多边形无结构网格.通过纯对流和纯扩散算例对模型进行检验和验证,结果表明,模型能够较好地模拟物质输运的对流扩散问题.应用模型模拟瓯江河口的盐度输运,通过计算值与实测值对比,进一步检验模型.     

积分节点迎风偏移的节点积分无单元Galerkin方法

建立求解稳态对流-扩散方程的-种稳定、高效的无单元Galerkin方法.该方法计算积分时采用基于局部Taylor展开的节点积分,并根据对流占优的程度对积分节点进行自适应迎风偏移.与传统的使用稳定化的无单元Galerkin方法相比,该方法是-种不依赖于背景网格积分的纯无网格方法,具有更好的稳定性和较高的计算效率,其程序实施更为简便.     

二维对流扩散方程的格子Boltzmann方法

给出了二维对流扩散方程的格子Boltzmann方法,用对流扩散方程中的对流系数和扩散系数确定了局哉平衡分布函数Chapman-Enskog展开中的可调参数,并对该方法进行了讨论. 关键词：     

改进的节块展开方法求解圆柱几何对流扩散方程

为高效求解球床高温气冷堆物理-热工耦合问题,发展改进节块展开法求解圆柱几何下的对流扩散方程.针对圆柱几何和对流扩散方程的特殊性,采用三阶多项式和指数函数作为r向横向积分方程的展开函数,在节块展开法的框架下高效求解对流扩散方程.数值验证表明,改进的节块展开方法具有固有的迎风特性,在使用粗网节块时依然能保持稳定性和较高的计算精度.     

求解对流扩散方程的紧致修正方法

提出了求解对流扩散方程的紧致修正方法,该方法是在低阶离散格式的源项中,引入紧致修正项,从而构造高阶紧致修正格式,并进行求解.采用紧致修正方法对典型的对流扩散方程进行计算.结果表明,紧致修正方法虽然与二阶经典差分方法建立在相同的结点数上,但紧致修正方法的精度与紧致方法的精度相同,均具有四阶精度.所以紧致修正方法可以在少网...     

基于微可压液体的统一对流扩散型流体力学方程组

为了表达上的方便及求解格式的统一,通常采用统一的方程形式来表达连续方程,动量方程、能量方程、湍动能方程和耗散方程等.除了连续方程外,其他方程都可以写成对流扩散方程的形式,由于没有扩散项,连续方程比较特别,也相对不便处理.在微可压液体区,通过合理的数学推导,不作任何近似、假定与简化,本文得到一套全新的连续方程形式.该新方程以压力为未知变量,是对流扩散型的,使得所有的流体动力学方程组都具有完全统一的方程形式.     

统一的对流扩散型可压缩流体力学方程与解法

流体力学的动量方程、能量方程、湍动能方程和耗散方程都具有对流扩散方程的形式,但连续方程却不是对流扩散型的。对于可压缩问题,本文通过合理的数学推导,不作任何近似、假定与简化,得到一个全新的连续方程形式．该连续方程以压力为未知变量,并具有对流扩散型形式,使得所有的流体动力学方程组都具有完全统一的方程形式,给出了这种三维对流扩散方程组的有限精确差分计算格式。对流体力学的进一步发展具有一定意义．     

对流作用下枝晶形貌演化的数值模拟和实验研究

合金凝固过程中存在于枝晶尖端液相区的强制对流和自然对流均能改变溶质扩散层厚度,从而会对枝晶形貌产生较大影响.在元胞自动机模型基础上,耦合液体流动方程、热传导方程和溶质对流扩散方程,建立了新的计算微观组织演化的数值模型,并利用该模型研究了强制对流和自然对流对枝晶生长的影响.三维数值模拟结果再现了强制对流作用下等轴枝晶的生长过程,揭示了强制对流对枝晶生长速率和尖端半径的影响特点.同时利用该模型模拟了NH₄Cl-H₂O溶液定向凝固过程中自然对流对柱状晶生长的影响,并采用相应的实验进行验证.模拟结果与实验结果符合良好,从而证明该模型是可靠的,可推广到实际合金系中. 关键词： 元胞自动机 对流 4Cl-H₂O溶液')" href="#">NH₄Cl-H₂O溶液 定向凝固     

反常扩散与分数阶对流-扩散方程

反常扩散现象在自然界和社会系统中广泛存在.考虑了扩散过程的时间相关和时空相关性，用非局域性的处理方法，在传统的二阶对流 扩散方程基础上，得到了分数阶对流 扩散方程，以此方程来描述反常扩散.在此方程中，弥散项和对时间的导数为分数阶导数所代替.由此分数阶对流 扩散方程，对传统的费克扩散定律进行推广，得到了广义的分数费克扩散定律，分数费克扩散定律说明某时刻空间中某点的流量不仅与其领域内的浓度梯度有关，而且与整个空间中其他不同点的粒子浓度、浓度变化的历史，甚至初始时刻的浓度有关.讨论了方程的解——分数稳定分布，并由此说明了扩散运动的平均平方位移是运移时间的非线性函数. 关键词： 扩散 分数阶微积分 稳定分布（Lévy分布） 费克扩散定律     

求解对流扩散方程的四种差分格式的比较

利用对流扩散方程，在边界和参数存在随机扰动的情况下，考察四种差分格式的优劣，为求 解对流扩散方程提供一种可靠的差分格式，并得到通过空间加密网格的方法可以控制边界、 参数随机影响的结论. 关键词： 对流扩散方程 差分格式 随机扰动     

三维定常问题的高阶紧致差分格式向非定常问题的推广

将已经建立的求解三维定常对流扩散方程的高阶紧致差分格式直接推广到三维非定常对流扩散方程的数值求解,时间导数项利用二阶向后欧拉差分公式,所得到的高阶隐式紧致差分格式时间为二阶精度,空间为四阶精度,并且是无条件稳定的.数值实验结果验证了本文方法的精确性和稳健性.     

SUPG-Stabilized Virtual Element Method for Optimal Control Problem Governed by a Convection Dominated Diffusion Equation

In this paper, the streamline upwind/Petrov Galerkin (SUPG) stabilized virtual element method (VEM) for optimal control problem governed by a convection dominated diffusion equation is investigated. The virtual element discrete scheme is constructed based on the first-optimize-then-discretize strategy and SUPG stabilized virtual element approximation of the state equation and adjoint state equation. An a priori error estimate is derived for both the state, adjoint state, and the control. Numerical experiments are carried out to illustrate the theoretical findings.     

数值模拟水环境污染的一种L∞稳定的分步杂交方法

本文提出了一种数值求解对流扩散方程的分步杂交方法。在不规则的三角形网格上,采用迎风离散格式或改型特征线方法处理对流算子;采用集中质量的有限元方法处理扩散算子。详细分析了这种算法的稳定性同题,在数学上严格证明了在满足①Δt≤min((2d)/v,(d²)/(3K)),其中d是三角形网格中最短垂线的长度,V和K分别为流场中的最大速度和扩散系数。②所有三角形的内角θ≤π/2的条件下,整个计算格式是L_∞稳定的,从而保证了在海洋环境和水质的数值模拟中海水的盐度、污染物的浓度和核电站冷却水系统中的超温不会出现负值。应用非线性的对流扩散方程对此方法的精度和收敛性进行了检验。通过数值解与精确解的比较,表明本方法的数值耗散很小,用改型特征线方法处理对流算子较迎风离散格式有更高的精度;两种处理对流算子的方法都没有伪振荡现象发生。本方法由于具有算法简单、L_∞稳定、计算网格灵活等优点,可推广使用于实际的海洋环境(潮波、海流、海洋污染)、港口和海湾的数值模拟以及不可压粘性流和对流传热同题的数值计算。     

三维对流扩散方程非等距网格上的四阶紧致格式及其多重网格方法

提出了数值求解三维变系数对流扩散方程非等距网格上的四阶精度19点紧致差分格式,为了提高求解效率,采用多重网格方法求解高精度格式所形成的大型代数方程组。数值实验结果表明本文方法对于不同的网格雷诺数问题,在精确性、稳定性和减少计算工作量方面均明显优于7点中心差分格式。     

Simulations of complex fluids by mixed lattice Boltzmann—finite difference methods

We present the numerical results of simulations of complex fluids under shear flow. We employ a mixed approach which combines the lattice Boltzmann method for solving the Navier–Stokes equation and a finite difference scheme for the convection–diffusion equation. The evolution in time of shear banding phenomenon is studied. This is allowed by the presented numerical model which takes into account the evolution of local structures and their effect on fluid flow.     

A charge-conservative approach for simulating electrohydrodynamic two-phase flows using volume-of-fluid

In the present study we propose a charge-conservative scheme to solve two-phase electrohydrodynamic (EHD) problems using the volume-of-fluid (VOF) method. EHD problems are usually simplified by assuming that the fluids involved are purely dielectric (insulators) or purely conducting. Gases can be considered as perfect insulators but pure dielectric liquids do not exist in nature and insulating liquids have to be approximated using the “Taylor–Melcher leaky dielectric model” [1], [2] in which a leakage of charge through the liquid due to ohmic conduction is allowed. It is also a customary assumption to neglect the convection of charge against the ohmic conduction. The scheme proposed in this article can deal with any EHD problem since it does not rely on any of the above simplifications. An unrestricted EHD solver requires not only to incorporate electric forces in the Navier–Stokes equations, but also to consider the charge migration due to both conduction and convection in the electric charge conservation equation [3]. The conducting or insulating nature of the fluids arise on their own as a result of their electric and fluid mechanical properties. The EHD solver has been built as an extension to Gerris, a free software solver for the solution of incompressible fluid motion using an adaptive VOF method on octree meshes developed by Popinet [4], [5].     

DIMENSIONAL METHOD TO SOLVE THE DIFFUSION CONVECTION EQUATION OF SOLAR COSMIC RAYS

Physical processes of the propagation of the solar cosmic rays in the interplanetary space include the diffusion in interplanetary disordered magnetic fields and the convection in solar winds. Dimensional method can be applied to solve those equations convertible into Bessel equation, the results obtained are identical with those solved by the commonly used separate variable method. In order to derive an analytic solution to the diffusion convection equation in an unbounded, uniform medium, two dimensionless parameters reflecting the diffusion and convection characteristics of the particles are introduced. In the diffusion dominated case, the solution is similar in form to the diffusion of a source moving with the convection velocity and is modified by another convection term, which can be expanded into a power series of the convection parameter with coefficients composed of the generalized hypergeometric function series of the diffusion parameter. This solution has a clear physical meaning, and can suitably be used in the discussion of the rise phase characteristics of the solar cosmic rays from medium to high energies （E_p≥10¹ MeV）.     

Multigrid method based on the transformation-free HOC scheme on nonuniform grids for 2D convection diffusion problems

In this paper, a multigrid method based on the high order compact (HOC) difference scheme on nonuniform grids, which has been proposed by Kalita et al. [J.C. Kalita, A.K. Dass, D.C. Dalal, A transformation-free HOC scheme for steady convection–diffusion on non-uniform grids, Int. J. Numer. Methods Fluids 44 (2004) 33–53], is proposed to solve the two-dimensional (2D) convection diffusion equation. The HOC scheme is not involved in any grid transformation to map the nonuniform grids to uniform grids, consequently, the multigrid method is brand-new for solving the discrete system arising from the difference equation on nonuniform grids. The corresponding multigrid projection and interpolation operators are constructed by the area ratio. Some boundary layer and local singularity problems are used to demonstrate the superiority of the present method. Numerical results show that the multigrid method with the HOC scheme on nonuniform grids almost gets as equally efficient convergence rate as on uniform grids and the computed solution on nonuniform grids retains fourth order accuracy while on uniform grids just gets very poor solution for very steep boundary layer or high local singularity problems. The present method is also applied to solve the 2D incompressible Navier–Stokes equations using the stream function–vorticity formulation and the numerical solutions of the lid-driven cavity flow problem are obtained and compared with solutions available in the literature.