Convecção mista em cavidades com fontes de calor aletadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Sugars Electrooxidation at Glassy Carbon Electrode Decorate with Multi-Walled Carbon Nanotubes with Nickel Oxy-Hydroxide

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Influence of nanoparticle size and gap size on nucleate boiling using HFE7100

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Development and application of a screen-printed electrode modified with MWCNT for the electrocatalytic detection of thiram

Interest in the electronic properties of carbon nanotubes has increased in recent years. These materials can be used in the development of electrochemical sensors for the measurement and monitoring of analytes of environmental interest, such as pharmaceuticals, dyes, and pesticides. This work describes the use of homemade screen-printed electrodes modified with multi-walled carbon nanotubes (MWCNT) for the electrochemical detection of the fungicide thiram. The electrochemical characteristics of the proposed system were evaluated using cyclic voltammetry, with investigation of the electrochemical behavior of the sensor in the presence of the analyte, and estimation of electrochemical parameters including the diffusion coefficient, electron transfer coefficient (α), and number of electrons transferred in the catalytic electro-oxidation. The sensor response was optimized using amperometry. The best sensor performance was obtained in 0.1 mol L-1 phosphate buffer solution at pH 8.0, where a detection limit of 7.9 x 10-6 mol L-1 was achieved. Finally, in order to improve the sensitivity of the sensor, square wave voltammetry (SWV) was used for thiram quantification, instead of amperometry. Using SWV, a response range for thiram from 9.9 x 10-6 to 9.1 x 10-5 mol L-1 was obtained, with a sensitivity of 30948 µA mol L-1, and limits of detection and quantification of 1.6 x 10-6 and 5.4 x 10-6 mol L-1, respectively. The applicability of this efficient new alternative methodology for thiram detection was demonstrated using analyses of enriched soil samples.

Modelo distribuído aplicado à análise de evaporadores do tipo tubo aletado

Pós-graduação em Engenharia Mecânica - FEIS

Flow Boiling Characteristics for R1234ze(E) in 1.0 and 2.2 mm Circular Channels

Experimental flow boiling heat transfer results are presented for horizontal 1.0 and 2.2 mm I. D. (internal diameter) stainless steel tubes for tests with R1234ze(E), a new refrigerant developed as a substitute for R134a with a much lower global warming potential (GWP). The experiments were performed for these two tube diameters in order to investigate a possible transition between macro and microscale flow boiling behavior. The experimental campaign includes mass velocities ranging from 50 to 1500 kg/m(2) s, heat fluxes from 10 to 300 kW/m(2), exit saturation temperatures of 25, 31 and 35 degrees C, vapor qualities from 0.05 to 0.99 and heated lengths of 180 mm and 361 mm. Flow pattern characterization was performed using high speed videos. Heat transfer coefficient, critical heat flux and flow pattern data were obtained. R1234ze(E) demonstrated similar thermal performance to R134a data when running at similar conditions. [DOI: 10.1115/1.4004933]

CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

An investigation into complex inorganic materials with Mössbauer spectroscopy

