Commit c0952bdb authored by Jameson Graef Rollins's avatar Jameson Graef Rollins

use gwinc.finesse in substrate thermal calculations

Arm cavity finesse is already calculated in precompIFO, so just use that
instead of re-calculating.
parent 106863c7
Pipeline #77276 passed with stages
in 59 seconds
...@@ -21,10 +21,8 @@ def carrierdensity_adiabatic(f, ifo): ...@@ -21,10 +21,8 @@ def carrierdensity_adiabatic(f, ifo):
diffHole = ifo.Materials.Substrate.HoleDiffusion diffHole = ifo.Materials.Substrate.HoleDiffusion
cdDens = ifo.Materials.Substrate.CarrierDensity cdDens = ifo.Materials.Substrate.CarrierDensity
r0 = ifo.Optics.ITM.BeamRadius/np.sqrt(2) r0 = ifo.Optics.ITM.BeamRadius/np.sqrt(2)
T = ifo.Optics.ITM.Transmittance
L = ifo.Infrastructure.Length L = ifo.Infrastructure.Length
Finesse = 2*pi/T gPhase = ifo.gwinc.finesse*2/pi
gPhase = 2*Finesse/pi
psdElec = 4*H*gammaElec**2*cdDens*diffElec/(pi*r0**4*Omega**2) # units are meters psdElec = 4*H*gammaElec**2*cdDens*diffElec/(pi*r0**4*Omega**2) # units are meters
psdHole = 4*H*gammaHole**2*cdDens*diffHole/(pi*r0**4*Omega**2) # units are meters psdHole = 4*H*gammaHole**2*cdDens*diffHole/(pi*r0**4*Omega**2) # units are meters
...@@ -53,12 +51,7 @@ def carrierdensity_exact(f, ifo): ...@@ -53,12 +51,7 @@ def carrierdensity_exact(f, ifo):
cdDens = ifo.Materials.Substrate.CarrierDensity cdDens = ifo.Materials.Substrate.CarrierDensity
gammaElec = ifo.Materials.Substrate.ElectronIndexGamma gammaElec = ifo.Materials.Substrate.ElectronIndexGamma
gammaHole = ifo.Materials.Substrate.HoleIndexGamma gammaHole = ifo.Materials.Substrate.HoleIndexGamma
gPhase = ifo.gwinc.finesse*2/pi
T = ifo.Optics.ITM.Transmittance
FSR = c/(2*L) # in Hz
Finesse = 2*pi/T
cavPole = FSR/(2*Finesse) # in Hz
gPhase = 2*Finesse/pi
omega = 2*pi*f omega = 2*pi*f
...@@ -99,9 +92,7 @@ def thermorefractiveITM_adiabatic(f, ifo): ...@@ -99,9 +92,7 @@ def thermorefractiveITM_adiabatic(f, ifo):
kBT = scipy.constants.k * Temp kBT = scipy.constants.k * Temp
r0 = ifo.Optics.ITM.BeamRadius/np.sqrt(2) r0 = ifo.Optics.ITM.BeamRadius/np.sqrt(2)
L = ifo.Infrastructure.Length L = ifo.Infrastructure.Length
T = ifo.Optics.ITM.Transmittance gPhase = ifo.gwinc.finesse*2/pi
Finesse = 2*pi/T
gPhase = 2*Finesse/pi
psd = 4*H*beta**2*kappa*kBT*Temp/(pi*r0**4*Omega**2*(rho*C)**2) # units are meters psd = 4*H*beta**2*kappa*kBT*Temp/(pi*r0**4*Omega**2*(rho*C)**2) # units are meters
psdMeters = 2*psd # two ITMs psdMeters = 2*psd # two ITMs
...@@ -127,13 +118,8 @@ def thermorefractiveITM_exact(f, ifo): ...@@ -127,13 +118,8 @@ def thermorefractiveITM_exact(f, ifo):
beta = ifo.Materials.Substrate.dndT beta = ifo.Materials.Substrate.dndT
C = ifo.Materials.Substrate.MassCM C = ifo.Materials.Substrate.MassCM
kappa = ifo.Materials.Substrate.MassKappa kappa = ifo.Materials.Substrate.MassKappa
T = ifo.Optics.ITM.Transmittance gPhase = ifo.gwinc.finesse*2/pi
FSR = c/(2*L) # in Hz
Finesse = 2*pi/T
#cavPole = FSR/(2*Finesse) # in Hz
#gPhase = 2*Finesse/pi * (1 + (f/cavPole)**2)**(-1/2)
gPhase = 2*Finesse/pi
omega = 2*pi*f omega = 2*pi*f
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment