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

talk about the most suitable way to discretize diffusion term in UEqn

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 16, 2014, 04:30
Talking talk about the most suitable way to discretize diffusion term in UEqn
  #1
Member
 
xuhe-openfoam
Join Date: Aug 2013
Location: DaLian,china
Posts: 82
Rep Power: 3
bieshuxuhe is on a distinguished road
hi,
I often see two kinds of UEqn in "Foam of multiphase flow"(e.g. interFoam ,MRFInterFoam ......) ! these are :
the first is
Code:
   surfaceScalarField muEff
    (
        "muEff",
        mixture.muf()
      + fvc::interpolate(rho*turbulence->nut())
    );

    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(mixture.rhoPhi(), U)
      - fvm::laplacian(muEff, U)
      - (fvc::grad(U) & fvc::grad(muEff))
    //- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
    );
this kind of UEqn often used for porous or MRF ...

the second is
Code:
    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      + turbulence->divDevRhoReff(rho, U)
    );

to make divDevRhoReff() more clear ,look down (for kEpsilon ):

00195 tmp<fvVectorMatrix> kEpsilon::divDevRhoReff
00196 (
00197     const volScalarField& rho,
00198     volVectorField& U
00199 ) const
00200 {
00201     volScalarField muEff("muEff", rho*nuEff());
00202 
00203     return
00204     (
00205       - fvm::laplacian(muEff, U)
00206       - fvc::div(muEff*dev(T(fvc::grad(U))))
00207     );
00208 }
this kind of UEqn often used for interFoam, InterDyMFoam ...
I think the red part in the code is the main difference ! That is the diffution term ! could you talk about how to make a choice between the two kinds of UEqn ?
for example : if we take the first UEqn for InterDyMFoam , will it be unsuitable for some reason ?
bieshuxuhe is offline   Reply With Quote

Old   April 16, 2014, 10:10
Default
  #2
Member
 
xuhe-openfoam
Join Date: Aug 2013
Location: DaLian,china
Posts: 82
Rep Power: 3
bieshuxuhe is on a distinguished road
you could discuss whatever you want about the discretization of diffution term .
bieshuxuhe is offline   Reply With Quote

Reply

Tags
diffusion term

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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
Adding an explicit source term in the UEqn() Anne Lincke OpenFOAM Running, Solving & CFD 0 June 25, 2013 12:19
Source term problem in UEqn of simpleFoam fisch OpenFOAM Programming & Development 1 June 17, 2011 10:57
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 15:45
secondary diffusion term Steve Moore FLUENT 9 September 7, 2005 23:14


All times are GMT -4. The time now is 19:06.