Foam::error::printStack after adding New Drag model

surajkvs · February 26, 2018, 08:40

I have added one new drag model and compiled it successfully.
After that my new solver also I compiled.
Than I run the default case 'Goldschmidt' with the new drag model it run succesfully in MPPICFoam solver.
But when I am running my own case it gets stop somewhere between 0.1 to 0.2 sec of simulation.and shows following error.

.................................................. ...........
#0 Foam::error:

rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4 pow in "/lib/x86_64-linux-gnu/libm.so.6"
#5 Foam::EmmsAtdDragForce<Foam::KinematicCloud<Foam:: Cloud<Foam::MPPICParcel<Foam::KinematicParcel<Foam :

article> > > > >::calcCoupled(Foam::MPPICParcel<Foam::KinematicPa rcel<Foam:

article> > const&, double, double, double, double) const at ??:?
#6 Foam::ParticleForceList<Foam::KinematicCloud<Foam: :Cloud<Foam::MPPICParcel<Foam::KinematicParcel<Foa m:

article> > > > >::calcCoupled(Foam::MPPICParcel<Foam::KinematicPa rcel<Foam:

article> > const&, double, double, double, double) const at ??:?
#7 ? at ??:?
#8 ? at ??:?
#9 ? at ??:?
#10 ? at ??:?
#11 ? at ??:?
#12 ? at ??:?
#13 ? at ??:?
#14 ? at ??:?
#15 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#16 ? at ??:?
Floating point exception (core dumped)
.................................................. ....................................

Can anybody please help me.I have read few PrintStack problems but unable to find my solution.

anuragm · February 26, 2018, 09:54

Usually happens when there is division by zero or similar error. I'd suggest recompiling your code using the following flag in your Make/options

Code:

 -DFULLDEBUG -g -O0

This will help you to trace the source of your error. Perhaps you will need to ensure your model is valid in extreme cases by adding a small term to the denominator.

surajkvs · March 10, 2018, 12:11

Thanks a lot Anurag for reply,Actually it worked when I changed the mesh.

February 26, 2018, 08:40	Foam::error::printStack after adding New Drag model	#1
surajkvs New Member SURAJ Join Date: Jun 2017 Location: IIT KANPUR,India Posts: 15 Rep Power: 8	I have added one new drag model and compiled it successfully. After that my new solver also I compiled. Than I run the default case 'Goldschmidt' with the new drag model it run succesfully in MPPICFoam solver. But when I am running my own case it gets stop somewhere between 0.1 to 0.2 sec of simulation.and shows following error. .................................................. ........... #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 ? in "/lib/x86_64-linux-gnu/libm.so.6" #4 pow in "/lib/x86_64-linux-gnu/libm.so.6" #5 Foam::EmmsAtdDragForce<Foam::KinematicCloud<Foam:: Cloud<Foam::MPPICParcel<Foam::KinematicParcel<Foam :article> > > > >::calcCoupled(Foam::MPPICParcel<Foam::KinematicPa rcel<Foam:article> > const&, double, double, double, double) const at ??:? #6 Foam::ParticleForceList<Foam::KinematicCloud<Foam: :Cloud<Foam::MPPICParcel<Foam::KinematicParcel<Foa m:article> > > > >::calcCoupled(Foam::MPPICParcel<Foam::KinematicPa rcel<Foam:article> > const&, double, double, double, double) const at ??:? #7 ? at ??:? #8 ? at ??:? #9 ? at ??:? #10 ? at ??:? #11 ? at ??:? #12 ? at ??:? #13 ? at ??:? #14 ? at ??:? #15 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #16 ? at ??:? Floating point exception (core dumped) .................................................. .................................... Can anybody please help me.I have read few PrintStack problems but unable to find my solution.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Adding New phaseCompressible Turbulence Model	scram	OpenFOAM Programming & Development	34	March 25, 2020 08:08
Adding a custom turbulence model in OpenFOAM 3.0	Sylvain	OpenFOAM Programming & Development	3	July 18, 2016 08:15
Airfoil lift and drag using k-kl-omega turbulence model	hylleman	OpenFOAM Running, Solving & CFD	6	June 17, 2016 15:10
How to compile a new drag force in Lagrangian model	lilinmin	OpenFOAM Programming & Development	4	June 9, 2016 04:15
multiphase model and drag law	Yasmail AKARIOUH	FLUENT	0	April 29, 2008 07:44

February 26, 2018, 09:54		#2
anuragm Member Anurag Join Date: Aug 2014 Location: Germany Posts: 57 Rep Power: 11	Usually happens when there is division by zero or similar error. I'd suggest recompiling your code using the following flag in your Make/options Code: -DFULLDEBUG -g -O0 This will help you to trace the source of your error. Perhaps you will need to ensure your model is valid in extreme cases by adding a small term to the denominator. surajkvs likes this.

March 10, 2018, 12:11		#3
surajkvs New Member SURAJ Join Date: Jun 2017 Location: IIT KANPUR,India Posts: 15 Rep Power: 8	Thanks a lot Anurag for reply,Actually it worked when I changed the mesh.