SpheresPacking
The function
SpheresPacking(Fluid$,m_dot, d, phi, A_fr, L, T_i, T_o, P, Cor: h, DELTAp, A_s, DELTAp_c, alpha, f, j_H)
returns the fluid side behavior of a packed matrix of spheres.
Inputs:
Fluid$ - string variable indicating a specific fluid in EES database or in the Solid-Liquid_Props fluids list
m_dot - mass flow rate of fluid [kg/s] or [lbm/min]
d - diameter of the spheres [m] or [ft]
phi - porosity of the bed [-]
A_fr - the frontal area of the matrix exposed to fluid flow [m^2] or [ft^2]
L - length/depth of stacked screens [m] or [ft]
T_i - inlet temperature of the fluid [C], [K], [F], [R]
T_o - outlet temperature of the fluid [C], [K], [F], [R]
P - pressure of the fluid at the inlet [Pa], [kPa], [bar], [MPa], [atm], or [psia]
Cor - correlation to use (1 = Kays and London)
Ouputs:
h - heat transfer coefficient [W/m^2-K] or [Btu/hr-ft^2-R]
DELTAp - total pressure drop [Pa], [kPa], [bar], [MPa], [atm], or [psia]
A_s - surface area [m^2] or [ft^2]
DELTAp_c - core-only pressure drop [Pa], [kPa], [bar], [MPa], [atm], or [psia]
alpha - surface area per volume [m^2/m^3] or [ft^2/ft^3]
f - friction factor [-]
j_H - Stanton j-factor [-]
Notes: This function determines the fluid properties at the stated pressure and temperature, and the Reynolds number from the sphere diameter. The non-dimensional function PackedSpheres_ND is called to determine the values for the Colburn j and friction factor. The Colburn j factor is interpreted by PackedSpheres and returned as a heat transfer coefficient. The porosity for a theoretically perfect spherical matrix, i.e., all the spheres are exactly the same diameter and are packed as close as possible, is 0.30 but porosity usually ranges from 0.32 to 0.37. The core pressure drop is determined based on frictional considerations. The total pressure drop also includes fluid acceleration.
Example:
$UnitSystem SI Mass J K Pa Radian
$VarInfo alpha units=m^2/m^3
$VarInfo A_s units=m^2
$VarInfo DELTAp units=Pa
$VarInfo DELTAp_c units=Pa
$VarInfo h units=W/m^2-K
Fluid$='air'
m_dot=0.1 [kg/s]
d=0.002 [m]
phi=0.36
A_fr=0.01 [m^2]
L=0.01 [m]
T_i=300 [K]
T_o=320 [K]
P=1e6 [Pa]
Cor=1
Call spherespacking(Fluid$,m_dot, d, phi, A_fr, L, T_i, T_o, P, Cor: h, DELTAp, A_s, DELTAp_c, alpha, f, j_H)
{Solution:
alpha=1920 [m^2/m^3]
A_s=0.192 [m^2]
DELTAp=932.4 [Pa]
DELTAp_c=929.9 [Pa]
f=0.5079
h=974.9 [W/m^2-K]
j_H=0.02818
}