This article provides a step-by-step tutorial for the full thermal comfort assessment in a meeting room using CFD simulation, including models for convective heat transfer, radiation heat transfer, and age of air.
This tutorial teaches how to:
The typical SimScale workflow will be followed:
First of all, click the button below. It will copy the tutorial project containing the geometry into your own Workbench.
The following picture demonstrates what should be visible after importing the tutorial project.
You can notice an Open inner region operation performed on the geometry called Meeting Room. As the input CAD model included the room walls, occupants, and furniture, the operation was used to create geometry for the volume of air, which represents the negative of the room and everything inside.
The original geometry without the flow volume extraction is included for comparison, under the name Meeting Room (Original). If you want to learn the details of the operation, you can refer to the documentation page on Flow Volume Extraction.
You can visualize the internal cavities of the extracted volume region by changing the render mode to translucent surfaces. Do this by using the top bar at the viewer.
From now on, make sure you are working on the geometry with the flow volume extraction operation, not the original one presented for reference.
Topological entity sets are groups of faces created at this point to be used in assignments for boundary conditions and other concepts, during the simulation setup. They can be found at the right panel:
A number of sets are already provided in the project, but one set for the window is still missing. The picture below shows how to add one:
In the pop-up dialog that appears, name the set Window and click ‘Create new set‘.
Be sure that you have created all geometry operations before adding topological entity sets. If you do it the other way round, you will loose the sets.
Now we can start with the simulation setup. Follow the steps presented in the picture below to create a new simulation:
The simulation library window appears:
Here you can select the analysis type you need.
Perform the following changes:
Do not forget to click the check mark button at the top to save the changes.
In order to have an overview, the following is a description of the simulation model and conditions:
You can explore the corresponding faces for each condition by clicking the topological entity sets at the right panel.
The above situation has some convection induced by the ventilation system and the window is considered to be closed. You will also find a version without forced ventilation by the ducts and the window considered to be open later in the tutorial.
Select the ‘Model’ tree element to specify the scalar transport properties and gravitational acceleration. For the LMA (local mean age of air) model, the following parameters are used:
To define and assign a material, please click on ‘+’ next to ‘Materials’. Doing so, the SimScale material library will pop up:
Now we need to set up the boundary conditions presented in Figure 9. Note that there are two versions of the setup:
However, we will start with the settings necessary for both scenarios, later on you can decide which scenario you want to simulate.
Internal walls are assigned as a wall boundary condition, and as they are facing the inside of the building, they are considered adiabatic.
Create a wall boundary condition and assign the Internal wall set, have a look at Figure 3 to see where to find the topological entity sets. Set up the parameters as shown in the reference picture:
For the occupants, a wall boundary condition is also used, this time with a heat flux power source to model the metabolic heat generation rate. Create a Wall boundary condition and assign the Occupants set. Set up the parameters as shown in the picture:
For the furniture, a wall boundary condition is also used, with adiabatic thermal behavior. Create a Wall boundary condition and assign the Furniture set. Set up the parameters as shown in the picture:
For the external wall, a wall boundary condition is also used, with a thermal wall model and convection to the exterior. Create a Wall boundary condition and assign the External Wall set. Set up the parameters as shown in the picture:
The following articles provide additional information about thermal wall modelling:
You can add a natural convection setup to simulate a model without air conditioning, where the window is open for natural cooling, or a forced convection setup to simulate the scenario when the window is closed and the air conditioning is on.
A flow velocity inlet boundary condition is created according to the following picture:
Now the setup for the velocity inlet boundary condition will pop up:
Please modify the following:
For the flow outlet boundary condition, follow the same procedure, but select a Pressure outlet. Leave all values as default and assign the Outlet Set (or the outlet surfaces), as shown in the image:
For the window, a wall boundary condition is also used, but this time with a layer wall thermal model, convection to the exterior temperature, and an external radiation source to model sunlight. Create a Wall boundary condition and assign the Window set. Set up the parameters as shown in the image:
For the natural convection boundary condition, select a Natural convection inlet/outlet. Assign the Ambient temperature as the expected temperature surrounding your model. Radiation effects on the inlet and outlet won’t be considered with the ‘Transparent’ radiative behavior applied.
For the window, a natural convection boundary condition is used, similar to the one applied to the inlet and outlet faces but this time define the Additional radiative source as 200 W/m2. This setup takes into consideration an external radiation source to model sunlight and the open boundary conditions allowing the air to circulate through the surface in both directions:
Under advanced concepts, you can define further advanced physical conditions to set up your simulation. For this tutorial, you need to define a passive scalar source.
For the LMA aging model, a Volumetric passive scalar source is used, as shown in the image:
Assign the flow region to the concept and specify a Flux value of 1.
For thermal comfort simulations, these options might be interesting:
For the numeric solver parameters, the only change made is the addition of non-orthogonality correctors. This will improve the convergence for the tetrahedral mesh created by the Standard mesher algorithm. Setup the Number of non-orthogonal correctors as shown in the picture:
Regarding simulation control, we keep the setup to the default values.
Result control items are used to retrieve specific computations from the numerical solver. By using them, we can have a look at specific variables at specific regions by querying the computation and output of our quantities of interest.
In order to measure the average LMA at the outlet face, a Result control item is used. It is created as shown in the picture:
Assign the Outlet set and check that the setup coincides with the picture:
Thermal comfort parameters are also queried in the result control items, under Field calculations. Create the concept as shown in the picture:
Setup the parameters as shown in the picture:
A description of the computed quantities and resulting fields can be found on this SimScale documentation page: Thermal Comfort Parameters.
In the mesh setup, all settings are left as default. You do not need to click Generate, as the mesh will be computed as part of the simulation run.
Now that the simulation setup is complete, a new Simulation Run can be created to perform the computation. In the following picture, the whole tree setup is shown, and the item to create the simulation run is highlighted:
In the pop-up window, press Start to immediately begin the computation run.
The computation takes around one and a half hours to complete. If you can’t wait to see the results, at the end of the article there is a link to the completed version of the project.
We will use the integrated post-processor to visualize the temperature inside the meeting room, evaluate the thermal comfort of the room by visualizing the Percentage Mean Vote (PMV), and the local mean age of the air. Before starting to post-process your results, make sure that there are no predefined filters, you are at the last timestep of your simulation, and you hide the ‘Flow region’ in the MESH dialog.
Visualize the temperature by selecting ‘Temperature’ as the Coloring for your parts. Change the units of temperature to Celcius by clicking the units in the legend and selecting ‘°C’ and change the maximum temperature visualized to ’40 °C’. Hide the walls of the room by selecting them and right-clicking on your mouse to select ‘Hide selection’. The temperature distribution inside the room can be seen in the figure below:
From the figure above, we can see that the bodies emit heat and that the temperature increases as it is further from the inlet. The room is unevenly heated with the floor being slightly cool to neutral as we progress towards the ceiling where it is slightly warm.
Next, we will focus on visualizing the thermal comfort parameters, which are the predicted mean vote (PMV) and predicted percentage dissatisfied (PPD). The thermal comfort parameters will be visualized by using the cutting plane filter. Follow the steps below to create a cutting plane:
The predicted mean vote (PMV) thermal comfort index distribution is displayed in the following picture:
The value of the index is clipped to the recommended range of [-3, 3]. This way, we can visualize the areas that are below and above them, with blue and red colored regions respectively.
You can access the Local Mean Age (LMA) by going to ‘T1’ under Area average 1 in the simulation run, which calculates the average quantities at the outlet. The Local Mean Age (LMA) was computed at 348.97 seconds and can be seen in the figure below:
The distribution of the LMA across the room can be visualized using cutting planes, which is similar to creating the cutting plane for the Predicted Mean Vote (PMV). Add the velocity vectors by sliding the slider besides Vectors and change the Coloring to black Solid color and change the Scale factor to ‘0.5’. The settings for the cutting plane in the y-axis and the final distribution for the LMA can be observed as below:
Blue regions show areas with a low age of air and red regions show areas with high age of air. Red areas have poor ventilation reach which is at the top of the meeting room.
Analyze and explore your results with the SimScale post-processor. Have a look at our post-processing guide to learn how to use the post-processor.
Congratulations! You finished the tutorial!
Last updated: January 8th, 2021
We appreciate and value your feedback.