diff --git a/src/classes/@ElemPerf/updateupstreamvel.m b/src/classes/@ElemPerf/updateupstreamvel.m index cea27813be3c133eb4a8a63e1a287f98aec141c3..778d6a987b2d67be70bb69e0db2a507983f37028 100644 --- a/src/classes/@ElemPerf/updateupstreamvel.m +++ b/src/classes/@ElemPerf/updateupstreamvel.m @@ -46,7 +46,7 @@ function updateupstreamvel(self, PrevOpRot) % area ratio of previous rotor contraction vec_radiusRatio = sqrt(vec_areaRatio); - vcRadius = PrevOpRot.Rot.radius * vec_radiusRatio(1, end); % vena contracta radius + vcRadius = PrevOpRot.Rot.Bl.y(end) * vec_radiusRatio(1, end); % vena contracta radius vcRadIdx = find(self.Rot.Bl.y <= vcRadius, 1, 'last'); % tangential velocity at the previous rotor disk