Next Article in Journal
Non-Symmetric Gyroscope Skewed Pyramids
Next Article in Special Issue
Validation of a Test Platform to Qualify Miniaturized Electric Propulsion Systems
Previous Article in Journal
A Test-Bed For Measuring UAS Servo Reliability
Previous Article in Special Issue
Solar Module Integrated Converters as Power Generator in Small Spacecrafts: Design and Verification Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Numerical Analysis of the Primary Composite Structure of a CubeSat

1
National Council of Science and Technology (CONACYT)—National Center for Aeronautics Technologies—Center for Engineering and Industrial Development (CENTA-CIDESI), Colón, Querétaro 76270, Mexico
2
National Center for Aeronautics Technologies (CENTA), Colón, Querétaro 76270, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Aerospace 2019, 6(9), 97; https://doi.org/10.3390/aerospace6090097
Submission received: 7 June 2019 / Revised: 12 August 2019 / Accepted: 26 August 2019 / Published: 4 September 2019
(This article belongs to the Special Issue Verification Approaches for Nano- and Micro-Satellites)

Abstract

:
A thermal computational analysis for the composite structure of a CubeSat is presented. The main purpose of this investigation is to study the thermal performance of carbon fibre/epoxy resin composite materials with Zinc Oxide nanoparticles in order to be used in the panels of the primary structure of a CubeSat. The radiative heat fluxes over each composite panel are computed according to the orbit trajectory and they are utilized as boundary conditions for the analysis. The direct solar, albedo and Earth infrared radiation fluxes are considered in this study. The model implementation, including the computation of the orthotropic thermal conductivity of the composite material is presented. The thermal simulations were performed for three different orbit inclination angles: the selected mission ( β = 57 ), the worst hot ( β = 90 ) and the worst cold ( β = 0 ). The temperature ranges in the electronic boards are analyzed in order to show that are into the operating limits of each electronic component.

1. Introduction

