Reply To: Additional Force Field
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Additional Force Field › Reply To: Additional Force Field
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.