CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Error while compiling icoLagrangian Foam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 13, 2009, 03:28
Default hi all, Thanks Bernhard I w
  #21
Member
 
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 8
sachin is on a distinguished road
hi all,
Thanks Bernhard
I would like to add thermal effect for drag force calculation and eventually to the momentum change in particle....so would like to know where are the calculation for smoment_ is done so that i can add these term...smoment_ is actually defined just above tsource calculation where it is used then where are the calculations for smoment_ done

thanks in advance
sachin
sachin is offline   Reply With Quote

Old   September 29, 2009, 10:35
Default Parallel processing hangs-icoLagrangianFoam
  #22
New Member
 
Suraj Deshpande
Join Date: Mar 2009
Location: Madison, WI, USA
Posts: 18
Rep Power: 8
suraj is on a distinguished road
Hello All,
I am Suraj. I am trying to integrate Lagrangian particle tracking with interFoam (OpenFOAM-1.5). The code icoLagrangianFoam (written for OpenFOAM-1.5) was really very helpful.
Now, I am ready with my code but it does not work very well in parallel .

I added some more physical variables apart from U (ie, pd and rho) to the code and changed the equations of motion. What the code does is that it runs in parallel for some time steps and after that it just does not do anything, as if it is in an infinite loop. It also does not respond to the changes made to controlDict. Is that an MPI issue? Is it common for particle codes (I did not see decomposeParDict in any of Lagrangian code tutorials)?

It would be very helpful to have your suggestion on this.

Another thing I saw in the code was the time step passed to the the HardBallParticle::updateProperties(...) function (ie the line updateProperties(deltaT,data,cellI,face()); ) . I think it should be dt, as defined in HardBallParticle::move(...) function.


Thanks,
Regards,
Suraj Deshpande
suraj is offline   Reply With Quote

Old   September 29, 2009, 18:01
Default
  #23
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by suraj View Post
Hello All,
I am Suraj. I am trying to integrate Lagrangian particle tracking with interFoam (OpenFOAM-1.5). The code icoLagrangianFoam (written for OpenFOAM-1.5) was really very helpful.
Now, I am ready with my code but it does not work very well in parallel .

I added some more physical variables apart from U (ie, pd and rho) to the code and changed the equations of motion. What the code does is that it runs in parallel for some time steps and after that it just does not do anything, as if it is in an infinite loop. It also does not respond to the changes made to controlDict. Is that an MPI issue? Is it common for particle codes (I did not see decomposeParDict in any of Lagrangian code tutorials)?

It would be very helpful to have your suggestion on this.
At first: which version of OF are you using and which version of the solver?

One word of warning: I'm not planing to port the solver to OF 1.6 as starting with 1.5 there are example solvers that use lagrangian particles in the distribution. The 1.5-version of icoLagrangianFoam was just a straight port from the 1.4.1-version and didn't use the new particle classes.

Could you be more specific on the topic "hangs": what is the last output (that would help to tell in which phase the problem occurs)?

Have you ever run a debug-version of the code? Some problems come to light when running one (if you do any serious development you should have one)

Does this always happen or only with certain decompositions? Try some different decompositions

Quote:
Originally Posted by suraj View Post
Another thing I saw in the code was the time step passed to the the HardBallParticle::updateProperties(...) function (ie the line updateProperties(deltaT,data,cellI,face()); ) . I think it should be dt, as defined in HardBallParticle::move(...) function.
That may be right. I'll have a look

Bernhard
gschaider is offline   Reply With Quote

Old   October 2, 2009, 18:23
Default
  #24
New Member
 
Suraj Deshpande
Join Date: Mar 2009
Location: Madison, WI, USA
Posts: 18
Rep Power: 8
suraj is on a distinguished road
Dear Bernhard,

I am using OpenFOAM-1.5. I am making changes to the standard interFoam solver that came with it. I have not used the debug version of the code.

The code stops while moving the particles. To be more precise, it stops at the statement: this->move(td) in IncompressibleCloud.C .

I have seen this happen with larger (about half a million cells) cases; smaller 2D cases do not stop. Also, interestingly, the case hangs at different times for different number of decompositions (it hangs much earlier for cases decomposed into large number of subdomains eg 32 and does not hang at all for single processor cases, no matter how big the case is). This is what leads me into thinking if this is a problem with inter-processor communication in my code.

This is what I see in the log:

Courant Number mean: 0.0022608 max: 0.166463
Moving Particles
529 Particles moved. 0 walls hit. 0 particles left the model. 1 particles injected.
Time = 4.0672

MULES: Solving for gamma
Liquid phase volume fraction = 0.517803 Min(gamma) = -2.17137e-09 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.517803 Min(gamma) = -1.04872e-10 Max(gamma) = 1
DICPCG: Solving for pd, Initial residual = 0.00124773, Final residual = 4.05203e-05,
No Iterations 3
DICPCG: Solving for pd, Initial residual = 0.000151898, Final residual = 6.90998e-06
, No Iterations 6
DICPCG: Solving for pd, Initial residual = 3.09272e-05, Final residual = 9.78173e-08
, No Iterations 43
time step continuity errors : sum local = 3.48316e-10, global = -7.79987e-11, cumulat
ive = 3.38557e-09
ExecutionTime = 538.01 s ClockTime = 541 s

Courant Number mean: 0.00226081 max: 0.167813
Moving Particles


...... hangs here

Also, there are times when the the number of particles crossing processor patch goes astronomically high (3-4Million) for number of droplets as small as 150. I am not sure if this an issue too...

Thanks,
Suraj

Last edited by suraj; October 2, 2009 at 19:57.
suraj is offline   Reply With Quote

Old   October 9, 2009, 13:00
Default
  #25
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by suraj View Post
Dear Bernhard,

I am using OpenFOAM-1.5. I am making changes to the standard interFoam solver that came with it. I have not used the debug version of the code.

The code stops while moving the particles. To be more precise, it stops at the statement: this->move(td) in IncompressibleCloud.C .

I have seen this happen with larger (about half a million cells) cases; smaller 2D cases do not stop. Also, interestingly, the case hangs at different times for different number of decompositions (it hangs much earlier for cases decomposed into large number of subdomains eg 32 and does not hang at all for single processor cases, no matter how big the case is). This is what leads me into thinking if this is a problem with inter-processor communication in my code.

This is what I see in the log:

Courant Number mean: 0.0022608 max: 0.166463
Moving Particles
529 Particles moved. 0 walls hit. 0 particles left the model. 1 particles injected.
Time = 4.0672

MULES: Solving for gamma
Liquid phase volume fraction = 0.517803 Min(gamma) = -2.17137e-09 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.517803 Min(gamma) = -1.04872e-10 Max(gamma) = 1
DICPCG: Solving for pd, Initial residual = 0.00124773, Final residual = 4.05203e-05,
No Iterations 3
DICPCG: Solving for pd, Initial residual = 0.000151898, Final residual = 6.90998e-06
, No Iterations 6
DICPCG: Solving for pd, Initial residual = 3.09272e-05, Final residual = 9.78173e-08
, No Iterations 43
time step continuity errors : sum local = 3.48316e-10, global = -7.79987e-11, cumulat
ive = 3.38557e-09
ExecutionTime = 538.01 s ClockTime = 541 s

Courant Number mean: 0.00226081 max: 0.167813
Moving Particles


...... hangs here

Also, there are times when the the number of particles crossing processor patch goes astronomically high (3-4Million) for number of droplets as small as 150. I am not sure if this an issue too...
But sometimes particle cross to other processors sucessfully? No idea. Sprinkle the HardballParticle.C with "

Pout << something << endl;

statements (especially the stuff concerning parallelization) to find what might be the problem

Bernhard
gschaider is offline   Reply With Quote

Old   October 9, 2009, 13:05
Default
  #26
New Member
 
Suraj Deshpande
Join Date: Mar 2009
Location: Madison, WI, USA
Posts: 18
Rep Power: 8
suraj is on a distinguished road
Thanks Bernhard,
I will try this and let you know if I see something interesting.

Thanks again,
Regards,
Suraj Deshpande
suraj is offline   Reply With Quote

Old   October 13, 2009, 20:46
Default
  #27
New Member
 
Suraj Deshpande
Join Date: Mar 2009
Location: Madison, WI, USA
Posts: 18
Rep Power: 8
suraj is on a distinguished road
Hi Bernhard,
I used the Pout statements as you suggested. I think I have a problem in Cloud<HardBallParticle>::move(TrackingData& td) of Cloud.C. The code stops at the statement "combineReduce(allNTrans, combineNsTransPs());" .
I will try to see what this function does.
Any suggestions?

Thanks,
Suraj Deshpande
suraj is offline   Reply With Quote

