Simulación numérica mediante analogía termoeléctrica de problemas de transmisión de calor en distintas geometrías con propiedades de los materiales dependientes de la temperatura

Abstract

Para la resolución de problemas transmisión de calor, generalmente, una de las aproximaciones más extendidas es suponer que las propiedades de los materiales son constantes con la temperatura. Sin embargo, para modelizar y simular este tipo de problemas con mayor precisión, podría resultar interesante considerar la densidad, la conductividad térmica y el calor específico de los materiales variables con la temperatura. Por ello, en la tesis, se planteó como objetivo general la implementación de problemas de transmisión de calor considerando la dependencia de las propiedades de los materiales con la temperatura. En este sentido, se propusieron ejemplos y ensayos con distintos materiales y distintas geometrías para valorar cuantitativa y cualitativamente la reducción de errores y, con ello, la mayor exactitud de cálculo en la resolución numérica de los procesos de transmisión de calor. La metodología que se siguió contó con varias fases diferenciadas. En primer lugar, se desarrollaron los modelos matemáticos en distintas geometrías (pared plana, cilíndrica y esférica) para la resolución de los problemas de transmisión de calor, considerando en todos los casos modelos con las propiedades de los materiales dependientes de la temperatura. En segundo lugar, se diseñaron los modelos en red equivalentes a los modelos matemáticos desarrollados anteriormente para cada una de las geometrías. De esta forma, se definieron los circuitos equivalentes para cada una de las celdas asociadas a cada geometría. A continuación, se programaron los modelos en red en Matlab y Ngspice y las distintas posibilidades para la representación gráfica de la solución de los problemas. Antes de comenzar las simulaciones de problemas reales, se validó el software mediante un ensayo de laboratorio comparando los resultados obtenidos experimentalmente con los resultados arrojados por el programa. Para ello, se tomaron tres botellas de uso común en el consumo de agua de diferentes materiales, con un doble objetivo: en primer lugar, validar el modelo propuesto con datos experimentales y, en segundo lugar, estudiar la variación de temperatura en el interior de la botella, ya que ésta podría conllevar un aumento en la migración de sustancias al líquido en materiales plásticos. Los resultados obtenidos en el experimento nos llevaron a concluir que el software implementado seguía fielmente la realidad de lo sucedido en el laboratorio, por tanto, se validó de forma favorable. Tras este experimento, se simularon problemas con el software para las diferentes geometrías. Se realizaron 4 simulaciones distintas: tanque cilíndrico con etanol líquido y acero inoxidable AISI 316; tanque esférico con 〖CO〗_2 líquido y acero inoxidable AISI 304; botella cilíndrica con amoniaco y aluminio Al 7075 y puerta de aluminio Al 6061. Para todos estos ejemplos, se comparó el efecto existente en la consideración o no de las propiedades de los materiales dependientes de la temperatura y el error cometido al considerarlas constantes. Como conclusiones de la tesis puede afirmarse, por un lado, que la variabilidad de las propiedades de los materiales: densidad, conductividad térmica y calor específico, debe tenerse en cuenta al simular problemas de transmisión de calor donde se requiera una mayor precisión en los resultados y, por otro lado, que el programa desarrollado en Matlab, permite simular problemas con las geometrías: pared plana, cilíndrica y esférica considerando las propiedades dependientes de la temperatura, siguiendo una ecuación de grado dos como máximo, minimizando, con ello, el error de cálculo respecto al que se cometería considerando las propiedades constantes.

In solving heat transfer problems, one of the most common approaches is to assume that material properties are constant with temperature. However, in order to model and simulate this type of problems more accurately, it could be interesting to consider the density, thermal conductivity and specific heat of materials varying with temperature. Therefore, in the thesis, the general objective was the implementation of heat transfer problems considering the dependence of material properties on temperature. In this sense, examples and tests with different materials and different geometries were proposed in order to quantitatively and qualitatively evaluate the reduction of errors and, thus, the highest calculation accuracy in the numerical resolution of heat transfer processes. The methodology followed consisted of several distinct phases. First, mathematical models were developed for different geometries (flat wall, cylindrical and spherical) to solve the heat transfer problems, considering in all cases models with temperature-dependent material properties. Secondly, the network models equivalent to the mathematical models previously developed for each of the geometries were designed. In this way, the equivalent circuits were defined for each of the cells associated with each geometry. Then, the network models were programmed in Matlab and Ngspice and the different possibilities for the graphical representation of the solution of the problems were programmed. Before starting the simulations of real problems, the software was validated by means of a laboratory test comparing the results obtained experimentally with the results obtained by the program. For this purpose, three bottles of different materials commonly used in water consumption were taken, with a double objective: firstly, to validate the proposed model with experimental data and, secondly, to study the temperature variation inside the bottle, since this could lead to an increase in the migration of substances to the liquid in plastic materials. The results obtained in the experiment led us to conclude that the software implemented faithfully followed the reality of what happened in the laboratory and was therefore validated favourably After this experiment, problems were simulated with the software for the different geometries. Four different simulations were performed: cylindrical tank with liquid ethanol and AISI 316 stainless steel; spherical tank with liquid CO2 and AISI 304 stainless steel; cylindrical bottle with ammonia and Al 7075 aluminum and Al 6061 aluminum door. For all these examples, the existing effect of considering or not considering temperature-dependent material properties and the error made when considering them constant were compared. As conclusions of the thesis it can be stated, on the one hand, that the variability of the properties of the materials: density, thermal conductivity and specific heat, must be taken into account when simulating heat transfer problems where greater precision in the results is required and, on the other hand, that the program developed in Matlab, allows simulating problems with the geometries: flat, cylindrical and spherical wall considering the properties dependent on temperature, following an equation of degree two at most, thus minimizing the calculation error with respect to that which would be committed considering the constant properties.