I can't find the error in this code. Please help with the dimension. Error in line 42

Question

Titas Chattopadhyay 2021-12-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1614685-i-can-t-find-the-error-in-this-code-please-help-with-the-dimension-error-in-line-42

评论： Titas Chattopadhyay 2021-12-21

clc
clear all
close all
%% Initial and Boundary conditions %%
g = 9.81;
J = 1.2./1000;
h = 0.075;
lambda = 0.075;
u_star = sqrt(g.*J.*h);
ny = 11; nz = 11; nt = 10; nit = 10;
%% Pressure term %%
P = (9.*u_star.^2)./h;
nu = 10^(-6)./(u_star.*h);
kappa = 0.4;  
XXx = 0:0.1:1;
YYx = (XXx.*(1-XXx))./(1+12.*0.19.*XXx.^2.*(1-XXx));
FF = trapz(XXx,YYx);
nu_t = kappa.*FF;
%% Grids %%
dy = lambda./(ny-1);
dz = h./(nz-1);
dt = 0.01;
y = 0:0.1:1;
z = 0:0.1:1;
f = zeros(nz,ny);
g = zeros(nz,ny);
fn = zeros(nz,ny);
gn = zeros(nz,ny);
[Z,Y] = meshgrid(z,y);
%% Main calculation %%
fn = f; gn = g ; jeta = 1:nz-1; eta = 1:ny-1;
for i = 2:nz-1
    for j = 2:ny-1
        f(i,j) = ((dt.*g.*J) + fn(i+1,j+1))./(1 + dt.*(2.*pi.*1i.*jeta).*fn(i,j) + (4.*dt.*pi.^2).*((nu+nu_t).*(jeta.^2 + eta.^2)+(nu_t.*jeta.*eta)) + (dt.*nu_t.*2.*pi.*1i.*jeta.*gn(i+1,j+1))) ;
        
        g(i,j) = ((dt.*g.*P) + gn(i+1,j+1))./(1 + dt.*(2.*pi.*1i.*eta).*gn(i,j) + (4.*dt.*pi.^2).*((nu+nu_t).*(jeta.^2 + eta.^2)+(nu_t.*jeta.*eta)) + (dt.*nu_t.*2.*pi.*1i.*eta.*fn(i+1,j+1))) ;
        
        
    end
end
Arrays have incompatible sizes for this operation.

2 个评论
显示无隐藏无

G A 2021-12-20

f(i,j) is a scalar, g = zeros(nz,ny) is an 11x11 array

Titas Chattopadhyay 2021-12-20

Thank you for your answer.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

VBBV 2021-12-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1614685-i-can-t-find-the-error-in-this-code-please-help-with-the-dimension-error-in-line-42#answer_858895

在 MATLAB Online 中打开

clc
clear all
close all
%%Initial and Boundary conditions %%
g = 9.81;
J = 1.2./1000;
h = 0.075;
lambda = 0.075;
u_star = sqrt(g.*J.*h);
ny = 11; nz = 11; nt = 10; nit = 10;
%%Pressure term %%
P = (9.*u_star.^2)./h;
nu = 10^(-6)./(u_star.*h);
kappa = 0.4;  
XXx = 0:0.1:1;
YYx = (XXx.*(1-XXx))./(1+12.*0.19.*XXx.^2.*(1-XXx));
FF = trapz(XXx,YYx);
nu_t = kappa.*FF
%%Grids %%
dy = lambda./(ny-1);
dz = h./(nz-1);
dt = 0.01;
y = 0:0.1:1;
z = 0:0.1:1;
f = zeros(nz,ny);
g = zeros(nz,ny);
fn = zeros(nz,ny);
gn = zeros(nz,ny);
[Z,Y] = meshgrid(z,y);
%%Main calculation %%
fn = f; gn = g ; jeta = 1:nz-1; eta = 1:ny-1;
for i = 2:nz-1
  for j = 2:ny-1
      f(i,j) = ((dt.*g(i,j).*J) + fn(i+1,j+1))./(1 + dt.*(2.*pi.*1i.*jeta(j)).*fn(i,j) + (4.*dt.*pi.^2).*((nu+nu_t).*(jeta(j).^2 + eta(j).^2)+(nu_t.*jeta(j).*eta(j))) + (dt.*nu_t.*2.*pi.*1i.*jeta(j).*gn(i+1,j+1))) ;
      g(i,j) = ((dt.*g(i,j).*P) + gn(i+1,j+1))./(1 + dt.*(2.*pi.*1i.*eta(j)).*gn(i,j) + (4.*dt.*pi.^2).*((nu+nu_t).*(jeta(j).^2 + eta(j).^2)+(nu_t.*jeta(j).*eta(j))) + (dt.*nu_t.*2.*pi.*1i.*eta(j).*fn(i+1,j+1))) ;    
  end
end