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

TwoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 6, 2008, 01:53
Default Hi everyone I'm just starting
  #1
New Member
 
Sara Schairer
Join Date: Mar 2009
Posts: 26
Rep Power: 7
sara is on a distinguished road
Hi everyone
I'm just starting with OpenFOAM. I want to simulate oil mist filtration and I chose the solver twoPhaseEulerFoam.
After changing some of the files and filling in my parameters I run the calculation. Then I get this error message:

DILUPBiCG: Solving for alpha, Initial residual = 1, Final residual = 7.5185e-14, No Iterations 3
Dispersed phase volume fraction = 5.09482e-08 Min(alpha) = 0 Max(alpha) = 0.0005236
DILUPBiCG: Solving for alpha, Initial residual = 0.00235546, Final residual = 1.54916e-13, No Iterations 2
Dispersed phase volume fraction = 5.09482e-08 Min(alpha) = -3.26858e-18 Max(alpha) = 0.0005236
#0 Foam::error::printStack(Foam:stream&) in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb7f9d420]
#3 Foam::DICPreconditioner::calcReciprocalD(Foam::Fie ld<double>&, Foam::lduMatrix const&) in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::DICSmoother::DICSmoother(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<foam::field,> const&, Foam::FieldField<foam::field,> const&, Foam::UPtrList<foam::lduinterfacefield> const&) in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix::smoother::addsymMatrixConstructor ToTable<foam::dicsmoother>::Ne w(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<foam::field,> const&, Foam::FieldField<foam::field,> const&, Foam::UPtrList<foam::lduinterfacefield> const&) in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::lduMatrix::smoother::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<foam::field,> const&, Foam::FieldField<foam::field,> const&, Foam::UPtrList<foam::lduinterfacefield> const&, Foam::Istream&) in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::GAMGSolver::initVcycle(Foam::PtrList<foam::f ield<double> >&, Foam::PtrList<foam::field<double> >&, Foam::PtrList<foam::ldumatrix::smoother>&) const in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#8 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#9 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/schairer/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libfiniteVolume.so"
#10 main in "/home/schairer/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/twoPhaseEul erFoam"
#11 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#12 Foam::regIOobject::readIfModified() in "/home/schairer/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/twoPhaseEul erFoam"
Floating point exception

Does anybody know what this is about? I didn't change anything in the directories mentioned in there.
I'd be grateful for some help.
Cheers, Sara
sara is offline   Reply With Quote

Old   November 6, 2008, 02:01
Default Hi Sara, The listing says: .
  #2
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 644
Rep Power: 10
dmoroian is on a distinguished road
Hi Sara,
The listing says:
Quote:
...
#7 Foam::GAMGSolver::...
...
which I think complains about your pressure field that could not be resolved (GAMG is usually used for solving the pressure equation).
So it seems there is something wrong with your pressure settings (boundary conditions, initial conditions).

I hope this is helpful,
Dragos
dmoroian is offline   Reply With Quote

Old   November 6, 2008, 19:26
Default Hi Drago Maybe it was. I chan
  #3
New Member
 
Sara Schairer
Join Date: Mar 2009
Posts: 26
Rep Power: 7
sara is on a distinguished road
Hi Drago
Maybe it was. I changed alpha, too. Because it was set to 0 in the internal field and I thought that maybe somewhere during the calculation it tried to divide by zero.
In the p file I only changed the expression of the value from 1e+05 to 100000.
And I changed zeroGradient for the boundaries to fixedValue...
At least it works now, though I'm not sure if it's doing what I want it to do ;)
Anyway, thanks a lot for the quick answer.
Cheers, Sara
sara 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
TwoPhaseEulerFoam Documentation paul OpenFOAM Running, Solving & CFD 48 January 24, 2012 09:29
Bug in twoPhaseEulerFoam alberto OpenFOAM Bugs 2 May 20, 2008 21:25
TwoPhaseEulerFoam Bug alondono OpenFOAM Bugs 1 February 19, 2008 20:01
TwoPhaseEulerFOAM application hemph OpenFOAM Bugs 0 November 16, 2006 07:27
TwoPhaseEulerFoam newbee OpenFOAM 0 March 27, 2006 07:41


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