|
[Sponsors] |
August 27, 2009, 19:52 |
multiphaseInterFoam for RAS turbulence model
|
#1 |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
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 |
|
August 28, 2009, 02:39 |
|
#2 |
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 22 |
You are getting a "Floating exception" so it looks like you have not initialized the turbulence fields appropriately.
H |
|
August 28, 2009, 06:36 |
|
#3 |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
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). -------------------------------------------------------------------------- |
|
August 31, 2009, 02:43 |
|
#4 |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
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 Last edited by chiven; September 1, 2009 at 19:21. |
|
April 7, 2011, 07:05 |
|
#5 |
New Member
Sarah Köhler
Join Date: Jun 2010
Location: Leoben
Posts: 11
Rep Power: 16 |
Hey Guys,
did you figured out how to run multiphaseInterFoam with RASModel? I get almost the same message! Thanks in advance sarah |
|
August 22, 2011, 22:25 |
|
#6 |
Member
Join Date: Mar 2009
Location: adelaide, SA, Australia
Posts: 32
Rep Power: 17 |
Hi All,
I got the same problem. Any updated solution to this turbulent initial condition for multiphase dambreak problems?
__________________
mali |
|
August 23, 2011, 02:38 |
|
#7 |
Member
Join Date: Mar 2009
Location: adelaide, SA, Australia
Posts: 32
Rep Power: 17 |
Hi All,
After many attempts, I found that by setting the boundaryfiled of the /nut to zeroGradient, the simulation run without any error. Cheers
__________________
mali |
|
August 23, 2011, 02:52 |
|
#8 |
New Member
Sarah Köhler
Join Date: Jun 2010
Location: Leoben
Posts: 11
Rep Power: 16 |
Hey mali,
nice, thanks a lot! Do you use zeroGradient just for the walls, or for in- and Outlet as well? Greetings sarah |
|
December 6, 2017, 02:08 |
|
#9 |
Senior Member
Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 16 |
Hi,
Would you please someone explain to me what the following part of multiphaseInterFoam mean? Code:
forAllIter(PtrDictionary<phase>, phases_, iter) { phase& alpha = iter(); ... Elham |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
SimpleFoam case with SpalartAllmaras turbulence model implemented | nedved | OpenFOAM Running, Solving & CFD | 2 | November 30, 2014 22:43 |
Eul-Eul flow, k-e-kp-ep-Theta Turbulence model | us | FLUENT | 5 | April 5, 2011 02:29 |
KOmega Turbulence model from wwwopenFOAMWikinet | philippose | OpenFOAM Running, Solving & CFD | 30 | August 4, 2010 10:26 |
validating turbulence model on flow around a car | Pedro | CFX | 1 | February 20, 2008 16:32 |
SSG Reynolds Turbulence Model | Georges | CFX | 1 | February 28, 2007 16:15 |