[garcfd]

I just spotted a couple of clues in the source code:

MOMENTUM_EXCHANGE_FORCE_H
//Momentum Exchange Wen2014
return (f1*(c-pVel) + f2*(c+pVel));

//Momentum Exchange Ladd1994
return ((f1 + f2)*c);

[garcfd]

Thanks – can you point me to the paper (or code) which describes how the momentum exchange method was implemented (please)?

[garcfd]

also here…

wall reflection method.pdf

[garcfd]

I put some details of the approach here in case anyone wants to give it a try:
https://docs.google.com/presentation/d/1aXyArj2_nsMdsfY9GIQVvt6nI69F3QBoqEMKDnpeoEU

[garcfd]

Here is the full 3D animation of the Ferrari Dino case (just using bounce-back).
It was 100 million cells, ran on runpod GPU for 50k iterations, took 42 minutes.
Ferrari Dino (2)

[garcfd]

Here is an animation from some time ago – this one was (I think) 2D by the look of it.
Also just using bounce-back in this case, which can be seen by the low wall velocity.
Ferrari Dino animation

[garcfd]

porsche animation

[garcfd]

Hi Adrian, thanks for your comments. I’m not sure about my ability to code OpenLB myself, (my code is Fortran) but I would be happy to work with (explain the method) to someone else who is familiar with coding in OpenLB.
Regards Giles.

[garcfd]

Hi Mathias, that is very useful and interesting to know.

what I was wondering is what settings are required to get it to just run
since when I increase the velocity to say 10 m/sec it crashed, so expect
that have to change some solver settings for higher speeds? Thank Giles.

[Author]

Posts