Thermal Energy Storage for Solar Power Production

Solar energy is the most abundant persistent energy resource. It is also an intermittent one available for only a fraction of each day while the demand for electric power never ceases. To produce a significant amount of power at the utility scale, electricity generated from solar energy must be dispatchable and able to be supplied in response to variations in demand. This requires energy storage that serves to decouple the intermittent solar resource from the load and enables around-the-clock power production from solar energy. Practically, solar energy storage technologies must be efficient as any energy loss results in an increase in the amount of required collection hardware, the largest cost in a solar electric power system. Storing solar energy as heat has been shown to be an efficient, scalable, and relatively low-cost approach to providing dispatchable solar electricity. Concentrating solar power systems that include thermal energy storage (TES) use mirrors to focus sunlight onto a heat exchanger where it is converted to thermal energy that is carried away by a heat transfer fluid and used to drive a conventional thermal power cycle (e.g., steam power plant), or stored for later use. Several approaches to TES have been developed and can generally be categorized as either thermophysical (wherein energy is stored in a hot fluid or solid medium or by causing a phase change that can later be reversed to release heat) or thermochemical (in which energy is stored in chemical bonds requiring two or more reversible chemical reactions).

Design and Preliminary Results of an Atmospheric-Pressure Model Gas Turbine Combustor Utilizing Varying CO2 Doping Concentration in CH4 to Emulate Biogas Combustion

Utilization of biogas can provide a source of renewable energy in both heat and power generation. Combustion of biogas in land-based gas turbines for power generation is a promising approach to reducing greenhouse gases and US dependence on foreign-source fossil fuels. Biogas is a byproduct from the decomposition of organic matter and consists primarily of CH4 and large amounts of CO2. The focus of this research was to design a combustion device and investigate the effects of increasing levels of CO2 addition to the combustion of pure CH4 with air. Using an atmospheric-pressure, swirl-stabilized dump combustor, emissions data and flame stability limitations were measured and analyzed. In particular, CO2, CO, and NOx emissions were the main focus of the combustion products. Additionally, the occurrence of lean blowout and combustion pressure oscillations, which impose significant limitations in operation ranges for actual gas turbines, was observed. Preliminary kinetic and equilibrium modeling was performed using Cantera and CEA for the CH4/CO2/Air combustion systems to analyze the effect of CO2 upon adiabatic flame temperature and emission levels. The numerical and experimental results show similar dependence of emissions on equivalence ratio, CO2 addition, inlet air temperature, and combustor residence time. (C) 2014 Elsevier Ltd. All rights reserved.

Performance characteristics and design recommendations for biomass-burning stoves using earthen construction materials

This report provides an analysis of the thermal performance and emissions characteristics of improved biomass stoves constructed using earthen materials. Commonly referred to as mud stoves, this type of improved stove incorporates high clay content soil with an organic binder in the construction of its combustion chamber and body. When large quantities of the mud material are used to construct the stove body, the stove does not offer significant improvements in fuel economy or air quality relative to traditional open fire cooking. This is partly because a significant amount of heat is absorbed by the mass of the stove reducing combustion efficiency and heat transfer to the cook pot. An analysis of the thermal and mechanical properties of stove materials was also performed. A material mixture containing a one‐to‐one ratio by volume of high content clay soil and straw was found to have thermal properties comparable to fired ceramics used in more advanced improved stove designs. Feedback from mud stove users in Mauritania and Mali, West Africa was also collected during implementation. Suggestions for stove design improvements were developed based on this information and the data collected in the performance, emissions, and material properties analysis. Design suggestions include reducing stove height to accommodate user cooking preferences and limiting overall stove mass to reduce heat loss to the stove body.

Experimental results on gravity driven fully condensing flows in vertical tubes, their agreement with theory, and their differences with shear driven flows' boundary-condition sensitivities

This doctoral thesis presents the experimental results along with a suitable synthesis with computational/theoretical results towards development of a reliable heat transfer correlation for a specific annular condensation flow regime inside a vertical tube. For fully condensing flows of pure vapor (FC-72) inside a vertical cylindrical tube of 6.6 mm diameter and 0.7 m length, the experimental measurements are shown to yield values of average heat transfer co-efficient, and approximate length of full condensation. The experimental conditions cover: mass flux G over a range of 2.9 kg/m2-s ≤ G ≤ 87.7 kg/m2-s, temperature difference ∆T (saturation temperature at the inlet pressure minus the mean condensing surface temperature) of 5 ºC to 45 ºC, and cases for which the length of full condensation xFC is in the range of 0 < xFC < 0.7 m. The range of flow conditions over which there is good agreement (within 15%) with the theory and its modeling assumptions has been identified. Additionally, the ranges of flow conditions for which there are significant discrepancies (between 15 -30% and greater than 30%) with theory have also been identified. The paper also refers to a brief set of key experimental results with regard to sensitivity of the flow to time-varying or quasi-steady (i.e. steady in the mean) impositions of pressure at both the inlet and the outlet. The experimental results support the updated theoretical/computational results that gravity dominated condensing flows do not allow such elliptic impositions.

