CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

particle-laden flow simulation using DPMFoam: turbulent dispersion is not heppening

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 8, 2021, 22:00
Default
  #21
New Member
 
Shuo Mi
Join Date: Nov 2020
Location: London
Posts: 19
Rep Power: 5
Msure is on a distinguished road
Quote:
Originally Posted by atul1018 View Post
Hello


-In the original reference experimental setup of Fessler & Eaton, it is clearly written that the setup is held vertically (See 2.1, first paragraph first sentence: https://citeseerx.ist.psu.edu/viewdo...=rep1&type=pdf)


-I also had downloaded the presentation from that seminar but in my slides the gravity is showed to set in +x direction (see attachement).


-Yes, you should compare the particle velocity profiles at measurement locations.


Best Regards
Atul Jaiswal
Hello,

I think you are right about the direction of gravity, then I tried all kinds of boundary conditions, meshes, and dispersion model. They didn't work. The results are the same as yours - particles not right but the fluid velocity is right compared with the experiment. Now I will try another experiment which is particle dynamical movement in a ventilated room to validate the DEM model of OpenFOAM, to see if the phenomenon we got in OpenFOAM in the backward-case is individual one. Thanks.

Regards,
Shuo
Msure is offline   Reply With Quote

Old   January 11, 2021, 05:10
Default
  #22
Senior Member
 
Join Date: Jun 2020
Posts: 100
Rep Power: 6
atul1018 is on a distinguished road
Hello Shuo


So far, i couldn't find any problem in my case but unfortunately particle dispersion is not occurring as expected.


Anyway good luck for new case and keep me updated about your progress.


Best Regards
Atul
atul1018 is offline   Reply With Quote

Old   January 13, 2021, 08:06
Default
  #23
New Member
 
Shuo Mi
Join Date: Nov 2020
Location: London
Posts: 19
Rep Power: 5
Msure is on a distinguished road
Quote:
Originally Posted by atul1018 View Post
Hello Shuo


So far, i couldn't find any problem in my case but unfortunately particle dispersion is not occurring as expected.


Anyway good luck for new case and keep me updated about your progress.


Best Regards
Atul
Hello Atul,

When I tested another case, the results are still not very good. What I compared is the particle concentration in a ventilated room. I doubt I calculated the concentration wrong. I calculated it in this way:
Quote:
template<class CloudType>
inline const Foam::tmp<Foam::volScalarField>
Foam::KinematicCloud<CloudType>::concentration() const
{
tmp<volScalarField> tconcentration
(
new volScalarField
(
IOobject
(
this->name() + ":concentration",
this->db().time().timeName(),
this->db(),
IOobject::NO_READ,
IOobject::NO_WRITE,
false
),
mesh_,
dimensionedScalar(dimless, Zero)
//extrapolatedCalculatedFvPatchScalarField::typeName
)
);


scalarField& concentration = tconcentration.ref().primitiveFieldRef();

for (const parcelType& p : *this)
{
const label celli = p.cell();

concentration[celli] += p.nParticle()*p.mass();
}

concentration /= mesh_.V();
//concentration.primitiveFieldRef() /= 1;//mesh_.V();
//concentration.correctBoundaryConditions();

return tconcentration;
}
and also a similar way to calculate the particle velocity in that backward facing case. So how did you calculate particle velocity in the y-direction in backward case. I want to know if my statistical approach for particles are wrong.
Thanks!

Regard,
Shuo
Msure is offline   Reply With Quote

Old   January 13, 2021, 08:28
Default
  #24
Senior Member
 
Join Date: Jun 2020
Posts: 100
Rep Power: 6
atul1018 is on a distinguished road
Hello


As i don't know the details about your ventilation case, I cant really say about your code regarding concentration calculation.



Regarding BFS case, I extracted particle ID, velocity and position at certain time (say 0.07 sec after insertion of particles). The extracted data are evaluated at slice around each measurement location (x/H=2,5,7,9,12) keeping the measurement location at center of slice. Now you have particle distribution along y direction in each slice. Now you can plot particle velocity profile for the particle in that slice.


I extracted particle data exactly the same way done by Greifzu et.al. You should read the paper of Greifzu et.al.: https://www.tandfonline.com/doi/pdf/...0.2015.1104266



Best Regards
Atul
YANGLU1021 likes this.
atul1018 is offline   Reply With Quote

Old   May 10, 2021, 09:32
Default
  #25
New Member
 
YANG Lu
Join Date: Feb 2021
Location: Usa
Posts: 4
Rep Power: 5
YANGLU1021 is on a distinguished road
Send a message via ICQ to YANGLU1021 Send a message via AIM to YANGLU1021 Send a message via Yahoo to YANGLU1021 Send a message via Skype™ to YANGLU1021
Hi, atul1018. Did you solve the problem? I have met the same problem recently. I also found that the particles in the Y-axis direction do not settle down. If you have found the solution, could you please share your method? Thanks a lot.
YANGLU1021 is offline   Reply With Quote

Old   May 11, 2021, 23:04
Default
  #26
New Member
 
YANG Lu
Join Date: Feb 2021
Location: Usa
Posts: 4
Rep Power: 5
YANGLU1021 is on a distinguished road
Send a message via ICQ to YANGLU1021 Send a message via AIM to YANGLU1021 Send a message via Yahoo to YANGLU1021 Send a message via Skype™ to YANGLU1021
Quote:
Originally Posted by atul1018 View Post
Hello


As i don't know the details about your ventilation case, I cant really say about your code regarding concentration calculation.



Regarding BFS case, I extracted particle ID, velocity and position at certain time (say 0.07 sec after insertion of particles). The extracted data are evaluated at slice around each measurement location (x/H=2,5,7,9,12) keeping the measurement location at center of slice. Now you have particle distribution along y direction in each slice. Now you can plot particle velocity profile for the particle in that slice.


I extracted particle data exactly the same way done by Greifzu et.al. You should read the paper of Greifzu et.al.: https://www.tandfonline.com/doi/pdf/...0.2015.1104266



Best Regards
Atul
Hi, atul1018. Did you solve the problem? I have met the same problem recently. I also found that the particles in the Y-axis direction do not settle down. If you have found the solution, could you please share your method? Thanks a lot.
YANGLU1021 is offline   Reply With Quote

Old   May 12, 2021, 04:41
Default
  #27
Senior Member
 
Join Date: Jun 2020
Posts: 100
Rep Power: 6
atul1018 is on a distinguished road
Hey YANGLU1021
Quote:
Did you solve the problem? I have met the same problem recently. I also found that the particles in the Y-axis direction do not settle down. If you have found the solution, could you please share your method? Thanks a lot.
-I tried with all I could but I could not able to get better results for particles. Let me know if you make any progress.
atul1018 is offline   Reply With Quote

Old   May 13, 2021, 22:20
Default
  #28
New Member
 
YANG Lu
Join Date: Feb 2021
Location: Usa
Posts: 4
Rep Power: 5
YANGLU1021 is on a distinguished road
Send a message via ICQ to YANGLU1021 Send a message via AIM to YANGLU1021 Send a message via Yahoo to YANGLU1021 Send a message via Skype™ to YANGLU1021
Quote:
Originally Posted by atul1018 View Post
Hey YANGLU1021

-I tried with all I could but I could not able to get better results for particles. Let me know if you make any progress.
Thanks for your reply. I will try some ways to see if this situation can be improved. I will tell you if I get any progress. Thanks and have a good day!
YANGLU1021 is offline   Reply With Quote

Reply

Tags
dpmfoam, komegasst, particle-laden flow, rans model

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Will the results of steady state solver and transient solver be same? carye OpenFOAM Running, Solving & CFD 9 December 28, 2019 05:21
Multiphase flow - incorrect velocity on inlet Mike_Tom CFX 6 September 29, 2016 01:27
Turbulent dispersion in Lagrangian particle tracking cvh CFX 0 October 29, 2015 12:19
Beginner level doubts with turbulent flow simulation setup Sagun OpenFOAM 1 October 22, 2012 09:56
Turbulent Dispersion in Particle Tracking ariel CFX 2 April 22, 2008 22:02


All times are GMT -4. The time now is 11:29.