CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

An additional pressure term in two phase flow

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

LinkBack Thread Tools Display Modes
Old   April 10, 2012, 12:27
Default An additional pressure term in two phase flow
Jeong Kim
Join Date: Feb 2010
Posts: 42
Rep Power: 9
enoch is on a distinguished road
I've been struggling with the below problem, but it is still unsolved. I started from modiftication of pEqn.H in bubbleFoam. There is no comilation error, but when it runs, it blows off due to high Courant number.

The momentum equations for phase1 and phase 2 are defined,
=-grad(p)/rhoa - p*grad(alpha)/rhoa + ...
=-grad(p)/rhob - p*grad(beta)/rhob + ...

There are additional pressure terms in the momenturm euqations for each phase: p*grad(alpha)/rhoa, p*grad(beta)/rhob

The face fluxes for alpha phase is defined.
phia = phia* - 1/(rhoa*Aa).f*S.f*grad(p).f - 1/(rhoa*A).f*S.f*(grad(alpha)/alpha).f*p

phib = ...

where phia* = (Ha/Aa).f*S.f + phiDraga, phib*=...
phiDraga = (explict-term drag) + (g&S.f)/Aa.f, phiDragb=...

In the pEqn.H,
Coefficents, Dp and Dpp for the lagrangian and divergence are defined in line 42 and 54, respectively.

Solve fvm::laplacian(Dp, p) + fvm::div(Dpp, p)== fvc::div(phi)

Update phia, phib, Ua and Ub in line 75, 78, 86 and 91, respectively.

Did I make any mistakes when addint the additional pressure terms in pEqn.H?

Below is a screen when the code run a case file. source and case files are linked below too.

Courant Number mean: 1.18108e+17 max: 6.18938e+19
Time = 0.01
Max Ur Courant Number = 1.53306e+21
DILUPBiCG: Solving for alpha, Initial residual = 0.064977, Final residual = 0.306707, No Iterations 1001
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
#1 Foam::sigFpe::sigHandler(int) in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
#2 in "/lib/x86_64-linux-gnu/"
#3 Foam:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
#4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
#5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam:ILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
#8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/lib/"
in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/bin/bubbleFoamMod"
in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/bin/bubbleFoamMod"
#11 __libc_start_main in "/lib/x86_64-linux-gnu/"
in "/home/OpenFOAM/openfoam210/platforms/linux64GccDPOpt/bin/bubbleFoamMod"
Floating point exception
enoch is offline   Reply With Quote


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
total pressure change in a pipe flow Im FLUENT 5 February 23, 2017 09:55
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Source term of multiphase flow in CFX44 youngan CFX 0 September 28, 2003 22:27
pressure gradient term in low speed flow Atit Koonsrisuk Main CFD Forum 2 January 10, 2002 11:52
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 21:31

All times are GMT -4. The time now is 21:03.