CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

what "If" condition means in rebound

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 12, 2014, 10:18
Unhappy what "If" condition means in rebound
  #1
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 12
brbbhatti is on a distinguished road
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

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

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<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
brbbhatti is offline   Reply With Quote

Reply

Tags
printstack, rebound, splash

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
UDF for degassing boundary condition peaker007 Fluent UDF and Scheme Programming 5 November 23, 2015 13:55
Internal flow operating condition? kookguy FLUENT 2 June 26, 2014 01:15
conjugate boundary condition Daniel_Khazaei OpenFOAM Programming & Development 0 December 31, 2013 14:11
Flow simulation in wind tunnel mukut OpenFOAM 15 September 11, 2013 23:43
Bai model droplet-wall impact / elastic rebound crevoise STAR-CD 3 November 15, 2012 15:06


All times are GMT -4. The time now is 00:46.