View Issue Details

IDProjectCategoryView StatusLast Update
0000383OpenFOAMBugpublic2012-01-09 11:56
Reporteruser322Assigned Touser2 
PriorityhighSeveritycrashReproducibilityalways
Status resolvedResolutionno change required 
Platformx86_64OSUbuntuOS Version10.10
Summary0000383: InterDyMFoam with AMI crashes in parallel
DescriptionI made following project:
2D multiphase with rotating impeller based with Arbitrary Mesh Interface (AMI).
solver: InterDymFoam

Everything is going good in single processor calculation. But when I change to parallel computation program crashes in first iteration. I have test it on one processor with 4 cores. Also i check it on other computer with Suse linux 11. The situation was the same.

Is there any special boundary condition parameters or solution methods wich should be chosen to parallel computation
Steps To ReproduceI run following commands:
setFields
decomposePar
mpirun -np 4 interDyMFoam -parallel > log &
Additional Informationthis is a log from computation:
Build : 2.1.x-eb976ba31c36
Exec : interDyMFoam -parallel
Date : Jan 08 2012
Time : 15:11:59
Host : "tomekGRANT"
PID : 10143
Case : /home/tomek/OpenFOAM/test
nProcs : 4
Slaves :
3
(
"tomekGRANT.10144"
"tomekGRANT.10145"
"tomekGRANT.10146"
)

Pstream initialized with:
    floatTransfer : 0
    nProcsSimpleSum : 0
    commsType : nonBlocking
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone WIRNIK_ZONE
Reading field p_rgh

Reading field alpha1

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 152 source faces and 88 target faces
AMI: Patch source weights min/max/average = 1, 1.00023, 1.00009
AMI: Patch target weights min/max/average = 1.00018, 1.00028, 1.00023
Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu 0.09;
    C1 1.44;
    C2 1.92;
    sigmaEps 1.3;
}


Reading g
Calculating field g.h


PIMPLE: Operating solver in PISO mode

time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0 max: 0
Time = 0.0005

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.0005 transformation: ((0 0 0) (0.999969 (0 0 0.0078539)))
AMI: Creating addressing and weights between 152 source faces and 88 target faces
AMI: Patch source weights min/max/average = 1, 1.00023, 1.00009
AMI: Patch target weights min/max/average = 1.00019, 1.00028, 1.00023
Execution time for mesh.update() = 0.33 s
MULES: Solving for alpha1
Liquid phase volume fraction = 0.314942 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.314942 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.314942 Min(alpha1) = 0 Max(alpha1) = 1
[2] #0 Foam::error::printStack(Foam::Ostream&)[3] #0 Foam::error::printStack(Foam::Ostream&) in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigHandler(int) in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #2 in "/lib/libc.so.6"
[2] #3 Foam::inv(Foam::Field<Foam::SymmTensor<double> >&, Foam::UList<Foam::SymmTensor<double> > const&) in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2 in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #4 in "/lib/libc.so.6"
[3] #3 Foam::inv(Foam::Field<Foam::SymmTensor<double> >&, Foam::UList<Foam::SymmTensor<double> > const&)
[2] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[2] #5 in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #4
[2] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[2] #6
[3] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[3] #5
[3] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[3] #6
[2] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[2] #7
[2] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[2] #8 __libc_start_main
[3] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[3] #7 in "/lib/libc.so.6"
[2] #9

[2] in "/home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[tomekGRANT:10145] *** Process received signal ***
[tomekGRANT:10145] Signal: Floating point exception (8)
[tomekGRANT:10145] Signal code: (-6)
[tomekGRANT:10145] Failing at address: 0x3e8000027a1
[tomekGRANT:10145] [ 0] /lib/libc.so.6(+0x33c20) [0x7f99b9c71c20]
[tomekGRANT:10145] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7f99b9c71ba5]
[tomekGRANT:10145] [ 2] /lib/libc.so.6(+0x33c20) [0x7f99b9c71c20]
[tomekGRANT:10145] [ 3] /home/tomek/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam3invERNS_5FieldINS_10SymmTensorIdEEEERKNS_5UListIS2_EE+0x546) [0x7f99bad0de96]
[tomekGRANT:10145] [ 4] interDyMFoam() [0x45a5c2]
[tomekGRANT:10145] [ 5] interDyMFoam() [0x498170]
[tomekGRANT:10145] [ 6] interDyMFoam() [0x4988f0]
[tomekGRANT:10145] [ 7] interDyMFoam() [0x42b4b5]
[tomekGRANT:10145] [ 8] /lib/libc.so.6(__libc_start_main+0xfe) [0x7f99b9c5cd8e]
[tomekGRANT:10145] [ 9] interDyMFoam() [0x4211a9]
[tomekGRANT:10145] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 10145 on node tomekGRANT exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------


In attachemnt my case files

Best Regards.

Tomasz
TagsNo tags attached.

Activities

user322

2012-01-08 14:26

 

test.zip (1,806,961 bytes)

user2

2012-01-09 11:56

  ~0000912

> Is there any special boundary condition parameters or solution methods wich should be chosen to parallel computation

There are no special settings required when running AMI cases in parallel

The problem here arises due to your mesh: it is not 2-D - if you zoom in, you can see that the front and back planes are not aligned. After generating a new mesh by extruding the 'TOP' patch, all works well in serial and parallel modes.

Issue History

Date Modified Username Field Change
2012-01-08 14:26 user322 New Issue
2012-01-08 14:26 user322 File Added: test.zip
2012-01-09 11:56 user2 Note Added: 0000912
2012-01-09 11:56 user2 Status new => resolved
2012-01-09 11:56 user2 Fixed in Version => 2.1.x
2012-01-09 11:56 user2 Resolution open => no change required
2012-01-09 11:56 user2 Assigned To => user2