SimScale CAE Forum

CFD simulation/ dust distribution at pedestrian level


#1

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.

    with trees

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
https://www.simscale.com/workbench/?pid=1301496307486336272&mi=spec%3Adad7914b-b839-4571-9d97-00135ac164fc%2Cservice%3ASIMULATION%2Cstrategy%3A3

thank you for your help,


#2

Hi @Mo_Elsebaei - very nice project idea! :slight_smile:

  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 - for the calculation of the concentration of dust you could use the approach from @tcakir here: Pollution exhaust in a garage

  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

Source for the porosity calculation: https://www.cfd-online.com/Forums/main/13996-porosity-trees.html


#3

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.


#4

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


#5

thank you roy, unfortunately many trail of refinements ended by failing because out of memory error
Capture
do you think i should select a smaller area for the study?


#6

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


#7

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… ?)


#8

Hi @Mo_Elsebaei!

Could you find additional information about the coefficients first before entering them? Here is a doc entry for porous media: https://www.simscale.com/docs/simulation-setup/model/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


#9

hi @Mo_Elsebaei

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.

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

Cheers,
Rohit.


#10

thank you so much :slight_smile: