Reply To: Additional Force Field

[Henderson_C]

Hello Adrian,

Thank you for your reply. I have add the force field but I still don’t observe any convective transport (Change in velocity). My approach was based on the Poiseuille Flow 2D Test case (Flowtype forced):

/// Code snippet
typedef D2Q9<FORCE, POROSITY, VELOCITY_NUMERATOR, VELOCITY_DENOMINATOR, CONTACT_DETECTION> NSDESCRIPTOR;
// Prepare Lattice
T Ly = converter.getLatticeLength( charL/2. );
std::vector<T> poiseuilleForce( 2,T() );
poiseuilleForce[0] = 2.*converter.getLatticeViscosity()* converter.getCharLatticeVelocity() / ( Ly*Ly ); // This rule came from fully developed flow for parallel plate
AnalyticalConst2D<T,T> force( poiseuilleForce );

// Initialize force
NSlattice.defineField<FORCE>(bulkIndicator, force);
NSlattice.defineField<FORCE>(superGeometry.getMaterialIndicator({2,5}), force);
// End of the code

I have tried the forcedBGK dynamics and it worked but I want to use the porousparticleBGK dynamics. Any suggestions?

Regards,