沿全尺寸模型水平轴 MHK 涡轮机（DOE RM1）叶片性能表征计算（Matlab代码实现）

           目录

💥1 概述

📚2 运行结果

🎉3 参考文献

👨‍💻4 Matlab代码

---

💥1 概述

该存储库包括 MATLAB 脚本和输入文件示例，用于沿全尺寸模型水平轴 MHK 涡轮机叶片计算 3D AOA、CL 和 CD。该脚本是专门为DOE RM1涡轮机设计的性能表征而开发的，DOE RM1涡轮机设计是由美国能源部和美国国家实验室开发的两叶水平轴MHK涡轮机。MATLAB命令从涡轮机的 RANS 仿真中获取输入，并报告变量来表征叶片性能。

📚2 运行结果

 

 

🎉3 参考文献

[1]汪辉. 基于Magnus效应的共水平轴海流涡轮机的性能研究[D].华北电力大学,2018.

👨‍💻4 Matlab代码

主函数部分代码：

clc;
clear all;
close all;

%=====User Inputs==============
%Operating conditions
roh=1025;     %Density in [kg/m^3]
v_inf=1.9;    %Free stream velocity in [m/sec]
v_local=1.4;  %1D Upstream velocity in [m/sec]
omega=1.2;    %Turbine Anguler velocity in [rad/sec]

%Blade geometry specifications
pitch_angle=0;

file_geometry='DOE_ref_model1_geometry_28sections.txt';
fid_geometry=fopen(file_geometry,'r');

for j=1
    fgetl(fid_geometry);
end
C_geometry = textscan(fid_geometry, '%f %f %f %f %f');
fclose(fid_geometry);

radius=C_geometry{1};
twist=C_geometry{2};
chord=C_geometry{3};
pitch(1:size(radius),1)=pitch_angle;
thickness=C_geometry{5};

line_to_skip=39;
line_to_skip_center_of_pressure=4;
%=====End User Inputs==============

%===== Reading/Sorting Streamwise and Perpendicular Forces========%
%=====  on Each Blade Section Extracted from ANSYS FLUENT ========%

%=====Forces in the streamwise direction==============
file='fy_seawater_KW.txt';
fid=fopen(file,'r');

for j=1:line_to_skip
    fgetl(fid);
end
C_y = textscan(fid, '%s %f %f %f %f %f %f');
fclose(fid);
temp_y=C_y{4};     

%Data rearrangment from root to the tip of the blade
F_y(1:2,1)=temp_y(1:2,1);       %blade0 and 1
F_y(3,1)=temp_y(13,1);          %blade2
F_y(4:10,1)=temp_y(22:28,1);    %blade3  to 9    
F_y(11:20,1)=temp_y(3:12,1);    %blade10 to 19
F_y(21:28,1)=temp_y(14:21,1);   %blade20 to 27
clear temp_x temp_y;

%=====Forces perpendicular to streamwise direction==============
file='fx_seawater_KW.txt';
fid=fopen(file,'r');

for j=1:line_to_skip
    fgetl(fid);
end

转载：https://blog.csdn.net/m0_73907476/article/details/128528645