[rafal]

Hi Hrv,

I have chosen a second option:
Extending the cylindrical sliding interface in vertical directions – this works very good (no errors, nice distribution of parameters). However I showed only a simplified case. On the top of my mixer i have something that will not "like" moving mesh topology so this solution is not suitable for me.

I concentrated on one sliding interface out of both sliders. This is something i am working on now and it is causing problems indeed.

Now I am diving in the code to understand better what is happening. I have still several things to check. If i manage to make it work i will post a solution. it may take me however 2-3 weeks (from Tuesday i'm on holiday in hot country not in Rainchester ).

Thanks for suggestions.

rafal

[hjasak]

Hehe, Rainchester: I know exactly what you mean.

To be honest with you, I need encouragement and/or company to get this finished and it's quite a hard job. I am trying to assemble a group of people who can either join in or contribute in other ways to get this work restarted. For all those out there who need sliding interfaces, complex topo changes and parallelism, now is the time to join in.

In any case, give me a shout when you come back to the real world and we can look for a way to take this further.

Have a good holiday,

Hrv

[auvinen]

Dear FOAM community,

I'm having a terrible time producing anything with icoDyMFoam. I cannot even reproduce the example cases delivered in this forum (kupplung, mixer3D). Okay, mixer2D ran fine. I do need some help for I'm from the Fortran community and a newbie to OpenFOAM.

1.) Case Kupplung runs until Time = 0.00573, but not beyond no matter what I do. The dramatic dying moment is captured here:

Mean and max Courant Numbers = 1.83329e-05 0.000287003
Time = 0.005729

ICCG: Solving for pcorr, Initial residual = 1, Final residual = 9.17333e-09, No Iterations 87
ICCG: Solving for pcorr, Initial residual = 0.0285899, Final residual = 6.106e-09, No Iterations 66
time step continuity errors : sum local = 1.40674e-20, global = 8.82553e-23, cumulative = 4.29905e-14
BICCG: Solving for Ux, Initial residual = 0.00324129, Final residual = 1.1653e-10, No Iterations 1
BICCG: Solving for Uy, Initial residual = 0.00317451, Final residual = 1.16643e-10, No Iterations 1
BICCG: Solving for Uz, Initial residual = 0.000524492, Final residual = 1.53121e-11, No Iterations 1
ICCG: Solving for p, Initial residual = 0.00444813, Final residual = 9.78411e-08, No Iterations 39
ICCG: Solving for p, Initial residual = 7.57075e-06, Final residual = 8.85904e-08, No Iterations 6
time step continuity errors : sum local = 3.72149e-14, global = -9.21735e-16, cumulative = 4.20687e-14
ICCG: Solving for p, Initial residual = 1.35521e-07, Final residual = 7.30238e-08, No Iterations 1
ICCG: Solving for p, Initial residual = 7.307e-08, Final residual = 7.307e-08, No Iterations 0
time step continuity errors : sum local = 3.06951e-14, global = -5.79509e-16, cumulative = 4.14892e-14
ICCG: Solving for p, Initial residual = 7.30763e-08, Final residual = 7.30763e-08, No Iterations 0
ICCG: Solving for p, Initial residual = 7.30763e-08, Final residual = 7.30763e-08, No Iterations 0
time step continuity errors : sum local = 3.06978e-14, global = -5.79044e-16, cumulative = 4.09102e-14
ICCG: Solving for p, Initial residual = 7.30763e-08, Final residual = 7.30763e-08, No Iterations 0
ICCG: Solving for p, Initial residual = 7.30763e-08, Final residual = 7.30763e-08, No Iterations 0
time step continuity errors : sum local = 3.06978e-14, global = -5.79043e-16, cumulative = 4.03312e-14
ExecutionTime = 27.52 s ClockTime = 28 s

Mean and max Courant Numbers = 1.83367e-05 0.000287053
Time = 0.00573

ICCG: Solving for pcorr: solution singularity
ICCG: Solving for pcorr: solution singularity
time step continuity errors : sum local = nan, global = nan, cumulative = nan
BICCG: Solving for Ux: solution singularity
BICCG: Solving for Uy: solution singularity
BICCG: Solving for Uz: solution singularity
ICCG: Solving for p: solution singularity
ICCG: Solving for p: solution singularity
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ICCG: Solving for p: solution singularity
ICCG: Solving for p: solution singularity
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ICCG: Solving for p: solution singularity
ICCG: Solving for p: solution singularity
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ICCG: Solving for p: solution singularity
ICCG: Solving for p: solution singularity
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ExecutionTime = 28.34 s ClockTime = 29 s

