ANÁLISIS Y SIMULACIÓN DE TRANSFERENCIA DE CALOR EN ÓRBITA DE UN CUBESAT USANDO iOS

R.C. CABRIALES*, L.A. REYES*, C. CHÁVEZ*, D. COBOS*,
PATRICIA C. ZAMBRANO-ROBLEDO*

CIENCIA UANL / AÑO 23, No.103 septiembre-octubre 2020

DOI: https://doi.org/10.29105/cienciauanl23.103-1

RESUMEN

La estructura de un CubeSat permite realizar distintas misiones en órbita terrestre baja, siendo de relevancia sus condiciones de operación en la prevención y control de fenómenos térmicos. Dependiendo de su órbita e inclinación, el satélite estará sujeto a diferentes efectos de transferencia de calor. En este trabajo se presenta el desarrollo de un modelo matemático en la evaluación de las condiciones críticas de transferencia de calor en los estados estable y transitorio de un CubeSat. Asimismo, se desarrolló una aplicación en iOS que simula los resultados obtenidos mediante la herramienta SceneKit para producir gráficos 3D.

Palabras clave: CubeSat, radiación, transferencia de calor, iOS, SceneKit.

ABSTRACT

CubeSat´s structure allows to perform space missions in low Earth orbit. Therefore, it is necessary to consider at which conditions it will be exposed, focusing on the prevention and control of thermal phenomena. Depending on its orbit and inclination, the satellite will be exposed to different heat transfer effects. In this research, temperature gradient simulations were performed using a mathematical model in order to obtain the heat transfer critical conditions steady and transient states. An iOS application was then created to simulate the results using the tool SceneKit to produce 3D graphics.

Keywords: CubeSat, radiation, heat transfer, iOS, SceneKit.

En las últimas décadas, la industria aeroespacial ha tenido un significativo interés en el desarrollo e investigación de satélites cada vez más pequeños para reducir su peso y costo. El estándar CubeSat fue creado por la Stanford University y por la California Polytechnic State University en 1999, éste especifica que una unidad de 1U estándar es un cubo de 10 cm (10x10x10 cm3) con una masa de hasta 1.33 kg (Poghosyan y Golkar, 2017). Un 1U CubeSat podría servir como un satélite independiente o podría combinarse para construir una nave espacial más grande. Por ejemplo, un 3U CubeSat tendrá un factor de forma similar a tres 1U CubeSats combinados. Una de las principales ventajas de esta estandarización es permitir a los fabricantes de vehículos de lanzamiento adoptar un sistema de despliegue común independiente del fabricante CubeSat. Dada la naturaleza exitosa del CubeSat, como las unidades 1U y 3U, se ha creado un estándar para mayores CubeSat (6U, 12U y 27U). Por lo general, los satélites pequeños se clasifican basándose sólo en su masa, sin embargo, en el caso del estándar CubeSat también se considera el volumen. La figura 1 proporciona una clasificación generalmente aceptada para satélites pequeños junto con una comparación con el estándar CubeSat.

La composición relativamente simple de un CubeSat permite realizar cualquier misión asignada en órbita terrestre baja, por esta razón, es necesario considerar las diferentes condiciones a las que estará expuesto, prestando atención a la prevención y control de los fenómenos térmicos. En los CubeSat el riesgo radica principalmente en la operación de los componentes electrónicos, así como en las baterías si son expuestas a temperaturas extremas, sobre todo por radiación solar, por la radiación reflejada de la superficie de la Tierra (conocida como albedo), así como la radiación infrarroja emitida por la Tierra.

El albedo es la fracción de la luz solar reflejada desde cualquier planeta al espacio, aumenta con la latitud, varía con la inclinación de la órbita y depende tanto de la superficie local como de las propiedades atmosféricas (Wen y Mudawar, 2004). El factor de albedo promedio y la intensidad de la radiación solar por planeta se muestran en la tabla I.

La actividad solar varía diariamente, y es la variación entre máximos y mínimos solares la que a menudo dicta el diseño del satélite. Los CubeSats están diseñados para permanecer en órbita por un corto periodo de tiempo y por lo tanto el promedio de flujo solar durante un lapso prolongado es suficiente para el análisis de radiación.

La transferencia de calor por radiación ocurre en forma de ondas electromagnéticas que pueden viajar por el vacío, la ley de Stefan-Boltzmann se utiliza para obtener la razón máxima de la radiación de un cuerpo negro que se puede emitir desde una superficie a una temperatura.

donde la constante de Stefan-Boltzmann es σ = 5.67×10⁻⁸ W/(m² K⁴), A es el superficie de emisión y T la temperatura absoluta.

Existen dos propiedades importantes relativas a la radiación sobre una superficie. La primera es la emisividad ε que es la relación entre la radiación emitida por la superficie real y la emitida por un cuerpo negro. La segunda es la absortividad α, que es la fracción de la energía de radiación incidente absorbida por la superficie. Así, la radiación térmica emitida por una superficie real puede ser obtenida cómo:

La radiación térmica absorbida por una superficie real se obtiene:

