From f4c7f38f0cf30f8b93032f4e255fa17e9a821f16 Mon Sep 17 00:00:00 2001 From: andy Date: Mon, 12 Nov 2012 14:45:59 +0000 Subject: [PATCH] ENH: Do not shut-down pyrolysis reactions in the presence of a liquid film --- .../reactingOneDim/reactingOneDim.C | 30 ------------------- 1 file changed, 30 deletions(-) diff --git a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C index 5c36e4ad53..21d2c5c0ec 100644 --- a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C +++ b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C @@ -684,36 +684,6 @@ void reactingOneDim::preEvolveRegion() { solidChemistry_->setCellReacting(cellI, true); } - - // De-activate reactions if pyrolysis region coupled to (valid) film - if (filmCoupled_) - { - const volScalarField& filmDelta = filmDeltaPtr_(); - - forAll(intCoupledPatchIDs_, i) - { - const label patchI = intCoupledPatchIDs_[i]; - const scalarField& filmDeltap = filmDelta.boundaryField()[patchI]; - - forAll(filmDeltap, faceI) - { - const scalar filmDelta0 = filmDeltap[faceI]; - if (filmDelta0 > reactionDeltaMin_) - { - const labelList& cells = boundaryFaceCells_[faceI]; - - // TODO: only limit cell adjacent to film? - //solidChemistry_->setCellNoReacting(cells[0]) - - // Propagate flag through 1-D region - forAll(cells, k) - { - solidChemistry_->setCellReacting(cells[k], false); - } - } - } - } - } }