|
[Sponsors] |
Different momentum equation in different versions of settlingFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 15, 2015, 14:09 |
Different momentum equation in different versions of settlingFoam
|
#1 |
New Member
Shenan
Join Date: Nov 2014
Posts: 5
Rep Power: 12 |
The momentum equation for settlingFoam is different in OpenFOAM 2.2.1 and FOAM-extend 3.1. The term
Code:
-(fvc::grad(U) & fvc::grad(muEff)) In OF 2.2, the momentum equation is: Code:
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(phi, U) + fvc::div ( (Alpha/(scalar(1.001) - Alpha))*(sqr(rhoc)/rho)*Vdj*Vdj, "div(phiVdj,Vdj)" ) - fvm::laplacian(muEff, U, "laplacian(muEff,U)") - (fvc::grad(U) & fvc::grad(muEff)) ); Code:
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(phi, U) + fvc::div ( (Alpha/(scalar(1.001) - Alpha))*(sqr(rhoc)/rho)*Vdj*Vdj, "div(phiVdj,Vdj)" ) - fvm::laplacian(mu, U, "laplacian(muEff,U)") ); In OF 2.3.1, settlingFoam has been replaced with driftFluxFoam. Its momentum equation is: Code:
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) + fvc::div(UdmModel.tauDm()) + turbulence->divDevRhoReff(U) == fvOptions(rho, U) ); Code:
tmp<fvVectorMatrix> kEpsilon::divDevRhoReff ( const volScalarField& rho, volVectorField& U ) const { volScalarField muEff("muEff", rho*nuEff()); return ( - fvm::laplacian(muEff, U) - fvc::div(muEff*dev(T(fvc::grad(U)))) ); } Code:
- fvm::laplacian(muEff, U, "laplacian(muEff,U)") - (fvc::grad(U) & fvc::grad(muEff)) Code:
turbulence->divDevRhoReff(U) If so, please can you tell me the relationship between Code:
(fvc::grad(U) & fvc::grad(muEff)) Code:
fvc::div(muEff*dev(T(fvc::grad(U)))) Thank you. Last edited by Shenan; April 15, 2015 at 14:15. Reason: ease of reading code |
|
Tags |
drift flux model, settlingfoam |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
solving the momentum equation in UEqn. | callahance | OpenFOAM | 2 | October 18, 2012 10:38 |
Coefficients discretized momentum equation | michujo | Main CFD Forum | 4 | June 20, 2012 02:33 |
Question on the discretization of momentum equation in icoFoam | MPJ | OpenFOAM | 3 | October 4, 2011 10:44 |
What is the Boussinesq Term in Momentum Equation | CFDtoy | Main CFD Forum | 0 | August 11, 2008 10:56 |
momentum equation | cfp | CFX | 0 | July 8, 2002 05:48 |