CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   pimpleDyMFoam failed in tutorial cases (2.2.x) (https://www.cfd-online.com/Forums/openfoam-solving/117343-pimpledymfoam-failed-tutorial-cases-2-2-x.html)

waku2005 May 6, 2013 22:05

pimpleDyMFoam failed in tutorial cases (2.2.x)
 
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. :confused:

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
  :


sharonyue May 7, 2013 00:02

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 May 7, 2013 02:08

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 May 13, 2013 04:28

Hi All

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

Regards
waku2005


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