非定常对流扩散方程的高精度多重网格方法

由已有的求解定常对流扩散方程的高阶紧致差分格式出发,直接推导出了数值求解非定常对流扩散方程的一种高阶隐式紧致差分格式,其时间为二阶精度,空间为四阶精度,并且是无条件稳定的。为了加快传统迭代法在求解隐格式时在每一个时间步上的迭代收敛速度,采用了多重网格加速技术。数值实验结果验证了本文方法的高阶精度、高效性及高稳定性。     

非定常NS/Boussinesq方程的高精度多重网格方法

提出了一种求解二维非定常不可压Navier—Stokes/Boussinesq方程的高精度全隐式紧致差分格式,为了提高隐格式的求解效率,在每一个时间步上,采用多重网格的全近似格式(FAS)加速其迭代收敛过程,其主要特点是既适于线性问题的求解又适用于非线性问题的求解。作为方法精确性和可靠性的验证,对方腔内部的自然对流问题进行了数值模拟。取Pr=0.71,在最大网格等分数为128×128网格上,Ra数最大算到107,所得结果与已有文献结果吻合的很好。     

二维对流反应方程的高精度多重网格方法

利用一阶偏导数项的四阶紧致差分算子,直接推导出了数值求解二维对流反应方程的一种新的高精度紧致差分格式。为了提高差分方程的求解效率,采用多重网格加速技术,建立了与之相适应的多重网格V循环算法。数值实验结果验证了本文方法的精确性和可靠性。     

构造定常对流扩散方程高精度紧致差分格式的新方法

以一维定常对流扩散方程的高精度差分格式为基础，提出了一种构造二维定常对扩散方程高精度紧致差分格式的新方法，并给出数值例子。     

二维半线性扩散反应方程的高精度全隐格式及其多重网格方法

针对二维非定常半线性扩散反应方程，空间导数项采用四阶紧致差分公式离散，时间导数项采用四阶向后Euler公式进行离散，提出一种无条件稳定的高精度五层全隐格式.格式截断误差为O（τ⁴+τ²h²+h⁴），即时间和空间均具有四阶精度.对于第一、二、三时间层采用Crank-Nicolson方法进行离散，并采用Richardson外推公式将启动层时间精度外推到四阶.建立适用于该格式的多重网格方法，加快在每个时间层上迭代求解代数方程组的收敛速度，提高计算效率.最后通过数值实验验证格式的精确性和稳定性以及多重网格方法的高效性.     

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

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

基于Voronoi网格的扩散方程差分格式

在Voronoi网格上利用一种基于回路积分法的有限体积法构造扩散方程的的差分格式.在这种特殊的网格上离散扩散方程比通常在四边形网格上离散的格式要简单,不会引进角点未知量,提高了对网格边上的流的离散精度,及差分格式整体精度.这种Voronoi网格上的扩散计算也可以与单元中心流体力学计算耦合.数值算例表明这种格式比四边形网格上的格式精度高,且能更好的应对网格扭曲情形.     

解二维对流扩散方程的块ADI方法

构造了解二维对流扩散方程的块交替方向隐式(ADI)方法。块ADI方法绝对稳定,并且具有计算和通讯的局部化特点,适合在分布存储的大规模并行计算系统上应用。文中给出了工作站机群系统上并行计算的数值例子     

一维非线性对流占优扩散方程的变网格特征差分方法

针对一维非线性对流占优扩散方程,提出了一类变网格特征差分格式,该格式能够根据解的梯度变化及时对计算网格进行调整.与均匀网格格式相比,给出的变网格特征差分格式对于对流占优扩散问题有着更好的计算效果.     

解二维时间相关光子扩散方程的交替方向隐式法

研究了应用交替方向隐式法求解非均匀生物组织中的二维时间相关光子扩散方程,并同解析解以及Monte Carlo模拟的结果作了比较,数值结果表明,该方法具有较好的稳定性和较高的精度.     

适应强扭曲网格的扩散方程时空二阶精度计算

以全局支撑算子方法为基础,通过引入面通量,构造了具有局部模板点的时空二阶精度格式。对于大变形扭曲网格,格式采用法向修正技术和合理的单元角体积计算方法,可以保持通量的精确性。算例表明该格式在非凸网格上能够精确获得线性解; 在非光滑网格上可以达到时空二阶精度; 能够较好地保持对称性; 并适合于三维非结构网格上的求解。     

基于时间的驱动方腔流的高精度多重网格方法数值模拟

提出了一种求解二维非定常不可压缩Navier-Stokes方程组的全隐二阶时间推进和空间四阶差分紧致格式,在每一个时间步上,采用多重网格的全近似格式(FAS)加速其迭代收敛过程。利用该方法对驱动方腔流动问题进行了直接数值模拟,结果显示对于Re≤5000,本文在粗网格上(64×64等分和128×128等分)即可得到较为准确的定常层流解;对于Re=7500和10000,由于更多二次涡的出现,本文在256×256等分网格上同样可得到与前人的结果相吻合的非定常周期性解。     

解三维扩散方程的积分内插法格式

研究三维扩散方程的数值模拟.在非正规六面体网格上,使用积分内插法建立扩散方程差分格式,涉及到27个相邻网格,适用于大变形网格上带间断系数的拟线性扩散方程的计算.叙述差分格式的建立,推导通量流和网格顶点温度的计算公式,给出了数值试验结果.     

基于强度传输方程的多重网格算法

提出一种新的切趾函数,用此切趾函数对过零单边干涉图加权,使具有相同光程差的两点光强与旋转因子的乘积之和为它们的平均值,减小了由于计算过程中零光程差点附近的数据被利用两次造成的误差.研究结果表明:与Mertz提出的切趾函数相比,本文提出的切趾函数对非对称性较严重的过零单边干涉图数据处理有更好的加权效果,能够有效减小光谱失真,同时计算效率得到了一定程度的提高,可广泛应用到傅里叶变换光谱仪中的过零单边干涉图处理中.     

基于强度传输方程的多重网格算法

为了从直接测量得到的强度图像中恢复相位信息,研究了基于强度传输方程的多重网格算法.从最粗层开始计算,给定初值,迭代出一个解,将此解作为最细层的初值,然后在最细层计算出一个近似解.进一步计算其残差,并将残差限制到较粗网格层求解,直至最粗层,然后逐层修正细网格层的解.利用循环在粗细不同的网格层来消除不同频率的误差分量,得到相位的精确解.仿真实验和真实实验表明,多重网格算法能够较好地从强度图像恢复物体的真实相位.     

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

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

基于亏量方程的多重网格解法

基于亏量方程提出了一种生成多重网格插值公式的新方法,新插值公式充分利用了粗网格的信息,因而具有更高的精度.对Poisson方程,各向异性方程,双调和方程,甚至三维问题的数值试验表明,新插值公式改进了多重网格法的渐近收敛速度,节省了存储空间及计算时间.     

随机对流扩散方程的数值仿真

本文针对对流一扩散随机过程在随机输入（即随机输运和源项）,作用下进行数值仿真。我们先将对流扩散随机微分方程中的随机函数采用有限项截断的多项式浑沌展开（Polynomial Chaos Expansion）展开,再由Galerkin映射法得到求解浑沌展开系数的确定性方程组。这是一个在物理空间包含多尺度解的大方程组。为此我...