View Issue Details

Jump to NotesJump to History
IDProjectCategoryView StatusDate SubmittedLast Update
0001344OpenFOAMBugpublic2014-07-09 23:102016-02-03 22:05
Reporteralundilong Assigned Tohenry  
PrioritynormalSeveritymajorReproducibilityalways
Status closedResolutionno change required 
PlatformLinuxOSUbuntuOS Version10.04
Summary0001344: redundant divU in solver of interPhaseChangeFoam
DescriptionI think this divU term at line of 105 is redundant because it is already accounted for in tphiAlphaCorr(), where both flux from convection term and compression term are included, at line of 103. My derivation is attached.
Please correct me if I am wrong. Thanks!
 
 98 MULES::explicitSolve
 99 (
100 geometricOneField(),
101 alpha1,
102 phi,
103 tphiAlphaCorr(),
104 vDotvmcAlphal,
105 (divU*alpha1 + vDotcAlphal)(),
106 1,
107 0
108 );

this is tphiAlphaCorr
==============================================
 47 tmp<surfaceScalarField> tphiAlphaCorr
 48 (
 49 fvc::flux
 50 (
 51 phi,
 52 alpha1,
 53 alphaScheme
 54 )
 55 + fvc::flux
 56 (
 57 -fvc::flux(-phir, alpha2, alpharScheme),
 58 alpha1,
 59 alpharScheme
 60 )
 61 );
TagsNo tags attached.

Activities

alundilong

2014-07-09 23:10

reporter  

questionOnAlphaEqn.pdf (110,998 bytes)

henry

2016-02-03 22:05

manager   ~0005912

The divU is necessary to ensure boundedness of the phase-fraction during phase-change.

Issue History

Date Modified Username Field Change
2014-07-09 23:10 alundilong New Issue
2014-07-09 23:10 alundilong File Added: questionOnAlphaEqn.pdf
2016-02-03 22:05 henry Note Added: 0005912
2016-02-03 22:05 henry Status new => closed
2016-02-03 22:05 henry Assigned To => henry
2016-02-03 22:05 henry Resolution open => no change required