Energy Flows in Thermal Liquid Networks

Upwind Energy Scheme

Thermal energy advection should always match the direction of flow. To satisfy this condition, the Thermal Liquid domain employs a modeling framework known as the upwind scheme. This scheme works by matching the specific internal energy at a port to that just upwind of the port:

If the flow at a port is directed into the component, the specific internal energy at that port is equal to that carried by the incoming fluid. In the figure, the specific internal energy at port A of component2等于在港口的组件吗1.
If the fluid flow at a port is directed out of the component, the specific internal energy at that port is equal to that in the component’s internal fluid volume. In the figure, the specific internal energy at port A of component2is equal to that in the internal fluid volume of the same component.

The upwind scheme describes the advection thermal energy flow rate using a conditional expression that, for port A of some component, reads:

$Φ_{A}^{A d} = {\begin{matrix} {\dot{m}}_{A} u_{A}, & {\dot{m}}_{A} > 0 \\ {\dot{m}}_{A} u_{I}, & {\dot{m}}_{A} \leq 0 \end{matrix},$

where:

Φ_A^Adis the thermal energy flow rate due to advection through port A.
${\dot{m}}_{A}$ is the mass flow rate through port A.
u_Ais the specific internal energy upwind of port A.
u_I具体的内部能源内部fluid volume.

Without numerical smoothing, this expression introduces a number of computational challenges. It adds a slope discontinuity to the thermal energy flux at zero mass flow rates that makes the model more prone to simulation errors. It also adds a jump discontinuity to the specific internal energy and makes its value ill-defined during flow reversals.

Upwind Energy Scheme without Smoothing

Numerical Smoothing

The Thermal Liquid domain smoothes the numerical discontinuities introduced by the upwind energy scheme by adding to the thermal energy flow rate a thermal conduction term. At port A of some component:

$Φ_{A}^{T h} = {\begin{matrix} {\dot{m}}_{A} u_{A} + G (u_{A} - u_{I}), & {\dot{m}}_{A} > 0 \\ {\dot{m}}_{A} u_{I} + G (u_{A} - u_{I}), & {\dot{m}}_{A} \leq 0 \end{matrix},$

where:

Φ_A^This the thermal energy flow rate due to advection through port A and conduction between port A and the internal fluid volume.
Gis a thermal conductance coefficient computed from the fluid properties and component geometry.

Rewriting the modified expression in a more convenient form:

$Φ_{A}^{T h} = {\begin{matrix} ({\dot{m}}_{A} + G) u_{A} - G u_{I}, & {\dot{m}}_{A} > 0 \\ G u_{A} + ({\dot{m}}_{A} - G) u_{I}, & {\dot{m}}_{A} \leq 0 \end{matrix}$

Usingminandmaxfunctions to collapse the conditional expression:

$Φ_{A}^{T h} = [\max ({\dot{m}}_{A}, 0) + G] u_{A} + [\min ({\dot{m}}_{A}, 0) - G] u_{I}$

Approximating the collapsed expression as a numerically smooth expression:

$Φ_{A}^{T h} \approx \frac{{\dot{m}}_{A} + \sqrt{{\dot{m}}_{A}^{2} + 4 G^{2}}}{2} u_{A} + \frac{{\dot{m}}_{A} - \sqrt{{\dot{m}}_{A}^{2} + 4 G^{2}}}{2} u_{I}$

The end expression removes the slope discontinuity from the thermal energy flow rate curve and the jump discontinuity from the specific internal energy curve. The thermal conductance coefficient determines the amount of smoothing applied in both cases:

$G = \frac{k S}{c_{v} L},$

where:

kis the thermal conductivity defined in theThermal Liquid Settings (TL)block.
Sis the flow cross-sectional area at the port considered.
c_vis the specific heat defined in theThermal Liquid Settings (TL)block.
Lis a characteristic distance between the port and the component’s fluid volume.

Upwind Energy Scheme with Smoothing

Total Energy Flow Rate

The Through variable of the Thermal Liquid domain is the total energy flow rate. This variable is defined in terms of the smoothed upwind thermal energy flow rate as:

$Φ_{A} = Φ_{A}^{T h} + \frac{p}{ρ},$

where:

Φ_Ais the total energy flow rate through port A.
Φ_A^This the thermal energy flow rate through port A computed from the smoothed upwind energy scheme.
pis the absolute pressure at port A.
ρis the fluid density at port A.

The kinetic energy contribution to the total energy flow rate is assumed negligible in Thermal Liquid networks and is ignored in this domain.