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

pimpleDyMFoam moving mesh error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 14, 2017, 14:28
Default pimpleDyMFoam moving mesh error
  #1
New Member
 
Pouria Taghikhani
Join Date: Sep 2016
Posts: 11
Rep Power: 2
Taghi is on a distinguished road
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:
Code:
 runApplication moveDynamicMesh -checkAMI
(as openfoam tutorials suggested)

but the following message apears:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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.
Attached Images
File Type: png section.png (4.3 KB, 18 views)

Last edited by Taghi; May 15, 2017 at 04:25.
Taghi is offline   Reply With Quote

Old   May 16, 2017, 18:07
Default
  #2
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,801
Blog Entries: 6
Rep Power: 32
Tobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

for your moveDynamicMesh application you should call the following:

Code:
moveDynamicMesh -noFunctionObjects -checkAMI
After that, check if your mesh is correct moving.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   May 17, 2017, 04:28
Default
  #3
New Member
 
Pouria Taghikhani
Join Date: Sep 2016
Posts: 11
Rep Power: 2
Taghi is on a distinguished road
Thanks. It worked. But I still have problems with pimplyDyMFoam
Taghi is offline   Reply With Quote

Old   May 17, 2017, 06:18
Default
  #4
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,801
Blog Entries: 6
Rep Power: 32
Tobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Is your AMI moving correctly and are the BC correct? You have a floating point exception (division by zero). Maybe you missed some

Code:
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.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   June 21, 2017, 08:00
Default
  #5
New Member
 
Pouria Taghikhani
Join Date: Sep 2016
Posts: 11
Rep Power: 2
Taghi is on a distinguished road
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?
Taghi is offline   Reply With Quote

Old   June 21, 2017, 08:06
Default
  #6
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,801
Blog Entries: 6
Rep Power: 32
Tobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
I think not, but you can use ACMI and directly include the boundary into the moving zone.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Reply

Tags
moving mesh, pimpledymfoam, propellers

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
GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM 300 October 29, 2014 19:00
How to install CGNS under windows xp? lzgwhy Main CFD Forum 1 January 11, 2011 19:44
groovyBC sega OpenFOAM Running, Solving & CFD 12 February 17, 2010 10:30
How to get the max value of the whole field waynezw0618 OpenFOAM Running, Solving & CFD 4 June 17, 2008 05:07
Compiling problems with hello worldC fw407 OpenFOAM Installation 21 January 6, 2008 18:38


All times are GMT -4. The time now is 00:13.