Development of one-dimensional and two-dimensional computational tools that simulate steady internal condensing flows in terrestrial and zero-gravity environments

This dissertation presents an effective quasi one-dimensional (1-D) computational simulation tool and a full two-dimensional (2-D) computational simulation methodology for steady annular/stratified internal condensing flows of pure vapor. These simulation tools are used to investigate internal condensing flows in both gravity as well as shear driven environments. Through accurate numerical simulations of the full two dimensional governing equations, results for laminar/laminar condensing flows inside mm-scale ducts are presented. The methodology has been developed using MATLAB/COMSOL platform and is currently capable of simulating film-wise condensation for steady (and unsteady flows). Moreover, a novel 1-D solution technique, capable of simulating condensing flows inside rectangular and circular ducts with different thermal boundary conditions is also presented. The results obtained from the 2-D scientific tool and 1-D engineering tool, are validated and synthesized with experimental results for gravity dominated flows inside vertical tube and inclined channel; and, also, for shear/pressure driven flows inside horizontal channels. Furthermore, these simulation tools are employed to demonstrate key differences of physics between gravity dominated and shear/pressure driven flows. A transition map that distinguishes shear driven, gravity driven, and “mixed” driven flow zones within the non-dimensional parameter space that govern these duct flows is presented along with the film thickness and heat transfer correlations that are valid in these zones. It has also been shown that internal condensing flows in a micro-meter scale duct experiences shear driven flow, even in different gravitational environments. The full 2-D steady computational tool has been employed to investigate the length of annularity. The result for a shear driven flow in a horizontal channel shows that in absence of any noise or pressure fluctuation at the inlet, the onset of non-annularity is partly due to insufficient shear at the liquid-vapor interface. This result is being further corroborated/investigated by R. R. Naik with the help of the unsteady simulation tool. The condensing flow results and flow physics understanding developed through these simulation tools will be instrumental in reliable design of modern micro-scale and spacebased thermal systems.

Numerical study on the curling and warping of hardened rigid pavement slabs

In-service hardened concrete pavement suffers from environmental loadings caused by curling and warping of the slab. Traditionally, these loadings are computed on the basis of treating the slab as an elastic material, and of evaluating separately the curling and warping components. This dissertation simulates temperature distribution and moisture distribution through the slabs by use of a developed numerical model that couples the heat transfer and moisture transport. The computation of environmental loadings treats the slab as an elastic-viscous material, which considers the relaxation behavior and Pickett effect of the concrete. The heat transfer model considers the impacts of solar radiation, wind speed, air temperature, pavement slab albedo, etc. on the pavement temperature distribution. This dissertation assesses the difference between documented models that aim to predict pavement temperature, highlighting their pros and cons. The moisture transport model is unique for the documented models; it mimics the wetting and drying events occurring at the slab surface. These events are estimated by a proposed statistical algorithm, which is verified by field rainfall data. Analysis of the predicted results examines on the roles of the local air RH (relative humidity), wind speed, rainy pattern in the moisture distribution through the slab. The findings reveal that seasonal air RH plays a decisive role on the slab‘s moisture distribution; but wind speed and its daily variation, daily RH variation, and seasonal rainfall pattern plays only a secondary role. This dissertation sheds light on the computation of environmental loadings that in-service pavement slabs suffer from. Analysis of the computed stresses centers on the stress relaxation near the surface, stress evolution after the curing ends, and the impact of construction season on the stress‘s magnitude. An unexpected finding is that the total environmental loadings at the cyclically-stable state divert from the thermal stresses. At such a state, the total stress at the daytime is roughly equal to the thermal stress; whereas the total stress during the nighttime is far greater than the thermal stress. An explanation for this phenomenon is that during the night hours, the decline of the slab‘s near-surface temperature leads to a drop of the near-surface RH. This RH drop results in contraction therein and develops additional tensile stresses. The dissertation thus argues that estimating the environmental loadings by solely computing the thermally-induced stresses may reach delusive results. It recommends that the total environmental loadings of in-service slabs should be estimated by a sophisticated model coupling both moisture component and temperature component.

