以喷管雷诺数，当地大气压及射流有效温度计算所需质量流量，静温

要解决的问题是：对于给定的Re及管道出口得到的Tr，应该设定质量流量为多少？

对于给定的Re

为了得到质量流量，需要确定喷管出口处μ的值，这由当地静温T0决定。

因此，需要确定T0的值。但目前我们测量得到的是Tr的值。

对于空气而言，

其中，T*为总温，T0为静温。热电阻测试得到的温度为有效温度，记为Tr，恢复系数记为r，则：

需要确定的参数就是r和Ma

恢复系数r一般与测温传感器的结构和安装方式（与气流垂直还是平行）有关。对于裸露的热电偶传感器，r的数据在0.6~0.9范围内。对于本试验所采用的热电阻而言，查阅相关资料获得其温度恢复系数与Ma的关系如下表所示。

拟合公式如下：

那么

需要确定马赫数

其中，密度又与当地静压和静温（absolute temperature）有关（注:甲方提供的PPT中将密度与总温静压联系在一起，应该是翻译错误）

P0为管道出口位置静压，一般与环境大气压一致。南京当地大气压为101166Pa，空气的R为287。

因此：

其中，Re, Tr,P0,D,γ,R为已知数，Ma,m,T0,μ为未知数。四个方程，四个未知数，理论上可联立求解。

首先写计算公式的头文件

#pragma once
#include <cmath>
extern double Pi;
extern double e;
extern double P0;
extern double k;
extern double D;
extern double Re;
extern double Tr;
double M_Re(double Re, double mu,double D)
{double M;M = (Re * Pi * D * mu) / 4;return M;
};
double Mu(double T0) //
{double Temp1 = (288.15 + 110.4) / (T0 + 110.4);double Temp2 = pow((T0 / 288.5), 1.5);double mu = Temp1 * Temp2 * 0.000017894;return mu;
};
double T0(double Ma, double Tr)
{double Temp = 1 +( 0.2 * (0.6495 * pow(e, (0.082 * Ma)) - 90.09 * pow(e, (-68.23 * Ma))) * pow(Ma, 2));double T0 = Tr / Temp;return T0;
};
double Ma(double M, double T0,double D)
{double Temp1 = 4 * M * sqrt(287 * T0);double Temp2 = Pi * P0 * D * D * sqrt(k);double Ma = Temp1 / Temp2;return Ma;
};

再写主程序

#include <iostream>
#include "Re.h"
#include <cmath>
using namespace std;
double Pi = 3.1415926;
double e = 2.718281828;
double P0 = 101300;//当地大气压
double k = 1.4;
double D = 0.03;/*0.01   0.02   0.03*///特征长度
double Re = 130000;//Re数
double Tr = 273.15 + 32.9;//测量有效温度
double M;//所需质量流量
double T;//静温
double V ;//射流速度
int i = 0;/*20mm_0.0115    0.020     0.029   0.0375*//*10mm_0.0055    0.010     0.014   0.018*//*30mm_0.0175    0.0305     0.0438   0.0565*/int main(){for (double m = 0.0565; m < 0.06; m = m + 0.0000001){for (double t = 273.15; t < 335.15; t = t + 0.01){V = Ma(m, t,D);T = T0(V, Tr);M = M_Re(Re, Mu(T), D);if (fabs(T - t) <= 0.01 && fabs(M - m) <= 0.0000001){cout << M<< "    " << Re << "    " << Tr << "    " << T << "    " << V << endl;};};};cout << "Calculate complete" << endl;return 0;
};

以喷管雷诺数，当地大气压及射流有效温度计算所需质量流量，静温相关推荐

喷管流动的守恒型CFD解法及激波捕捉（附完整代码）
入门CFD,主要参考书目<计算流体力学基础及其应用>(John D.Anderson 著,吴颂平等译) 实现了第 7.6 节激波捕捉的代码,采用的是 MacCormack 方法, ...
拟一维喷管流动的数值解——全亚声速等熵喷管流动的非守恒型CFD解法（MacCormack方法）
一.Matlab代码片 %全亚声速等熵喷管流动非守恒型麦考马克方法数值求解 clear; %清理内存变量 clc; %清理工作窗中的所有显示内容 r=1.4; %比热比 L=3; %喷管长度 i=3 ...
拟一维喷管流动的数值解——亚声速-超声速等熵喷管流动的非守恒型CFD解法（MacCormack方法）
一.Matlab代码片 %亚声速-超声速等熵喷管流动非守恒型麦考马克方法数值求解 clear; %清理内存变量 clc; %清理工作窗中的所有显示内容 r=1.4; %比热比 L=3; %喷管长度 ...
拟一维喷管流动的数值解——亚声速-超声速等熵喷管流动的守恒型CFD解法（MacCormack方法）
一.Matlab代码片 %亚声速-超声速等熵喷管流动守恒形CFD解法 MacCormack方法 clear; %清理内存变量 clc; %清理工作窗中的所有显示内容 r=1.4; %比热比 L=3; ...
Matlab 编程《计算流体力学基础及应用（约翰D安德森）》全亚声速等熵喷管流动CFD解法拉瓦尔喷管非守恒形式方程解法
Matlab 编程 <计算流体力学基础及应用(约翰D安德森)> 全亚声速等熵喷管流动CFD解法拉瓦尔喷管非守恒形式方程解法问题之全亚声速等熵喷管流动的CFD解法初始化参数迭代过 ...
流体力学概念总结(涉及所有重点)
表面力:又称面积力,是毗邻流体或其它物体,作用在隔离体表面上的直接施加的接触力.它的大小与作用面积成比例. 剪力.拉力.压力质量力:是指作用于隔离体内每一流体质点上的力,它的大小与质量成正比. ...
[渝粤教育] 上海交通大学流体力学参考资料
教育 -流体力学-章节资料考试资料-上海交通大学[] 单元作业单元测验 1.[单选题]一个物体微元的变形大小与所受切应力相关,这个物体可以是: A.刚体 B.弹性体 C.流体 D.无参考资料[ ] ...
【渝粤题库】陕西师范大学600007 化工基础作业（专升本）
陕西师范大学内部题库教育 (yuyueshool) 编制 <化工基础>作业一填空 1．化工产品种类繁多,生产工艺流程千差万别,具体内容是动量传递.质量传递.热量传递和 ,但其 ...
ROCKET PROPULSION ELEMENTS——DEFINITIONS AND FUNDAMENTALS笔记
火箭推进理论包括力学.热力学.化学单位换算磅力 lbflbflbf 磅,磅质量 lbmlbmlbm,lblblb,poundspoundspounds,lbslbslbs 1slug=1lbf1f ...

以喷管雷诺数，当地大气压及射流有效温度计算所需质量流量，静温

以喷管雷诺数，当地大气压及射流有效温度计算所需质量流量，静温相关推荐

最新文章

热门文章