[vorsorge1988]

Hello Foamers,

I want to implement the wallLubricationForce in the twoPhaseEulerFoam solver. The compiling was succesfull, but if i want to solve the bubbleColumn-Case i got the following error:

Courant Number mean: 0.0596327 max: 0.205145
Max Ur Courant Number = 0.128547
Time = 0.002

MULES: Solving for alpha1
MULES: Solving for alpha1
Dispersed phase volume fraction = 0.0001 Min(alpha1) = 0 Max(alpha1) = 0.5
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/tpef_wall"
#5
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/tpef_wall"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/tpef_wall"

I think the errorsource is the walldistance (yw) in my liftDragCoeff-File which i get from wallDist.H. I tried a lot, but my knowlegde in C++ is nearly zero. ^^

volVectorField Ur(U2 - U1);
volScalarField magUr(mag(Ur) + residualSlip);
volVectorField nw = wallDistReflection(mesh).n();
volScalarField yw = wallDist(mesh).y();

volScalarField Cwall(0.5*CwlEo*d*(1.0/sqrt(yw)-1.0/sqrt(Dpipe-yw)));

volScalarField WallForceMag(CwlEo*alpha1*dichte*sqrt(mag(Ur-(Ur & nw)*nw)));

volVectorField WallForce(WallForceMag*(-nw));

Thank you for your help and ideas

Best regards Daniel