def Parameters(vwind, Tlist, mf, del_z):
A = Constants()[0]
hgam = hgam_fun(vwind, Tlist, del_z)
hga = hga_fun(Tlist, del_z)
haab = haab_fun(Tlist, del_z)
hf = hf_fun(mf, Tlist, del_z)
hr1 = sigma / (1 / ep + 1 / eg - 1)
Vol_correct = del_z / Lc
# print(hgam, hga, haab, hf, hr1)
b = A / Lc
B = Cg * rhog * Vg * Vol_correct
C = hgam * b * del_z
D = hr1 * b * del_z
E = hga * b * del_z
F = tau * alpha_sg * b * del_z
H = Ca * rhoa * Va * Vol_correct
I = haab * b * del_z
J = alpha_sab * b * del_z
K = Cab * rhoab * Vab * Vol_correct
L = ki * b * del_z / Li
M = np.pi * din * del_z * hf * 8 #remember plate number
N = Ci * rhoi * Vi * Vol_correct
O = Cf * rhof * (np.pi / 4 * din**2) * del_z * 8 #remember plate number
P = mf * Cf
Q = sigma * b * del_z
return [B, C, D, E, F, H, I, J, K, L, M, N, O, P, Q]
def Integration_step(G, vwind, mf, Tlist, Tf_prev, del_t, del_z):
A, alpha_sab, alpha_sg, Ca, Cab, Cf, Cg, Ci, din, eg, eg, ep, g, ka, kf, \
ki, Lc, Li, Ltube, rhoa, rhoab, rhof, rhog, rhoi, sigma, tau, Va, Vab,\
Vg, Vi = Constants()
Tg, Ta, Tab, Tf, Ti, Tam = Tlist
B, C, D, E, F, H, I, J, K, L, M, N, O, P, Q = Parameters(
vwind, Tlist, mf, del_z)
Tsky = 0.0552 * Tam**1.5
G = G / 2.9
dTgdt = C/B*(Tam-Tg) + D/B*(Tab**4-Tg**4) + E/B*(Ta-Tg) + F/B*G \
+ Q/B*(Tsky**4 - Tg**4)
dTadt = E / H * (Tg - Ta) + I / H * (Tab - Ta)
dTabdt = J/K*G + D/K*(Tg**4-Tab**4) + I/K*(Ta-Tab) + L/K*(Ti-Tab) \
+ M/K*(Tf-Tab)
dTfdt = M / O * (Tab - Tf) - P / O * (Tf - Tf_prev)
dTidt = L / N * (Tab - Ti) + C / N * (Tam - Ti)
Tg += del_t * dTgdt
Ta += del_t * dTadt
Tab += del_t * dTabdt
Tf += del_t * dTfdt
Ti += del_t * dTidt
return [Tg, Ta, Tab, Tf, Ti]
# -*- coding: utf-8 -*-
"""
Created on Thu Feb 16 17:02:12 2017
@author: Wian
"""
import numpy as np
from Constants_main import Constants
A, alpha_sab, alpha_sg, Ca, Cab, Cf, Cg, Ci, din, eg, eg, ep, g, ka, kf, ki, Lc, Li, Ltube, rhoa, rhoab, rhof, rhog, rhoi, sigma, tau, Va, Vab, Vg, Vi = Constants(
)
def hgam_fun(v_wind, Tlist, del_z): #hiam has same function
Tg, Ta, Tab, Tf, Ti, Tam = Tlist
mu = 1.849E-5
Re = rhoa * v_wind * Lc / mu
Pr = 0.7202
if Re > 5E5:
Nu = 0.037 * Re**0.8 * Pr**(1 / 3)
else:
Nu = 0.664 * Re**0.5 * Pr**(1 / 3)
return Nu * ka / Lc
def hf_fun(mf, Tlist, del_z):
Tg, Ta, Tab, Tf, Ti, Tam = Tlist
muf = 8.9E-4
u = mf / (np.pi / 4 * din**2) / rhof / 8
Re = rhof * u * din / muf
from Integration_main import Integration_step
from Data_Experimental import Import_data
from Constants_main import Constants
import matplotlib.pyplot as plt
import numpy as np
A, alpha_sab, alpha_sg, Ca, Cab, Cf, Cg, Ci, din, eg, eg, ep, g, ka, kf, ki, \
Lc, Li, Ltube, rhoa, rhoab, rhof, rhog, rhoi, sigma, tau, Va, Vab, Vg, Vi \
= Constants()
Nodes = 1
Times, Tin, Tout, G, Air_Temp = Import_data()
Tin = np.array(Tin) + 273.15
Tout = np.array(Tout) + 273.15
def mass_flow(t):
if t < 34*60:
return 0.01296
elif t < 55*60:
return 0.0327
else:
return 0.007788
for Nodes in [20]:
del_z = Ltube/Nodes
Temps = np.array([35, 35, 25, 37, 35]) + 273.15
iterations = 2000
tspan = np.linspace(0, Times[-1], iterations)
del_t = (tspan[1] - tspan[0])