Deletion of CD39 inhibits development of metastatic liver cancer

New vessel formation and tumor infiltrating lymphocytes (TIL) influence host responses to malignant tissues. Extracellular adenosine-mediated pathways promote both vascular endothelial cell proliferation and inhibit cytotoxic T cells, thereby potentiating cancer growth. CD39 is the dominant ectonucleotidase of vascular and T regulatory cells and has the potential to generate high levels of adenosine locally. We have previously shown that deletion of Cd39 results in angiogenic failure and T regulatory cell dysfunction with loss of immune suppressive functions. Aim: Investigate impact of CD39 upon development of hepatic metastases. Methods and Results: We studied the development of metastatic liver deposits following portal vein infusion of 1.5x105 melanoma B16/F10 cells, with luciferase expression, in wild type and Cd39-null C57BL/6 mice (n=24). Tumor formation in liver was directly examined and animals imaged at days 7-17 after tumor cell implantation. As predicted, the formation of hepatic malignant foci was markedly suppressed in Cd39-null mice, at all time points examined. To test whether the major impact of Cd39-deletion was upon neovasculature formation or immune responsiveness, adoptive transfer experiments were conducted. Bone marrow transplants (BMT) from Cd39-null or wild type BL/6 mice were placed in lethally irradiated control and/or null mice, in a crossover manner (total n=24 for each group, respectively). Eight weeks postadoptive transfer, melanoma cells were infused via portal vein as before and tumor growth studied. The Cd39-null mice that received wild type BMT mirrored the wild type phenotype with progressive tumor growth observed (n=8 per time point; p=0.015). In contrast, metastases were significantly inhibited in both number and size and ultimately became necrotic in the wild type mice that had received Cd39-null BMT. Conclusions: Bone marrow derived cells mediate the major inhibitory effects of CD39 deletion on tumor growth. Pharmacological inhibition of CD39 may find utility as an adjunct therapy in the management of hepatic malignancy.

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).

Estimation of vertical groundwater fluxes into a streambed through continuous temperature profile monitoring and the relationship of groundwater fluxes to coaster brook trout spawning habitat

We hypothesized that the spatial distribution of groundwater inflows through river bottom sediments is a critical factor associated with the selection of coaster brook trout (a life history variant of Salvelinus fontinalis,) spawning sites. An 80-m reach of the Salmon Trout River, in the Huron Mountains of the upper peninsula of Michigan, was selected to test the hypothesis based on long-term documentation of coaster brook trout spawning at this site. Throughout this site, the river is relatively similar along its length with regard to stream channel and substrate features. A monitoring well system consisting of an array of 27 wells was installed to measure subsurface temperatures underneath the riverbed over a 13-month period. The monitoring well locations were separated into areas where spawning has and has not been observed. Over 200,000 total temperature measurements were collected from 5 depths within each of the 27 monitoring wells. Temperatures within the substrate at the spawning area were generally cooler and less variable than river temperatures. Substrate temperatures in the non-spawning area were generally warmer, more variable, and closely tracked temporal variations in river temperatures. Temperature data were inverted to obtain subsurface groundwater velocities using a numerical approximation of the heat transfer equation. Approximately 45,000 estimates of groundwater velocities were obtained. Estimated velocities in the spawning and non-spawning areas confirmed that groundwater velocities in the spawning area were primarily in the upward direction, and were generally greater in magnitude than velocities in the non-spawning area. In the non-spawning area there was a greater occurrence of velocities in the downward direction, and velocity estimates were generally lesser in magnitude than in the spawning area. Both the temperature and velocity results confirm the hypothesis that spawning sites correspond to areas of significant groundwater influx to the river bed.

THERMAL CHARACTERIZATION OF COMBUSTION CHAMBER COMPONENTS IN A GASOLINE TURBOCHARGED DIRECT INJECTION (GTDI) ENGINE

