基于四叉树网格下的浸入边界-格子Boltzmann方法

提出了一种基于非均匀四叉树网格下的浸入边界-格子Boltzmann方法。在不同层网格界面上,为了保证物理量的连续性,需要在时空方向上进行插值。由于四叉树同层网格在空间上步长相等、且相邻时间层上时间步长相等,所以在时空方向上采用平均值插值,这样做的优点是时空方向插值既可以达到二阶精度,又可以简化计算过程,节省资源,对任意边界加密下的网格,格子Boltzmann方法的实现比较容易。为了充分利用均匀笛卡尔网格的优势,物面边界的处理采用了速度修正法,与传统浸入边界-格子Boltzmann方法中的直接力法、动量交换法相比,无滑移边界条件得到了较好的保证。同时,把大涡模型加入到浸入边界-格子Boltzmann方法中,实现了在四叉树网格数据结构、边界处理技术、大涡模拟几种模型相结合下,绕障碍物的较大雷诺数流动的模拟。通过不可压缩粘性流中圆柱绕流算例验证,结果与其它方法结果吻合良好。

基于CV Newmark-β法的桥梁风致振动FSI数值模拟

为求解桥梁断面风致振动问题,首先介绍了两类数值微分方程解法,然后以Ansys Fluent为计算平台,通过嵌入自定义函数的方法实现了流线型桥梁断面的流固耦合数值模拟。通过理论推导发现,通过以常规的Newmark-β法嵌入UDF来驱动桥梁断面附近网格做刚体运动建立起的流固耦合计算模型计算得到的位移,与Fluent程序中网格更新的实际位移不一致,因而提出了一种修正速度的Newmark-β法以消除这种误差效应,并且建立了相应的桥梁断面流固耦合计算模型。针对某具体桥梁断面分别采用以常规Newmark-β法和修正速度的Newmark-β法建立的流固耦合计算模型,进行了低风速下和高风速下的桥梁断面风致振动数值模拟。研究表明:断面小振幅运动下不同计算模型获得的位移时程曲线基本吻合,以常规的Newmark-β法计算获得的位移与网格运动的真实位移之间的误差较小,低风速下这种算法与网格真实运动之间的位移不匹配效应产生的误差可以忽略;高风速桥梁断面大振幅颤振下不同计算模型获得的位移时程曲线差距较大,以常规的Newmark-β法建立起的流固耦合计算模型计算获得的位移与网格真实运动之间的位移不一致效应造成的误差不可忽略;低风速和高风速下以修正速度的Newmark-β法建立起的流固耦合计算模型获得的位移均与程序中网格真实运动位移一致。进行桥梁断面风致振动CFD数值模拟颤振问题时,应尤其注意处理这种位移误差效应,以建立合理的流固耦合计算模型。

带自由面流体运动的拉格朗日有限元分析

介绍了拉格朗日描述下的非定常不可压缩带自由面流体运动的有限元分析。在时间积分中采用速度修正法,采用速度修正分步法不仅可以让速度和压力均采用同阶插值函数,而且可以使算法变得更简单。用四边形单元对所求区域加以划分,并借助于Galerkin加权余量法导出相应的有限元方程组。通过实例计算,阐明了用该方法来分析带自由面流体运动的有效性及其实用性。

输油管线升力的数值模拟与实验研究

在波浪场中 ,利用 PIV装置进行水平输油管线升力及其周围流场的同步测量。同时采用速度压力修正法求解 N- S方程和泊松压力方程 ,得到了雷诺数较低条件下的水平输油管线周围的流场分布。