Mathematical model of thermal shields for long-term stability optical resonators

Modern experiments aiming at tests of fundamental physics, like measuring gravitational waves or testing Lorentz Invariance with unprecedented accuracy, require thermal environments that are highly stable over long times. To achieve such a stability, the experiment including typically an optical resonator is nested in a thermal enclosure, which passively attenuates external temperature fluctuations to acceptable levels. These thermal shields are usually designed using tedious numerical simulations or with simple analytical models. In this paper, we propose an accurate analytical method to estimate the performance of passive thermal shields in the frequency domain, which allows for fast evaluation and optimization. The model analysis has also unveil interesting properties of the shields, such as dips in the transfer function for some frequencies under certain combinations of materials and geometries. We validate the results by comparing them to numerical simulations performed with commercial software based on finite element methods.


Introduction
The ever increasing accuracy of current and future fundamental physics experiments on ground and in space demands very stable and controlled environments to meliorate spurious effects such as temperature fluctuations.In this paper, we present a method that allows the rapid assessment of the performance of passive thermal shields in the frequency domain.This is done with analytical models instead of time consuming and precision limited numerical simulations.The frequency domain is of special interest since stability requirements are typically given in such domain.
One of the key applications of this method is the design of thermal isolations for optical resonators (ORs), which are used as frequency references in highly sensitive interferometers, in experiments testing fundamental physics, and as optical frequency standards, see e.g.[1][2][3][4][5], [6][7][8][9][10][11][12][13][14], and [15][16][17][18], respectively.ORs consist of a spacer and two high-reflectivity mirrors.The distance between the mirrors defines the resonance frequency, which is directly linked to the length of the spacer.Consequently, temperature fluctuations result in frequency fluctuations, i.e., δ ν ν where ν is the optical frequency, ℓ is the spacer length, α is the spacer coefficient of thermal expansion (CTE) and T is its temperature.Fluctuations of the variables are expressed as δ .The spacers of ORs are made of materials with very low CTEs, e.g.ultra low expansion glass (ULE) or Zerodur with CTEs about 10 −8 K −1 at room temperature [19,20].In addition, the ORs are often driven at temperatures close to their CTE zero crossing (or at the CTE minimum) [15,21].
In case these temperatures cannot be achieved, different techniques have been developed to tune the temperature of the zero crossing [22,23].Nevertheless, the accuracy of current optical experiments is already so high that these low CTEs do not suffice and temperature fluctuations need to be strongly attenuated.The damping of temperature fluctuations at frequencies in the milli-Hertz and sub-milli-Hertz regime by several orders of magnitude demands a particularly careful thermal shield design.This low frequency range is crucial in different missions and experiments since the expected science signals are in the same frequency band.This is especially important for the space-based gravitational wave detector eLISA [1] and foreseen missions testing the fundamental principles of general relativity such as the Lorentz Invariance with Michelson-Morley experiments and a Kennedy-Thorndike experiment like STAR (SpaceTime Asymmetry Research [11]), mini-STAR, and BOOST (BOOst Symmetry Test [12]), for earlier Kennedy-Thorndike experiments see, e.g., [6,7].
The thermal shield design is driven by the required temperature stability of the OR, i.e., the required frequency stability for a given spacer CTE, and the temperature variations at the outer most thermal shield layer.For eLISA the required temperature stability (if a cavity is used for the laser pre-stabilization) is µK Hz −1/2 in the mili-Hertz band.The thermal environment in eLISA will be extremely stable [24] and the required attenuation will not be very demanding.However, the main problem is for laboratory-based demonstration experiments [2,25] where µK Hz −1/2 temperature stability in the sub-milli-Hertz range is needed.For the Kennedy Thorndike experiments the allowed temperature fluctuations are typically 0.1 µK at the orbital period (∼90 minutes) and the expected stability in the satellites around 1 K and ∼10 mK can be achieved in the payload by active temperature control.This poses a stringent requirement on the thermal shields: about five orders of magnitude attenuation at ∼0.2 mHz.For the Gravity Recovery and Climate Experiment Follow-On (GRACE-FO) interferometer similar figures are required [5].The thermal shield performance is usually estimated with the help of tedious numerical simulations [16,24] or with rather simplified analytical models.However, a method to estimate its transfer function yielding results close to the numerical simulations is to our best knowledge still lacking.In this paper, we provide such a tool, using thermal shields for optical resonators as an example.
The paper is organized as follows: First the general mathematical model is introduced, which is applied to the cases of spherical and cylindrically shaped shields in Sec. 3.These results are compared with numerical simulations based on finite element methods (FEM) in Sec. 4. Sec. 5 summarizes our conclusions.Supplementary models and mathematical details are provided in the appendix.

