Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.

Convective Cold Pools over the Atlas Mountains and Their Influence on the Saharan Heat Low

The West African Monsoon (WAM) and its representation in numerical models are strongly influenced by the Saharan Heat Low (SHL), a low-pressure system driven by radiative heating over the central Sahara and ventilated by the cold and moist inflow from adjacent oceans. It has recently been shown that a significant part of the southerly moisture flux into the SHL originates from convective cold pools over the Sahel. These density currents driven by evaporation of rain are largely absent in models with parameterized convection. This crucial issue has been hypothesized to contribute to the inability of many climate models to reproduce the variability of the WAM. Here, the role of convective cold pools approaching the SHL from the Atlas Mountains, which are a strong orographic trigger for deep convection in Northwest Africa, is analyzed. Knowledge about the frequency of these events, as well as their impact on large-scale dynamics, is required to understand their contribution to the variability of the SHL and to known model uncertainties. The first aspect is addressed through the development of an objective and automated method for the generation of multi-year climatologies not available before. The algorithm combines freely available standard surface observations with satellite microwave data. Representativeness of stations and influence of their spatial density are addressed by comparison to a satellite-only climatology. Applying this algorithm to data from automated weather stations and manned synoptic stations in and south of the Atlas Mountains reveals the frequent occurrence. On the order of 6 events per month are detected from May to September when the SHL is in its northernmost position. The events tend to cluster into several-days long convectively active periods, often with strong events on consecutive days. This study is the first to diagnose dynamical impacts of such periods on the SHL, based on simulations of two example cases using the Weather Research and Forecast (WRF) model at convection-permitting resolution. Sensitivity experiments with artificially removed cold pools as well as different resolutions and parameterizations are conducted. Results indicate increases in surface pressure of more than 1 hPa and significant moisture transports into the desert over several days. This moisture affects radiative heating and thus the energy balance of the SHL. Even though cold pool events north of the SHL are less frequent when compared to their Sahelian counterparts, it is shown that they gain importance due to their temporal clustering on synoptic timescale. Together with studies focusing on the Sahel, this work emphasizes the need for improved parameterization schemes for deep convection in order to produce more reliable climate projections for the WAM.

Heat fluxes between the Guadalquivir river and the Gulf of Cádiz Continental Shelf

An 18-year time series of daily sea surface temperature of Gulf of Cadiz and an 18-month time series of temperature collected in the vicinity of the Guadalquivir estuary mouth have been analyzed to investigate the heat exchange between the estuary and the adjacent continental shelf. The first time identifies a continental shelf area where seasonal thermal oscillation signal (amplitudes and phase) changes abruptly. In order to explain this anomaly, the second data set allows a description of thermal fluctuations in a wide range of frequencies and an estimation of the upstream heat budget of the Guadalquivir estuary. Results show that high frequency thermal signal, diurnal and semidiurnal, and water flux signal through Guadalquivir mouth, mainly semidiurnal, apparently interact randomly to give a small exchange of thermal energy at high frequency. There is no trace, at the estuary's mouth, of daily heat exchanges with intertidal mudflats probably because it tends to cancel on daily time scales. Results also show that fluctuations of estimated air-sea fluxes force fluctuations of temperature in a quite homogeneous estuarine, with a delay of 20 days. The along-channel thermal energy gradient reaches magnitudes of 300-400 J m-4 near the mouth during the summer and winter and drives the estuary-shelf exchange of thermal energy at seasonal scale. Particularly, the thermal heat imported by the estuary from the shelf area during late fall-winter-early spring of 2008/2009 is balanced by the thermal heat that the estuary exports to the shelf area during late spring-summer of 2008. In summary, Guadalquivir river removes/imports excess of thermal energy towards/from the continental shelf seasonally, as a mechanism to accommodate excess of heat from one side respect to the other side.

A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.

Characterization of commercial double-walled carbon nanotube material : composition, structure, and heat capacity

