pisoLagrangianFoam
 

Hi there!

I am trying to upgrade the pisoLagrangianFoam but I get an error in the momentum predictor because of units mismatch.
I tried adding the particles' density in order to solve the problem but the problem still persists.
My question is how to change the momentum equation in order to take into account the particles influence.
Thanks a lot in advance.
Regards,

Alex

This is the error I get (3217 is the density of the particles but apparently it is not taking into account the units):

--> FOAM FATAL ERROR:
incompatible dimensions for operation
[U[0 1 -2 0 0 0 0] ] == [(smom|3217)[1 -2 -2 0 0 0 0] ]

And this is the code


// Pressure-velocity PISO corrector
{
// Momentum predictor

fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
+ turbulence->divDevReff(U)
==
particles.momentumSource()/particles.rhop()
);

I solved it some time ago.
Thanks anyway!