From 521f581df99482e578d9e1c947bc10c9ac9a790b Mon Sep 17 00:00:00 2001
From: sp89hili <sp89hili@pmd196.csi.tu-darmstadt.de>
Date: Tue, 18 Feb 2025 17:53:33 +0100
Subject: [PATCH] Pu_offset outlet area by a according to Schneider
---
@BearImp/calcCap_puchtler2025.m | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/@BearImp/calcCap_puchtler2025.m b/@BearImp/calcCap_puchtler2025.m
index 33c2328..f9c7709 100644
--- a/@BearImp/calcCap_puchtler2025.m
+++ b/@BearImp/calcCap_puchtler2025.m
@@ -30,13 +30,13 @@ function C_out = calcCap_puchtler2025(s,R_RE,R_R,R_ZL,R_ZR,B,R_L,epsilon_r,alpha
Theta_1l = @(phi) atan((+B_RhL - Delta_x)./(Delta_z*cos(phi) - sign(R_L)*sqrt(R_ZL^2-Delta_z^2*sin(phi).^2)));
Theta_1r = @(phi) atan((-B_RhR - Delta_x)./(Delta_z*cos(phi) - sign(R_L)*sqrt(R_ZR^2-Delta_z^2*sin(phi).^2)));
if method.outlet == 'air'
- Theta_0lIn = @(phi) max(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) - max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),phi.*(phi>0)),Theta_1l(phi));
- Theta_0rIn = @(phi) min(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) + max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),phi.*(phi>0)),Theta_1r(phi));
+ Theta_0lIn = @(phi) max(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) - max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),(a/R_RE+phi).*(phi>0)),Theta_1l(phi));
+ Theta_0rIn = @(phi) min(atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) + max(a*sqrt(max(1/R_RE^2-phi.^2/b^2,0)),(a/R_RE+phi).*(phi>0)),Theta_1r(phi));
Theta_0lOut = @(phi) atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) - a*sqrt(max(1/R_RE^2-phi.^2/b^2,0));
Theta_0rOut = @(phi) atan((R_RE*sin(alpha))./(-Delta_z*cos(phi) + sign(R_L)*sqrt(R_L^2 -Delta_z^2*sin(phi).^2))) + a*sqrt(max(1/R_RE^2-phi.^2/b^2,0));
- phi_0 = a/R_RE/sqrt(1+a^2/b^2);
- phi_1l = fzero(@(phi) Theta_0l(phi) - phi - Theta_1l(phi),pi/4);
- phi_1r = fzero(@(phi) Theta_0r(phi) + phi - Theta_1r(phi),pi/4);
+ phi_0 = 0;%a/R_RE/sqrt(1+a^2/b^2);
+ phi_1l = fzero(@(phi) Theta_0l(phi) - (a/R_RE+phi) - Theta_1l(phi),pi/4);
+ phi_1r = fzero(@(phi) Theta_0r(phi) + (a/R_RE+phi) - Theta_1r(phi),pi/4);
% fprintf('phi_0 = %d | phi_1l = %d | phi_1r = %d \n',phi_0,phi_1l,phi_1r)
end
Theta_2l = @(phi) atan((+B/2 - Delta_x)./(Delta_z*cos(phi) - sign(R_L)*sqrt(R_ZL^2-Delta_z^2*sin(phi).^2)));
--
GitLab