Donde a G se le conoce como irradiación; otro parámetro usado es el factor de forma F que representa la razón del flujo neto entre sus superficies debido a su inclinación, en el caso de un CubeSat alrededor de la Tierra, debido a sus diferencias de tamaño, el factor de forma se aproxima a 1. De igual manera, entre el Sol y el satélite, debido a su distancia, también se aproxima a 1 por lo que no será necesario considerarlo en el modelo 1 (Incropera et al., 2007; Duffie y Beckman, 2006).

En el caso de la irradiancia solar, ésta se ha obtenido experimentalmente a partir de satélites patrocinados por agencias espaciales, principalmente la NASA. Siendo uno de los más conocidos el Experimento Solar de Radiación y Clima (SORCE), que se lanzó el 25 de enero de 2003 y ha recibido mediciones precisas de radiación solar; para este trabajo se usaron los datos del instrumento abordo llamado TIM (Total Irradiance Monitor). La nave espacial SORCE fue lanzada en una órbita de 645 km y 40 grados y es operada por el Laboratorio de Física Atmosférica y Espacial en Colorado, EE. UU.

La absorción del Sol al satélite puede ser calculada como el producto de la absortividad solar αꚃ, el área expuesta del satélite e irradiación solar Gꚃ. La absorción del CubeSat por el albedo puede ser calculada con la irradiación solar Gꚃ multiplicada por el factor del albedo fᴀ. La irradiación infrarroja emitida por la superficie de la Tierra puede ser calculada con base en la energía solar incidente Iꚍ sobre el área transversal de la Tierra considerando un radio R:

Esta radiación se extiende sobre la superficie esférica de la Tierra, por lo que la energía emitida por unidad de superficie puede ser calculada (Cisneros, et al., 2017) como:

Los parámetros previamente mencionados se deberán considerar en el desarrollo del modelo térmico. La simulación se realizó a través de Swift, un lenguaje de programación multiparadigma creado por Apple para el desarrollo de aplicaciones para iOS y macOS. Swift es un lenguaje de desarrollo rápido de propósito general, usa el compilador LLVM incluido en Xcode que es el entorno de desarrollo para macOS. En 2015 este lenguaje pasó a código abierto (Apple Inc., 2018). Este lenguaje es derivado de C con conceptos de programación similares a lenguajes de scripting como Ruby o Phyton (Swift, 2016).

La visualización gráfica se elaboró usando SceneKit, un conjunto de librerías y códigos que permite trabajar de acuerdo con una metodología que proporciona los recursos requeridos en una escena 3D, este conjunto maneja automáticamente el código necesario para un renderizado de OpenGL (Open Graphics Library) (Apple Inc., 2018), es decir, la especificación multilenguaje y multiplataforma para escribir aplicaciones que generan gráficos 2D y 3D.

CONTRIBUCIÓN

En este trabajo se utilizaron modelos de transferencia de calor por radiación en estado estable y transitorio para elaborar una aplicación móvil en iOS que simule las condiciones térmicas críticas de la órbita de un CubeSat utilizando herramientas tecnológicas que sean accesibles y fáciles de entender para todo público. Los resultados obtenidos fueron comparados con datos experimentales presentados en la bibliografía y la aplicación ha sido publicada en la appstore como aplicación gratuita con el nombre de CubeSat en la liga: https://itunes.apple.com/us/app/cubesat/id1358767556?l=es&ls=1 &mt=8

MODELO MATEMÁTICO

En condiciones estacionarias sobre el CubeSat se tiene la siguiente relación:

Donde las energías de radiación involucradas son:

La emisión del satélite al espacio:

Donde es la absortividad solar (adimensional),  la absortividad por radiación ir de la Tierra (adimensional),   la emisividad por radiación ir de la Tierra (adimensional),  la irradiancia solar (W/m²), fᴀ el factor de albedo (adimensional), ET la energía ir emitida por la Tierra (W/m²), la temperatura del espacio (K), el área expuesta del satélite hacia el Sol (m²), el área expuesta del satélite por el albedo (m²),  el área expuesta del satélite hacia la Tierra (m²) y el área expuesta del satélite hacia el espacio (m²).

De la ley de Kirchhoff de la radiación térmica se tiene la siguiente igualdad:

Considerando condiciones extremas donde una cara A del CubeSat se encuentra de forma perpendicular al Sol, otra cara se encuentra enfocada hacia la Tierra y las demás se encuentran dirigidas al espacio, podemos escribir:

Dividiendo los dos lados de la igualdad entre A:

Las condiciones en el estado transitorio en el CubeSat pueden ser representadas por:

Donde es la temperatura del satélite en el tiempo,  la temperatura del satélite en estado estable, siendo el límite para el estado transitorio; Cp es el calor especifico del satélite (J/(kg,K) y m la masa del satélite (kg).

Además de las ecuaciones de estado estable se conoce que:

La ecuación (19) en notación simple puede ser extendida a los diferentes materiales que presenta el satélite en la superficie, siendo Cτ el calor específico total considerando distintos materiales, en este trabajo i tiene el valor de dos componentes: aluminio y la celda solar. Además, descomponiendo en fracciones parciales e integrando:

En esta expresión T no se puede evaluar directamente; siendo T el cambio de temperatura con respecto al tiempo t, su valor puede ser encontrado iterativamente debido a que la expresión del lado izquierdo de la igualdad está compuesta por cuatro constantes y el tiempo como única variable, mientras la expresión del lado derecho presenta dos constantes Tο y , donde Tο es la temperatura inicial obtenida de los datos experimentales del CubeSat CP3
(Friedel y McKibbon, 2011) y se obtiene de (15). La variable T de la expresión del lado derecho se puede calcular para un determinado valor de t desarrollando iteraciones en el valor de T hasta encontrar el valor que satisfaga la igualdad. Mediante el desarrollo de un algoritmo se puede variar el valor de t desde un tiempo cero
hasta un tiempo definido, dependiendo del tiempo que completa media órbita para el lado iluminado o lado oscuro según sea el caso. En esta ecuación, representa el valor límite que podría alcanzar T, presentando un valor diferente para el lado iluminado en comparación al lado oscuro debido a la ausencia del término   en este último caso.

RESULTADOS

Para la validación del modelo se utilizaron datos experimentales recopilados del CubeSat CP3 (1U) lanzado en 2007 con una órbita de inclinación de 97.3 grados y un periodo de 99.14 minutos (Friedel y McKibbon, 2011). En el CubeSat CP3, 42% del área está cubierta con paneles solares con una emisividad ε = 0.825, absortividad α = 0.805 y = 1600 J/(kg K), el resto de la superficie puede considerarse 58% de aluminio 7075 con ε=0.045, α=0.25, =960 J/(kg K), además, considerando que el estándar CubeSat recomienda que la estructura de aluminio no debe pesar más de 30% de la masa total, 0.3 kg de aluminio, 0.03 kg de celda solar y una temperatura del espacio de 2.7 K (Poghosyan y Golkar, 2017; Wertz y Larson, 2005). La energía emitida por la Tierra también debe corregirse debido al albedo de (5).

Para un =35%, que puede variar con el ángulo beta, y con base en el modelo matemático, se obtuvieron dos temperaturas del estado estable según la irradiancia solar, una para el lado iluminado (=1355.827 W/m²) y otra para el lado oscuro (=0 W/m²), con estas temperaturas se realizaron los cálculos del estado transitorio con condiciones iniciales de 245 K (Friedel y McKibbon, 2011), obteniendo la variación de temperatura sobre su órbita. Los resultados del algoritmo para la órbita se muestran en la figura 2, la línea negra representa los valores del modelo y la línea azul los valores experimentales del CP3, es importante observar que el traslape presente en los valores mínimos es de esperarse debido a que como temperatura inicial del modelo se está utilizando la misma temperatura inicial del satélite CP3. Se puede resaltar que ésta se ajusta correctamente al comportamiento observado.

La simulación y la visualización que se realizó mediante SceneKit se muestra en la figura 3.

CONCLUSIONES

Se desarrolló una aplicación que proporciona un marco para comprender los efectos térmicos en un CubeSat accesible tanto para investigadores como para estudiantes, la aplicación que se desarrolló proporciona distribuciones de temperatura y capacidades de modelado. Esta herramienta muestra la evolución temporal del perfil de temperatura del CubeSat. La aplicación está diseñada como una herramienta eficiente para el análisis térmico de nanosatélites, la cual complementará a las ya existentes herramientas de investigación.

 

*Universidad Autónoma de Nuevo León.
Contacto: roberto.cabrialesgm@uanl.edu.mx

 

REFERENCIAS

Apple Inc. (2018). https://swift.org
Apple Inc. (2018). https://developer.apple.com/documentation/scenekit
Chang-Da, W., y Mudawar, I. (2004). Emissivity characteristics of roughened aluminum alloy surfaces and assessment of multispectral radiation thermometry (MRT) emissivity models. International Journal of Heat and Mass Transfer. 47:3591-3605.
Cisneros, J.C., Martínez, N.S., Mendoza, U.A., et al. (2017). Analysis of the temperature of a 1U CubeSat due to radiation in space, Advances in Structural Engineering and Mechanics, Seoul, Korea.
Duffie, J.A., y Beckman, W.A. (2006). Solar Engineering of Thermal Processes. Wiley and Sons, Hoboken, New Jersey, third edition.
Friedel, J., y McKibbon, S. (2011). Thermal Analysis of the CubeSat CP3 Satellite, Senior Project, California Polytechnic State University, Aerospace Engineering Department.
Incropera, F., DeWitt, D., Bergman, T., et al. (2007). Fundamentals of Heat and Mass Transfer. 6th Edition. John Wiley & Sons Inc. USA
Lattner, C. (2016). Swift: Opportunities for Language and Compiler Research IBM PL Day 2016, Yorktown Heights, NY
Poghosyan, A., y Golkar, A. (2017). CubeSat evolution: Analyzing CubeSat capabilities for conducting science missions, Progress in Aerospace Sciences. Skolkovo Institute of Science and Technology. 88:59-83.
Wertz, J., y Larson J., W. (2010). Space Mission Analysis and Design. Space Technology Library. California, USA.