Old   October 16, 2009, 05:46
Default
  #28
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by suraj View Post
Hi Bernhard,
I used the Pout statements as you suggested. I think I have a problem in Cloud<HardBallParticle>::move(TrackingData& td) of Cloud.C. The code stops at the statement "combineReduce(allNTrans, combineNsTransPs());" .
I will try to see what this function does.
Any suggestions?
This is hardcore (that part should work in any case). Do ALL processors reach that statement?

Which version of 1.5 are you using (unpatched, 1.5.x or something in between)? (I think there was a fix sometime it affected parallel running of dieselFoam which might have been the same problem)

Bernhard
gschaider is offline   Reply With Quote

Old   August 20, 2010, 06:06
Default
  #29
Senior Member
 
xinguang cui
Join Date: Mar 2009
Posts: 116
Rep Power: 8
flying is on a distinguished road
Hello OpenFoamer:

I wish all of you are enjoying the life now.

I am implementing a model in the icoLagrangianFoam which can solve the sub-micro particles. I need to generate the random number. I try to use the class random().scalar01() in the HardBallParticle::updateProperties, but it is pity that there is some error as following:

HardBallParticle.C:181: error: request for member 'scalar01' in 'random()', which is of non-class type 'long int'

who can give me some hint on it?

Last edited by flying; August 20, 2010 at 08:49.
flying is offline   Reply With Quote

Old   August 20, 2010, 08:47
Default
  #30
Senior Member
 
xinguang cui
Join Date: Mar 2009
Posts: 116
Rep Power: 8
flying is on a distinguished road
Hello OpenFoamer:

I wish all of you are enjoying the life now.

I am implementing a model in the icoLagrangianFoam which can solve the sub-micro particles. I need to generate the random number. I try to use the class random().scalar01() in the HardBallParticle::updateProperties, but it is pity that there is some error as following:

HardBallParticle.C:181: error: request for member 'scalar01' in 'random()', which is of non-class type 'long int'

who can give me some hint on it?
flying is offline   Reply With Quote

Old   February 6, 2013, 02:28
Default
  #31
Senior Member
 
Jie
Join Date: Jan 2010
Location: Australia
Posts: 119
Rep Power: 7
jiejie is on a distinguished road
Quote:
Originally Posted by flying View Post
Hello OpenFoamer:

I wish all of you are enjoying the life now.

I am implementing a model in the icoLagrangianFoam which can solve the sub-micro particles. I need to generate the random number. I try to use the class random().scalar01() in the HardBallParticle::updateProperties, but it is pity that there is some error as following:

HardBallParticle.C:181: error: request for member 'scalar01' in 'random()', which is of non-class type 'long int'

who can give me some hint on it?
I tried to include the Random.H header file, but it doesn't really help either???

Any help will be appraciated.
jiejie is offline   Reply With Quote

Old   February 6, 2013, 05:33
Default
  #32
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by jiejie View Post
I tried to include the Random.H header file, but it doesn't really help either???

Any help will be appraciated.
Case is important in C++. "random" is a system-function (that returns a long int - thus the error). "Random" is a class in Foam.

BTE: the way you call it will probably always give you the same number - but I would have to look at the implementation to be 100% sure
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   February 6, 2013, 09:45
Default
  #33
Senior Member
 
Jie
Join Date: Jan 2010
Location: Australia
Posts: 119
Rep Power: 7
jiejie is on a distinguished road
Quote:
Originally Posted by gschaider View Post
Case is important in C++. "random" is a system-function (that returns a long int - thus the error). "Random" is a class in Foam.

BTE: the way you call it will probably always give you the same number - but I would have to look at the implementation to be 100% sure
Thanks gschaider. Hope to hear from you soon.
jiejie is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Errors when Compiling UDF: error C2040/error C2099 Julian K. FLUENT 1 December 21, 2008 01:23
Compiling MRFSimple foam error mrangitschdowcom OpenFOAM Installation 0 September 15, 2008 15:09
FOAM FATAL IO ERROR msrinath80 OpenFOAM Running, Solving & CFD 4 July 30, 2008 10:06
Adding collision to icoLagrangian foam alan_fergus OpenFOAM Running, Solving & CFD 0 December 12, 2006 08:55
FOAM FATAL IO ERROR sita OpenFOAM Running, Solving & CFD 2 August 23, 2005 04:37


All times are GMT -4. The time now is 17:06.