Skip to content

Commit

Permalink
Added correction factor on heat transfer models
Browse files Browse the repository at this point in the history
  • Loading branch information
Giovanni Mangola committed Feb 11, 2020
1 parent 690799d commit 2565c84
Showing 1 changed file with 27 additions and 26 deletions.
53 changes: 27 additions & 26 deletions ThermoPower/Thermal.mo
View file
Expand Up @@ -959,7 +959,7 @@ The swapping is performed if the counterCurrent parameter is true (default value
if not useAverageTemperature then T[j+1]
else if not adaptiveAverageTemperature then (T[j] + T[j + 1])/2
else (T[j]+T[j+1])/2 + (T[j+1]-T[j])/2*exp(-w_wnom/sigma);
Qw[j] = (Tw[j] - Tvol[j])*omega*l*gamma*Nt;
Qw[j] = (Tw[j] - Tvol[j])*omega*l*kc*gamma*Nt;
end for;

annotation (
Expand Down Expand Up @@ -1025,10 +1025,10 @@ The swapping is performed if the counterCurrent parameter is true (default value
equation
Hv*Qv = Hw*Qw "Energy balance on coarser grid";
if Nw >= Nv then
Qw = omega*lw*gamma*(Tw - G*Tv)
Qw = omega*lw*kc*gamma*(Tw - G*Tv)
"Convective heat transfer on finer grid";
else
Qv = omega*lv*gamma*(Tv - G*Tw)
Qv = omega*lv*kc*gamma*(Tv - G*Tw)
"Convective heat transfer on finer grid";
end if;
for j in 1:Nv loop
Expand Down Expand Up @@ -1094,7 +1094,7 @@ the global nominal thermal conductance UA is given instead of the nominal specif
if not useAverageTemperature then T[j+1]
else if not adaptiveAverageTemperature then (T[j] + T[j + 1])/2
else (T[j]+T[j+1])/2 + (T[j+1]-T[j])/2*exp(-w_wnom/sigma);
Qw[j] = (Tw[j] - Tvol[j])*gamma*omega*l*Nt;
Qw[j] = (Tw[j] - Tvol[j])*kc*gamma*omega*l*Nt;
end for;
annotation (
Diagram(coordinateSystem(preserveAspectRatio=false, extent={{-100,-100},{
Expand Down Expand Up @@ -1175,7 +1175,7 @@ the global nominal thermal conductance UAnom is given instead of the nominal spe
for j in 1:Nw loop
Tvol[j] = if useAverageTemperature then (T[j] + T[j + 1])/2 else T[j + 1];
gamma_vol[j] = if useAverageTemperature then (gamma[j] + gamma[j+1])/2 else gamma[j+1];
Qw[j] = (Tw[j] - Tvol[j])*omega*l*gamma_vol[j]*Nt;
Qw[j] = (Tw[j] - Tvol[j])*kc*omega*l*gamma_vol[j]*Nt;
end for;
annotation (
Diagram(coordinateSystem(preserveAspectRatio=false, extent={{-100,-100},{
Expand Down Expand Up @@ -1282,37 +1282,37 @@ the global nominal thermal conductance UAnom is given instead of the nominal spe

for j in 1:Nw loop
if noEvent((h[j] < hl and h[j + 1] < hl) or (h[j] > hv and h[j + 1]> hv)) then // 1-phase liquid or vapour
Qw[j] = (Tw[j] - Tvol[j])*omega*l*Nt*((gamma1ph[j] + gamma1ph[j+1])/2);
Qw[j] = (Tw[j] - Tvol[j])*omega*l*Nt*((gamma1ph[j] + gamma1ph[j+1])/2)*kc;
state[j] = 1;
alpha_l[j] = 0;
alpha_v[j] = 0;
elseif noEvent((h[j] < hl and h[j + 1] >= hl and h[j + 1] <= hv)) then // liquid --> 2-phase
Qw[j] = alpha_l[j]* (Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*((gamma1ph[j] + gamma_bubble)/2) +
(1 - alpha_l[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph;
Qw[j] = (alpha_l[j]* (Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*((gamma1ph[j] + gamma_bubble)/2) +
(1 - alpha_l[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph)*kc;
state[j] = 2;
alpha_l[j] = (hl - h[j])/(h[j + 1] - h[j]);
alpha_v[j] = 0;
elseif noEvent(h[j] >= hl and h[j] <= hv and h[j + 1] >= hl and h[j + 1]<= hv) then // 2-phase
Qw[j] = (Tw[j] - Ts)*omega*l*Nt*gamma2ph;
Qw[j] = (Tw[j] - Ts)*omega*l*Nt*gamma2ph*kc;
state[j] = 3;
alpha_l[j] = 0;
alpha_v[j] = 0;
elseif noEvent(h[j] >= hl and h[j] <= hv and h[j + 1] > hv) then // 2-phase --> vapour
//Qw[j] = alpha_v[j]*(Tw[j] - (T[j + 1] + Ts)/2)*omega*l*(gamma1ph[j] + gamma1ph[j+1])/2 + (1 - alpha_v[j])*(Tw[j] - Ts)*omega*l*gamma2ph;
Qw[j] = alpha_v[j]* (Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*(gamma_dew + gamma1ph[j+1])/2 +
(1 - alpha_v[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph;
Qw[j] = (alpha_v[j]* (Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*(gamma_dew + gamma1ph[j+1])/2 +
(1 - alpha_v[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph)*kc;
state[j] = 4;
alpha_l[j] = 0;
alpha_v[j] = (h[j + 1] - hv)/(h[j + 1] - h[j]);
elseif noEvent(h[j] >= hl and h[j] <= hv and h[j + 1] < hl) then // 2-phase --> liquid
Qw[j] = alpha_l[j]* (Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*(gamma_bubble + gamma1ph[j+1])/2 +
(1 - alpha_l[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph;
Qw[j] = (alpha_l[j]* (Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*(gamma_bubble + gamma1ph[j+1])/2 +
(1 - alpha_l[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph)*kc;
state[j] = 5;
alpha_l[j] = (hl - h[j + 1])/(h[j] - h[j + 1]);
alpha_v[j] = 0;
else // if noEvent(h[j] > hv and h[j + 1] <= hv and h[j + 1] >= hl) then // vapour --> 2-phase
Qw[j] = alpha_v[j]* (Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*(gamma1ph[j] + gamma_dew)/2 +
(1 - alpha_v[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph;
Qw[j] = (alpha_v[j]* (Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*(gamma1ph[j] + gamma_dew)/2 +
(1 - alpha_v[j])*(Tw[j] - Ts)*omega*l*Nt*gamma2ph)*kc;
state[j] = 6;
alpha_l[j] = 0;
alpha_v[j] = (h[j] - hv)/(h[j] - h[j + 1]);
Expand Down Expand Up @@ -1397,41 +1397,41 @@ the global nominal thermal conductance UAnom is given instead of the nominal spe

for j in 1:Nw loop
if noEvent(h[j] < hl and h[j + 1] < hl) then // liquid
Qw[j] = (Tw[j] - Tvol[j])*omega*l*Nt*gamma_liq;
Qw[j] = (Tw[j] - Tvol[j])*omega*l*Nt*gamma_liq*kc;
state[j] = 0;
alpha_l[j] = 1;
alpha_v[j] = 0;
elseif noEvent(h[j] > hv and h[j + 1]> hv) then // vapour
Qw[j] = (Tw[j] - Tvol[j])*omega*l*Nt*gamma_vap;
Qw[j] = (Tw[j] - Tvol[j])*omega*l*Nt*gamma_vap*kc;
state[j] = 1;
alpha_l[j] = 0;
alpha_v[j] = 1;
elseif noEvent((h[j] < hl and h[j + 1] >= hl and h[j + 1] <= hv)) then // liquid --> 2-phase
Qw[j] = alpha_l[j] *(Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*gamma_liq +
(1 - alpha_l[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph;
Qw[j] = (alpha_l[j] *(Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*gamma_liq +
(1 - alpha_l[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph)*kc;
state[j] = 2;
alpha_l[j] = (hl - h[j])/(h[j + 1] - h[j]);
alpha_v[j] = 0;
elseif noEvent(h[j] >= hl and h[j] <= hv and h[j + 1] >= hl and h[j + 1]<= hv) then // 2-phase
Qw[j] = (Tw[j] - Ts)*omega*l*Nt*gamma_2ph;
Qw[j] = (Tw[j] - Ts)*omega*l*Nt*gamma_2ph*kc;
state[j] = 3;
alpha_l[j] = 0;
alpha_v[j] = 0;
elseif noEvent(h[j] >= hl and h[j] <= hv and h[j + 1] > hv) then // 2-phase --> vapour
Qw[j] = alpha_v[j] *(Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*gamma_vap +
(1 - alpha_v[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph;
Qw[j] = (alpha_v[j] *(Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*gamma_vap +
(1 - alpha_v[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph)*kc;
state[j] = 4;
alpha_l[j] = 0;
alpha_v[j] = (h[j + 1] - hv)/(h[j + 1] - h[j]);
elseif noEvent(h[j] >= hl and h[j] <= hv and h[j + 1] < hl) then // 2-phase --> liquid
Qw[j] = alpha_l[j] *(Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*gamma_liq +
(1 - alpha_l[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph;
Qw[j] = (alpha_l[j] *(Tw[j] - (T[j + 1] + Ts)/2)*omega*l*Nt*gamma_liq +
(1 - alpha_l[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph)*kc;
state[j] = 5;
alpha_l[j] = (hl - h[j + 1])/(h[j] - h[j + 1]);
alpha_v[j] = 0;
else // if noEvent(h[j] > hv and h[j + 1] <= hv and h[j + 1] >= hl) then // vapour --> 2-phase
Qw[j] = alpha_v[j] *(Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*gamma_vap +
(1 - alpha_v[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph;
Qw[j] = (alpha_v[j] *(Tw[j] - (T[j] + Ts)/2)*omega*l*Nt*gamma_vap +
(1 - alpha_v[j])*(Tw[j] - Ts) *omega*l*Nt*gamma_2ph)*kc;
state[j] = 6;
alpha_l[j] = 0;
alpha_v[j] = (h[j] - hv)/(h[j] - h[j + 1]);
Expand Down Expand Up @@ -2115,6 +2115,7 @@ This package contains models to compute the material properties needed to model
parameter SI.MassFlowRate wnom "Nominal mass flow rate (single tube)"
annotation(Dialog(enable=false, tab = "Set by Flow1D model"));
final parameter SI.Length l=L/(Nw) "Length of a single volume";
parameter Modelica.SIunits.PerUnit kc = 1 "Correction factor for heat transfer";

Medium.Temperature T[Nf] "Temperatures at the fluid side nodes";
Medium.Temperature Tw[Nw] "Temperatures of the wall volumes";
Expand Down

0 comments on commit 2565c84

Please sign in to comment.