modelBellDelaware

BellDelaware: Single phase (oil/water) correlation for flow in the shell side of a STHX

Extends from TRANSFORM.Fluid.ClosureRelations.HeatTransfer.Models.DistributedPipe_1D.PartialHeatTransfer_setT (Base model).

Information

The Lyon-Martinelli heat transfer correlation is for flow of liquid metal in circular tubes with contant heat flux along the tube wall.

The correlation was taken from the following textbook:

M. M. Wakil

Nuclear Heat Transport 1993

eq. 10-2 pg. 268

and is reported valid for:

  • liquid metal coolant
  • flow in circular tubes
  • constant heat flux along tube wall

Parameters

TypeNameDefaultDescription
BooleanisGasfalsetrue if Medium is a gas
Realnp0Gas specific exponential correction factor (e.g., air = 0; N2 = 0.12)
BooleantoggleEndChannelfalse=true for end channel; false for center channel
BooleantoggleStaggeredtruetrue = staggered grid type; false = in-line
SI.Lengthd_BDiameter of holes in baffles
SI.Lengthd_oOuter diameter of tubes
SI.LengthD_iInside shell diameter
SI.LengthD_lBaffle Diameter
SI.LengthDBTube bundle diameter
SI.LengthHHeight of baffle cut
SI.Lengths1Tube to tube pitch parallel to baffel edge
SI.Lengths2Tube to tube pitch perpindicular to baffel edge
SI.LengthSBaffle spacing between baffles
SI.LengthS_EBaffle spacing between the heat exchanger sheets and adjacent baffles
SI.Lengthe1Space between tubes and shell
Realnes# of shortest connections connecting neighboring tubes
Realn_W# of tubes in both the upper and lower window
Realn_TTotal # of tubes (including blind and support)
Realn_MR# of main resistances in cross flow path
Realn_s# of pairs of sealing strips
Internal Interface
RealnParallel (from PartialHeatTransfer_setT)1Number of parallel components
IntegernHT (from PartialHeatTransfer_setT)1Number of heat transfer segments
IntegerflagIdeal (from PartialHeatTransfer_setT)0Flag for models to handle ideal boundary
Internal Interface › Initialization
SI.Temperature[nHT]Ts_start (from PartialHeatTransfer_setT)
Advanced
SI.ReynoldsNumberRe_lam (from PartialHeatTransfer_setT)2300Laminar transition Reynolds number
SI.ReynoldsNumberRe_turb (from PartialHeatTransfer_setT)4000Turbulent transition Reynolds number
Units.NonDim[nHT]CF (from PartialHeatTransfer_setT)fill(1.0, nHT)Correction Factor: Q = CF*alpha*A*dT

Connectors

TypeNameDefaultDescription
HeatAndMassTransfer.Interfaces.HeatPort_Flow[nHT]heatPorts (from PartialHeatTransfer_setT)

Components

TypeNameDefaultDescription
Medium.ThermodynamicState[nHT]states (from PartialHeatTransfer_setT)Thermodynamic state of fluid
SI.Velocity[nHT]vs (from PartialHeatTransfer_setT)Fluid Velocity
SI.Diameter[nHT]dimensions (from PartialHeatTransfer_setT)Characteristic dimension (e.g. hydraulic diameter)
SI.Area[nHT]crossAreas (from PartialHeatTransfer_setT)Cross sectional flow area
SI.Length[nHT]dlengths (from PartialHeatTransfer_setT)Characteristic length of heat transfer segment
SI.Height[nHT]roughnesses (from PartialHeatTransfer_setT)Average height of surface asperities
SI.Area[nHT]surfaceAreas (from PartialHeatTransfer_setT)Surface area for heat transfer
SI.Temperature[nHT]Ts_fluid (from PartialHeatTransfer_setT)Medium.temperature(states)Fluid temperature
SI.Temperature[nHT]Ts_wall (from PartialHeatTransfer_setT)heatPorts.TWall temperature
SI.MassFlowRate[nHT]m_flows (from PartialHeatTransfer_setT)Fluid mass flow rate
SI.ReynoldsNumber[nHT]Res (from PartialHeatTransfer_setT)Reynolds number
SI.PrandtlNumber[nHT]Prs (from PartialHeatTransfer_setT)Prandtl number
SI.CoefficientOfHeatTransfer[nHT]alphas (from PartialHeatTransfer_setT)Coefficient of heat transfer
SI.NusseltNumber[nHT]Nus (from PartialHeatTransfer_setT)Nusselt number
SI.HeatFlowRate[nHT]Q_flows (from PartialHeatTransfer_setT)heatPorts.Q_flow/nParallelHeat flow rate
Realas1/d_o
Realbs2/d_o
Realc((a/2)^2 + b^2)^(0.5)
SI.Lengthe(if toggleStaggered then (if b >= 0.5*(2*a + 1)^(0.5) then (a - 1)*d_o else (c - 1)*d_o) else (a - 1)*d_o)
SI.LengthL_E2*e1 + e*nes
SI.AreaA_Eif toggleEndChannel then S*L_E else S_E*L_E
SI.AreaA_fif toggleEndChannel then D_i*S_E else D_i*S
SI.AreaA_Bif e < (D_i - DB) then S*(D_i - DB - e) else 0
RealR_BA_B/A_E
RealR_Sn_s/n_MR
SI.ReynoldsNumberRes_psil
Realbeta
Realf_B
Realf_w
Realf_P
RealPrs_w
SI.NusseltNumberNus_llam
SI.NusseltNumberNus_lturb
SI.NusseltNumberNus_l0
SI.NusseltNumberNus_bundle0
SI.NusseltNumberNus_bundle
Realgamma2*Modelica.Math.acos(1 - 2*H/D_l)*180/pi
SI.AreaA_GTB(n_T - n_W/2)*pi*(d_B^2 - d_o^2)/4
SI.AreaA_GSBpi/4*(D_i^2 - D_l^2)*(360 - gamma)/360
SI.AreaA_SGA_GTB + A_GSB
RealR_LA_SG/A_E
RealR_Gn_W/n_T
Realf_Lif A_SG < 1e-6 then 1 else 0.4*A_GTB/A_SG + (1 - 0.4*A_GTB/A_SG)*exp(-1.5*R_L)
Realf_G1 - R_G + 0.524*R_G^(0.32)
Realpsiif b >= 1 then 1 - pi/(4*a) else 1 - pi/(4*a*b)
Realf_Aif toggleStaggered then 1 + 2/(3*b) else 1 + 0.7*(b/a - 0.3)/(psi^(1.5)*(b/a + 0.7)^2)
Realf_N1
SI.Lengthllpi/2*d_o
SI.Velocityw
Medium.ThermodynamicStatestates_w
SI.PressurePsMedium.pressure(states)