[kewsinger]

Hey there,
i try to simulate a bloodpump. Well, quite easy model of axial flow with rotating part using pimpleDyMFoam. The mesh works and rotates but the case doesn't run, producing following error:

++++++++++++++++++++++++++++++++++++++++

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : simpleFoam
Date   : Nov 20 2014
Time   : 10:21:16
Host   : "iz-lvp4-42.HS-Karlsruhe.DE"
PID    : 11582
Case   : /home/ADS/kuma1031/OpenFOAM/kuma1031-2.3.0/run/SHM_v19_rotate_merge/merge_v02/merged
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

Reading field p

AMI: Creating addressing and weights between 13777 source faces and 13652 target faces
AMI: Patch source sum(weights) min/max/average = 0.501124, 1.36466, 0.998511
AMI: Patch target sum(weights) min/max/average = 0, 2.62764, 0.985836
Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.1 average: 0.1
#0  Foam::error::printStack(Foam::Ostream&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib64/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvPatchField>(Foam::FieldField<Foam::fvPatchField, double>&, Foam::FieldField<Foam::fvPatchField, double> const&, Foam::FieldField<Foam::fvPatchField, double> const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#6  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::average<double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::average<double>(Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::average<double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::average<double>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10  Foam::bound(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensioned<double> const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#11  Foam::incompressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#12  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kEpsilon>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#13  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#14  
 in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleFoam"
#15  __libc_start_main in "/lib64/libc.so.6"
#16  
 in "/usr/lib64/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception

+++++++++++++++++++++++++++++++++++++++++++++

I just can't find my mistakes. Maybe anybody has some ideas to fix this!
The whole case is added as a dropbox link
https://www.dropbox.com/s/23nuj0pi7kubzzv/merged.zip

Thanks in advance,
Mario

[kewsinger]

maybe this is somehow connected to this error occurring while moveDynamicMesh command?!

Create time

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone rotating
Writing VTK files with weights of AMI patches.

Time = 0.0005
solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.0005 transformation: ((0 0 0) (0.995953 (0 0 0.0898785)))
AMI: Creating addressing and weights between 13777 source faces and 13652 target faces
AMI: Patch source sum(weights) min/max/average = 0.500287, 1.36597, 0.998531
AMI: Patch target sum(weights) min/max/average = 0, 1.58001, 0.985922
--> FOAM Warning : 
    From function solidBodyMotionFvMesh::update()
    in file solidBodyMotionFvMesh/solidBodyMotionFvMesh.C at line 203
    Did not find volVectorField U. Not updating U boundary conditions.
    Point usage OK.
    Upper triangular ordering OK.
    Topological cell zip-up check OK.
    Face vertices OK.
    Face-face connectivity OK.
    Mesh topology OK.
    Boundary openness (1.03605e-16 1.9261e-16 -1.84562e-16) OK.
    Max cell openness = 3.2531e-16 OK.
    Max aspect ratio = 21.1884 OK.
    Minimum face area = 5.11228e-09. Maximum face area = 1.96121e-05.  Face area magnitudes OK.
    Min volume = 6.98259e-13. Max volume = 5.90534e-08.  Total volume = 0.000474058.  Cell volumes OK.
    Mesh non-orthogonality Max: 64.947 average: 11.9615
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.74385 OK.
    Mesh geometry OK.
Mesh OK.
Calculating AMI weights between owner patch: AMI1 and neighbour patch: AMI2
ExecutionTime = 7.72 s  ClockTime = 8 s

what actually goes wrong here?

[alexeym]

Hi,

if your moveDynamicMesh somehow interacts with values of epsilon or k, then maybe. According to your log, error happens during construction of turbulence model, it is FPE, and it happens during division operation. Usually it's division by zero.

If you take a look at kEpsilon.C (constructor part):

...
    bound(k_, kMin_);
    bound(epsilon_, epsilonMin_);

    nut_ = Cmu_*sqr(k_)/epsilon_;
    nut_.correctBoundaryConditions();

    printCoeffs();
...

So first "bound" was performed, second seems to be skipped and your simulation halts with error either during calculation of nut_ or nut_.correctBoundaryConditions() call.

Also checkMesh in attached case:

Checking topology...
 ****Problem with boundary patch 0 named ader of type wall. The patch should start on face no 691406 and the patch specifies 718835.
Possibly consecutive patches have this same problem. Suppressing future warnings.
 ***Boundary definition is in error.

[kewsinger]

thanks Alexey for your input! So do you think the problem is the messed patch face no? What would be an idea to repair it?