[发明专利]一种基于蒸汽流场计算液滴运动相变参数的方法有效
申请号: | 201910036741.8 | 申请日: | 2019-01-15 |
公开(公告)号: | CN109740281B | 公开(公告)日: | 2022-06-21 |
发明(设计)人: | 赵富龙;周娅;卢瑞博;谭思超;余霖;何宇豪;孟涛;曾陈 | 申请(专利权)人: | 哈尔滨工程大学 |
主分类号: | G06F30/28 | 分类号: | G06F30/28;G06F113/08 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 150001 黑龙江省哈尔滨市南岗区*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 蒸汽 计算 运动 相变 参数 方法 | ||
1.一种基于蒸汽流场计算液滴运动相变参数的方法,其特征在于,包括以下步骤:
步骤1:构建目标结构的几何模型,并进行网格划分,且针对局部位置区域进行网格加密;导出网格文件;
步骤2:将所述步骤1导出的网格文件导入流体计算软件,并计算蒸汽的流场参数信息,且网格节点的空间坐标与该网格节点对应的蒸汽流场参数信息之间呈一一对应关系;导出该蒸汽流场参数信息文件;
步骤3:将所述步骤2导出的蒸汽流场参数信息文件导入C++程序,建立KD搜索算法获取目标液滴位置周围最近的M个网格节点的空间坐标;
步骤4:获取所述最近的M个网格节点所对应的蒸汽流场参数信息,建立快速插值方法以计算目标液滴位置处的蒸汽流场参数信息;其中,M=1或4;
步骤4.1:当目标液滴的空间位置坐标点与最近的网格节点重合时,根据该重合的网格节点的空间坐标,获取对应的蒸汽流场参数信息,作为该目标液滴所在位置处的的蒸汽流场参数信息;
步骤4.2:若目标液滴的空间位置坐标点与最近的网格节点均不重合,则根据已经找出的离目标液滴位置xi距离最近的M个网格点,获取该M个网格节点的空间坐标所对应的蒸汽流场参数信息,并采用距离反比插值方法得到目标液滴位置xi处的蒸汽流场参数信息Fi,插值格式为:
式中:lj表示位置矢量xi和第j个网格点的欧式距离;
步骤5:利用步骤4获取的目标液滴位置处的蒸汽流场参数信息,求解液滴运动相变模型,以计算液滴运动相变参数;
步骤5.1:将目标液滴位置处的蒸汽流场参数信息在C++程序中传递给液滴运动相变模型;
步骤5.2:获取t时刻液滴的运动相变参数,以计算t+Δt时刻液滴的运动相变参数,计算使用的液滴运动相变模型包括以下方程中的至少一个:
其中,位于方程左侧的待求量包括:液滴的位置矢量x,液滴的旋转速度ω,液滴的速度v,液滴的半径r,液滴的温度T;
位于方程右侧的已知物理量包括:
(1)蒸汽流场参数信息,包括:蒸汽密度ρg,蒸汽流动速度u,蒸汽的温度Tg,蒸汽的压力Pg,蒸汽的动力粘度μg,流场旋度Ω,其中
(2)t时刻液滴的运动相变参数,包括:加速度a、速度v、位置坐标、温度Td、压力pd、密度ρd、粘度μd、定压比热容cp、汽化潜热γ、摩尔质量M;
(3)根据经验关系式计算得到的常数,包括:重力加速度g;气体常数R;蒸发冷凝系数α,对流换热系数h;
(4)根据蒸汽和液滴参数求得的相关系数,包括:转矩系数CM,曵力系数CD,升力系数CMa,升力系数CSa;
(5)归一化系数λ1~λ5,计算表达式分别为:λ1=-15ρg/16πρd、λ2=3ρf/(8ρdr+4ρgr)、λ3=3ρg(4ρd+2ρg)、λ4=[1.615(μd+2μg/3)2/(μd+μg)2(μfρg)0.5]/(ρdπr/3+ρgπr/6)、λ5=2(ρd-ρg)/(2ρd+ρg)。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学,未经哈尔滨工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910036741.8/1.html,转载请声明来源钻瓜专利网。