Mathematical model of the thermal shields
The basic assumption of the thermal shields' mathematical model is that heat is transferred only by radiation and conduction.In Sec.2.1, the analytical transfer function of the thermal shields is derived, when the latter is negligible.In Sec.2.2, the conductive links between the shields are included in the model, too.

Radiative heat transfer
The radiative heat transfer between two gray bodies (from body i to body j) is defined as [26]: where σ is the Stefan-Boltzmann constant (= 5.67 × 10 −8 W m −2 K −4 ), A j is the area, T j is the temperature in Kelvin of the respective layers and β i j is a term including the view factors and emissivities, ε -see Sec. 3. The temperature change of the layer j due to the heat exchange given by Eq. ( 2) is where m j and c j are the mass and the specific heat of the layer j, respectively.We assume that the thermal shields have high thermal conductivity and are relatively thin.Consequently, the temperature distribution in the shields homogenizes rapidly even in the cases, where point heat sources are present -see Sec. 3. Hence, we do not consider a spatial temperature distribution.Combining Eqs. ( 2) and (3) results in the differential equation that has been linearized in order to find an analytical solution.T 0 is the average temperature of the bodies.Typically, temperature changes during the experiments are small compared to T 0 and, thus, the linearization is justified.The transfer function of the thermal shield (from layer i to j) is obtained after applying the Fourier transform to Eq. ( 4) and taking the ratio between both temperatures: which corresponds to a first-order low-pass filter with a cut-off angular frequency of Analogously to the electrical case, the cut-off angular frequency is split in a thermal resistance (in units of K W −1 ) and a thermal capacitance (in units of J K −1 ): C j = m j c j (8) and the product is the time constant of the system, τ i j = θ i j C j = ω −1 c .Once the transfer function of one shield (formed by 2 layers) has been determined the behavior of N concentric thermal shields (N + 1 layers) can be calculated.In order to obtain a compact expression, it is assumed that all the layers are exactly the same (in material and in size), i.e., they all have the same time constant τ i j = τ.The validity of this assumption is discussed in Sec. 3, where we show that it holds for thin layers that are sufficiently close to each other.In Appendix A, the transfer function for N layers is derived and the result is This is depicted in Fig. 1 for different N. The left panel shows, how the layers interact with each other causing some of the poles of the transfer function, cf.Eq. ( 9), of the system to appear at lower frequencies than the poles of each individual layer.For this reason, the system exhibits a stronger temperature damping at frequencies around the cut-off angular frequency than if the layers are assumed uncoupled, i.e., However, for ω 10ω c the response of both systems (coupled or uncoupled) is the same and Eq. ( 10) can be used instead of Eq. ( 9). 1. Left: transfer functions using the coupled case Eq. ( 9) (solid lines) and the uncoupled case Eq. (10) (dashed lines) for N = 3, 4 and 5.The cut-off angular frequency of each layer is ω c (cf. Eq. 6).Notice that some of the poles of the exact solution (cross marks) are at lower frequencies than ω c (C: coupled filters.U: uncoupled filters), which improves the damping for frequencies around ω c .For ω 10ω c the results are the same.The bottom plot shows the ratio between the coupled and uncoupled solutions.Right: | H(ω)| as a function of the frequency and the number of layers (in logarithmic scale).
In order to calculate the insulator transfer function, cf.Appendix A, the cut-off angular frequency needs to be well known to avoid errors in | H(ω)|.The errors in the cut-off frequency are due to the assumption that all the time constants of the shields, τ, are the same.The ratio between the transfer function considering a −10% relative error in ω c and the actual transfer function is shown in Fig. 2. The errors only become significant when large number of layers are used and for ω ≫ ω c .This will be used in Sec. 3 to asses the validity of the assumption that τ i j = τ.For completeness, we present here the transfer function in the case when the OR is included.The OR acts as an extra low-pass filter.However, its thermal resistance, θ OR , and thermal capacitance, C OR , can be, in general, different from the thermal shield one, i.e., τ = τ OR .The transfer function is -see Appendix A for details, where τ OR = θ OR C OR .If τ OR ≪ τ and C OR ≪ C, Eq. ( 9) is recovered and the effect of the resonator is negligible.Here N is the number of shields without including the resonator.However, in some cases C OR ≥ C since optical resonators are typically bulky bodies made of Zerodur or ULE.The thermal resistance (between the last layer of the thermal shield and the resonator) and thermal capacitance of the resonator are: where it is important to notice that θ N OR does not depend on the resonator's emissivity (assum- ing ε N ≪ ε OR and that r OR ≃ r N , where r OR is a representative radius of the optical resonator), ε OR , and for this reason it will usually be similar to θ .If such condition is met and the thermal capacitance of the shield and the resonator are also similar, Eq. ( 9) with an extra layer can be used instead of Eq. ( 11).