Mössbauer Spektroskopie ist ein unverzichtbares Instrument für die Bestimmung von Oxidationszuständen und für die Analyse von lokalen Ordnungsphänomenen von Mössbauer aktiven Atomen. Weil es sich um eine lokale Methode handelt können sowohl kristalline als auch amorphe Materialien untersucht werden. Die Kombination von lokaler Prüfung mit Mössbauer Spektroskopie und globaler Untersuchung z.B. mit Röntgendiffraktometrie ermöglicht die Studie von Ordnungseffekten von statistisch besetzten Positionen in einer geordneten Matrix. Das wurde hier eingesetzt um die lokale Umgebung in zwei Serien von Heuslerverbindungen, Co2-xFe1+xSi and Co2Mn1-xFexAl zu untersuchen. Für die Co2Mn1-xFexAl Serie wurde eine L21 geordnete Phase in einer insgesamt B2 geordneten Probe detektiert. Ein Wechsel von der AlCu2Mn zu der CuHg2Ti Struktur wurde für die Co2-xFe1+xSi Proben gefunden. Die Transformation von einem Glas zu einem keramischen Material wurde mit 119Sn Mössbauer Spektroskopie untersucht. Die höhere Ordnung in der Keramik wurde von einer kleiner werdenden Mössbauerlinienbreite begleitet. Demzufolge geben die Modifikationen der Sn Umgebungen klar die Transformation des gesamten Materials wieder. Ist die lokale Umgebung von unregelmäßig auftretenden Atomen in einer amorphen Matrix von Interesse, sind lokal prüfende Methoden die zuverlässigsten Methoden die zur Verfügung stehen. In dieser Arbeit wurde 119Sn Mössbauer Spektroskopie eingesetzt um die Oxidationszustände, die lokalen Umgebungen und relativen Intensitäten von Zinn Atomen in einer Silikatmatrix zu bestimmen. Modifikationen dieser Parameter als Funktion von Prozess bestimmenden Parametern wie der Sauerstoffpartialdruck, die Temperatur, die Behandlungsdauer und der Abkühlprozess genauso wie der SnO2 Gehalt sind von Interesse, weil durch Reduktions- und Diffusionsprozesse Änderungen des Koordinations- und des Oxidationszustands der Zinnatome auftreten. Da diese Änderungen in der Glasmatrix verursachen, die das fertige Produkt im industriellen Fertigungsprozess ruinieren können sind diese feinen Veränderungen sehr wichtig. Wenigstens zwei Mössbauerlinien korrespondierend mit zwei verschiedenen Umgebungen für Sn2+ und Sn4+ sind für eine Analyse mit ausreichender Qualität notwendig. Durch Vergleich von den bestimmten Hyperfein Parametern mit den Parametern von Modelsubstanzen werden lokale Umgebungen der Zinnatome entworfen. Für Sn2+ werden zwei auf einer trigonalen Pyramide basierende Umgebungen mit variierender Anzahl von bindenden und nicht-bindenden Sauerstoffatomen formuliert. Für Sn4+ wurde eine tetraedrische und eine oktaedrische Umgebung postuliert. Die relativen Intensitäten der vier Mössbauerlinien wurden um ein Diffusions- und Reaktionsmodell zu entwickeln und um einen Satz von Diffusions- und Transferkoeffizienten zu bestimmen eingesetzt. Die bestimmten Diffusionskoeffizienten stimmen mit den Literaturdaten überein. Der Massentransferkoeffizient ist kleiner als der bestimmte Wert, aber immer noch in der gleichen Größenordnung. Im Gegensatz zu den Erwartungen ist der präsentierte Diffusionskoeffizient für Sn4+ bestimmt als der von Sn2+. Das wiederum kann durch Berücksichtigung von Elektronhoppingprozessen erklärt werden.

Computational fluids domain reduction to a simplified fluid network

The primary goal of this project is to demonstrate the practical use of data mining algorithms to cluster a solved steady-state computational fluids simulation (CFD) flow domain into a simplified lumped-parameter network. A commercial-quality code, “cfdMine” was created using a volume-weighted k-means clustering that that can accomplish the clustering of a 20 million cell CFD domain on a single CPU in several hours or less. Additionally agglomeration and k-means Mahalanobis were added as optional post-processing steps to further enhance the separation of the clusters. The resultant nodal network is considered a reduced-order model and can be solved transiently at a very minimal computational cost. The reduced order network is then instantiated in the commercial thermal solver MuSES to perform transient conjugate heat transfer using convection predicted using a lumped network (based on steady-state CFD). When inserting the lumped nodal network into a MuSES model, the potential for developing a “localized heat transfer coefficient” is shown to be an improvement over existing techniques. Also, it was found that the use of the clustering created a new flow visualization technique. Finally, fixing clusters near equipment newly demonstrates a capability to track temperatures near specific objects (such as equipment in vehicles).

A Study of Radon in Air and Water in Maine Schools

