【摘 要】
：
本文在Newton迭代方法的基础上，对高阶精度间断Galerkin方法／有限体积(DG／FV)混合格式的时间隐式格式进行了研究。Newton迭代法的优势在于收敛效率高效，并且定常和非定常问题能够统一处理，对于非定常问题无需引入双时间步策略。为了避免大型矩阵的求逆，我们采用一步Gallss-Seidel迭代和Matrix-free技术消去残值Jacobi矩阵的上、下三角矩阵，从而只需计算和存储对角(
【机 构】
：
中国空气动力研究与发展中心空气动力学国家重点实验室,四川绵阳 621000 中国空气动力研究与发展
论文部分内容阅读
本文在Newton迭代方法的基础上，对高阶精度间断Galerkin方法／有限体积(DG／FV)混合格式的时间隐式格式进行了研究。Newton迭代法的优势在于收敛效率高效，并且定常和非定常问题能够统一处理，对于非定常问题无需引入双时间步策略。为了避免大型矩阵的求逆，我们采用一步Gallss-Seidel迭代和Matrix-free技术消去残值Jacobi矩阵的上、下三角矩阵，从而只需计算和存储对角(块)矩阵。对角(块)矩阵采用数值方法计算。利用该方法，数值模拟了Bump绕流和NACA0012翼型绕流。计算结果表明，与显式的Runge—Kutta时间格式相比，隐式格式所需的迭代步数和CPU时间均在很大程度上得到减少，计算效率能够提高1～2个量级。
其他文献
A selective, rapid and sensitive ultra performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) method was developed for the determination and pharmacokinetic study of azithromycin in h
本文理论研究井孔声电波场与地层参数的内在联系。首先，推导出低频斯通利波诱导的井孔轴向电场与声压比值的表达式。然后，通过数值计算和简化处理，获得了诱导电场和声场相位差异与地层参数的近似关系式。在此基础上，提出利用声电效应测井全波资料计算地层渗透率的方法。
本文基于结构对接网格／有限体积方法，采用可压缩大涡模拟方法求解N-S方程，对Mach=2.46，Re=1e6的圆柱超音速流场进行数值模拟，计算采用Schumann亚格子模型湍动能模式对小尺度涡作用进行模拟，数值模拟出了超音速底部流动的分离涡、膨胀波与再压缩激波等现象，从速度型以及压力分布看，采用的大涡模拟方法获取了较为可靠的数值模拟结果，为超音速底部流动的数值模拟提供了准确的数值模拟程序。
本文对将有源噪声控制技术应用于Kang等提出的通风隔声窗上进行了实验研究。实验结果表明，当噪声源为200至600 Hz白噪声时，按照本文所给出的实验配置，可以在误差传声器处实现全频带15 dB左右的降噪。在房间内的观察点处，在大部分频带范围内，有7至8 dB的降噪量。总声压级的降噪量在10 dB左右。这表明将有源噪声控制技术应用于此型通风隔声窗中，能够在低频取得较好的降噪效果。最后，分析了降噪不均
本文采用层流+经验转捩准则+湍流模式方法对高超声速边界层转捩进行预测。采用f(Reθ,Me)类经验转捩准则耦合SA、SST湍流模式，对超声速平板边界层、2°攻角尖锥裙流动和类X-51高超声速飞行器前体流动进行了数值模拟研究。给出了摩阻系数、热流系数和湍动能、涡粘性系数的计算结果，所采用的层流+经验转捩准则+湍流模式方法能够较为正确地预测高超声速边界层转捩。所采用的湍流模式本身的性能对边界层转捩的预
本文采用五阶wENO格式直接模拟了初始Taylor微尺度雷诺数为72～153，初始湍流马赫数为0.3～1.5的可压缩各向同性湍流，目的是研究可压缩均匀各向同性湍流的湍动能衰减规律。研究表明，可压缩均匀各向同性湍流湍动能存在一个幂指数形式的衰减，相应的衰减指数n=1.4～1.6，接近不可压缩均匀各向同性湍流统计平衡态前的弛豫过程，也就是湍动能初期衰减指数的理论预测值1.43；这表明压缩性对可压缩各向
本文对高阶格式涉及到的带宽优化和非线性格式技术进行了分析研究。带宽优化方面总结了具有代表性的优化方法，并提出类似Adams&Sharff的方法；应用该方法我们提出了四阶和五阶精度的线性带宽优化格式；研究指出WeirS&Manins的三阶格式存在不稳定逆耗散问题。非线性格式技术方面，我们分析了加权技术中的光滑度量构造，指出五阶WENO存在的一个精度分析错误，提出了两种基于变差的光滑度量；文中还发展了
雷诺数影响准确大型飞机气动特性的最重要的因素之一。针对高雷诺数流动特点，在WCNS高精度平台工程版本的基础上，建立了相应的数值模拟软件。软件中实现了基于壁面函数法的自动壁面处理，现有工作完成了自动壁面处理方法针对WCNS-E5的参数标定，并在湍流平板上对计算方法进行了初步的测试。
采用时间分裂法对化学源项与流动项进行解耦处理，建立了一种数值模拟化学非平衡流动的有限体积新算法，从线性模型方程出发可以证明具有时间和空间二阶精度。该算法有效解决超声速非平衡流动数值模拟时遇到的刚性问题，还具有公式推导简单易于软件工程化的优点。方法应用于超声速钝头体在可燃气体中诱导的振荡燃烧流动模拟，不但很好的分辨出流场内部的复杂结构，而且计算得到的振荡频率与文献一致，验证了方法的精度和适用性。在此
松耦合算法在数值模拟流固耦合问题时具有应用简单和计算高效的优势，但是需要解决流固界面上的信息传递算法。本文发现导致松耦合算法在界面上精度降低的原因是流体力学和固体力学方程性质和计算方法不同，理论上同时满足速度和位移相等条件是不可能的。速度相等条件体现了动量守恒，位移相等条件保证了能量守恒，数值实验表明位移相等条件得到的计算结果精度高。目前Farhat等人通过引入结构预测器概念构造出的二阶时间精度松