Computational simulations of latent heat thermal energy storage systems - with innovative and first-principles based simulation for the underlying unsteady melting (and solidification) processes

This thesis develops an effective modeling and simulation procedure for a specific thermal energy storage system commonly used and recommended for various applications (such as an auxiliary energy storage system for solar heating based Rankine cycle power plant). This thermal energy storage system transfers heat from a hot fluid (termed as heat transfer fluid - HTF) flowing in a tube to the surrounding phase change material (PCM). Through unsteady melting or freezing process, the PCM absorbs or releases thermal energy in the form of latent heat. Both scientific and engineering information is obtained by the proposed first-principle based modeling and simulation procedure. On the scientific side, the approach accurately tracks the moving melt-front (modeled as a sharp liquid-solid interface) and provides all necessary information about the time-varying heat-flow rates, temperature profiles, stored thermal energy, etc. On the engineering side, the proposed approach is unique in its ability to accurately solve – both individually and collectively – all the conjugate unsteady heat transfer problems for each of the components of the thermal storage system. This yields critical system level information on the various time-varying effectiveness and efficiency parameters for the thermal storage system.

Sensible and latent heat fluxes in the atmospheric surface layer over the equatorial Atlantic (Appendix)

During a four weeks anchoring station of R.V. ,,Meteor" on the equator at 30° W longitude, vertical profiles of wind, temperature, and humidity were measured by means of a meteorological buoy carrying a mast of 10 m height. After eliminating periods of instrumental failure, 18 days are available for the investigation of the diurnal variations of the meteorological parameters and 9 days for the investigation of the vertical heat fluxes. The diurnal variations of the above mentioned quantities are caused essentially by two periodic processes: the 24-hourly changing solar energy supply and the 12-hourly oscillation of air pressure, which both originate in the daily rotation of the earth. While the temperature of the water and of the near water layers of the air show a 24 hours period in their diurnal course, the wind speed, as a consequence of the pressure wave, has a 12 hours period, which is also observable in evaporation and, consequently, in the water vapor content of the surface layer. Concerning the temperature, a weak dependence of the daily amplitude on height was determined. Further investigation of the profiles yields relations between the vertical gradients of wind, temperature, and water vapor and the wind speed, the difference between sea and air of temperature and water vapor, respectively, thus giving a contribution to the problem of parameterizing the vertical fluxes. Mean profile coefficients for the encountered stabilities, which were slightly unstable, are presented, and correction terms are given due to the fact that the conditions at the very surface are not sufficiently represented by measuring in a water depth of 20 cm and assuming water vapor saturation. This is especially true for the water vapor content, where the relation between the gradient and the air-sea difference suggests a reduction of relative humidity to appr. 96% at the very surface, if the gradients are high. This effect may result in an overestimation of the water vapor flux, if a ,,bulk"-formula is used. Finally sensible and latent heat fluxes are computed by means of a gradient-formula. The influence of stability on the transfer process is taken into account. As the air-sea temperature differences are small, sensible heat plays no important role in that region, but latent heat shows several interesting features. Within the measuring period of 18 days, a regular variation by a factor of ten is observed. Unperiodic short term variations are superposed by periodic diurnal variations. The mean diurnal course shows a 12-hours period caused by the vertical wind speed gradient superposed by a 24-hours period due to the changing stabilities. Mean values within the measuring period are 276 ly/day for latent heat and 9.41y/day for sensible heat.

Ultra high temperature latent heat energy storage and thermophotovoltaic energy conversion

A conceptual energy storage system design that utilizes ultra high temperature phase change materials is presented. In this system, the energy is stored in the form of latent heat and converted to electricity upon demand by TPV (thermophotovoltaic) cells. Silicon is considered in this study as PCM (phase change material) due to its extremely high latent heat (1800 J/g or 500 Wh/kg), melting point (1410 C), thermal conductivity (~25 W/mK), low cost (less than $2/kg or $4/kWh) and abundance on earth. The proposed system enables an enormous thermal energy storage density of ~1 MWh/m3, which is 10e20 times higher than that of lead-acid batteries, 2e6 times than that of Li-ion batteries and 5e10 times than that of the current state of the art LHTES systems utilized in CSP (concentrated solar power) applications. The discharge efficiency of the system is ultimately determined by the TPV converter, which theoretically can exceed 50%. However, realistic discharge efficiencies utilizing single junction TPV cells are in the range of 20e45%, depending on the semiconductor bandgap and quality, and the photon recycling efficiency. This concept has the potential to achieve output electric energy densities in the range of 200-450 kWhe/m3, which is comparable to the best performing state of the art Lithium-ion batteries.

Economic efficiency of mobile latent heat storages

In a pilot project an optimized mobile latent heat storage based on a system available on the market has been tested at Fraunhofer Institute for Environmental, Safety and Energy Technology. Initially trials were conducted with the aim of optimizing the process of charging and discharging. A specifically constructed test rig at the incineration trials centre at the institute allowed charging and discharging procedures of the mobile latent heat storage with adjustable parameters. In addition an evaluation model was constructed to further optimize the heat exchanger systems. In conclusion the prototype of the mobile latent heat storage was tested in practical operation. The economic and technical feasibility of heat transportation was shown if not utilized waste heat is available. © 2014 The Authors.

Multiple solutions and numerical analysis to the dynamic and stationary models coupling a delayed energy balance model involving latent heat and discontinuous albedo with a deep ocean

We study a climatologically important interaction of two of the main components of the geophysical system by adding an energy balance model for the averaged atmospheric temperature as dynamic boundary condition to a diagnostic ocean model having an additional spatial dimension. In this work, we give deeper insight than previous papers in the literature, mainly with respect to the 1990 pioneering model by Watts and Morantine. We are taking into consideration the latent heat for the two phase ocean as well as a possible delayed term. Non-uniqueness for the initial boundary value problem, uniqueness under a non-degeneracy condition and the existence of multiple stationary solutions are proved here. These multiplicity results suggest that an S-shaped bifurcation diagram should be expected to occur in this class of models generalizing previous energy balance models. The numerical method applied to the model is based on a finite volume scheme with nonlinear weighted essentially non-oscillatory reconstruction and Runge–Kutta total variation diminishing for time integration.

Mathematical simulation parameters for drying of cassava starch pellets

Based on experimental tests, it was obtained the equations for drying, equilibrium moisture content, latent heat of vaporization of water contained in the product and the equation of specific heat of cassava starch pellets, essential parameters for realizing modeling and mathematical simulation of mechanical drying of cassava starch for a new technique proposed, consisting of preformed by pelleting and subsequent artificial drying of starch pellets. Drying tests were conducted in an experimental chamber by varying the air temperature, relative humidity, air velocity and product load. The specific heat of starch was determined by differential scanning calorimetry. The generated equations were validated through regression analysis, finding an appropriate correlation of the data, which indicates that by using these equations, can accurately model and simulate the drying process of cassava starch pellets.

An assessment of the latent and sensible heat flux on the simulated regional climate over Southwestern South Atlantic Ocean

A Regional Climate Model (RegCM3) 10-year (1990-1999) simulation over southwestern South Atlantic Ocean (SAO) is evaluated to assess the mean climatology and the simulation errors of turbulent fluxes over the sea. Moreover, the relationship between these fluxes and the rainfall over some cyclogenetic areas is also analyzed. The RegCM3 results are validated using some reanalyses datasets (ERA40, R2, GPCP and WHOI). The summer and winter spatial patterns of latent and sensible heat fluxes simulated by the RegCM3 are in agreement with the reanalyses (WHOI, R2 and ERA40). They show large latent heat fluxes exchange in the subtropical SAO and at higher latitudes in the warm waters of Brazil Current. In particular, the magnitude of RegCM3 latent heat fluxes is similar to the WHOI, which is probably related to two factors: (a) small specific humidity bias, and (b) the RegCM3 flux algorithm. In contrast, the RegCM3 presents large overestimation of sensible heat flux, though it simulates well their spatial pattern. This simulation error is associated with the RegCM3 underestimation of the 2-m air temperature. In southwestern SAO, in three known cyclogenetic areas, the reanalyses and the RegCM3 show the existence of different physical mechanisms that control the annual cycles of latent/sensible heating and rainfall. It is shown that over the eastern coast of Uruguay (35A degrees-43A degrees S) and the southeastern coast of Argentina (44A degrees-52A degrees S) the sea-air moisture and heat exchange play an important role to control the annual cycle of precipitation. This does not happen on the south/southeastern coast of Brazil.

Influência da temperatura do ar de secagem no calor latente de vaporização de água em feijão macassar (Vigna unguiculata (L.) Walp.), variedade sempre-verde

In order to determine the energy needed to artificially dry an agricultural product the latent heat of vaporization of moisture in the product, H, must be known. Generally, the expressions for H reported in the literature are of the form H = h(T)f(M), where h(T) is the latent heat of vaporization of free water, and f(M) is a function of the equilibrium moisture content, M, which is a simplification. In this article, a more general expression for the latent heat of vaporization, namely H = g(M,T), is used to determine H for cowpea, always-green variety. For this purpose, a computer program was developed which automatically fits about 500 functions, with one or two independent variables, imbedded in its library to experimental data. The program uses nonlinear regression, and classifies the best functions according to the least reduced chi-squared. A set of executed statistical tests shows that the generalized expression for H used in this work produces better results of H for cowpea than other equations found in literature.

Modeling of heat transfer in fluidized-bed coating of cylinders

A numerical model of heat transfer in fluidized-bed coating of solid cylinders is presented. By defining suitable dimensionless parameters, the governing equations and its associated initial and boundary conditions are discretized using the method of orthogonal collocation and the resulting ordinary differential equations simultaneously solved for the dimensionless coating thickness and wall temperatures. Parametric Studies showed that the dimensionless coating thickness and wall temperature depend on the relative heat capacities of the polymer powder and object, the latent heat of fusion and the size of the cylinder. Model predictions for the coating thickness and wall temperature compare reasonably well with numerical predictions and experimental coating data in the literature and with our own coating experiments using copper cylinders immersed in nylon-11 and polyethylene powders. (C) 2001 Elsevier Science Ltd. All rights reserved.

