Fill out the form to download

Required field
Required field
Not a valid email address
Required field

Documentation

Advanced Tutorial: Smoke Propagation From a Chimney

This article provides a step-by-step tutorial for smoke propagation in open air, making use of the passive scalar transport formulation.

passive scalar analysis of a power plant with streamlines
Figure 1: Smoke propagating off the chimneys of a power plant.

Overview

This smoke propagation simulation tutorial teaches how to:

  • Set up and run a convective heat transfer simulation with passive species;
  • Assign boundary conditions, material, and other properties to the simulation;
  • Set up walls with a given wall roughness;
  • Define an Atmospheric Boundary Layer (ABL) profile for the inlet;
  • Mesh with the SimScale standard meshing algorithm;

We are following the typical SimScale workflow:

  1. Preparing the CAD model for the simulation.
  2. Setting up the simulation.
  3. Creating the mesh.
  4. Run the simulation and analyze the results.

1. Prepare the CAD Model and Select the Analysis Type

First of all click, the button below. It will copy the smoke propagation project containing the geometry into your workbench.

The following picture shows how the workbench looks after importing the project.

smoke propagation initial geometry
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.

enclosure external aerodynamics size
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:

creating a new simulation
Figure 4: Creating a simulation for a geometry.

Hitting the ‘Create simulation’ button leads to the following options:

creating a convective heat transfer simulation
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:

adding passive species to convective heat transfer
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 2 Passive 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:

choosing a material to the parts
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.

gravity set up cfd simulations
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:

boundary conditions pollutant dispersion 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:

creating a boundary condition
Figure 10: From the list of boundary conditions, choose ‘wall’.
  1. After hitting the ‘+ button’ next to boundary conditions, a drop-down menu will appear, where one can choose between different boundary conditions;
  2. Select the ‘Wall’ boundary condition.
slip wall external aerodynamics
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.

setting up a pressure outlet boundary condition wind tunnel
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.

passive scalar transport inlet
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:

velocity inlet with passive species
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.

atmospheric boundary layer set up cfd
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:

$$u = \frac{u^{*}}{K} \cdot ln \left (\frac{z + z_{0}}{z_{0}} \right )$$

Where

  • K‘ is the von Karman constant, equal to 0.41;
  • z0‘ is the aerodynamic roughness length. In this tutorial, we will consider it 0.03 meters;
  • z‘ is the height in which the velocity is calculated;
  • u*‘ is the friction velocity, given by:

$$u^{*} = K \cdot \frac{u_{ref}} {ln \left (\frac{z_{ref} + z_{0}} {z_{0}} \right )}$$

The turbulent kinetic energy (k) formulation is shown below. \(c_\mu\) is the turbulent viscosity constant, equal to 0.09:

$$k = \frac{{u^{*}}^{2}} {\sqrt{c_{\mu}}}$$

Finally, for the specific dissipation rate \(\omega\), the formula is:

$$\omega = \frac{u^{*}}{K \cdot c_{\mu }} \cdot \frac{1}{z + z_{0}}$$


For more information about the setup parameters, please visit this page.

We will use formula inputs for velocity, k and \(\omega\). Based on the equations above:

Velocity = (0.41*5/log((10+0.03)/0.03))/0.41*log((z+0.03)/0.03)

k: ((0.41*5)/(log((10+0.03)/(0.03)))^2/(0.09)^0.5)

\(\omega\): (0.41*5/log((10+0.03)/0.03))/(0.41*(0.09)^0.5)*1/(z+0.03)

Follow the same procedure as presented in figure 10 but choose a custom boundary condition:

formula input for velocity
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:

setting up an atmospheric boundary layer for velocity
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:

custom boundary condition for atmospheric boundary layer
Figure 18: Fully set up custom boundary condition for the inlet. The settings different from default are highlighted.
  • Pressure type: ‘Zero gradient’;
  • (T) temperature: ‘Fixed value’ at ‘19.85 °C‘;
  • \(\alpha_t\) Turb. thermal diffusivity: ‘Zero gradient‘.
  • Passive scalars: ‘Fixed value’ at ‘0’.

e. Building Walls

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.

applying wall roughness in cfd simulations
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.

terrain wall roughness cfd simulation
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}\):

velocity global settings initial condition
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.

simulation control settings
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:

setting result controls to assess convergence
Figure 23: This result control monitors the average levels for all parameters at the outlet.
  1. Expand the Result control tab;
  2. Click on the ‘+ button‘ next to Surface data;
  3. Set Write interval to 5. This should be enough to see the convergence patterns;
  4. 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.

standard mesh external flow set up
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.

local element size refinement in the standard mesher
Figure 25: Local element size refinement applied to a topological entity set.
  1. Click on the ‘+ button’ next to Refinements.
  2. Limit the Element size to 0.3 m:
  3. 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

boundary layer refinement standard algorithm
Figure 26: Inflate boundary layer refinement set up. Layers should never be added to inlets and outlets.
  1. Create an Inflate boundary layer refinement.
  2. 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:

region refinement creation for cfd simulations
Figure 27: Creating a region refinement.
  1. First, create a Region refinement.
  2. 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:
region refinement to capture wake effects
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.

simscale starting a simulation
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.

accessing post processing environment
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):

passive scalar result control
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. For example, velocity, temperature, and passive species levels at the outlet.

5.2 Post-Processing Pictures

The following picture shows a slice of the geometry, created with a cutting plane filter. Velocity vectors are plotted onto the plane, highlighting the atmospheric boundary layer at the inlet.

atmospheric boundary layer generated at the inlet
Figure 32: Highlighting the velocity profile generated at the inlet.

The Isovolumes filter is particularly helpful to show how the passive scalars are spreading in the environment. In the figure below, we are showing the regions where T2 is between 0.01 and 1.

passive scalar transport showing smoke propagation
Figure 33: Using the isovolumes filter for T2 (second passive species). It’s possible to see how it is spreading through 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!

References

[1] 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: October 21st, 2020

Contents
Data Privacy