From 767f4b0bebc88ae2660b30109ef339f3d903c1f1 Mon Sep 17 00:00:00 2001
From: Thomas Lambert <t.lambert@uliege.be>
Date: Fri, 19 May 2023 11:41:59 +0200
Subject: [PATCH] refact(indvel): more logical place for nBlades

---
 src/configs/caradonna1981.m | 2 +-
 src/solvers/indvel.m        | 6 +++---
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/configs/caradonna1981.m b/src/configs/caradonna1981.m
index 8c49f3b..cdc0bd4 100644
--- a/src/configs/caradonna1981.m
+++ b/src/configs/caradonna1981.m
@@ -94,7 +94,7 @@ Airfoil.polarType = 'file'; % Type of polar to use ('file', 'polynomial')
 Airfoil.coordFile = 'airfoil_data/naca0012.dat';
 
 % If Airfoil.polarType == 'file'
-Airfoil.polarFile = 'airfoil_data/NACA_0012-Re_2e5-1e7.mat';
+Airfoil.polarFile = 'airfoil_data/NACA_0012-Re_1e5-1e7.mat';
 Airfoil.extrapMethod = 'viterna'; % Polar extrapol. ('none', 'spline', 'Viterna')
 
 % ==================================================================================================
diff --git a/src/solvers/indvel.m b/src/solvers/indvel.m
index b7d0423..98e53ff 100644
--- a/src/solvers/indvel.m
+++ b/src/solvers/indvel.m
@@ -89,8 +89,8 @@ function indvel(OpRot, Mod)
         dD = coeff2force(cd, OpRot, relVel);
 
         % Thrust and torque contributions
-        dFa = dL .* cos(phi) - dD .* sin(phi);
-        dFu = dL .* sin(phi) + dD .* cos(phi);
+        dFa = OpRot.Rot.nBlades * (dL .* cos(phi) - dD .* sin(phi));
+        dFu = OpRot.Rot.nBlades * (dL .* sin(phi) + dD .* cos(phi));
 
         % Final system to solve for v and u
         vw = OpRot.Op.speed + dFa ./ dmdot ./ K_T;
@@ -125,7 +125,7 @@ function force = coeff2force(coeff, OpRot, speed)
     % COEFF2FORCE Convert aerodynamic coefficient to force
 
     rho = OpRot.Op.Flow.rho;
-    nondim = OpRot.Rot.nBlades .* (0.5 * rho * OpRot.Rot.Bl.area .* speed.^2);
+    nondim = 0.5 * rho * 2 * pi * OpRot.Rot.Bl.area .* speed.^2;
 
     force = coeff .* nondim;
 
-- 
GitLab