平板边界层内的流速分布实验

最新推荐文章于 2023-02-15 10:01:46 发布
最新推荐文章于 2023-02-15 10:01:46 发布
阅读量2.4k 收藏 5
点赞数 5
分类专栏： 杂货铺 文章标签： matlab 编程语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/dasdkhaksd/article/details/115013265
版权

目录

(1) 代码

clc;clear;

Rg = 287; %气体常数（J/(kg·K)）
g = 9.8; %重力加速度（m/s^2）
waterDensity = 1000; %水的密度（kg/m^3）

%对已知实验参数的输入
atmosphericPressure = input('请输入大气压：(Pa)');
temperature = input('请输入温度：(℃)');
KVC = input('请输入流体运动粘性系数：(10^-7m^2/s)'); %Kinematic viscosity coefficient
K = input('请输入风洞系数：');
gearFactor = input('请输入微压计所选的档位系数：');
inletVacuum = input('请输入风洞进口处真空度：(mmH2O)');
staticPressure0 = input('请输入流场静压：(mmC2H5OH)');

%第一截面实验数据输入
x10 = input('请输入第一截面距前缘的距离：(mm)');
y10 = input('请输入第一截面边界层内距离：(mm)');
micromanometerP10 = input('请输入第一截面用微压计测得的读数数据：(mmC2H5OH)');
transmitterP10 = input('请输入第一截面用压差变送器测得的读数数据：(Pa)'); %Differential pressure transmitter

%第二截面实验数据输入
x20 = input('请输入第二截面距前缘的距离：(mm)');
y20 = input('请输入第二截面边界层内距离的数据：(mm)');
micromanometerP20 = input('请输入第二截面用微压计测得的读数数据：(mmC2H5OH)');
transmitterP20 = input('请输入第二截面用压差变送器测得的读数数据：(Pa)');

%实验过程中通过计算出来的、不变的数据
airDensity = atmosphericPressure / (Rg * (273.15 + temperature)); %空气密度
staticPressure = atmosphericPressure - gearFactor * staticPressure0 * waterDensity * g * 0.001; %流场静压
x1 = x10 * 0.001;
y1 = y10 * 0.001;
x2 = x20 * 0.001;
y2 = y20 * 0.001;
velocityInf = sqrt((0.001 * (2 * K * g * waterDensity * inletVacuum)) / airDensity);
fprintf('=>风洞进口流速为：%.3f m/s\n',velocityInf);
micromanometerP1 = atmosphericPressure - micromanometerP10 * gearFactor * waterDensity * g * 0.001;
transmitterP1 = atmosphericPressure - transmitterP10;
micromanometerP2 = atmosphericPressure - micromanometerP20 * gearFactor * waterDensity * g * 0.001;
transmitterP2 = atmosphericPressure - transmitterP20;

%对第一截面数据进行处理
fprintf('>>>>>>>>>>**第一截面数据计算**<<<<<<<<<<\n');
Re1 = (velocityInf * x1) / (KVC * 1e-7);
fprintf('=>雷诺数 Re= %.3f\n',Re1);
fprintf('=>微压计实际读数计算出来的绝对压强：(Pa)\n');
disp(micromanometerP1);
micromanometerVelocity1 = zeros(1,length(micromanometerP1));
for i=1:length(micromanometerP1)
    micromanometerVelocity1(i) = sqrt(2 * (micromanometerP1(i) - staticPressure) / airDensity);
end
fprintf('=>用微压计测量的数据计算得到的边界层内的流速：（m/s）\n');
disp(micromanometerVelocity1);
transmitterVelocity1 = zeros(1,length(transmitterP1));
for i=1:length(transmitterP1)
    transmitterVelocity1(i) = sqrt(2 * (transmitterP1(i) - staticPressure) / airDensity);
end
fprintf('=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）\n');
disp(transmitterVelocity1);
micromanometerVelocityRatio1 = micromanometerVelocity1 / velocityInf;
fprintf('=>用微压计测量的数据计算得到的边界层内的速度比：\n');
disp(micromanometerVelocityRatio1);
transmitterVelocityRatio1 = transmitterVelocity1 / velocityInf;
fprintf('=>用压差变送器测量的数据计算得到的边界层内的速度比：\n');
disp(transmitterVelocityRatio1);

%对第二截面数据进行处理
fprintf('>>>>>>>>>>**第二截面数据计算**<<<<<<<<<<\n');
Re2 = (velocityInf * x2) / (KVC * 1e-7);
fprintf('=>雷诺数 Re= %.3f\n',Re2);
fprintf('=>微压计实际读数计算出来的绝对压强：(Pa)\n');
disp(micromanometerP2);
micromanometerVelocity2 = zeros(1,length(micromanometerP2));
for i = 1:length(micromanometerP2)
    micromanometerVelocity2(i) = sqrt(2 * (micromanometerP2(i) - staticPressure) / airDensity);
end
fprintf('=>用微压计测量的数据计算得到的边界层内的流速：（m/s）\n');
disp(micromanometerVelocity2);
transmitterVelocity2 = zeros(1,length(transmitterP2));
for i = 1:length(transmitterP2)
    transmitterVelocity2(i) = sqrt(2 * (transmitterP2(i) - staticPressure) / airDensity);
end
fprintf('=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）\n');
disp(transmitterVelocity2);
micromanometerVelocityRatio2 = micromanometerVelocity2 / velocityInf;
fprintf('=>用微压计测量的数据计算得到的边界层内的速度比：\n');
disp(micromanometerVelocityRatio2);
transmitterVelocityRatio2 = transmitterVelocity2 / velocityInf;
fprintf('=>用压差变送器测量的数据计算得到的边界层内的速度比：\n');
disp(transmitterVelocityRatio2);

%数据可视化
figure; %第一截面
hold on;grid on;
plot(micromanometerVelocityRatio1,y1,'p','markersize',10); %微压计测量所得数据点
plot(transmitterVelocityRatio1,y1,'*','markersize',10); %压差变送器测量所得数据点
micromanometerT1 = 0.5 : 0.01 : 0.95;
micromanometerU1 = spline(micromanometerVelocityRatio1,y1,micromanometerT1);
micromanometerBLT1 = spline(micromanometerVelocityRatio1,y1,0.93); %Boundary layer thickness
fprintf('=>用微压计测量得出的边界层厚度为：%.3f mm\n',micromanometerBLT1 * 1000);
plot(micromanometerT1,micromanometerU1,'-.','linewidth',2); %微压计测量拟合曲线。
transmitterT1 = 0.5 : 0.01 : 0.95;
transmitterU1 = spline(transmitterVelocityRatio1,y1,transmitterT1);
transmitterBLT1 = spline(transmitterVelocityRatio1,y1,0.93);
fprintf('=>用压差变送器测量得出的边界层厚度为：%.3f mm\n',transmitterBLT1 * 1000);
plot(transmitterT1,transmitterU1,'--','linewidth',2); %压差变送器测量拟合曲线

micromanometerY1 = 0 : 0.00001 : transmitterBLT1;
micromanometerX1 = ((3 * micromanometerY1 - (micromanometerY1.^3) / (transmitterBLT1^2))...
    / (2 * transmitterBLT1));
plot(micromanometerX1,micromanometerY1,'-'); %层流理论曲线

transmitterY1 = 0 : 0.00001 : transmitterBLT1;
transmitterX1 = (transmitterY1 / transmitterBLT1).^(1/7);
plot(transmitterX1,transmitterY1,'-'); %紊流理论曲线

title('第一截面边界层速度分布曲线');
xlabel('速度比');
ylabel('离平板壁面距离y（m）');
axis tight;
legend('微压计测量所得数据点','压差变送器测量所得数据点','微压计测量拟合曲线','压差变送器测量拟合曲线',...
    '层流理论曲线','紊流理论曲线','location','northwest');

figure; %第二截面
hold on;grid on;
plot(micromanometerVelocityRatio2,y2,'p','markersize',10); %微压计测量所得数据点
plot(transmitterVelocityRatio2,y2,'*','markersize',10); %压差变送器测量所得数据点
micromanometerT2 = 0.5 : 0.01 : 0.95;
micromanometerU2 = spline(micromanometerVelocityRatio2,y2,micromanometerT2);
micromanometerBLT2 = spline(micromanometerVelocityRatio2,y2,0.93);
fprintf('=>用微压计测量得出的边界层厚度为：%.3f mm\n',micromanometerBLT2 * 1000);
plot(micromanometerT2,micromanometerU2,'-.','linewidth',2); %微压计测量拟合曲线
transmitterT2 = 0.5 : 0.01 : 0.95;
transmitterU2 = spline(transmitterVelocityRatio2,y2,transmitterT2);
transmitterBLT2 = spline(transmitterVelocityRatio2,y2,0.93);
fprintf('=>用压差变送器测量得出的边界层厚度为：%.3f mm\n',transmitterBLT2 * 1000);
plot(transmitterT2,transmitterU2,'--','linewidth',2); %压差变送器测量拟合曲线

micromanometerY2 = 0:0.00001 : transmitterBLT2;
micromanometerX2 = ((3 * micromanometerY2 - (micromanometerY2.^3) / (transmitterBLT2^2)))...
    / (2*transmitterBLT2);
plot(micromanometerX2,micromanometerY2,'-'); %层流理论曲线

transmitterY2 = 0 : 0.00001 : transmitterBLT2;
transmitterX2 = (transmitterY2 / transmitterBLT2).^(1/7);
plot(transmitterX2,transmitterY2,'-'); %紊流理论曲线

title('第二截面边界层速度分布曲线');
xlabel('速度比');
ylabel('离平板壁面距离y（m）');
axis tight;
legend('微压计测量所得数据点','压差变送器测量所得数据点','微压计测量拟合曲线','压差变送器测量拟合曲线',...
    '层流理论曲线','紊流理论曲线','location','northwest');

%流态分析
fprintf('>>>>>>>>>>**流态分析**<<<<<<<<<<\n');
fprintf('一、对第一截面：\n');
fprintf('1.用雷诺数判定流态：\n');
if micromanometerBLT1 >= 3e6
    fprintf('=>（用微压计测量所得的数据分析）流动为紊流\n');
else
    fprintf('=>（用微压计测量所得的数据分析）流动为层流\n');
end
if transmitterBLT1 >= 3e6
    fprintf('=>（用压差变送器测量所得的数据分析）流动为紊流\n');
else
    fprintf('=>（用压差变送器测量所得的数据分析）流动为层流\n');
end
fprintf('2.用边界层厚度计算值与理论值比较判定流态：\n');
if abs(micromanometerBLT1 - 5 * sqrt((KVC * 1e-7 * x1) / velocityInf)) < 0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为层流\n');
elseif abs(micromanometerBLT1 - 0.37 * x1 * ((KVC * 1e-7) / (velocityInf * x1))^0.2)<0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为紊流\n');
else
    fprintf('=>测量值与理论值相差较大，无法判定流态！\n');
end
if abs(transmitterBLT1 - 5 * sqrt((KVC * 1e-7 * x1) / velocityInf)) < 0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为层流\n');
elseif abs(transmitterBLT1 - 0.37 * x1 * ((KVC * 1e-7) / (velocityInf*x1))^0.2) < 0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为紊流\n');
else
    fprintf('=>测量值与理论值相差较大，无法判定流态！\n');
end
fprintf('二、对第二截面：\n');
fprintf('1.用雷诺数判定流态：\n');
if micromanometerBLT2 >= 3e6
    fprintf('=>（用微压计测量所得的数据分析）流动为紊流\n');
else
    fprintf('=>（用微压计测量所得的数据分析）流动为层流\n');
end
if transmitterBLT2 >= 3e6
    fprintf('=>（用压差变送器测量所得的数据分析）流动为紊流\n');
else
    fprintf('=>（用压差变送器测量所得的数据分析）流动为层流\n');
end
fprintf('2.用边界层厚度计算值与理论值比较判定流态：\n');
if abs(micromanometerBLT2 - 5 * sqrt((KVC * 1e-7 * x2) / velocityInf)) < 0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为层流\n');
elseif abs(micromanometerBLT2 - 0.37 * x2 * ((KVC * 1e-7) / (velocityInf * x2))^0.2) < 0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为紊流\n');
else
    fprintf('=>测量值与理论值相差较大，无法判定流态！\n');