A purified commercial double-walled carbon nanotube (DWCNT) sample was investigated by transmission electron microscopy (TEM), thermogravimetry (TG), and Raman spectroscopy. Moreover, the heat capacity of the DWCNT sample was determined by temperature-modulated differential scanning calorimetry in the range of temperature between -50 and 290 °C. The main thermo-oxidation characterized by TG occurred at 474 °C with the loss of 90 wt% of the sample. Thermo-oxidation of the sample was also investigated by high-resolution TG, which indicated that a fraction rich in carbon nanotube represents more than 80 wt% of the material. Other carbonaceous fractions rich in amorphous coating and graphitic particles were identified by the deconvolution procedure applied to the derivative of TG curve. Complementary structural data were provided by TEM and Raman studies. The information obtained allows the optimization of composites based on this nanomaterial with reliable characteristics.

A quasi stationary approach to drying kinetics of cylindrical food particulates in a heat pump asisted fluid bed dryer

Experiments were undertaken to study drying kinetics of moist cylindrical shaped food particulates during fluidised bed drying. Cylindrical particles were prepared from Green beans with three different length:diameter ratios, 3:1, 2:1 and 1:1. A batch fluidised bed dryer connected to a heat pump system was used for the experimentation. A Heat pump and fluid bed combination was used to increase overall energy efficiency and achieve higher drying rates. Drying kinetics, were evaluated with non-dimensional moisture at three different drying temperatures of 30, 40 and 50o C. Numerous mathematical models can be used to calculate drying kinetics ranging from analytical models with simplified assumptions to empirical models built by regression using experimental data. Empirical models are commonly used for various food materials due to their simpler approach. However problems in accuracy, limits the applications of empirical models. Some limitations of empirical models could be reduced by using semi-empirical models based on heat and mass transfer of the drying operation. One such method is the quasi-stationary approach. In this study, a modified quasi-stationary approach was used to model drying kinetics of the cylindrical food particles at three drying temperatures.

Preventing physical activity induced heat illness in school settings

The climatic conditions of tropical and subtropical regions within Australia present, at times, extreme risk of physical activity induced heat illness. Many administrators and teachers in school settings are aware of the general risks of heat related illness. In the absence of reliable information applied at the local level, there is a risk that inappropriate decisions may be made concerning school events that incorporate opportunities to be physically active. Such events may be prematurely cancelled resulting in the loss of necessary time for physical activity. Under high or extremely high risk conditions however, the absence of appropriate modifications or continuation could place the health of students, staff and other parties at risk. School staff and other key stakeholders should understand the mechanisms of escalating risk and be supported to undertake action to reduce the level of risk through appropriate policies, procedures, resources and action plans.

High-flux ceramic membranes with a nanomesh of metal oxide nanofibers

Traditional ceramic separation membranes, which are fabricated by applying colloidal suspensions of metal hydroxides to porous supports, tend to suffer from pinholes and cracks that seriously affect their quality. Other intrinsic problems for these membranes include dramatic losses of flux when the pore sizes are reduced to enhance selectivity and dead-end pores that make no contribution to filtration. In this work, we propose a new strategy for addressing these problems by constructing a hierarchically structured separation layer on a porous substrate using large titanate nanofibers and smaller boehmite nanofibers. The nanofibers are able to divide large voids into smaller ones without forming dead-end pores and with the minimum reduction of the total void volume. The separation layer of nanofibers has a porosity of over 70% of its volume, whereas the separation layer in conventional ceramic membranes has a porosity below 36% and inevitably includes dead-end pores that make no contribution to the flux. This radical change in membrane texture greatly enhances membrane performance. The resulting membranes were able to filter out 95.3% of 60-nm particles from a 0.01 wt % latex while maintaining a relatively high flux of between 800 and 1000 L/m2·h, under a low driving pressure (20 kPa). Such flow rates are orders of magnitude greater than those of conventional membranes with equal selectivity. Moreover, the flux was stable at approximately 800 L/m2·h with a selectivity of more than 95%, even after six repeated runs of filtration and calcination. Use of different supports, either porous glass or porous alumina, had no substantial effect on the performance of the membranes; thus, it is possible to construct the membranes from a variety of supports without compromising functionality. The Darcy equation satisfactorily describes the correlation between the filtration flux and the structural parameters of the new membranes. The assembly of nanofiber meshes to combine high flux with excellent selectivity is an exciting new direction in membrane fabrication.