Conductive heat transfer
Thermal shield layers need, of course, mechanical supports between them.In this section, the effect of such a support structure on the performance of the thermal shields (2 layers) is analyzed.The supports or thermal links are supposed to have a very low thermal conductivity and, therefore, the lumped model approximation is not adequate.It is not guaranteed that the temperature is uniform along the supports [26][27][28].Consequently, the spatial gradient term of the Fourier heat transfer equation needs to be included.In this case, the transfer function of a shield is -see Appendix B for details: where κ s , c s , ρ s , A s and ℓ are the conductivity, the specific heat, the density, the cross section and the length of the supports.Note that the supports are assumed to have the same length and the cross-section is the sum of all of them.For the derivation of the transfer function it has been assumed that no radiative heat transfer occurs between the supports and the thermal shields.Therefore, the model is accurate if the supports are covered with high reflective coating.In Sec.3.2, a detailed numerical analysis is shown and in Sec. 4 they are compared with FEM simulations.

Concentric spheres and concentric cylinders thermal shields
Basically two geometries are used as thermal shields: concentric hollow spheres (or cubes) and concentric hollow cylinders (or rectangular cuboids).The cut-off angular frequency ω c of a thermal shield defines completely the attenuation of temperature fluctuations if we consider solely radiation.It is given in Eq. ( 6), where β depends on the geometry of the shields and the emissivity of the material [26]: , for spheres ( 16) where r is the radius.The view factors for finite cylinders, cubes and cuboids are more complicated than the ones given in Eqs. ( 16) and ( 17), however, the approximation of cubes to spheres and finite cylinders and cuboids to infinite cylinders yields accurate results -see Sec.3.2.If the shields are assumed to have approximately equal size (r i ≈ r j ) and high reflectivity materials are used (ε ≪ 1), Eqs. ( 16) and ( 17) reduce to β ≃ 2/ε and the cut-off angular frequency of a thermal shield is where h is the thickness of the layers.In Sec.2.1, we assumed that τ is the same for each shield.Equation (18) shows that this assumption is valid for the aforementioned conditions.However, in Eq. ( 18) two simplifications have been used: (i) the volume of the hollow sphere (or cylinder) has been approximated to 4πr 2 h (and similarly for the cylinder) and, (ii) the radius of all layers has been considered the same (r i = r j ).Both of these assumptions cause errors in the calculated cut-off angular frequency.To keep them below 10% -see Fig. 2, the ratio between the inner and outer radii of the layer i has to be r i,in /r i,out ≥0.9 for spheres (the layer thickness has to be h ≤ 0.1r out ) and r i,in /r i,out ≥0.8 for cylinders -see Fig. 3 (top); and the ratio between the radius of two consecutive layers has to be r i+1 /r i ≥0.9 for spheres and r i+1 /r i ≥0.8 for cylinderssee Fig. 3 (bottom).Both conditions, thin layers and small distance between the layers, are usually met when designing and constructing thermal shields.Fig. 3. Top: relative error in ω c due to approximating the volume of a layer to 4πr 2 h (or 2πrhℓ for cylinders).Bottom: relative error in ω c due to the assumption that the distance between two consecutive layers is zero (r i = r j ) when it is not.
The only assumption in Sec.2.1, which has not yet been justified, is the one concerning the temperature homogenization of the thermal shields.This matters, when a point-like heat source is present.A homogeneous temperature distribution can be assumed if the time constant of the shield itself is much smaller than the time constant of the radiative heat transfer, i.e., τ ≫ τ shield .The time constant of the shield is (for the spherical case) and considering thin layers (h = 0.1r): which is much smaller than one for high reflectivity and high conductivity materials.

