需要一些指导,解决这个问题的速度和压力变量

1视图(30天)

显示旧的评论

Kumaresh Kumaresh 2022年8月29日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1789390-need-some-guidance-to-solve-this-problem-for-velocity-and-pressure-variables

评论道: Kumaresh Kumaresh2022年8月29日

你好所有的,

在我的问题,我需要解决以下方程计算变量速度(从情商。(1和2))和压力(从情商。(3))。

∇。 G =年代 (1)

∇× G - - - - - - G ×∇(ln (F1 + F2 * G )) = 0 (2)

∇ ^ 2 ( P ^ 2)= - (F1 + F2 * G )S - G ∙∇(F1 + F2 * G ) (3)

在那里, G , G = 质量速度矢量及其大小,P =压力和S =源项。

下面这是我的代码和查询:

1)定义连续性Eq。(1)时,我面临一个错误,

%错误:数组索引必须是正整数或逻辑值。

2)没有情商。(1),获得的速度和压力值并不像预期的那样好。

我不确定我的方法是如何找到合适的速度和压力。

有人请检查代码和判断我的方法,进一步帮助我。

谢谢你！

                          清晰;clc;关闭所有;
                         
                          L = 0.25;W = 5;
                         
                          Nx = 10;纽约= 200;
                         
                          x = linspace (0 L Nx);
                         
                          y = linspace (0 W纽约);
                         
                          (X, Y) = meshgrid (X, Y);
                         
                          dx = x / Nx;dy = y /纽约;
                         
                          rhoo = 0.8;
                         
                          为什么= 0.3;μ= 0.8 e-5;dia = 0.003;
                         
                          P = 101325。*的(纽约,Nx);
                         
                          V = 0(纽约、Nx);
                         
                          Nx Gvec = 0(纽约);
                         
                          步伐= 1;
                         
                          Vx =梯度(X,步伐);
                         
                          v =梯度(Y,步伐);
                         
                          Vmag =√Vx。^ 2 + v。^ 2);
                         
                          %边界条件
                         
                          V (: 1) = 0;
                         
                          :V(纽约)= 0;
                         
                          V (1:) = (Ny-1:);
                         
                          V (:, Nx) = (:, Nx-1);
                         
                          V(纽约,1)= 0;
                         
                          V (1,1) = 0;
                         
                          S = [0 0 0.1458 0.2058 0.2161 0.2171 0 0 0 0];
                         
                          f1 = (300 * ((1-por) ^ 2) * P *μ)/ (rhoo。* (dia。^ 2) * (^ 3));
                         
                          f2 = (3.5。* (1-por)。* P) / (rhoo * dia。* (^ 3));
                         
                          F1 = F1。*的(纽约,Nx);
                         
                          F2 = F2。*的(纽约,Nx);
                         
                          Gvec = rhoo。* V;
                         
                          Gmag = rhoo。* Vmag;
                         
                          %的差异Eq (Gvec) = 0%。(1)
                         
                          %错误:数组索引必须是正整数或逻辑值。
                         
                          d1 =日志(F1 + (F2。* Gmag));
                         
                          d2 =梯度(d1);
                         
                          A =旋度(Gvec X);
                         
                          B =旋度(d2, Gvec);
                         
                          V = A - B% Eq。(2)
                         
                             V =
                             200×10
                            
                              0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139 0.0139

                          d3 = F1 + (F2。* Gmag);
                         
                          d4 =梯度(d3);
                         
                          P =梯度(梯度((P . ^ 2))) +散度(d3, Gvec) + d3。*% Eq。(3)
                         
                             P =
                             200×10
                            
                              1.0 e + 08 * 0 0 0 0 0.0000 - 0.0000 0 0 1.5155 1.0178 1.4367 1.5086 1.0178 0.0000 0.0000 1.4367 1.5086 1.5155 0 0 0 0 0 0 0.0000 - 0.0000 0 0 1.5155 1.0178 1.4367 1.5086 1.0178 0.0000 0.0000 1.4367 1.5086 1.5155 0 0 0 0 0 0 0.0000 - 0.0000 0 0 1.5155 1.0178 1.4367 1.5086 1.0178 0.0000 0.0000 1.4367 1.5086 1.5155 0 0 0 0 0 0 0.0000 - 0.0000 0 0 1.5155 1.0178 1.4367 1.5086 1.0178 0.0000 0.0000 1.4367 1.5086 1.5155 0 0 0 0 0 0 0.0000 - 0.0000 0 0 1.5155 1.0178 1.4367 1.5086 1.0178 - 1.4367 1.5086 - 1.5155 0.0000 - 0.0000 0 0