
[Sponsors] 
February 20, 2012, 06:27 
UEqn without pressure term

#1 
Member
Join Date: Jun 2011
Posts: 35
Rep Power: 7 
Hello,
I have a momentum transport equation which is not pressuredriven, in addition to the continuity equation. Instead it has gravity and drag terms as the driving force. I can include it into OpenFOAM as: div(phi,U) == gravity + drag However, I can not include the continuity equation with this one. The solver does not compile when I try fvm::div(phi) or fvc::div(phi) I know that in simple algorithm, the continuity equation is included in the pressure equation. But i don't have to solve for pressure, it is not anywhere in the equations. How can I do it? 

February 21, 2012, 08:25 

#2 
Member
Join Date: Jun 2011
Posts: 35
Rep Power: 7 
I need to include the equation shown in attachment, which is a mass conservation equation, where "a" is phase volume fraction in a multiphase flow.
I defined "a" as a volScalarField. When I do at as: Code:
fvScalarMatrix alphaEqn ( fvm::ddt(alpha) + fvm::div(alpha,U) ); solve(alphaEqn); error: no matching function for call to ‘div(Foam::volScalarField&, Foam::volVectorField&)’ I suppose, I can not use two volume fields as argument in fvm::div Then I switch to Code:
fvScalarMatrix alphaEqn ( fvm::ddt(alpha) + fvc::div(phi) ); solve(alphaEqn); This time solver converges, but when I solve my case, it gives floating point exception if my internal field value for "a" is different than the voundary condition value. 

February 21, 2012, 15:47 

#4 
Member
Join Date: Jun 2011
Posts: 35
Rep Power: 7 
I modified it as
linearInterpolate(alpha*U) & mesh.Sf() in order to get phase flux. I am not sure if this makes sense. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Pulsatile pressure inlet with pressure outlet  a.lynchy  FLUENT  3  March 23, 2012 14:45 
pressure transient term?  Romek  CFX  4  April 19, 2011 21:19 
Pressure work term in turbulent K.E. equation  lost.identity  Main CFD Forum  0  March 8, 2011 13:21 
UDF to define or adjust pressure???  engahmed  FLUENT  0  July 6, 2010 17:19 
Pressure enigma.  Sara  Main CFD Forum  22  July 30, 2009 04:46 