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

Problems with compressible twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 26, 2012, 08:36
Default Problems with compressible twoPhaseEulerFoam
  #1
Member
 
Sören Sander
Join Date: Sep 2012
Location: Germany
Posts: 33
Rep Power: 4
Sören Sander is on a distinguished road
Hi here,

I imported my own mesh to the testcase of the compressible tPEF, adopted the BC and started a case. Anyways it always throws the same error:

Courant Number mean: 6.63685e-09 max: 0.00185723
Max Ur Courant Number = 0.00185723
deltaT = 1.2e-05
Time = 1.2e-05

DILUPBiCG: Solving for alpha1, Initial residual = 1, Final residual = 4.67011e-17, No Iterations 1
Dispersed phase volume fraction = 1 Min(alpha1) = 0 Max(alpha1) = 1
DILUPBiCG: Solving for alpha1, Initial residual = 0.0621701, Final residual = 2.18343e-18, No Iterations 1
Dispersed phase volume fraction = 0.999999 Min(alpha1) = 0 Max(alpha1) = 1
DILUPBiCG: Solving for epsilon, Initial residual = 0.0171299, Final residual = 9.50906e-06, No Iterations 87
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 8.80289e-06, No Iterations 166
DILUPBiCG: Solving for T1, Initial residual = 1, Final residual = 2.2903e-14, No Iterations 2
DILUPBiCG: Solving for T2, Initial residual = 1, Final residual = 2.08691e-11, No Iterations 1
#0 Foam::error:rintStack(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::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam"
#5
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam"
Floating point exception
Sören Sander is offline   Reply With Quote

Old   November 1, 2012, 04:58
Default
  #2
Member
 
Sören Sander
Join Date: Sep 2012
Location: Germany
Posts: 33
Rep Power: 4
Sören Sander is on a distinguished road
I found that the problem was with the inlet boundary condition of U2. I tried to use a fixedValue of uniform (-1 0 0) instead of $internalField. If I change the BC to $interalField again and the solver runs without problems. However I want to specify an inlet velocity for phase 2. How can I do this?
Sören Sander is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Blow of compressible solver while using K-epsilon model in openfoam Amit Mathur OpenFOAM 16 October 6, 2013 11:09
Compressible code problems praveen OpenFOAM Programming & Development 0 April 11, 2011 00:00
Fixed grid methods for compressible viscous flow liujmljm Main CFD Forum 1 November 7, 2010 18:54
compressible flow maria teresa FLUENT 1 September 7, 2007 16:58
urgent help needed with 2d compressible flow James FLUENT 2 June 20, 2007 04:22


All times are GMT -4. The time now is 17:53.