A single unit (U) CubeSat is a satellite with 10 × 10 × 10 cm dimensions and mass up to 1.33 kg. These space platforms have two main advantages, lower cost, time of design and manufacturing development compared with the traditional systems [1]. Nevertheless, weight restrictions are one of the most critical issues during the design process of a CubeSat. Traditionally, CubeSats have been manufactured using aluminum alloys, the inherent reduction in weight utilizing composites offers an increasing on the payload. However, if any new material is intended to be used in space applications, it must be evaluated in order to validate its thermal and mechanical performance in the environmental and operating conditions to which it will be exposed in space [2]. Several efforts have been developed in order to design composite materials for the structure of a nanosatellite [3,4,5].
Some other efforts to develop nanosatellites using composite materials include the 3D printed polymer based CubeSats in which additive manufacturing technologies are utilized in order to have rapid prototyping and embed part of the electronics into the 3D printed structure during the manufacturing process [6,7]. A successful implementation of composites into a secondary CubeSat structure was developed in the SwissCube project [8], which is the first satellite entirely built in Switzerland. The Swiss CubeSat was manufactured using an Al alloy frame with glass-epoxy composite panels cured by autoclave. Main mechanical loads are supported by the Al frame and composite panels are used as a secondary structure. Thermal loads are restricted in the glass-fiber panels due to the presence of solar cells bonded at their top surface. The SwissCube was launched in 2013 and still continues orbiting the Earth.
Currently, the incorporation of nanoparticles as polymer reinforcement has shown huge potential for providing enhancement performances for different applications at very low filler content (lower than 1 wt %) [9,10]. In recent works, they allow an increase of the Glass transition temperature ( T g ) of polymers and/or a shift of their continuous working temperature towards higher values [11]. Composite materials added with different kind of nanoparticles can be turn into a solution to implement composites for aerospace applications, specially, when the environmental conditions require high performance behavior from the thermal, mechanical or UV resistant points of view.
The CubeSat design development requires predicting the temperature ranges of the CubeSat structure at the space environmental conditions in order to evaluate whether or not the electronic devices and payload remain into their operational and safety limits [1]. In Reference [12], the authors presented a thermal analysis for the STEP Cube Lab (Cube Laboratory for Space Technology Experimental Projects), a CubeSat framework for educational purposes. The analysis were performed using SINDA/FLUINT Thermal Desktop tool. The radiation fluxes and exchange factors were computed using RadCAD [13]. The temperature behavior as function of time is reported for the external panels and the main internal components of the CubeSat, showing that all the electronic systems satisfy the temperature requirements.
A thermal analysis for a 3U CubeSat was carried out in Reference [14]. Based on the orbit mechanics results, the authors computed the direct solar and albedo radiative heat fluxes as function of the CubeSat position along its orbit. The transient thermal model was implemented in ANSYS 13 software and, minimum and maximum temperatures reached by the CubeSat are presented considering only the aluminum structure and a multi-layer insulation throughout the aluminum structure. A similar but more complete analysis was done by Reference [15] for the TRIO-CINEMA (Triplet Ionospheric Observatory - CubeSat for Ion, Neutral, Electron, MAgnetic fields). The authors took into account the optical properties of the main components of the CubeSat to compute the radiation heat fluxes. Using the computational results, a passive thermal control was design in order to comply the temperature requirements for each electronic device. For the case of composite structures, in Reference [3], the authors compared the thermal performance of aluminum and composite materials for satellite structures. They also presented the possibility to use composite materials with thermal conductivity larger than the aluminum. As the authors commented, this can help to remove the excess of heat from the electronic devices reducing the thermal control cost. Their results showed that the use of this high conductivity composite materials can significantly decrease the temperature gradients along the satellite structure.
In the scope of the AEM-CONACYT (Mexican Aerospace Agency) project entitled, “Design and characterization of composite materials for CubeSat structure”, the research objectives are the following: (1) To investigate a family of low-cost composite materials that can handle space environments and (2) To evaluate the mechanical-thermal and radiation performance of the down-selected composite materials in a CubeSat prototype. At this point, it is important to comment that actual space composite systems are based in cyanate-ester resins. Cyanate-ester resins have as their main advantage a high glass transition temperature ( T g = 400 °C), low dielectric constant, thermal stability and low flammability, which makes it an excellent candidate for electronic and spatial applications [16]. Nonetheless, the cost of cyanate-ester resin is more than 6 times (2000 USD/gal) [17] compared with a low flammable FAR 25 [18] certified epoxy resin (300 USD/gal) [19]. Recent studies have compared both systems in terms of flame retardancy and thermal properties [20,21]. Conclusions have established that epoxy systems added with fire retardant and/or ceramic nanofillers improve their thermal behavior and can be considered as a candidate for spatial applications [22]. Having in mind those references, an epoxy resin was selected considering the following aspects: (1) mechanical and thermal properties, (2) compatibility with carbon fibres and (3) cost-benefit.
At the actual stage of the project, the family of promising composite materials has been developed and tested at coupon level. The purpose of this study is to communicate the thermal performance of the down-selected composite material considering the CubeSat structure in a possible space environment.
In the present work, a thermal analysis is implemented to evaluate the performance of the proposed CubeSat composite structure. The design configuration of our CubeSat employs the composite panels as main structure, because no Al main frame is considered. A Carbon fibre/Epoxy resin composite with Zinc Oxide nanoparticles was used for the CubeSat panels. Additionally, we propose this Out-of-Autoclave (OoA) manufacturing process, such as Resin Transfer Molding (RTM), due to the novelty trends respecting environmental agreements such as Horizon 2020 and Europe Horizon [23]. The principal reason for using a composite material is the reduction of weight with respect to the classic aluminum structure used in most of the CubeSats that are currently operating. During the first stage of the project it was found that the addition of nanoparticles increase the resistance to UV radiation with respect to the original composite material [24], this effect has been found for different composite and coatings materials [25,26]. Also, the characterization of this composite demonstrated that the presence of the Zinc Oxide nanoparticles (ZnO) into the polymeric matrix slightly improve its thermal properties [27].

2. Materials and Methods

2.1. Conceptual Design

The proposed primary CubeSat structure consists in six quasi-isotropic composite panels as shown in Figure 1a. The conceptual model of the structure is simplified with respect to a final design of a CubeSat. However, at this stage of the design process, only the evaluation of the composite material with nanoparticles was performed under the severe thermal conditions along a Low-Earth-Orbit. With this in mind, all the other external components such as solar cells are not considered in this study. The main objective is to demonstrate that the composite panels can provide proper thermal operating conditions for the electronic devices during the satellite mission. In Figure 1a, the panels nomenclature with respect to the coordinate reference frame is shown. The internal electronic devices configuration of the satellite is presented in Figure 1b. Only the relevant electronic devices from the thermal point of view were taken into account for the analysis (Electric Power Supply EPS, GPS, communication boards and battery pack, see Table 1).
The thermal conditions to which the structure of the satellite will be exposed are strongly dependent on the orbital parameters. Then, a specific mission was defined to estimate these orbital parameters and consequently the time-dependent radiative heat fluxes over the composite panels of the CubeSat. An orbit inclination angle β = 57 and an altitude of 250 km (Low-Earth-Orbit) were used as input parameters. From this information and assuming a circular orbit, the period, angular radius, sunlight and eclipse stages were determined.
A nadir-pointing scenario was selected (see Figure 1c) to evaluate the thermal performance of the composite panels. Nadir-pointing is a case in which the same face of the CubeSat is always pointing to the centre of the Earth. This scenario allows for the four faces exposed to direct solar heat flux to exhibit a sinusoidal behavior of radiation [28]. The proposed composite material for the primary structure panels is a quasi-isotropic [0/ 45 2 / 0 ¯ 3 ] S carbon fibre/epoxy resin laminate with 3 mm in thickness. Resin Transfer Molding (RTM) process was implemented to manufacture them. Before the RTM process, the epoxy resin was added with 0.25 wt % of Zinc Oxide nanoparticles. The Zinc Oxide nanoparticles decreased the degradation of the material when it is subjected to UV radiation. However, this effect is outside of the scope of the present investigation and is therefore not quantified. In order to characterize the thermal properties of the composite material, thermogravimetric (TGA) and transient thermal response analysis were performed and the thermal conductivity, diffusivity and heat capacity were calculated [27]. In Table 2, the thermal properties of the epoxy resin/ZnO nanoparticles and the carbon fibre are shown, the properties of the carbon fibre were taken from Reference [29].
The appropriate thermal design of the CubeSat can allow to maintain the on-board equipment within their operating temperature limits. These operating limits are shown in Table 1 for the most relevant electronic devices. As it was commented, the motivation to change from the classical aluminum to a composite primary structure is the weight reduction that can be achieved. In this sense, the weight of the proposed composite structure is 226 g. A comparable structure with solid-walls configuration and manufactured in Al 5052-H32 sheet with a thickness of 0.080’ (2.032 mm) weighs around 324 g. This means that using the composite panels, a 30 percent reduction of the weight of the CubeSat structure can be obtained.