Design and construction of precision heat fluxmeters

The possibility of local elastic instabilities is considered in a first¿order structural phase transition, typically a thermoelastic martensitic transformation, with associated interfacial and volumic strain energy. They appear, for instance, as the result of shape change accommodation by simultaneous growth of different crystallographic variants. The treatment is phenomenological and deals with growth in both thermoelastic equilibrium and in nonequilibrium conditions produced by the elastic instability. Scaling of the transformed fraction curves against temperature is predicted only in the case of purely thermoelastic growth. The role of the transformation latent heat on the relaxation kinetics is also considered, and it is shown that it tends to increase the characteristic relaxation times as adiabatic conditions are approached, by keeping the system closer to a constant temperature. The analysis also reveals that the energy dissipated in the relaxation process has a double origin: release of elastic energy Wi and entropy production Si. The latter is shown to depend on both temperature rate and thermal conduction in the system.

Giant heat dissipation at the low temperature reversible-irreversible transition in Gd5Ge4

The heat exchanged at the low-temperature first-order magnetostructural transition is directly measured in Gd5Ge4 . Results show that the origin and the temperature dependence of the heat exchanged varies with the reversible/irreversible character of the first-order transition. In the reversible regime, the heat exchanged by the sample is mostly due to the latent heat at the transition and decreases with decreasing temperature, while in the irreversible regime, the heat is irreversibly dissipated and increases strongly with decreasing temperature, reaching a value of 237 J/kg at 4 K.

Moisture adsorption characteristics of ginger slices

The moisture adsorption characteristics of dried ginger slices was studied to determine the effect of storage conditions on moisture adsorption for the purpose of shelf life prediction, selection of appropriate packaging materials, evaluate the goodness-of-fit of sorption models, and determine the thermodynamics of moisture adsorption for application in drying. There was a highly significant effect (p < 0.05) of water activity (a w), temperature, and pre-treatment on the equilibrium moisture content (EMC) of the dried ginger slices. At constant a w, the EMC decreased as temperature increased. The EMC of all samples increased as the a w increased at constant temperature. The sorbed moisture of the unpeeled ginger slices was higher than the peeled while those of unblanched samples were higher than the blanched. Henderson equation allows more accurate predictions about the isotherms with the lowest %RMS, and therefore, it describes best the adsorption data followed by GAB, Oswin, and Halsey models in that order. The monolayer moisture generally decreased with temperature for all samples. The isosteric heat decreased with moisture content approaching the asymptotic value or the latent heat of vaporization of pure water (∆Hst = 0) while the entropy of sorption was observed to increase with moisture content.

Giant heat dissipation at the low temperature reversible-irreversible transition in Gd5Ge4

The heat exchanged at the low-temperature first-order magnetostructural transition is directly measured in Gd5Ge4 . Results show that the origin and the temperature dependence of the heat exchanged varies with the reversible/irreversible character of the first-order transition. In the reversible regime, the heat exchanged by the sample is mostly due to the latent heat at the transition and decreases with decreasing temperature, while in the irreversible regime, the heat is irreversibly dissipated and increases strongly with decreasing temperature, reaching a value of 237 J/kg at 4 K.

Study of the Oceanic Heat Budget Components over the Arabian Sea during the Formation and Evolution of Super Cyclone, Gonu

Oceans play a vital role in the global climate system. They absorb the incoming solar energy and redistribute the energy through horizontal and vertical transports. In this context it is important to investigate the variation of heat budget components during the formation of a low-pressure system. In 2007, the monsoon onset was on 28th May. A well- marked low-pressure area was formed in the eastern Arabian Sea after the onset and it further developed into a cyclone. We have analysed the heat budget components during different stages of the cyclone. The data used for the computation of heat budget components is Objectively Analyzed air-sea flux data obtained from WHOI (Woods Hole Oceanographic Institution) project. Its horizontal resolution is 1° × 1°. Over the low-pressure area, the latent heat flux was 180 Wm−2. It increased to a maximum value of 210 Wm−2 on 1st June 2007, on which the system was intensified into a cyclone (Gonu) with latent heat flux values ranging from 200 to 250 Wm−2. It sharply decreased after the passage of cyclone. The high value of latent heat flux is attributed to the latent heat release due to the cyclone by the formation of clouds. Long wave radiation flux is decreased sharply from 100 Wm−2 to 30 Wm−2 when the low-pressure system intensified into a cyclone. The decrease in long wave radiation flux is due to the presence of clouds. Net heat flux also decreases sharply to −200 Wm−2 on 1st June 2007. After the passage, the flux value increased to normal value (150 Wm−2) within one day. A sharp increase in the sensible heat flux value (20 Wm−2) is observed on 1st June 2007 and it decreased there- after. Short wave radiation flux decreased from 300 Wm−2 to 90 Wm−2 during the intensification on 1st June 2007. Over this region, short wave radiation flux sharply increased to higher value soon after the passage of the cyclone.