【技术实现步骤摘要】
包含任意流固界面非定常不可压流动数值模拟方法及设备
本公开涉及计算流体力学领域，尤其涉及一种适用于包含任意流固界面的非定常不可压流动的快速准确数值模拟方法及计算设备。
技术介绍
借助计算流体力学(computationalfluiddynamics，CFD)方法模拟围绕物体的非定常流动是许多工程研究设计工作中的重要环节。为提高计算模拟结果对实际设计应用的参考价值，所模拟的物体外形往往需要与真实情况大致吻合，因此这些模拟问题中的流固界面往往较为复杂。在计算机算力得到极大提升的今天，使用传统的计算流体力学方法仍需要首先构建生成贴合复杂流固界面的计算网格，而后采用更为繁琐的控制方程离散方法实现数值模拟。随着仿真性要求的提升，所处理的流固界面具有更高的几何复杂性，以上两个问题导致采用传统CFD方法的数值模拟方案所消耗的计算资源也会相应地提升。而且，根据不同问题中的流固界面，需要构建不同的网格，采用不同的控制方程离散方法，这导致数值模拟方案具有对流动模拟对象的特殊性，而不具备一般普适性。采用传统CFD方法来模拟复杂流固界面的非...
【技术保护点】
1.一种包含任意流固界面非定常不可压流动数值模拟方法，其特征在于，包括：/n使用笛卡尔网格识别流固界面，包括：/n使用笛卡尔网格离散描述模型问题的计算域，使用曲线坐标离散描述流固界面；以及/n描述界面-网格关系，明确界面附近网格点所处空间域的信息，标识为对应网格点的空间域旗标；以及/n基于空间域旗标，明确所有网格点相对流固界面的位置信息，标识为对应网格点的网格线旗标；以及/n施加流固界面边界条件，包括：/n根据网格线旗标，进行信息探测操作，调整流固界面附近网格点的空间离散模板，以在空间离散控制方程时植入当地流固界面边界条件；以及/n对空间离散处理后的流动控制方程进行时间推进。/n
【技术特征摘要】
1.一种包含任意流固界面非定常不可压流动数值模拟方法，其特征在于，包括：
使用笛卡尔网格识别流固界面，包括：
使用笛卡尔网格离散描述模型问题的计算域，使用曲线坐标离散描述流固界面；以及
描述界面-网格关系，明确界面附近网格点所处空间域的信息，标识为对应网格点的空间域旗标；以及
基于空间域旗标，明确所有网格点相对流固界面的位置信息，标识为对应网格点的网格线旗标；以及
施加流固界面边界条件，包括：
根据网格线旗标，进行信息探测操作，调整流固界面附近网格点的空间离散模板，以在空间离散控制方程时植入当地流固界面边界条件；以及
对空间离散处理后的流动控制方程进行时间推进。
2.如权利要求1所述的数值模拟方法，其特征在于，根据网格线旗标，进行信息探测操作，调整流固界面附近网格点的空间离散模板，包括：
当某个网格点处于流体域，存在包含在该网格点空间离散模板中的相关网格点不处于流体域的情况时，采用信息重构算法进行信息探测。
3.如权利要求2所述的数值模拟方法，其特征在于，所述采用信息重构算法进行信息探测包括：
根据位于流体域的该网格点的流场信息以及交点的流场信息，沿网格线方向进行线性外推插值得到该相关网格点的虚拟流场信息；
所述交点为该网格点与该相关网格点之间的网格线与所述流固界面的交点。
4.如权利要求3所述的数值模拟方法，其特征在于，所述线性外推插值采用的数学公式为：vG＝αBvB+αAvA+αEvE；
vG为正邻点或负邻点的虚拟流场信息；vB为网格线与所述流固界面交点的流场信息；vA为与正邻点或负邻点对应的该网格点的流场信息；vE为与该网格点相邻的另一个位于流体域的网格点的流场信息；αB,αA和αE均为插值系数；
当Δb≥Δg时，αE＝0；
当Δb<Δg时，αB＝2，
Δb为该网格点与所述交点之间沿网格线的距离，Δg为所述交点与正邻点或负邻点之间沿网格线的距离，Δe为该网格点与相邻流体域网格点之间沿网格线的距离。
5.如权利要求1所述的数值模拟方法，其特征在于，所述描述界面-网格关系，明确界面附近网格点所处空间域的信息，标识为对应网格点的空间域旗标，包括：
识别距离某一流固界面离散点最近的网格点并构建以该网格点为中心的网格点集合；
计算所述网格点集合中任意相邻两个网格点之间的网格线与流固界面的交点；
识别出与流固界面有交点的网格线对应的两个网格点；计算以该流固界面离散点为起点，分别以所述两个网格点为终点的两个单位向量；
将该流固界面离散点处的法向向量分别与所述两个单位向量进行点积，确认点...
【专利技术属性】
技术研发人员：潘天宇，李秋实，苏冠廷，郑孟宗，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。