The transfer coefficient of radon from water to air was investigated in schools. Kitchens, bathrooms and locker rooms were studied for seven schools in Maine. Simulations were done in water-use rooms where radon in air detectors were in place. Quantities measured were radon in water (270-24500 F) and air (0-80 q), volume of water used, emissivities (0.01-0.99) and ventilation rates (0.012-0.066A). Variation throughout the room of the radon concentration was found. Values calculated for the transfer coefficient for kitchens and baths were ranged from 9.6 x to 2.0 x The transfer coefficient was calculated using these parameters and was also measured using concentrations of radon in water and air. This provides a means by which radon in air can be estimated using the transfer coefficient and the concentration in the water in other schools and it can be used to estimate the dose caused by radon released from water use. This project was partially funded by the United States Environmental Protection Agency (grant #X828l2 101-0) and by the State of Maine (grant #10A500178). These are the first measurements of this type to be done in schools in the United States.

Thermal food processing optimization: Algorithms and software

The algorithms and graphic user interface software package ?OPT-PROx? are developed to meet food engineering needs related to canned food thermal processing simulation and optimization. The adaptive random search algorithm and its modification coupled with penalty function?s approach, and the finite difference methods with cubic spline approximation are utilized by ?OPT-PROx? package (http://tomakechoice. com/optprox/index.html). The diversity of thermal food processing optimization problems with different objectives and required constraints are solvable by developed software. The geometries supported by the ?OPT-PROx? are the following: (1) cylinder, (2) rectangle, (3) sphere. The mean square error minimization principle is utilized in order to estimate the heat transfer coefficient of food to be heated under optimal condition. The developed user friendly dialogue and used numerical procedures makes the ?OPT-PROx? software useful to food scientists in research and education, as well as to engineers involved in optimization of thermal food processing.

Laminar counterflow parallel-plate heat exchangers: Exact and approximate solutions

Multilayered, counterflow, parallel-plate heat exchangers are analyzed numerically and theoretically. The analysis, carried out for constant property fluids, considers a hydrodynamically developed laminar flow and neglects longitudinal conduction both in the fluid and in the plates. The solution for the temperature field involves eigenfunction expansions that can be solved in terms of Whittaker functions using standard symbolic algebra packages, leading to analytical expressions that provide the eigenvalues numerically. It is seen that the approximate solution obtained by retaining the first two modes in the eigenfunction expansion provides an accurate representation for the temperature away from the entrance regions, specially for long heat exchangers, thereby enabling simplified expressions for the wall and bulk temperatures, local heat-transfer rate, overall heat-transfer coefficient, and outlet bulk temperatures. The agreement between the numerical and theoretical results suggests the possibility of using the analytical solutions presented herein as benchmark problems for computational heat-transfer codes.

High power beam dump project for the accelerator prototype LIPAc: Cooling design and analysis

En el campo de la fusión nuclear y desarrollándose en paralelo a ITER (International Thermonuclear Experimental Reactor), el proyecto IFMIF (International Fusion Material Irradiation Facility) se enmarca dentro de las actividades complementarias encaminadas a solucionar las barreras tecnológicas que aún plantea la fusión. En concreto IFMIF es una instalación de irradiación cuya misión es caracterizar materiales resistentes a condiciones extremas como las esperadas en los futuros reactores de fusión como DEMO (DEMOnstration power plant). Consiste de dos aceleradores de deuterones que proporcionan un haz de 125 mA y 40 MeV cada uno, que al colisionar con un blanco de litio producen un flujo neutrónico intenso (1017 neutrones/s) con un espectro similar al de los neutrones de fusión [1], [2]. Dicho flujo neutrónico es empleado para irradiar los diferentes materiales candidatos a ser empleados en reactores de fusión, y las muestras son posteriormente examinadas en la llamada instalación de post-irradiación. Como primer paso en tan ambicioso proyecto, una fase de validación y diseño llamada IFMIFEVEDA (Engineering Validation and Engineering Design Activities) se encuentra actualmente en desarrollo. Una de las actividades contempladas en esta fase es la construcción y operación de una acelarador prototipo llamado LIPAc (Linear IFMIF Prototype Accelerator). Se trata de un acelerador de deuterones de alta intensidad idéntico a la parte de baja energía de los aceleradores de IFMIF. Los componentes del LIPAc, que será instalado en Japón, son suministrados por diferentes países europeos. El acelerador proporcionará un haz continuo de deuterones de 9 MeV con una potencia de 1.125 MW que tras ser caracterizado con diversos instrumentos deberá pararse de forma segura. Para ello se requiere un sistema denominado bloque de parada (Beam Dump en inglés) que absorba la energía del haz y la transfiera a un sumidero de calor. España tiene el compromiso de suministrar este componente y CIEMAT (Centro de Investigaciones Energéticas Medioambientales y Tecnológicas) es responsable de dicha tarea. La pieza central del bloque de parada, donde se para el haz de iones, es un cono de cobre con un ángulo de 3.5o, 2.5 m de longitud y 5 mm de espesor. Dicha pieza está refrigerada por agua que fluye en su superficie externa por el canal que se forma entre el cono de cobre y otra pieza concéntrica con éste. Este es el marco en que se desarrolla la presente tesis, cuyo objeto es el diseño del sistema de refrigeración del bloque de parada del LIPAc. El diseño se ha realizado utilizando un modelo simplificado unidimensional. Se han obtenido los parámetros del agua (presión, caudal, pérdida de carga) y la geometría requerida en el canal de refrigeración (anchura, rugosidad) para garantizar la correcta refrigeración del bloque de parada. Se ha comprobado que el diseño permite variaciones del haz respecto a la situación nominal siendo el flujo crítico calorífico al menos 2 veces superior al nominal. Se han realizado asimismo simulaciones fluidodinámicas 3D con ANSYS-CFX en aquellas zonas del canal de refrigeración que lo requieren. El bloque de parada se activará como consecuencia de la interacción del haz de partículas lo que impide cualquier cambio o reparación una vez comenzada la operación del acelerador. Por ello el diseño ha de ser muy robusto y todas las hipótesis utilizadas en la realización de éste deben ser cuidadosamente comprobadas. Gran parte del esfuerzo de la tesis se centra en la estimación del coeficiente de transferencia de calor que es determinante en los resultados obtenidos, y que se emplea además como condición de contorno en los cálculos mecánicos. Para ello por un lado se han buscado correlaciones cuyo rango de aplicabilidad sea adecuado para las condiciones del bloque de parada (canal anular, diferencias de temperatura agua-pared de decenas de grados). En un segundo paso se han comparado los coeficientes de película obtenidos a partir de la correlación seleccionada (Petukhov-Gnielinski) con los que se deducen de simulaciones fluidodinámicas, obteniendo resultados satisfactorios. Por último se ha realizado una validación experimental utilizando un prototipo y un circuito hidráulico que proporciona un flujo de agua con los parámetros requeridos en el bloque de parada. Tras varios intentos y mejoras en el experimento se han obtenido los coeficientes de película para distintos caudales y potencias de calentamiento. Teniendo en cuenta la incertidumbre de las medidas, los valores experimentales concuerdan razonablemente bien (en el rango de 15%) con los deducidos de las correlaciones. Por motivos radiológicos es necesario controlar la calidad del agua de refrigeración y minimizar la corrosión del cobre. Tras un estudio bibliográfico se identificaron los parámetros del agua más adecuados (conductividad, pH y concentración de oxígeno disuelto). Como parte de la tesis se ha realizado asimismo un estudio de la corrosión del circuito de refrigeración del bloque de parada con el doble fin de determinar si puede poner en riesgo la integridad del componente, y de obtener una estimación de la velocidad de corrosión para dimensionar el sistema de purificación del agua. Se ha utilizado el código TRACT (TRansport and ACTivation code) adaptándalo al caso del bloque de parada, para lo cual se trabajó con el responsable (Panos Karditsas) del código en Culham (UKAEA). Los resultados confirman que la corrosión del cobre en las condiciones seleccionadas no supone un problema. La Tesis se encuentra estructurada de la siguiente manera: En el primer capítulo se realiza una introducción de los proyectos IFMIF y LIPAc dentro de los cuales se enmarca esta Tesis. Además se describe el bloque de parada, siendo el diseño del sistema de rerigeración de éste el principal objetivo de la Tesis. En el segundo y tercer capítulo se realiza un resumen de la base teórica así como de las diferentes herramientas empleadas en el diseño del sistema de refrigeración. El capítulo cuarto presenta los resultados del relativos al sistema de refrigeración. Tanto los obtenidos del estudio unidimensional, como los obtenidos de las simulaciones fluidodinámicas 3D mediante el empleo del código ANSYS-CFX. En el quinto capítulo se presentan los resultados referentes al análisis de corrosión del circuito de refrigeración del bloque de parada. El capítulo seis se centra en la descripción del montaje experimental para la obtención de los valores de pérdida de carga y coeficiente de transferencia del calor. Asimismo se presentan los resultados obtenidos en dichos experimentos. Finalmente encontramos un capítulo de apéndices en el que se describen una serie de experimentos llevados a cabo como pasos intermedios en la obtención del resultado experimental del coeficiente de película. También se presenta el código informático empleado para el análisis unidimensional del sistema de refrigeración del bloque de parada llamado CHICA (Cooling and Heating Interaction and Corrosion Analysis). ABSTRACT In the nuclear fusion field running in parallel to ITER (International Thermonuclear Experimental Reactor) as one of the complementary activities headed towards solving the technological barriers, IFMIF (International Fusion Material Irradiation Facility) project aims to provide an irradiation facility to qualify advanced materials resistant to extreme conditions like the ones expected in future fusion reactors like DEMO (DEMOnstration Power Plant). IFMIF consists of two constant wave deuteron accelerators delivering a 125 mA and 40 MeV beam each that will collide on a lithium target producing an intense neutron fluence (1017 neutrons/s) with a similar spectra to that of fusion neutrons [1], [2]. This neutron flux is employed to irradiate the different material candidates to be employed in the future fusion reactors, and the samples examined after irradiation at the so called post-irradiative facilities. As a first step in such an ambitious project, an engineering validation and engineering design activity phase called IFMIF-EVEDA (Engineering Validation and Engineering Design Activities) is presently going on. One of the activities consists on the construction and operation of an accelerator prototype named LIPAc (Linear IFMIF Prototype Accelerator). It is a high intensity deuteron accelerator identical to the low energy part of the IFMIF accelerators. The LIPAc components, which will be installed in Japan, are delivered by different european countries. The accelerator supplies a 9 MeV constant wave beam of deuterons with a power of 1.125 MW, which after being characterized by different instruments has to be stopped safely. For such task a beam dump to absorb the beam energy and take it to a heat sink is needed. Spain has the compromise of delivering such device and CIEMAT (Centro de Investigaciones Energéticas Medioambientales y Tecnológicas) is responsible for such task. The central piece of the beam dump, where the ion beam is stopped, is a copper cone with an angle of 3.5o, 2.5 m long and 5 mm width. This part is cooled by water flowing on its external surface through the channel formed between the copper cone and a concentric piece with the latter. The thesis is developed in this realm, and its objective is designing the LIPAc beam dump cooling system. The design has been performed employing a simplified one dimensional model. The water parameters (pressure, flow, pressure loss) and the required annular channel geometry (width, rugoisty) have been obtained guaranteeing the correct cooling of the beam dump. It has been checked that the cooling design allows variations of the the beam with respect to the nominal position, being the CHF (Critical Heat Flux) at least twice times higher than the nominal deposited heat flux. 3D fluid dynamic simulations employing ANSYS-CFX code in the beam dump cooling channel sections which require a more thorough study have also been performed. The beam dump will activateasaconsequenceofthe deuteron beam interaction, making impossible any change or maintenance task once the accelerator operation has started. Hence the design has to be very robust and all the hypotheses employed in the design mustbecarefully checked. Most of the work in the thesis is concentrated in estimating the heat transfer coefficient which is decisive in the obtained results, and is also employed as boundary condition in the mechanical analysis. For such task, correlations which applicability range is the adequate for the beam dump conditions (annular channel, water-surface temperature differences of tens of degrees) have been compiled. In a second step the heat transfer coefficients obtained from the selected correlation (Petukhov- Gnielinski) have been compared with the ones deduced from the 3D fluid dynamic simulations, obtaining satisfactory results. Finally an experimental validation has been performed employing a prototype and a hydraulic circuit that supplies a flow with the requested parameters in the beam dump. After several tries and improvements in the experiment, the heat transfer coefficients for different flows and heating powers have been obtained. Considering the uncertainty in the measurements the experimental values agree reasonably well (in the order of 15%) with the ones obtained from the correlations. Due to radiological reasons the quality of the cooling water must be controlled, hence minimizing the copper corrosion. After performing a bibligraphic study the most adequate water parameters were identified (conductivity, pH and dissolved oxygen concentration). As part of this thesis a corrosion study of the beam dump cooling circuit has been performed with the double aim of determining if corrosion can pose a risk for the copper beam dump , and obtaining an estimation of the corrosion velocitytodimension the water purification system. TRACT code(TRansport and ACTivation) has been employed for such study adapting the code for the beam dump case. For such study a collaboration with the code responsible (Panos Karditsas) at Culham (UKAEA) was established. The work developed in this thesis has supposed the publication of three articles in JCR journals (”Journal of Nuclear Materials” y ”Fusion Engineering and Design”), as well as presentations in more than four conferences and relevant meetings.

Global air-sea flux of CO2: An estimate based on measurements of sea–air pCO2 difference

Approximately 250,000 measurements made for the pCO2 difference between surface water and the marine atmosphere, ΔpCO2, have been assembled for the global oceans. Observations made in the equatorial Pacific during El Nino events have been excluded from the data set. These observations are mapped on the global 4° × 5° grid for a single virtual calendar year (chosen arbitrarily to be 1990) representing a non-El Nino year. Monthly global distributions of ΔpCO2 have been constructed using an interpolation method based on a lateral advection–diffusion transport equation. The net flux of CO2 across the sea surface has been computed using ΔpCO2 distributions and CO2 gas transfer coefficients across sea surface. The annual net uptake flux of CO2 by the global oceans thus estimated ranges from 0.60 to 1.34 Gt-C⋅yr−1 depending on different formulations used for wind speed dependence on the gas transfer coefficient. These estimates are subject to an error of up to 75% resulting from the numerical interpolation method used to estimate the distribution of ΔpCO2 over the global oceans. Temperate and polar oceans of the both hemispheres are the major sinks for atmospheric CO2, whereas the equatorial oceans are the major sources for CO2. The Atlantic Ocean is the most important CO2 sink, providing about 60% of the global ocean uptake, while the Pacific Ocean is neutral because of its equatorial source flux being balanced by the sink flux of the temperate oceans. The Indian and Southern Oceans take up about 20% each.

Mathematical modeling of drying process of unripe banana slices

Unripe banana flour (UBF) production employs bananas not submitted to maturation process, is an interesting alternative to minimize the fruit loss reduction related to inappropriate handling or fast ripening. The UBF is considered as a functional ingredient improving glycemic and plasma insulin levels in blood, have also shown efficacy on the control of satiety, insulin resistance. The aim of this work was to study the drying process of unripe banana slabs (Musa cavendishii, Nanicão) developing a transient drying model through mathematical modeling with simultaneous moisture and heat transfer. The raw material characterization was performed and afterwards the drying process was conducted at 40 ºC, 50 ºC e 60 ºC, the product temperature was recorded using thermocouples, the air velocity inside the chamber was 4 m·s-1. With the experimental data was possible to validate the diffusion model based on the Fick\'s second law and Fourier. For this purpose, the sorption isotherms were measured and fitted to the GAB model estimating the equilibrium moisture content (Xe), 1.76 [g H2O/100g d.b.] at 60 ºC and 10 % of relative humidity (RH), the thermophysical properties (k, Cp, ?) were also measured to be used in the model. Five cases were contemplated: i) Constant thermophysical properties; ii) Variable properties; iii) Mass (hm), heat transfer (h) coefficient and effective diffusivity (De) estimation 134 W·m-2·K-1, 4.91x10-5 m-2·s-1 and 3.278?10-10 m·s-2 at 60 ºC, respectively; iv) Variable De, it presented a third order polynomial behavior as function of moisture content; v) The shrinkage had an effect on the mathematical model, especially in the 3 first hours of process, the thickness experienced a contraction of about (30.34 ± 1.29) % out of the initial thickness, finding two decreasing drying rate periods (DDR I and DDR II), 3.28x10-10 m·s-2 and 1.77x10-10 m·s-2, respectively. COMSOL Multiphysics simulations were possible to perform through the heat and mass transfer coefficient estimated by the mathematical modeling.