2.2. Heat Radiation Sources

The environmental conditions to which the CubeSat is exposed are mainly dictated by the orbital parameters and the external radiation sources as the direct solar ( q s ), albedo ( q a ) and Earth infrared ( q i ) radiation. It is important to consider the specific translation and rotation motion of the CubeSat along its orbit, since this effect causes that the satellite receives heat on one side and dissipate heat through the remaining side, promoting a temperature gradient between the panels of the CubeSat.
The computation of the radiative heat fluxes is the most critical issue in the thermal analysis for a nanosatellite. The direct solar radiation depends on the solar constant ( C s ) defined as the intensity of sunlight radiation perpendicular to a surface at the mean distance from Earth to the Sun. For this analysis an average solar constant equal to 1367 W/m2 was used, since a circular orbit assumption for the CubeSat was established. Due to the cubic shape of the nanosatellite, there is no coupling between faces, it means that the rays reflected by one face of the CubeSat never falls upon other. The direct solar flux is computed as the projection of the solar vector to the corresponding face normal:
q s = C s cos γ i
where γ i is the angle between solar vector and the normal vector to the i CubeSat face. The albedo radiation flux is quite more difficult to compute, since it depends on many parameters. Some guidelines for the selection of the parameters used to compute albedo heat flux can be found in Reference [30]. Following a simple formulation presented in Reference [31], the albedo flux can be computed as function of the albedo reflectivity coefficient (a), the sun-earth-satellite angle ( θ ) and the so-called view factor ( F E ) as:
q a = C s a ( cos 0.9 θ ) e F E
The view factor of each face can be easily estimated as:
F E = r 2.1 ( sin ρ / 2 ) e
where r = R / ( R + h ) , being R the Earth radius and h the mission altitude and the exponent e is computed as function of r through a polynomial expression that can be found in Reference [31]. This formulation is an alternative approach to the chart presented in Reference [32] for the albedo heat flux and the view factor.
For the Earth infrared heat flux an average value E B B = 240 W/m2 can be used multiplied by the corresponding view factor of the i face.
Additionally to the external radiative heat fluxes, the face radiation to the space must be calculated. The standard environmental conditions used for a spacecraft thermal analysis considers a temperature of 5 K. The CubeSat will radiates energy to the space and it can be computed through the Stephan-Boltzmann radiation law. Finally, the internal power generation from the electronic devices is modelled as a volumetric heat source in each one of the electronic boards of the CubeSat. In Table 1, the heat dissipation from the most relevant electronic devices is listed.
A python script was developed in order to compute the orbital parameters and radiative heat fluxes as function of the mission definition. In Figure 2, the total radiative heat fluxes ( q = q s + q a + q i ) on each CubeSat panel are shown for one complete orbit. As can be observed, the fluxes exhibit almost a sinusoidal behavior. It is also illustrated the eclipse period of the nanosatellite orbit. Due to the inclination angle, the panel 6 (P6) receives the largest amount of direct solar radiation and the nadir (P2) the largest amount of albedo and Earth infrared radiation. The large difference between the incident energy over each composite panel will produce temperature gradients that can affect the integrity of the electronic devices and compromise the correct operation of the electronic boards. From the solution of the heat transfer model it can be determined if the composite material with nanoparticles will be able to maintain the temperature range into the safety conditions of the electronic systems in the CubeSat.

2.3. Model Implementation

