[brbbhatti]

Hello Everyone,

I am trying to implement splash model in PatchInteraction Model by using the thermosurfacefilm model equations. I have created a new class in kinematic cloud and a new submodel Splash. I get to compiling the library and solver with new model successfully but i have a problem when i run the simulation.

It went to endless loop in the start. After studying the rebound sub model in StandardWallinteractionModel i come to the point that there is "if" condition in "Itrebound" which prevents the model from going into a continuous loop because if i remove that condition it goes to continuous loop too.

So its means i have to write some if condition for splash in StandardWallinteractionModel as it is written for the rebound model in it.
I have written some if condition similar to rebound for splash model and now it does not go to the continuous loop but give Printstack error.

Can anyone please help me to understand what does that if condition do in rebound

case PatchInteractionModel<CloudType>::itRebound:

{
keepParticle = true;
active = true;
vector nw;
vector Up;
this->owner().patchData(p, pp, trackFraction, tetIs, nw, Up);
// Calculate motion relative to patch velocity
U -= Up;
scalar Un = U & nw;
vector Ut = U - Un*nw;

if (Un > 0)
{
U -= (1.0 + e_)*Un*nw;
}

U -= mu_*Ut;
// Return velocity to global space
U += Up;
break;
}

the condition i have written for splash is

if (Uj > 0)
{
U -= (1.0 + e_)*Un*nj;
}

And how can i write the correct condition for splash phenomena to work ??.

And second this is the error now i am getting after writing this condition and running the simulation

It is splash and its not working is an info command

Time = 1

Solving cloud reactingCloud1

--> Cloud: reactingCloud1
Added 5 new parcels

It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
It is splash and its not working
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM_new.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM_new.so"
#2 in "/lib64/libc.so.6"
#3 Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMul tiphaseParcel<Foam::ReactingParcel<Foam::ThermoPar cel<Foam::KinematicParcel<Foam:article> > > > > >:atchData_splash(Foam::ReactingMultiphaseParcel <Foam::ReactingParcel<Foam::ThermoParcel<Foam::Kin ematicParcel<Foam:article> > > > const&, Foam:olyPatch const&, Foam::tetIndices const&, double, double, double, double, bool&) in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/lib/liblagrangianIntermediate_new_splash.so"
#4 Foam::StandardWallInteraction<Foam::KinematicCloud <Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam:: ReactingParcel<Foam::ThermoParcel<Foam::KinematicP arcel<Foam:article> > > > > > >::correct(Foam::ReactingMultiphaseParcel<Foam::Re actingParcel<Foam::ThermoParcel<Foam::KinematicPar cel<Foam:article> > > >&, Foam:olyPatch const&, bool&, double, Foam::tetIndices const&) in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/lib/liblagrangianIntermediate_new_splash.so"
#5
in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/bin/steadyParcelSolver_h_neu_1_splash"
#6
in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/bin/steadyParcelSolver_h_neu_1_splash"
#7
in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/bin/steadyParcelSolver_h_neu_1_splash"
#8
at steadyParcelSolver_h_neu_1_splash.C:0
#9
in "/home/bhb8fe/OpenFOAM/bhb8fe-2.1.0/platforms/linux64GccDPOpt/bin/steadyParcelSolver_h_neu_1_splash"
#10 __libc_start_main in "/lib64/libc.so.6"
#11
at /home/abuild/rpmbuild/BUILD/glibc-2.14.1/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception

any idea or help will be nice.

Regards,

Bilal