CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

multiphase turbulance case floating error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 7, 2014, 09:41
Smile multiphase turbulance case floating error
  #1
Senior Member
 
harshawardhank
Join Date: Mar 2014
Posts: 209
Rep Power: 13
harsha_kulkarni is on a distinguished road
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
Attached Images
File Type: jpg 0_new.rpm.jpg (29.4 KB, 19 views)
harsha_kulkarni is offline   Reply With Quote

Old   October 9, 2014, 02:28
Default
  #2
Senior Member
 
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 109
Rep Power: 17
rcastilla is on a distinguished road
Check where is k diverging. It can be a mesh problem...
rcastilla is offline   Reply With Quote

Old   October 11, 2014, 04:36
Default
  #3
Senior Member
 
harshawardhank
Join Date: Mar 2014
Posts: 209
Rep Power: 13
harsha_kulkarni is on a distinguished road
but for singlePhase same model working for given k. when I'm including multiphase it is producing error.
harsha_kulkarni is offline   Reply With Quote

Old   February 18, 2016, 05:06
Default
  #4
Member
 
Sebastian W.
Join Date: Nov 2012
Location: Saxony, Germany
Posts: 43
Rep Power: 13
nero235 is on a distinguished road
Send a message via ICQ to nero235
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
eg.:
Code:
Gauss linear limited 0.5
Regards, Sebastian
nero235 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 10:07.