Skip to content
Snippets Groups Projects
Commit 6cc6a0e0 authored by andy's avatar andy
Browse files

ENH: Updated porous zone treatment in reactingParcelFoam and LTSReactingParcelFoam

parent 6ce8c61d
Branches
Tags
No related merge requests found
......@@ -64,7 +64,6 @@ int main(int argc, char *argv[])
#include "createRadiationModel.H"
#include "createClouds.H"
#include "createExplicitSources.H"
#include "createPorousZones.H"
#include "initContinuityErrs.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
......
fvVectorMatrix UEqn
(
//pZones.ddt(rho, U)
fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
......@@ -14,8 +13,6 @@
sources.constrain(UEqn);
pZones.addResistance(UEqn);
if (pimple.momentumPredictor())
{
solve(UEqn == -fvc::grad(p));
......
Info<< "Creating porous zones" << nl << endl;
IOporosityModelList pZones(mesh);
......@@ -65,7 +65,6 @@ int main(int argc, char *argv[])
#include "createRadiationModel.H"
#include "createClouds.H"
#include "createExplicitSources.H"
#include "createPorousZones.H"
#include "initContinuityErrs.H"
#include "readTimeControls.H"
#include "compressibleCourantNo.H"
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment