一维Zakharov-Rubenchik方程的有效紧致差分格式

偏微分方程的有限差分法是科学计算中的一种有效方法,采用经典的一阶和二阶有限差分格式对方程进行数值求解,要想得到较高精度的近似解是不容易的,一种合理的方法是设计高阶紧致差分格式.为了研究一维Zakharov-Rubenchik方程的有效紧致差分格式及其数值计算.针对一般形式的Zakharov-Rubenchik方程,提出了一种半隐式紧致有限差分格式,该格式克服了传统差分格式效率低、精确度不足的缺点,并在离散层次上保持了质量和能量的守恒性.最后,通过数值算例验证了该格式的精确程度及守恒性,并对几种不同差分格式的误差和计算耗时进行了比较,数值结果表明了半隐式紧致差分格式的高阶收敛性及有效性.     

三维双调和方程的高阶紧致差分格式及其多重网格方法

提出一类求解三维双调和方程的高精度紧致差分格式.该类格式是以泊松方程的高精度格式为基础的四阶精度19点紧致差分格式和六阶精度27点紧致差分格式.采用多重网格方法求解由高精度紧致差分格式所形成的代数方程组,并与低精度方法进行比较.讨论多重网格方法中不同松驰算子的迭代收敛效果.数值实验结果验证四阶紧致差分格式和六阶紧致差分格式的精度以及多重网格方法的可靠性和高效性.     

Cattaneo模型的紧致有限差分法

紧致差分格式是一种高精度的有限差分方法.本文给出了Cattaneo模型的四阶紧致差分格式,通过对具体算例进行数值模拟,和二阶差分格式比较,验证了紧致差分方法的精确性和有效性.     

求解不可压流体Navier-Stokes方程的四阶精度有限容积紧致格式

采用四阶精度的有限容积紧致格式在交错网格上对二维非定常不可压流体的Navier-Stokes方程中的对流项和扩散项进行离散．压力项则由压力Poisson方程求得，并给出了新的压力Poisson方程的四阶精度有限容积紧致格式的离散表达式．用低存贮的三阶Runge-Kutta方法对Navier-Stokes方程进行时间推进．Fourier分析表明，有限容积紧致格式比一般的有限容积非紧致格式有更高的分辨率．最后以Taylor涡为例，得到了很好的结果．     

2维Maxwell方程的局部1维高阶紧致格式

将算子分裂方法与高阶紧致差分方法相结合,构造了2维Maxwell方程的局部1维紧致时域有限差分格式.该格式在时间方向和空间方向分别具有1阶和4阶收敛精度,并且具有计算效率高、无条件稳定的优点.数值实验表明:新构造的格式是能量守恒、高效率的.     

基于紧致交错差分格式的二维声波及黏滞声波方程数值模拟

针对常规网格差分难以适用于地震波场数值模拟中复杂介质的问题,首次将紧致交错有限差分格式应用于黏滞声波方程的数值模拟研究并同声波方程的数值模拟进行了模拟精度、频散关系和稳定性分析等方面的比较。理论研究结果表明：当差分精度相同时,紧致交错网格所需节点数要少于常规的中心差分和交错差分格式,计算效率更高;同常规的交错差分与中心差分格式相比,紧致差分的截断误差更小,数值频散也更低,能够适用于粗网格计算;差分精度相同情况下时进行数值模拟,紧致交错格式所需要的时间网格更小,稳定性条件也更为严格;紧致交错差分格式在完全匹配层(perfectly matched layer, PML)条件下,能够对边界反射进行有效吸收。最后,对均匀、水平层状介质以及Marmousi模型进行了黏滞声波方程的数值模拟和波场特征分析,实验结果证明了该方法对于复杂介质的数值模拟的适应性和有效性,并具有较高的模拟精度及计算效率。     

一种非均匀网格上的高精度紧致差分格式

提出了一种形式简单、网格剖分灵活、具有一定通用性的非均匀网格上的三点四阶紧致差分格式,对格式的截断误差进行了分析.采用文中提出的格式对Burgers方程和对流方程进行数值求解,并与均匀网格上的三点四阶紧致差分格式所得数值解对比,结果证明本文提出的格式对于大梯度问题的数值模拟有更高的精度.     

二维泊松方程非均匀网格上的高精度紧致差分格式

提出了数值求解二维泊松方程基于非均匀网格的高阶紧致差分格式,通过选取合适的网格分布参数求解具有边界层的数值算例,空间可以达到四阶精度.并与均匀网格上的计算结果进行比较,充分验证了本文非均匀网格高精度紧致格式的精确性和优越性.     

一类变系数非稳态对流扩散方程的四阶紧致差分格式

针对一类变系数非稳态对流扩散问题,构造了一种四阶Runge-Kutta高阶紧致有限差分格式.该格式具有时空四阶收敛精度,即O(h4,4τ),而且构造方法简单、易推广应用到其他问题.最后给出数值算例验证了所提出方法在求解非齐次对流扩散问题上的有效性和可靠性.     

一种求解RLW方程的紧致差分格式

利用紧致有限差分方法进行空间离散,龙格库塔方法进行时间离散,建立了一种求解RLW方程的数值格式,较好地解决了对空间与时间混合导数项的离散问题,并在空间和时间上都保持了高阶精度.所得数值结果证实了该数值格式具有较高的精度.     

Higher-order Accurate and High-resolution Implicit Upwind Finite Volume Scheme for Solving Euler/Reynolds-averaged Navier-Stokes Equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｅｍｅｒｇｅｎｃｅｏｆＣＦＤａｂｏｕｔ15ｙｅａｒｓａｇｏｐｒｏｖｉｄｅｄａｍａｊｏｒｉｍｐｅｔｕｓｔｏｓｏｌｖｅｔｈｅＥｕｌｅｒａｎｄＮａｖｉｅｒＳｔｏｋｅｓｅｑｕａｔｉｏｎｓｇｏｖｅｒｎｉｎｇｔｈｅｆｌｏｗｆｉｅｌｄｉｎｅｘｔｅｒｎａｌａｎｄｉｎｔｅｒｎａｌｆｌｏｗｓ.ＭａｊｏｒｐｒｏｇｒｅｓｓｈａｓｂｅｅｎｍａｄｅｉｎｔｈｅｆｉｅｌｄｏｆＣＦＤａｌｇｏｒｉｔｈｍｄｅｖｅｌｏｐｍｅｎｔ,ｓｕｃｈａｓｍｕｌｔｉｇｒｉｄａｃｃｅｌｅｒａｔｉｏｎｔｅｃｈｎｉｑｕｅｓ,ｄｉ…     

四阶紧致格式有限体积法湍流大涡模拟

为准确预测不可压复杂湍流,提出了一种应用于大涡模拟的高精度有限体积法。该方法在非交错网格上数值求解大涡模拟方程,空间离散采用有限体四阶紧致格式,时间推进用四阶Runge-Kutta法,压力速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动和振荡平板上方的流动,证实了该方法具有近四阶精度;并在此基础上,采用动力Smagorinsky亚格子应力模式,成功地完成了充分发展槽道湍流的大涡模拟计算,所得结果与直接数值模拟结果吻合良好。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

一维对流扩散方程第三边值问题的紧有限体积格式

针对一维常系数对流扩散方程第三边值问题提出一种紧有限体积格式,该格式形成的线性代数方程组具有三对角性质,可以使用追赶法求解.用能量估计法证明了格式按照离散L2范数、H1半范数和最大模范数均具有4阶收敛精度.数值算例验证了理论分析的正确性,并说明了格式的有效性.     

三对角四阶紧致差分格式的优化和初步应用

 提高数值解的精度和分辨率,有助于更精确地求解日趋复杂的工程问题。本文依据差分格式的伪波数应该在尽可能大的波数范围内接近物理波数的思想,构造了满足四阶精度的具有高分辨率的三对角紧致差分格式。一方面,它可以与近些年发展的求解(循环)三对角方程组的高效算法相结合,以更高的分辨率、更小的计算量来计算一阶导数;另一方面,与传统格式相比,该格式的最大精确求解波数可以达到2.5761,大于传统格式的1.13097。因此,优化格式更适合模拟小尺度波动。数值计算结果表明：(1) 虽然优化格式仍然是四阶精度,但要比传统四阶紧致差分格式的计算误差小,尤其对于小尺度波动,优化格式的计算误差会更小;(2) 对于行波问题,优化格式能够更加准确地模拟波动的传播行为,其优势也更加明显。理论分析和数值算例的比较结果均表明,优化的紧致差分格式更适合求解小尺度波动问题。     

三维不可压缩N-S方程的紧致有限差分和Fourier谱方法

采用高精度紧致有限差分－Fourie谱杂交的方法直接数值模拟了三维不可压缩的NavierStokes方程,该算法的时间离散采三阶精度混合显隐分裂格式,空间离散则结合Fourie谱方法及高精度紧天才有限差分逼近,该方法与普通的有限差分格式相比,具有很高的逼的精度及波数分辨率,针对三维平面槽道流的情况,应用该算法,直接数值模拟了三维T－S波在平面槽道流的传播问题,计算结果与流动稳定性分析结果吻合一致。     

二阶抛物方程的高精度差分格式

利用广义差分法构造二阶抛物方程的广义差分格式,讨论它的相容性、稳定性和收敛性．对经典差分格式、紧致差分格式、广义差分格式进行比较,并用数值例子验证理论结果     

高阶紧致格式结合外推技巧求解对流扩散方程

将Padé 逼近方法、Richardson外推技巧、两重网格算法与紧致差分方法相结合,对线性对流扩散问题构造了一种新的高阶紧致差分格式,给出了格式的实现过程和误差估计,分析了稳定性。最后给出数值算例说明算法的有效性。     

基于非稳态对流扩散方程高阶紧致差分的煤自燃数值模拟

为提高求解煤自燃非稳态对流扩散方程的有限差分法精度,构造高阶紧致差分格式,采用Richardson外推法进行求解;并针对煤自燃多场耦合模型的求解复杂度,提出一种具有承袭性的解耦策略。利用本文的解法和策略,对宁夏枣泉矿煤进行煤自燃多场耦合数值模拟,模拟计算的结果与相同实验条件下的实验结果相符。     

紧致差分格式的分辨率与精度的实例比较与讨论

通过比较先前建立的4阶最优紧致差分格式,以及传统的6阶和8阶紧致差分格式,来研究精度和分辨率之间的关系,主要比较了空间离散格式的有效波数范围、实际数值计算精度、以及对小尺度波动的模拟能力.数值试验结果表明:(1)这3种格式的计算精度都可以达到理论精度,并且此时精度越高,误差越小;(2)对于小尺度波动,最优4阶紧致格式比6阶和8阶紧致格式具有更高的分辨率;(3)对于行波问题,最优4阶紧致格式能够更加准确地模拟波动的传播行为.理论分析和数值算例的比较结果均表明,数值格式的精度和分辨率并不能互相替代,而是要根据计算问题的需要选择具有合适的精度和分辨率的数值格式.