The following picture shows how the workbench looks after importing the project.
Figure 2: Workbench after importing the tutorial project.
In this example, we want to have multiple faces on the ground to assign different terrain roughness. Therefore, we modeled the air domain in a CAD software and split the ground surface.
Did you know?
For simulations involving external flow, it’s important that the flow region is large. This is mostly necessary to avoid interferences from the boundaries. Furthermore, the domain should be longer downstream from the geometry, to capture the wake effects.
Assuming “L” is the longest edge of any building in the geometry, use the following rule of thumb to create a flow domain.
Figure 3: Minimum enclosure size. The enclosure should be longer in the downstream direction.
1.2 Create the Simulation
To set up a new simulation, select the geometry in the left-hand side panel:
Figure 4: Creating a simulation for a geometry.
Hitting the ‘Create simulation’ button leads to the following options:
Figure 5: Analysis types available in SimScale.
Choose ‘Convective heat transfer’ for analysis type and ‘Create the simulation’.
Now you can see the default simulation tree, like presented in the following picture:
Figure 6: Global settings for the convective heat transfer simulation.
In the global simulation set up, please change the Turbulence model to ‘k-omega SST’. This model switches in between k-omega and k-epsilon automatically, therefore it takes advantage of both models.
Also, add 2Passive species. These will help to define the sources of smoke from the chimneys and differentiate the concentration of them and clean air.
2. Set up the Simulation
In this part of the tutorial we will define the physics of the model we want the solver to calculate.
2.1 Define a Material
This simulation will use ‘Air’ as fluid material. Therefore click on the ‘+ button’ next to Materials. In doing so, a library of fluid materials opens:
Figure 7: Fluid materials library.
Select ‘Air’ and click ‘Apply’. In the window that opens, you can leave the default values and assign them to the entire flow region.
2.2 Model
In the next step, we define gravity according to the global coordinate system. Therefore click on Model in the simulation tree. For this geometry, gravity direction is negative z-axis. Hence, input ‘-9.81’ in the z-direction. Leave the Diffusion coefficients and Turbulent Schmidt number as default.
Figure 8: Configuring gravity in the model tab. The direction is based on the orientation cube in the bottom right corner.
2.2 Assign the Boundary Conditions
The following picture summarizes the boundary conditions applied in this simulation:
Figure 9: Overview of the boundary conditions for this tutorial.
a. Walls – Slip
Follow the instructions presented in the picture below to add a new boundary condition:
Figure 10: From the list of boundary conditions, choose ‘wall’.
After hitting the ‘+ button’ next to boundary conditions, a drop-down menu will appear, where one can choose between different boundary conditions;
Select the ‘Wall’ boundary condition.
Figure 11: Assigning a slip wall condition to the top and sides of the enclosure.
Change (U) velocity to ‘Slip’ and (T) Temperature to ‘Adiabatic’. Now assign the top, left and right walls to this condition.
b. Pressure Outlet
Create a new boundary condition, but this time choose ‘Pressure outlet’. The outlet is further away from the chimneys.
Figure 12: Pressure outlet boundary condition.
c. Velocity Inlet – Chimneys
Create yet another boundary condition. Select ‘Velocity inlet’. Let’s first set up the inlet for the tall chimney.
Figure 13: Setting up a velocity inlet with passive scalars.
Smoke will be released at ‘8 m/s’ in the z direction (Uz) and at
A (T) Temperature of ‘450 °C’.
The first passive species is therefore linked to the inlet flow from the tall chimney. That is why we select ‘1’ for the first Value within Passive scalars.
Set up another velocity inlet, this time for the short chimney:
Figure 14: Velocity inlet for the short chimney. This one is linked to the second passive species.
Velocity will be ‘5 m/s’ in the positive z-direction.
The temperature is ‘300 °C’.
This inlet will receive a Value of ‘1’ for the second passive species.
d. Velocity Inlet – ABL
In this step, we will set up an Atmospheric Boundary Layer (ABL). It consists of a pre-defined velocity profile, with value zero on the ground, gradually increasing with height. In this tutorial, \(u_{ref}\) is 5 m/s and \(h_{ref}\) is 10 m.
Figure 15: Atmospheric boundary layer showing the velocity profile at the inlet.
Refer to these equations for the formulation of the logarithmic profile:
This logarithmic profile has the following formulation:
Follow the same procedure as presented in figure 10 but choose a custom boundary condition:
Figure 16: Steps to input a formula for a parameter.
To input a formula for velocity, select a ‘Fixed value’, and then click on the highlighted button to open further options:
Figure 17: In this case, the velocity in the x-direction is a function of height.
Afterward, two tabs will be available: Formula and Table. Switch to the formula input, and paste the velocity formula posted above in the x-direction. Hit ‘Apply’ to save.
Repeat the same process for k and \(\omega\). The other parameters should be as follows:
Figure 18: Fully set up custom boundary condition for the inlet. The settings different from default are highlighted.
Solid walls will receive a no-slip condition. Furthermore, we will specify a surface roughness for the building walls, as well as for the terrains. For more information about the surface roughness model, check the relevant documentation.
Create a Custom boundary condition.
Figure 19: With this configuration, the building walls receive a uniform wall roughness of 0.001 meters.
Building walls will receive a Roughness height of ‘0.001 m’ with a
Roughness constant of ‘0.5’.
To save time, assign this boundary condition to a pre-created topological entity set named ‘Building walls‘:
f. Terrains
The terrains will receive roughness as well. The custom boundary condition is the same as above. The only difference is the Roughness height:
Terrain 1: Roughness height of 0.01 m;
Terrain 2: Roughness height of 0.0075 m;
Lastly, for terrain 3: Roughness height of 0.005 m.
You can also use the pre-saved topological entity sets to correctly assign the terrains.
Figure 20: Setting up the wall roughness for a terrain.
Repeat the same process for terrain 2 and 3.
Note
Check out this page, if you are interested in other boundary conditions available in SimScale.
2.3 Initial Conditions
To enhance convergence, please initialize the velocity field with ‘5 m/s‘ in the x-direction. This velocity represents \(u_{ref}\):
Figure 21: Initializing the global velocity field with the reference wind velocity value.
2.4 Numerics and Simulation Control
The default settings for Numerics are optimized based on the analysis type. Therefore, they work well for most of the simulations and require no changes. However, if you are a simulation expert, you have full control over the settings and may change them as you like.
Access the Simulation control. Here we can define the number of iterations for the simulations as well as the write interval for specifying how many intermediate results to safe.
Figure 22: Adjusting the number of iterations and writing controls for the simulation.
In Simulation control, change End time and Write interval to ‘600’.
Note
For more information about simulation control and numerics, check the corresponding documentation pages.
2.5 Result Control
By setting result controls, you can observe the convergence behavior of several parameters of interest. Hence it is an important indicator to evaluate the quality and trustability of the results.
For this simulation, set up an Area average result control on the outlet. Please proceed as below:
Figure 23: This result control monitors the average levels for all parameters at the outlet.
Expand the Result control tab;
Click on the ‘+ button‘ next to Surface data;
Set Write interval to 5. This should be enough to see the convergence patterns;
Assign it to the outlet face.
3. Mesh
To create the mesh, we recommend using the standard algorithm. This algorithm is a good choice in general, as it is quite automated and delivers good results for most of the geometries.
Figure 24: Main settings for the standard algorithm.
Disable Automatic boundary layers. The layers will be specified manually later. To improve the resolution around the areas of interest, we will set up three refinements.
3.1 Local Element Size Refinement
Click on the ‘+ button’ next to Refinements. The first refinement is a Local element size. This refinement limits the maximum edge length of cells in chosen entities. Hence, it’s useful to improve the resolution in specific regions.
Figure 25: Local element size refinement applied to a topological entity set.
Click on the ‘+ button’ next to Refinements.
Limit the Element size to 0.3 m:
To save time, assign it to a previously created topological entity set named Local element size refinement. This set consists of the solid walls of the buildings and the smoke inlets. you can find it in the geometry panel on the right side of the workbench.
3.2 Inflate Boundary Layer Refinement
Figure 26: Inflate boundary layer refinement set up. Layers should never be added to inlets and outlets.
Create an Inflate boundary layer refinement.
Keep the settings as default and assign it to these 4 topological entity sets:
Terrain 1, 2 and 3
Building Walls
3.3 Region Refinement
To improve the resolution of the region around the chimneys, let’s create a region refinement. Please proceed as below:
Figure 27: Creating a region refinement.
First, create a Region refinement.
Leave 1 meter as Maximum edge length and click on the ‘+ button’ to create a geometry primitive. Choose a Cartesian box and give it the following coordinates:
Figure 28: Cartesian box created for the region refinement. The box should be longer downstream to capture the wake.
Finally, assign this cartesian box to the region refinement. This step finalizes the set up of the simulation.
4. Start the Smoke Propagation Simulation
Click on the ‘+ button‘ next to Simulation runs and ‘Start‘ the process.
Figure 29: Starting a new simulation run.
While it’s running, you can access the intermediate results by clicking on Solution fields. They’re updated as the iterations go by!
The entire smoke propagation simulation takes around 100 minutes to finish. In the reference project, which is linked below, we calculated a few more iterations, to assure complete convergence.
Figure 30: To access the post-processing environment, click on either of these two buttons.
5. Post-Processing
5.1 Result Controls
Once the run is finished for your smoke propagation simulation, open the area average at the outlet result control. For steady-state simulations, the parameters are expected to converge and stabilize as iterations go on. Let’s inspect the average level of the first passive scalar (identified as T1):
Figure 31: Passive scalar level at the outlet, showing good convergence. This is the scalar released by the tall chimney.
For this simulation, other parameters are also useful to assess convergence, such as velocity, temperature, and passive species levels at the outlet.
5.2 Particle Trace
The first post-processor filter that we will use is the ‘Particle Trace’ filter. This filter is very useful to know how the chimneys affect the airflow in the fluid domain. Please follow the steps below to visualize the particle traces:
Make sure that you are at the last timestep of your simulation and that there are no predefined filters in the post-processor.
Figure 32: Ensure that you are at the last timestep of the simulation and that there are no predefined filters.
Click the ‘Add Filter’ button in the Filters panel and choose ‘Particle Trace’. This will take you to the settings of the particle trace.
Figure 33: Click the ‘Add Filter’ button and choose‘Particle Trace’
Next, place the seeds of the streamlines by clicking on the icon next to Pick Position. Place this at the inlet and align it with the height of the taller chimney. If you want to make sure the height of the seeds is appropriate, you can select the top and side enclosure walls, and hide them by right-clicking and choosing ‘Hide selection’.
Figure 34: Click the button beside Pick Position to choose the origin point of the particle traces and align it with the height of the chimney
Adjust the visualization of the streamlines. Increase the number of origin points by changing the #Seeds vertically to ’20’. Reduce the Spacing to ‘0.5’ and Size to ‘0.1’.
Figure 35: Increase the #Seeds vertically to ’20’ and reduce the Spacing to ‘0.5’ and Size to ‘0.1’
For the second chimney, you can repeat the same steps outlined above. As a result, the particle traces will be similar to the figure below:
Figure 36: Particle traces aligned with the chimneys with the flow rising due to gas coming out of of the chimney
Based on Figure 36, the gas coming out of the chimneys causes the flow to rise and creates a higher velocity.
5.3 Cutting Plane
Next, let’s utilize a cutting plane filter to observe how the smoke propagation, or gas coming out of the chimneys, affects the temperature of the surroundings. Create a cutting plane by following the steps below:
Click the ‘Add Filter’ button in the Filters panel and choose ‘Cutting Plane’
Under Position, change the cutting plane placement by typing in the following coordinates: ’50, 0.85, 25′. The Orientation of the cutting plane should be in the ‘Y’ axis. Furthermore, please enable the velocity vectors with the on/off toggle next to Vectors. To improve the visualization, change the vectors’ Coloring to a black ‘Solid color’. Reduce the Scale factor to ‘0.01’ and the Spacing to ‘0.005’.
Figure 37: Cutting plane in the ‘Y’ axis with velocity vectors
From the figure above, we can see that hot gas comes out of the chimney, affecting the air temperature downstream. The buoyancy effects are also noticeable, with the hot exhaust air quickly rising towards the top of the enclosure.
5.4 Iso Volume
The final filter that will be showcased is the Iso volume filter which is very useful to observe the propagation of the gases. Please follow the steps below to set an iso volume filter:
Click the ‘Add Filter’ button in the Filters panel and choose ‘Iso Volume’
Change the Iso scalar and Coloring to ‘Passive scalar 1’ and the Min. iso value range to ‘0.1’ and the Max iso value to ‘1’. With these settings, we are going to highlight the volumes which consist of 10% to 100% of exhaust gases
Finally, change the Opacity to ‘0.5’.
Figure 38: Isovolume in the range of ‘0.1’ to ‘1’ for the gas coming out of the taller chimney
Follow the steps above for the gas coming out of the other chimney and changing the Iso scalar and Coloring to ‘Passive scalar 2’.
Figure 39: Isovolumes of passive scalar coming from both chimneys in the range of ‘0.1’ to ‘1’.
From Figure 39, you can see the smoke propagation results as the concentration of passive scalars is highest at the outlet of the chimney and decreases as it goes into the environment.
Analyze your smoke propagation 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 smoke propagation tutorial!
Yang, Yi, et al. “Consistent Inflow Boundary Conditions for Modelling the Neutral Equilibrium Atmospheric Boundary Layer for the SST k-ω Model.” Wind and Structures, vol. 24, no. 5, 2017, pp. 465–480., DOI:10.12989/was.2017.24.5.465.
Note
If you have questions or suggestions, please reach out either via the forum or contact us directly.
Last updated: January 8th, 2021
Did this article solve your issue?
How can we do better?
We appreciate and value your feedback.
Contents
Cookie Policy
SimScale uses cookies to improve your user experience. By using this website you consent to our cookie policy.
This website uses cookies so that we can provide you with the best user experience possible. Cookie information is stored in your browser and performs functions such as recognising you when you return to our website and helping our team to understand which sections of the website you find most interesting and useful.
Strictly Necessary Cookies
Strictly Necessary Cookie should be enabled at all times so that we can save your preferences for cookie settings.
If you disable this cookie, we will not be able to save your preferences. This means that every time you visit this website you will need to enable or disable cookies again.