# CFD simulation/ dust distribution at pedestrian level

**URL:** https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094
**Category:** Project Support
**Created:** [August 12, 2019, 7:57pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094 "2019-08-12T19:57:37Z")
**Posts on this page:** 17
**Page:** 1

<div class="post-metadata">

### Author: ![Mo\_Elsebaei](https://www.simscale.com/forum/letter_avatar/mo_elsebaei/32/5_5575768a8748004e209b776fc1b2916d.png) [@Mo\_Elsebaei](https://www.simscale.com/forum/u/Mo_Elsebaei)
#### Post date: [August 12, 2019, 7:57pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/1 "2019-08-12T19:57:37Z")

</div>

Hello everyone,

- Overview about the project:  
my project is to investigate the dust flow and concentration at different points in outdoor space of small community before and after adding different types of trees as buffer zone around the community.  
the dust is presented as PM10 which is any particles less than or equal to 10 micrometers in diameter , these particles are one of the major pollutants in the air and have many bad effects on human health.  
first i want to clarify that PM10 are solids but once they are lifted by the air they are converted from a static solid-like state to a dynamic fluid-like state under the phenomena of fluidization.  
the study area surrounded by desert from three sides, which i will consider as the source of the dust in my simulation.  
 ![edited](https://www.simscale.com/forum/uploads/default/original/3X/a/4/a456abc2dfe37e86e049bac8de09f510e674b756.jpeg)  
with trees  
 ![2](https://www.simscale.com/forum/uploads/default/original/3X/7/f/7fc051c6bb4bf372d30f4a06e1d040586d8f5960.jpeg)

i will simulate the windest day in three months before and after adding three types of plants where i will define the wind speed and direction in that day as will as the dust flux value as measured in the field in microgram/m2/s.

my questions are:  
1- is incompressible analysis with k-omega SST Turbulence model/ Steady state / 1 Passive species, the right simulation?  
2- what is the right way to define dust as passive scalar and the 3 surrounding area as the source of dust  
3- how can i define the trees as porous media  
4- is the program capable to present the dust concentration values and distribution in the results ?

i did a test simulation but all the results were showing only wind  
flow

the link for the project

> **[SimScale Login](https://www.simscale.com/workbench/?pid=1301496307486336272&mi=spec%3Adad7914b-b839-4571-9d97-00135ac164fc%2Cservice%3ASIMULATION%2Cstrategy%3A3)**
>
> Log in to your SimScale account. SimScale makes high-fidelity engineering simulation technically and economically accessible. At any scale. In the cloud.

thank you for your help,

---

<div class="post-metadata">

### Author: ![jousefm](https://www.simscale.com/forum/user_avatar/www.simscale.com/jousefm/32/34908_2.png) [@jousefm](https://www.simscale.com/forum/u/jousefm)
#### Post date: [August 12, 2019, 9:22pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/2 "2019-08-12T21:22:54Z")

</div>

Hi @Mo_Elsebaei - very nice project idea! 🙂

1. Sounds good but it depends if you want to see transient structural behavior on the buildings or is steady-steate fine for you? Did you already decide? For the turbulence model maybe @1318980 or @ananthu_ajit can comment as I have not tested k-omega SST myself, all I know is that k-epsilon underpredicts in the strong wind and wake regions, maybe they can comment on that.

2. Maybe this approach here helps you: **[Passive Scalar Sources](https://www.simscale.com/docs/content/simulation/model/advancedConcepts/passivescalarSources/passivescalarSources.html)** - for the calculation of the concentration of dust you could use the approach from @tcakir here: **[Pollution exhaust in a garage - #4 by tcakir](https://www.simscale.com/forum/t/pollution-exhaust-in-a-garage/85652/4)**

3. You can define the porosity ( P ) of a (coniferous) tree by

P = 1-\frac{Dc}{4\*(LAI/h)}

where

- LAI = Leaf Area Index (LAI)
- Dc = Needle diameter
- h = Height of the canopy

* * *

1. Not if you want to use the approach mentioned in 2) I would say as you need the calculator filter which is not yet implemented (but on the roadmap)

Hope I could help you a bit!

Jousef

* * *

**Relevant articles:** Pekin, B.; Macfarlane, C. Measurement of Crown Cover and Leaf Area Index Using Digital Cover Photography and Its Application to Remote Sensing. _Remote Sens._ **2009** , _1_ , 1298-1320.doi:[10.3390/rs1041298](http://dx.doi.org/10.3390/rs1041298)

**Source for the porosity calculation:** [Porosity of Trees -- CFD Online Discussion Forums](https://www.cfd-online.com/Forums/main/13996-porosity-trees.html)

---

<div class="post-metadata">

### Author: ![roy\_g](https://www.simscale.com/forum/user_avatar/www.simscale.com/roy_g/32/25462_2.png) [@roy\_g](https://www.simscale.com/forum/u/roy_g)
#### Post date: [August 13, 2019, 10:42am UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/3 "2019-08-13T10:42:56Z")

</div>

I’ve never done anything like this so not much input other than pay attention to your mesh, especially those trees. Use mesh refinements otherwise you will end up with a massive number of elements. Those tiny features (the trees) are not ideal.

---

<div class="post-metadata">

### Author: ![Mo\_Elsebaei](https://www.simscale.com/forum/letter_avatar/mo_elsebaei/32/5_5575768a8748004e209b776fc1b2916d.png) [@Mo\_Elsebaei](https://www.simscale.com/forum/u/Mo_Elsebaei)
#### Post date: [August 13, 2019, 11:37am UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/4 "2019-08-13T11:37:41Z")

</div>

hi @jousefm, thank you for your fast answer!  
i’ve been reading and looking at all the possible related topics and projects but i couldn’t make even any general idea of how to start my project… it may seems easy to you but for me unfortunately i don’t find it user friendly program… so i will highly highly appreciate more detailed or step by step answer,  
1- my concern is only the outdoor pedestrian comfort level, so i only need to see how the air circulates dust between the building and on the streets. that being said, is transient behavior will give me more accurate result? i sow a project that used Laminar/ Transient, but i really don’t no what the differences.  
2- so i have to define areas surrounding the study area which are the sources of the dust as geometry primitives and then select them as passive scalar source and define the flux, when i did so in a test simulation the highest value of T was 1 where the flux was set to 50/s, so how this is possible? i think i missed something  
3- can you please explain how can i define porosity? and which model support this? (Darcy forchheimer medium, Fixed coeff… ?)

thank you

---

<div class="post-metadata">

### Author: ![Mo\_Elsebaei](https://www.simscale.com/forum/letter_avatar/mo_elsebaei/32/5_5575768a8748004e209b776fc1b2916d.png) [@Mo\_Elsebaei](https://www.simscale.com/forum/u/Mo_Elsebaei)
#### Post date: [August 13, 2019, 12:25pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/5 "2019-08-13T12:25:37Z")

</div>

thank you roy, unfortunately many trail of refinements ended by failing because out of memory error

 ![Capture](https://www.simscale.com/forum/uploads/default/original/3X/2/e/2e130c2b8350d5884a6b6e3ffa6fcdbdafd9e5e1.jpeg)  
do you think i should select a smaller area for the study?

---

<div class="post-metadata">

### Author: ![roy\_g](https://www.simscale.com/forum/user_avatar/www.simscale.com/roy_g/32/25462_2.png) [@roy\_g](https://www.simscale.com/forum/u/roy_g)
#### Post date: [August 13, 2019, 1:24pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/6 "2019-08-13T13:24:38Z")

</div>

> [@Mo\_Elsebaei](#):
>
> refinements ended by failing because out of memory er

You could do, and you can simplify your model so it has less small features especially small rounded features

---

<div class="post-metadata">

### Author: ![Mo\_Elsebaei](https://www.simscale.com/forum/letter_avatar/mo_elsebaei/32/5_5575768a8748004e209b776fc1b2916d.png) [@Mo\_Elsebaei](https://www.simscale.com/forum/u/Mo_Elsebaei)
#### Post date: [August 16, 2019, 2:00pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/7 "2019-08-16T14:00:16Z")

</div>

> [@jousefm](#):
>
> 1. You can define the porosity ( P ) of a (coniferous) tree by
> 
> P=1−Dc4∗(LAI/h)

hello @jousefm, can you please tell me after i calculate the porosity , where to put the value in the simulation? (Darcy forchheimer medium, Fixed coeff… ?)

---

<div class="post-metadata">

### Author: ![jousefm](https://www.simscale.com/forum/user_avatar/www.simscale.com/jousefm/32/34908_2.png) [@jousefm](https://www.simscale.com/forum/u/jousefm)
#### Post date: [August 16, 2019, 4:41pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/8 "2019-08-16T16:41:57Z")

</div>

Hi @Mo_Elsebaei!

Could you find additional information about the coefficients first before entering them? Here is a doc entry for porous media: **[Porosity & Porous Media | Advanced Concepts | SimScale](https://www.simscale.com/docs/simulation-setup/advanced-concepts/porous-media/)**.

@varsey, could you help our user here to clarify if anything is not clear? The use of the porous media module is straightforward but we would have to know the coefficients as previously mentioned. I am not sure what you could find in the literature.

Cheers,

Jousef

---

<div class="post-metadata">

### Author: ![ROHIT\_SR](https://www.simscale.com/forum/user_avatar/www.simscale.com/rohit_sr/32/25850_2.png) [@ROHIT\_SR](https://www.simscale.com/forum/u/ROHIT_SR)
#### Post date: [August 17, 2019, 5:38am UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/9 "2019-08-17T05:38:24Z")

</div>

hi @Mo_Elsebaei

> [@Mo\_Elsebaei](#):
>
> 3- how can i define the trees as porous media

if you want to apply porous media for the tree, first decided porosity for the tree (in term of percentage).

when velocity flow through the tree it will change its nature and become nonlinear and if we consider the steady-state condition than you can only use Darcy-Forchheimer equation for this case.

you can apply porosity and particles diameter only in Darcy-Forchheimer equation (pressure-flux equation) and calculated the value of Darcy coefficients and Forchheimer coefficients.

 ![25%20AM](https://www.simscale.com/forum/uploads/default/original/3X/2/d/2dbec4e2461607ce477009c5f0455c5262fa34b7.jpeg)

follow above formula and used your value for calculation.  
if you still face any problem contact me at any time.😉

Cheers,  
Rohit.

---

<div class="post-metadata">

### Author: ![Mo\_Elsebaei](https://www.simscale.com/forum/letter_avatar/mo_elsebaei/32/5_5575768a8748004e209b776fc1b2916d.png) [@Mo\_Elsebaei](https://www.simscale.com/forum/u/Mo_Elsebaei)
#### Post date: [August 17, 2019, 2:35pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/10 "2019-08-17T14:35:10Z")

</div>

thank you so much 🙂

---

<div class="post-metadata">

### Author: ![Mo\_Elsebaei](https://www.simscale.com/forum/letter_avatar/mo_elsebaei/32/5_5575768a8748004e209b776fc1b2916d.png) [@Mo\_Elsebaei](https://www.simscale.com/forum/u/Mo_Elsebaei)
#### Post date: [August 30, 2019, 11:11am UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/11 "2019-08-30T11:11:50Z")

</div>

@ROHIT_SR , Good morning! i’ve been looking in the literature to find an easy way to calculate the porosity of a tree because the information needed to solve this (P=1−Dc4∗(LAI/h)) is not available for the plants that i will use in the simulation (not for most of the plants in fact!)  
so i found an article that have used a program “SideLook” to calculate the porosity, the program depends on high resolution digital photo of the tree and converts it to black and white image by defining the blue spectrum then it defines the characteristic based on the ratio of white and black pixels. (Velarde, Jacinto Garrido et al. “Estimation of Optical Porosity or Canopy Structure of Two Species of Tree with Hemispherical and Vertical Images.” (2018))  
and the below describes in details the process of defining the vegetation characteristic using the program ([https://reader.elsevier.com/reader/sd/pii/S0367253004700614?token=9F37A8F421208239FE7C3BE8DC5B3A5B4F5D32DB890515AA568B207B567E74F1A9BE87C9A28BA7BF7CFF0CE80FD6BFF9](https://reader.elsevier.com/reader/sd/pii/S0367253004700614?token=9F37A8F421208239FE7C3BE8DC5B3A5B4F5D32DB890515AA568B207B567E74F1A9BE87C9A28BA7BF7CFF0CE80FD6BFF9)).

i wanted to share with you this information in case you are interested in this topic and i also want to ask you what we mean by particles diameter (for Darcy-Forchheimer equation) in the case of a tree? and after i calculate both of Darcy coefficients and Forchheimer coefficients, are they the same for x, y, z or it depends on the porosity of the material in each direction? and what unit vector is used for?

 ![Unit%20vector](https://www.simscale.com/forum/uploads/default/original/3X/8/1/81b1fa571711eb46e7b8fcbbd917510109822fbb.jpeg)

Thank you, i highly appreciate your help

---

<div class="post-metadata">

### Author: ![ROHIT\_SR](https://www.simscale.com/forum/user_avatar/www.simscale.com/rohit_sr/32/25850_2.png) [@ROHIT\_SR](https://www.simscale.com/forum/u/ROHIT_SR)
#### Post date: [August 30, 2019, 11:49am UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/12 "2019-08-30T11:49:32Z")

</div>

hi @Mo_Elsebaei,

> [@Mo\_Elsebaei](#):
>
> i wanted to share with you this information in case you are interested in this topic

Always interested …!  
happy to see your progress because you done lots of research and trying to understand from basic.👏

> [@Mo\_Elsebaei](#):
>
> what we mean by particles diameter

Particles dia. means its spacing between two solid pieces. it was different from material to material. In your case ( Tree) I think the gap between two branches but it was in a circle every where. so I would suggest you, don’t apply porous media for the tree( its too complicated for this simulation.)

> [@Mo\_Elsebaei](#):
>
> both of Darcy coefficients and Forchheimer coefficients, are they the same for x, y, z

yes. it was same for all direction.  
NOTE:- it was also dependent on your B.C. beacuse Darcy coefficients & Forchheimer coefficients was different.  
[read my post for more details](https://www.simscale.com/forum/t/porous-simulation-project-by-rohit-sr/86149/2)

> [@Mo\_Elsebaei](#):
>
> and what unit vector is used for?

Unit vector applied in principle direction ( fluid flow direction). used your geometry x, y and z value for e1 and e2.  
if you have any other doubt, contact me at any time. 📧

Best,  
Rohit.

---

<div class="post-metadata">

### Author: ![Mo\_Elsebaei](https://www.simscale.com/forum/letter_avatar/mo_elsebaei/32/5_5575768a8748004e209b776fc1b2916d.png) [@Mo\_Elsebaei](https://www.simscale.com/forum/u/Mo_Elsebaei)
#### Post date: [August 30, 2019, 12:20pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/13 "2019-08-30T12:20:05Z")

</div>

thank you so much for your support!! 😃

> [@ROHIT\_SR](#):
>
> Particles dia. means its spacing between two solid pieces. it was different from material to material. In your case ( Tree) I think the gap between two branches but it was in a circle every where. so I would suggest you, don’t apply porous media for the tree( its too complicated for this simulation.)

so it is actually refers to the void… well, SideLook also can estimate the average diameter of the holes within the tree with reference to a distance the user defines at the image. i will use this value as particle diameter. i know that the result will not be accurate 100% but at least within an acceptable range 🙂  
thank you so much again

---

<div class="post-metadata">

### Author: ![ROHIT\_SR](https://www.simscale.com/forum/user_avatar/www.simscale.com/rohit_sr/32/25850_2.png) [@ROHIT\_SR](https://www.simscale.com/forum/u/ROHIT_SR)
#### Post date: [August 30, 2019, 1:23pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/14 "2019-08-30T13:23:26Z")

</div>

> [@Mo\_Elsebaei](#):
>
> SideLook also can estimate the average diameter of the holes within the tree with reference to a distance the user defines at the image.

interesting to see your post-processing image. 👍

> [@Mo\_Elsebaei](#):
>
> i know that the result will not be accurate 100%

yes, you are right.

best,  
Rohit.

---

<div class="post-metadata">

### Author: ![mkayed](https://www.simscale.com/forum/letter_avatar/mkayed/32/5_5575768a8748004e209b776fc1b2916d.png) [@mkayed](https://www.simscale.com/forum/u/mkayed)
#### Post date: [September 14, 2019, 4:43pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/15 "2019-09-14T16:43:40Z")

</div>

Hello!  
i have been looking for the formula that you have used to calculate the inertial resistance forchheimer coefficient in the literature ; there is several formulas but none of them have “3.5” , can you please give me the reference of this formula?

Best regards,

---

<div class="post-metadata">

### Author: ![ROHIT\_SR](https://www.simscale.com/forum/user_avatar/www.simscale.com/rohit_sr/32/25850_2.png) [@ROHIT\_SR](https://www.simscale.com/forum/u/ROHIT_SR)
#### Post date: [September 14, 2019, 5:57pm UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/16 "2019-09-14T17:57:34Z")

</div>

hi @mkayed

> [@mkayed](#):
>
> there is several formulas but none of them have “3.5” ,

yes, you are right. no one mention 3.5.  
but,  
if you know the equation of Ergun than you can understand all term involve in inertial resistance (forchheimer coefficient). in Ergun equation 1.7 was mention & if you solve Ergun equation than you get 3.5 at the end.

follow my post for more information:

> **[Porous Media Test by ROHIT\_SR | SimScale](https://www.simscale.com/projects/ROHIT_SR/porous/)**
>
> The aim of the project was to implement a porous media option and try to understand the basic concept behind through CFD setup. for this case, incompressible flow simulation type was used. In this case, the pressure drop across the two was calculated...

reference,

> **[Ergun equation](https://en.wikipedia.org/wiki/Ergun_equation)**
>
> The Ergun equation, derived by the Turkish chemical engineer Sabri Ergun in 1952, expresses the friction factor in a packed column as a function of the modified Reynolds number.
> where 
>   
>     
>       
>         
> f
>           
> p
>           
>         
>       
>     
> {\\displaystyle f\_{p}}
>   
> and 
>   
>     
>       
> G
>         
> r
>           
> p
>           
>         
>       
>     
> {\\displaystyle Gr\_{p}}
>   
> are defined as

if you face any problem contact me any time. 😉

cheers!  
Rohit.

---

<div class="post-metadata">

### Author: ![mkayed](https://www.simscale.com/forum/letter_avatar/mkayed/32/5_5575768a8748004e209b776fc1b2916d.png) [@mkayed](https://www.simscale.com/forum/u/mkayed)
#### Post date: [September 15, 2019, 3:58am UTC](https://www.simscale.com/forum/t/cfd-simulation-dust-distribution-at-pedestrian-level/86094/17 "2019-09-15T03:58:11Z")

</div>

thanks!
