CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   twoPhaseEulerFoam : bubble column modeling (http://www.cfd-online.com/Forums/openfoam-solving/103070-twophaseeulerfoam-bubble-column-modeling.html)

chittipo June 10, 2012 08:30

twoPhaseEulerFoam : bubble column modeling
 
Dear members,

I am trying to model the rectangular bubble column using twoPhaseEulerFoam. The column is filled with water and the air bubbles are injected from the bottom of the tank.

Can somebody helping for this case, I am getting the following error message

"-------------------------------------------------------------------------------
Selecting dragModel for phase a: SchillerNaumann
Selecting dragModel for phase b: SchillerNaumann
dragPhase is a
Selecting viscosityModel Syamlal
Selecting conductivityModel HrenyaSinclair
Selecting radialModel Gidaspow
Selecting granularPressureModel Lun
Selecting frictionalStressModel JohnsonJackson
Courant Number mean: 0 max: 0

PIMPLE: Operating solver in PISO mode


Starting time loop

Reading/calculating field UaMean

Reading/calculating field UbMean

Reading/calculating field alphaMean

Reading/calculating field pMean

fieldAverage: starting averaging at time 0

Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
Time = 0.002

DILUPBiCG: Solving for alpha, Initial residual = 0, Final residual = 0, No Iterations 0
Dispersed phase volume fraction = 0 Min(alpha) = 0 Max(alpha) = 1
DILUPBiCG: Solving for alpha, Initial residual = 0, Final residual = 0, No Iterations 0
Dispersed phase volume fraction = 0 Min(alpha) = 0 Max(alpha) = 1
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam"
#9 __libc_start_main in "/lib/libc.so.6"
#10
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam"
Floating point exception

-----------------------------------------------------------------------------------------------"

enoch June 10, 2012 15:15

why is Co zero?
 
Quote:

Originally Posted by chittipo (Post 365672)
Dear members,

I am trying to model the rectangular bubble column using twoPhaseEulerFoam. The column is filled with water and the air bubbles are injected from the bottom of the tank.

Can somebody helping for this case, I am getting the following error message

"-------------------------------------------------------------------------------
Selecting dragModel for phase a: SchillerNaumann
Selecting dragModel for phase b: SchillerNaumann
dragPhase is a
Selecting viscosityModel Syamlal
Selecting conductivityModel HrenyaSinclair
Selecting radialModel Gidaspow
Selecting granularPressureModel Lun
Selecting frictionalStressModel JohnsonJackson
Courant Number mean: 0 max: 0

PIMPLE: Operating solver in PISO mode


Starting time loop

Reading/calculating field UaMean

Reading/calculating field UbMean

Reading/calculating field alphaMean

Reading/calculating field pMean

fieldAverage: starting averaging at time 0

Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
Time = 0.002

DILUPBiCG: Solving for alpha, Initial residual = 0, Final residual = 0, No Iterations 0
Dispersed phase volume fraction = 0 Min(alpha) = 0 Max(alpha) = 1
DILUPBiCG: Solving for alpha, Initial residual = 0, Final residual = 0, No Iterations 0
Dispersed phase volume fraction = 0 Min(alpha) = 0 Max(alpha) = 1
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam"
#9 __libc_start_main in "/lib/libc.so.6"
#10
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam"
Floating point exception

-----------------------------------------------------------------------------------------------"

Firstly, I will check material properties and BCs to figure out why Co is zero.

chittipo June 11, 2012 06:12

It works now. Thanks for the reply.


All times are GMT -4. The time now is 18:56.