基于多网格的 CFD 求解盖驱动腔问题附matlab代码

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

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

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

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

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

⛄ 内容介绍

实现了多种时间、空间数值格式以及多种湍流模式与嵌补网格技术的有机结合.开发了基于嵌补网格技术,适合复杂外型亚跨超音速流场模拟的计算软件.经多个典型算例的验证,证明软件及所使用的方法具有良好的精度和可靠性.对复杂流场具有很强的模拟能力.

⛄ 部分代码

%% Intermediate Velocity

function [u_tmp, v_tmp] = Predictor(nx, ny, u, v, Re, dx, dy)

    u_tmp = u; v_tmp = v;

    % x - momentum

    for i = 2 : nx+1

        for j = 2 : ny

                ue = (u(i,j) + u(i,j+1))/2; 

                uw = (u(i,j) + u(i,j-1))/2; 

                us = (u(i,j) + u(i-1,j))/2; vs = (v(i-1,j) + v(i-1,j+1))/2;

                un = (u(i,j) + u(i+1,j))/2; vn = (v(i,j) + v(i,j+1))/2;

                convective = -(ue*ue-uw*uw)/dx - (un*vn-us*vs)/dy;

                viscous = 1/Re*( (u(i+1,j) - 2*u(i,j) + u(i-1,j))/(dy*dy) + ...

                                 (u(i,j+1) - 2*u(i,j) + u(i,j-1))/(dx*dx) );

                u_tmp(i,j) = viscous + convective;

        end

    end

    % y - momentum

    for i = 2 : nx

        for j = 2 : ny+1

                ue = (u(i,j)   + u(i+1,j))/2;   ve = (v(i,j) + v(i,j+1))/2;

                uw = (u(i,j-1) + u(i+1,j-1))/2; vw = (v(i,j) + v(i,j-1))/2;

                                                vs = (v(i,j) + v(i-1,j))/2;

                                                vn = (v(i,j) + v(i+1,j))/2;

                convective = -(ue*ve - uw*vw)/dx - (vn*vn - vs*vs)/dy;

                viscous = 1/Re*( (v(i+1,j) - 2*v(i,j) + v(i-1,j))/(dy*dy) + ...

                                 (v(i,j+1) - 2*v(i,j) + v(i,j-1))/(dx*dx) );

                v_tmp(i,j) = viscous + convective;

        end

    end

end

⛄ 运行结果

⛄ 参考文献

[1] 宋慎义,王彦棡,刘冰,等.基于GPU的非结构网格CFD求解器的设计与优化[J].科研信息化技术与应用, 2012, 3(1):9.DOI:CNKI:SUN:KYXH.0.2012-01-005.

[2] ZHONG Ying,钟英,ZHANG Li-Lun,等.并行JFNK在结构多块网格CFD隐式求解中的应用[C]//全国高性能计算学术年会.中国计算机学会, 2013.

[3] 晏畅.基于MATLAB和CFD数据库的流场可视化的实现[J].河南科技, 2019(002):000.

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除