zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

multiphase MRF solvers: added Coriolis contribution to the virtual-mass force

Browse Source
This commit is contained in:
Henry
2012-03-21 17:35:54 +00:00
parent 503eefc325
commit 78522ae3c6
10 changed files with 19 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/DDtU.H
View File

@ -3,9 +3,11 @@
fvc::ddt(U1) fvc::ddt(U1)
+ fvc::div(phi1, U1) + fvc::div(phi1, U1)
- fvc::div(phi1)*U1; - fvc::div(phi1)*U1;
mrfZones.addCoriolis(U1, DDtU1);
DDtU2 = DDtU2 =
fvc::ddt(U2) fvc::ddt(U2)
+ fvc::div(phi2, U2) + fvc::div(phi2, U2)
- fvc::div(phi2)*U2; - fvc::div(phi2)*U2;
mrfZones.addCoriolis(U2, DDtU2);
} }

4
applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/UEqns.H
View File

@ -45,7 +45,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
- fvm::Sp(dragCoeff/rho1, U1) - fvm::Sp(dragCoeff/rho1, U1)
- alpha1*alpha2/rho1*(liftForce - Cvm*rho2*DDtU2) - alpha1*alpha2/rho1*(liftForce - Cvm*rho2*DDtU2)
); );
mrfZones.addCoriolis(alpha1, U1Eqn); mrfZones.addCoriolis(alpha1*(1 + Cvm*rho2*alpha2/rho1), U1Eqn);
U1Eqn.relax(); U1Eqn.relax();
} }
@ -76,7 +76,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
- fvm::Sp(dragCoeff/rho2, U2) - fvm::Sp(dragCoeff/rho2, U2)
+ alpha1*alpha2/rho2*(liftForce + Cvm*rho2*DDtU1) + alpha1*alpha2/rho2*(liftForce + Cvm*rho2*DDtU1)
); );
mrfZones.addCoriolis(alpha2, U2Eqn); mrfZones.addCoriolis(alpha2*(1 + Cvm*rho2*alpha1/rho2), U2Eqn);
U2Eqn.relax(); U2Eqn.relax();
} }
} }

2
applications/solvers/multiphase/multiphaseEulerFoam/DDtU.H
View File

@ -6,4 +6,6 @@ forAllIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
fvc::ddt(phase.U()) fvc::ddt(phase.U())
+ fvc::div(phase.phi(), phase.U()) + fvc::div(phase.phi(), phase.U())
- fvc::div(phase.phi())*phase.U(); - fvc::div(phase.phi())*phase.U();
mrfZones.addCoriolis(phase.U(), phase.DDtU());
} }

6
applications/solvers/multiphase/multiphaseEulerFoam/UEqns.H
View File

@ -49,7 +49,11 @@ forAllIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
) )
) )
); );
mrfZones.addCoriolis(alpha, UEqns[phasei]); mrfZones.addCoriolis
(
alpha*(1 + (1/phase.rho())*fluid.Cvm(phase)),
UEqns[phasei]
);
UEqns[phasei].relax(); UEqns[phasei].relax();
phasei++; phasei++;

2
applications/solvers/multiphase/twoPhaseEulerFoam/DDtU.H
View File

@ -3,9 +3,11 @@
fvc::ddt(U1) fvc::ddt(U1)
+ fvc::div(phi1, U1) + fvc::div(phi1, U1)
- fvc::div(phi1)*U1; - fvc::div(phi1)*U1;
mrfZones.addCoriolis(U1, DDtU1);
DDtU2 = DDtU2 =
fvc::ddt(U2) fvc::ddt(U2)
+ fvc::div(phi2, U2) + fvc::div(phi2, U2)
- fvc::div(phi2)*U2; - fvc::div(phi2)*U2;
mrfZones.addCoriolis(U2, DDtU2);
} }

4
applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
View File

@ -53,7 +53,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
//+ alpha2/rho1*K*U2 // Explicit drag transfered to p-equation //+ alpha2/rho1*K*U2 // Explicit drag transfered to p-equation
- alpha2/rho1*(liftCoeff - Cvm*rho2*DDtU2) - alpha2/rho1*(liftCoeff - Cvm*rho2*DDtU2)
); );
mrfZones.addCoriolis(U1Eqn); mrfZones.addCoriolis(scalar(1) + Cvm*rho2*alpha2/rho1, U1Eqn);
U1Eqn.relax(); U1Eqn.relax();
} }
@ -93,7 +93,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
//+ alpha1/rho2*K*U1 // Explicit drag transfered to p-equation //+ alpha1/rho2*K*U1 // Explicit drag transfered to p-equation
+ alpha1/rho2*(liftCoeff + Cvm*rho2*DDtU1) + alpha1/rho2*(liftCoeff + Cvm*rho2*DDtU1)
); );
mrfZones.addCoriolis(U2Eqn); mrfZones.addCoriolis(scalar(1) + Cvm*rho2*alpha1/rho2, U2Eqn);
U2Eqn.relax(); U2Eqn.relax();
} }
} }

2
src/finiteVolume/cfdTools/general/MRF/MRFZone.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License

2
src/finiteVolume/cfdTools/general/MRF/MRFZone.H
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License

2
src/finiteVolume/cfdTools/general/MRF/MRFZones.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License

2
src/finiteVolume/cfdTools/general/MRF/MRFZones.H
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License