[mo_na]

Hi guys,

here is a link to the case I am working on: http://www.cfd-online.com/Forums/ope...ixed-body.html

I thought it would be better to open a new thread for the problem I am facing now, because probably its unrelated to my other problem.

I had two mesh parts with AMI1 and AMI2 as cyclicAMI patches both defined as patches in the blockMeshDict. I didn't manage to run snappyHexMesh like this so I used

stitchMesh AMI1 AMI2

After I was able to use snappyHexMesh and added entries for AMI1 and AMI2 for all my files in the 0 folder. I used

    type            cyclicAMI;
    transform       rotational;
    rotationAxis    (0 1 0);
    rotationCentre  (0 0 0);

for everything except pointDisplacement, where I used "slip" (not sure if thats right).

When I try to run the solver (ihDyMFoam) then I get the following error:

Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: sixDoFRigidBodyMotion
Selecting sixDoFSolver CrankNicolson
Translational constraint tensor (0 0 0 0 0 0 0 0 0)
Rotational constraint tensor (0 0 0 0 1 0 0 0 0)

PIMPLE: Operating solver in PISO mode

Reading field porosityIndex

Porosity NOT activated

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian


--> FOAM FATAL ERROR: 
Not implemented

    From function cyclicAMI::coupledFvPatchField<Type>::snGrad()
    in file /home/thesis/wilhelm/OpenFOAM/OpenFOAM-3.0.1/src/finiteVolume/lnInclude/coupledFvPatchField.H at line 153.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::coupledFvPatchField<double>::snGrad() const at ??:?
#3  Foam::fv::gaussGrad<double>::correctBoundaryConditions(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) at ??:?
#4  Foam::fv::gaussGrad<double>::calcGrad(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const at ??:?
#5  Foam::fv::gradScheme<double>::grad(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const at ??:?
#6  Foam::tmp<Foam::GeometricField<Foam::outerProduct<Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) at ??:?
#7  Foam::interfaceProperties::calculateK() at ??:?
#8  Foam::interfaceProperties::interfaceProperties(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::IOdictionary const&) at ??:?
#9  Foam::immiscibleIncompressibleTwoPhaseMixture::immiscibleIncompressibleTwoPhaseMixture(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#10  ? at ??:?
#11  __libc_start_main in "/lib64/libc.so.6"
#12  ? at /home/abuild/rpmbuild/BUILD/glibc-2.18/csu/../sysdeps/x86_64/start.S:125

Am I doing something wrong or could this be a bug?
Thanks for your help everyone!

Cheers,
Mona

[mo_na]

Ok apparently it wasn't right to use stitchMesh in this case. I separated my case in 2 parts now and then used mergeMeshes after snapping in the geometry in one of them (I explained it in the other thread) .
Now I don't get that error message anymore.