[DanielRoman]

Hello everybody. I have modified the shallowWaterFoam solver in order to add some restrictive stresses and address a land slide flow model which momentum equation is quite similar to these from the shallow water problem. I compiled the solver and run simple a case (the same squareBump), it begins calculating but then after some iterations I got this error:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
 #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  ? at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? at ??:?

My momentum equation reads:

if (pimple.momentumPredictor())
            {
                solve(hUEqn == (-1.0/rho)*(h*kap*rho*gTz*soli*fvc::div(Tsd)
                               +sign(graUyx)*h*rho*gTz*soli*fvc::div(Tss)
                               +sign(graUxy)*h*rho*gTz*soli*fvc::div(Tsi)
                               +h*rho*gTz*soli*tan(phib)*(GUx+GUy)
                               -nuf*mu*h*divgUT
                               -rho*gTrans*h)
                     );
            
                // Constrain the momentum to be in the geometry if 3D geometry
                if (mesh.nGeometricD() == 3)
                {
                    hU -= (gHat & hU)*gHat;
                    hU.correctBoundaryConditions();
                }
            }

I guess that something is happening with a negative square root argument, due to the sigFpe Error but if someone has a clear idea I will appreciate any suggestions. Thanks in advance!



Regards,


Daniel.


PS: I attached couple files. One is the "log" file and the other is the solver landSlideFoam.C