双时间隐式方法求解非定常绕流的相关问题

空间采用三阶精度的WeightedNND(WNND)格式,时间方向则采用具有二阶精度的双时间步隐式方法求解二维非定常Euler方程,数值模拟了NACA0012翼型俯仰强迫振动流场,讨论了内迭代步数、内迭代收敛判据以及物理时间步长对计算结果、计算效率与计算CPU时间的影响,并将计算结果与实验值进行了比较。     

测试投稿，大家不要管哈~

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

高超声速钝锥体俯仰阻尼导数影响因素分析

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

泰氟隆烧蚀产物对电子密度的影响

采用WNND格式,对有泰氟隆烧蚀产物引射的化学非平衡NS方程进行了数值模拟。采用7组元纯空气化学反应和19组元、28种反应的空气—泰氟隆化学反应系统,对照计算了壁面有泰氟隆烧蚀产物引射和纯空气绕流两种流场,研究了泰氟隆烧蚀产物对电子密度的影响。     

一类非线性类薛定谔方程的有限差分格式

给出了一类非线性类薛定谔方程的有限差分格式,推导了这类非线性类薛定谔方程的积分形式,并得到积分形式下系统的2个守恒定理.通过对积分方程离散化,并对系统的守恒量进行有限差分,得到了这类非线性类薛定谔方程的有限差分格式.     

五阶HWCNS格式在低速复杂流场中的应用

采用五阶精度显式HWCNS格式求解雷诺平均NS方程,利用多块对接结构网格技术,对30P-30N多段翼型进行网格收敛性研究,在不考虑转捩的情况下采用SA一方程湍流模型研究HWCNS格式与二阶精度MUSCL格式对该翼型压力分布和典型站位速度型的影响,并与实验结果进行对比分析。采用HWCNS格式和SA一方程湍流模型模拟梯形翼高升力构型低速复杂流场,通过对总体气动特性和压力分布的分析,探讨五阶精度显式HWCNS格式在低速复杂外形流动中的应用能力。     

椭圆型方程 Chebyshev 拟谱区域分解格式的等价形式

本文对一维、二维椭圆型方程建立了Ｃｈｅｂｙｓｈｅｖ拟谱区域分解格式,对这种拟谱区域分解格式给出了一种等价的广义变分形式。     

运送有毒物泄漏造成水体污染的扩散模型及求解方法研究

运送有毒物发生交通事故,导致有毒物泄漏引发水体污染,会造成严重影响。在分析污染事故的特点、危害的基础上,根据水力学及污染物的迁移转化规律,建立了有毒物泄漏水体后的转移弥散模型,即将控制水流的圣维南方程和描述有毒物迁移的对流弥散方程进行耦合。运用低数值耗散的半离散中心迎风格式求解方程,其结果可以预测污染河段长度和延续时间。最后,通过实例进行了模拟计算。     

非结构网格上解二维Hamilton-Jacobi方程的自适应算法

提出了一种非结构网格上解Hamilton Jacobi方程的自适应局部加密方法。讨论了加密方法关键的一步:光滑度指标的确定,通过采取无振荡插值,得到了无振荡数值格式,典型算例表明,该方法只需增加少量的点和计算量却能够提高计算精度和对间断的分辨能力。     

用动网格法计算理想平板的颤振导数

提出了一种数值模拟振动的理想平板绕流场 ,并提取其气动导数的方法。时间相关的不可压N S方程采用Projetion 2格式解耦 ,关于中间速度的动量方程的时间和空间离散采用二阶半隐格式 ,压力Possion方程迭代用多层网格法加速收敛。分别计算平板作竖向强迫振动或扭转强迫振动的气动力 ,用动网格法考虑平板和气流的耦合作用。由计算得到的气动力用最小二乘法确定 8个气动导数。计算结果和理想平板的Theodorsen理论值有较好的一致性     

振动矩形机翼非线性绕流数值研究

发展了可用于模拟包含运动边界的非定常流动的数值方法。该方法采用非结构动网格技术并构造了一种MUSCL类型有限体积格式,求解考虑了动网格效应的三维非定常Euler方程。非结构网格的变形运动由弹簧近似技术实现。用上述方法计算了作正弦振动的矩形机翼非线性绕流,计算结果与实验结果吻合,得到了流体运动与固体运动耦合的非线性流动特性。计算表明,该数值方法具有较好的精度和可靠性,完全可用于气动弹性或其他流固耦合问题的数值模拟。     

用高阶高精度 WENO 格式求解二维激波附面层干扰流场

采用时间相关法求解二维 Navier-Stokes方程 ,数值模拟二维平板层流附面层与激波干扰流场 ,给出了物面压力分布和应力分布。计算中 ,对流项空间导数的差分离散采用高阶高精度 WENO格式 ,时间方向采用具有 TVD性质的 Runge-Kutta方法 ,粘性项采用二阶中心差分。所得压力分布和应力分布与国外实验结果吻合较好 ,计算实践表明高阶 WENO格式具有优异的性能 ,应用前景广阔     

非定常超音速流动数值计算方法比较

为适用于非定常流场计算 ,在保证时、空二阶精度的前提下 ,通过引入变系数残值光顺技术 ,提高了四步Runge Kutta方法的计算效率 ;并针对球锥简谐振动过程 ,数值比较Runge Kutta多步法和“双时间法”的计算精度。     

球锥体高超声速绕流的俯仰阻尼导数的数值计算

采用Navier Stokes方程描述球锥体外形非定常振动流场 ,在Etkin理论下给出俯仰阻尼导数的计算公式。对定常流场的计算采用ADI形式的NND格式 ,对非定常流场的计算采用四步Runge Kutta方法 ,并引入变系数残值光顺技术加速收敛。将俯仰阻尼导数计算结果与实验及工程计算进行了比较 ,并数值研究了其随振动中心的变化规律。     

超声速燃烧流的双时间步计算方法研究

发展了模拟非定常超声速燃烧流场的隐式双时间步方法。采用时间分裂法对流动和反应进行解耦处理,流动方程组由双时间步方法求解,内迭代过程采用LU-SGS方法;反应源项方程组通过隐式二阶梯形公式求解。分析了时间分裂格式和时间步长对计算结果的影响,结果显示:一阶时间精度的分裂格式会略微高估化学效率,应该采用二阶时间精度的分裂格式;时间步长的选取对计算结果影响显著,为了保证解耦算法的计算精度,时间步长应足够小以能够较准确捕捉到主导各种输运过程的大尺度涡团的非定常行为。     

一类广义Liénard方程极限环的不存在性

利用Cauchy问题解的唯一性,研究了一类二阶非线性微分方程极限环的不存在性,得到了在奇点唯一的条件下极限环的不存在的充分条件,所得结果运用和发展了以往的方法。