|
[Sponsors] |
July 7, 2022, 12:01 |
Particles captured on walls
|
#1 |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Good morning everyone. First of all, thanks in advance to whoever helps me.
I've been struggling with a model for a while now. To be objective, I'm trying to get the collision rate of 1micron diameter solid particles into a filter geometry on a simple tube. SST turbulence, turbulent dispersion force, fully coupled. To do so, I have set the filter as a Wall with restitution coefficients of 0. From the inlet, 1000 particles run with zero slip velocity. When I run it as steady state, I get the number of captured particles on wall and that's it. However, my simulation never gave me a decent convergence (monitor wouldn't stabilize, residuals wouldn't go below 10^-3). So I studied a lot and improved it using many tips from this forum (thanks for that). I gave up running as steady state and after running it as transient, adaptive, 3-5 coef. loops, I got the timestep of 1.5e-6; courant of 0 and max courant of 0.06; residuals below 10^-6; flat monitor graph. Ran again as transient using this timestep and everything is beautiful. (steady state reached, I believe?) But when I run it as transient, I can't get the output I have been working with on steady state: "Entered domain: 1000; Captured on walls: xxxx; Left domain: xxxx". First of all, I only get 1 particle per iteration entering the domain. Secondly, if I understood correctly, for the 8s residence time I have in this model, for 1.5e-6 timestep, I'd need to run the transient model for 5 millions accumulated timesteps for it to end? Or should I just get the 1.5e-6 timestep and put it into the steady state model? Because when I do that, I dont get the same graphs - my residuals won't go as low, and my monitor won't stabilize. I apologize if that's too big a post, but I tried to give all the info that would help understand my situation. Best regards! Last edited by lgtmelo; July 7, 2022 at 12:13. Reason: including picture of enclosure |
|
July 7, 2022, 17:56 |
|
#2 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
Thanks for the detailed post. It is very useful when people post what they are doing and other details as then we can understand what you are trying to achieve.
Some comments: Steady state versus transient - It is quite likely that the flow through this assembly is transient as you will get the flow moving about behind the bluff filter particles. So I would expect that a transient simulation is required. But I would not expect a transient simulation to settle down into a steady state solution. Are you sure you have run the transient solution long enough to see what happens? Slow particle rate: This is set by the mass flow rate of particles and particle diameter. There are also some particle modelling parameters which affect it as well. So you can increase the number of particles by adjusting these parameters. Fully coupled versus One way coupling: Are you sure you need this to be fully coupled? If you can run it with one way coupling it will be MUCH easier and a hundred times less CPU time.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
July 7, 2022, 18:10 |
|
#3 | ||
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
thank you for answering Glenn. in my learning curve here, I think your posts from 10+ years ago were already very helpful.
Quote:
Quote:
i am definitely NOT sure I need to be fully coupled. I just read that it is a more reliable model, since each fluid will "act" into the other, instead of just behaving as if they were independent. |
|||
July 7, 2022, 18:22 |
|
#4 | |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
Quote:
Your chart shows it decaying to a flat line. But it is just velocity at a point. The transient stuff is likely to be just behind a bit of your filter. Right in the middle of the filter actually. So you are probably not looking at a transient bit. Hold off on the transient particles per second for now, I am not convinced you need transient yet. Yes, fully coupled is a more "reliable" model in that it includes more physics. But the price of reliability is a slower simulation, and in my experience transient particle tracking models can be 100 times slower to run than one way coupled. This obviously depends on exactly what you are modelling, but the point is fully coupled simulations are MUCH slower. So only use fully coupled if you really need it. In your case, I understand you have 1um particles (meaning they are really small), they appear to be bacteria (meaning the density of them is likely to be about the same as the water fluid) and I do not see anything which generates a big slip velocity between the bacteria and fluid - this all suggests a one way coupled model will be fine.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
||
July 7, 2022, 18:26 |
|
#5 | |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Quote:
thanks! ps: i think the oldest post from you i was able to find was from around 2005 haha. and YES, all of them were very helpful to me. Last edited by lgtmelo; July 7, 2022 at 18:36. Reason: question already answered :D |
||
July 7, 2022, 18:43 |
|
#6 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
Transient or steady state is another issue entirely. I would:
* Have a look at your transient final results and see if you can spot where it is transient. * Run it steady state and include the equation residuals in the results file. Then you can look in the post-processor and see where the high residuals are and this is the problem area. This will help decide if it is transient or steady state. Also, I would just run it SS and see if it is different. If no difference then SS should be fine. Finally - note that turbulent dispersion on particle tracking with coupled particle tracking means a steady state result cannot exist (look this up in the doco). So your lack of convergence might be just the turbulent dispersion model. If this is the case then this will be fixed by going to one way coupled particles (if this assumption is valid).
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
July 7, 2022, 18:46 |
|
#7 | |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Quote:
|
||
July 7, 2022, 19:00 |
|
#8 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
You can still have turbulent dispersion in a one way coupled model. It is just that the random jiggling of the particles then does not affect the flow field, so the flow will converge fine.
Also consider whether you actually have turbulent dispersion occurring. The flow in many filters is laminar, so adding turbulent dispersion is not appropriate. Finally - be aware that the particle tracking model does not have complete physics of whether particles impact the surface or not. For instance the particles are modelled as points and it is assumed they hit the walls when that point passes through a wall. But the actual particle has a diameter, so will hit the wall before its centroid point does. This is not taken into account. There are other more subtle effects not taken into account as well.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
July 7, 2022, 19:14 |
|
#9 |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
disregard the next sentence: thats weird, now solver wont show the "particles source change rate" graph, or mention particles entering/exiting domain or captured on walls. is this because of the one way coupling or did I mess something up?
"One-way coupled particles are solved at the end of the simulation" got it! Last edited by lgtmelo; July 7, 2022 at 19:19. Reason: "One-way coupled particles are solved at the end of the simulation" |
|
July 7, 2022, 19:16 |
|
#10 | |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Quote:
"One-way coupled particles are solved at the end of the simulation" got it! Last edited by lgtmelo; July 7, 2022 at 19:19. Reason: "One-way coupled particles are solved at the end of the simulation" |
||
July 7, 2022, 19:18 |
|
#11 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
One way coupling does not need to converge the flow and particle coupling (which is what particle source change rate is). It just iterates the fluids then solves the particles once only. No convergence required.
It will show the particles entering/exiting the domain, just in a different part of the output file. Look in the end section.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
July 7, 2022, 19:21 |
|
#12 | |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Quote:
thanks |
||
July 7, 2022, 19:26 |
|
#13 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
Turbulence model and particle tracking are totally different things. You choose a turbulent model which accurately models the flow, and then the particles do their thing in the flow. There is no direct interaction between the turbulence model and the particles (except when you activate the turbulent dispersion model, and even then the interaction is very diffuse).
I repeat my previous comment: Are you sure this flow is turbulent? Most filters are laminar.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
July 7, 2022, 19:31 |
|
#14 | |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Quote:
I ran it as laminar initially, but solver suggested it was turbulent. also the reynolds it pointed was >5k. |
||
July 7, 2022, 19:38 |
|
#15 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
The solver message about turbulence model and Reynolds Numbers are guides only. You should calculate the Reynolds number of the flow yourself - in your case it looks like a tube, so the Reynolds Number from the tube diameter and liquid properties and velocity. Then compare this to known data (eg the Moody Diagram) to work out whether it is laminar or turbulent. And then you select laminar or turbulent based on that.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
July 7, 2022, 19:40 |
|
#16 | |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Quote:
pic of the infill: |
||
July 7, 2022, 19:53 |
|
#17 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
Then do a quick calc to see if the filter flow is turbulent: You say the diameter is ~3mm, and look at your simulations to see a typical velocity in them (or estimate the velocity if you can). Then you can get a representative Reynolds Number for the flow in the filter.
(But I know the answer already from experience: If the main flow is laminar then the flow in this infill is going to be laminar as well. This means you can turn off the turbulence model and just model it as laminar, and the turbulent dispersion force option is not available. This will greatly simplify your model.)
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
July 7, 2022, 19:55 |
|
#18 | |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
Quote:
|
||
July 7, 2022, 20:42 |
|
#19 |
Member
Join Date: May 2022
Posts: 61
Rep Power: 3 |
quick update while its still running: so far it seems weird. i am considering trying a higher timestep (~7.5s = residence time) . currently on auto, it is 7.76235E-02.
pics: |
|
July 7, 2022, 20:54 |
|
#20 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,706
Rep Power: 143 |
The turbulence model puts extra dissipation in the model which stabilises it. The k-e model will have lots of additional dissipation, the SST model not so much, but still some. When you remove that dissipation the simulation is more accurate, but more likely to go transient. But this transient flow is likely to be real.
This FAQ might help now: https://www.cfd-online.com/Wiki/Ansy...gence_criteria
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
DEM Particles protruding through walls | connor.dio12 | STAR-CCM+ | 1 | March 2, 2023 10:29 |
dsmcFoam setup | hherbol | OpenFOAM Pre-Processing | 1 | November 19, 2021 01:52 |
UDF for deleting particles in DPM | imanmirzaii | Fluent UDF and Scheme Programming | 12 | November 25, 2020 19:27 |
Boundary Conditions k-omega-SST with slip walls | shock77 | OpenFOAM Running, Solving & CFD | 6 | October 23, 2020 16:57 |
[DPM-UDF] Re-injecting escaping particles at different position | CeesH | Fluent UDF and Scheme Programming | 7 | May 13, 2020 10:34 |