From aaf8819c92066c6ecabb67c0be198c644f08be6c Mon Sep 17 00:00:00 2001 From: Henry Weller <http://cfd.direct> Date: Tue, 16 Aug 2016 11:30:17 +0100 Subject: [PATCH] tetPtI -> tetPti --- .../polyMeshTetDecomposition/tetIndices.C | 8 +- .../polyMeshTetDecomposition/tetIndices.H | 2 +- .../polyMeshTetDecomposition/tetIndicesI.H | 6 +- .../wallBoundedParticleTemplates.C | 5 +- .../clouds/Templates/DSMCCloud/DSMCCloud.C | 4 +- .../clouds/Templates/DSMCCloud/DSMCCloud.H | 2 +- .../parcels/Templates/DSMCParcel/DSMCParcel.H | 2 +- .../Templates/DSMCParcel/DSMCParcelI.H | 4 +- .../parcels/derived/dsmcParcel/dsmcParcel.C | 4 +- .../basic/indexedParticle/indexedParticle.H | 4 +- src/lagrangian/basic/particle/particle.C | 10 +- src/lagrangian/basic/particle/particle.H | 16 +-- src/lagrangian/basic/particle/particleI.H | 135 +++++++++--------- src/lagrangian/basic/particle/particleIO.C | 6 +- .../basic/particle/particleTemplates.C | 41 +++--- .../basic/passiveParticle/passiveParticle.H | 4 +- .../CollidingParcel/CollidingParcel.H | 4 +- .../CollidingParcel/CollidingParcelI.H | 8 +- .../KinematicParcel/KinematicParcel.H | 4 +- .../KinematicParcel/KinematicParcelI.H | 8 +- .../Templates/MPPICParcel/MPPICParcel.H | 4 +- .../Templates/MPPICParcel/MPPICParcelI.H | 8 +- .../ReactingMultiphaseParcel.H | 4 +- .../ReactingMultiphaseParcelI.H | 8 +- .../Templates/ReactingParcel/ReactingParcel.H | 4 +- .../ReactingParcel/ReactingParcelI.H | 8 +- .../Templates/ThermoParcel/ThermoParcel.H | 4 +- .../Templates/ThermoParcel/ThermoParcelI.H | 8 +- .../CellZoneInjection/CellZoneInjection.C | 4 +- .../CellZoneInjection/CellZoneInjection.H | 2 +- .../ConeInjection/ConeInjection.C | 4 +- .../ConeInjection/ConeInjection.H | 2 +- .../ConeNozzleInjection/ConeNozzleInjection.C | 14 +- .../ConeNozzleInjection/ConeNozzleInjection.H | 4 +- .../FieldActivatedInjection.C | 4 +- .../FieldActivatedInjection.H | 2 +- .../InflationInjection/InflationInjection.C | 4 +- .../InflationInjection/InflationInjection.H | 2 +- .../InjectionModel/InjectionModel.C | 36 +++-- .../InjectionModel/InjectionModel.H | 4 +- .../KinematicLookupTableInjection.C | 4 +- .../KinematicLookupTableInjection.H | 2 +- .../ManualInjection/ManualInjection.C | 4 +- .../ManualInjection/ManualInjection.H | 2 +- .../InjectionModel/NoInjection/NoInjection.H | 2 +- .../PatchFlowRateInjection.C | 4 +- .../PatchFlowRateInjection.H | 2 +- .../PatchInjection/PatchInjection.C | 4 +- .../PatchInjection/PatchInjection.H | 2 +- .../PatchInjection/patchInjectionBase.C | 8 +- .../PatchInjection/patchInjectionBase.H | 2 +- .../SurfaceFilmModel/SurfaceFilmModel.C | 4 +- .../ReactingLookupTableInjection.C | 4 +- .../ReactingLookupTableInjection.H | 2 +- .../ReactingMultiphaseLookupTableInjection.C | 4 +- .../ReactingMultiphaseLookupTableInjection.H | 2 +- .../ThermoLookupTableInjection.C | 4 +- .../ThermoLookupTableInjection.H | 2 +- .../molecule/molecule/molecule.H | 2 +- .../molecule/molecule/moleculeI.H | 4 +- src/lagrangian/solidParticle/solidParticle.H | 2 +- src/lagrangian/solidParticle/solidParticleI.H | 4 +- .../Templates/SprayParcel/SprayParcel.H | 4 +- .../Templates/SprayParcel/SprayParcelI.H | 8 +- 64 files changed, 250 insertions(+), 239 deletions(-) diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.C b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.C index f380512362f..8e0d4f3581a 100644 --- a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.C +++ b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.C @@ -34,7 +34,7 @@ Foam::tetIndices::tetIndices() faceBasePtI_(-1), facePtAI_(-1), facePtBI_(-1), - tetPtI_(-1) + tetPti_(-1) {} @@ -53,7 +53,7 @@ Foam::tetIndices::tetIndices faceBasePtI_(faceBasePtI), facePtAI_(facePtAI), facePtBI_(facePtBI), - tetPtI_(tetPtI) + tetPti_(tetPtI) {} @@ -70,7 +70,7 @@ Foam::tetIndices::tetIndices faceBasePtI_(-1), facePtAI_(-1), facePtBI_(-1), - tetPtI_(tetPtI) + tetPti_(tetPtI) { const faceList& pFaces = mesh.faces(); const labelList& pOwner = mesh.faceOwner(); @@ -81,7 +81,7 @@ Foam::tetIndices::tetIndices faceBasePtI_ = mesh.tetBasePtIs()[facei_]; - label facePtI = (tetPtI_ + faceBasePtI_) % f.size(); + label facePtI = (tetPti_ + faceBasePtI_) % f.size(); label otherFacePtI = f.fcIndex(facePtI); if (own) diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.H b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.H index 9c4bdcad5bc..0c80210e636 100644 --- a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.H +++ b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.H @@ -104,7 +104,7 @@ class tetIndices //- Point on the face, *relative to the base point*, which // characterises this tet on the face. - label tetPtI_; + label tetPti_; public: diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndicesI.H b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndicesI.H index 7a30eea0999..1024595e4f3 100644 --- a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndicesI.H +++ b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndicesI.H @@ -23,8 +23,6 @@ License \*---------------------------------------------------------------------------*/ -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // Foam::label Foam::tetIndices::cell() const @@ -59,7 +57,7 @@ Foam::label Foam::tetIndices::facePtB() const Foam::label Foam::tetIndices::tetPt() const { - return tetPtI_; + return tetPti_; } @@ -185,7 +183,7 @@ Foam::label& Foam::tetIndices::facePtB() Foam::label& Foam::tetIndices::tetPt() { - return tetPtI_; + return tetPti_; } diff --git a/src/functionObjects/field/wallBoundedStreamLine/wallBoundedParticleTemplates.C b/src/functionObjects/field/wallBoundedStreamLine/wallBoundedParticleTemplates.C index ddbe21109a7..7ee23cac58d 100644 --- a/src/functionObjects/field/wallBoundedStreamLine/wallBoundedParticleTemplates.C +++ b/src/functionObjects/field/wallBoundedStreamLine/wallBoundedParticleTemplates.C @@ -34,7 +34,6 @@ void Foam::wallBoundedParticle::patchInteraction const scalar trackFraction ) { -// typedef TrackData::CloudType cloudType; typedef typename TrackData::cloudType::particleType particleType; particleType& p = static_cast<particleType&>(*this); @@ -45,7 +44,7 @@ void Foam::wallBoundedParticle::patchInteraction label origFacei = facei_; label patchi = patch(facei_); - // No action taken for tetPtI_ for tetFacei_ here, handled by + // No action taken for tetPti_ for tetFacei_ here, handled by // patch interaction call or later during processor transfer. @@ -175,7 +174,7 @@ Foam::scalar Foam::wallBoundedParticle::trackToEdge // Check angle to nbrCell tet. Is it in the direction of the // endposition? I.e. since volume of nbr tet is positive the // tracking direction should be into the tet. - tetIndices nbrTi(nbrCelli, tetFacei_, tetPtI_, mesh_); + tetIndices nbrTi(nbrCelli, tetFacei_, tetPti_, mesh_); if ((nbrTi.faceTri(mesh_).normal() & (endPosition-position())) < 0) { // Change into nbrCell. No need to change tetFace, tetPt. diff --git a/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.C b/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.C index 02a900a722c..6da1c1aa8e8 100644 --- a/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.C +++ b/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.C @@ -468,7 +468,7 @@ void Foam::DSMCCloud<ParcelType>::addNewParcel const scalar Ei, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId ) { @@ -480,7 +480,7 @@ void Foam::DSMCCloud<ParcelType>::addNewParcel Ei, celli, tetFacei, - tetPtI, + tetPti, typeId ); diff --git a/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.H b/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.H index 6c4a936ac7a..d7ea2be9b41 100644 --- a/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.H +++ b/src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.H @@ -455,7 +455,7 @@ public: const scalar Ei, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId ); diff --git a/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcel.H b/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcel.H index c9b81509d74..0435546342e 100644 --- a/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcel.H +++ b/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcel.H @@ -183,7 +183,7 @@ public: const scalar Ei, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId ); diff --git a/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcelI.H b/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcelI.H index fb41a6cee33..5fb82c5cdad 100644 --- a/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcelI.H +++ b/src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcelI.H @@ -60,11 +60,11 @@ inline Foam::DSMCParcel<ParcelType>::DSMCParcel const scalar Ei, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId ) : - ParcelType(mesh, position, celli, tetFacei, tetPtI), + ParcelType(mesh, position, celli, tetFacei, tetPti), U_(U), Ei_(Ei), typeId_(typeId) diff --git a/src/lagrangian/DSMC/parcels/derived/dsmcParcel/dsmcParcel.C b/src/lagrangian/DSMC/parcels/derived/dsmcParcel/dsmcParcel.C index ea58c2f59b9..e0ec951fccd 100644 --- a/src/lagrangian/DSMC/parcels/derived/dsmcParcel/dsmcParcel.C +++ b/src/lagrangian/DSMC/parcels/derived/dsmcParcel/dsmcParcel.C @@ -44,7 +44,7 @@ Foam::dsmcParcel::dsmcParcel const scalar Ei, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId ) : @@ -56,7 +56,7 @@ Foam::dsmcParcel::dsmcParcel Ei, celli, tetFacei, - tetPtI, + tetPti, typeId ) {} diff --git a/src/lagrangian/basic/indexedParticle/indexedParticle.H b/src/lagrangian/basic/indexedParticle/indexedParticle.H index 9c3a017a24b..75ed0783e42 100644 --- a/src/lagrangian/basic/indexedParticle/indexedParticle.H +++ b/src/lagrangian/basic/indexedParticle/indexedParticle.H @@ -68,11 +68,11 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label index = 0 ) : - particle(mesh, position, celli, tetFacei, tetPtI), + particle(mesh, position, celli, tetFacei, tetPti), index_(index) {} diff --git a/src/lagrangian/basic/particle/particle.C b/src/lagrangian/basic/particle/particle.C index fa1fc78420f..88de1e7dda7 100644 --- a/src/lagrangian/basic/particle/particle.C +++ b/src/lagrangian/basic/particle/particle.C @@ -50,7 +50,7 @@ Foam::particle::particle const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : mesh_(mesh), @@ -59,7 +59,7 @@ Foam::particle::particle facei_(-1), stepFraction_(0.0), tetFacei_(tetFacei), - tetPtI_(tetPtI), + tetPti_(tetPti), origProc_(Pstream::myProcNo()), origId_(getNewParticleID()) {} @@ -79,7 +79,7 @@ Foam::particle::particle facei_(-1), stepFraction_(0.0), tetFacei_(-1), - tetPtI_(-1), + tetPti_(-1), origProc_(Pstream::myProcNo()), origId_(getNewParticleID()) { @@ -98,7 +98,7 @@ Foam::particle::particle(const particle& p) facei_(p.facei_), stepFraction_(p.stepFraction_), tetFacei_(p.tetFacei_), - tetPtI_(p.tetPtI_), + tetPti_(p.tetPti_), origProc_(p.origProc_), origId_(p.origId_) {} @@ -112,7 +112,7 @@ Foam::particle::particle(const particle& p, const polyMesh& mesh) facei_(p.facei_), stepFraction_(p.stepFraction_), tetFacei_(p.tetFacei_), - tetPtI_(p.tetPtI_), + tetPti_(p.tetPti_), origProc_(p.origProc_), origId_(p.origId_) {} diff --git a/src/lagrangian/basic/particle/particle.H b/src/lagrangian/basic/particle/particle.H index d6f1855edbd..7da868fa2f3 100644 --- a/src/lagrangian/basic/particle/particle.H +++ b/src/lagrangian/basic/particle/particle.H @@ -156,7 +156,7 @@ protected: //- Index of the point on the face that defines the decomposed // tet that the particle is in. Relative to the face base // point. - label tetPtI_; + label tetPti_; //- Originating processor id label origProc_; @@ -189,7 +189,7 @@ protected: const label tetPlaneBasePtI, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const scalar tol ) const; @@ -203,7 +203,7 @@ protected: const label tetPlaneBasePtI, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const scalar tol ) const; @@ -211,12 +211,12 @@ protected: inline void tetNeighbour(label triI); //- Cross the from the given face across the given edge of the - // given cell to find the resulting face and tetPtI + // given cell to find the resulting face and tetPti inline void crossEdgeConnectedFace ( const label& celli, label& tetFacei, - label& tetPtI, + label& tetPti, const edge& e ); @@ -316,7 +316,7 @@ public: DefinePropertyList ( "(Px Py Pz) celli facei stepFraction " - "tetFacei tetPtI origProc origId" + "tetFacei tetPti origProc origId" ); //- Cumulative particle counter - used to provode unique ID @@ -343,10 +343,10 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); - //- Construct from components, tetFacei_ and tetPtI_ are not + //- Construct from components, tetFacei_ and tetPti_ are not // supplied so they will be deduced by a search particle ( diff --git a/src/lagrangian/basic/particle/particleI.H b/src/lagrangian/basic/particle/particleI.H index a31aadd57ed..5df1a22e92e 100644 --- a/src/lagrangian/basic/particle/particleI.H +++ b/src/lagrangian/basic/particle/particleI.H @@ -53,7 +53,7 @@ inline void Foam::particle::findTris tetPlaneBasePtIs[i], celli_, tetFacei_, - tetPtI_, + tetPti_, tol ); @@ -74,7 +74,7 @@ inline Foam::scalar Foam::particle::tetLambda const label tetPlaneBasePtI, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const scalar tol ) const { @@ -91,7 +91,7 @@ inline Foam::scalar Foam::particle::tetLambda tetPlaneBasePtI, celli, tetFacei, - tetPtI, + tetPti, tol ); } @@ -146,7 +146,7 @@ inline Foam::scalar Foam::particle::movingTetLambda const label tetPlaneBasePtI, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const scalar tol ) const { @@ -169,7 +169,7 @@ inline Foam::scalar Foam::particle::movingTetLambda vector n0 = Zero; { - tetIndices tetIs(celli, tetFacei, tetPtI, mesh_); + tetIndices tetIs(celli, tetFacei, tetPti, mesh_); // Tet at timestep start tetPointRef tet00 = tetIs.oldTet(mesh_); @@ -337,7 +337,7 @@ inline void Foam::particle::tetNeighbour(label triI) << abort(FatalError); } - label facePtI = (tetPtI_ + tetBasePtI) % f.size(); + label facePtI = (tetPti_ + tetBasePtI) % f.size(); label otherFacePtI = f.fcIndex(facePtI); switch (triI) @@ -348,7 +348,7 @@ inline void Foam::particle::tetNeighbour(label triI) // neighbour cell over tetFacei // Modification of celli_ will happen by other indexing, - // tetFacei_ and tetPtI don't change. + // tetFacei_ and tetPti don't change. break; } @@ -358,7 +358,7 @@ inline void Foam::particle::tetNeighbour(label triI) ( celli_, tetFacei_, - tetPtI_, + tetPti_, Foam::edge(f[facePtI], f[otherFacePtI]) ); @@ -368,9 +368,9 @@ inline void Foam::particle::tetNeighbour(label triI) { if (own) { - if (tetPtI_ < f.size() - 2) + if (tetPti_ < f.size() - 2) { - tetPtI_ = f.fcIndex(tetPtI_); + tetPti_ = f.fcIndex(tetPti_); } else { @@ -378,16 +378,16 @@ inline void Foam::particle::tetNeighbour(label triI) ( celli_, tetFacei_, - tetPtI_, + tetPti_, Foam::edge(f[tetBasePtI], f[otherFacePtI]) ); } } else { - if (tetPtI_ > 1) + if (tetPti_ > 1) { - tetPtI_ = f.rcIndex(tetPtI_); + tetPti_ = f.rcIndex(tetPti_); } else { @@ -395,7 +395,7 @@ inline void Foam::particle::tetNeighbour(label triI) ( celli_, tetFacei_, - tetPtI_, + tetPti_, Foam::edge(f[tetBasePtI], f[facePtI]) ); } @@ -407,9 +407,9 @@ inline void Foam::particle::tetNeighbour(label triI) { if (own) { - if (tetPtI_ > 1) + if (tetPti_ > 1) { - tetPtI_ = f.rcIndex(tetPtI_); + tetPti_ = f.rcIndex(tetPti_); } else { @@ -417,16 +417,16 @@ inline void Foam::particle::tetNeighbour(label triI) ( celli_, tetFacei_, - tetPtI_, + tetPti_, Foam::edge(f[tetBasePtI], f[facePtI]) ); } } else { - if (tetPtI_ < f.size() - 2) + if (tetPti_ < f.size() - 2) { - tetPtI_ = f.fcIndex(tetPtI_); + tetPti_ = f.fcIndex(tetPti_); } else { @@ -434,7 +434,7 @@ inline void Foam::particle::tetNeighbour(label triI) ( celli_, tetFacei_, - tetPtI_, + tetPti_, Foam::edge(f[tetBasePtI], f[otherFacePtI]) ); } @@ -458,7 +458,7 @@ inline void Foam::particle::crossEdgeConnectedFace ( const label& celli, label& tetFacei, - label& tetPtI, + label& tetPti, const edge& e ) { @@ -542,17 +542,17 @@ inline void Foam::particle::crossEdgeConnectedFace { // The point is the base point, so this is first tet // in the face circulation - tetPtI = 1; + tetPti = 1; } else if (eIndex == otherFace.size() - 1) { // The point is the last before the base point, so // this is the last tet in the face circulation - tetPtI = otherFace.size() - 2; + tetPti = otherFace.size() - 2; } else { - tetPtI = eIndex; + tetPti = eIndex; } break; @@ -621,19 +621,19 @@ inline Foam::label& Foam::particle::tetFace() inline Foam::label Foam::particle::tetPt() const { - return tetPtI_; + return tetPti_; } inline Foam::label& Foam::particle::tetPt() { - return tetPtI_; + return tetPti_; } inline Foam::tetIndices Foam::particle::currentTetIndices() const { - return tetIndices(celli_, tetFacei_, tetPtI_, mesh_); + return tetIndices(celli_, tetFacei_, tetPti_, mesh_); } @@ -676,21 +676,22 @@ inline void Foam::particle::initCellFacePt() position_, celli_, tetFacei_, - tetPtI_ + tetPti_ ); - if (celli_ == -1) + if (debug && celli_ == -1) { - FatalErrorInFunction - << "cell, tetFace and tetPt search failure at position " - << position_ << abort(FatalError); + WarningInFunction + << "cell, tetFace and tetPt search failure for position " + << position_ + << endl; } } else { - mesh_.findTetFacePt(celli_, position_, tetFacei_, tetPtI_); + mesh_.findTetFacePt(celli_, position_, tetFacei_, tetPti_); - if (tetFacei_ == -1 || tetPtI_ == -1) + if (tetFacei_ == -1 || tetPti_ == -1) { label oldCelli = celli_; @@ -699,43 +700,45 @@ inline void Foam::particle::initCellFacePt() position_, celli_, tetFacei_, - tetPtI_ + tetPti_ ); - if (celli_ == -1 || tetFacei_ == -1 || tetPtI_ == -1) + if (celli_ == -1 || tetFacei_ == -1 || tetPti_ == -1) { - // The particle has entered this function with a cell - // number, but hasn't been able to find a cell to - // occupy. + // The particle has entered this function with a cell number, + // but hasn't been able to find a cell to occupy. if (!mesh_.pointInCellBB(position_, oldCelli, 0.1)) { - // If the position is not inside the (slightly - // extended) bound-box of the cell that it thought - // it should be in, then this is considered an - // error. - - FatalErrorInFunction - << "position " << position_ << nl - << " for requested cell " << oldCelli << nl - << " If this is a restart or " - "reconstruction/decomposition etc. it is likely that" - " the write precision is not sufficient.\n" - " Either increase 'writePrecision' or " - "set 'writeFormat' to 'binary'" - << abort(FatalError); + // If the position is not inside the (slightly extended) + // bound-box of the cell that it thought it should be in, + // then this is considered an error. + + if (debug) + { + WarningInFunction + << "position " << position_ + << " not in specified cell " << oldCelli + << " with centre " << mesh_.cellCentres()[oldCelli] + << endl; + } + + celli_ = -1; + tetFacei_ = -1; + tetPti_ = -1; + + return; } - // The position is in the (slightly extended) - // bound-box of the cell. This situation may arise - // because the face decomposition of the cell is not - // the same as when the particle acquired the cell - // index. For example, it has been read into a mesh - // that has made a different face base-point decision - // for a boundary face and now this particle is in a - // position that is not in the mesh. Gradually move - // the particle towards the centre of the cell that it - // thought that it was in. + // The position is in the (slightly extended) bound-box of the + // cell. This situation may arise because the face + // decomposition of the cell is not the same as when the + // particle acquired the cell index. For example, it has been + // read into a mesh that has made a different face base-point + // decision for a boundary face and now this particle is in a + // position that is not in the mesh. Gradually move the + // particle towards the centre of the cell that it thought that + // it was in. celli_ = oldCelli; @@ -756,7 +759,7 @@ inline void Foam::particle::initCellFacePt() celli_, newPosition, tetFacei_, - tetPtI_ + tetPti_ ); iterNo++; @@ -777,7 +780,7 @@ inline void Foam::particle::initCellFacePt() << " to " << newPosition << " in cell " << celli_ << " tetFace " << tetFacei_ - << " tetPt " << tetPtI_ << nl + << " tetPt " << tetPti_ << nl << " (A fraction of " << 1.0 - mag(cC - newPosition)/mag(cC - position_) << " of the distance to the cell centre)" @@ -797,7 +800,7 @@ inline void Foam::particle::initCellFacePt() << " Found" << " cell " << celli_ << " tetFace " << tetFacei_ - << " tetPt " << tetPtI_ << nl + << " tetPt " << tetPti_ << nl << " This is a different cell to that which was supplied" << " (" << oldCelli << ")." << nl << endl; diff --git a/src/lagrangian/basic/particle/particleIO.C b/src/lagrangian/basic/particle/particleIO.C index 1c4b2bf9ad6..6dbb45c7b5c 100644 --- a/src/lagrangian/basic/particle/particleIO.C +++ b/src/lagrangian/basic/particle/particleIO.C @@ -51,7 +51,7 @@ Foam::particle::particle(const polyMesh& mesh, Istream& is, bool readFields) facei_(-1), stepFraction_(0.0), tetFacei_(-1), - tetPtI_(-1), + tetPti_(-1), origProc_(Pstream::myProcNo()), origId_(-1) { @@ -64,7 +64,7 @@ Foam::particle::particle(const polyMesh& mesh, Istream& is, bool readFields) is >> facei_ >> stepFraction_ >> tetFacei_ - >> tetPtI_ + >> tetPti_ >> origProc_ >> origId_; } @@ -111,7 +111,7 @@ Foam::Ostream& Foam::operator<<(Ostream& os, const particle& p) << token::SPACE << p.facei_ << token::SPACE << p.stepFraction_ << token::SPACE << p.tetFacei_ - << token::SPACE << p.tetPtI_ + << token::SPACE << p.tetPti_ << token::SPACE << p.origProc_ << token::SPACE << p.origId_; } diff --git a/src/lagrangian/basic/particle/particleTemplates.C b/src/lagrangian/basic/particle/particleTemplates.C index 876f765debd..15e1340a271 100644 --- a/src/lagrangian/basic/particle/particleTemplates.C +++ b/src/lagrangian/basic/particle/particleTemplates.C @@ -88,7 +88,7 @@ void Foam::particle::correctAfterParallelTransfer tetFacei_ = facei_ + ppp.start(); // Faces either side of a coupled patch have matched base indices, - // tetPtI is specified relative to the base point, already and + // tetPti is specified relative to the base point, already and // opposite circulation directions by design, so if the vertices // are: // source: @@ -108,7 +108,7 @@ void Foam::particle::correctAfterParallelTransfer // This relationship can be verified for other points and sizes of // face. - tetPtI_ = mesh_.faces()[tetFacei_].size() - 1 - tetPtI_; + tetPti_ = mesh_.faces()[tetFacei_].size() - 1 - tetPti_; // Reset the face index for the next tracking operation if (stepFraction_ > (1.0 - SMALL)) @@ -226,7 +226,7 @@ Foam::scalar Foam::particle::trackToFace // Pout<< "stepFraction " << stepFraction_ << nl // << "celli " << celli_ << nl // << "tetFacei " << tetFacei_ << nl - // << "tetPtI " << tetPtI_ + // << "tetPti " << tetPti_ // << endl; scalar trackFraction = 0.0; @@ -245,10 +245,10 @@ Foam::scalar Foam::particle::trackToFace // + pA and pB are the remaining points on the face, such that // the circulation, {basePt, pA, pB} produces a positive // normal by the right-hand rule. pA and pB are chosen from - // tetPtI_ do accomplish this depending if the cell owns the - // face, tetPtI_ is the vertex that characterises the tet, and + // tetPti_ do accomplish this depending if the cell owns the + // face, tetPti_ is the vertex that characterises the tet, and // is the first vertex on the tet when circulating around the - // face. Therefore, the same tetPtI represents the same face + // face. Therefore, the same tetPti represents the same face // triangle for both the owner and neighbour cell. // // Each tet has its four triangles represented in the same order: @@ -309,7 +309,7 @@ Foam::scalar Foam::particle::trackToFace label basePtI = f[tetBasePtI]; - label facePtI = (tetPtI_ + tetBasePtI) % f.size(); + label facePtI = (tetPti_ + tetBasePtI) % f.size(); label otherFacePtI = f.fcIndex(facePtI); label fPtAI = -1; @@ -430,7 +430,7 @@ Foam::scalar Foam::particle::trackToFace tetPlaneBasePtIs[tI], celli_, tetFacei_, - tetPtI_, + tetPti_, lambdaDistanceTolerance ); @@ -457,7 +457,7 @@ Foam::scalar Foam::particle::trackToFace tetBasePtI, fPtAI, fPtBI, - tetPtI_ + tetPti_ ); } else if (triI > 0) @@ -470,7 +470,7 @@ Foam::scalar Foam::particle::trackToFace // << " " << celli_ // << " " << facei_ // << " " << tetFacei_ - // << " " << tetPtI_ + // << " " << tetPti_ // << " " << triI // << " " << lambdaMin // << " " << trackFraction @@ -479,7 +479,7 @@ Foam::scalar Foam::particle::trackToFace // Pout<< "# Tracking loop tet " // << origId_ << " " << origProc_<< nl // << "# face: " << tetFacei_ << nl - // << "# tetPtI: " << tetPtI_ << nl + // << "# tetPti: " << tetPti_ << nl // << "# tetBasePtI: " << mesh_.tetBasePtIs()[tetFacei_] << nl // << "# tet.mag(): " << tet.mag() << nl // << "# tet.quality(): " << tet.quality() @@ -555,7 +555,7 @@ Foam::scalar Foam::particle::trackToFace label origFacei = facei_; label patchi = patch(facei_); - // No action taken for tetPtI_ for tetFacei_ here, handled by + // No action taken for tetPti_ for tetFacei_ here, handled by // patch interaction call or later during processor transfer. if @@ -1000,8 +1000,8 @@ void Foam::particle::hitCyclicPatch tetFacei_ = facei_; - // See note in correctAfterParallelTransfer for tetPtI_ addressing. - tetPtI_ = mesh_.faces()[tetFacei_].size() - 1 - tetPtI_; + // See note in correctAfterParallelTransfer for tetPti_ addressing. + tetPti_ = mesh_.faces()[tetFacei_].size() - 1 - tetPti_; const cyclicPolyPatch& receiveCpp = cpp.neighbPatch(); label patchFacei = receiveCpp.whichFace(facei_); @@ -1053,10 +1053,15 @@ void Foam::particle::hitCyclicAMIPatch if (patchFacei < 0) { - FatalErrorInFunction + // If the patch face of the particle is not known assume that + // the particle is lost and to be deleted + td.keepParticle = false; + + WarningInFunction << "Particle lost across " << cyclicAMIPolyPatch::typeName << " patches " << cpp.name() << " and " << receiveCpp.name() - << " at position " << position_ << abort(FatalError); + << " at position " << position_ + << endl; } // Convert face index into global numbering @@ -1066,8 +1071,8 @@ void Foam::particle::hitCyclicAMIPatch tetFacei_ = facei_; - // See note in correctAfterParallelTransfer for tetPtI_ addressing. - tetPtI_ = mesh_.faces()[tetFacei_].size() - 1 - tetPtI_; + // See note in correctAfterParallelTransfer for tetPti_ addressing. + tetPti_ = mesh_.faces()[tetFacei_].size() - 1 - tetPti_; // Now the particle is on the receiving side diff --git a/src/lagrangian/basic/passiveParticle/passiveParticle.H b/src/lagrangian/basic/passiveParticle/passiveParticle.H index 572692566d7..1c8f7be0f7f 100644 --- a/src/lagrangian/basic/passiveParticle/passiveParticle.H +++ b/src/lagrangian/basic/passiveParticle/passiveParticle.H @@ -64,10 +64,10 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - particle(mesh, position, celli, tetFacei, tetPtI) + particle(mesh, position, celli, tetFacei, tetPti) {} //- Construct from components, with searching for tetFace and diff --git a/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.H b/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.H index b581ef6eb9d..113f4ac51bd 100644 --- a/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcel.H @@ -171,7 +171,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); //- Construct from components @@ -181,7 +181,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcelI.H b/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcelI.H index a08b546441c..c40c4e5d54e 100644 --- a/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcelI.H +++ b/src/lagrangian/intermediate/parcels/Templates/CollidingParcel/CollidingParcelI.H @@ -66,10 +66,10 @@ inline Foam::CollidingParcel<ParcelType>::CollidingParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - ParcelType(owner, position, celli, tetFacei, tetPtI), + ParcelType(owner, position, celli, tetFacei, tetPti), f_(Zero), angularMomentum_(Zero), torque_(Zero), @@ -84,7 +84,7 @@ inline Foam::CollidingParcel<ParcelType>::CollidingParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, @@ -102,7 +102,7 @@ inline Foam::CollidingParcel<ParcelType>::CollidingParcel position, celli, tetFacei, - tetPtI, + tetPti, typeId, nParticle0, d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H index f5e57a39282..f3ebe37a843 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H @@ -328,7 +328,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); //- Construct from components @@ -338,7 +338,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H index 5f71e7b655d..42e6d144090 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H @@ -76,10 +76,10 @@ inline Foam::KinematicParcel<ParcelType>::KinematicParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - ParcelType(owner, position, celli, tetFacei, tetPtI), + ParcelType(owner, position, celli, tetFacei, tetPti), active_(true), typeId_(-1), nParticle_(0), @@ -103,7 +103,7 @@ inline Foam::KinematicParcel<ParcelType>::KinematicParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, @@ -112,7 +112,7 @@ inline Foam::KinematicParcel<ParcelType>::KinematicParcel const constantProperties& constProps ) : - ParcelType(owner, position, celli, tetFacei, tetPtI), + ParcelType(owner, position, celli, tetFacei, tetPti), active_(true), typeId_(typeId), nParticle_(nParticle0), diff --git a/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.H b/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.H index 313bd60025c..c02d97e1e0c 100644 --- a/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcel.H @@ -189,7 +189,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); //- Construct from components @@ -199,7 +199,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcelI.H b/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcelI.H index 3ee0b680544..1af5110019e 100644 --- a/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcelI.H +++ b/src/lagrangian/intermediate/parcels/Templates/MPPICParcel/MPPICParcelI.H @@ -32,10 +32,10 @@ inline Foam::MPPICParcel<ParcelType>::MPPICParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - ParcelType(owner, position, celli, tetFacei, tetPtI), + ParcelType(owner, position, celli, tetFacei, tetPti), UCorrect_(Zero) {} @@ -47,7 +47,7 @@ inline Foam::MPPICParcel<ParcelType>::MPPICParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, @@ -63,7 +63,7 @@ inline Foam::MPPICParcel<ParcelType>::MPPICParcel position, celli, tetFacei, - tetPtI, + tetPti, typeId, nParticle0, d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.H b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.H index 2f76f9ef7cb..417f085e713 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.H @@ -278,7 +278,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); @@ -289,7 +289,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcelI.H b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcelI.H index 005cafb6322..c64e29b262f 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcelI.H +++ b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcelI.H @@ -71,10 +71,10 @@ inline Foam::ReactingMultiphaseParcel<ParcelType>::ReactingMultiphaseParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - ParcelType(mesh, position, celli, tetFacei, tetPtI), + ParcelType(mesh, position, celli, tetFacei, tetPti), YGas_(0), YLiquid_(0), YSolid_(0), @@ -89,7 +89,7 @@ inline Foam::ReactingMultiphaseParcel<ParcelType>::ReactingMultiphaseParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, @@ -111,7 +111,7 @@ inline Foam::ReactingMultiphaseParcel<ParcelType>::ReactingMultiphaseParcel position, celli, tetFacei, - tetPtI, + tetPti, typeId, nParticle0, d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H index d48fa45e944..81db6586f3c 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H @@ -232,7 +232,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); //- Construct from components @@ -242,7 +242,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcelI.H b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcelI.H index 78fea81bc53..a6359c40062 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcelI.H +++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcelI.H @@ -66,10 +66,10 @@ inline Foam::ReactingParcel<ParcelType>::ReactingParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - ParcelType(mesh, position, celli, tetFacei, tetPtI), + ParcelType(mesh, position, celli, tetFacei, tetPti), mass0_(0.0), Y_(0), pc_(0.0) @@ -83,7 +83,7 @@ inline Foam::ReactingParcel<ParcelType>::ReactingParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, @@ -102,7 +102,7 @@ inline Foam::ReactingParcel<ParcelType>::ReactingParcel position, celli, tetFacei, - tetPtI, + tetPti, typeId, nParticle0, d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H index fd8af64a7b8..84b6e9c5e66 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H @@ -288,7 +288,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); //- Construct from components @@ -298,7 +298,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcelI.H b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcelI.H index b4e80ad19a9..9b2f365b447 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcelI.H +++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcelI.H @@ -77,10 +77,10 @@ inline Foam::ThermoParcel<ParcelType>::ThermoParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - ParcelType(mesh, position, celli, tetFacei, tetPtI), + ParcelType(mesh, position, celli, tetFacei, tetPti), T_(0.0), Cp_(0.0), Tc_(0.0), @@ -95,7 +95,7 @@ inline Foam::ThermoParcel<ParcelType>::ThermoParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, @@ -113,7 +113,7 @@ inline Foam::ThermoParcel<ParcelType>::ThermoParcel position, celli, tetFacei, - tetPtI, + tetPti, typeId, nParticle0, d0, diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.C index 8336f8c88fa..07ec1e9b26d 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.C @@ -320,13 +320,13 @@ void Foam::CellZoneInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { position = positions_[parcelI]; cellOwner = injectorCells_[parcelI]; tetFacei = injectorTetFaces_[parcelI]; - tetPtI = injectorTetPts_[parcelI]; + tetPti = injectorTetPts_[parcelI]; } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.H index 9213037ee9f..be21d3f5ece 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/CellZoneInjection/CellZoneInjection.H @@ -157,7 +157,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.C index 8e979312267..88c30c6e0e2 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.C @@ -241,7 +241,7 @@ void Foam::ConeInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { const label i = parcelI % positionAxis_.size(); @@ -249,7 +249,7 @@ void Foam::ConeInjection<CloudType>::setPositionAndCell position = positionAxis_[i].first(); cellOwner = injectorCells_[i]; tetFacei = injectorTetFaces_[i]; - tetPtI = injectorTetPts_[i]; + tetPti = injectorTetPts_[i]; } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.H index a482be6af8a..9cd2a829f55 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeInjection/ConeInjection.H @@ -172,7 +172,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C index 392f901f463..30e9540c9fd 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.C @@ -49,7 +49,7 @@ void Foam::ConeNozzleInjection<CloudType>::setInjectionMethod() ( injectorCell_, tetFacei_, - tetPtI_, + tetPti_, position_, false ); @@ -111,7 +111,7 @@ Foam::ConeNozzleInjection<CloudType>::ConeNozzleInjection position_(this->coeffDict().lookup("position")), injectorCell_(-1), tetFacei_(-1), - tetPtI_(-1), + tetPti_(-1), direction_(this->coeffDict().lookup("direction")), parcelsPerSecond_ ( @@ -215,7 +215,7 @@ Foam::ConeNozzleInjection<CloudType>::ConeNozzleInjection position_(im.position_), injectorCell_(im.injectorCell_), tetFacei_(im.tetFacei_), - tetPtI_(im.tetPtI_), + tetPti_(im.tetPti_), direction_(im.direction_), parcelsPerSecond_(im.parcelsPerSecond_), flowRateProfile_(im.flowRateProfile_), @@ -252,7 +252,7 @@ void Foam::ConeNozzleInjection<CloudType>::updateMesh() ( injectorCell_, tetFacei_, - tetPtI_, + tetPti_, position_ ); } @@ -316,7 +316,7 @@ void Foam::ConeNozzleInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { cachedRandom& rndGen = this->owner().rndGen(); @@ -331,7 +331,7 @@ void Foam::ConeNozzleInjection<CloudType>::setPositionAndCell position = position_; cellOwner = injectorCell_; tetFacei = tetFacei_; - tetPtI = tetPtI_; + tetPti = tetPti_; break; } @@ -346,7 +346,7 @@ void Foam::ConeNozzleInjection<CloudType>::setPositionAndCell ( cellOwner, tetFacei, - tetPtI, + tetPti, position, false ); diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.H index 5d7e34f1180..1ccce9208e8 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ConeNozzleInjection/ConeNozzleInjection.H @@ -131,7 +131,7 @@ private: label tetFacei_; //- Index of tet point for injector cell - label tetPtI_; + label tetPti_; //- Injector direction [] vector direction_; @@ -244,7 +244,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.C index eb93b942d99..94c25b01c8e 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.C @@ -201,13 +201,13 @@ void Foam::FieldActivatedInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { position = positions_[parcelI]; cellOwner = injectorCells_[parcelI]; tetFacei = injectorTetFaces_[parcelI]; - tetPtI = injectorTetPts_[parcelI]; + tetPti = injectorTetPts_[parcelI]; } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.H index c55ebc829b6..85974580629 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/FieldActivatedInjection/FieldActivatedInjection.H @@ -174,7 +174,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.C index 8951c6c6536..ca0ac40eba3 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.C @@ -446,7 +446,7 @@ void Foam::InflationInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { position = newParticles_[parcelI].first().first(); @@ -455,7 +455,7 @@ void Foam::InflationInjection<CloudType>::setPositionAndCell ( cellOwner, tetFacei, - tetPtI, + tetPti, position, false ); diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.H index 7f0bc21c2e6..e2248282d58 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InflationInjection/InflationInjection.H @@ -170,7 +170,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C index 77c18d476e4..3042707ca67 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C @@ -73,7 +73,7 @@ bool Foam::InjectionModel<CloudType>::prepareForNextTimeStep } else { - // injection should have started, but not sufficient volume to + // Injection should have started, but not sufficient volume to // produce (at least) 1 parcel - hold value of timeStep0_ validInjection = false; } @@ -93,7 +93,7 @@ bool Foam::InjectionModel<CloudType>::findCellAtPosition ( label& celli, label& tetFacei, - label& tetPtI, + label& tetPti, vector& position, bool errorOnNotFound ) @@ -107,7 +107,7 @@ bool Foam::InjectionModel<CloudType>::findCellAtPosition position, celli, tetFacei, - tetPtI + tetPti ); label proci = -1; @@ -125,7 +125,7 @@ bool Foam::InjectionModel<CloudType>::findCellAtPosition { celli = -1; tetFacei = -1; - tetPtI = -1; + tetPti = -1; } // Last chance - find nearest cell and try that one - the point is @@ -138,7 +138,15 @@ bool Foam::InjectionModel<CloudType>::findCellAtPosition { position += SMALL*(cellCentres[celli] - position); - if (this->owner().mesh().pointInCell(position, celli)) + this->owner().mesh().findCellFacePt + ( + position, + celli, + tetFacei, + tetPti + ); + + if (celli > 0) { proci = Pstream::myProcNo(); } @@ -150,7 +158,7 @@ bool Foam::InjectionModel<CloudType>::findCellAtPosition { celli = -1; tetFacei = -1; - tetPtI = -1; + tetPti = -1; } } @@ -377,9 +385,7 @@ Foam::InjectionModel<CloudType>::~InjectionModel() template<class CloudType> void Foam::InjectionModel<CloudType>::updateMesh() -{ - // do nothing -} +{} template<class CloudType> @@ -443,7 +449,7 @@ void Foam::InjectionModel<CloudType>::inject(TrackData& td) // tetFace and tetPt label celli = -1; label tetFacei = -1; - label tetPtI = -1; + label tetPti = -1; vector pos = Zero; @@ -455,7 +461,7 @@ void Foam::InjectionModel<CloudType>::inject(TrackData& td) pos, celli, tetFacei, - tetPtI + tetPti ); if (celli > -1) @@ -468,7 +474,7 @@ void Foam::InjectionModel<CloudType>::inject(TrackData& td) // Create a new parcel parcelType* pPtr = - new parcelType(mesh, pos, celli, tetFacei, tetPtI); + new parcelType(mesh, pos, celli, tetFacei, tetPti); // Check/set new parcel thermo properties cloud.setParcelThermoProperties(*pPtr, dt); @@ -563,7 +569,7 @@ void Foam::InjectionModel<CloudType>::injectSteadyState // tetFace and tetPt label celli = -1; label tetFacei = -1; - label tetPtI = -1; + label tetPti = -1; vector pos = Zero; @@ -575,7 +581,7 @@ void Foam::InjectionModel<CloudType>::injectSteadyState pos, celli, tetFacei, - tetPtI + tetPti ); if (celli > -1) @@ -585,7 +591,7 @@ void Foam::InjectionModel<CloudType>::injectSteadyState // Create a new parcel parcelType* pPtr = - new parcelType(mesh, pos, celli, tetFacei, tetPtI); + new parcelType(mesh, pos, celli, tetFacei, tetPti); // Check/set new parcel thermo properties cloud.setParcelThermoProperties(*pPtr, 0.0); diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.H index 9e219bca101..e1901c69d7b 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.H @@ -160,7 +160,7 @@ protected: ( label& celli, label& tetFacei, - label& tetPtI, + label& tetPti, vector& position, bool errorOnNotFound = true ); @@ -320,7 +320,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) = 0; //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.C index 3708cb89dc3..d081be437ff 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.C @@ -176,7 +176,7 @@ void Foam::KinematicLookupTableInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { label injectorI = 0; @@ -193,7 +193,7 @@ void Foam::KinematicLookupTableInjection<CloudType>::setPositionAndCell position = injectors_[injectorI].x(); cellOwner = injectorCells_[injectorI]; tetFacei = injectorTetFaces_[injectorI]; - tetPtI = injectorTetPts_[injectorI]; + tetPti = injectorTetPts_[injectorI]; } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.H index 622ac451c03..565a6b221a6 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.H @@ -160,7 +160,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.C index d19010de6f4..21fc2471f66 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.C @@ -206,13 +206,13 @@ void Foam::ManualInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { position = positions_[parcelI]; cellOwner = injectorCells_[parcelI]; tetFacei = injectorTetFaces_[parcelI]; - tetPtI = injectorTetPts_[parcelI]; + tetPti = injectorTetPts_[parcelI]; } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.H index 01ecf0d5c6a..d8789a96463 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/ManualInjection/ManualInjection.H @@ -151,7 +151,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/NoInjection/NoInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/NoInjection/NoInjection.H index ae8a082ba9c..9aaf8793875 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/NoInjection/NoInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/NoInjection/NoInjection.H @@ -106,7 +106,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); virtual void setProperties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C index 4703885b37a..01873a3091d 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C @@ -221,7 +221,7 @@ void Foam::PatchFlowRateInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { patchInjectionBase::setPositionAndCell @@ -231,7 +231,7 @@ void Foam::PatchFlowRateInjection<CloudType>::setPositionAndCell position, cellOwner, tetFacei, - tetPtI + tetPti ); } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.H index df767d2fd66..9df1aac0b99 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.H @@ -157,7 +157,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); virtual void setProperties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.C index ef10b640cf9..b73c4c0d749 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.C @@ -176,7 +176,7 @@ void Foam::PatchInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { patchInjectionBase::setPositionAndCell @@ -186,7 +186,7 @@ void Foam::PatchInjection<CloudType>::setPositionAndCell position, cellOwner, tetFacei, - tetPtI + tetPti ); } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.H index b1c0825ecf8..61ce2c971db 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.H @@ -152,7 +152,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); virtual void setProperties diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.C index a83812ab010..082c2a73273 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.C @@ -153,7 +153,7 @@ void Foam::patchInjectionBase::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { scalar areaFraction = rnd.globalPosition(scalar(0), patchArea_); @@ -209,13 +209,13 @@ void Foam::patchInjectionBase::setPositionAndCell // the base point on the face as the tetPt. The tracking will pick // the cell consistent with the motion in the first tracking step tetFacei = mesh.cells()[cellOwner][0]; - tetPtI = 1; + tetPti = 1; } else { cellOwner = -1; tetFacei = -1; - tetPtI = -1; + tetPti = -1; // Dummy position position = pTraits<vector>::max; @@ -225,7 +225,7 @@ void Foam::patchInjectionBase::setPositionAndCell { cellOwner = -1; tetFacei = -1; - tetPtI = -1; + tetPti = -1; // Dummy position position = pTraits<vector>::max; diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.H b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.H index e6fdd27c9c8..47c003b60b7 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.H +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/patchInjectionBase.H @@ -121,7 +121,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); }; diff --git a/src/lagrangian/intermediate/submodels/Kinematic/SurfaceFilmModel/SurfaceFilmModel/SurfaceFilmModel.C b/src/lagrangian/intermediate/submodels/Kinematic/SurfaceFilmModel/SurfaceFilmModel/SurfaceFilmModel.C index d7be65546ce..c5016097725 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/SurfaceFilmModel/SurfaceFilmModel/SurfaceFilmModel.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/SurfaceFilmModel/SurfaceFilmModel/SurfaceFilmModel.C @@ -152,7 +152,7 @@ void Foam::SurfaceFilmModel<CloudType>::inject(TrackData& td) // point as the tetPt. The tracking will pick the cell // consistent with the motion in the first tracking step. const label tetFacei = this->owner().mesh().cells()[celli][0]; - const label tetPtI = 1; + const label tetPti = 1; // const point& pos = this->owner().mesh().C()[celli]; @@ -172,7 +172,7 @@ void Foam::SurfaceFilmModel<CloudType>::inject(TrackData& td) pos, celli, tetFacei, - tetPtI + tetPti ); // Check/set new parcel thermo properties diff --git a/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.C b/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.C index f0b33169cdb..40290f9a851 100644 --- a/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.C +++ b/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.C @@ -175,7 +175,7 @@ void Foam::ReactingLookupTableInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { label injectorI = 0; @@ -192,7 +192,7 @@ void Foam::ReactingLookupTableInjection<CloudType>::setPositionAndCell position = injectors_[injectorI].x(); cellOwner = injectorCells_[injectorI]; tetFacei = injectorTetFaces_[injectorI]; - tetPtI = injectorTetPts_[injectorI]; + tetPti = injectorTetPts_[injectorI]; } diff --git a/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.H b/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.H index 4d18307b48e..100046e5f82 100644 --- a/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.H +++ b/src/lagrangian/intermediate/submodels/Reacting/InjectionModel/ReactingLookupTableInjection/ReactingLookupTableInjection.H @@ -160,7 +160,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.C b/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.C index 9db4d1cd2d0..7361402d21b 100644 --- a/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.C +++ b/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.C @@ -181,7 +181,7 @@ void Foam::ReactingMultiphaseLookupTableInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { label injectorI = 0; @@ -198,7 +198,7 @@ void Foam::ReactingMultiphaseLookupTableInjection<CloudType>::setPositionAndCell position = injectors_[injectorI].x(); cellOwner = injectorCells_[injectorI]; tetFacei = injectorTetFaces_[injectorI]; - tetPtI = injectorTetPts_[injectorI]; + tetPti = injectorTetPts_[injectorI]; } diff --git a/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.H b/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.H index 9a01e6314cb..36276e54872 100644 --- a/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.H +++ b/src/lagrangian/intermediate/submodels/ReactingMultiphase/InjectionModel/ReactingMultiphaseLookupTableInjection/ReactingMultiphaseLookupTableInjection.H @@ -162,7 +162,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.C b/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.C index 35e14f2682e..b93166d4d4a 100644 --- a/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.C +++ b/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.C @@ -176,7 +176,7 @@ void Foam::ThermoLookupTableInjection<CloudType>::setPositionAndCell vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ) { label injectorI = 0; @@ -193,7 +193,7 @@ void Foam::ThermoLookupTableInjection<CloudType>::setPositionAndCell position = injectors_[injectorI].x(); cellOwner = injectorCells_[injectorI]; tetFacei = injectorTetFaces_[injectorI]; - tetPtI = injectorTetPts_[injectorI]; + tetPti = injectorTetPts_[injectorI]; } diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.H b/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.H index 4fb4b10dcde..bb68a49ff4f 100644 --- a/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.H +++ b/src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.H @@ -160,7 +160,7 @@ public: vector& position, label& cellOwner, label& tetFacei, - label& tetPtI + label& tetPti ); //- Set the parcel properties diff --git a/src/lagrangian/molecularDynamics/molecule/molecule/molecule.H b/src/lagrangian/molecularDynamics/molecule/molecule/molecule.H index a54c8744952..d2e38e5c296 100644 --- a/src/lagrangian/molecularDynamics/molecule/molecule/molecule.H +++ b/src/lagrangian/molecularDynamics/molecule/molecule/molecule.H @@ -242,7 +242,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const tensor& Q, const vector& v, const vector& a, diff --git a/src/lagrangian/molecularDynamics/molecule/molecule/moleculeI.H b/src/lagrangian/molecularDynamics/molecule/molecule/moleculeI.H index c9c4077949a..5c3058d0d4e 100644 --- a/src/lagrangian/molecularDynamics/molecule/molecule/moleculeI.H +++ b/src/lagrangian/molecularDynamics/molecule/molecule/moleculeI.H @@ -223,7 +223,7 @@ inline Foam::molecule::molecule const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const tensor& Q, const vector& v, const vector& a, @@ -236,7 +236,7 @@ inline Foam::molecule::molecule ) : - particle(mesh, position, celli, tetFacei, tetPtI), + particle(mesh, position, celli, tetFacei, tetPti), Q_(Q), v_(v), a_(a), diff --git a/src/lagrangian/solidParticle/solidParticle.H b/src/lagrangian/solidParticle/solidParticle.H index 72af01cf281..1a0258a48fd 100644 --- a/src/lagrangian/solidParticle/solidParticle.H +++ b/src/lagrangian/solidParticle/solidParticle.H @@ -133,7 +133,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const scalar d, const vector& U ); diff --git a/src/lagrangian/solidParticle/solidParticleI.H b/src/lagrangian/solidParticle/solidParticleI.H index f1074b7298b..c1b6fdd4189 100644 --- a/src/lagrangian/solidParticle/solidParticleI.H +++ b/src/lagrangian/solidParticle/solidParticleI.H @@ -48,12 +48,12 @@ inline Foam::solidParticle::solidParticle const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const scalar d, const vector& U ) : - particle(mesh, position, celli, tetFacei, tetPtI), + particle(mesh, position, celli, tetFacei, tetPti), d_(d), U_(U) {} diff --git a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.H b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.H index 056c2b66690..4f3c6246901 100644 --- a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.H +++ b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.H @@ -192,7 +192,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ); //- Construct from components @@ -202,7 +202,7 @@ public: const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, diff --git a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcelI.H b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcelI.H index 526341eaf2e..d0f7391b826 100644 --- a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcelI.H +++ b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcelI.H @@ -110,10 +110,10 @@ inline Foam::SprayParcel<ParcelType>::SprayParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI + const label tetPti ) : - ParcelType(mesh, position, celli, tetFacei, tetPtI), + ParcelType(mesh, position, celli, tetFacei, tetPti), d0_(this->d()), position0_(position), sigma_(0.0), @@ -137,7 +137,7 @@ inline Foam::SprayParcel<ParcelType>::SprayParcel const vector& position, const label celli, const label tetFacei, - const label tetPtI, + const label tetPti, const label typeId, const scalar nParticle0, const scalar d0, @@ -165,7 +165,7 @@ inline Foam::SprayParcel<ParcelType>::SprayParcel position, celli, tetFacei, - tetPtI, + tetPti, typeId, nParticle0, d0, -- GitLab