CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

wallLubricationForce in twoPhaseEulerFoamSolver

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   July 25, 2013, 10:48
Default wallLubricationForce in twoPhaseEulerFoamSolver
New Member
Daniel Wittig
Join Date: Oct 2012
Posts: 1
Rep Power: 0
vorsorge1988 is on a distinguished road
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/"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#2 in "/lib/x86_64-linux-gnu/"
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/"
#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"
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/tpef_wall"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/"
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
vorsorge1988 is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

All times are GMT -4. The time now is 13:39.