A time dependent heat transfer model was implemented in ANSYS 19 [33]. The mathematical model is based on the energy conservation equation [1]:
T t = α 2 T + q ˙ ρ c p
where α = k / ρ c p is the thermal diffusivity, k is the thermal conductivity, ρ is the density, c p is the heat capacity of the material and q ˙ is the heat dissipation from the internal electronic devices. Due to the orthotropic nature of the composite materials, the thermal conductivity must be evaluated in each direction. The orthotropic thermal conductivity must be referenced by a local coordinate system ( x , y , z ) on each face of the CubeSat at which the z -axis is always on the thickness direction of the composite panel. For the in-plane thermal conductivity ( k x , y ) and the heat capacity, the simple rule of mixtures was used:
k x , y = k f V f + k m ( 1 V f )
where k f and k m are the thermal conductivities of the carbon fibre and epoxy matrix with ZnO nanoparticles, respectively, and V f is the fibre volume fraction. For the computation of the thermal conductivity in the thickness direction of the composite laminate ( k z ), the Clayton model was used [34]:
k z = K m 4 ( 1 V f ) 2 K f K m 1 2 + 4 K f K m ( 1 V f ) K f K m 1 2
The surface to surface radiation in the interior of the CubeSat was also modeled using the S2S model implemented in ANSYS 19. This model allows to simulate thermal radiation exchange between surfaces forming an enclosure. The medium that fills the space between the surfaces is non-participating, it means that the medium inside the CubeSat does not emit, absorb or scatter radiative energy. Therefore, the amount of radiation received and emitted by each surface is only defined by their optical properties. The radiation model S2S from ANSYS assumes the surfaces to be gray and diffuse, then, the emissivity ( ϵ ) is equal to the absorptivity ( α ). In the computational simulations the emissivity for the composite panels was taken as 0.81 [35] and for the PCBs (FR4) equal to 0.6 [15].
A hexahedral mesh was built to discretize the CubeSat structure including the composite panels and the internal electronic boards. The electronic boards were model using the FR4 material properties as in Reference [12]. In Figure 3 a total mesh with 52,400 elements is shown. For the composite panels a mesh with 10 elements along the thickness was used in order to capture the effect due to the orthotropic conductivity of the composite material. The total radiative heat fluxes were applied as boundary conditions on each composite panel (see Figure 3). The solution of the model was done for five orbits in order to minimize the initial conditions effect and compute the periodic behavior of the temperature along the orbit of the CubeSat.

3. Results and Discussion

The average temperature as function of time was computed for each panel. The simulations were performed for the defined mission with an inclination angle of 57 . Additionally, in order to show the feasibility for using this composite material for any orbital parameters, the radiation fluxes were computed for the worst cold and hot cases ( β = 0 , 90 ) and the thermal analysis were also carried out for these both conditions. In Figure 4, the temperature field in each panel for the case of β = 57 at t = 19,620 s (the coldest instant of time during eclipse period of the orbit) is shown. As can be observed, the panel 5 is the coldest one, since no albedo and Earth infrared radiation fluxes fall on this panel along the orbit. The temperature fields on the CubeSat composite panels is shown in Figure 5 for an instant of time during the sunlight period of the orbit. Due to the inclination angle of the orbit, the maximum temperature appears in P6, while in P1 and P3 the temperature falls below 263 K (−10 °C). At this point, it is important to discuss the results in terms of the temperatures ranges at which the electronic devices into the CubeSat can operate (Table 1). The average temperature in each panel as function of time is plotted in Figure 6 during five complete orbits. The minimum and maximum temperatures occurs periodically at the eclipse and sunlight stages of the orbit. The results show that the average temperature in each panel never falls below the lower limit of operation of the CubeSat electronic systems. Due to the selected orbit, P1 is the coldest one, since this panel is not exposed to the direct solar radiation along the whole orbit. On the other hand, panel P6 is the hottest reaching a temperature upon 340 K (67 °C).
The electronic boards temperature fields at t = 19,620 s and t = 22,856 s are shown in Figure 7. It can be observed that the temperature in each electronic board for the hottest and coldest instants along the orbit, depends mostly on the temperature gradient that appears over the external composite panels. In order to analyze quantitatively the temperature in each electronic board, the average temperature of each one as function of time is plotted in Figure 8. As can be seen, the temperature in the electronic boards follows almost the same behavior as the temperature in the composite panels. The temperature for each one of the electronic boards at any instant of time are into the range of their safety operating limits, showing the feasibility to use the composite panels as the primary structure for a CubeSat that perform the selected mission.
To compare the thermal behavior between the composite material and an aluminum structure, the simulation was carried out with the properties of the aluminum 5052-H32 sheet with a thickness of 0.080’ (2.032 mm) black paint, the thermal properties are listed in Table 2 and the emissivity of this alloy was taken as 0.89 [15]. The average temperature in each panel as function of time is plotted in Figure 9. As can be observed, the instantaneous average temperatures in all the aluminum panels are almost the same and the temperature range is similar to the exhibited in the case of the composite structure.
It is important to highlight that from the thermal point of view, the aluminum structure is better than the composite one, since its high thermal conductivity allows to reduce the temperature gradients in the CubeSat, however, the composite structure complies with the thermal requirements and the reduction of weight due to the use of a composite material makes it an attractive option for applications where the payload needs to be maximized.
The thermal validation for the conceptual design of the composite structure requires the evaluation of the behavior at the worst hot and cold orbit conditions. The coldest orbit occurs when the CubeSat spends the longest time eclipsed by the Earth ( β = 0 ) and the electronic devices operate at the minimum power. The hottest orbit case occurs when the eclipse time is minimized ( β = 90 ). For both cases the radiation fluxes were computed and incorporated in the heat transfer model. The average temperatures in each panel as function of time are plotted in Figure 10 and Figure 11 for worst hot and cold orbit cases, respectively. For the worst hot case, the panel 5 receives perpendicularly the direct solar heat flux along the entire orbit reaching temperatures around 350 K (77 °C), on the other hand, P1 does not receive any radiative heat flux and its minimum temperature is about 263 K (−10 °C), the large temperature difference between the hottest and coldest composite panels can promote thermo-mechanical stresses due to the thermal expansion of the material that must be studied during the structural analysis of the CubeSat. Nonetheless, as a quasi-isotropic configuration was chosen, it can be assumed that their thermal expansion will be the same in the in-plane principal axes of the panels (0, 90, ±45). For the thermal analysis purposes, the important issue is that the P5 reaches a temperature that is very close to the operating limits of the electronic devices. For the worst cold case (Figure 11), the temperature of the panels present a periodic behavior in a range of temperatures from −245 to 325 K. The minimum temperature in the composite panels for the worst coldest case is lower than the operating limits of batteries (−20 °C) and close to the communication board limit (−30 °C). In order to wrap up the analysis, the average temperature in the electronic boards for the worst hot and cold cases are plotted in Figure 12 and Figure 13. For the worst hot case the temperature range limits for all the electronic boards are satisfied. The average temperature in the electronic boards for the worst cold case has a range about 70 K during a complete orbit (see Figure 13). The issue in this particular case is that the temperature in the electronic boards falls below 255 K (−18 °C), compromising the operation of nanosatellite batteries. Then, for an orbit with an inclination angle close to zero, inner heaters will be necessary to ensure the correct behavior of all the CubeSat electronic systems.