Selection of the material
Given a requirement for the attenuation of temperature fluctuations, the thermal shields can be optimized in either mass or volume.In this section, this is briefly analyzed.For two thermal shields of different materials, A and B, having the same cut-off angular frequency the ratio between the thickness of the shields is (assuming both materials have the same emissivity, which can be easily done by coating the material) and the mass ratio is If we want to minimize the volume for a fixed mass of the layer (m A = m B ), a material with a low specific heat should be chosen following Eq.( 22).In contrast, if the mass should be minimized for a given radius of the layer (r A = r B ), a material with a high specific heat should be selected.

Numerical analysis
In this section, a numerical analysis of the equations shown in Sec.2.1 and Sec.2.2 is performed.First, Eq. ( 9) is used to calculate the required number of shields, N, for a given temperature stability requirement.If the attenuation is needed for ω ≥ 10ω c , the shields can be assumed uncoupled -see Fig. 1, and N is: where ⌈•⌉ is the ceiling function.Figure 4 shows the required number of shields for different attenuation levels: from 10 −7 to 10 −4 and for two frequencies: ω/2π=0.1 mHz and 1 mHz.The shields are assumed to be made of gold coated aluminum (ε=0.03) with a thickness of 0.1 mm, 0.5 mm and 1.5 mm, which correspond to cut-off angular frequencies of 60 µHz, 12 µHz and 4 µHz, respectively.Next, we design the supports between the layers such that they do not degrade the attenuation of the temperature fluctuations.To do so, we employ Eq. ( 14). Figure 5 (top) shows the attenuation at the radiative cut-off angular frequency -cf.Eq. ( 18), for two materials as a function of the length of the supports ℓ and the total cross section A s of the supports.For the calculations, the layer thickness has been set to h=0.5 mm and as material we chose aluminum.The distance between the layers ℓ is kept such that r i+1 /r i = 0.9 (the distance between the layers is 0.1r i ).The total cross-section of the supports has been constraint to A j /10 (the supports only cover 10% of the thermal shields area).For larger areas the model is not accurate since β is no longer of the form given in Eq. ( 16) -see Sec. 4.However, such large areas for the supports are unusual for typical designs of thermal shields.The properties of the support materials are given in Table 1.The materials are chosen due to their low thermal conductivity.
Figure 5 shows how the behavior of the system changes at a certain length ℓ min : For fixed ℓ (> ℓ min ), increasing the cross-section of the supports improves the attenuation of temperature fluctuations and for ℓ < ℓ min diminishes it.This turning point, ℓ min , occurs for τ s 2τ, where Table 1.Properties of the support materials used for the numerical evaluation in Fig. 5.
Ultem 1000 0.122 1280 2000 Macor 1.46 2520 790 and τ = θC.Therefore, the respective minimum length of the supports is: which depends only on the materials of the shields and supports as well as the thickness of the layers h.For instance, for h=0.5 mm we have ℓ min = 35.5 mm for Ultem 1000 and ℓ min = 150 mm for Macor.It is important to remark that once ℓ > ℓ min the support cross-sections do not have any constraints provided A s /A j < 1/10, which we assumed in the derivation.Note that it is not necessary that the shield layers are at that distance ℓ min , because the thermal supports can be routed conveniently between the shields allowing still a compact designs as long as the radiative heat transport between the supports and the layers is negligible.The latter assumption of our model can be met by using supports with high reflectivity coatings.This at first glance counter-intuitive change in behavior occurs since for a certain length increasing the area of the support causes an increase of the thermal capacitance, which compensates the reduction of the thermal resistance due to the larger cross-section of the support.The bottom left plot in Fig. 5 shows the attenuation at the radiative cut-off angular frequency ω c for a given A s as a function of the length of the supports for Ultem 1000 and Macor.The minimum length corresponds to that where the attenuation is 1/ √ 2, i.e., the length where the cut-off angular frequency of the thermal shield with supports does not change with respect to the one without supports.The minimum length for Ultem is much shorter than the one for Macor and also the attenuation drops much faster when increasing the length.However, in the asymptotic limit the attenuation for Macor is stronger than the Ultem one: The right bottom plot shows the transfer function for three scenarios: (i) ℓ = ℓ min /2, (ii) ℓ = 2ℓ min and (iii) radiation only.The rest of the parameters are kept the same.The change of behavior in the attenuation is again clear.It is also interesting that combining the thermal shields and the supports (radiation and conduction) adequately a dip in the transfer function can be obtained.Such behavior could be used to design thermal shields, where temperature fluctuations have to be attenuated in a very narrow frequency band.Bottom left: the attenuation of temperature fluctuations at ω c for different support lengths.We set in this calculations r i+1 /r i = 0.9 and A s /A j = 0.1.If ℓ < ℓ min the attenuation is degraded around ω c compared to a system without the supports (the attenuation at ω c is 1/ √ 2 -horizontal dashed line).Once ℓ > ℓ min , increasing the length of the supports improves the damping of temperature fluctuations.This improvement reaches an asymptotic value not too far from ℓ min .Bottom right: transfer functions as a function of the frequency for two scenarios with different support lengths and the radiative case only.
The results shown in Fig. 5 are for one thermal shield consisting of two layers.When nesting N shields the transfer function including the supports can be approximated by: This assumption of uncoupled layers yields errors similar to those shown in Fig. 1.The derivation of the solution for the coupled case is given in Appendix B.

