Fluid Flow Loss Analysis Using Real Gas Properties

AI Thread Summary
The discussion focuses on modeling fluid flow and losses in a component using real gas properties, specifically for refrigerants like R134A. The author emphasizes the inadequacy of existing analyses that assume ideal gas behavior, which does not account for density changes and associated losses. They present a 1D analysis method that involves adjusting mass flow rates to maintain specific enthalpy stagnation at the inlet and outlet. Key calculations include determining fluid properties at both the inlet and exit, as well as the impact of drag forces on internal energy. The conversation highlights the need for clarity on system specifics and suggests that a mechanical energy balance may be more appropriate than a thermal energy balance for this analysis.
ScottAllenRauch
Messages
1
Reaction score
0
I am trying to model fluid flow and losses through a component (e.g., pipe) using REAL GAS properties since 1) I have access to RefProp from NIST, and 2) I am dealing with refrigerants, which are far from ideal.

I have seen nowhere an analysis of fluid flow losses (e.g., drag, friction, incidence, pressure drop) using real gas properties. In every discussion and book I have examined, somewhere in the analysis the gas is assumed to be ideal or perfect or steam, or no velocity change is allowed (which happens due to density change, if nothing else).

However, I have access to real gas properties (as I am sure many of us do), and I would greatly appreciate someone with a firmer grasp on reality than I have double-checking my 1D analysis.

The procedure in general is:
Given inlet fluid conditions, and pressure out, change mass flow rate until h0_out = h0_in (i.e., specific enthalpy stagnation at component exit = specific enthalpy stagnation at component inlet).

The way I figure it, all the energy to accelerate the fluid (due to the change in density) and overcome the losses (shear forces due to viscosity) must come from the internal energy, resulting in a reduction in fluid temperature.

I know this is long, but all the formulas and results (including units) are given, so I hope it is at least clear. Thanks in advance.
Code:
Inputs
      Fluid                 = R134A

      T_in_C                Temperature at component inlet
                            = 5.0000  C

      Cd                    Drag coefficient
                            = 0.50000

      Area                  Characteristic area
                            = 0.10000  m²

      PressureRatio         p_out / p_in
                            = 0.90000

Calculations       ƒ()  =>  fluid function RefProp from NIST
  Calculations Not Iterated
    Fluid Properties at Component Inlet

      T_in                  Temperature at component inlet
                            = T_in_C converted to K
                            = 278.15  K

      p_in                  Pressure at component inlet
                            = Saturation pressure @ (T_in - 0.01 K) to ensure pure vapor = ƒ(T)
                            = 349.54  kPa

      rho_in                Density at component inlet
                            = ƒ(P_in,T_in)
                            = 17.124  kg/m³

      s_in                  Specific entropy at component inlet
                            = ƒ(P_in,T_in)
                            = 1.7245  m²/(s²·K)

      h_in                  Specific enthalpy at component inlet
                            = ƒ(P_in,T_in)
                            = 401.50  m²/s²

      u_in                  Specific internal energy at component inlet
                            = ƒ(P_in,T_in)
                            = 381.08  m²/s²

      u_in_reality_check    u_in calculated differently
                            = h_in - p_in / rho_in
                            = 381.08  m²/s²

    Fluid Properties at Component Exit

      p_out                 Pressure at component exit
                            = p_in · PressureRatio
                            = 314.58  kPa

      u_out_s               Specific internal energy at component exit isentropic
                            = ƒ(P_out,s_in)
                            = 379.05  m²/s²

  Iterated Calculations       Solve for mdot

      mdot                  Mass flow rate. Note circular iterative self-reference
                            = mdot · h0_target / h0_out_calculated
                            = 3.9505  kg/s

      rho_out               Density at component exit. Uses u_out calculated below.
                            = ƒ(P_out,u_out)
                            = 15.366  kg/m³      Vdot_in               Volumetric flow rate at component inlet
                            = mdot / rho_in
                            = 0.23069  m³/s

      c_in                  Relative fluid velocity at component inlet
                            = Vdot_in / Area
                            = 2.3069  m/s

      ke_in                 Specific kinetic energy at component inlet
                            = c_in² / 2
                            = 2.6610  m²/s²

      h0_in                 Specific enthalpy stagnation at component inlet
                            = h_in + ke_in
                            = 404.16  m²/s²

      h0_out_target         Specific enthalpy stagnation at component exit
                            = h0_in
                            = 404.16  m²/s²

      Vdot_out              Volumetric flow rate at component exit
                            = mdot / ρout
                            = 0.25709  m³/s

      c_out                 Relative fluid velocity at component exit
                            = Vdot_out / Area
                            = 2.5709  m/s

      ke_out                Specific kinetic energy at component exit
                            = c_out² / 2
                            = 3.3047  m²/s²

      F_drag                Drag force
                            = Cd · rho_in · c_in² · area / 2
                            = 2.2784  kg·m/s²

      Pwr_drag              Drag power
                            = F_drag · c_in
                            = 5.2561  kg·m²/s³

      Delta_u               Δ specific internal energy due to drag
                            = Pwr_drag / mdot
                            = 1.3305  m²/s²

      u_out                 Specific internal energy at component exit
                            = u_out_s + Delta_u
                            = 380.38  m²/s²

      h0_out_calculated
                            = u_out + p_out / rho_out + ke_out
                            = 404.16  m²/s²/code]
 
Last edited:
Engineering news on Phys.org
In order to answer the OP's questions, we need to know more about the specific system being analyzed, including geometry and imposed operating conditions. Offhand, it looks like the OP is trying to use the thermal energy balance when he should be using the mechanical energy balance equation.
 
How did you find PF?: Via Google search Hi, I have a vessel I 3D printed to investigate single bubble rise. The vessel has a 4 mm gap separated by acrylic panels. This is essentially my viewing chamber where I can record the bubble motion. The vessel is open to atmosphere. The bubble generation mechanism is composed of a syringe pump and glass capillary tube (Internal Diameter of 0.45 mm). I connect a 1/4” air line hose from the syringe to the capillary The bubble is formed at the tip...
Thread 'Physics of Stretch: What pressure does a band apply on a cylinder?'
Scenario 1 (figure 1) A continuous loop of elastic material is stretched around two metal bars. The top bar is attached to a load cell that reads force. The lower bar can be moved downwards to stretch the elastic material. The lower bar is moved downwards until the two bars are 1190mm apart, stretching the elastic material. The bars are 5mm thick, so the total internal loop length is 1200mm (1190mm + 5mm + 5mm). At this level of stretch, the load cell reads 45N tensile force. Key numbers...
I'd like to create a thread with links to 3-D Printer resources, including printers and software package suggestions. My motivations are selfish, as I have a 3-D printed project that I'm working on, and I'd like to buy a simple printer and use low cost software to make the first prototype. There are some previous threads about 3-D printing like this: https://www.physicsforums.com/threads/are-3d-printers-easy-to-use-yet.917489/ but none that address the overall topic (unless I've missed...
Back
Top