CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Rotor-stator case with non-overlapping zones (http://www.cfd-online.com/Forums/openfoam-solving/109382-rotor-stator-case-non-overlapping-zones.html)

jrrygg November 16, 2012 06:10

Rotor-stator case with non-overlapping zones
 
Hi all,

I am trying to run a case with a rotor-stator interface. I am using OpenFOAM 1.6-ext, and set up a case with cyclicGgi boundary conditions.

I tried to run it (with interDyMFoam), but it gives me the following error:

Code:

ubuntu@ubuntu:~/OpenFOAM/ubuntu-1.6-ext/run/MD5-OF-GGI-2$ interDyMFoam
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM Extend Project: Open source CFD        |
|  \\    /  O peration    | Version:  1.6-ext                              |
|  \\  /    A nd          | Web:      www.extend-project.de                |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 1.6-ext
Exec  : interDyMFoam
Date  : Nov 16 2012
Time  : 10:00:17
Host  : ubuntu
PID    : 5605
Case  : /home/ubuntu/OpenFOAM/ubuntu-1.6-ext/run/MD5-OF-GGI-2
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh mixerGgiFvMesh
void mixerGgiFvMesh::addZonesAndModifiers() : Zones and modifiers already present.  Skipping.
Mixer mesh:
    origin: (0 0 0)
    axis  : (0 1 0)
    rpm  : 13

Reading g
Reading field pd

Reading field alpha1

Reading field U

Reading/calculating face flux field phi

Initializing the GGI interpolator between master/shadow patches: statif/rotif

    From function void GGIInterpolation<MasterPatch, SlavePatch>::rescaleWeightingFactors() const
    in file /build/buildd/openfoam-1.6-ext-1.6.0/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/GGIInterpolationWeights.C at line 535
    Uncovered faces found.  On master: 981 on slave: 839
Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Killed

I have some zones that are not overlapping, as I have simplified my mesh to cover only 1/3 of 360 degrees. Is this what is causing it to crash? Does GGI support non-overlapping zones?

Picture of the setup:
https://dl.dropbox.com/u/2820596/Pel...bined-mesh.png


All times are GMT -4. The time now is 11:18.