4. Conclusions

The thermal analysis for the primary composite structure of a CubeSat was carried out. The main idea was to demonstrate from the thermal point of view, the feasibility for using a carbon based composite material with ZnO nanoparticles for the primary structure of a CubeSat. A computational analysis was performed solving the transient energy conservation equation implemented in ANSYS 19. The orthotropic nature of the composite material was included in the solution using Clayton model for the thermal conductivity calculation. The direct solar, albedo and Earth infrared radiative heat fluxes were computed as function of the orbit parameters and the position of each composite panel along the orbit trajectory. For the selected orbit inclination angle ( β = 57 ) and the worst hot orbit ( β = 90 ), the temperature ranges of the composite panels and the internal boards met with the safety operating limits of the electronic systems. However, for the worst cold orbit ( β = 0 ) the minimum temperature reached is very close to the lowest operational limit of the battery pack (−18 °C) and additional thermal control system must be incorporated to the CubeSat design. It is important to note that the hottest temperatures and temperature gradients between the panels reached by the CubeSat are higher than the temperatures for the aluminum structures computed (see Figure 9 and reported in the literature [1,12,15]. The main reason is the lower thermal conductivity of the composite material that reduces the heat transfer in the structure of the CubeSat and increases the temperature gradients between the panels. Although, with the results presented in this work, it has been verified that the thermal design of the CubeSat structure using the proposed material satisfies the temperature requirement of on-board electronic equipment. In addition to the thermal analysis presented in this work, static and dynamic structural analyzes are required to complement the feasibility study of the use of the proposed composite material for the main structure of a CubeSat. Specifically, it must carry out vibrational analysis to validate the mechanical performance for the material and establish the compliance of the CubeSat Standard [2].

Author Contributions

Conceptualization, S.P. and M.T.; methodology, S.P., M.T., S.L.; software, S.P.; formal analysis, S.P., M.T. and S.L.; investigation, M.T., S.P.; resources, M.T.; writing—original draft preparation, S.P., M.T. and S.L.; writing—review and editing, M.T. and S.P.; visualization, S.P.; project administration, M.T.; funding acquisition, M.T.

Funding

This research was funded by CONACYT through the Research Fund for the Development of Space Activities (FIDAE), grant number 275783.

Acknowledgments

Saul Piedra and Mauricio Torres, convey their special appreciation to the “CONACYT Research Fellow Program (Catedras CONACYT)”. All of the authors convey their appreciation to Miguel Vergara for his technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Corpino, S.; Caldera, M.; Nichele, F.; Masoero, M.; Viola, N. Thermal design and analysis of a nanosatellite in low earth orbit. Acta Astronaut. 2015, 115, 247–261. [Google Scholar] [CrossRef] [Green Version]
  2. Mehrparvar, A. Cubesat Design Specifications Rev. 13; California Polytechnic State University: San Luis Obispo, CA, USA; Stanford University: Stanford, CA, USA, 2014. [Google Scholar]
  3. Niaki, K.S.; Anvari, A.; Farhani, F. Aluminum and composite materials for satellite structures—A comparison of thermal performance. Mater. Sci. Res. India 2007, 4, 25–34. [Google Scholar] [CrossRef]
  4. Condruz, M.R.; Voicu, L.R.; Puscasu, C.; Vintila, I.S.; Sima, M.; Deaconu, M.; Dragasanu, L. Composite material designs for lightweight space packaging structures. INCAS Bull. 2018, 10, 13–25. [Google Scholar] [CrossRef]
  5. Ampatzoglou, A.; Baltopoulos, A.; Kotzakolios, A.; Kostopoulos, V. Qualification of Composite Structure for Cubesat Picosatellites as a Demonstration for Small Satellite Elements. IJASAR 2014, 1, 1–10. [Google Scholar] [CrossRef]
  6. Shemelya, C.; De la Rosa, A.; Torrado, A.R.; Yu, K.; Domanowski, J.; Bonacuse, P.J.; Martin, R.E.; Juhasz, M.; Hurwitz, F.; Wicker, R.B.; et al. Anisotropy of thermal conductivity in 3D printed polymer matrix composites for space based cube satellites. Addit. Manuf. 2017, 16, 186–196. [Google Scholar] [CrossRef]
  7. Marshall, W.M.; Zemba, M.; Shemelya, C.; Wicker, R.; Espalin, D.; MacDonald, E.; Keif, C.; Kwas, A. Using Additive Manufacturing to Print a CubeSat Propulsion System. In Proceedings of the 51st AIAA/SAE/ASEE Joint Propulsion Conference, Orlando, FL, USA, 27–29 July 2015. [Google Scholar]
  8. Borgeaud, M.; Scheidegger, N.; Noca, M.; Roerhlisberger, G.; Jordan, F.; Choueiri, T.; Steiner, N. SwissCube: The First Entirely-Built Swiss Student Satellite with an Earth Observation Payload. In Small Satellite Missions for Earth Observation; Sandau, R., Roeser, H.P., Valenzuela, A., Eds.; Springer: Berlin, Germany, 2010; pp. 207–213. ISBN 978-3-642-03501-2. [Google Scholar]
  9. Hussain, F.; Hojjati, M.; Okamoto, M.; Gorga, R.E. Review article: Polymer-matrix Nanocomposites, Processing, Manufacturing, and Application: An Overview. J. Compos. Mater. 2006, 40, 1511–1575. [Google Scholar] [CrossRef]
  10. Paul, D.R.; Robeson, L.M. Polymer nanotechnology: Nanocomposites. Polymer 2008, 49, 3187–3204. [Google Scholar] [CrossRef] [Green Version]
  11. Gonzalez-Garcia, P.; Ramirez-Aguilar, R.; Torres, M.; Franco-Urquiza, E.A.; May-Crespo, J.; Camacho, N. Mechanical and thermal behavior dependence on graphite and oxidized graphite content in polyester composites. Polymer 2018, 153, 9–16. [Google Scholar] [CrossRef]
  12. Kang, S.J.; Oh, H.U. On-Orbit Thermal Design and Validation of 1U Standardized CubeSat of STEP Cube Lab. Int. J. Aerosp. Eng. 2016, 2016, 4213189. [Google Scholar] [CrossRef]
  13. Panczak, T.D.; Ring, S.G. RadCAD®: Next Generation Thermal Radiation Analyzer. SAE Int. 1997. [Google Scholar] [CrossRef]
  14. Chandrasekaran, V.; Subramanian, E.R. Transient thermal analysis of a nanosatellite in low earth orbit. In Proceedings of The Eighth International Conference on Engineering Computational Technology; Topping, B.H.V., Ed.; Civil-Comp Press: Stirlingshire, UK, 2012. [Google Scholar]
  15. Yoo, J.; Jin, H.; Seon, J.; Jeong, Y.H.; Glaser, D.; Lee, D.H.; Lin, R.P. Thermal Analysis of TRIO-CINEMA Mission. J. Astron. Space Sci. 2012, 29, 23–31. [Google Scholar] [CrossRef] [Green Version]
  16. Kandelbauer, A. Cyanate Esters. In Handbook of Thermoset Plastics, 3rd ed.; Dodiuk, H., Goodman, S.H., Eds.; William Andrew Publishing: Boston, MA, USA, 2014; pp. 425–457. [Google Scholar]
  17. Toray Advanced Composites. Toray EX-1510 Cyanate Resin, Two Part System Data Sheet. Available online: https://www.toraytac.com/product-explorer/products/ntro/EX-1510 (accessed on 8 August 2019).
  18. FAR 25.853—Appendix F Part II. Fire Test to Aircraft Material—Flammability of Seat Cushions; Federal Aviation Administration, Department of Transportation: Washington, DC, USA, 2018.
  19. Sika-AXSON Corporation. EPOLAM 2500, Epoxy System FAR 25 Approved. Available online: https://advanced-resins.sika.com (accessed on 10 September 2018).
  20. Nakamura, S.; Fuji, T.; Matsukawa, S.; Katagiri, M.; Fukuyama, H. Specific heat, thermal conductivity, and magnetic susceptibility of cyanate-ester resins—An alternative to commonly used epoxy resins. Cryogenics 2018, 95, 76–81. [Google Scholar] [CrossRef]
  21. Toldy, A.; Szlancsik, A.; Szolnoki, B. Reactive flame retardancy of cyanate ester/epoxy resin blends and their carbon fibre reinforced composites. Polym. Degrad. Stab. 2016, 128, 29–38. [Google Scholar] [CrossRef] [Green Version]
  22. Chen, J.; Liu, B.; Yan, L. Transport of thermal energy in epoxy matrix composites reinforced with a hybrid carbon nanofiller. Results Phys. 2019, 14, 102363. [Google Scholar] [CrossRef]
  23. European Commission Directorate-General for Energy. Horizon 2020, EU Energy Trends to 2030. Climate Action DG and Mobility and Transport DG, 2020 Climate and Energy Package; Publications Office of the European Union: Luxembourg, 2010. [Google Scholar]
  24. Torres, M.; Piedra, S.; Ledesma, S.; Perez, R.; Franco, J.A. Thermal and mechanical analysis of a Carbon-Epoxy laminates for CubeSat primary structures. In Proceedings of the 5th International Conference on Mechanics of Composites (MechComp 2019), Lisbon, Portugal, 1–4 July 2019. [Google Scholar]
  25. Patterson, B.A.; Sodano, H.A. Enhanced Interfacial Strength and UV Shielding of Aramid Fiber Composites through ZnO Nanoparticle Sizing. ACS Appl. Mater. Interfaces 2016, 8, 33963–33971. [Google Scholar] [CrossRef] [PubMed]
  26. Salla, J.; Pandey, K.K.; Srinivas, K. Improvement of UV resistance of wood surfaces by using ZnO nanoparticles. Polym. Degrad. Stab. 2012, 97, 592–596. [Google Scholar] [CrossRef]
  27. Camacho, N.; May-Crespo, J.F.; Rojas-Trigos, J.B.; Mondragon-Rodriguez, G.C.; Martinez, K.; Marin, E. Thermal properties and degradation kinetics of epoxy-γ-alumina and epoxy-zinc oxide composites. Polym. Test. 2019. Submitted. [Google Scholar]
  28. Sanchez-San Juan, S.; Gonzalez-Llorente, J.; Hurtado-Velasco, R. Comparison of the Incident Solar Energy and Battery Storage in a 3U CubeSat Satellite for Different Orientation Scenarios. J. Aerosp. Technol. Manag. 2016, 8, 91–102. [Google Scholar] [CrossRef] [Green Version]
  29. HEXCEL Corporation. HexTow AS4 Carbon fibre Datasheet; Hexcel Corporation: Stamford, CT, USA, 2016. [Google Scholar]
  30. Anderson, B.J.; Justus, C.J.; Batts, G.W. Guidelines for the Selection of Near-Earth Thermal Environment Parameters for Spacecraft Design; NASA/TM-2001-211221; NASA: Huntsville, AL, USA, 2001.
  31. Jaques, L. Thermal Design of the OUFTI-1 Nanosatellite. Master’s Thesis, University of Liege, Liege, Belgium, 2009. [Google Scholar]
  32. Gilmore, D. Spacecraft Thermal Control Handbook, Volume 1: Fundamental Technologies; Aerospace Press: Reston, VA, USA, 2002. [Google Scholar]
  33. ANSYS Mechanical User’s Guide. Available online: https://www.sharcnet.ca/Software/Ansys/18.2.2/en-us/help/ai_sinfo/mech_intro.html (accessed on 25 January 2019).
  34. Dasgupta, A.; Agarwal, R.K. Orthotropic thermal conductivity of plain-weave fabric composites using a homogenization technique. J. Compos. Mater. 1992, 26, 2736–2758. [Google Scholar] [CrossRef]
  35. Hubbard, J.A.; Brown, A.L.; Dodd, A.B.; Gomez-Vasquez, S.; Ramirez, C.J. Carbon Fiber Composite Characterization in Adverse Thermal Environments; SAND2011-2833; Sandia National Laboratories: Albuquerque, NM, USA, 2011.
Figure 1. Conceptual design of the CubeSat primary structure, (a) external composite panels (b) internal electronic devices arrangement and (c) panels configuration and Nadir-pointing scenario scheme.
Figure 1. Conceptual design of the CubeSat primary structure, (a) external composite panels (b) internal electronic devices arrangement and (c) panels configuration and Nadir-pointing scenario scheme.
Aerospace 06 00097 g001
Figure 2. Radiative heat fluxes over each CubeSat panel as function of time along one orbit ( β = 57 ).
Figure 2. Radiative heat fluxes over each CubeSat panel as function of time along one orbit ( β = 57 ).
Aerospace 06 00097 g002
Figure 3. Computational mesh and radiative heat fluxes over panel 2 as function of time, black, red and green curves correspond to the direct solar, albedo and Earth infrared radiative heat fluxes over P2 ( β = 57 ).
Figure 3. Computational mesh and radiative heat fluxes over panel 2 as function of time, black, red and green curves correspond to the direct solar, albedo and Earth infrared radiative heat fluxes over P2 ( β = 57 ).
Aerospace 06 00097 g003
Figure 4. Temperature field on the surface of the CubeSat panels at t = 19,620 s ( β = 57 ).
Figure 4. Temperature field on the surface of the CubeSat panels at t = 19,620 s ( β = 57 ).
Aerospace 06 00097 g004
Figure 5. Temperature field on the surface of the CubeSat panels at t = 22,856 s ( β = 57 ).
Figure 5. Temperature field on the surface of the CubeSat panels at t = 22,856 s ( β = 57 ).
Aerospace 06 00097 g005
Figure 6. Average temperature in each composite panel as function of time ( β = 57 ).
Figure 6. Average temperature in each composite panel as function of time ( β = 57 ).
Aerospace 06 00097 g006
Figure 7. Temperature field at the surface of the CubeSat panels and electronic boards, (a) t = 19,620 s, (b) t = 22,856 s ( β = 57 ).
Figure 7. Temperature field at the surface of the CubeSat panels and electronic boards, (a) t = 19,620 s, (b) t = 22,856 s ( β = 57 ).
Aerospace 06 00097 g007
Figure 8. Average temperature in the electronic boards as function of time ( β = 57 ).
Figure 8. Average temperature in the electronic boards as function of time ( β = 57 ).
Aerospace 06 00097 g008
Figure 9. Average temperature in each aluminum panel as function of time ( β = 57 ).
Figure 9. Average temperature in each aluminum panel as function of time ( β = 57 ).
Aerospace 06 00097 g009
Figure 10. Average temperature in each composite panel as function of time ( β = 90 ).
Figure 10. Average temperature in each composite panel as function of time ( β = 90 ).
Aerospace 06 00097 g010
Figure 11. Average temperature in each composite panel as function of time ( β = 90 ).
Figure 11. Average temperature in each composite panel as function of time ( β = 90 ).
Aerospace 06 00097 g011
Figure 12. Average temperature in the electronic boards as function of time ( β = 90 ).
Figure 12. Average temperature in the electronic boards as function of time ( β = 90 ).
Aerospace 06 00097 g012
Figure 13. Average temperature in the electronic boards as function of time ( β = 0 ).
Figure 13. Average temperature in the electronic boards as function of time ( β = 0 ).
Aerospace 06 00097 g013
Table 1. Operating temperature limits for the main on-board systems [12].
Table 1. Operating temperature limits for the main on-board systems [12].
UnitPower Dissipation (W)Operating Temperature Range (°C)
EPS board0.1−4085
Battery pack0.6−206060
GPS board0.27−4085
Communication board0.27−3060
Table 2. Thermal properties of the materials [19,27].
Table 2. Thermal properties of the materials [19,27].
MaterialThermal Conductivity (W/m K)Heat Capacity (J/kg K)Density (kg/m3)
Epolam 2500/ZnO 0.25 wt %0.418901210
Epolam 25000.3517601210
AS4 Carbon Fibre6.8311301790
Al 5052-H321388802680

Share and Cite

MDPI and ACS Style

Piedra, S.; Torres, M.; Ledesma, S. Thermal Numerical Analysis of the Primary Composite Structure of a CubeSat. Aerospace 2019, 6, 97. https://doi.org/10.3390/aerospace6090097

AMA Style

Piedra S, Torres M, Ledesma S. Thermal Numerical Analysis of the Primary Composite Structure of a CubeSat. Aerospace. 2019; 6(9):97. https://doi.org/10.3390/aerospace6090097

Chicago/Turabian Style

Piedra, Saul, Mauricio Torres, and Saul Ledesma. 2019. "Thermal Numerical Analysis of the Primary Composite Structure of a CubeSat" Aerospace 6, no. 9: 97. https://doi.org/10.3390/aerospace6090097

APA Style

Piedra, S., Torres, M., & Ledesma, S. (2019). Thermal Numerical Analysis of the Primary Composite Structure of a CubeSat. Aerospace, 6(9), 97. https://doi.org/10.3390/aerospace6090097

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop