Skip to content
Snippets Groups Projects
Commit b69d4267 authored by Christopher Wipf's avatar Christopher Wipf
Browse files

documentation fixes

parent 85c771da
No related branches found
No related tags found
No related merge requests found
Pipeline #121721 passed
...@@ -61,8 +61,6 @@ def substrate_carrierdensity(f, materials, wBeam, exact=False): ...@@ -61,8 +61,6 @@ def substrate_carrierdensity(f, materials, wBeam, exact=False):
def substrate_thermorefractive(f, materials, wBeam, exact=False): def substrate_thermorefractive(f, materials, wBeam, exact=False):
"""Substrate thermal displacement noise spectrum from thermorefractive fluctuations """Substrate thermal displacement noise spectrum from thermorefractive fluctuations
For semiconductor substrates.
:f: frequency array in Hz :f: frequency array in Hz
:materials: gwinc optic materials structure :materials: gwinc optic materials structure
:wBeam: beam radius (at 1 / e^2 power) :wBeam: beam radius (at 1 / e^2 power)
......
...@@ -143,11 +143,13 @@ def wireGeometry(r, N, RibbonThickness=None, TaperedEndRadius=None, **kwargs): ...@@ -143,11 +143,13 @@ def wireGeometry(r, N, RibbonThickness=None, TaperedEndRadius=None, **kwargs):
"""Compute wire geometry-dependent factors """Compute wire geometry-dependent factors
r is the wire radius, or ribbon width. r is the wire radius, or ribbon width.
N is the number of wires or ribbons
RibbonThickness must be set when ribbons are used, or RibbonThickness must be set when ribbons are used, or
TaperedEndRadius when tapered fibers are used. TaperedEndRadius when tapered fibers are used.
Other kwargs are ignored.
Returns the cross-sectional area and moment of inertia, Returns cross-sectional areas (central and end) and moment of inertia,
and the modified surface to volume ratios (vertical and horizontal). and modified surface to volume ratios (vertical and horizontal).
""" """
# Usual case: round wire/fiber # Usual case: round wire/fiber
...@@ -203,12 +205,22 @@ def wireTELoss(w, tension, xsectEnd, xII, Temp, alpha, beta, rho, C, K, Y, ...@@ -203,12 +205,22 @@ def wireTELoss(w, tension, xsectEnd, xII, Temp, alpha, beta, rho, C, K, Y,
"""Thermoelastic calculations for wires """Thermoelastic calculations for wires
Repeated for upper and lower joint of each stage. Repeated for upper and lower joint of each stage.
w = angular frequency
tension = weight supported per wire
xsectEnd = cross sectional area of wire end
xII = cross sectional moment of inertia
Temp = temperature
alpha = coeff of thermal expansion alpha = coeff of thermal expansion
beta = temp dependence of Young's modulus beta = temp dependence of Young's modulus
rho = mass density rho = mass density
C = heat capacity C = heat capacity
K = thermal conductivity W/(m K) K = thermal conductivity W/(m K)
Y = Young's modulus Y = Young's modulus
RibbonThickness must be set when ribbons are used
Other kwargs are ignored
Returns the loss angle associated with thermoelastic damping
(wire horizontal)
""" """
# horizontal TE time constant, wires # horizontal TE time constant, wires
...@@ -237,6 +249,9 @@ def bladeTELoss(w, t, Temp, alpha, beta, rho, C, K, Y): ...@@ -237,6 +249,9 @@ def bladeTELoss(w, t, Temp, alpha, beta, rho, C, K, Y):
"""Thermoelastic calculations for blades """Thermoelastic calculations for blades
Invoked for upper joint only (there is no lower blade) Invoked for upper joint only (there is no lower blade)
w = angular frequency
t = blade thickness
Temp = temperature
alpha = coeff of thermal expansion alpha = coeff of thermal expansion
beta = temp dependence of Young's modulus beta = temp dependence of Young's modulus
rho = mass density rho = mass density
...@@ -244,6 +259,9 @@ def bladeTELoss(w, t, Temp, alpha, beta, rho, C, K, Y): ...@@ -244,6 +259,9 @@ def bladeTELoss(w, t, Temp, alpha, beta, rho, C, K, Y):
K = thermal conductivity W/(m K) K = thermal conductivity W/(m K)
Y = Young's modulus Y = Young's modulus
Returns the loss angle associated with thermoelastic damping
(blade vertical)
""" """
# vertical TE time constant, blades # vertical TE time constant, blades
tau = (rho * C * t**2) / (K * pi**2) tau = (rho * C * t**2) / (K * pi**2)
...@@ -274,6 +292,17 @@ def bladeTELoss(w, t, Temp, alpha, beta, rho, C, K, Y): ...@@ -274,6 +292,17 @@ def bladeTELoss(w, t, Temp, alpha, beta, rho, C, K, Y):
def continuumWireKh(w, N, length, tension, xsect, xII, rho, Y, phi): def continuumWireKh(w, N, length, tension, xsect, xII, rho, Y, phi):
"""Horizontal spring constant, including violin modes """Horizontal spring constant, including violin modes
w = angular frequency
N = number of wires
length = wire length
tension = weight supported per wire
xsect = wire cross sectional area
xII = cross sectional moment of inertia
rho = mass density
Y = Young's modulus
phi = loss angle
Returns the spring constant (wire horizontal)
""" """
Y = Y * (1 + 1j * phi) # complex Young's modulus Y = Y * (1 + 1j * phi) # complex Young's modulus
...@@ -309,6 +338,18 @@ def continuumWireKv(w, N, length, xsect, xsectEnd, rho, Y, phi, ...@@ -309,6 +338,18 @@ def continuumWireKv(w, N, length, xsect, xsectEnd, rho, Y, phi,
TaperedEndLength=None, **kwargs): TaperedEndLength=None, **kwargs):
"""Vertical spring constant, including bounce mode. """Vertical spring constant, including bounce mode.
w = angular frequency
N = number of wires
length = wire length
xsect = wire cross sectional area
xsectEnd = cross sectional area of wire end
rho = mass density
Y = Young's modulus
phi = loss angle
TaperedEndLength must be set when tapered fibers are used
Other kwargs are ignored
Returns the spring constant (wire vertical)
""" """
Y = Y * (1 + 1j * phi) # complex Young's modulus Y = Y * (1 + 1j * phi) # complex Young's modulus
k = sqrt(rho / Y) * w k = sqrt(rho / Y) * w
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment