songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ENH: Adding pressure work terms for compressibleInterFoam

Browse Source 
Removing temporary fix from compressibleInterDyFOAM to re-calculte Uf with morphing mesh
This commit is contained in:
sergio 2018-05-18 01:27:13 -07:00
parent 647e1c892d
commit e0912a058c
2 changed files with 6 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
applications/solvers/multiphase/compressibleInterFoam/TEqn.H
View File

@ -1,16 +1,15 @@
{
fvScalarMatrix TEqn
(
fvm::ddt(rho, T) + fvm::div(rhoPhi, T)
- fvm::Sp(contErr, T)
fvm::ddt(rho, T) + fvm::div(rhoPhi, T) - fvm::Sp(contErr, T)
- fvm::laplacian(turbulence.alphaEff(), T)
+ (
divU*p
+ fvc::ddt(rho, K) + fvc::div(rhoPhi, K)
(divU*p)() - contErr/rho*p
+ (fvc::ddt(rho, K) + fvc::div(rhoPhi, K))() - contErr*K
)
*(
alpha1/mixture.thermo1().Cv()
+ alpha2/mixture.thermo2().Cv()
alpha1()/mixture.thermo1().Cv()()
+ alpha2()/mixture.thermo2().Cv()()
)
==
fvOptions(rho, T)

6
applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
View File

@ -119,12 +119,8 @@ int main(int argc, char *argv[])
ghf = (g & mesh.Cf()) - ghRef;
}
if ((mesh.changing() && correctPhi) || mesh.topoChanging())
if ((mesh.changing() && correctPhi))
{
// Calculate absolute flux from the mapped surface velocity
// Note: temporary fix until mapped Uf is assessed
Uf = fvc::interpolate(U);
// Calculate absolute flux from the mapped surface velocity
phi = mesh.Sf() & Uf;