View Issue Details

IDProjectCategoryView StatusLast Update
0003324OpenFOAMBugpublic2019-08-09 14:57
Reporterhosseinheydari2013Assigned Towill 
PriorityurgentSeveritymajorReproducibilityalways
Status closedResolutionfixed 
PlatformLINUXOSRed Hat Enterprise Linux ServerOS Version7.6 (Maipo)
Product Version6 
Fixed in Version7 
Summary0003324: spray generation inconsistency
DescriptionI am trying to develop my own solver and I noticed something strange and then I was curious to know whether I am doing something wrong or is it a bug in OF. And as it turns out, lagrangian spray generation for the serial and parallel runs in OF are not consistent with each other at all!!. I set up a test with following conditions:

serial case:
OF6
solver: sprayFoam
testCase: original aachenBomb (OpenFOAM-6/tutorials/lagrangian/sprayFoam/aachenBomb)

parallel case:
OF6
solver: sprayFoam (mpirun -np 4 sprayFoam -parallel)
testCase: original aachenBomb (OpenFOAM-6/tutorials/lagrangian/sprayFoam/aachenBomb)

The parallel spray is tilted to one side. I assume there should be an issue with the parallelization of the injection algorithm in OF. I attached two animations of the spray top view which shows the discrepancy between them. Also note that this error is not due to having multiple processors injecting lagrangian particles because I decomposed my domain in a way that only one processor would inject the lagrangian particles.

 
This is really annoying because suppose somebody develops a solver in serial version first. Then he wants to proceed further by parallelizing the solver and there will be no possibility to generate an identical testcase for both versions to compare. Moreover, I mentioned my OS is RedHat which is for a HPC center however I tested this on our own local cluster running openSUSE 13.2 (Harlequin) and the same issue persists.
Steps To ReproduceAs I described in Description it is the original tutorial of OF. One may find it in "OpenFOAM-6/tutorials/lagrangian/sprayFoam/aachenBomb".
TagsNo tags attached.

Activities

hosseinheydari2013

2019-08-09 10:35

reporter  

parallel.gif (613,823 bytes)
serial.gif (654,923 bytes)

will

2019-08-09 14:02

manager   ~0010663

Can you upload the decomposeParDict you used, please?

hosseinheydari2013

2019-08-09 14:11

reporter   ~0010664

Hiiiii will,

Thanks for your reply. Attached you may find it.

decomposeParDict (1,154 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

numberOfSubdomains 4;

method         simple;

simpleCoeffs
{
    n               ( 1 4 1 );
    delta           0.001;
}

hierarchicalCoeffs
{
    n               (8 1 1 );
    delta           0.001;
    order           xyz;
}

manualCoeffs
{
    dataFile        "";
}

distributed     no;

roots           ( );

// ************************************************************************* //
decomposeParDict (1,154 bytes)

will

2019-08-09 14:57

manager   ~0010666

Reproduced the problem in version 6. This has been fixed in version 7 and dev. I don't know the exact reason for the issue in version 6, but the cone and cone nozzle injectors had a number of problems with them and were were combined and significantly overhauled by the following commits.

https://github.com/OpenFOAM/OpenFOAM-dev/commit/06d8f79814f2ed2e8c476bbdfc47a98093a7529d
https://github.com/OpenFOAM/OpenFOAM-dev/commit/b785c988c6abc1580da24b2ed39d482791ab0c0c

Issue History

Date Modified Username Field Change
2019-08-09 10:35 hosseinheydari2013 New Issue
2019-08-09 10:35 hosseinheydari2013 File Added: parallel.gif
2019-08-09 10:35 hosseinheydari2013 File Added: serial.gif
2019-08-09 14:02 will Note Added: 0010663
2019-08-09 14:11 hosseinheydari2013 File Added: decomposeParDict
2019-08-09 14:11 hosseinheydari2013 Note Added: 0010664
2019-08-09 14:57 will Assigned To => will
2019-08-09 14:57 will Status new => closed
2019-08-09 14:57 will Resolution open => fixed
2019-08-09 14:57 will Fixed in Version => 7
2019-08-09 14:57 will Note Added: 0010666