• OpenLB - Open source lattice Boltzmann code

September 6, 2023 at 11:20 am #7728

Keymaster

Ok, you will need to use dynamics that include a forcing scheme. This is not the case for PorousParticleBGK:


template<typename T, typename DESCRIPTOR, typename MOMENTA=momenta::BulkTuple>
using PorousParticleBGKdynamics = dynamics::Tuple<
  T, DESCRIPTOR,
  MOMENTA,
  equilibria::SecondOrder,
  collision::PorousParticle<collision::BGK,false>,
  dynamics::ExposePorousParticleMomenta
>;

As this dynamics already use a combination rule (dynamics::ExposePorousParticleMomenta) you can not directly insert one of the existing forcing schemes in src/dynamics/forcing.h. Depending on what exactly you want to model you will need to adapt your own combination rule to both expose the porous particle momenta and apply the forcing. But first I would recommend to check the literature / verify your model w.r.t. the specific forcing scheme you want to use.

Once you have selected a forcing scheme you can combine the two rules dynamics::ExposePorousParticleMomenta and forcing::* into a new one for your situation.

Reply To: Additional Force Field