# what "If" condition means in rebound

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

August 12, 2014, 10:18
what "If" condition means in rebound
#1
Member

bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 11
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.

Quote:
 case PatchInteractionModel::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

Quote:
 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

Quote:
 It is splash and its not working is an info command
Quote:
 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 > > > > >:atchData_splash(Foam::ReactingMultiphaseParcel > > > 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 > > > > > >::correct(Foam::ReactingMultiphaseParcel > > >&, 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

 Tags printstack, rebound, splash