matlab大气湍流处理,rans湍流模型
一内容介绍传统雷诺平均模拟湍流模型大多基于布西奈斯克湍流涡粘本构模型，该本构模型是在研究半无限大平板流构型中得到的。
2部分代码%{
紊流-项目2:紊流的雷诺平均模拟数值模拟
通过一个渠道。
数值模型：
-长度混合(0等式)
——TKE(一个方程式)
-k-(2个方程)
问题描述：
从完全发展的飞机的直接数值模拟中提取的数据
湍流通道流。网址：托罗贾http://号。DMT。芬欧汇川。es/通道/数据/
Re _ tau=180 Re _ tau=550 Re _ tau=950 Re _ tau=2000
%}
清晰；色度控制中心帮助mainadd路径(gen路径(pwd))；
%数字保存路径
创建目录数字；addpath( .\数字\ );fpath=strcat(pwd， \ figures \ )；
savefigure=false
% - %
%问题配置
% - %
global Re lm n nfrec sigma _ k c CD c _ mu c _ ep1 c _ ep2 sigma _ EP h
Re=180% Re_tau
EPS=1e-3；%右侧(RHS)误差-静态稳定
nfrec=1000%绘图频率
h=1；%通道的一半（无量纲h==h^*=1)
n=圆形(Re/4)；网格节点的百分比
CFL=1；与粘性项相关的% Corant数
a=1.5节点的%指数放大器
mplot=真实百分比网格图？
x0=0；%向左
xf=h；%右界
itmax=2e6%最大迭代次数
% - %
% DNS数据
% - %
% y/h y U U v w -Om _ z Om _ x
% om _ y om _ z uv uw VW pr PS psto p
load(strcat(Re ，num2str(Re)， .prof’))；
第一节.d(:)=eval(Re ，num 2 str(Re)))；
% y/h y diss IP product p-应变 p-diff t-diff
% v-diff bal tp-kbal
load(strcat(Re ，num2str(Re)， bal。kbal’)；
第二节.d(:)=eval(Re ，num2str(Re)， _ bal )；
% - %
err=1；err 1=err err 2=err err 3=erit=1；
y(:1)=mesh1D(n，a，x0，xf，MP lot)；
% y(:1)=d(1).d(:1)；%复制网格域名服务器(Domain Name Server)
% n=长度(y)；%复制网格域名服务器(Domain Name Server)
yplus=y.*Re ./h；
deltay(:1)=y(2:end)-y(1:end-1)；
deltay max=max(deltay)；
mu=deltay(1:结束)。/deltay max；
% dt=CFL*min(deltay .^2*re)；
% - %
%初始化
% - %
函数初始化百分比
% u=15 .*y.^(0.25)；%初始化：混合TKE k-
% k=8 .*y.^(0.25).*预期(-2.5).* y)；%初始化：TKE k-
% ep=0.75 .*y.^(0.25).*exp(-5 .* y)；%初始化：k-
ep=500 .*y.^(0.25).*exp(-5 .* y)；%初始化：k-
% ep=750*y.^(0.25).*exp(-5 .* y)；%初始化：k-
精确初始化百分比- 1
y_dns=d(1).d(:1)；
u_dns=d(1).d(:3)；
k_dns=0.5 .*(d(1)).d(:4).^2 dd(:5).^2 dd(:6).^2);
ep_dns=-Re .*d(2).d(:3)；
u=spline(y_dns，u_dns，y)；
k=spline(y_dns，k_dns，y)；
% ep=spline(y_dns，ep_dns，y)；
精确初始化百分比- 2
% load(kepsilon_2000.mat )
% u=spline(y_kepsilon，u_kepsilon，y)；
% k=spline(y_kepsilon，k_kepsilon，y)；
% ep=spline(y_kepsilon，ep_kepsilon，y)；
% - %
来自域名服务器(Domain Name Server)数据的%湍流粘度：
dudy_dns=difx(u_dns，y _ DNS)；
tau12_dns=d(1).d(:11)；
nut_dns=-tau12_dns ./dudy _ DNS；
nut _ DNS=nut _ DNS(1:end-1)；%删除错误的值
% - %
nut _ spline=spline(y _ DNS(1:end-1)，nut_dns，y)；
dt=CFL .*(1/Re abs(nut_spline(1:end-1)).^(-1).*增量.^2;
dt=[dt；dt(end)]；
nyplus1=find(yplus 1，1， last )；
% dt=min(CFL .*(1/Re abs(nut_spline(1:end-1)).^(-1).*增量.^2);固定步长百分比
% - %
%显示配置
显示(-配置-)
fprintf(y=%.3e \n ，h)
fprintf(n=%d \n ，n)
fprintf(dy_min=%.3e \n ，min(deltay))
fprintf(dy_max=%.3e \n ，deltaymax)
fprintf(dt_min=%.3e \n ，min(dt))
fprintf(CFL=%.3e \n ，CFL)
fprintf(Re=%.3e \n ，Re)
fprintf(eps=%.3e \n ，eps)
fprintf(Nodes y 1=%d \n ，nyplus1)
显示(-)
如果nyplus1 2
注意了！y 1中的节点数小于2。)
目标
%%混合长度模型
混合；
%% TKE模型
动能
% %
% kepsilon% Chien
% kepsilonRM % Chien fmu Rodi Mansour
%开普西龙2%长野田川
kepsilon3% Launder-Sharmar
% kepsilon4% Lam-Bremhorst
% % k
% komega%还没有
%%加载绘图结果
npoint=1；每个n点的绘图百分比
负载_结果；
%%图混合，TKE和k-
普罗塔卢；
普罗托尔克
plotallep
plotallnut
%%以墙为单位绘制（点)
dnsplot1=20
dnsplot 2=8；
plotalluplus
plotallkplus
plotallepplus
plotallnutplus
%%以墙为单位绘制（线)
dnsplot1=16
dnsplot 2=8；
plotalluplus2
plotallkplus2
plotallepplus2
plotallnutplus2
%%绘制粘性应力与雷诺应力的关系图
普罗陶；
%%地块生产与耗散
plotalllep_prod_dis_plus
%%速度-亏损定律
绘图速度
%%仅绘制k-模型
% plotkep _ uplus
% plotkep _ kplus
% plotkep _ epplus
% plotkep _ nutplus
3运行结果
四参考文献[1]郭永涛，徐弘一。基于方、矩形环管直接数值模拟结果的各向异性雷诺平均模拟湍流模型构建[J]的缩写.复旦学报（自然科学版), 2019, 58(01):5-17.
[2]王翔宇106 . RANS/莱斯混合方法在湍流精细数值模拟中的应用与改进d .西北工业大学。
郑重声明：本文由网友发布，不代表盛行IT的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。