end
if abs(transmitterBLT2 - 5 * sqrt((KVC * 1e-7 * x2) / velocityInf)) < 0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为层流\n');
elseif abs(transmitterBLT2 - 0.37 * x2 * ((KVC * 1e-7) / (velocityInf * x2))^0.2) < 0.1
    fprintf('=>(用微压计测量所得的数据分析)流动为紊流\n');
else
    fprintf('=>测量值与理论值相差较大，无法判定流态！\n');
end

(2) 运行结果

输入：
请输入大气压：(Pa)98700
请输入温度：(℃)26
请输入流体运动粘性系数：(10^-7m^2/s)8.982
请输入风洞系数：1.08
请输入微压计所选的档位系数：0.2
请输入风洞进口处真空度：(mmH2O)22.8
请输入流场静压：(mmC2H5OH)118
请输入第一截面距前缘的距离：(mm)150
请输入第一截面边界层内距离：(mm)[0.45 1.45 2.45 3.45 4.45 5.45 6.45 6.95 7.45]
请输入第一截面用微压计测得的读数数据：(mmC2H5OH)[77.5 55 38.8 26 15 10.5 7.5 6.8 6.4]
请输入第一截面用压差变送器测得的读数数据：(Pa)[164.23 118.72 85.559 56.076 34.772 23.870 21.223 20.668 20.202]
请输入第二截面距前缘的距离：(mm)250
请输入第二截面边界层内距离的数据：(mm)[0.45 1.45 2.45 3.45 4.45 5.45 6.45 7.45 8.45 9.45]
请输入第二截面用微压计测得的读数数据：(mmC2H5OH)[83 64 51 39.8 29 21 14 11 9.8 8]
请输入第二截面用压差变送器测得的读数数据：(Pa)[173.85 129.31 104.47 79.882 55.702 40.633 28.892 22.750 19.701 18.167]
输出：
=>风洞进口流速为：20.490 m/s
>>>>>>>>>>**第一截面数据计算**<<<<<<<<<<
=>雷诺数 Re= 3421782.658
=>微压计实际读数计算出来的绝对压强：(Pa)
   1.0e+04 *
  列 1 至 5
    9.8548    9.8592    9.8624    9.8649    9.8671
  列 6 至 9
    9.8679    9.8685    9.8687    9.8687
=>用微压计测量的数据计算得到的边界层内的流速：（m/s）
  列 1 至 5
   11.7516   14.6568   16.4336   17.7118   18.7408
  列 6 至 9
   19.1458   19.4111   19.4725   19.5075
=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）
  列 1 至 5
   10.8004   13.9937   15.9222   17.4588   18.4898
  列 6 至 9
   18.9958   19.1166   19.1418   19.1630
=>用微压计测量的数据计算得到的边界层内的速度比：
  列 1 至 5
    0.5735    0.7153    0.8020    0.8644    0.9146
  列 6 至 9
    0.9344    0.9474    0.9504    0.9521
=>用压差变送器测量的数据计算得到的边界层内的速度比：
  列 1 至 5
    0.5271    0.6830    0.7771    0.8521    0.9024
  列 6 至 9
    0.9271    0.9330    0.9342    0.9353
>>>>>>>>>>**第二截面数据计算**<<<<<<<<<<
=>雷诺数 Re= 5702971.096
=>微压计实际读数计算出来的绝对压强：(Pa)
   1.0e+04 *
  列 1 至 5
    9.8537    9.8575    9.8600    9.8622    9.8643
  列 6 至 10
    9.8659    9.8673    9.8678    9.8681    9.8684
=>用微压计测量的数据计算得到的边界层内的流速：（m/s）
  列 1 至 5
   10.9246   13.5696   15.1150   16.3295   17.4207
  列 6 至 10
   18.1868   18.8316   19.1013   19.2081   19.3672
=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）
  列 1 至 5
    9.9957   13.3192   14.8532   16.2294   17.4774
  列 6 至 10
   18.2120   18.7644   19.0470   19.1857   19.2551
=>用微压计测量的数据计算得到的边界层内的速度比：
  列 1 至 5
    0.5332    0.6623    0.7377    0.7970    0.8502
  列 6 至 10
    0.8876    0.9191    0.9322    0.9375    0.9452
=>用压差变送器测量的数据计算得到的边界层内的速度比：
  列 1 至 5
    0.4878    0.6500    0.7249    0.7921    0.8530
  列 6 至 10
    0.8888    0.9158    0.9296    0.9364    0.9398
=>用微压计测量得出的边界层厚度为：5.211 mm
=>用压差变送器测量得出的边界层厚度为：5.736 mm
=>用微压计测量得出的边界层厚度为：7.133 mm
=>用压差变送器测量得出的边界层厚度为：7.491 mm
>>>>>>>>>>**流态分析**<<<<<<<<<<
一、对第一截面：
1.用雷诺数判定流态：
=>（用微压计测量所得的数据分析）流动为层流
=>（用压差变送器测量所得的数据分析）流动为层流
2.用边界层厚度计算值与理论值比较判定流态：
=>(用微压计测量所得的数据分析)流动为层流
=>(用微压计测量所得的数据分析)流动为层流
二、对第二截面：
1.用雷诺数判定流态：
=>（用微压计测量所得的数据分析）流动为层流
=>（用压差变送器测量所得的数据分析）流动为层流
2.用边界层厚度计算值与理论值比较判定流态：
=>(用微压计测量所得的数据分析)流动为层流
=>(用微压计测量所得的数据分析)流动为层流

在这里插入图片描述
在这里插入图片描述

参 考 资 料 来 源 参考资料来源 参考资料来源

  1. 《流体力学》.张鸣远.高等教育出版社.
  2. 《热与流体实验教程》.王小丹.孟倩.张可.吴青平.唐上朝.西安交通大学出版社.

源 码 作 者 ： A i d e n   L e e 源码作者：Aiden\ Lee 源码作者：Aiden Lee

博 客 创 作 ： A i d e n   L e e 博客创作：Aiden\ Lee 博客创作：Aiden Lee

特别声明：文章仅供学习参考，转载请注明出处，严禁盗用！

