From 6940748b095a01af3f0b2540701bd1a766d9bdd4 Mon Sep 17 00:00:00 2001
From: "Mads M. Pedersen" <mmpe@dtu.dk>
Date: Wed, 19 Mar 2025 07:47:05 +0000
Subject: [PATCH] Set ctlim = 0.96 in Nygaard_2022. which corresponds to
 limiting the term in...

---
 py_wake/literature/turbopark.py | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/py_wake/literature/turbopark.py b/py_wake/literature/turbopark.py
index f5b416de1..ca4a6564d 100644
--- a/py_wake/literature/turbopark.py
+++ b/py_wake/literature/turbopark.py
@@ -27,7 +27,11 @@ class Nygaard_2022(PropagateDownwind):
         wake_deficitModel = TurboGaussianDeficit(
             ct2a=ct2a_mom1d,
             groundModel=Mirror(),
-            rotorAvgModel=GaussianOverlapAvgModel())
+            rotorAvgModel=GaussianOverlapAvgModel(),
+            # ctlim=0.96 corresponds to limiting the term in the sqrt of epsilon to 3 as in
+            # the matlab code from Ørsted (https://github.com/OrstedRD/TurbOPark/blob/main/TurbOPark.m#L197)
+            ctlim=0.96
+        )
 
         # Ørsted scales the deficit with respect to the ambient wind speed of the downstream turbine:
         wake_deficitModel.WS_key = 'WS_jlk'
-- 
GitLab