[Taghi]

Hi all
I am trying to simulate a ducted propeller with pimpleDyMFoam solver. There is no problem with propeller alone. But when I try to add the duct and solve it again, solver stops.
I am now trying to figure the problem by using:

 runApplication moveDynamicMesh -checkAMI

(as openfoam tutorials suggested)

but the following message apears:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*   Windows 32 and 64 bit porting by blueCAPE: http://www.bluecape.com.pt   *\
|  Based on Windows porting (2.0.x v4) by Symscape: http://www.symscape.com   |
\*---------------------------------------------------------------------------*/
Build  : 4.x-ed69f631ce88
Exec   : C:/PROGRA~1/BLUECF~1/OpenFOAM-4.x/platforms/mingw_w64GccDPInt32Opt/bin/moveDynamicMesh.exe -checkAMI
Date   : May 14 2017
Time   : 21:54:08
Host   : "DESKTOP-PU3N81D"
PID    : 5112
Case   : C:/Users/Pouria/Desktop/propeller400
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing 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 innerCylinderSmall
Writing VTK files with weights of AMI patches.


PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

forces Rotor:
--> FOAM Warning : 
    From function Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const
    in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 794
    Cannot find any patch or group names matching propellerRotor
    Not including porosity effects
forces DAS:
--> FOAM Warning : 
    From function Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const
    in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 794
    Cannot find any patch or group names matching propellerDuct
    Not including porosity effects
Time = 5e-006
PIMPLE: iteration 1
--> FOAM Warning : 
    From function virtual bool Foam::solidBodyMotionFvMesh::update()
    in file solidBodyMotionFvMesh/solidBodyMotionFvMesh.C at line 228
    Did not find volVectorField U Not updating Uboundary conditions.
PIMPLE: iteration 2
    Point usage OK.
    Upper triangular ordering OK.
    Topological cell zip-up check OK.
    Face vertices OK.
  <<Number of duplicate (not baffle) faces found: 2. This might indicate a problem.
  <<Number of faces with non-consecutive shared points: 3. This might indicate a problem.
    Mesh topology OK.
    Boundary openness (3.72253e-016 -9.13045e-017 9.01232e-017) OK.
    Max cell openness = 4.35445e-016 OK.
    Max aspect ratio = 11.8457 OK.
    Minimum face area = 3.57658e-009. Maximum face area = 0.00294544.  Face area magnitudes OK.
    Min volume = 9.05772e-011. Max volume = 0.000110349.  Total volume = 1.31244.  Cell volumes OK.
    Mesh non-orthogonality Max: 64.9606 average: 9.12233
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.94105 OK.
    Mesh geometry OK.
Mesh OK.
Calculating AMI weights between owner patch: AMI1 and neighbour patch: AMI2
AMI: Creating addressing and weights between 45130 source faces and 46232 target faces
AMI: Patch source sum(weights) min/max/average = 0.0734142, 2.13295, 0.999935
AMI: Patch target sum(weights) min/max/average = 0, 1.84503, 0.999232
ExecutionTime = 39.753 s  ClockTime = 40 s



--> FOAM FATAL ERROR: 
Could not find U, p

    From function void Foam::functionObjects::forces::initialise()
    in file forces/forces.C at line 186.

FOAM exiting

and this is pimpleDyMFoam error:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : pimpleDyMFoam -parallel
Date   : May 15 2017
Time   : 12:40:46
Host   : "pouria"
PID    : 3167
Case   : /home/pouria/Desktop/propeller400
nProcs : 4
Slaves : 
3
(
"pouria.3168"
"pouria.3169"
"pouria.3170"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing 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 innerCylinderSmall

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading field p

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 45146 source faces and 46392 target faces
AMI: Patch source sum(weights) min/max/average = 0.565692, 2.10459, 1.00072
AMI: Patch target sum(weights) min/max/average = 0, 1.84453, 0.998807
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.06 average: 0.06
bounding epsilon, min: 0 max: 0.0495 average: 0.0495
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}

No MRF models present

Reading/calculating face velocity Uf

No finite volume options present

Courant Number mean: 6.98835e-07 max: 0.00055208
forces propellerTip:
    Not including porosity effects
forces DAS:
    Not including porosity effects

Starting time loop

Courant Number mean: 6.98835e-07 max: 0.00055208
deltaT = 1e-05
Time = 1e-05

AMI: Creating addressing and weights between 45146 source faces and 46392 target faces
AMI: Patch source sum(weights) min/max/average = 0.00638659, 2.10717, 1.00071
AMI: Patch target sum(weights) min/max/average = 0, 1.84552, 0.998866
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0302785, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.000126242, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0305505, No Iterations 1
[1] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[1] #1  Foam::sigFpe::sigHandler(int) at ??:?
[1] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[1] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[1] #5  ? at ??:?
[1] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7  ? at ??:?
[pouria:03168] *** Process received signal ***
[pouria:03168] Signal: Floating point exception (8)
[pouria:03168] Signal code:  (-6)
[pouria:03168] Failing at address: 0x3e800000c60
[pouria:03168] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f15358094b0]
[pouria:03168] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f1535809428]
[pouria:03168] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f15358094b0]
[pouria:03168] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE+0x114)[0x7f1536c46264]
[pouria:03168] [ 4] pimpleDyMFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_11dimensionedIdEERKS8_+0x1d7)[0x474ba7]
[pouria:03168] [ 5] pimpleDyMFoam[0x4259ac]
[pouria:03168] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f15357f4830]
[pouria:03168] [ 7] pimpleDyMFoam[0x4286c9]
[pouria:03168] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 3168 on node pouria exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

by the way. I am uploading a section of my project geometry. The purple block shows AMI region.

thank you.

[Tobi]

Hi,

for your moveDynamicMesh application you should call the following:

moveDynamicMesh -noFunctionObjects -checkAMI

After that, check if your mesh is correct moving.

[Taghi]

Thanks. It worked. But I still have problems with pimplyDyMFoam

[Tobi]

Is your AMI moving correctly and are the BC correct? You have a floating point exception (division by zero). Maybe you missed some

value uniform xxx;

in a BC. Sometimes it happen that this causes the problem. However, your information are not sufficient to give any other advice. The only point is, that the solver give up after you solve your pressure equation. You can also try to deactivate the momentum predictor but there are a lot of things that can cause your problem.

[Taghi]

Hi again. It took a long time but I think I figured were the problem is. The AMI is too close to Boundaries. Is there any way that I can rotate my propeller without AMI. MRF of SRF?

[Tobi]

I think not, but you can use ACMI and directly include the boundary into the moving zone.