[waku2005]

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

[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

[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
  :

[sharonyue]

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.

[waku2005]

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]

Hi All

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

Regards
waku2005