Reply To: Additional Force Field
Ok, you will need to use dynamics that include a forcing scheme. This is not the case for
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
forcing::* into a new one for your situation.