SIMPLE-Algorithmus Berechnung der Geschwindigkeitskorrektur

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

PaulM

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 18.09.2014, 14:32 Titel: SIMPLE-Algorithmus Berechnung der Geschwindigkeitskorrektur

Hallo liebe Forum- User,

ich schreibe ein Programm für den SIMPLE- Algorithmus in Matlab und bin noch neu in Matlab. Hierfür muss ich die Geschwindigkeiten berechnen, aber es taucht ein Fehler auf, den ich mir nicht erklären kann:

"Error using /
Matrix dimensions must agree.

Error in StroemungfiniteVolumen (line 103)
vstern(i,j)=(av_E*vstern(i+1,j)+av_W*vstern(i- 1,j)+ av_N*vstern(i,j+1) + av_S*vstern(i,j-1)+a_P0*vstern_alt(i,j))/av_P
+ (pstern(i-1,j)-pstern(i,j))*dx/av_P"

Der dazugehörige Code lautet:

Code:

for i=2:ni+1
for j=2:nj+1

%Massenflüsse für Impulsgleichung in x-Richtung
Fu_w=0.5*(Fu(i,j)+Fu(i-1,j));
Fu_e=0.5*(Fu(i+1,j)+Fu(i,j));
Fu_s=0.5*(Fv(i,j)+Fv(i-1,j));
Fu_n=0.5*(Fv(i,j+1)+Fv(i-1,j+1));
%Massenflüsse für Impulsgleichung in y-Richtung
Fv_w=0.5*(Fu(i,j)+Fu(i,j-1));
Fv_e=0.5*(Fu(i+1,j)+Fu(i+1,j-1));
Fv_s=0.5*(Fv(i,j-1)+Fv(i,j));
Fv_n=0.5*(Fv(i,j)+Fv(i,j+1));
%Koeffizienten für x-Impuls
au_E=D_e+max(-Fu_e,0);
au_W=D_w+max(Fu_w,0);
au_N=D_n+max(-Fu_n,0);
au_S=D_s+max(Fu_s,0);
a_P0=rho*dx*dy/dt;
au_P=au_E+au_W+au_N+au_S+a_P0;
%diskretisierte Bestimmungsgleichung für ustern
line 95: ustern(i,j)=(au_E*ustern(i+1,j)+au_W*ustern(i-1,j)+au_N*ustern(i,j+1)+au_S*ustern(i,j-1)+a_P0*ustern_alt(i,j))/au_P+(pstern(i-1,j)-pstern(i,j))*dy/au_P;
%Koeffizienten für y-Impuls
av_E=D_e+max(-Fv,0);
av_W=D_w+max(Fv_w,0);
av_N=D_n+max(-Fv_n,0);
av_S=D_s+max(Fv_s,0);
av_P=av_E+av_W+av_N+av_S+a_P0;
%diskretiserte Bestimmungsgleichung für vstern
line 103: vstern(i,j)=(av_E*vstern(i+1,j)+av_W*vstern(i-1,j)+av_N*vstern(i,j+1)+av_S*vstern(i,j-1)+a_P0*vstern_alt(i,j))/av_P + (pstern(i,j-1)-pstern(i,j))*dx/av_P;
end
end

Funktion ohne Link?

Ich verstehe nicht, wieso die Matrix Dimensionen nicht übereinstimmen, da in Zeile 95 fast das gleiche steht und ich dort keine Fehlermeldung bekomme.

Vielen Dank im Vorraus

PaulM

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 18.09.2014, 22:15 Titel:

Ich habe vergessen, die Definitionen der verschiedenen Matrizen hinzuzufügen, daher:

Code:

%Abmessungen des Berechnungsgitters
Lx=3;
Ly=0.2;
cv_x=150; %Anzahl Kontrollvolumen in x-Richtung
cv_y=10; %Anzahl Kontrollvolumen in y-Richtung
dx=Lx/cv_x; %Breite eines KV
dy=Ly/cv_y; %Höhe eines KV
ni=cv_y+2;
nj=cv_x+2;

%Erstellung der Matrizen für die gesuchten Größen
u=zeros(ni+1,nj+1);
v=zeros(ni+1,nj+1);
p=zeros(ni,nj);
ustern=u;
vstern=v;
pstern=p;
ustrich=u;
vstrich=v;
pstrich=p;

%Erstellung der Matrizen für die Massenflüsse über KV-Wände
%und Diffusionsleitfähigkeiten an KV-Wänden
Fu=zeros(ni+1,nj+1);
Fv=Fu;

Funktion ohne Link?


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.