Parallel using icoLagrangianFoam
 

Hey Foamers:

It is very my pleasure to discuss the parallel running of icoLagrangianFoam.

When I use icoLagrangianFoam, at first I use the command' decomposePar' and there are two folders of 'processor0' and 'processor1' to be created. And then I used the command of ''mpriun -np icoLagrangianFoam'.

The case is run, but there is no result to be saved in the folder of 'processor1'and'processor0', and the result is saved under the folder of case.

I feel a little confused whether the icoLangrangianFoam could be run paralleled or the result is saved under the folder of case other in the 'processor*`. Would anyone like to help me for the question?

Thanks and best wishes!

Try again with
mpriun -np 2 icoLagrangianFoam -parallel
(what you did was run two serial runs in parallel

Bernhard

Sincerely thanks! And it is work now. But for running paralleled with two processors is even slower than one processor.

For small cases that is not surprising, because the communication needs more time than the actual calculation

Thanks a lot

Long time no see. And wish all foamers are fine.

There are big bugs on "icoLagrangianFoam" in the case of paralleling computing. There is no problem for the case are computed using one process. If the mesh is decomposed in the z direction for cavity, the particle could be moved, but the result is different with case of using one process, it is the same situation for the y direction. The worst thing is that the particles can not be moved if it is decomposed in x direction after the computation starts some steps. It is absolutely there is some problem for the particle transfered between processors.

Does anyone also find the bug?
It is very sad I know few about the part of transfer field of particles. Who would like to give me some advice on it?

Thanks and best wishes!

Before we start talking: which version? (OpenFOAM and icoLagrangianFoam)

I just had a quick look at the latest version and noticed that in method HardBallParticle::hitProcessorPatch it reads
td.switchProcessor=false;
where in my opinion it should be
td.switchProcessor=true;
but I'm not 100% sure whether this will fix it and I havn't got the time to test it

Bernhard

It is icoLangrangianFoam for OpenFoam 1.5.

In the HardBallParticle.C, there is no lines in the member function of HardBallParticle::hitProcessorPatch. Maybe I should add these lines. I will try it. Thanks a lot.

There are two such methods. One is empty. The other not.

Hi Bernhard,

Yes, this solves the problem.
However, if I set useMomentumSource to 1, i.e. two-way coupling of solid and fluid phase, a extremely high value of Co will be achieved. This obviously not physical. This may be dedicated to the suddenly appearance of particles in a new processor.

Any idea to solve this?

Thanks.

Roro

Does this happen only in parallel runs? (I doubt it)

I think it IS physical: particles with a momentum appear out of nowhere (OK. THAT is unphysical). What should the fluid do? Suffer silently (that would be unphysical) or try to incorporate the additional momentum (that is what you're seeing)

Solutions all have to do with the injector:
- inject less (volume fraction of the particles should be well below 10% for the solver to be valid)
- inject slower
- write a different injector (on the patch, one that injects with no relative velocity to the fluid)

Bernhard

Hi, Bernard,

Yes, it's not the parallel problem. Particles can smoothly pass the processorPatch.

Thanks.

Roro

problems in parallel
 

Hi All,

I downloaded the icoLagrangianFoam from (http://openfoam-extend.svn.sourcefor...agrangianFoam/) and everything compiled fine in OF-1.5-dev

I tried to run icoLagrangianFoam in parallel and it failed, giving the following error messages:

I thought it might be a recent update of my kernel, but a quick parallel run of interFoam on the damBreak case finished without any problems.

I ran a parallel run using
Im sure its something I'm doing wrong. Any help?

Dan

Hi Dan!

Could you try the same with the rhoTurbTwinParcelFoam? Just to make sure whether the problem is with the solver or with the library?

Bernhard

same error...
 

Thanks for the reply, I tried the rhoTurbTwinParcelFoam solver and it compiled fine. When i ran the test case received the error:

the code works fine in serial.

I then went to the lagrangian folder in the $FOAM_SRC and performed and "svn update" and then recompiled the lagrangian files with the ./Allwmake script. Still the same error. Is there another folder I should update? Should I do a complete update of 1.5-dev?

Dan

update and recompile did not help
 

I just wanted to make sure it wasn't something else that was out of date and needed to be updated. I updated all of of-1.5-dev through svn and recompiled the changes, tried to run parallel again with icoLagrangianFoam and same error messages.

I'm working on 64bit ubuntu 9.10

Dan

This seems to be a general problem with the 1.5-dev (I'm getting it too). The nature of the stack-trace indicates that somewhere memory is freed twice somewhere

Bernhard

Hi All,

Does anyone know how to get around the failure with resulting stack trace problem when running icoLagrangianFoam in parallel?.
I am also getting this with OF 1.5-dev. It works fine on one processor, but as soon as I try to run it in parallel I also get this error.

I also noted that icoLagrangianFoam is not included in OF 1.6-ext. I tried to move it across but when I tried to compile it with (wmake icoLagrangianFoam ) it complains about a missing CintDefs.H file.

Thanks in advance (again),

Paul.

Hi Paul!

icoLagrangianFoam should be in 1.6-ext, but not in the solver directory but in $FOAM_TUTORIALS/lagrangian

Bernhard

Hi Bernhard,

I compiled the icoLagrangianFoam files in $FOAM_TUTORIALS/lagrangian of OF 1.6-ext and everything now works.

Thanks again,

Paul.