# Problems with compressible twoPhaseEulerFoam

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

 October 26, 2012, 08:36 Problems with compressible twoPhaseEulerFoam #1 Member   Sören Sander Join Date: Sep 2012 Location: Germany Posts: 43 Rep Power: 7 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 const&, Foam::UList const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::tmp > Foam:perator/(Foam::dimensioned const&, Foam::tmp > 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

 November 1, 2012, 04:58 #2 Member   Sören Sander Join Date: Sep 2012 Location: Germany Posts: 43 Rep Power: 7 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?

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Amit Mathur OpenFOAM 16 October 6, 2013 11:09 praveen OpenFOAM Programming & Development 0 April 11, 2011 00:00 liujmljm Main CFD Forum 1 November 7, 2010 18:54 maria teresa FLUENT 1 September 7, 2007 16:58 James FLUENT 2 June 20, 2007 04:22

All times are GMT -4. The time now is 22:41.