Analytical solutions and FEM simulations results
In this section, we compare the results and properties found in the previous section with the results of FEM simulations performed with commercial software (Comsol).Figure 6 shows the transfer function, if we take only radiation into account.The model considered here consists of four aluminum layers (N=3) with ε=0.03 and h=1.5 mm.The radii of the shields are from the inner shield to the outer: 100 mm, 111.1 mm, 123.9 mm and 137.9 mm, respectively i.e., the ratio r i+1 /r i has been kept at approximately 0.9.The maximum mesh element size is 31 mm and the minimum one is 5.5 mm.The mesh settings are valid for the rest of the simulations with finer mesh when including the conductive thermal links.We use a sweep sine with frequencies from 1 µHz to 50 µHz as the boundary condition of the outermost layer of the thermal shield.
The transfer function is calculated as the ratio of the Fourier transformed temperature at the innermost and the outermost layer of thermal shield.Figure 6 shows that the numerical simulations are in very good agreement with the analytical results calculated using Eq. ( 9).The discrepancy between the simulations and the model is close to the expected one due to the fact of assuming τ i = τ i+1 when r i+1 /r i =0.9 -see Fig. 2 and Fig. 3.The simulations were carried out for different cross-sections with an oscillating temperature at ∼ ω c as the boundary condition in the outer aluminum layer.Three configurations have been simulated for a given ℓ=70 mm and different ℓ min : (i) ℓ min = 1.75ℓ=40 mm (blue trace), (ii) ℓ min = 0.58ℓ=120 mm (red trace) and (iii) ℓ min = 0.2ℓ = 350 mm.The results from the simulations ("×" marks) and from the analytical equations (dashed lines) are also in good agreement.The right panel in Fig. 7 shows the simulated and analytical transfer functions for different parameters, when the supports are included.The model considers only one shield and the properties have been chosen to simulate different expected behaviors.They are summarized in Table 2. Cubes have been used instead of spheres and, therefore, A j needs to be calculated accordingly since it appears in θ in Eq. ( 14).The results from the simulations and the analytical calculations are in very well agreement, too.Table 2. Properties of the supports used in Fig. 7 (right panel): ρ s =1200 kg m 3 , c s =7200 J kg −1 K −1 , ℓ = 0.053 m; properties of the layers of the thermal shield: aluminum with h=1.5 mm, ε=0.03,A 2 =0.06 m 2 , r 2 =0.1 m and r 1 =0.153 m (2r is the length of cube's edge).

Model
12.2 10 −3 (iv) 12.2 6.2 × 10 −3 .FEM simulations and analytical results.Left: temperature fluctuations attenuation at ω c for ℓ > ℓ min (blue) and ℓ < ℓ min . (red and black).ℓ =70 mm for all the simulations.Right: transfer functions for different models -see Table 2, and the ratios between the simulations and the analytical transfer functions, which indicate the good agreement between them.The simulations failed for frequencies higher than 100 µHz due to numerical resolution issues.
Finally, simulations have been performed to quantify two limitations of the analytical model and its assumptions: (i) the ratio of the cross-section of the supports and the area of the shields was assumed to be sufficiently small -see Sec.3.2 and, (ii) radiative heat transfer between the supports and the shields was neglected.If (i) is not met, the view factors used for the derivation in Sec.3.2 are not valid anymore.Figure 8 (left) shows the comparison between the simulations (solid lines) and the model (dashed lines).The results are in good agreement for ratios smaller than ∼0.15.The disrepancy at the dip for A s /A j = 0.15 is mainly due to the limited frequency resolution of the estimated transfer function from the simulations.For larger ratios (black trace), the errors become significant.For instance, for A s /A j = 0.66 the theoretical transfer function and the one obtained by FEM simulations differ by a factor of ∼ 10 at high frequencies.However, such large ratios are rarely used in thermal shield designs.The right panel in Fig. 8 compares the simulations when including radiative heat transfer between the supports and the shield layers.The simulations (solid and dashed lines) and the model (dashdotted lines) are in well agreement for small ε s .for different values of the supports cross sections and no radiation between the supports and the shields, ε s =0.The model agrees with the simulations when the total area of the supports is significantly smaller than the area of the shield (A s /A j 0.15) as previously stated.The bottom plot shows the ratio between the simulations and the analytical transfer functions.The discrepancies are significant for A s /A j =0.66 at high frequencies.Right: FEM simulations (solid and dashed lines) for different supports emissivities and A s /A j ratios and the theoretical transfer functions (dash-dotted lines) where ε s = 0.The difference between the simulations and the model appears for large emissivity values.The bottom plot shows the ratios between the numerical simulations transfer functions and the analytical model.

Conclusion
Highly sensitive interferometers and experiments testing fundamental physics require an isolation against external and environmental influences.In particular, the attenuation of external temperature fluctuations by several orders of magnitude at low frequencies (sub-milli-Hertz and milli-Hertz) using thermal shields is crucial.Such thermal shields are typically designed by numerical FEM simulations, which are time consuming and precision limited.Otherwise, simple analytical "toy" models that do not reflect all the features of the system are employed.However, an accurate analytical method is to our knowledge still lacking.We derived in this paper a fast and accurate method to calculate the performance of thermal shields in the frequency domain.This analysis has also yielded interesting results, such as that for some configurations a dip in the transfer function can be achieved in a narrow frequency range, and that the cross-section of the supports can be increased without jeopardizing the shields performance as far as they reach a critical minimum length.The former property is of interest for experiments with a narrow and well-known frequency range, such as BOOST.The latter can be useful for space experiments to design mechanically stable shields without loosing thermal insulating capability.The analytical model has been validated by comparing it with the results obtained by FEM simulations.This method allows for quick investigations of the design of thermal shields to optimize its performance, which later can be assessed and fine-tuned by FEM simulations.

Fig. 2 .
Fig. 2. Ratio between |H(ω, ω c )| and |H(ω, 0.9ω c )| as a function of N and angular fre- quency.The relative error of the cut-off angular frequency is −10%.The errors in the transfer function are tolerable even for large N if the relative error in ω c is kept smaller than −10%.

Fig. 4 .
Fig. 4. Number of shields needed as a function of the required attenuation for 0.1 mHz and 1 mHz and different shield thickness, h.

Fig. 5 .
Fig. 5. Transfer function for different support lengths and cross-sections for two materials and h=0.5 mm.Top left: Ultem 1000.Top right: Macor.The vertical dashed lines indicate ℓ min .Bottom left: the attenuation of temperature fluctuations at ω c for different support lengths.We set in this calculations r i+1 /r i = 0.9 and A s /A j = 0.1.If ℓ < ℓ min the attenuation is degraded around ω c compared to a system without the supports (the attenuation at ω c is 1/ √ 2 -horizontal dashed line).Once ℓ > ℓ min , increasing the length of the supports improves the damping of temperature fluctuations.This improvement reaches an asymptotic value not too far from ℓ min .Bottom right: transfer functions as a function of the frequency for two scenarios with different support lengths and the radiative case only.

Fig. 6 .
Fig. 6.Left: FEM simulation and theoretical results -cf.Eq. (9) for the transfer function of a thermal shield consisting of four aluminum layers (N=3) considering only radiation.The bottom panel shows the ratio between the results.Right: thermal shield model used for the calculations and simulations.The OR (12.5 cm) is not included in the simulations.The mesh is shown in the left figure inset.

Figure 7 (
Figure7(left panel) shows the simulations performed to evaluate the behavior found in Sec.3.2, when conductive elements are included in the model.Different configurations have been simulated with supports (70 mm in length) used to connect the thermal shield layers.The properties of the shields and the supports have been tweaked to achieve different ℓ min values.The simulations were carried out for different cross-sections with an oscillating temperature at ∼ ω c as the boundary condition in the outer aluminum layer.Three configurations have been simulated for a given ℓ=70 mm and different ℓ min : (i) ℓ min = 1.75ℓ=40 mm (blue trace), (ii) ℓ min = 0.58ℓ=120 mm (red trace) and (iii) ℓ min = 0.2ℓ = 350 mm.The results from the simulations ("×" marks) and from the analytical equations (dashed lines) are also in good agreement.The right panel in Fig.7shows the simulated and analytical transfer functions for different parameters, when the supports are included.The model considers only one shield and

Fig. 8 .
Fig.8.Left: FEM simulations (solid lines) and theoretical transfer functions (dashed lines) for different values of the supports cross sections and no radiation between the supports and the shields, ε s =0.The model agrees with the simulations when the total area of the supports is significantly smaller than the area of the shield (A s /A j 0.15) as previously stated.The bottom plot shows the ratio between the simulations and the analytical transfer functions.The discrepancies are significant for A s /A j =0.66 at high frequencies.Right: FEM simulations (solid and dashed lines) for different supports emissivities and A s /A j ratios and the theoretical transfer functions (dash-dotted lines) where ε s = 0.The difference between the simulations and the model appears for large emissivity values.The bottom plot shows the ratios between the numerical simulations transfer functions and the analytical model.