基于非结构自适应网格的复合有限体积法

利用文献[1]中将Lax-Wendroff格式和Lax-Friedrichs格式整体复合作用构成二维无结构网格上的复合型有限体积法,同时利用Delaunay方法,根据流场流动特性变化的梯度值为指示器对网格进行加密和粗化,实现自适应,并将此方法应用到二维浅水波方程的求解上,进行了二维部分溃坝,倾斜水跃的数值实验.结果表明,该方法是一个计算稳定、能适应复杂的求解域、能很好地捕捉激波、且计算速度快的算法.     

无结构网格有限体积方法及在热对流中的应用研究

首次将无结构三角网格的有限体积方法和压强连接半隐式算法相结合 ,用于求解非平行壁管道中的热对流问题 .并由此分析了化学汽相淀积薄膜生长的均匀性问题 .计算结果对于分析一类管道中热和动量输运现象均有普遍指导意义     

同位非结构和结构网格摄动有限体积法（PFV）求解二维Navier-Stokes方程组

根据NS方程组的一阶迎风和二阶中心有限体积（UFV和CFV）格式,导出NS方程组迎风和中心摄动有限体积(UPFV和CPFV)格式．该格式通过把控制体界面质量通量摄动展开成网格间距的幂级数,并由守恒方程本身求得幂级数系数而获得．迎风和中心摄动有限体积格式使用了与一阶迎风和二阶中心格式相同的基点数和相同的表达形式,宜于计算机编程．顶盖驱动方腔流和驻点流标量输运的数值实验证明,迎风PFV格式比一阶UFV、二阶CFV格式有更高的精度,更高的分辨率．尤其是良好的鲁棒特性．对顶盖驱动方腔流,在Re数从102到104范围内,亚松弛系数可在0.3～0.8任取,收敛性能良好．     

大涡模拟的代数模型及其应用

大涡模拟是近十多年来发展起来的湍流研究的重要手段。本文介绍了大涡模拟的基本思想,提出了大涡模拟的代数模型,并详细地讨论了模型的建立。并将它应用于平直槽道和弯曲槽道中湍流的数值模拟,计算的初步结果表明本模型是可行的,并有望在进一步改善后可用于较复杂湍流的数值模拟。     

基于非结构自适应网格的二维Euler方程数值求解方法研究

提出了一种基于非结构自适应网格的二维Euler方程的数值解法.采用有限体积法进行空间离散,通量计算采用Jamson中心格式,使得它适用于任意多边形计算单元.为了得到定常解,采用一种显式的四步Runge-Kutta迭代方法对时间进行积分.根据流场参数的变化梯度确定加密边,由加密准则进行自适应网格剖分,然后得到分布合理的加密过后的网格.求解二维Euler方程,对NACA0012翼型进行了数值模拟,通过对自适应前后的数值解的对比,说明所建立的方法是正确的.     

非结构网格上一类满足局部极值原理的三阶精度有限体积方法

对二维标量双曲型守恒律方程,发展了一类满足局部极值原理的非结构网格有限体积格式.其构造思想是,以单调数值通量为基础,通过应用基于最小二乘法的二次重构和极值限制器,使数值解满足局部极值原理.为保证数值解在光滑区域达到三阶精度,该格式可结合局部光滑探测器使用.本文从理论上分析了格式的稳定性条件,数值实验验证了格式的精度和对间断的分辨能力.     

基于线性化Navier-Stokes方程二维水槽内单涡到双涡的数值模拟

建立了Navier-Stokes方程的预估-校正有限差分方法,在此基础上求得了二维水槽内部单涡到双涡的数值解,所得结果与前人的数值结果和解析解吻合很好.数值模拟结果表明,自由振动运动中自由面波高因粘性作用会发生衰减,且Reynolds数越大衰减越缓慢.在短时间内倾斜加速度激励下对于不同Reynolds数会出现一定周期的单涡.经过长时间的倾斜激励,水槽内涡场由单涡变化成双涡,而且只在较低的Reynolds数条件下出现双涡.     

超音速探测器-刚性盘-缝-带型降落伞系统的大涡模拟研究

研究了Mach数为2时,流场不同块结构自适应网格加密精度对探测器-刚性盘-缝-带型降落伞系统的气动减速性能以及流场结构特性的影响.对于非定常可压缩流体流动,采用了兼顾激波与湍流的WENO(weighted essen-tially non-oscillatory)和TCD(tuned center difference...     

不可压缩流体三维Rayleigh-Taylor不稳定性的大涡模拟

对不可压缩流体三维Rayleigh-Taylor不稳定性问题建立被动标量输运模型,用大涡模拟方法计算了正弦初始扰动和随机初始扰动下不稳定性发展各个阶段的瞬时速度场和标量场,以及混合过程中计算尺度和亚格子尺度上的平均湍流脉动能、平均剪切应力和被动标量通量;分析了界面形状、被动标量浓度分布的演化规律及气泡、尖钉速度和混合层宽度随时间的变化规律,计算结果与其他数值模拟和实验结果相吻合,验证了大涡模拟方法应用于该问题的可行性.     

填充烧结铜球的T型管中流体混合的大涡模拟

在FLUENT软件平台上,运用大涡模拟湍流模型及Smagorinsky-Lilly亚格子尺度模型,对填充有烧结铜球多孔介质的T型管道内冷热流体混合过程的流动与传热情况进行了数值计算,与未填充多孔介质时混合区域内的平均温度和温度波动、平均速度和速度波动等数据进行了对比,并对温度波动进行了功率谱密度分析．数值结果表明,多孔介质可有效削弱T型通道流体混合区域内的温度和速度波动,有效降低1 Hz至10 Hz频域中的温度波动的功率谱密度．     

Poisson方程差分格式迭代解法中一些最优参数的有理拟合方法(英文)

本文针对二维Poisson方程五点和九点差分格式,导出了求解这些格式的SOR方法中最优松弛因子与区域剖分数的有理拟合公式,给出了Jacobi结合Chebyshev加速方法中Jacobi迭代矩阵谱半径的有理拟合公式.实际计算表明这些公式计算效果良好.     

变厚度簿板弯曲问题的任意网格差分解法

本文提出一种求解变厚度薄板弯曲问题的任意网格差分格式,可适应不同边界,各种荷载和复杂形状板.计算实例表明,该方法具有格式简单、通用性强、计算精度高,计算量少等特点.     

一类伪双曲型方程的特征-差分方法

给出了一类伪双曲型方程的特征－差分格式,得到位移ｕ和速度ｕ／ｔ的差分解和最优ｈ＾１模及ｌ＾２模误差估计,并对计算中遇到的离散点会落在区域外这一问题,给出了具体的解决方法。     

Sobolev方程基于POD的降阶外推差分算法

用奇异值分解和特征投影分解(proper orthogonal decomposition,简记POD)方法建立Sobolev方程的一种降阶外推有限差分算法,并给出误差估计.最后用数值例子,验证基于POD方法降阶外推有限差分算法的可行性和有效性.     

抛物方程基于POD方法的降阶外推差分算法

用奇值分解和特征投影分解(Proper Orthogonal Decomposition,简记POD)方法去建立抛物方程的一种降阶外推有限差分算法,并给出误差估计.最后用数值例子验证这种基于POD方法降阶外推有限差分算法的可行性和有效性.     

自然对流换热问题基于混合元法的差分格式及其数值模拟

研究自然对流换热问题,通过对于空间变量采用有限元离散而对于时间变量用差分离散,导出一种基于混合有限元法的最低阶的差分格式,这种格式可以同时求出流体的速度、温度和压力的数值解,并给出了模拟方腔流的自然换热的数值例子。