【摘 要】
：
最近，相关研究者提出了一种两步四阶时间离散方法(LW4)。 这种方法不同于传统的Runge-Kutta方法(RK)的地方在于，这种方法是Lax-Wendroff离散思想的衍生，是独立于Runge-Kutta方法之外的一种时间离散方法体系。但是目前该方法还没被用在有限差分领域。本文采用邓小刚等人提出的WCNS有限差分离散方案，实现了利用有限差分方法来执行四阶时间离散方法。相较于半离散匹配Runge-
【机 构】
：
北京应用物理与计算数学研究所,北京100094
论文部分内容阅读
最近，相关研究者提出了一种两步四阶时间离散方法(LW4)。 这种方法不同于传统的Runge-Kutta方法(RK)的地方在于，这种方法是Lax-Wendroff离散思想的衍生，是独立于Runge-Kutta方法之外的一种时间离散方法体系。但是目前该方法还没被用在有限差分领域。本文采用邓小刚等人提出的WCNS有限差分离散方案，实现了利用有限差分方法来执行四阶时间离散方法。相较于半离散匹配Runge-Kutta方法方案，新方法显示出了更符合守恒律方程本质的一些特点。
其他文献
针对超燃冲压发动机膨胀循环系统内低比转速小流量煤油离心泵，开展了数值计算与设计优化研究。利用Isight平台集成SolidWords、Ansys等软件实现数值计算自动化，并采用试验设计方法(DOE)建立样本数据库，进行多轮迭代优化。在此基础上进行了设计参数对离心泵操作特性的敏感性分析，获得了性能提升幅度较大的煤油离心泵设计方案。本文的研究为进一步开展关于煤油涡轮泵性能研究的数值计算与实验研究提供了
为了发展适用于高维精细化气动优化设计的更高效的代理优化算法，本文通过引入低可信度数据和梯度信息，发展了基于梯度增强分层kriging模型的变可信度代理优化方法。该方法的思想是：将NS方程疏网格计算的结果作为低可信度数据，将NS方程密网格计算的结果作为高可信度数据，同时采用Adjoint方法计算梯度，建立梯度增强的分层kriging模型，采用EI加点准则加入新的样本点，更新代理模型，直到优化收敛。本
针对有复杂状态方程的可压缩多相流动的高精度数值模拟，本文给出了一种基于特征重构和界面压缩的高精度数值算法。首先，在保证压力平衡的基础上，我们推导出了一种能保持五方程模型压力平衡的一种平均状态的表达式。借助于此平均状态，我们通过进行高精度局部特征重构来求解网格单元边界的左右状态。
对于含有移动界面管道瞬变流问题的计算，不论是特征线法还是其他网格方法，都不得不使用复杂的界面跟踪和数值插值技术。这不仅影响了计算的效率和稳定性，也可能因插值误差的积累而降低计算精度。针对以上问题，本文提出了一种分析含移动界面瞬变流的拉格朗日模型，给出了相应的无网格粒子计算方法，并分析了方法的数值精度。通过与典型工况下的实验结果对比，验证了本文所提出的拉格朗日粒子模型的有效性。
本文研究了一种基于自适应笛卡尔网格的相容守恒浸入边界法，主要解决两方面问题。一是固体颗粒在不可压缩粘性流体中的沉降或上浮问题。二是工程实际中具有几何构型复杂和多区域特点的传热问题。本文基于八叉树数据结构的各向同性自适应方法，采用两层次自适应网格策略。
本文提出了低雷诺数流场中弹性结构运动的浸没边界-格子波尔兹曼建模方法，其中，流场和弹性结构分别采用三维的格子波尔兹曼法和空间弹簧模型描述，而在流固运动界面上运用浸没边界法。然后运用该方法分析剪切流内纤维结构的动态运动问题，讨论了正弦周期剪切流条件下弹性纤维结构的运动规律，结果表明该方法可以较精确模拟复杂流动条件下流体-结构耦合系统的运动学和动力学特性，可以成为传统CFD建模技术的有力替代手段。
本文针对TVB限制器的改进进行了研究.TVB限制器广泛应用于求解各种间断问题，是一种抑制非物理性振荡有效的手段，但是TVB限制器中基于网格特征尺度的常数M的选取制约了TVB限制器的应用范围.本文利用TV随时间的变化率，将问题在时域上的划分为不同的物理阶段，让M按照规律分段选取不同的常数.更进一步设计了M的调节因子，使M能够跟随计算结果做出快速调整.为了研究和比较TVB限制器和改进的TVB限制器，应
有限体积法物理意义明确、适应性强，在工程问题中得到了广泛的应用，因此研究高阶精度有限体积法不但具有理论意义更具有应用价值。在传统有限体积法的基础上，构造适应任意质量网格的精确高阶精度有限体积法。本文方法基于网格单元及部分邻居单元的物理量均值，通过奇异值分解进行加权最小二乘求解，重构出网格单元内物理量分布的高阶逼近多项式，以此对Navier-Stokes方程的无粘项和粘性项进行高阶精度空间离散。以圆
尽管以二阶精度计算格式为基础的计算流体力学(CFD)软件已在复杂外形的复杂流动数值模拟中取得了巨大成功，但是二阶精度计算格式具有较大的数值耗散和色散，因此对湍流、分离等非常复杂的流动现象难以给出精细的流场结构。对于湍流的直接数值模拟(DNS)、大涡模拟(LES)和计算气动声学(CAA)，采用高阶精度格式是必要的。
Acker等人指出对于加权基本无振荡(WENO)格式，适当地提高相对光滑模板的权值比重，比改善极值点处的计算精度更有意义。本文首先分析了几种三阶WENO格式(包括WENO-JS、WENO-M、WENO-Z)的权值函数在相对光滑区的分布规律，进一步根据Acker等人的思想提出了三阶WENO格式新的映射函数。该函数形式简单，在保证格式基本无振荡的前提下，提高了非线性权值在相对光滑区趋于理想权值的比例，