This report is a PhD dissertation proposal to study the in-cylinder temperature and heat flux distributions within a gasoline turbocharged direct injection (GTDI) engine. Recent regulations requiring automotive manufacturers to increase the fuel efficiency of their vehicles has led to great technological achievements in internal combustion engines. These achievements have increased the power density of gasoline engines dramatically in the last two decades. Engine technologies such as variable valve timing (VVT), direct injection (DI), and turbocharging have significantly improved engine power-to-weight and power-to-displacement ratios. A popular trend for increasing vehicle fuel economy in recent years has been to downsize the engine and add VVT, DI, and turbocharging technologies so that a lighter more efficient engine can replace a larger, heavier one. With the added power density, thermal management of the engine becomes a more important issue. Engine components are being pushed to their temperature limits. Therefore it has become increasingly important to have a greater understanding of the parameters that affect in-cylinder temperatures and heat transfer. The proposed research will analyze the effects of engine speed, load, relative air-fuel ratio (AFR), and exhaust gas recirculation (EGR) on both in-cylinder and global temperature and heat transfer distributions. Additionally, the effect of knocking combustion and fuel spray impingement will be investigated. The proposed research will be conducted on a 3.5 L six cylinder GTDI engine. The research engine will be instrumented with a large number of sensors to measure in-cylinder temperatures and pressures, as well as, the temperature, pressure, and flow rates of energy streams into and out of the engine. One of the goals of this research is to create a model that will predict the energy distribution to the crankshaft, exhaust, and cooling system based on normalized values for engine speed, load, AFR, and EGR. The results could be used to aid in the engine design phase for turbocharger and cooling system sizing. Additionally, the data collected can be used for validation of engine simulation models, since in-cylinder temperature and heat flux data is not readily available in the literature..

THERMAL CHARACTERIZATION OF A GASOLINE TURBOCHARGED DIRECT INJECTION (GTDI) ENGINE UTILIZING LEAN OPERATION AND EXHAUST GAS RECIRCULATION (EGR)

The push for improved fuel economy and reduced emissions has led to great achievements in engine performance and control. These achievements have increased the efficiency and power density of gasoline engines dramatically in the last two decades. With the added power density, thermal management of the engine has become increasingly important. Therefore it is critical to have accurate temperature and heat transfer models as well as data to validate them. With the recent adoption of the 2025 Corporate Average Fuel Economy(CAFE) standard, there has been a push to improve the thermal efficiency of internal combustion engines even further. Lean and dilute combustion regimes along with waste heat recovery systems are being explored as options for improving efficiency. In order to understand how these technologies will impact engine performance and each other, this research sought to analyze the engine from both a 1st law energy balance perspective, as well as from a 2nd law exergy analysis. This research also provided insights into the effects of various parameters on in-cylinder temperatures and heat transfer as well as provides data for validation of other models. It was found that the engine load was the dominant factor for the energy distribution, with higher loads resulting in lower coolant heat transfer and higher brake work and exhaust energy. From an exergy perspective, the exhaust system provided the best waste heat recovery potential due to its significantly higher temperatures compared to the cooling circuit. EGR and lean combustion both resulted in lower combustion chamber and exhaust temperatures; however, in most cases the increased flow rates resulted in a net increase in the energy in the exhaust. The exhaust exergy, on the other hand, was either increased or decreased depending on the location in the exhaust system and the other operating conditions. The effects of dilution from lean operation and EGR were compared using a dilution ratio, and the results showed that lean operation resulted in a larger increase in efficiency than the same amount of dilution with EGR. Finally, a method for identifying fuel spray impingement from piston surface temperature measurements was found. Note: The material contained in this section is planned for submission as part of a journal article and/or conference paper in the future.

Influence of Mechanical and Thermal Boundary Conditions on Stabilizing/Destabilizing Mechanisms in Evaporating Liquid Films

Liquid films, evaporating or non-evaporating, are ubiquitous in nature and technology. The dynamics of evaporating liquid films is a study applicable in several industries such as water recovery, heat exchangers, crystal growth, drug design etc. The theory describing the dynamics of liquid films crosses several fields such as engineering, mathematics, material science, biophysics and volcanology to name a few. Interfacial instabilities typically manifest by the undulation of an interface from a presumed flat state or by the onset of a secondary flow state from a primary quiescent state or both. To study the instabilities affecting liquid films, an evaporating/non-evaporating Newtonian liquid film is subject to a perturbation. Numerical analysis is conducted on configurations of such liquid films being heated on solid surfaces in order to examine the various stabilizing and destabilizing mechanisms that can cause the formation of different convective structures. These convective structures have implications towards heat transfer that occurs via this process. Certain aspects of this research topic have not received attention, as will be obvious from the literature review. Static, horizontal liquid films on solid surfaces are examined for their resistance to long wave type instabilities via linear stability analysis, method of normal modes and finite difference methods. The spatiotemporal evolution equation, available in literature, describing the time evolution of a liquid film heated on a solid surface, is utilized to analyze various stabilizing/destabilizing mechanisms affecting evaporating and non-evaporating liquid films. The impact of these mechanisms on the film stability and structure for both buoyant and non-buoyant films will be examined by the variation of mechanical and thermal boundary conditions. Films evaporating in zero gravity are studied using the evolution equation. It is found that films that are stable to long wave type instabilities in terrestrial gravity are prone to destabilization via long wave instabilities in zero gravity.

