From 681efab1495d4959d01aec01db1397c10457cacf Mon Sep 17 00:00:00 2001 From: Ilan Upfal <103550915+iupfal@users.noreply.github.com> Date: Mon, 9 Mar 2026 07:09:41 -0400 Subject: [PATCH 1/2] fix tangential induction --- MITRotor/TangentialInduction.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/MITRotor/TangentialInduction.py b/MITRotor/TangentialInduction.py index cd493ec..d44369c 100644 --- a/MITRotor/TangentialInduction.py +++ b/MITRotor/TangentialInduction.py @@ -57,7 +57,7 @@ def __call__( eff_yaw = calc_eff_yaw(yaw, tilt) aprime = ( np.clip(aero_props.C_tau_corr, -2, 2) - / (4 * np.maximum(geom.mu_mesh, 0.1) ** 2 * tsr * (1 - aero_props.an) * np.cos(eff_yaw)) + / (4 * np.maximum(geom.mu_mesh, 0.1) * tsr * (1 - aero_props.an) * np.cos(eff_yaw)) ) return aprime From 2a78178948a82b8ae826032634f27bdd70d9a46f Mon Sep 17 00:00:00 2001 From: Ilan Upfal <103550915+iupfal@users.noreply.github.com> Date: Mon, 9 Mar 2026 07:13:19 -0400 Subject: [PATCH 2/2] increase iterations --- MITRotor/BEMSolver.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/MITRotor/BEMSolver.py b/MITRotor/BEMSolver.py index daea9bc..559e657 100644 --- a/MITRotor/BEMSolver.py +++ b/MITRotor/BEMSolver.py @@ -130,7 +130,7 @@ def Ctprime(self, grid: Literal["sector", "annulus", "rotor"] = "rotor"): return average(self.geom, Ctprime, grid=grid) -@adaptivefixedpointiteration(max_iter=500, relaxations=[0.25, 0.5, 0.96]) +@adaptivefixedpointiteration(max_iter=1000, relaxations=[0.25, 0.5, 0.96]) class BEM: """ A generic BEM class which facilitates dependency injection for various models.