Mean and max Courant Numbers = nan 0.000287104
Time = 0.005731

Foam::error::printStack(Foam:stream&)
Foam::sigSegv::sigSegvHandler(int)
/lib64/libc.so.6 [0x30ca22ee80]
Foam::Field<double>::map(Foam::UList<double> const&, Foam::FieldMapper const&)
Foam::Field<double>::autoMap(Foam::FieldMapper const&)
void Foam::MapGeometricFields<double,>(Foam::fvMeshMapp er const&)
Foam::fvMesh::mapFields(Foam::mapPolyMesh const&)
Foam::fvMesh::updateMesh(Foam::mapPolyMesh const&)
Foam::polyTopoChanger::changeMesh()
Foam::mixerFvMesh::update()
icoDyMFoam [0x41497b]
__libc_start_main
__gxx_personality_v0
Segmentation fault

___________________________________________

2.) Mixer3D on the other hand meets its fatal end at the first time step. Here is the restart from latestTime:

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.3 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : icoDyMFoam . mixer3D
Root : /home/mauvinen/OpenFOAM/mauvinen-1.3/run/cases
Case : mixer3D
Nprocs : 1
Create time

Create mesh

Selecting dynamicFvMesh mixerFvMesh
void mixerFvMesh::addZonesAndModifiers() : Zones and modifiers already present. Skipping.
Mixer mesh:
origin: (0 0 0)
axis: (0 0 1)
rpm: 10
Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

ICCG: Solving for pcorr, Initial residual = 1, Final residual = 6.59728e-09, No Iterations 73
ICCG: Solving for pcorr, Initial residual = 0.000170241, Final residual = 8.33034e-09, No Iterations 46
time step continuity errors : sum local = 5.20285e-14, global = 4.66565e-22, cumulative = 4.66565e-22

Starting time loop

Mean and max Courant Numbers = 0.000225328 0.00300178
deltaT = 0.000595238
Time = 2.0006

ICCG: Solving for pcorr, Initial residual = 1, Final residual = 6.66178e-09, No Iterations 73
ICCG: Solving for pcorr, Initial residual = 0.000397392, Final residual = 9.09398e-09, No Iterations 48
time step continuity errors : sum local = 3.84912e-12, global = 7.8478e-22, cumulative = 1.25134e-21
BICCG: Solving for Ux, Initial residual = 0.000881017, Final residual = 3.12e-09, No Iterations 1
BICCG: Solving for Uy, Initial residual = 0.000881708, Final residual = 3.04676e-09, No Iterations 1
BICCG: Solving for Uz, Initial residual = 0.00101462, Final residual = 4.41322e-09, No Iterations 1
ICCG: Solving for p, Initial residual = 0.990379, Final residual = 9.22838e-07, No Iterations 67
ICCG: Solving for p, Initial residual = 0.0157626, Final residual = 6.4386e-07, No Iterations 48
time step continuity errors : sum local = 1.3987e-10, global = 2.30073e-21, cumulative = 3.55208e-21


--> FOAM FATAL ERROR : Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.

From function adjustPhi(surfaceScalarField& phi, const volVectorField& U,const volScalarField& p
in file cfdTools/general/adjustPhi/adjustPhi.C at line 108.

FOAM exiting

___________________________________________

Any help on this matter would be highly appreciated. I'm studying fairly complex pump simulations with our own code and I'd like to replicate some of our cases with OpenFOAM. If only I could get past these problems ...

- mikko

[hjasak]

Have a look through the forum: the last sourceforge release is riddled with bugs and you need a patched up version.

Hrv

[clo]

Hi all, can someone send me the mixer3D case that works with icoDyMFoam? As Mikko saied the case crash and I have the latest version...
thank you!

[deepblue17]

With the latest commit of the overlapGgi boundary condition in the OpenFOAM-1.6-ext git repository, this problem is now finally solved for me.

Thanks,
Oliver