wibibaboo
关注
  • 5
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
两平行板间流动流体的速度分布：粘性流-matlab开发
05-29
完全由 Sreetam Bhaduri 创建。
边界层应用程序：这个 MATLAB 应用程序提供了一个 GUI 来研究平板上流动的层流边界层问题。-matlab开发
06-01
这是 Virtual Thermal/Fluid Lab 系列中的第一个 MATLAB 应用程序。 此 MATLAB 应用程序允许您： 1. 可视化边界层2. 研究边界层厚度随自由流速度的增长3. 可视化流线和速度剖面4. 学习如何用 TDMA 数值求解边界层问题5.查看GUI源代码，看看它是如何创建的 查看有关使用 MATLAB 进行虚拟流体力学和传热实验室的网络研讨会： http://www.mathworks.com/videos/teaching-fluid-mechanics-and-heat-transfer-with-interactive-matlab-apps-81962.html
平板上的Blasius方程和速度分布图：此代码解决了三阶非线性Blasius方程并在边界层内绘制了速度向量-matlab开发
05-31
此代码求解了三阶非线性 Blasius 方程并绘制了边界层内的速度向量。 可以进一步改进以在平板上绘制流线
边界层
05-13 2295
边界层，又称附面层是一个流体力学名词，表示流体中紧接着管壁或其他固定表面的部份。边界层是由黏滞力产生的效应，和雷诺数Re有关。 一般提到的边界层是指速度的边界层。在边界层外，流体的速度接近定值，不随位置而变化。在边界层内，在固定表面上流速为0，距固定表面越远，速度会趋近一定值。 边界层厚度：严格而言，边界层区和主流区之间无明显边界，通常以速度达到主流区的0.99U作为边界层的外缘。由边界层
【CFD学习】yplus，y+，边界层理论，壁面第一层网格
热门推荐
dabing1212的博客
07-29 7万+
谈谈本人对计算流体力学中y+的理解，欢迎批评指正！ 先思考几个问题： Q:什么时候需要考虑y+？ A:粘性计算（求解N-S方程）才涉及y+，无粘计算（求解Euler方程）不考虑边界层，不涉及y+的讨论 Q:为什么需要考虑y+?（理论基础） A：粘性计算时，壁面附近是流场中可以预先知道的高梯度区域，即边界层区域。即使边界层不是流场中速度梯度最大的区域，但至少是已知的大速度梯度区域。因此，应当预先考虑：需要用多少个网格，才能分辨出边界层呢？ 要回答这个问题，不仅需要预估边界层厚度，还需要大致..
matlab计算边界层厚度,边界层厚度.PPT
weixin_30134145的博客
03-25 1636
边界层厚度* 第三节 边界层的概念 由于N-S方程的非线性特征，使得问题的求解非常困难。在许多情况下，需要根据流动的特点对方程进行不同程度的简化。求解粘性不可压缩流体方程组的近似解时，可分为两种情况： 1、小雷诺数情形，此时粘性力较惯性力大得多(雷诺数代表惯性力与粘性力之比)，可以全部或部分地忽略惯性力，得到简化的线性方程。 2、大雷诺数情形，此时惯性力较粘性力大得多...
c语言Dvv实验报告,电渗实验报告
weixin_39921689的博客
05-19 1868
篇一：实验35电渗实验35 电渗一、目的①用电渗法测定SiO2对水的?电势②观察电渗现象，了解电渗法实验技术概要。二、基本原理电渗是胶体常见的电动现象的一种。早在1809年,就观察到在电场作用下，水能通过多孔沙土或粘土隔膜的现象(图Ⅱ.97)。这种现象是胶体常见的电动现象的一种。多孔固体在与液体接触的界面处因吸附离子或本身电离而带电荷，分散介质则带相反的电 荷。在外电场的作用下,介质将通过多孔固体...
传热学期末复习
zhuyujie0216的博客
12-20 9772
文章目录一、热量传递的三种基本方式1. 导热的特点2. 对流（热对流Convection)的概念2.1 对流传热的分类2.2 对流传热系数2.3 ==对流换热==的特点3. 热辐射的特点二、（总）传热过程三、温度场、稳态的概念四、导热基本定律五、热导率六、导热微分方程七、热扩散率 a八、定解条件九、热阻十、圆筒壁的导热热阻 一、热量传递的三种基本方式 热传导=导热、热对流、热辐射 【区别】 导热+热对流 =（表面）对流传热：流体经过固体壁 第五六七章 总传热过程：高温流体→固体壁→低温流体 第
边界层气象学期末复习笔记
m0_60394632的博客
02-15 4302
边界层气象学期末复习笔记，参考王咏薇老师的课堂ppt
绕圆柱湍流流场及其发展的计算* (2013年)
05-08
为研究绕圆柱湍流流场及其发展变化，应用RNG肛￡模型和有限体积法求解非定常雷诺平均纳维尔-斯托克斯(URANS)方程，以二阶隐式时间离散和二阶迎风空间离散格式对控制方程进行离散，采用SIMPLEC算法耦合求解平均速度场和压力场．在Re=103～10 7范围内对流场进行了数值模拟计算，得到了速度场、涡量场及其湍流脉动参数在下游的发展变化规律。计算与分析结果表明：在高雷诺数条件下，回流区减小，尾流区变窄，尾流在主流方向的影响区域变小，涡量向下游的衰减更慢，湍流动能和湍流耗散率在下游的影响区域也减小．
计算流体力学平板边界层流动换热数值计算源代码
09-15
在这个压缩包中，我们关注的是平板边界层流动换热的数值计算，这是一个在CFD领域常见的基础问题，对于理解和研究流体流动与传热有重要意义。 平板边界层流动通常发生在流体接近固体表面时，如空气或水流动经过平直...
波流边界层内外流速剖面的数值模拟
02-02
波流边界层内外流速剖面的数值模拟，张弛，李诚，准确模拟波流相互作用下的流场结构是河口海岸动力学研究与工程实践的一个重要课题。在传统波流边界层数学模型的基础上，引入了波
论文研究 - 矩形明渠横向剪切层及其流速分布
05-22
关于速度如何在壁边界层中横向变化还不是很了解。 本文根据Shiono＆Knight提出的深度平均动量方程，给出了矩形明渠中横向速度分布的解析解。 对于给定的粗糙通道，在壁面剪切层中获得的横向速度分布与横向涡流粘度...
混流泵内三维边界层计算.rar
05-05
边界层则是指流体在固体表面附近形成的一层薄流体区域，其流速从固壁的零速度逐渐过渡到主流区的速度。在混流泵中，边界层的行为直接影响泵的效率、流动损失和湍流特性。 计算混流泵内的三维边界层通常需要以下步骤...
边界层理论讲义.pdf
02-24
边界层理论讲义.pdf
工程流体力学笔记暂记34 (平板层流边界层的方程)
ResumeProject的博客
11-26 3046
粘性流体绕物体的流动中最简单的情况就是绕平板的运动 因为只有摩擦阻力 没有压差阻力 边界层动量积分方程的应用 根据理想流体区的欧拉积分方程 我们可以确定边界层外边界的压强为一常数 所以可得到边界层内所有点的压强均是相等的 将所得带入边界层动量积分方程 对其进行简化 所得方程对层流和湍流均是适用的 我们假设速度的方程为多项式形式 当然次数越高 方程的描述越精确、但其对应的计算就更为复杂 这里我们用二次的 a0 a1 a2 为待定系数 通过边界条件获得 第三个边界条件是由理想流体区的速度分布规律确定的
工程流体力学笔记暂记35 (平板层流边界层和平板混合边界层)
ResumeProject的博客
11-25 2845
湍流边界层比层流边界层复杂的多 目前还没有精确的解析解、 就基于边界层动量积分方程进行近似求解 同平版层流边界层的计算相类似 求解需补充速度和切应力的啷个关系式 下面利用圆管中的湍流流动的结果来给出平板湍流边界层速度和切应力关系式 普朗特认为演平板的边界层与管道流动的边界层没有明显区别 对于充分发展的管流可以认为边界层的厚度就是管道的半径r0 当流体绕过平板流动的时候，当层流边界层只占板的很小一部分时可以近似的看成整个版都是湍流 但是当层流边界层和湍流边界层相接近时就应该按混合边界层来计算
MATLAB相对声压分布彩图,热噪声试验中平板温度边界层对近壁面声压分布影响数值分析...
weixin_31600439的博客
03-22 301
声波入射到非均匀介质区域会发生散射，导致入射方向声压的衰减。设不均匀散射区的密度、压缩系数和声速分别为${\rho _{\rm{e}}}\left( r \right)$、${\kappa _{\rm{e}}}\left( r \right)$和${c_{\rm{e}}}\left( r \right)$，其中${c_{\rm{e}}}\left( r \right) = \dfrac{1}{{\...
matlab求解外掠平板层流边界层流动和换热
最新发布
11-15
matlab是一种强大的数学软件，可以用来求解外掠平板的层流边界层流动和换热问题。首先，我们可以通过编写数学模型来描述外掠平板的流动和换热特性，包括质量守恒、动量守恒和能量守恒方程。然后，我们可以利用matlab中的偏微分方程求解器来求解这些方程，得到外掠平板的流场和温度场的分布。 在matlab中，我们可以利用有限元法或有限差分法求解偏微分方程，得到外掠平板上的速度场和温度场的分布。通过求解得到的结果，我们可以分析外掠平板的边界层流动特性，包括速度剖面、边界层厚度和摩阻系数等。同时，我们还可以计算外掠平板的换热系数和温度分布，分析其换热性能。 除此之外，matlab还可以进行参数化研究和优化设计，通过改变外掠平板的几何形状和流动条件，分析其对边界层流动和换热性能的影响。通过matlab的求解和分析，我们可以更深入地理解外掠平板的流动和换热特性，为工程设计和优化提供重要参考。 总之，利用matlab求解外掠平板的层流边界层流动和换热问题，可以帮助我们深入理解其流动特性和换热性能，为工程设计和优化提供重要参考。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值