duct flow理论进口速度分布

我们都知道对于泊肃叶流动，其解析解为速度沿y方向抛物线分布，如果拓展到三维，并且spanwize方向不再是无穷大，而是具有两个侧壁时，流动的速度分布由如下公式确定。

我们采用matlab来绘制出这个速度型，

%% write by rqli 2023/3/10
%%求解duct flow理论速度分布，仅有x方向，且公式为无量纲量。
format long
L = 8;W = 2;H = 1;
gridNumY = 20;gridNumZ = 40;
U = zeros(gridNumY,gridNumZ);
y = linspace(H/(2*gridNumY),H-H/(2*gridNumY),gridNumY);
z = linspace(W/(2*gridNumZ),W-W/(2*gridNumZ),gridNumZ);

for i = 1:gridNumY
    for j = 1:gridNumZ
        U(i,j) = duct(y(i),z(j));
    end
end

[Y,Z] = meshgrid(y,z);
surf(Y,Z,U')

% a = duct(0.1,0.1);
function ux = duct(y,z)
u1 = y-y^2;
u2 = 0;
for i = 1:2:50
    u2 = u2+(8*cosh(i*pi*(z-1))*sin(i*pi*y))/((i*pi)^3)/cosh(i*pi*1);
end
ux = u1-u2;
end```
对于有限体积法，数据存储在网格格心处，这里我们只给出了非边界处的速度，至于边界处的速度由无滑移边界条件，设置为0。这里并没有给出。
![duct flow速度分布](https://img-blog.csdnimg.cn/3f09b7c086704362ae89e18bcd109d90.png#pic_center)
over