Large droplet terminal velocity

[jbarr444]

Hi @slayoo,

Happy new year. I have been running into a persistent and annoying error: ValueError: Radii can be interpolated up to 0.006 m... I imagine this is because drops are getting larger than whatever table of values for terminal velocity is used can support. I'm currently using a gravitational kernel with perfect collection efficiency, so droplets are growing rather fast. Is there any way to change the terminal velocity scheme? I am not so interested in super realistic results right now, just playing with the mathematics of the system.

Thanks,
Jason

[slayoo]

H @jbarr444, happy new year!

One option could be to use the SimpleGeometric collision kernel which does not take into account the terminal velocity (it allows a simple multiplicative constant to be passed as ctor argument).

Another option could be to use the RogersYau formulation of the terminal velocity (as exemplified, e.g., here).

HTH,
Sylwester