Fill out the form to download

Required field
Required field
Not a valid email address
Required field
Required field
  • Set up your own cloud-native simulation in minutes.

  • Documentation

    Validation Case: Underfloor Air Distribution

    The underfloor air distribution validation case belongs to fluid dynamics. This test case aims to validate the following parameter:

    • Convective heat transfer

    The validation is performed against experimental data obtained by Zhang and Chen\(^1\), who analyzed the particle transport and distribution in ventilated rooms.

    It is possible to use their results because it was concluded that those particles were passive scalars and, therefore, had no impact on the fluid behavior or the room temperature.

    Geometry

    The geometry consists of a room with an underfloor air distribution (UFAD) ventilation system. The heat sources used are four human simulators and six ceiling lights.

    Cool air is fed to the room through the inlets on the floor, while the humans and the lamps heat the surrounding air, creating convection currents. Figure 1 gives an overview of the geometry and room dimensions.

    underfloor air distribution validation geometry
    Figure 1: Dimensions of the room geometry. Adapted from [Zhang and Chen]\(^1\).

    Although the figure adapted from Zhang and Chen\(^1\) has a scale; the only details provided by the paper were the dimensions of the room. Based on that and other schematic figures provided in the reference paper, it was possible to estimate closely the dimensions of the other elements.

    The geometry used in the present air distribution validation case is shown below:

    underfloor air distribution validation geometry used
    Figure 2: CAD model used in this air distribution validation, based on schematics presented by [Zhang and Chen]\(^1\).

    Table 1 gives further details on the dimensions of each element:

    ElementLength in x-Direction \([m]\)Length in y-Direction \([m]\)Length in z-Direction \([m]\)
    Room4.84.22.4
    Human Simulator0.380.90.38
    Inlets0.250.25
    Exhaust0.250.25
    Lamps1.50.1
    Table 1: Dimensions of room elements, in meters

    Analysis Type and Mesh

    Tool Type: OpenFOAM®

    Analysis Type: Steady state, incompressible convective heat transfer analysis

    Turbulence Model: k-epsilon

    Mesh and Element Types: The mesh used in this validation case was created in SimScale with the standard algorithm. Appropriate refinements on the lamps, human simulators, inlets, and exhaust surfaces were applied. Additionally, boundary layer refinements were applied on all surfaces, to ensure the correct results of the thermal boundary layer.

    Table 2 presents metrics and details of the resulting mesh:

    DomainMesh TypeCellsElement Type
    RoomStandard33406213D tetrahedral/hexahedral
    Table 2: Overview of the mesh generated for the room geometry

    Find below, in Figure 3, the standard mesh generated in SimScale.

    standard mesh for a room used in simscale
    Figure 3: Resulting standard mesh. Local refinements around the outlet and lamps are visible.

    Using a mesh clip to inspect the interior cells, it’s possible to observe the transition from tetrahedral cells to hexahedral cells. The boundary layers are also visible:

    standard mesher hybrid tetrahedral hexahedral cells
    Figure 4: Highlighting the transition between tetrahedral and hexahedral cells around one of the human simulators along with boundary layer refinement

    Simulation Setup for Underfloor Air Distribution

    In the reference paper, the authors provided the power generation for the human simulators as a whole. For the air distribution simulation, the total power was distributed according to each face of the simulators. Additionally, the Boussinesq approximation (incompressible flow) was used for this simulation as no significant temperature differences would be present.

    Material:

    • Viscosity model: Newtonian
    • \((\nu)\) Kinematic viscosity: 1.5295e-5 \(m^2/s\)
    • \((\rho)\) Density: 1.1965 \(kg/m^3\)
    • Thermal expansion coefficient: 0.00343 \(\frac {1} {K}\)
    • \((T_0)\) Reference temperature: 298.15 \(K\)
    • \((Pr_{lam})\) Laminar Prandtl number: 0.713
    • \((Pr_{t})\) Turb. Prandtl number: 0.85
    • Specific heat: 1004 \(\frac {J} {kg.K}\)

    Boundary Conditions:

    The boundary conditions will be defined based on the nomenclature and orientation of Figure 1.

    In the table below, the configuration for both velocity and pressure are given at each of the boundaries:

    ElementBoundary TypeBoundary Condition Description
    InletVelocity inlet0.0472 \(m^3/s\) per inlet, at 293 \(K\)
    ExhaustPressure outletGauge pressure fixed at 0 \(Pa\)
    Human simulators (top face)WallNo-slip condition with a turbulent heat flux. The power heat source is defined as 9.6 \(W\) per face.
    Human simulators (side faces)WallNo-slip condition with a turbulent heat flux. The power heat source is defined as 22.6 \(W\) per face.
    LampsWallNo-slip condition with a turbulent heat flux. The power heat source is defined as 64 \(W\) per lamp.
    Wall: positive x-directionWallNo-slip condition, with a fixed temperature of 297.7 \(K\)
    Wall: negative x-directionWallNo-slip condition, with a fixed temperature of 298 \(K\)
    Ceiling: positive y-directionWallNo-slip condition, with a fixed temperature of 298.7 \(K\)
    Floor: negative y-directionWallNo-slip condition, with a fixed temperature of 297 \(K\)
    Wall: positive z-directionWallNo-slip condition, with a fixed temperature of 298.5 \(K\)
    Wall: negative z-directionWallNo-slip condition, with a fixed temperature of 298.3 \(K\)
    Table 3: Overview of the boundary conditions. Wall functions were used for the wall treatment.

    Numerics

    The following settings are different from the default:

    • Relaxation type: Manual
      • \((P)\) Modified pressure field: 0.7
      • \((U)\) Velocity equation: 0.3
      • \((T)\) Temperature equation: 0.5
      • \((k)\) Turb. kinetic energy equation: 0.7
      • \((\epsilon)\) Dissipation rate equation: 0.7
    • Divergence scheme
      • div(phi,U): Gauss linear upwind limited grad

    Result Comparison

    In the reference study, the authors evaluated seven lines with probe points across the middle of the room \((z = 0)\). Each line contains up to 7 probe points, located at \(y\) = {0.1, 0.3, 0.6 1.1, 1.4, 1.7, 2.2} meters.

    The location of the probe points was obtained from Figure 5, provided in the reference paper:

    location of probe points underfloor velocity and temperature assessment
    Figure 5: Location of the velocity probe points, provided by Zhang and Chen [1].

    Velocity Profiles

    The lines of interest for this validation case were V1, V2, and V3. Find the comparison between the reference study and the SimScale results for velocity in Figure 6:

    velocity profile convective heat transfer validation
    Figure 6: Comparison of velocity profiles over lines V1, V2, and V3, between SimScale and the reference paper

    A comparison of the velocity profiles reveals that the SimScale results follow the same trend as the experimental results provided by the reference\(^1\). Since many dimensions had to be estimated and tested, discrepancies were expected.

    Additionally, the results provided by Zhang and Chen lack symmetry across the ZY plane, which should not be the case given that the geometry is symmetric. These errors in the experimental results add to the discrepancies between SimScale and the reference values. The most important thing is that the profiles follow the same trend.

    Temperature Profiles

    A similar comparison was done for temperatures. The temperatures were normalized using equation 1:

    $$T_{norm} = \frac {T-T_{inlet}} {T_{exhaust}-T_{inlet}}\tag{1}$$

    In SimScale’s results, \(T_{exhaust}\) was found to be 299.15 \(K\). In Figure 7, we compare the simulation results to values from the reference paper.

    temperature profile convective heat transfer validation
    Figure 7: Comparison of temperature profiles over lines V1, V2, and V3, between SimScale and the reference paper.

    A comparison of the temperature profiles reveals that the SimScale results follow the same trend as the experimental results, even more closely than the velocity ones. As explained above, disparities were to be expected.

    Having such close results for the temperature indicates that the difference in the velocity profile comes from the estimated size of the inlet sizes, which impacts velocity and not temperature.

    air distribution buoyancy driven flow temperature density gradients
    Figure 8: Buoyancy-driven flow observed around the occupants

    Last updated: September 4th, 2023

    Contents