From 7a1c873262b7301dae16dd4b7cef1e6451a57613 Mon Sep 17 00:00:00 2001
From: Christopher Wipf <wipf@ligo.mit.edu>
Date: Mon, 27 Aug 2018 16:54:04 -0700
Subject: [PATCH] de-mystify a constant in the SUS TE calculation
---
gwinc/suspension.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/gwinc/suspension.py b/gwinc/suspension.py
index 424f7483..101bf63e 100644
--- a/gwinc/suspension.py
+++ b/gwinc/suspension.py
@@ -234,7 +234,8 @@ def suspQuad(f, ifo, material='Silica'):
mu_h = 4 / r_w # surface to volume ratio, horizontal
mu_v = 2 / r_w # surface to volume ratio, vertical (wire)
- # horizontal TE time constant, wires ( WHAT IS THIS CONSTANT 7.37e-2? )
+ # horizontal TE time constant, wires
+ # The constant 7.37e-2 is 1/(4*q0^2) from eq 12, C. Zener 10.1103/PhysRev.53.90 (1938)
tau_h = 7.37e-2 * 4 * (rho_w * C_w * xsect) / (pi * K_w)
# vertical TE time constant, blades
--
GitLab