二维多介质欧拉自适应网格方法研究

要对聚能和流场的剧烈变化进行模拟，欧拉方法具有明显的优势。而在这些方面的研究中，所涉及的流场十分复杂，为达到所需的计算精度，必须采用很密的网格才能以较好的分辨率去模拟流场的剧烈变化部分和介质的界面，特别是大空间尺度流场局部细节的数值模拟，有些问题如果用统一网格计算，即使最快的计算机也不能提供足够的分辨率。所以目前计算机的内存和速度限制了整个计算区域的网格细分程度，对计算区域作局部的网格自适应细分、用大的动态空间分辨率划分流场是必要的和迫切的。     

欧拉流体力学数值方法中的界面计算

本文所提算法适用于二维和三维多介质流体力学两步欧拉数值方法中输运计算的混合网格(包括自由面网格)界面处理。在一个混合网格中，界面被近似地看作直线(二维)或平面(三维)。整个方法分为三步：(1)第一步，用混合网格周围的八个网格的介质面积份额(二维)或二十六个网格的介质体积份额(三维)确定界面的法线方向；第二步，用混合网格的本身的介质面积份额(二维)或体积份额(三维)确定界面的方程(位置)；第三步，用此直线方程求出通过网格边界的流以及下一时刻网格的面积份额(二维)或体积份额(三维)。最后给出了用此方法所做的一些数值计算及与SLIC算法的比较。     

三维空间二阶精度流体体积界面重构方法

给出三维空间网格模板含81个单元的最小二乘流体体积界面重构方法,并和Youngs方法及网格模板含125个单元的最小二乘流体体积界面重构方法进行比较.静态和动态的测试例子均表明:该方法能精确重构任意方向的平面界面,对C²光滑曲面它能达到二阶收敛精度.和网格模板含125个单元的最小二乘流体体积界面重构方法相比,在达到同样网格精度的条件下,减少了计算量,节省了计算时间,提高了计算效率.     

求解 Banach空间中的非线性方程的 修正的 Chebyshev迭代方法

本文给出了一个求解 Banach空间中的非线性方程的迭代方法 ,这一迭代方法实际上是对 Chebyshev迭代法的修正 ,它也是三阶收敛的 ,而且它对二次方程是四阶收敛的.     

二维高速碰撞问题欧拉数值模拟的混合网格计算

提出了适用于二维平面或轴对称多介质流体力学两步欧拉数值方法中输运计算的混合网格界面处理.在一个混合网格中,将界面近似看作直线.整个方法分为3步:①用混合网格周围的8个网格的介质面积份额确定界面的法线方向;②用混合网格的介质面积份额或体积份额确定界面的直线方程;③用此直线方程求出通过网格边界的流.给出了用此方法所做的测试、数值计算及与其它算法的比较.     

矩形网格上二阶精度流体体积界面追踪方法

欧拉方法模拟多流体问题的难点在于混合网格中界面的处理及物理量的确定。因此，界面处理技术是欧拉方法的主要研究内容之一。     

三维多介质欧拉有限体积格式的混合算法

对于多介质欧拉方法，混合网格物理量的计算是其难点和关键点之一。这里提出的方法是运用Yonugs界面重构技术确定出混合网格内物质的界面，界面确定后，混合网格内每一部分可能是非规则的四面体、五面体、六面体或七面体，采用对非规则区域适应性很强的有限体积法对每一部分分别进行计算。这种方法虽然比较复杂，但是它兼有拉氏方法的优点，因此计算出的混合网格内每一部分物质的物理量比较精确。     

三维多介质欧拉方法混合网格的有限体积计算格式

在采用Youngs界面重构技术的基础上，对三维欧拉方法混合网格的计算格式进行了研究。运用Youngs技术确定界面后，混合网格内每一部分物质一般不再是正规的六面体结构，可能是非规则的四面体、五面体、六面体或七面体。本文采用对非规则网格适应性很强的有限体积法对每一部分分别进行计算，给出了混合网格内每种物质的压力、速度、能量等的计算公式，比较有效地解决了混合网格的计算问题。