VVF12 - Fluid humidity calculation in a cylinder

Description

This test case examines the mixing of water vapor in air inside a cylinder and determines temperature and relative humidity of the mixture at the exit of the cylinder. Two streams of wet air with different values of relative humidity, inlet temperature, and mass flow rate are injected into the cylinder shown in the figure below. A high level of mixing is induced in the cylinder so that the fluid becomes homogeneously mixed before reaching the cylinder outlet.

Geometry

The geometry is made up of a cylinder with 750 mm length and 3 mm in diameter. One quadrant face is defined, and then extruded. The remaining 3 quadrants are formed by revolving the first quadrant. This results in a cylinder with 4 quadrants of equal dimensions shown in the following figure.

Simulation Model

Four separate meshes are defined on the inlet end of the cylinder, so that each quadrant has its own mesh, allowing different mass flow rate, humidity levels and inlet temperature to be assigned to each quadrant. The swept mesh is used to create a hexahedral elements mesh along the length of the cylinder. The two quadrants from opposite ends of the cylinder are selected as the source and the target faces. Mesh controls of 3 elements per edge per quadrant side are applied, and can be seen in the figure below. Mesh controls of 75 per edge are applied to the 5 edges running the length of the cylinder.

The fluid is modeled using incompressible air with the following properties:

• Mass density: ρ= 1.207 kg/m³
• Thermal conductivity: k = 0.0263 W/m·K
• Dynamic viscosity: µ = 1.85e-005 kg/m·s
• Specific heat at constant pressure: C_p= 1007 J/kg·K
• Coefficient of thermal expansion: β = 3.41e-006 m/K
• Gas constant: R = 287 J/kg·K

  Two wet air streams (A and B) enter at one end of the cylinder. Each stream has two inlet regions denoted by 1 and 2 in the figure below. The two regions of a single stream are diagonally opposite each other. This insures a good mixing of moist air, while keeping the geometry simple. Both streams enter with a swirl of 45 degree but in opposite directions.

  The following boundary conditions are applied:

  • Flow Boundary Condition: Inlet Flow on four surfaces of the cylinder
  • Flow Boundary Condition: Opening on the four surfaces at the opposite end of the cylinder
  

  Case	ΦA (%)	ṁA (kg/s)	TA (°C)	ΦB (%)	ṁB (kg/s)	TB (°C)	Mass ratio
  1	80	1.0e-4	13	20	5.0e-5	36	2/1
  2	60	6.05e-5	25	10	1.1e-4	42	6/11
  3	90	5.0e-5	17	10	5.0e-5	21	1/1
  4	60	1.3e-4	8	20	3.0e-5	40	13/3

  where:

  • T_A and T_B are the inlet temperature of the stream A and B respectively.
  • ṁ_A and ṁ_B are the mass flow rate of the stream A and B respectively.
  • Φ_A and Φ_B are the relative humidity of the stream A and B respectively (ratio of the mass of water vapor to the mass of saturated vapor).

The following solution options are set:

• Turbulence Model: Mixing Length
• Solution Type: Steady State

The default solver parameters are used.

Theory

The state of the homogeneous fluid mixture resulting from the mixing of two streams of wet air can be obtained theoretically with the use of a psychrometric charts as shown in the following figures. The psychrometric chart allows all the parameters of the air mixture to be determined knowing initial temperature, relative humidity, and mass flow rate of each stream of air at the inlet point [28].

Adiabatic mixing of two air stream is governed by the mass and energy balances:

• Mass of dry air is defined as:

• Mass of water vapor is defined as:

• Energy is defined as:

where:

• ṁ_A, ṁ_B, ṁ_M are the flow rate mass of the stream A, stream B, and mixture M respectively.
• ω_A, ω_B, and ω_M are the specific humidity of the stream A, stream B, and mixture M respectively (ratio of the mass of water vapor to the mass of dry air).
• h_A, h_B, and h_M are the enthalpy of the stream A, stream B, and mixture M respectively.

Eliminating ṁ_M gives the following:

From these balances, it can be shown that the state of the mixture M (T_M, ṁ_M, ω_M) lies on a straight line connecting states A (T_A, ṁ_A, ω_A) and B (T_B, ṁ_B, ω_B) on the psychrometric chart (figures below). The ratio of the distances AM and MB is equal to the ratio of mass flow rate ṁ_A and ṁ_B. First following figure demonstrates mixing process for the case 1. Air stream A, with flow rate of 0.0001 kg/sat 13 °C and relative humidity of 80 %, is mixed with air stream B, with flow rate of 0.00005 kg/sat 36 °C and relative humidity of 20%.

The mass ratio of the two streams of air is 2/1. Under these conditions, the point M, representing the mixture, is located on the line AB and has the following temperature:

From the chart, the relative humidity is, Φ_M, 49.5% and enthalpy is, h_M, 39.7 kJ/kg.

Results

The following table compares temperature and relative humidity values of the mixture, predicted by the thermal flow solver with the theoretical results calculated from Eq. VVC5 - Incompressible laminar heated flow over a cylinder and from psychrometric charts. The computed results are in agreement with the theoretical values.