|
|
1.INTRODUCTIONSpaceborne reflective optics are usually made of materials having a very low thermal expansion coefficient around the targeted operating temperature of the space mission. The low dilatation property guarantees a form stability of the optic though all the thermal transients, however some other effects can contribute to the optic deformation such as compaction and or thermomechanical change induced by ionizing radiations. To date, the effect of the ionizing radiations on spaceborne optic material is not comprehensively understood, hence, for conservative reasons, the metrological error budgets associated to this phenomenon often overestimate the impact of radiations on the material. Studies carried out in the past where almost exclusively investigating the effect of a homogeneously distributed dose-deposition on some properties of optical material. Practically, the dose-deposition expected through the optical material varies strongly through the depth and depends on the type of mission considered (LEO/GEO, duration of the mission, local shielding of the optic due to the superstructure of the satellite, etc…). For all these reasons, the determination of a phenomenological ageing law describing the evolution of the characteristics of the material as a function of the local deposited ionizing dose is a very valuable tool in order to define a realistic error budget tailored to their specific mission parameters. In this present work, we investigated the influence of the electrons on the change of coefficient of thermal expansion of ZERODUR®. In the first section we present the analytical approach we used in order to describe the impact of the ionizing radiations on the vertex radius of curvature of the optic. In the second section are presented the experimentation that we carried out in order to characterize the change in the thermal properties induced by the ionizing radiation. The third section describe the expected changes in the surface peak to valley and vertex radius of curvature under realistic electron irradiations. 2.INDUCED CURVATURE. ANALYTICAL CALCULATIONIn this section, we present a textbook 1D analytical models based originally on the bimetallic assumption that we considered for the calculation of the irradiation induced curvature. This approach could appear obsolete compared to the powerful FEM solution routinely available nowadays; however, the typical dose-depth profiles that we are investigating in this work are difficult to address numerically: Typically, a decrease by several orders of magnitude of the deposited dose over micrometric scales are expected. These profiles would require a unusually tight FEM meshing of the optic in order to be comprehensively described at the macroscopic scale: such simulations could be achieved at the expense of exaggeratedly long simulation times. 2.1Bimetallic approachPrior works investigating the effect of total ionizing irradiation on the change of curvature of optics were mostly based on the bimetallic approach. This model supposes that the optical device under consideration is composed of a homogeneously irradiated domain (front face of the optic) connected to a homogeneously shielded domain (back face of the optic). The geometrical determination of the irradiated domain and therefore of the estimated average irradiation level over this range were both arbitrary, these approximations biased the calculation of the curvature (and to the contrary, biased the estimation of the compaction law that had been derived using the bimetallic curvature expression). Because of this, we did not use directly the bimetallic description in this work. 2.2Continuous approachA continuous approach is described in Opto-structural Analysis: John W. Pepi [1]. The calculation of the Vertex Radius of Curvature (VRoC) of a rotation symmetric optic is given by: With being the second moment of area, and d the thickness of the optical element. And ε(z) being the local strain defined: Where ρ0 correspond to the local density of the sample taken at a reference temperature T0, z represents the position along the longitudinal axis of the cylindrical coordinate system, see Figure 1, and ρ(z) is the local density taken at a temperature T(z) with ρ(z) = ρ(T(z)). Conventionally, we can use the center of the discs in order to perform the integration. However, if we subtract the average contribution of the strain to the strain profile function, we can then perform the integration using any possible coordinate origin along the sample thickness (this is especially useful if the origin of the coordinate system matches with the z=0 origin of our dose-deposition profile). The averaged contribution of the strain is defined: Consequently, the effective integrated strain contribution can formally be defined: This continuous approach is extremely simple compared to e.g. the multi-metallic equation [3], in spite of its simplicity the results generated using this equation are very consistent with the multi-metallic approach and also with FEM results provided that radius/thickness ratio remains in the validity domain of the thin plate approximation. 2.3Strain modelThe local strain ε over the depth z of the optical device under consideration can be mathematically defined: The first term on the right hand side corresponds to the radiation induced compaction, with D(z) represents the total ionization dose (TID) depth profile and ρ the local density of the material. This contribution is not discussed in the present paper as the validity of the presently available compaction models is still debated. This phenomenon will be further discussed in future work [4]. The second term on the right hand side of the aforementioned equation corresponds to the thermal strain induced by the change of Coefficient of Thermal Expansion (CTE) combined with the temperature gradient over the depth (remember that we are working presently with one-dimensional models). Two scenarios will be discussed in the present paper, we will first investigate the change of sagitta induced by the homogeneous variation of the temperature close to the operating temperature of the optic, i.e. T(z) = T ≠ Top., in this case the contribution T(z) – Top. is constant and the observed deformation are solely due to the local variation of the CTE induced by the ionizing radiation: Second, we will further complexity the model by superimposing a temperature field onto the CTE field: And we will check the form stability of different optical materials at different operating temperature. 3.RADIATION INDUCED CTE VARIATIONSThe estimation of the sagitta induced by spatial (depth) CTE variation presupposes the knowledge of the phenomenological law binding the local dose absorbed to the CTE variations. These results are only partially available in the literature, e.g. only for one given temperature or correspond to the temperature averaged CTE variation over a T=0°C to 50°C range. Prior to our calculations, we had to quantify the impact of ionizing radiation on the thermal properties of Zerodur®. In this section are presented the irradiation setup as well as the analytical CTE measurement performed on our coupons for the estimation of the dose-dependent CTE variations. 3.1Sample geometryDilatometric measurements made in Schott AG metrological center typically require cylindrical samples 6 mm in diameter and 100 mm long. However the intended homogeneous irradiation with dose variation smaller than 10% of a 100 mm long sample in the radiation field of Gaussian shape could only be achieved at the expense of exaggeratedly low achievable dose rate. This would require a more time consuming and therefore costly irradiation run. In order to reduce the irradiation time we divided the samples in pairs of two 50 mm long rods that can be afterward stacked together into the dilatometer holder. With these smaller samples we could achieve the best compromise between radiation field homogeneity and dose rate. 3.2Irradiation procedureThe irradiation has been performed at the research linear accelerator of the Metrological Electron Accelerator Facility (MELAF) of the Physikalisch-Technische Bundesanstalt (PTB). The energy of the electrons used for the irradiation was 5 MeV. The electron beam is pulsed with a typical pulse width of 2.5 μs. The pulse repetition frequency used during the irradiation was 50 Hz. An Integrating Current Transformer (ICT, Bergoz Instrumentation) is installed at the end of the accelerator beamline for non-destructive charge measurement of the electron pulses for beam monitoring. At the end of the accelerator beam line the collimated electron beam passes through a 0.05 mm thick Ti foil vacuum window which scatters the beam. The diameter of this window is much larger than the beam width, thus all electrons detected by the ICT contributed to the radiation field. The scattering in the Ti foil results in a radial electron radiation field with a Gaussian lateral profile. The width σ of the field increases with decreasing energy and increasing distance to the source. The width of the lateral dose distribution as function of the distance was determined experimentally by means of an ionization chamber matrix (PTW OCTAVIUS Detector 1500) as well by MC simulation (see Figure 2). The samples were irradiated on the central axis of the radiation field at a distance of 1100 mm. At this distance the width of the field is large enough (σ = 92 mm for a 5 MeV electron beam) to achieve a homogenous irradiation, i.e. a dose variation smaller than 10% over the surface to be irradiated (radius of 43 mm). At this distance, the achievable mean dose rate amounts to about 103 Gy.min-1. The samples were irradiated with 103, 104, 105, and 106 Gy, respectively. For the different deposited doses, the same setup and beam setting was used. Only the duration of irradiation was different. In order to achieve a homogeneous dose distribution also in the depth of the samples, they were irradiated from two sides. For this purpose, the samples were rotated after one half of the irradiation time. 3.3Dose determinationThe absorbed dose to Zerodur® per incident electron in the current setup was determined by means of the Monte Carlo simulation Toolkit GEANT4 [5]. The number of electrons incident on the glass ceramic samples was measured in situ by means of a dedicated Faraday cup (FC) connected to a current integrator. The FC acting also as a sample holder (see Figure 3). Since accurate measurement of short intense charge pulses is challenging, the validity of this charge measurement was verified using a second independent method. For this purpose, the charge of the beam pulses in the beamline was measured inductively using a current transformer and the proportion that impinge on the cross-sectional area of the FC was determined. The back side of the FC is a metal block, thick enough to completely stop high-energy electrons up to 10 MeV. Furthermore, it serves to dissipate the heat, in order to avoid a heat up of the samples during the irradiation. The front side of the FC is made of graphite to minimize electron backscattering from the metal block and loss of secondary electrons. In addition, the electrical and heat conductivity of this material avoids electrostatic charging as well as heat up of thesamples. The high-energy electrons pass through the graphite and the glass ceramic samples, depositing in it a part of their energy (dose) and are then stopped in the metal block at the backside. The FC is electrically isolated from the ground by means of a wood wedge. The leakage current was negligible, i.e. all electrons collected by the FC are counted by the current integrator (Ortec 439). A resistor-capacitor circuit (120 Ohm, 15 μF) was mounted at the input of the current integrator to keep the current from the short high charge pulses low enough. This in situ measurement setup is especially valuable for long time irradiation run were the intensity of the electron beam could drift or if some other disruptions would occur during the irradiation. To verify the charge measurement method the charge of the beam pulses in the beamline was measured simultaneously using FC and by means of the ICT. The fraction f of electrons that impinge on the cross-sectional area of the FC is a function of the FC radius R, the distance z, and the width σ of the Gaussian shaped field: The measured fraction of electrons collected by the FC agrees well with the theoretical value f (see Figure 4). ICT and FC give the same result when the FC is directly in front of the exit window, i.e. f(z=0)=1. It therefore can be concluded that the charge measurement methods by means of FC and ICT agree (QFC(z=1100 mm)=f(z=1100 mm) x QICT). The estimation of the deposited dose has also been achieved using Geant4. Practically simulations have been performed using a fixed number of electrons being directly fired onto the target NSim. and the corresponding simulated dose-depth within ZERODUR® has been determined (DSim.) Using these two numerical values, we could determine the amount of electrons N required in order to experimentally achieve an average target TID of 〈D〉 in our samples: In this above formula V represents the volume of a CTE sample. We can reformulate this equation into: To this amount of electrons N corresponds a charge q defined: Where C represents the elementary charge of an electron. The charge q is the one that must be collected on the Faraday cup in order to achieve the target average dose 〈D〉 The irradiation at a given target average TID has been systematically performed in two phases. During the first phase we irradiated the samples until we reached one half of the target charge on the Faraday cup, then we flipped the samples (front to back) and the irradiation run were then finished, i.e. till we reach the total targeted charge on the FC. The corresponding dose depth profile is exemplary pictured for a TID of 105 Gy on Figure 5. 3.4Temperature monitoringThe heat up (as a corollary of the dose deposition rate) within the sample holder should be minimized. From previous Electron Spin Resonance investigations on glass ceramics and also from photoluminescence measurements is known that some defect healing could be thermally activated [6].Thus, by keeping the temperature close to the room temperature we avoid potential temperature activated effects at higher temperature. A pyrometer was used for in situ monitoring of the temperature of the sample holder. The defects induced by the ionizing radiation in the samples may also partially reduced by thermally induced relaxation and healing during the time span between the irradiation and the CTE measurement. In order to freeze the damaged state of the samples, we conserved them at low temperature (-17°C) directly after irradiation. The samples have been transported from PTB to the Schott Marienborn metrological center using a dedicated container that kept the samples at low temperature (≤ -17°C), and once in Marienborn, the samples have also been kept by ≤ -17°C prior to analyses. 3.5Analytical measureAs with most mirror materials, ZERODUR® averaged CTE value over the T=[0;50] °C temperature range is expressed, so-called CTE(0;50). This figure or merit is not fully relevant for those designing and operating space optics. Frequently, the embedded mirrors are designed to be optimal around a given operating temperature, for example the [0;20] °C range. It is important to determine the instant CTE value at each temperature as well as the change induced by the radiation on the mirror in order to characterize with precision the surface level effect of local temperature variation around the operating temperature. The instant CTE, CTE(T), corresponds to the first derivative of the Total Change of Length (TCL = Δl/l = (l(T) – l0)/l0 with l0 the length of the samples at a given constant temperature T0) with respect to the temperature T. In this work, the TCL measurements have been achieved by means of an optical dilatometer, the measurements are performed dynamically over a [-30;30] °C temperature range, see Figure 6. Due to the smallness of the changes observed, several iterations were required in order to resolve properly the CTE values. To this end we used different measurement protocols based different time-temperature cycles in order to determine the best accuracy the CTE. For the determination of the measurement protocols, we systematically used different pair of ZERODUR® rods in order to perform CTE characterization on as irradiated samples. An exemplary t-T program is picture on Figure 7. Previous work made in this field by (Paige L. Higby [6]) suggested that the TCL measurement could eventually heal defect the samples especially during the TCL measurement over the higher temperature (at that time 300 K). Being aware of this difficulty, we purposely performed some of our first TCL characterization starting from a low temperature plateau -10°C followed by a further decrease of the temperature down to -20°C (at -1.2 K.min-1). This phase was then followed by an increase up to 10°C (at +1.2 K.min-1) then again followed by cooling down to -20°C (at _+1.2 K.min.-1). We repeated this program with increasing progressively the upper temperature bound. This procedure has been designed in order to check whether any healing effect could be observed (we would then observe a systematic drift of the CTE over the T-t iteration). Interestingly, in spite of having cautiously investigated any thermal relaxation threshold, we could not observe any healing effect similar to those reported by P.L. Higby in her pioneering work. The thermal activated healing process having been set aside we performed the standard TCL investigation using the time-temperature program described in Figure 7. 3.6CTE variation lawThe different TCL curves have been fitted over a [-20;30] °C by means of a degree 5 polynomial. These polynomials have been thereafter derived with respected to the temperature and the local derivative of the TCL, i.e. the CTE: The CTEs values calculated for different averaged deposited doses have been graphically presented on Figure 8 below: Using these data, we proposed a power law estimation of the local CTE variation as a function of the TID, ΔCTE(Top.; D) = ATop. × DBTop.. Following this formalism, the overall CTE can then be written: The power law has already successfully be used in order to characterize compaction [7][8]. Moreover, this formulation naturally takes into account the fact that for zero dose, the induced changes must also be zero. The offset between the 106Gy and 105Gy instant CTE change appears to be linearly temperature dependent, in order to quantify this visual impression we plotted the difference between these two set of CTE data: This above relation has been plotted on Figure 9: The differences between the CTE for 106 and 105 Gy confirmed to be reasonably well described by a linear regression law. Based on this observation we arbitrarily decided to pack all the temperature dependency into the prefactor ATop. and to set the exponent BTop. = B as being independent from the temperature. Numerically we estimated that the prefactor B is quite close to 0.5 and we decided for the convenience of the numerical implementation to fix it equal to 0.5. Following this ansatz the CTE change law expressed in ppm.K-1 as a function of the temperature in °C and the deposited dose expressed in Gy can then be formally written: As the variations of Zerodur® CTE are by far smaller than initially assumed, only the data determined for an average deposited dose of 106 and 105 Gy could be used for the derivation of a phenomenological law. The authors are aware of the sparsity of these data, however with this prospective work we are now able to propose further irradiation points, e.g. 2 × 105Gy, in order to consolidate this phenomenological description. 4.SIMULATIONS4.1FEM validation of the 1D analytical modelThe continuous 1D description of the surface curvature is valid as long as we are working in the realm of the thin plate approximation, i.e. some drifts could be expected between the FEM measurement and the 1D analytical model if the diameter/thickness ratio is not large enough (typically ≥10:1). An aspect ratio of 10:1 is typical of many open-back lightweight space borne mirrors up to 1.5m diameter delivering 200 Hz free-free first eigenfrequency. Larger lightweight mirrors may go to somewhat chunkier aspect ratios. As a preamble to our analyses, we wanted to cross-check the agreement between the 1D plate model and the FEM predictions. To this end, we used the FEM simulations model that we considered for defining the DOE of a ZERODUR® compaction law, see [4]. The meshing of the structure has been made by means of ANSYS and the calculation has been made with ABAQUS on a axial symmetric system with only one fixed point. The finite element size in the z-direction were adapted to the discretization step of the dose-depth, i.e. 0,05 mm, whereas their radial extension was about ten times this value. The samples under consideration are assumed to be irradiated with a 0,5 MeV and 1,5 MeV e- flux. Due to the different e-penetration depth at these energies we considered two different geometries. Coupons having a 2 mm thickness are irradiated with 0,5 MeV whereas coupons having a 5 mm thickness are irradiated with 1,5 MeV. The dose-depth profiles along the z-axis (assuming a polar representation of the coordinate system) are translated into a strain profile by means of a compaction law. The final change of curvature of the optic under consideration is then calculated by relaxing the stress field within the material. These FEM data have then been compared to the analytical results generated by means of the 1D bending approach. The numerical calculations have been achieved by means of a Python script based on the SciPy library. We used the same sample geometry as those considered for the FEM for the calculation of the VRoC (solely determined by the relaxation of the irradiation induced strain profile along the length of the optic as well as for the Peak-to Valley (PV) surface change VRoC that requires the VRoC as well as the radius of the optic for the calculation of the SAG). The numerical agreement is excellent see Figure 10. 4.2Dose-depth profile over a typical spaceborne mirrorDose-depth profile are frequently formulated in terms of dose deposited by an omnidirectional particle flux at the center of a solid sphere as a function of the sphere radius, see Figure 11 and [9]. The ionizing dose deposition calculated in this way is overestimating the dose deposition that one would expect on a real embedded mirror. In a satellite system, the ZERODUR® mirror would be shielded from the omnidirectional flux by a tube assembly consisting of the outer stray light baffle, the telescope’s metering structure plus thermal shields. In order to estimate the shielding influence of the telescope tube assembly onto the primary mirror, we designated a reference tube assembly 2010 mm in length of aluminum (2.54 mm thick). In actual telescopes, the effective shielding thickness typically includes substantially more shielding mass than our reference case. The primary mirror assumes a 10 mm thick ZERODUR® facesheet on the front of support pocketed to lightweight. For simplicity, we neglected shielding from the secondary mirror assembly (see Figure 12). Moreover, since there was miniscule penetration of radiation beneath the facesheet, we legitimately neglected the dose-depth profile within the isogrid structure of the telescope. It is assumed that there will be a strong back, instruments and a spacecraft aft of the primary mirror; thus, irradiation coming from the back will be negligible. In this work. The exclusive effect of e-.for the characterization of the CTE variation laws has been considered. Consequently, in those subsequent sections, we only focus on the contribution of e- on the optic deformation. The particle flux distribution for a given orbit can be calculated by means of Spenvis [10]. The particle trajectory incident on the telescope assembly is assumed being omnidirectional, with the isotropic e- fluence coming from a 2π steradians upon the telescope then onto the mirror. In order to test the consistency of our calculation we cross checked the results given by Spenvis (MULASSIS module) with the ones generated with our Geant4 based simulation software. In our computed approach the e- are fired from the surface of a fictitious sphere centered toward our optic (the radius of the sphere is assumed large compared to our telescope design in order to mimic at best the omnidirectional flux), see Figure 13. Cross comparisons made with Spenvis on an unshielded system (i.e. without tube assembly) using the same electronic fluxes demonstrated consistency with the Geant4 simulations. See Figure 14 and Figure 15 for, respectively, for LEO and GEO orbits. The same simulations have been repeated taking into account this time the effect of the telescope tube assembly. In both cases, we observe a tube screening effect that reduces by ~15 to 20 times the magnitude of the dose deposited in the first millimeters of the mirror’s facesheet1. The radial dose profile is also extremely instructive, see Figure 16, we can see these that the radiation profile is almost flat, which means that we are very close to the working hypothesis of the 1D continuous model presented at the very beginning of this paper, where the dose distribution only depends on the depth. The deformation expected for those radius independent profile are expected to have a pure power contribution. The power deformation can easily be compensated by refocusing the optic. 4.31D estimated change of curvatureThe TID profile generated on the shielded system then have been used in the continuous 1D approach. We still considered a Ø=1m telescopes but because we are estimating bending, we take into account the isogrid structure of a 75 mm bulk mirror thickness (empirically found to be equivalent to 100 mm isogrid thickness). We first consider the telescope subject to a 10 years e- irradiation over two different orbits, namely LEO and GEO, with and without shielding. We next evaluate the effect of the TID on the local CTE and its consequence if the temperature of the mirror (isothermal assumption) is drifting by about ΔT = 5 K from the nominal operating temperature Top.: The results of this model are presented on the Figure 17 below: We can see that the impact of electrons is only marginal; moreover, one should keep in mind that these deformations are mostly defocus2 and can be easily mechanically corrected. Practically the reflective optic is also bound to the satellite structure; these connecting points are increasing the stiffness of the overall optic thus reducing further the deformations. In a second time, we showed the effect of a local temperature gradient of about δT = 0,75 K along the d=75 mm thickness of mirror onto the VRoC and PV (i.e 1 K over a 100 mm thick isogrid system), we performed these calculations taking also into account the change of CTE induced by the TID. In this configuration, the strain can be written: The calculation have been performed assuming two different materials, on one hand CO720® Cordierite material from Kyocera (unirradiated) and on the other hand ZERODUR® (unirradiated and irradiated assuming LEO/GEO orbits). The instant CTE of these two material are compared on Figure 18. The VRoC and the surface PV are provided on Figure 19. These calculations show that a gradient of temperature through the thickness has a by far more significant effect than a homogeneous temperature drift within an irradiation aged material. In this situation, the effect of the irradiation is almost negligible. This effect has been graphically represented on Figure 20. We can see that the CO720® is as good a Zerodur® very close to it optimal operating temperature (around 23 °C), but small variations of the operating temperature from the CTE null-crossing point already imply significant thermal strain and consequently some significant changes in VRoC. Practically, it is unusual and difficult to drive a telescope to isothermal as high as T ≈23 °C. Generally, steady-state gradients are present and acceptable in the baffles and structures. More complex are asymmetric transients introduced when the pointing vector changes with respect to the sun, or in the case of a earth observing system, for which solar view factors modulate with orbital phase. Ideally, the primary mirror is maintained under stable thermal conditions, not undergoing time-varying gradients. With time-varying gradients, performance is affected whether the material is high diffusivity, high CTE or low diffusivity, low CTE. These efforts include insulation, and in extreme cases, proportional heat radiated on the mirror. If a telescope primary mirror is to operate as warm as Room Temperature (~23-C), heat must be radiatively added to the back of the mirror. If the front of the telescope is looking at cold space, the face of the mirror will want to cool as the mirror and baffles look at this cold space, thus a gradient is present, with the principal thermal resistance being the reflective coating of the mirror. 5.CONCLUSIONIn this paper, we estimated the expected change of sagitta due to ionizing irradiation effect onto the CTE of Zerodur®. We first determined at the lab scale the impact of electrons onto the CTE for different TID in order to propose a phenomenological law. We then investigated the expected irradiation profile within our material over a LEO and a GEO orbit without and with baffle shielding. We then combined these dose-depth scenarios with our CTE ageing law in order to investigate by means of a continuous plate theory two scenario. In the first scenario we considered a worn Zerodur® facesheet under a spatially homogeneous temperature drift, the observed cambra is then due to the local change in the CTE caused by the irradiation field. In this approach, the change in surface PV is only very marginal. In the second scenario we introduced some more complexity by adding a thermal gradient through the thickness of the Zerodur® facesheet. These calculations showed once more that the effect of irradiation are only contributing marginally to the change observed in the VRoC. These calculations are all based on a 1D model whose validity has been previously cross-checked with FEM. However these results are assuming a free-free boundary conditions and are not taking into account any radial temperature distribution nor any complexity due to the isogrid structure of the reflective optic. In addition to this the CTE ageing law proposed in this work has been derived with a small dataset, the CTE being only marginally impacted by the ionization we only had two set of linear independent results which were significantly drifting from the baseline for feeding the fitting procedure. The generation of another set of data for a deposited dose of e.g. about 107 Gy would significantly help us in confirming the power law ansatz used for the description of the CTE changes. All these aforementioned aspects must be and will be further deepened in some future contributions. REFERENCESDöhring T., Hartmann, P., Lentes, F-T, Jedamzik, R, Davis M. J.,
in Proc. SPIE (7425),
(2009). Google Scholar
Pepi, J. W.,
in SPIE Press Bellingham,
(2018). Google Scholar
Vasudevan, M., Johnson, W., Appl. Sci. Res. B, 9 420
–430
(1961). https://doi.org/10.1007/BF02921840 Google Scholar
Carré, A., Hull, T., Valente, M., Krieg, J., Westerhoff, T., ICSO 2020 Paper session, 249 Google Scholar
Allison, J., Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 835 186
–225
(2016). https://doi.org/10.1016/j.nima.2016.06.125 Google Scholar
Higby, P. L., Friebele, E. J., Shaw, C. M., Rajaram M., Graham, E. K., Kinser, D. L., Wolff, E. G., J. Am, Ceram. Soc, 71 796
–802
(1988). https://doi.org/10.1111/jace.1988.71.issue-9 Google Scholar
Bourrieau, J., Roméro, M.,
“ESA Spacecraft Mater,”
In Space Environment, 275
–285
(1979). Google Scholar
Bourrieau, J.,
“Étude expérimentale sur l´irradiation par les particules à haute energie des miroirs utilisés dans le télescope radiomètre METEOSAT,”
Report N°2145/74/CG,
(1974). Google Scholar
Seltzer, S., National Bureau of Standards Technical Note 1116, U.S. Gov. Printing Office1980). Google Scholar
– “SPENVIS, Space Environment Information System,”
https://www.spenvis.oma.be/credits.php Google Scholar
|