multiphaseInterFoam for RAS turbulence model
 

I am trying to run the damBreak4phase case in RAS turbulence model, however, the RAS model is uninterpreted, and the kEpsilonCoeffs are not able to be read.
I also try to use multiphaseInterFoam to calculate my case with three phases, and the same problem is met.
The terminal prompts are attached.
Have a nice day!
Chiven

--------------------------------------------------------------------------------------------------------------
Create time
Create mesh for time = 0

Reading g
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
#0 _ZN4Foam5error10printStackERNS_7OstreamE-0x1466af0
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
#1 _ZN4Foam6sigFpe13sigFpeHandlerEi-0x140cae0
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 _ZNK4Foam14incompressible9RASModels33nutWallFuncti onFvPatchScalarField7calcNutEv-0x3267110
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libincompressibleRASModels.so"
#4 _ZN4Foam14incompressible9RASModels33nutWallFunctio nFvPatchScalarField12updateCoeffsEv-0x326aab0
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libincompressibleRASModels.so"
#5 Foam::fvPatchField<double>::evaluate(Foam::Pstream ::commsTypes) in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxIA64GccDPOpt/multiphaseInterFoam"
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxIA64GccDPOpt/multiphaseInterFoam"
#7 _ZN4Foam14incompressible9RASModels8kEpsilonC1ERKNS _14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldE NS_7volMeshEEERKNS3_IdNS_13fvsPatchFieldENS_11surf aceMeshEEERNS_14transportModelE-0x345ea40
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libincompressibleRASModels.so"
#8 _ZN4Foam14incompressible8RASModel31adddictionaryCo nstructorToTableINS0_9RASModels8kEpsilonEE3NewERKN S_14GeometricFieldINS_6VectorIdEENS_12fvPatchField ENS_7volMeshEEERKNS6_IdNS_13fvsPatchFieldENS_11sur faceMeshEEERNS_14transportModelE-0x3410a70
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libincompressibleRASModels.so"
#9 _ZN4Foam14incompressible8RASModel3NewERKNS_14Geome tricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volM eshEEERKNS2_IdNS_13fvsPatchFieldENS_11surfaceMeshE EERNS_14transportModelE-0x34d6870
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libincompressibleRASModels.so"
#10 _ZN4Foam14incompressible15turbulenceModel36addturb ulenceModelConstructorToTableINS0_8RASModelEE18New turbulenceModelERKNS_14GeometricFieldINS_6VectorId EENS_12fvPatchFieldENS_7volMeshEEERKNS5_IdNS_13fvs PatchFieldENS_11surfaceMeshEEERNS_14transportModel E-0x34c0310
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libincompressibleRASModels.so"
#11 _ZN4Foam14incompressible15turbulenceModel3NewERKNS _14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldE NS_7volMeshEEERKNS2_IdNS_13fvsPatchFieldENS_11surf aceMeshEEERNS_14transportModelE-0xe414a0
in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libincompressibleTurbulenceModel.so"
#12 main in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxIA64GccDPOpt/multiphaseInterFoam"
#13 __libc_start_main-0x10608e0
in "/lib/tls/libc.so.6.1"
#14 _start in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxIA64GccDPOpt/multiphaseInterFoam"
Floating exception

You are getting a "Floating exception" so it looks like you have not initialized the turbulence fields appropriately.

H

Hi, Dear Dr. Weller, thanks a lot for your friendly and precious comments. I am trying to improve the initialized turbulence fields. And I also try to use the laminar model to calculate my case, and when Time = 0.056, the similar error is met. Another comments?
yours sincerely,
Chiven