Bayesian latent trait modeling of migraine symptom data

Definition of disease phenotype is a necessary preliminary to research into genetic causes of a complex disease. Clinical diagnosis of migraine is currently based on diagnostic criteria developed by the International Headache Society. Previously, we examined the natural clustering of these diagnostic symptoms using latent class analysis (LCA) and found that a four-class model was preferred. However, the classes can be ordered such that all symptoms progressively intensify, suggesting that a single continuous variable representing disease severity may provide a better model. Here, we compare two models: item response theory and LCA, each constructed within a Bayesian context. A deviance information criterion is used to assess model fit. We phenotyped our population sample using these models, estimated heritability and conducted genome-wide linkage analysis using Merlin-qtl. LCA with four classes was again preferred. After transformation, phenotypic trait values derived from both models are highly correlated (correlation = 0.99) and consequently results from subsequent genetic analyses were similar. Heritability was estimated at 0.37, while multipoint linkage analysis produced genome-wide significant linkage to chromosome 7q31-q33 and suggestive linkage to chromosomes 1 and 2. We argue that such continuous measures are a powerful tool for identifying genes contributing to migraine susceptibility.

Latent variable models in statistical genetics

Understanding the complexities that are involved in the genetics of multifactorial diseases is still a monumental task. In addition to environmental factors that can influence the risk of disease, there is also a number of other complicating factors. Genetic variants associated with age of disease onset may be different from those variants associated with overall risk of disease, and variants may be located in positions that are not consistent with the traditional protein coding genetic paradigm. Latent Variable Models are well suited for the analysis of genetic data. A latent variable is one that we do not directly observe, but which is believed to exist or is included for computational or analytic convenience in a model. This thesis presents a mixture of methodological developments utilising latent variables, and results from case studies in genetic epidemiology and comparative genomics. Epidemiological studies have identified a number of environmental risk factors for appendicitis, but the disease aetiology of this oft thought useless vestige remains largely a mystery. The effects of smoking on other gastrointestinal disorders are well documented, and in light of this, the thesis investigates the association between smoking and appendicitis through the use of latent variables. By utilising data from a large Australian twin study questionnaire as both cohort and case-control, evidence is found for the association between tobacco smoking and appendicitis. Twin and family studies have also found evidence for the role of heredity in the risk of appendicitis. Results from previous studies are extended here to estimate the heritability of age-at-onset and account for the eect of smoking. This thesis presents a novel approach for performing a genome-wide variance components linkage analysis on transformed residuals from a Cox regression. This method finds evidence for a dierent subset of genes responsible for variation in age at onset than those associated with overall risk of appendicitis. Motivated by increasing evidence of functional activity in regions of the genome once thought of as evolutionary graveyards, this thesis develops a generalisation to the Bayesian multiple changepoint model on aligned DNA sequences for more than two species. This sensitive technique is applied to evaluating the distributions of evolutionary rates, with the finding that they are much more complex than previously apparent. We show strong evidence for at least 9 well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least 7 classes in an alignment of four mammals, including human. A pattern of enrichment and depletion of genic regions in the profiled segments suggests they are functionally significant, and most likely consist of various functional classes. Furthermore, a method of incorporating alignment characteristics representative of function such as GC content and type of mutation into the segmentation model is developed within this thesis. Evidence of fine-structured segmental variation is presented.