【物理应用】基于傅里叶伽辽金谱法二维纳维-斯托克斯附matlab代码

简介: 【物理应用】基于傅里叶伽辽金谱法二维纳维-斯托克斯附matlab代码

✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。

🍎个人主页:Matlab科研工作室

🍊个人信条:格物致知。

更多Matlab仿真内容点击👇

智能优化算法  神经网络预测雷达通信 无线传感器

信号处理图像处理路径规划元胞自动机无人机 电力系统

⛄ 内容介绍

【物理应用】基于傅里叶伽辽金谱法二维纳维-斯托克斯附matlab代码

⛄ 部分代码

function examples

%% customVortices

N=128;

Re=200;

te=50;

w0=customVortices(N);

FGM2DNS=FGM2D_NavierStokes(N,Re,w0,0,te,[]);

FGM2DNS=FGM2DNS.Preprocessing_FGM2D;

FGM2DNS=FGM2DNS.Main_FGM2D(0,1);


%% taylorVortex

N=64;

locx=0;

locy=0;

R=0.3;

Re=1000;

te=10;

w0=taylorVortex(N,locx,locy,R);

w_exact=singleTaylorVortexSol(N,locx,locy,R,Re,te);


dt=0.01;


FGM2DNS=FGM2D_NavierStokes(N,Re,w0,0,te,dt);

FGM2DNS=FGM2DNS.Preprocessing_FGM2D;

close

FGM2DNS=FGM2DNS.Main_FGM2D(0,1);


plot(FGM2DNS.x(N/2+1,N/2+1:end),w0(N/2+1,N/2+1:end))

hold on

plot(FGM2DNS.x(N/2+1,N/2+1:end),w_exact(N/2+1,N/2+1:end),'*',FGM2DNS.x(N/2+1,N/2+1:end),FGM2DNS.w(N/2+1,N/2+1:end))

grid on

xlabel('$$x$$',Interpreter='latex')

ylabel('$$\omega$$',Interpreter='latex')

title("Taylor vortex")

legend('Initial condition','Exact solution','Numerical solution')


%% corotatingVortexPair

locy=[0 0];

Re=1000;

te=3;

N=128;

d=pi/6;

R=[0.2 0.2];

locx=[-d/2 d/2];

w0=taylorVortex(N,locx,locy,R);

FGM2DNS=FGM2D_NavierStokes(N,Re,w0,0,te,[]);

FGM2DNS=FGM2DNS.Preprocessing_FGM2D;

close

FGM2DNS=FGM2DNS.Main_FGM2D(0,1);

contourf(FGM2DNS.x,FGM2DNS.y,FGM2DNS.w,200,'edgecolor','none')

axis equal

colorbar

title("$$R=$$"+R(1)+"$$\quad d=\frac{\pi}{6}$$","Interpreter","latex")

xlabel("Dimensionless $$x$$","Interpreter","latex")

ylabel("Dimensionless $$y$$","Interpreter","latex")


%% kelvinHelmholtzInstability

N=256;

Re=5000;

w0 = twoEqualOppositeMixingLayer(N,19,pi/3);

w0(:,linspace(1,256,6))=w0(:,linspace(1,256,6))*1.1;

FGM2DNS=FGM2D_NavierStokes(N,Re,w0,0,5,[]);

FGM2DNS=FGM2DNS.Preprocessing_FGM2D;

FGM2DNS=FGM2DNS.Main_FGM2D(0,1);


%% periodicArrayOfVortices

locy=[pi/4 pi/4 pi/4 pi/4 pi/4 -pi/4 -pi/4 -pi/4 -pi/4 -pi/4];

locx=[-pi -2*pi/4 0 2*pi/4 pi -pi -2*pi/4 0 2*pi/4 pi];

R=[0.3 0.3 0.3 0.3 0.3 -0.3 -0.3 -0.3 -0.3 -0.3];

Re=100;

N=128;


w0=taylorVortex(N,locx,locy,R);


FGM2DNS=FGM2D_NavierStokes(N,Re,w0,0,30,[]);

FGM2DNS=FGM2DNS.Preprocessing_FGM2D;

FGM2DNS=FGM2DNS.Main_FGM2D(0,1);


return

⛄ 运行结果

⛄ 参考文献


⛄ Matlab代码关注

❤️部分理论引用网络文献,若有侵权联系博主删除
❤️ 关注我领取海量matlab电子书和数学建模资料


相关文章
|
7天前
|
存储 编解码 算法
【多光谱滤波器阵列设计的最优球体填充】使用MSFA设计方法进行各种重建算法时,图像质量可以提高至多2 dB,并在光谱相似性方面实现了显著提升(Matlab代码实现)
【多光谱滤波器阵列设计的最优球体填充】使用MSFA设计方法进行各种重建算法时,图像质量可以提高至多2 dB,并在光谱相似性方面实现了显著提升(Matlab代码实现)
|
7天前
|
机器学习/深度学习 传感器 算法
【高创新】基于优化的自适应差分导纳算法的改进最大功率点跟踪研究(Matlab代码实现)
【高创新】基于优化的自适应差分导纳算法的改进最大功率点跟踪研究(Matlab代码实现)
86 14
|
7天前
|
机器学习/深度学习 算法
【概率Copula分类器】实现d维阿基米德Copula相关的函数、HACs相关的函数研究(Matlab代码实现)
【概率Copula分类器】实现d维阿基米德Copula相关的函数、HACs相关的函数研究(Matlab代码实现)
|
7天前
|
机器学习/深度学习 传感器 算法
【裂纹检测】检测和标记图片中的裂缝(Matlab代码实现)
【裂纹检测】检测和标记图片中的裂缝(Matlab代码实现)
|
7天前
|
存储 算法 安全
【多目标工程应用】基于MOGWO的地铁隧道上方基坑工程优化设计研究(Matlab代码实现)
【多目标工程应用】基于MOGWO的地铁隧道上方基坑工程优化设计研究(Matlab代码实现)
|
7天前
|
传感器 机器学习/深度学习 编解码
【电缆】中压电缆局部放电的传输模型研究(Matlab代码实现)
【电缆】中压电缆局部放电的传输模型研究(Matlab代码实现)
|
7天前
|
算法 调度 决策智能
【复现】同时考虑考虑孤岛与重构的配电网故障恢复运行策略(Matlab代码实现)
【复现】同时考虑考虑孤岛与重构的配电网故障恢复运行策略(Matlab代码实现)
|
7天前
|
机器学习/深度学习 运维 算法
【微电网多目标优化调度】多目标学习者行为优化算法MOLPB求解微电网多目标优化调度研究(Matlab代码实现)
【微电网多目标优化调度】多目标学习者行为优化算法MOLPB求解微电网多目标优化调度研究(Matlab代码实现)
|
9天前
|
算法 计算机视觉
【MPDR & SMI】失配广义夹角随输入信噪比变化趋势、输出信干噪比随输入信噪比变化趋势研究(Matlab代码实现)
【MPDR & SMI】失配广义夹角随输入信噪比变化趋势、输出信干噪比随输入信噪比变化趋势研究(Matlab代码实现)
|
9天前
|
编解码 人工智能 算法
【采用BPSK或GMSK的Turbo码】MSK、GMSK调制二比特差分解调、turbo+BPSK、turbo+GMSK研究(Matlab代码实现)
【采用BPSK或GMSK的Turbo码】MSK、GMSK调制二比特差分解调、turbo+BPSK、turbo+GMSK研究(Matlab代码实现)

热门文章

最新文章