---------------------------------------------------------------------------------------------------------------------------------
[7] #0 _ZN4Foam5error10printStackERNS_7OstreamE-0x18f5aa0
[7] in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
[7] #1 _ZN4Foam6sigFpe13sigFpeHandlerEi-0x140cae0
[7] in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
[7] #2 Uninterpreted:
[7] #3 _ZNK4Foam10GAMGSolver13scalingFactorERNS_5FieldIdE ERKS2_S5_S5_-0x1cbddde
[7] in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
[7] #4 _ZNK4Foam10GAMGSolver13scalingFactorERNS_5FieldIdE ERKNS_9lduMatrixES3_RKNS_10FieldFieldIS1_dEERKNS_8 UPtrListIKNS_17lduInterfaceFieldEEERKS2_h-0x1cbd370
[7] in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
[7] #5 _ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9ldu Matrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_RNS1 _IS8_EESD_h-0x1cb6b00
[7] in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
[7] #6 _ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h-0x1cb58a0
[7] in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so"
[7] #7 _ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE-0x2a75610
[7] in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libfiniteVolume.so"
[7] #8 main in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxIA64GccDPOpt/multiphaseInterFoam"
[7] #9 __libc_start_main-0x14ef870
[7] in "/lib/tls/libc.so.6.1"
[7] #10 _start in "/home/g2/e090012/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxIA64GccDPOpt/multiphaseInterFoam"
[jaeri02:26701] *** Process received signal ***
[jaeri02:26701] Signal: Floating point exception (8)
[jaeri02:26701] Signal code: (-6)
[jaeri02:26701] Failing at address: 0x2d330000684d
[jaeri02:26701] [ 0] [0xa0000000000107e0]
[jaeri02:26701] [ 1] [0xa000000000010641]
[jaeri02:26701] [ 2] /lib/tls/libc.so.6.1(gsignal-0x14c8360) [0x2000000002705190]
[jaeri02:26701] [ 3] /home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so(_ZN4Foam6sigFpe13sigFpeHandlerEi-0x140cad0) [0x20000000022c1eb0]
[jaeri02:26701] [ 4] [0xa0000000000107e0]
[jaeri02:26701] [ 5] /home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so(_ZNK4Foam10GAMGSolver13scalingFacto rERNS_5FieldIdEERKS2_S5_S5_-0x1cbddde) [0x2000000001f0f692]
[jaeri02:26701] [ 6] /home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so(_ZNK4Foam10GAMGSolver13scalingFacto rERNS_5FieldIdEERKNS_9lduMatrixES3_RKNS_10FieldFie ldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEE RKS2_h-0x1cbd370) [0x2000000001f10110]
[jaeri02:26701] [ 7] /home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_RNS1_IS8_EESD_h-0x1cb6b00) [0x2000000001f16990]
[jaeri02:26701] [ 8] /home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h-0x1cb58a0) [0x2000000001f17c00]
[jaeri02:26701] [ 9] /home/g2/e090012/OpenFOAM/OpenFOAM-1.6/lib/linuxIA64GccDPOpt/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERKNS _10dictionaryE-0x2a75610) [0x2000000001157ea0]
[jaeri02:26701] [10] multiphaseInterFoam [0x40000000000c8560]
[jaeri02:26701] [11] /lib/tls/libc.so.6.1(__libc_start_main-0x14ef870) [0x20000000026ddc50]
[jaeri02:26701] [12] multiphaseInterFoam [0x4000000000023fa0]
[jaeri02:26701] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 7 with PID 26701 on node jaeri02 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

Hello, Dear Dr. Weller, now the liminar model is OK for the damBreak4phase case and my case. For the RASModel, I try to initialize the turbulence fields again, like the U, k, and epsilon, but fail and the same error is met. The BC setting is show in follows. The k and epsilon are obtained from
k=1.5*(5%*Ux)^2
epsilon=(Cmut^2*k^1.5)/(0.07*Dc), Dc is diameter
The k and epsilon setting way is OK for the single-phase turbulence calculation.
I also have tried to use boxtocell to set the fields of U, k,and epsilon.
Now, I have no idea about this problem. Especially I am amazing why the kEpsilonCoeffs{Cmu 0.09; C1 1.44; C2 1.92; sigmaEps 1.3;} can NOT be read.
Thank you very much!
Best regards,
Chiven

Hey Guys,

did you figured out how to run multiphaseInterFoam with RASModel?
I get almost the same message!

Thanks in advance

sarah

Hi All,

I got the same problem. Any updated solution to this turbulent initial condition for multiphase dambreak problems?

Hi All,

After many attempts, I found that by setting the boundaryfiled of the /nut to zeroGradient, the simulation run without any error.

Cheers

Hey mali,

nice, thanks a lot!
Do you use zeroGradient just for the walls, or for in- and Outlet as well?

Greetings

sarah

Hi,

Would you please someone explain to me what the following part of multiphaseInterFoam mean?

Cheers,

Elham