|
[Sponsors] |
October 7, 2014, 09:41 |
multiphase turbulance case floating error
|
#1 |
Senior Member
harshawardhank
Join Date: Mar 2014
Posts: 209
Rep Power: 13 |
Hello friends,
I'm doing simulation of two opposed jet microextractor using OpenFOAM. from one jet water is comming and fro other toluene. I'm using multiphaseInterFOAM and k-epsilon for turbulance and I'm getting floating point error . I'm attaching the geometry. when I switch off the turbulance It woks. inlet velocity are given as 1m/s value of k-0.0015 value of epsilon -0.0454569 ================================================== ============= Courant Number mean: 67.1235 max: 2.55783e+09 Interface Courant Number mean: 3.33317 max: 3.18604e+08 Time = 5e-05 MULES: Solving for alpha.water water volume fraction, min, max = 7.64012e-06 -0.587206 1 MULES: Solving for alpha.toluene toluene volume fraction, min, max = 7.63143e-06 -0.00422416 1 MULES: Solving for alpha.air air volume fraction, min, max = -1.66329e+10 -1.31092e+34 1.63662e+34 Phase-sum volume fraction, min, max = -1.66329e+10 -1.31092e+34 1.63662e+34 MULES: Solving for alpha.water water volume fraction, min, max = 8.08753e-06 -50.491 61.9137 MULES: Solving for alpha.toluene toluene volume fraction, min, max = 8.07943e-06 -0.0364036 1 MULES: Solving for alpha.air air volume fraction, min, max = 8.37159e+16 -2.13494e+41 1.86964e+41 Phase-sum volume fraction, min, max = 8.37159e+16 -2.13494e+41 1.86964e+41 MULES: Solving for alpha.water water volume fraction, min, max = 8.53494e-06 -10980.1 17978.1 MULES: Solving for alpha.toluene toluene volume fraction, min, max = 8.52743e-06 -0.318669 1 MULES: Solving for alpha.air air volume fraction, min, max = 0 -1.55348e+49 1.02287e+49 Phase-sum volume fraction, min, max = 0 -1.55348e+49 1.02287e+49 MULES: Solving for alpha.water water volume fraction, min, max = 8.98242e-06 -2.66078e+07 2.61686e+07 MULES: Solving for alpha.toluene toluene volume fraction, min, max = 8.97543e-06 -3.42201 2.06948 MULES: Solving for alpha.air air volume fraction, min, max = 5.02697e+33 -6.86083e+57 4.84214e+57 Phase-sum volume fraction, min, max = 5.02697e+33 -6.86083e+57 4.84214e+57 GAMG: Solving for p_rgh, Initial residual = 6.89188e-33, Final residual = 6.89188e-33, No Iterations 0 time step continuity errors : sum local = 6.60636e+07, global = -0.0168154, cumulative = -0.0255993 GAMG: Solving for p_rgh, Initial residual = 7.30448e-33, Final residual = 7.30448e-33, No Iterations 0 time step continuity errors : sum local = 7.00186e+07, global = 0.0165987, cumulative = -0.00900067 smoothSolver: Solving for epsilon, Initial residual = 0.245263, Final residual = 1.77353e-12, No Iterations 1 bounding epsilon, min: -3.01345e+22 max: 1.8657e+32 average: 2.29564e+27 smoothSolver: Solving for k, Initial residual = 0.999858, Final residual = 7.32211e-06, No Iterations 1000 bounding k, min: -504.816 max: 6.00298e+22 average: 8.58603e+20 ExecutionTime = 649.02 s ClockTime = 674 s Courant Number mean: 5.34644e+07 max: 3.22865e+15 Interface Courant Number mean: 1.66442e+07 max: 3.22865e+15 Time = 6e-05 MULES: Solving for alpha.water water volume fraction, min, max = 9106.17 -1.36203e+16 1.325e+16 MULES: Solving for alpha.toluene toluene volume fraction, min, max = 9.42388e-06 -59600 167014 MULES: Solving for alpha.air air volume fraction, min, max = -4.31814e+42 -1.21498e+67 1.30926e+67 Phase-sum volume fraction, min, max = -4.31814e+42 -1.21498e+67 1.30926e+67 MULES: Solving for alpha.water water volume fraction, min, max = 5.0864e+08 -5.20075e+27 3.97036e+27 MULES: Solving for alpha.toluene toluene volume fraction, min, max = 9.87405e-06 -2.04296e+12 2.24217e+12 MULES: Solving for alpha.air air volume fraction, min, max = 4.63656e+51 -2.58722e+76 2.33644e+76 Phase-sum volume fraction, min, max = 4.63656e+51 -2.58722e+76 2.33644e+76 MULES: Solving for alpha.water water volume fraction, min, max = 5.88807e+14 -8.99561e+35 6.87842e+35 MULES: Solving for alpha.toluene toluene volume fraction, min, max = -0.00736009 -2.89893e+22 3.47649e+22 MULES: Solving for alpha.air air volume fraction, min, max = -1.99139e+61 -4.72101e+85 5.41195e+85 Phase-sum volume fraction, min, max = -1.99139e+61 -4.72101e+85 5.41195e+85 MULES: Solving for alpha.water water volume fraction, min, max = -1.80447e+26 -5.22978e+49 5.69986e+49 MULES: Solving for alpha.toluene toluene volume fraction, min, max = 1.12489e+21 -2.58103e+45 2.60982e+45 MULES: Solving for alpha.air air volume fraction, min, max = 0 -1.23016e+95 1.68637e+95 Phase-sum volume fraction, min, max = 0 -1.23016e+95 1.68637e+95 GAMG: Solving for p_rgh, Initial residual = 0.999984, Final residual = 8782.28, No Iterations 1000 time step continuity errors : sum local = 3.59764e+45, global = 8.32869e+31, cumulative = 8.32869e+31 GAMG: Solving for p_rgh, Initial residual = 1.77362e-66, Final residual = 1.77362e-66, No Iterations 0 time step continuity errors : sum local = 5.13781e+57, global = -2.39079e+41, cumulative = -2.39079e+41 [1] #0 [0] #0 [3] #0 [2] #0 Foam::error:rintStack(Foam::Ostream&)Foam::error :rintStack(Foam::Ostream&)Foam::error:rintStac k(Foam::Ostream&)Foam::error:rintStack(Foam::Ost ream&) at ??:? [0] #1 Foam::sigFpe::sigHandler(int) at ??:? [1] #1 Foam::sigFpe::sigHandler(int) at ??:? [3] #1 Foam::sigFpe::sigHandler(int) at ??:? [2] #1 Foam::sigFpe::sigHandler(int) at ??:? [2] #2 Uninterpreted: [2] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? [3] #2 Uninterpreted: [3] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? [0] #2 Uninterpreted: [0] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? [1] #2 Uninterpreted: [1] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? [2] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [0] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [1] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [3] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [0] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [3] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [2] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [1] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [2] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? [3] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? [0] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? [1] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? [3] #7 at ??:? [1] #7 at ??:? [0] #7 at ??:? [2] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::diction ary const&)Foam::fvMatrix<double>::solve(Foam::diction ary const&)Foam::fvMatrix<double>::solve(Foam::diction ary const&) at ??:? [0] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:? [1] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:? [2] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:? [3] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:? [1] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:? [0] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:? [2] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:? [0] #10 at ??:? [1] #10 at ??:? [3] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:? [2] #10 at ??:? [3] #10 ================================================== =============== can you tell me how to tackle this error for analysis purpose I also want velocity data which is also not obtained in log file |
|
October 9, 2014, 02:28 |
|
#2 |
Senior Member
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 109
Rep Power: 17 |
Check where is k diverging. It can be a mesh problem...
|
|
October 11, 2014, 04:36 |
|
#3 |
Senior Member
harshawardhank
Join Date: Mar 2014
Posts: 209
Rep Power: 13 |
but for singlePhase same model working for given k. when I'm including multiphase it is producing error.
|
|
February 18, 2016, 05:06 |
|
#4 |
Member
|
Late to the party but for multiphase flow you encounter a problem that you can't specify an constant y+ for your case (if U would be constant). The dimensionless wall distance y+ depends on the kinematic viscosity nu, which for multiphase flow is different in each phase. Say you have 3 phases: air, water and oil. You will have 3 different nu and therefore is y+ different for each phase, even if you have the same velocity in each phase.
This problem can be solved by using a limited scheme (instead of corrected) for the Code:
laplacianScheme, sndGradScheme Code:
Gauss linear limited 0.5 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[blockMesh] blockMesh with double grading. | spwater | OpenFOAM Meshing & Mesh Conversion | 92 | January 12, 2019 09:00 |
Compile calcMassFlowC | aurore | OpenFOAM Programming & Development | 13 | March 23, 2018 07:43 |
[OpenFOAM.org] Compile OF 2.3 on Mac OS X .... the patch | gschaider | OpenFOAM Installation | 225 | August 25, 2015 19:43 |
Errors in UDF | shashank312 | Fluent UDF and Scheme Programming | 6 | May 30, 2013 20:30 |
Free surface boudary conditions with SOLA-VOF | Fan | Main CFD Forum | 10 | September 9, 2006 12:24 |