一种非结构结构多层混合网格方法及其应用

对于粘性绕流的数值模拟,在自适应直角网格基础上,结合三角形非结构网格和结构化网格,利用其各自的优势和特点,提出一种生成混合杂交网格的思路和方法.在物面附近生成适合粘性流计算的大长宽比结构化网格,在远场分布自适应直角网格,快速离散计算空间.对于复杂的多体问题,采用三角形网格来连接各体网格,并运用网格合并的方法,保证各网格之间的光滑过渡与连接,提高网格质量.针对一些二维、三维外形的绕流问题,在上述网格基础上,采用B-L代数湍流模型和中心有限体积法,完成Navier Stokes和Euler方程数值模拟的对比计算,结果表明网格生成和流场计算是正确的.

一种中心型有限体积孔隙–裂隙渗流求解方法及其OpenMP并行化

为高效求解单相孔隙–裂隙渗流问题,发展一种基于任意网格的三维中心型有限体积渗流求解算法,并对其进行OpenMP并行化。该算法将压力置于单元中心处;使用串联弹簧模型在空间域离散;使用显式差分格式在时间域离散;使用动态松弛求解技术,逐个单元求解。算例研究表明,该算法与有限元相比具有类似的精度,但求解效率更高。OpenMP并行化使得该算法运算速度在CPUi7–3770上可提高至4.0倍,在CPUi7–4770上可提高至4.2倍;两台机器上的并行效率均高达50%以上。

多块网格技术及其在操纵面绕流计算中的应用

采用中心有限体积法求解Navier Stokes方程组 ,利用点对点搭接及面搭接多块网格技术成功地进行了多段翼型和存在外形间断 (剪刀差 )的带操纵面三维机翼的网格生成及粘性绕流流场数值计算 ,使用重叠面积加权插值及Ramshaw算法实现了面搭接的多块网格之间的信息传递并保证了通量守恒 ,并将流场计算结果与实验数据进行了对比 .

翼面结冰对翼身组合体气动特性影响研究

飞机翼面结冰会影响气动性能,危及飞行安全。针对翼身组合体的翼面结冰问题,设计3种不同的冰型:钝头体、双角体和尖头体,对结冰前后的绕流流场及气动特性,采用数值方法进行分析研究。利用多块网格技术,结合法向外推代数方法和求解椭圆型偏微分方程方法,生成计算网格,采用中心有限体积法、LU-SGS双时间隐式推进算法,结合B-L湍流模型,运用流场分区求解技术,完成绕流流场的N-S方程数值模拟。计算与试验数据进行对比分析,说明了所采用方法的可行性和正确性,并给出了不同冰型对气动特性的影响结果。

裂隙岩体注浆扩散范围及注浆量数值模拟

为预测浆液在裂隙中的扩散范围,建立二维正交裂隙网络宾汉浆液渗流模型,采用中心型有限体积法进行数值模拟。分别研究不同参数下,浆液在裂隙中的扩散范围,给出多因素影响下的注浆扩散半径公式和注浆量公式。研究表明,注浆压力、浆液黏度、浆液剪切强度、裂隙开度和裂隙粗糙度对浆液扩散半径起到主导作用;注浆压力、钻孔长度和裂隙开度对注浆量起到主导作用。通过该公式可定量预测任意参量下浆液扩散范围和注浆量。