[发明专利]一种弯曲激波压缩流场气流参数的快速确定方法有效

专利信息
申请号: 201610070343.4 申请日: 2016-02-01
公开(公告)号: CN105574294B 公开(公告)日: 2018-11-23
发明(设计)人: 王磊;张堃元;苏纬仪;金志光 申请(专利权)人: 南京航空航天大学
主分类号: G06F17/50 分类号: G06F17/50
代理公司: 南京瑞弘专利商标事务所(普通合伙) 32249 代理人: 黄成萍
地址: 210016*** 国省代码: 江苏;32
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 弯曲 激波 压缩 气流 参数 快速 确定 方法
【权利要求书】:

1.一种弯曲激波压缩流场气流参数的快速确定方法,其特征在于:分别独立计算压缩壁面发出的压缩波系与弯曲激波反射的膨胀波系造成的流动参数的变化,在计算过程中,将原流场中的压缩波系简化为汇聚于点S的Prandtl-Meyer压缩波,将原流场中的膨胀波系简化为经过点S的一道膨胀波,根据简化后的Prandtl-Meyer压缩波和膨胀波快速确定流场中的压缩壁面参数、弯曲激波形状、波后参数、流场内流线形状、流线上参数以及出口参数;

设自由来流的马赫数、压力和流动方向角分别为M、p和θ,压缩壁面的前缘点与弯曲激波的前缘点重合,记为点Ostart(0,0),压缩壁面的末端点记为Wend,弯曲激波的末端点记为Send;压缩壁面的形状为yw=f(xw),0≤xw≤L;根据压缩壁面形状计算点Ostart位置的流动方向角θstart,然后根据气体动力学中的斜激波关系计算点Ostart位置的马赫数Mstart和压力pstart

对原流场中其他参数采用如下方法确定:

(1)在原流场中,经过压缩壁面上点Wi(xwi,ywi)的压缩波与弯曲激波的交点记为Sj;从原流场中分离出点Wi之前的压缩波系并简化为汇聚于点Sj的Prandtl-Meyer压缩波,根据点Wi与点Ostart的位置关系计算Wi位置的流动方向角θwi,然后采用Prandtl-Meyer公式计算点Wi位置的马赫数Mwi和压力pwi;同时计算自由来流经过转折角等于点Wi位置流动方向角θwi时产生的斜激波之后的马赫数压力和流动方向角

对压缩壁面上从点Ostart开始的每一点Wi重复步骤(1),直至点Wi到达点Wend,就确定了原流场中的压缩波造成的压缩壁面的参数变化;

(2)从原流场中分离出点Sj之前的膨胀波系并简化为经过点Sj的一道膨胀波,该膨胀波与压缩壁面的交点记为Wj;该膨胀波导致的转折角δj为:

式中:k为气体比热容;根据下式计算点Wj的坐标(xwj,ywj):

μwi=arcsin(1/Mwi)

对压缩壁面上从点Ostart开始的每一点Wj重复步骤(2),直至点Wj到达点Wend,就确定了原流场中的膨胀波造成的压缩壁面的参数变化;

(3)计算点Wj位置经压缩波和膨胀波合成后流动相对于点Ostart的转折角δwj-correction

δwj-correction=(θwjstart)-2δj

采用Prandtl-Meyer公式,根据Ostart位置的马赫数Mstart和压力pstart计算点Wj位置经压缩波和膨胀波合成后的马赫数Mwj-correction和压力pwj-correction

对压缩壁面上从点Ostart开始的每一点Wj重复步骤(3),直至点Wj到达点Wend,就确定了原流场中的压缩波和膨胀波合成后造成的压缩壁面的参数变化;

(4)计算点Sj位置相对于自由来流的转折角δsj-correction

δsj-correction=θwij

根据转折角δsj-correction和气体动力学理论计算点Sj位置的马赫数Msj-correction、压力psj-correction和流动方向角θsj-correction

根据下式计算点Sj的坐标(xsj,ysj):

μws=arcsin(1/Mws)

对压缩壁面上从点Ostart开始的每一点Wj重复步骤(4),直至点Wj到达点Wend,就确定了弯曲激波的形状和波后气动参数;

(5)在原流场中,经过点Wk的压缩波与以点Sj为起点的流线的交点为R,点R的流动方向角θr-correction以及点R相对于点Sj的转折角δr-correction为:

θr-correction=θwk

δr-correction=(θwk-2δj)-(θwi-2δi)

采用Prandtl-Meyer公式计算点R位置的马赫数Mr-correction和压力pr-correction;根据下式计算点R的坐标(xr,yr):

μwr=arcsin(1/Mwr)

对压缩壁面上从点Wi开始的每一点Wj重复步骤(5),直至点Wj到达点Wend,就确定了以点Sj为起点的流线的形状和流线上气动参数;

(6)出口截面的马赫数、压力和总压取为压缩壁面、弯曲激波和流线末端参数的平均值:

式中:n为总的流线数量,Mw-end、pw-end和为点Wend位置经压缩波和膨胀波合成后的马赫数、压力和总压,Mr-end、pr-end和为第r条流线末端位置的马赫数、压力和总压,Ms-end、ps-end和为点Send位置的马赫数、压力和总压。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于南京航空航天大学,未经南京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201610070343.4/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top