Fan and Grille Boundary Conditions
Fans, pumps, or grilles (devices) can be represented using lumped curves implemented as boundary conditions. These simplifications also imply some assumptions. In particular, it is assumed that a given boundary can only be either an inlet or an outlet. Such a boundary should not be a mix of inlets/outlets, nor switch between them during a simulation.
Manufacturers usually provide curves that describe the static pressure as a function of flow rate for a fan. See Fan and Grille for the node settings.
Defining a Device at an Inlet
In this case, the device’s inlet is an external boundary, represented by the external circular boundary of the green domain on Figure 3-4. The device’s outlet is an interior face situated between the green and blue domains in Figure 3-4. The lumped curve gives the flow rate as a function of the pressure difference between the external boundary and the interior face. This boundary condition implementation follows the Pressure Boundary Condition for inlets with the Suppress backflow option:
(3-41)
The Grille boundary condition sets the following conditions:
(3-42)
Here, V0 is the flow rate across the boundary, pinput is the pressure at the device’s inlet, and Δpfan(V0) and Δpgrille(V0) are the static pressure functions of flow rate for the fan and the grille. Equation 3-41 and Equation 3-42 correspond to the compressible formulation. For incompressible flows, the term −(2/3)μ(∇ ⋅ u) vanishes. When a turbulence model with a transport equation for the turbulent kinetic energy is applied, the term −(2/3k appears on the left sides of Equation 3-41 and Equation 3-42. In such cases the turbulent kinetic energy, k, the turbulent relative fluctuations, ζ (for the v2-f turbulence model), and dissipation rate, ε, or specific dissipation rate, ω, must be specified on the downstream side. The turbulence conditions are specific to the design and operating conditions of the fan or grille. For the Fan condition, a reference velocity scale Uref is available in order to set default values according to Equation 3-40. For the Grille boundary condition the turbulence quantities on the downstream side are specified by defining a loss coefficient K = 2Δpgrille/ρ(u ⋅ n)2, from which a refraction coefficient is derived using Equation 3-55. Equation 3-48 through Equation 3-50 are then used to relate upstream and downstream turbulence quantities.
When an algebraic turbulence model, such as the Algebraic yPlus or L-VEL model, is used, the turbulence is only implicitly affected by the device through the change in the local Reynolds number.
The flow direction is enforced to be normal to the boundary by default. However, it is possible to define an arbitrary inflow velocity direction. In such case, Equation 3-41 and Equation 3-42 are modified similarly to Equation 3-36. In 2D, the thickness in the third direction, Dz, is used to define the flow rate. Fans are modeled as rectangles in this case.
For 3D and 2D axisymmetry when the swirl flow is included, it is possible to define a swirl on the downstream side of the fan. It sets the following boundary conditions:
Here, the swirl ratio csf is a positive number less than 1, defining the ratio of the rotation transferred from the fan to the flow, f is the number of revolutions per time for the fan, and rbp is the rotation axis base point. The swirl velocity is defined on the maximum inscribed circular region on the fan inlet boundary. Hence, it is only applicable on simply connected flat boundaries. For more complicated inlet topologies, a swirl velocity can be defined by specifying the inlet velocity vector.
Figure 3-4: A device at the inlet. The arrow represents the flow direction, the green circle represents the device (that should not be part of the model), and the blue cube represents the modeled domain with an inlet boundary condition described by a lumped curve for the attached device.
Defining a Device at an Outlet
In this case (see Figure 3-14), the fan’s inlet is the interior face situated between the blue (cube) and green (circle) domain while its outlet is an external boundary, here the circular boundary of the green domain. The lumped curve gives the flow rate as a function of the pressure difference between the interior face and the external boundary. This boundary condition implementation follows the Pressure Boundary Condition for outlets with the Suppress backflow option:
(3-43)
The Grille boundary condition sets the following conditions:
(3-44)
Here, V0 is the flow rate across the boundary, pexit is the pressure at the device outlet, and Δpfan(V0) and Δpgrille(V0) are the static pressure function of flow rate for the fan and the grille. Equation 3-43 and Equation 3-44 correspond to the compressible formulation. For incompressible flows, the term −(2/3)μ(∇ ⋅ u)n vanishes. When a turbulence model with a transport equation for the turbulence kinetic energy is applied, the term −(2/3kn appears on the left sides of Equation 3-43 and Equation 3-44. In 2D the thickness in the third direction, Dz, is used to define the flow rate. Fans are modeled as rectangles in this case.
Figure 3-5: A fan at the outlet. The arrow represents the flow direction, the green circle represents the fan (that should not be part of the model), and the blue cube represents the modeled domain with an outlet boundary condition described by a lumped curve for the attached fan.