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

pimpleDyMFoam failed in tutorial cases (2.2.x)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 6, 2013, 22:05
Default pimpleDyMFoam failed in tutorial cases (2.2.x)
  #1
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Dear All,

I'm using 2.2.x and recently experienced strange failures in using pimpleDyMFoam after "git pull".
In the four pimpleDyMFoam tutorials, here are my results.
1) propeller ... failed after 1 time step
2) mixerVesselAMI2D ... failed after 1 time step
3) movingCone ... success and results are OK
4) wingMotion .... success and results are OK

The 2 failed cases are both finished with Floating point exception after only 1 time step.

I think something is wrong in dynamic mesh process in particularly rotation of meshes.

Any information will be appreciated and thanks in advance !

waku2005


Code:
[waku@ensis10 mixerVesselAMI2D]$ pimpleDyMFoam 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.x-d2301c6b77ea
Exec   : pimpleDyMFoam
Date   : May 07 2013
Time   : 11:01:03
Host   : "ensis10"
PID    : 7243
Case   : /home/waku/OpenFOAM/waku-2.2.x/tutorials/incompressible/pimpleDyMFoam/mixerVesselAMI2D
nProcs : 1
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 rotor
Reading field p

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 96 source faces and 96 target faces
AMI: Patch source weights min/max/average = 1, 1, 1
AMI: Patch target weights min/max/average = 1, 1, 1
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Reading field rAU if present

Creating fintite volume options from fvOptions

Selecting finite volume options model type MRFSource
    Source: MRF1
    - applying source for all time
    - selecting cells using cellZone rotor
    - selected 1536 cell(s) with volume 0.000100459


PIMPLE: Operating solver in PISO mode


Starting time loop

Courant Number mean: 0 max: 0
deltaT = 0.00119048
Time = 0.00119048

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.00119048 transformation: ((0 0 0) (0.999999 (0 0 0.00124666)))
AMI: Creating addressing and weights between 96 source faces and 96 target faces
AMI: Patch source weights min/max/average = 0.0381209, 0.961918, 0.50002
AMI: Patch target weights min/max/average = 0, 1.00004, 0.50002
#0  Foam::error::printStack(Foam::Ostream&) in "/home/waku/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/waku/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  
 at sigaction.c:0
#3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/home/waku/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#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&) in "/home/waku/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
#5  
 in "/home/waku/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
#6  __libc_start_main in "/lib64/libc.so.6"
#7  
 in "/home/waku/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/pimpleDyMFoam"
Floating point exception
Code:
[waku@ensis10 OpenFOAM-2.2.x]$ git log
commit d2301c6b77ea0962b708dfd87628ec4818688f4d
Merge: 85a2812 39d7777
Author: Henry <Henry>
Date:   Mon May 6 21:36:45 2013 +0100

    Merge branch 'master' of github.com:OpenFOAM/OpenFOAM-2.2.x
  :
waku2005 is offline   Reply With Quote

Old   May 7, 2013, 00:02
Default
  #2
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 661
Rep Power: 8
sharonyue is on a distinguished road
Selecting finite volume options model type MRFSource
Source: MRF1
- applying source for all time
- selecting cells using cellZone rotor
- selected 1536 cell(s) with volume 0.000100459

Are you using pimpleDyMFoam?So in this two rotating case.and in dynamic mesh solver , it should not have MRF,rite? try to delete it.
sharonyue is offline   Reply With Quote

Old   May 7, 2013, 02:08
Default
  #3
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Hi

In the fvOptions, MRF zone is inactive and has no effects to AMI.
This fvOptions file is as provided by OF 2.2.x initial branch and never changed.

waku2005

// MRF source to enable comparison between AMI and MRF case
// Note: inactive by default
MRF1
{
type MRFSource;
active no; <<<<<<<<<<<<<<<
selectionMode cellZone;
cellZone rotor;

MRFSourceCoeffs
{
// Fixed patches (by default they 'move' with the MRF zone)
nonRotatingPatches ();

origin (0 0 0);
axis (0 0 1);
omega constant 104.72;
}
}
waku2005 is offline   Reply With Quote

Old   May 13, 2013, 04:28
Default
  #4
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Hi All

Recent update solved the problem and now get working fine in all tutorials.

Regards
waku2005
waku2005 is offline   Reply With Quote

Reply

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
Vessel tutorial problem hosseinhgf CFX 1 March 17, 2013 12:39
Exchanging the geomtry in the LES motorBike tutorial simpomann OpenFOAM 1 May 26, 2012 08:26
Tutorial Cases (Help) Rophys OpenFOAM Running, Solving & CFD 6 February 1, 2012 14:28
STAR-CD Tutorial shekhar aryal STAR-CD 4 March 22, 2010 04:25
Problem in running tutorial cases preetham OpenFOAM Installation 2 June 13, 2009 16:36


All times are GMT -4. The time now is 07:21.