From 3f322e8b909886de5188023ee6ca81cae776d465 Mon Sep 17 00:00:00 2001 From: Henry Date: Mon, 9 Mar 2015 10:40:51 +0000 Subject: [PATCH] functionObjects/forces: limit the bins to handle moving meshes This approach simply accumulates data outside the range of the bins into the first or last bin which may be OK for small motions. For larger motions it may be better if the bins are updated or operate in a coordinate system attached to the body for solid-body motion. Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1560 --- .../functionObjects/forces/forces/forces.C | 15 ++++++++------- 1 file changed, 8 insertions(+), 7 deletions(-) diff --git a/src/postProcessing/functionObjects/forces/forces/forces.C b/src/postProcessing/functionObjects/forces/forces/forces.C index fa8cfdcfc..1794a8b88 100644 --- a/src/postProcessing/functionObjects/forces/forces/forces.C +++ b/src/postProcessing/functionObjects/forces/forces/forces.C @@ -375,13 +375,14 @@ void Foam::forces::applyBins forAll(dd, i) { - label binI = floor(dd[i]/binDx_); - force_[0][binI] += fN[i]; - force_[1][binI] += fT[i]; - force_[2][binI] += fP[i]; - moment_[0][binI] += Md[i]^fN[i]; - moment_[1][binI] += Md[i]^fT[i]; - moment_[2][binI] += Md[i]^fP[i]; + label bini = min(max(floor(dd[i]/binDx_), 0), force_[0].size()); + + force_[0][bini] += fN[i]; + force_[1][bini] += fT[i]; + force_[2][bini] += fP[i]; + moment_[0][bini] += Md[i]^fN[i]; + moment_[1][bini] += Md[i]^fT[i]; + moment_[2][bini] += Md[i]^fP[i]; } } }