Chromosomal control of the dominant suppression of {\it ras\/}-mediated transformation

The role of tumor suppressor function in the multistep process of carcinogenesis was studied in the human teratocarcinoma cell line PA-1. Early passage PA-1 cells ($<$P100) are preneoplastic while late passage ($>$P100) PA-1 cells are spontaneously transformed. Previous work demonstrated a causal role for the N-ras oncogene in the neoplastic transformation of this cell line and the gene was cloned. A clonal cell line established at passage 40 has been shown to suppress the neoplastic transformation potential of the PA-1 N-ras oncogene in gene transfer experiments. This phenotype has been termed SRT+ for suppression of ras transformation. A clonal cell line established at passage 63 is neoplastically transformed by the N-ras in similar gene transfer experiments and is regarded as srt$-$. Somatic cell hybrids were formed between the SRT+ cell and two different N-ras transformed srt$-$ cells. The results indicate that five of the seven independent hybrid clones, and all 14 subclones, failed to form tumors in the nude mouse tumor assay. Chromosomal analysis of rare neoplastic segregants which arose from suppressed hybrid populations demonstrate that the general loss of chromosomes correlates with the reemergence of neoplastic transformation. Karyotype analyses demonstrate a statistically correlative loss of chromosomes 1, 4, 19, and to a lesser extent 11, 14, and 16. DNA hybridization analysis demonstrates a single copy of the intact N-ras oncogene in parental cells, suppressed hybrids, and neoplastically transformed hybrids. These results indicate that functional ras transformation suppression is a trans-dominant trait which may be controlled by sequences residing on particular chromosomes in the human genome. Furthermore, the suppression of ras transformation results from a unique step in the multistep process of carcinogenesis that is different from the induction of immortality. Thus, the neoplastic process of the PA-1 cell line involves at least three steps: (1) induction of immortality, (2) activation of the N-ras oncogene, and (3) loss of tumor suppressor function. ^

Regulation of sodium/glucose cotransporter expression in LLC-PK$\sb1$ cells

Expression of the Na$\sp+$/glucose cotransporter (SGLT1), a differentiated function of the pig kidney epithelial cell line LLC-PK$\sb1$ derived from proximal tubule, was further investigated. The differentiation inducer hexamethylene bisacetamide (HMBA) and IBMX, an inhibitor of cAMP phosphodiesterase, each stimulated a significant increase in Na$\sp+$/glucose cotransport activity, levels of the 75 kD cotransporter subunit and steady-state levels of the SGLT1 message. The action of HMBA is associated with involvement of polyamines and protein kinase C, and is synergistic with cAMP. We provide evidence that cAMP-elevating agents increase Na$\sp+$/glucose cotransporter expression, at least in part, via a post-transcriptional mechanism. Two molecular species of SGLT1 mRNA (3.9 kb and 2.2 kb) are transcribed from the same gene in LLC-PK$\sb1$ cells and differ only in the length of the 3$\sp\prime$ untranslated region (3$\sp\prime$ UTR). cAMP elevation differentially stabilized the 3.9 kb SGLT1 transcript from degradation but not the 22 kb species. UV-cross-linking and label transfer experiments indicated that cyclic AMP elevation was associated with formation of a 48 kD protein complex with a specific domain within the 3$\sp\prime$ UTR of SGLT1 mRNA. The binding was competitively inhibited by poly (U) and other U-rich RNA species such as c-fos ARE, and modulated by a protein kinase A-mediated phosphorylation/dephosphorylation mechanism. The binding site was mapped to a 120-nucleotide 3$\sp\prime$ UTR sequence which contains a uridine-rich region (URE). Our study provides the first demonstration that renal SGLT1 is post-transcriptionally regulated by a phosphorylation/dephosphorylation mechanism, and provides a deeper insight into gene regulation of this physiologically important cotransporter. ^

Three-dimensional Vortex Structures on Heated Micro-lattice in a Gas

This paper addresses the microscale heat transfer problem from heated lattice to the gas. A micro-device for enhanced heat transfer is presented and numerically investigated. Thermal creep induces 3-D vortex structures in the vicinity of the lattice. The gas flow is in the slip flow regime (Knudsen number Kn⩽0.1Kn⩽0.1). The simulations are performed using slip flow Navier–Stokes equations with boundary condition formulations proposed by Maxwell and Smoluchowski. In this study the wire thicknesses and distances of the heated lattice are varied. The surface geometrical properties alter significantly heat flux through the surface.