diff --git a/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.C b/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.C index 889e32abf1..03ab5d4285 100644 --- a/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.C +++ b/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.C @@ -251,6 +251,9 @@ void Foam::reconstruction::plicRDF::calcResidual const label celli = interfaceLabels_[i]; if (mag(normal_[celli]) == 0 || mag(interfaceNormal_[i]) == 0) { + normalResidual[i].celli = celli; + normalResidual[i].normalResidual = 0; + normalResidual[i].avgAngle = 0; continue; } diff --git a/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.H b/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.H index 6432b6632d..cbdb19d8cf 100644 --- a/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.H +++ b/src/transportModels/geometricVoF/reconstructionSchemes/plicSchemes/plicRDF/plicRDF.H @@ -81,9 +81,9 @@ class plicRDF //- residuals storage struct normalRes { - label celli; - scalar normalResidual; - scalar avgAngle; + label celli = {}; + scalar normalResidual = {}; + scalar avgAngle = {}; }; //- Reference to mesh