CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   AMI interDyMFoam for mixer nu problem (http://www.cfd-online.com/Forums/openfoam-programming-development/119526-ami-interdymfoam-mixer-nu-problem.html)

danny123 June 19, 2013 05:23

AMI interDyMFoam for mixer nu problem
 
Hello,

I have altered the mixerAMI tutorial case to have an interface between 2 phases. The calculation seems starting off ok, and build the 1st time folder (0.1 s) result out of 5 s. Then the program abborts having the following error message:

HTML Code:

dw@PC83U:~/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2$ mpirun -np 6 interDyMFoam -parallel > log &
[1] 3945
dw@PC83U:~/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2$ [1]
[1]
[1] --> FOAM FATAL IO ERROR:
[1] keyword nu is undefined in dictionary "/home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor1/constant/transportProperties"
[1]
[1] file: [2]
[2]
[2] --> FOAM FATAL IO ERROR:
[2] keyword nu is undefined in dictionary "/home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor2/constant/transportProperties"
[2]
[2] file: /home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor2/constant/transportProperties from line 0 to line 0[4] [5]
[5]
[5] --> FOAM FATAL IO ERROR:
[5] keyword nu is undefined in dictionary "/home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor5/constant/transportProperties"
[5]
[5] file: /home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor5/constant/transportProperties from line 0 to line 0.
[5]
[5]    From function dictionary::lookupEntry(const word&, bool, bool) const/home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor1/constant/transportProperties from line 0 to line 0.
[1]
[1]    From function dictionary::lookupEntry(const word&, bool, bool) const
[1]    in file db/dictionary/dictionary.C at line 400.
[1]
FOAM parallel run exiting
[1]
.
[2]
[2]    From function dictionary::lookupEntry(const word&, bool, bool) const
[2]    in file db/dictionary/dictionary.C at line 400.
[2]
FOAM parallel run exiting
[2]
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 5 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

[5]    in file db/dictionary/dictionary.C at line 400.
[5]
FOAM parallel run exiting
[5]
[0]
[0]
[0] --> FOAM FATAL IO ERROR:
[0] keyword nu is undefined in dictionary "/home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor0/constant/transportProperties"
[0]
[0] file: /home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor0/constant/transportProperties from line 20 to line 43.
[0]
[0]    From function dictionary::lookupEntry(const word&, bool, bool) const
[0]    in file db/dictionary/dictionary.C at line 400.
[0]
FOAM parallel run exiting
[0]
[3]
[3]
[3] --> FOAM FATAL IO ERROR:
[3] keyword nu is undefined in dictionary "/home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor3/constant/transportProperties"
[3]
[3] file: /home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor3/constant/transportProperties from line 0 to line 0.
[3]
[3]    From function dictionary::lookupEntry(const word&, bool, bool) const
[3]    in file db/dictionary/dictionary.C at line 400.
[3]
FOAM parallel run exiting
[3]

[4]
[4] --> FOAM FATAL IO ERROR:
[4] keyword nu is undefined in dictionary "/home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor4/constant/transportProperties"
[4]
[4] file: /home/dw/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2/processor4/constant/transportProperties from line 0 to line 0.
[4]
[4]    From function dictionary::lookupEntry(const word&, bool, bool) const
[4]    in file db/dictionary/dictionary.C at line 400.
[4]
FOAM parallel run exiting
[4]
--------------------------------------------------------------------------
mpirun has exited due to process rank 5 with PID 3951 on
node PC83U exiting improperly. There are two reasons this could occur:

1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
for all processes to call "init". By rule, if one process calls "init",
then ALL processes must call "init" prior to termination.

2. this process called "init", but exited without calling "finalize".
By rule, all processes that call "init" MUST call "finalize" prior to
exiting or it will be considered an "abnormal termination"

This may have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------
[PC83U:03945] 5 more processes have sent help message help-mpi-api.txt / mpi-abort
[PC83U:03945] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

[1]+  Exit 1                  mpirun -np 6 interDyMFoam -parallel > log
dw@PC83U:~/OpenFOAM/dw-2.1.x/run/MooneyMeter_AMI_interDymFoam_2$

The last iteration step looks like this:

HTML Code:

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0.00718959 max: 0.0288944
deltaT = 0.00863822
Time = 0.1

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.1 transformation: ((0 0 0) (0.999945 (0 0 0.0104718)))
AMI: Creating addressing and weights between 3834 source faces and 3834 target faces
AMI: Patch source weights min/max/average = 1, 1.00701, 1.00045
AMI: Patch target weights min/max/average = 1, 1.00435, 1.00039
Execution time for mesh.update() = 0.14 s
MULES: Solving for alpha1
Phase-1 volume fraction = 0.531065  Min(alpha1) = -2.5463e-34  Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.531065  Min(alpha1) = -1.57193e-33  Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.531065  Min(alpha1) = -1.57193e-33  Max(alpha1) = 1
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.000580287, Final residual = 7.26327e-06, No Iterations 7
smoothSolver:  Solving for Uy, Initial residual = 0.000580396, Final residual = 7.1926e-06, No Iterations 7
smoothSolver:  Solving for Uz, Initial residual = 0.0348143, Final residual = 6.43603e-06, No Iterations 17
GAMG:  Solving for p_rgh, Initial residual = 4.63679e-07, Final residual = 9.73762e-09, No Iterations 11
time step continuity errors : sum local = 6.73903e-08, global = 8.68606e-09, cumulative = 2.30308e-06
GAMG:  Solving for p_rgh, Initial residual = 2.225e-08, Final residual = 4.5959e-09, No Iterations 6
time step continuity errors : sum local = 1.80979e-07, global = 9.84622e-09, cumulative = 2.31293e-06
GAMG:  Solving for p_rgh, Initial residual = 2.65858e-08, Final residual = 7.77777e-09, No Iterations 1
time step continuity errors : sum local = 1.27862e-07, global = 9.97309e-09, cumulative = 2.3229e-06
GAMG:  Solving for p_rgh, Initial residual = 1.39376e-08, Final residual = 7.85448e-09, No Iterations 2
time step continuity errors : sum local = 1.23004e-07, global = 9.88666e-09, cumulative = 2.33279e-06
GAMG:  Solving for p_rgh, Initial residual = 1.6827e-08, Final residual = 2.94924e-09, No Iterations 1
time step continuity errors : sum local = 4.24369e-08, global = 9.88312e-09, cumulative = 2.34267e-06
GAMG:  Solving for p_rgh, Initial residual = 5.29784e-09, Final residual = 5.29784e-09, No Iterations 0
time step continuity errors : sum local = 7.37101e-08, global = 9.88904e-09, cumulative = 2.35256e-06
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.000406186, Final residual = 9.17091e-06, No Iterations 8
smoothSolver:  Solving for Uy, Initial residual = 0.000406245, Final residual = 9.11417e-06, No Iterations 8
smoothSolver:  Solving for Uz, Initial residual = 0.03295, Final residual = 7.27708e-06, No Iterations 17
GAMG:  Solving for p_rgh, Initial residual = 1.20815e-07, Final residual = 6.02623e-09, No Iterations 8
time step continuity errors : sum local = 1.17466e-07, global = 1.1788e-08, cumulative = 2.36435e-06
GAMG:  Solving for p_rgh, Initial residual = 4.85723e-08, Final residual = 9.21026e-09, No Iterations 6
time step continuity errors : sum local = 1.65825e-07, global = 1.0649e-08, cumulative = 2.375e-06
GAMG:  Solving for p_rgh, Initial residual = 1.4502e-07, Final residual = 8.6203e-09, No Iterations 7
time step continuity errors : sum local = 2.90993e-08, global = 1.05608e-08, cumulative = 2.38556e-06
GAMG:  Solving for p_rgh, Initial residual = 4.74253e-08, Final residual = 5.46874e-09, No Iterations 1
time step continuity errors : sum local = 2.40905e-08, global = 1.05746e-08, cumulative = 2.39613e-06
GAMG:  Solving for p_rgh, Initial residual = 1.16193e-08, Final residual = 4.55811e-09, No Iterations 2
time step continuity errors : sum local = 1.99323e-08, global = 1.05617e-08, cumulative = 2.40669e-06
GAMG:  Solving for p_rgh, Initial residual = 9.88582e-09, Final residual = 9.88582e-09, No Iterations 0
time step continuity errors : sum local = 4.26416e-08, global = 1.05606e-08, cumulative = 2.41726e-06
PIMPLE: iteration 3
smoothSolver:  Solving for Ux, Initial residual = 7.31042e-05, Final residual = 8.09888e-06, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 7.3104e-05, Final residual = 8.05193e-06, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.0126386, Final residual = 8.20881e-06, No Iterations 14
GAMG:  Solving for p_rgh, Initial residual = 7.15945e-08, Final residual = 7.47356e-09, No Iterations 1
time step continuity errors : sum local = 3.32439e-08, global = 1.059e-08, cumulative = 2.42785e-06
GAMG:  Solving for p_rgh, Initial residual = 1.86575e-08, Final residual = 4.7532e-09, No Iterations 1
time step continuity errors : sum local = 2.10096e-08, global = 1.0583e-08, cumulative = 2.43843e-06
GAMG:  Solving for p_rgh, Initial residual = 8.49776e-09, Final residual = 8.49776e-09, No Iterations 0
time step continuity errors : sum local = 3.53645e-08, global = 1.05813e-08, cumulative = 2.44901e-06
GAMG:  Solving for p_rgh, Initial residual = 1.04412e-08, Final residual = 3.86596e-09, No Iterations 1
time step continuity errors : sum local = 1.63832e-08, global = 1.05833e-08, cumulative = 2.45959e-06
GAMG:  Solving for p_rgh, Initial residual = 4.58133e-09, Final residual = 4.58133e-09, No Iterations 0
time step continuity errors : sum local = 1.87374e-08, global = 1.05827e-08, cumulative = 2.47018e-06
GAMG:  Solving for p_rgh, Initial residual = 5.12872e-09, Final residual = 5.12872e-09, No Iterations 0
time step continuity errors : sum local = 2.09763e-08, global = 1.05827e-08, cumulative = 2.48076e-06
PIMPLE: iteration 4
smoothSolver:  Solving for Ux, Initial residual = 3.00653e-05, Final residual = 7.61561e-06, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 3.00715e-05, Final residual = 7.61024e-06, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00640266, Final residual = 8.44997e-06, No Iterations 13
GAMG:  Solving for p_rgh, Initial residual = 2.59555e-08, Final residual = 4.07093e-09, No Iterations 1
time step continuity errors : sum local = 1.67074e-08, global = 1.05998e-08, cumulative = 2.49136e-06
GAMG:  Solving for p_rgh, Initial residual = 8.68281e-09, Final residual = 8.68281e-09, No Iterations 0
time step continuity errors : sum local = 3.60341e-08, global = 1.06025e-08, cumulative = 2.50196e-06
GAMG:  Solving for p_rgh, Initial residual = 1.17972e-08, Final residual = 3.7096e-09, No Iterations 1
time step continuity errors : sum local = 1.5422e-08, global = 1.05886e-08, cumulative = 2.51255e-06
GAMG:  Solving for p_rgh, Initial residual = 5.39397e-09, Final residual = 5.39397e-09, No Iterations 0
time step continuity errors : sum local = 2.22894e-08, global = 1.05884e-08, cumulative = 2.52314e-06
GAMG:  Solving for p_rgh, Initial residual = 6.74929e-09, Final residual = 6.74929e-09, No Iterations 0
time step continuity errors : sum local = 2.789e-08, global = 1.05884e-08, cumulative = 2.53373e-06
GAMG:  Solving for p_rgh, Initial residual = 7.69072e-09, Final residual = 7.69072e-09, No Iterations 0
time step continuity errors : sum local = 3.17803e-08, global = 1.05884e-08, cumulative = 2.54431e-06
PIMPLE: iteration 5
smoothSolver:  Solving for Ux, Initial residual = 1.45307e-05, Final residual = 7.77353e-06, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1.4495e-05, Final residual = 7.72975e-06, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.0039205, Final residual = 8.98137e-06, No Iterations 12
GAMG:  Solving for p_rgh, Initial residual = 1.855e-08, Final residual = 5.7775e-09, No Iterations 1
time step continuity errors : sum local = 2.3897e-08, global = 1.05728e-08, cumulative = 2.55489e-06
GAMG:  Solving for p_rgh, Initial residual = 1.00594e-08, Final residual = 3.79035e-09, No Iterations 1
time step continuity errors : sum local = 1.5562e-08, global = 1.05829e-08, cumulative = 2.56547e-06
GAMG:  Solving for p_rgh, Initial residual = 5.14613e-09, Final residual = 5.14613e-09, No Iterations 0
time step continuity errors : sum local = 2.10885e-08, global = 1.05838e-08, cumulative = 2.57605e-06
GAMG:  Solving for p_rgh, Initial residual = 6.09124e-09, Final residual = 6.09124e-09, No Iterations 0
time step continuity errors : sum local = 2.49615e-08, global = 1.05838e-08, cumulative = 2.58664e-06
GAMG:  Solving for p_rgh, Initial residual = 6.82151e-09, Final residual = 6.82151e-09, No Iterations 0
time step continuity errors : sum local = 2.79541e-08, global = 1.05838e-08, cumulative = 2.59722e-06
GAMG:  Solving for p_rgh, Initial residual = 7.36735e-09, Final residual = 7.36735e-09, No Iterations 0
time step continuity errors : sum local = 3.01909e-08, global = 1.05838e-08, cumulative = 2.60781e-06
PIMPLE: iteration 6
smoothSolver:  Solving for Ux, Initial residual = 8.0169e-06, Final residual = 9.08257e-07, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 7.99963e-06, Final residual = 8.9875e-07, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.00224408, Final residual = 8.402e-07, No Iterations 16
GAMG:  Solving for p_rgh, Initial residual = 1.59378e-08, Final residual = 4.26623e-09, No Iterations 1
time step continuity errors : sum local = 1.74918e-08, global = 1.0607e-08, cumulative = 2.61841e-06
GAMG:  Solving for p_rgh, Initial residual = 8.24514e-09, Final residual = 8.24514e-09, No Iterations 0
time step continuity errors : sum local = 3.42347e-08, global = 1.061e-08, cumulative = 2.62902e-06
GAMG:  Solving for p_rgh, Initial residual = 1.10524e-08, Final residual = 3.48832e-09, No Iterations 1
time step continuity errors : sum local = 1.44895e-08, global = 1.0589e-08, cumulative = 2.63961e-06
GAMG:  Solving for p_rgh, Initial residual = 5.03605e-09, Final residual = 5.03605e-09, No Iterations 0
time step continuity errors : sum local = 2.07288e-08, global = 1.05877e-08, cumulative = 2.6502e-06
GAMG:  Solving for p_rgh, Initial residual = 6.05585e-09, Final residual = 6.05585e-09, No Iterations 0
time step continuity errors : sum local = 2.49263e-08, global = 1.05877e-08, cumulative = 2.66079e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 6.71681e-09, Final residual = 6.71681e-09, No Iterations 0
time step continuity errors : sum local = 2.76469e-08, global = 1.05877e-08, cumulative = 2.67137e-06
ExecutionTime = 26.87 s  ClockTime = 27 s

I am a little confused. 1st: there is no transport property file within processor0 to 5 folders. There is, of course, in the parent folder. This looks like this:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.com                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

phase1
{
transportModel  CrossPowerLaw;

nu              nu [ 0 2 -1 0 0 0 0 ] 1;

CrossPowerLawCoeffs
{
    nu0            nu0 [ 0 2 -1 0 0 0 0 ] 310.553906;
    nuInf          nuInf [ 0 2 -1 0 0 0 0 ] 0;
    m              m [ 0 0 1 0 0 0 0 ] 2.26191333;
    n              n [ 0 0 0 0 0 0 0 ] 0.33808236;
}

rho            rho [ 1 -3 0 0 0 0 0 ] 803.792583;

}

phase2
{
    transportModel  Newtonian;
    nu              nu [ 0 2 -1 0 0 0 0 ] 1e-5;
    rho            rho [ 1 -3 0 0 0 0 0 ] 1;
}

sigma          sigma [ 1 0 -2 0 0 0 0 ] 0.05;


// ************************************************************************* //

One hint: the folder "forces" is empty, no data written except the header. This is, where it hangs?

Thanks for the help.

Regards,

Daniel

sharonyue June 19, 2013 22:48

Some times it happens to me if I wanne write the "forces". if you find a solution I would like to know it. Thanks.

danny123 June 20, 2013 02:58

Hello Sharonyue,

Thanks for the reply. You gave me a hint. I tried calculating without forces and it works!

Well, I need forces anyway, since is how to verify the simulation by tests. Do you have an idea, meaning can you specify "sometimes" a little more in detail?

Regards,

Daniel

PS:

Well, there is another issue: I have created the level using the setfields utility, since creating it manually did not work. However, the alpha1 file looks different to what I found in this forum, meaning one of your posts. setFields removes all the patches except the stator and the titel "boundaryfield" is also removed. Recopying the text back does not work, meaning I cannot decompose the case. Calculating seems ok. so far.

musahossein August 27, 2013 17:37

keyword nu is undefined in dictionary "home/..../transportProperties"
 
I have the same problem as of today when I introduced the following code in controlDict:

forces
{
type forces;
functionObjectLibs ("libforces.so"); //Lib to load
patches (leftWall rightWall lowerWall);
pName p;
UName U;
rhoName rhoInf;
rhoInf 998.2; //Reference density for fluid
nuInf 1e-06;
CofR (0 0 0); //Origin for moment calculations
outputControl timeStep;
outputInterval 1;
}

The error message looks like this:

keyword nu is undefined in dictionary "/home/...../transportProperties"
[7]
[7] file: /home/...../transportProperties from line 0 to line 0.

The referenced line is from 0 to line 0. Any suggestions for a possible fix?

danny123 August 28, 2013 11:56

Hello musahossein,

I did a couple trials using the MixerVesselAMI using fill level and, so far did not find any particular pattern, why "Forces" does not work. It seems not to make a difference, whether the initial fill is in the stator Region, the rotor or across the AMI Interface.

Regards,

Daniel

musahossein August 28, 2013 16:27

You might want to review the following post as a way to solve the problem. I will do it in the next few hours. If you are successful, please let me know.

http://www.cfd-online.com/Forums/ope....html#poststop

musahossein August 30, 2013 17:10

keyword nu is undefined in dictionary "home/..../transportProperties"
 
Quote:

Originally Posted by musahossein (Post 448607)
You might want to review the following post as a way to solve the problem. I will do it in the next few hours. If you are successful, please let me know.

http://www.cfd-online.com/Forums/ope....html#poststop

Ok here is what I have understood after researching the forum extensively and trying it out for my self. If you are going to call "libforces.so" in controldict, then you need to do the following in order to avoid crashing (such has keyword nu undefined etc type of messages):

In the TurbulenceProperties file do the following:

set simulationType RASModel or LESModel;

if simulationType is set to laminar, comment it out, and out one of the choices.

In the RasProperties file change the settings to:

RASModel (or LESModel) laminar;
turbulence off;
printCoeffs on;

It appears that libforces.so in OpenFOAM 2.2.1 cannot resolve forces involving 2 fluids as discussed here: http://www.cfd-online.com/Forums/ope....html#poststop.

I am running sloshingtank2d using interdymfoam. I dont know if this will help you case but I hope it may highlight the fact that the problem is triggered by libforces.so and nothing else.

wyldckat August 31, 2013 13:24

Greetings to all!

I've answered to Musaddeque on this topic here: http://www.cfd-online.com/Forums/ope...tml#post449068 - post #13

In essence, the problem is that the forces function object was designed for working with single-phase solvers, not with multiphase solvers. Read the indicated post for more information.

Best regards,
Bruno

danny123 September 6, 2013 02:34

Thanks Bruno, thanks musahossein,

This is a big help. You are completly right. I will reply on post449068


for further questions if I have any.


Regards,

Daniel


All times are GMT -4. The time now is 06:42.