MATLAb Tremine Code Druck Diffusion in TKE Budget

[Samoza]

Bye to all,

Ich versuche, alle Begriffe des Tke-Haushalts in einem einzigen Diagramm aufzuzeigen, um die Kurven mit den Daten aus Moser, Kim und Mansour (1999) zu vergleichen. die Formel des Tke-Haushalts wird im beigefügten Bild dargestellt:

https://imgur.com/d9mbt1

die betrachtete Strömung ist eine turbulente Strömung im Kanal. Ich konnte fast alle Begriffe (Produktion, Turbulenz Diffusion, viskose Diffusion, Dissipation Rate), erhalten konsistente Ergebnisse mit den Daten von Moser, Kim und Mansour wie im Bild gesehen:

https://imgur.com/9w9gfz5

Ich habe jedoch Schwierigkeiten, den Begriff Druckdiffusor korrekt aufzuzeigen.

verfügbare Daten enthalten:

- drei Matrizen 3d (256x128x128) `u_prime`, `v_prime` und `w_prime`, die die Schwankungen der Geschwindigkeitskomponenten entlang der Richtungen x, y bzw. z darstellen.
- eine 3d `p_prime` Matrix (256x128x128) die das Druckfeld darstellt.
- drei Matrizen 3d `x`, `y` und `z` (256x128x128) die räumlichen Koordinaten darstellen.

moser, kim und mansour data finden Sie auf der folgenden Website: dns Daten für turbulenten Kanal fließen bis zu Reibung Reynolds Anzahl 590 . die spezifische Datei, die die Daten zum Vergleich enthält, befindet sich auf dem Pfad `chan180/balances/chan180.kbal`.

hier ist mein Versuch:

cc;
klar;
alle schließen;


% Raumstatistik
% Importieren von Zellenzentren Koordinaten
ccx = importdata(fullfile('velfieldschannel_corsoturb', 'ccx'));
ccy = importdata(fullfile('velfieldschannel_corsoturb', 'ccy'));
ccz = importdata(fullfile('velfieldschannel_corsoturb', 'ccz'));


% Importgeschwindigkeitsfeld
Wähle einen Snapshot = 3; % wähle einen Snapshot (von 1 bis 22)
Dateipfad = fullfile('velfieldschannel_corsoturb', num2str(choosesnapshot), 'u.gz');
tempdir = fullfile(tempdir, 'tempextraction');
gunzip(filepath, tempdir);
u = Importdaten(fullfile(tempdir, 'u'));


% Importdruckfeld
pfield = importdata(fullfile('velfieldschannel_corsoturb', 'p'));


% Parameter
nx = 256;
= 128;
nz = 128;
Delta = 1;
u_tau = 1;
n = 1/180;
re = (delta * u_tau) / nu;
rho = 1


% Zuordnungswerte in Matrizen
u = Nullen (nx, ny, nz);
v = Null(nx, ny, nz);
w = Nullen (nx, ny, nz);
x = Null(nx, ny, nz);
y = Null(nx, ny, nz);
z = Null(nx, ny, nz);
p = Null(nx, ny, nz);
Index = 1;
für k = 1:nz
für j = 1:ny
für i = 1:nx
u(i, j, k) = u(index, 1);
v(i, j, k) = u(Index, 2.
w(i, j, k) = u(Index, 3);
x(i, j, k) = ccx(index);
y(i, j, k) = ccy(index);
z(i, j, k) = ccz(index);
p(i, j, k) = pfield(index);
Index = Index + 1;
Ende
Ende
Ende


% Berechnung der mittleren Profile für jede y
u_mean = mean(u, [1, 3];
v_mean = mittel(v) [1, 3];
w_mean = mittel(w) [1, 3];
p_mean = mittel(p) [1, 3];


% Berechnungsfelder
u_prime = u - u_mean;
v_prime = v - v_mean;
w_prime = w - w_mean;
p_prime = p - p_mean;


%extrahierende Positionsvektoren
x_direc = x, 1, 1;
y_direc = y(1, :, 1);
z_direc = z(1, 1, : );


% Berechnung y_plus
y_plus = (y_direc * u_tau) / nu;


% Berechnung unterschiedlicher Richtungen x, y und z
dx = diff(x_direc);
dy = diff(y_direc);
dz = diff(z_direc);


% Import von Moos, Kim und Mansour-Daten
filepath = fullfile('chan180', 'balances', 'chan180.kbal');
Dateiid = fopen(filepath, 'r));
chan180kbal = textcan(fileid, '%f%f%f%f%f%f%f%f%f%f%f%f', 'commentstyle', '#');
fclose(fileid);
y_plusmkm = chan180kbal{2};
dissipationmkm = chan180kbal{3};
pressdiffmkm = chan180kbal{6};


% Druckdiffusionsterm
mean_pu_prime = mean(p_prime .* u_prime, [1 3];
mid_pv_prime = mean(p_prime.* v_prime, [1 3];
(p _ prime. * w _ prime, [1 3];


dpu_dx = diff(mean_pu_prime) ./ dx;
dpv_dy = diff(mean_pu_prime) ./ dy;
dpw_dz = diff(mean_pu_prime) ./ dz;


pressdiff_pre = -(1/rho) * (dpu_dx + dpv_dy + dpw_dz);
pressdiff_pre = pressdiff_pre(1,:, 1);
Pressdiff = nan(1, 128),
pressdiff(1:127) = pressdiff_pre;




% Plotdruckdiffusionsbedingungen
Abbildung;
halten;
Grundstück(y_plus, pressdiff/re, 'displayname', '[imath]Warenbezeichnung[/imath]„;
Grundstück(y_plusmkm, pressdiffmkm, '-', 'displayname', '[imath]{mkm)[/imath]„;
Xlabel( '[imath]+[/imath]',' Dolmetscher ',' Latex ';
Xlim ([0 max(y_plusMKM)/2];
Plottitle = 'Druckdiffusionsterm';
Titel(plottitle, 'interpreter', 'latex');
Legende('show', 'location', 'northeast', 'interpreter', 'latex');
Netz auf;
Feld an;
abhalten;

die Ausgabe dieses Codes ist das folgende Bild:

https://imgur.com/twbq9cs

und wie man sieht, passen die beiden Kurven nicht.

da die Formel für Druckdiffusion im Bild gezeigt wird, wie würden Sie den Code schreiben, um den Begriff Druckdiffusion richtig zu erfassen?

für jede Klarstellung, zögern Sie nicht zu fragen, und ich werde versuchen, zu antworten.