From 8d2301c46f67add811b3325c2ba676afa099c5e9 Mon Sep 17 00:00:00 2001
From: Christopher Wipf <wipf@ligo.mit.edu>
Date: Tue, 29 Aug 2017 18:59:49 -0700
Subject: [PATCH] ifo submodule; some comments on material properties
---
__init__.py | 2 --
IFOModel_123_2000.py => ifo/Voyager.py | 18 +++++++++---------
ifo/__init__.py | 0
IFOModel.py => ifo/aLIGO.py | 0
4 files changed, 9 insertions(+), 11 deletions(-)
rename IFOModel_123_2000.py => ifo/Voyager.py (97%)
create mode 100644 ifo/__init__.py
rename IFOModel.py => ifo/aLIGO.py (100%)
diff --git a/__init__.py b/__init__.py
index 48b19c23..a137975d 100644
--- a/__init__.py
+++ b/__init__.py
@@ -1,3 +1 @@
from gwinc import gwinc
-from IFOModel import IFOModel
-from IFOModel_123_2000 import IFOModel_123_2000
diff --git a/IFOModel_123_2000.py b/ifo/Voyager.py
similarity index 97%
rename from IFOModel_123_2000.py
rename to ifo/Voyager.py
index df8db04b..92b2a8c7 100644
--- a/IFOModel_123_2000.py
+++ b/ifo/Voyager.py
@@ -1,6 +1,6 @@
from __future__ import division, print_function
from numpy import pi, NaN
-from util import SpotSizes
+from ..util import SpotSizes
import scipy.constants
import scipy.special
from scipy.io import loadmat
@@ -8,7 +8,7 @@ from scipy.io.matlab.mio5_params import mat_struct
import os
-def IFOModel_123_2000():
+def IFOModel():
"""IFOMODEL returns a structure describing an IFO for use in
benchmark programs and noise simulator. Part of the gwinc
package, which provides science-grounded figures of merit for
@@ -227,7 +227,7 @@ def IFOModel_123_2000():
ifo.Materials.Coating.CVhighn = 345.6*2250 # volume-specific heat capacity (J/K/m^3); http://journals.aps.org/prl/pdf/10.1103/PhysRevLett.96.055902
ifo.Materials.Coating.Alphahighn = 1e-9 # zero crossing at 123 K
ifo.Materials.Coating.Betahighn = 1.4e-4 # dn/dT
- ifo.Materials.Coating.ThermalDiffusivityhighn = 1 # W/m/K
+ ifo.Materials.Coating.ThermalDiffusivityhighn = 1 # W/m/K (this is a misnomer, meant to be thermal conductivity not diffusivity)
ifo.Materials.Coating.Phihighn = 3e-5 # just a guess (depends on prep)
ifo.Materials.Coating.Indexhighn = 3.5
@@ -238,7 +238,7 @@ def IFOModel_123_2000():
ifo.Materials.Coating.CVlown = 1.6412e6 # volume-specific heat capacity (J/K/m^3); Crooks et al, Fejer et al
ifo.Materials.Coating.Alphalown = 5.1e-7 # Fejer et al
ifo.Materials.Coating.Betalown = 8e-6 # dn/dT, (ref. 14)
- ifo.Materials.Coating.ThermalDiffusivitylown = 1.38 # Fejer et al
+ ifo.Materials.Coating.ThermalDiffusivitylown = 1.38 # Fejer et al (this is a misnomer, meant to be thermal conductivity not diffusivity)
ifo.Materials.Coating.Philown = 1e-4 # ?
# calculated for 123 K and 2000 nm following
@@ -251,12 +251,12 @@ def IFOModel_123_2000():
# Silicon @ 120K (http://www.ioffe.ru/SVA/NSM/Semicond/Si/index.html)
ifo.Materials.Substrate = mat_struct()
- # phi_sub = c2 * f^(MechLossExp)
- ifo.Materials.Substrate.c2 = 3e-13 # Coeff of freq dep. term for bulk loss (Lam & Douglass, 1981)
- ifo.Materials.Substrate.MechanicalLossExponent = 1 # Exponent for freq dependence of silicon loss
+ # phi_sub = c2 * f^(MechLossExp)
+ ifo.Materials.Substrate.c2 = 3e-13 # Coeff of freq dep. term for bulk loss (Lam & Douglass, 1981)
+ ifo.Materials.Substrate.MechanicalLossExponent = 1 # Exponent for freq dependence of silicon loss
ifo.Materials.Substrate.Alphas = 5.2e-12 # Surface loss limit ???
- ifo.Materials.Substrate.MirrorY = 155.8e9 # N/m^2; Youngs modulus (ioffe)
- ifo.Materials.Substrate.MirrorSigma = 0.27 # kg/m^3; Poisson ratio (ioffe)
+ ifo.Materials.Substrate.MirrorY = 155.8e9 # N/m^2; Youngs modulus (ioffe) -- what about anisotropy??
+ ifo.Materials.Substrate.MirrorSigma = 0.27 # kg/m^3; Poisson ratio (ioffe) -- what about anisotropy??
ifo.Materials.Substrate.MassDensity = 2329 # kg/m^3; (ioffe)
ifo.Materials.Substrate.MassAlpha = 1e-9 # 1/K; CTE = 0 @ 120 K
ifo.Materials.Substrate.MassCM = 0.3*1000 # J/kg/K; specific heat (ioffe @ 120K)
diff --git a/ifo/__init__.py b/ifo/__init__.py
new file mode 100644
index 00000000..e69de29b
diff --git a/IFOModel.py b/ifo/aLIGO.py
similarity index 100%
rename from IFOModel.py
rename to ifo/aLIGO.py
--
GitLab