diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C index 637bf15fa4..a6a3865ba1 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C @@ -39,6 +39,19 @@ void Foam::KinematicParcel::updateCellQuantities ) { rhoc_ = td.rhoInterp().interpolate(this->position(), cellI); + if (rhoc_ < SMALL) + { + WarningIn + ( + "void Foam::KinematicParcel::updateCellQuantities" + "(" + "TrackData&, " + "const scalar, " + "const label" + ")" + ) << "Density < " << SMALL << " in cell " << cellI << nl << endl; + } + Uc_ = td.UInterp().interpolate(this->position(), cellI); muc_ = td.muInterp().interpolate(this->position(), cellI); diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C index 4ff61b6602..8b59c5017a 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C @@ -40,6 +40,18 @@ void Foam::ReactingParcel::updateCellQuantities ThermoParcel::updateCellQuantities(td, dt, cellI); pc_ = td.pInterp().interpolate(this->position(), cellI); + if (pc_ < SMALL) + { + WarningIn + ( + "void Foam::ReactingParcel::updateCellQuantities" + "(" + "TrackData&, " + "const scalar, " + "const label" + ")" + ) << "Pressure < " << SMALL << " in cell " << cellI << nl << endl; + } } diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C index d8b999d0f5..b22a8e57e6 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C @@ -40,6 +40,19 @@ void Foam::ThermoParcel::updateCellQuantities KinematicParcel::updateCellQuantities(td, dt, cellI); Tc_ = td.TInterp().interpolate(this->position(), cellI); + if (Tc_ < SMALL) + { + WarningIn + ( + "void Foam::ThermoParcel::updateCellQuantities" + "(" + "TrackData&, " + "const scalar, " + "const label" + ")" + ) << "Temperature < " << SMALL << " in cell " << cellI << nl << endl; + } + cpc_ = td.cpInterp().interpolate(this->position(), cellI); }