CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Hi,I have a question about momentum equation in interFoam

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

Like Tree1Likes
  • 1 Post By fumiya

Reply
 
LinkBack Thread Tools Display Modes
Old   January 11, 2013, 10:19
Question Hi,I have a question about momentum equation in interFoam
  #1
New Member
 
anonymous
Join Date: Apr 2010
Posts: 14
Rep Power: 7
shchao is on a distinguished road
Dear Foamers:
I have read Henrik Rusche's phd thesis.Now i meet a question in the momentum equation of interFoam. Mr. Rusche defined an 'H' operator AH=As-An*phi where 'As' means source vector.Does As include the item like -grad(p)-g*x*grad(rho)...?Or What included in As?Thank you very much!
shchao is offline   Reply With Quote

Old   January 13, 2013, 06:59
Default
  #2
Senior Member
 
Fumiya Nozaki
Join Date: Jun 2010
Location: Yokohama, Japan
Posts: 188
Rep Power: 8
fumiya is on a distinguished road
Take a look at the interFoam/UEqn.H:

Code:
    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      + turbulence->divDevRhoReff(rho, U)
     ==
        sources(rho, U)
    );
In the case of UEqn.H(), 'H' operator is constructed from the above
fvVectorMatrix UEqn. In the above definition of UEqn, there is no
terms corresponding to grad(p) and g*x*grad(rho), so these terms
are not included in the H operator in the interFoam solver.

Fumiya
fumiya is offline   Reply With Quote

Old   January 14, 2013, 08:47
Default
  #3
New Member
 
anonymous
Join Date: Apr 2010
Posts: 14
Rep Power: 7
shchao is on a distinguished road
Here's the UEqn.H in of 2.1.1

surfaceScalarField muEff
(
"muEff",
twoPhaseProperties.muf()
+ fvc::interpolate(rho*turbulence->nut())
);
fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muEff, U)
- (fvc::grad(U) & fvc::grad(muEff))
//- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U)) ) & mesh.Sf()))
);
UEqn.relax();
if (pimple.momentumPredictor())
{
solve
(
UEqn
==
fvc::reconstruct
(
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
) * mesh.magSf()
)
);
}

The H operator is constructed by :
fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muEff, U)
- (fvc::grad(U) & fvc::grad(muEff))
//- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U)) ) & mesh.Sf()))
);

Right?
shchao is offline   Reply With Quote

Old   January 15, 2013, 08:26
Default
  #4
Senior Member
 
Fumiya Nozaki
Join Date: Jun 2010
Location: Yokohama, Japan
Posts: 188
Rep Power: 8
fumiya is on a distinguished road
I think so.

Fumiya
shchao likes this.
fumiya is offline   Reply With Quote

Old   January 28, 2013, 22:27
Default
  #5
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 420
Rep Power: 15
santiagomarquezd will become famous soon enough
Hi, the H operator gives which is defined in Hrv thesis for the implicit terms used in the fvVectorMatrix declaration.

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd 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
solving the momentum equation in UEqn. callahance OpenFOAM 2 October 18, 2012 09:38
momentum equation question mauz Main CFD Forum 0 March 26, 2010 12:34
Momentum equation used by mixture model George FLUENT 0 May 19, 2005 13:07
Problem with momentum equation Francisco Main CFD Forum 4 August 5, 2004 20:04
Nonstaggered Grids and Momentum Interpolations S. Wang Main CFD Forum 3 March 2, 2001 14:44


All times are GMT -4. The time now is 15:30.