一种TVD方法在三维非结构网格中的应用

以三维非结构网格的显式有限体积法为基础,采用了一种TVD方法求解三维Euler方程,使用Roe通量来计算网格单元边界处的守恒量通量.为了验证方法的可行性,用该方法模拟三维爆炸问题,得出的结论是我们的方法可行,稳定且有效.     

基于粘性笛卡尔网格的N-S方程计算

本文采用基于叉树数据结构的粘性笛卡尔网格求解N-S方程。为了使网格的分布更为合理,采用基于物体几何外形的方法对网格进行自适应加密。在物面处采用光顺、投影、分割加密的方法保证了物面处网格的正交性与贴体性,从而生成高质量粘性笛卡尔网格。将有限体积法应用于粘性笛卡尔网格求解N-S方程,采用了SA紊流模型。算例计算结果与实验数据符合良好。     

二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

翼—身—尾组合体绕流的Euler方程计算

本文采用重迭网格技术和Ｅｕｌｅｒ方程计算翼－身－尾组合体绕流。对翼－身与尾－身部分采用各自的Ｈ－Ｏ型网格。Ｅｕｌｅｒ方程求解采用Ｊａｍｅｓｏｎ有有限体积法，即中心差分挖和显式Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身－尾模型为例，计算的空气动力特性与实验符合较好。     

三维带操纵面机翼欧拉方程数值模拟

本文给出了用欧拉方程求解带操纵面机翼的数值结果，采用有限体积法，四步Ｒｏｎｇｅ－Ｋｕｔｔａ时间推进，结合不连续面通量守恒传递技术，对有操纵面偏转的三维机翼流场进行求解，计算结果与实验符合较好。     

一类旋转体的侧面积与体积公式

证明了满足定理1条件的封闭曲线Ls绕X轴旋转一周的侧面积A=2kπ@ls(其中ls为Ls的长度).满足定理2条件的封闭曲线Ls绕X轴旋转一周的体积V=2kπAs(其中As为图形S的面积).     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

一种新的纤维增强复合材料细观力学模型

研究了复合材料宏、细观特征之间的联系，将宏观复合材料体中的一点赋予了细现结构特征。基于细现结构周期性假设，建立了一种细观力学模型。模型中用高阶多项式函数模拟基体和增强相中细现位移场，通过对细现单元力学方程的分析与求解，建立了复合材料宏、细观力学变量之间的联系。该细现力学模型，不仅能用于复合材料宏观有效性能的预测及细现应力、应变场的分析，而且容易融入常规有限元法中，实现对复合材料结构的宏、细观一体化分析。以该细观力学模型为基础的计算结果与试验结果及理论计算值具有较好的一致性。