# what "If" condition means in rebound

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

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

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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post peaker007 Fluent UDF and Scheme Programming 5 November 23, 2015 13:55 kookguy FLUENT 2 June 26, 2014 00:15 Daniel_Khazaei OpenFOAM Programming & Development 0 December 31, 2013 14:11 mukut OpenFOAM 15 September 11, 2013 22:43 crevoise STAR-CD 3 November 15, 2012 15:06

All times are GMT -4. The time now is 19:53.