Data-Driven Dynamic Modeling of Coupled Thermal and Electric Outputs of Microturbines

Microturbines are among the most successfully commercialized distributed energy resources, especially when they are used for combined heat and power generation. However, the interrelated thermal and electrical system dynamic behaviors have not been fully investigated. This is technically challenging due to the complex thermo-fluid-mechanical energy conversion processes which introduce multiple time-scale dynamics and strong nonlinearity into the analysis. To tackle this problem, this paper proposes a simplified model which can predict the coupled thermal and electric output dynamics of microturbines. Considering the time-scale difference of various dynamic processes occuring within microturbines, the electromechanical subsystem is treated as a fast quasi-linear process while the thermo-mechanical subsystem is treated as a slow process with high nonlinearity. A three-stage subspace identification method is utilized to capture the dominant dynamics and predict the electric power output. For the thermo-mechanical process, a radial basis function model trained by the particle swarm optimization method is employed to handle the strong nonlinear characteristics. Experimental tests on a Capstone C30 microturbine show that the proposed modeling method can well capture the system dynamics and produce a good prediction of the coupled thermal and electric outputs in various operating modes.

Thomson scattering on non-thermal atmospheric pressure plasma jets

To characterize non-thermal atmospheric pressure plasmas experimentally, a large variety of methods and techniques is available, each having its own specific possibilities and limitations. A rewarding method to investigate these plasma sources is laser Thomson scattering. However, that is challenging. Non-thermal atmospheric pressure plasmas (gas temperatures close to room temperature and electron temperatures of a few eV) have usually small dimensions (below 1 mm) and a low degree of ionization (below 10^-4). Here an overview is presented of how Thomson scattering can be applied to such plasmas and used to measure directly spatially and temporally resolved the electron density and energy distribution. A general description of the scattering of photons and the guidelines for an experimental setup of this active diagnostic are provided. Special attention is given to the design concepts required to achieve the maximum signal photon flux with a minimum of unwanted signals. Recent results from the literature are also presented and discussed.

Development of a Rational Characterization Method for Iowa Fly Ash, HR-286, Annual Progress Report 2, 1987

The following report summarizes research activities on the project for the period December 1, 1986 to November 30, 1987. Research efforts for the second year deviated slightly from those described in the project proposal. By the end of the second year of testing, it was possible to begin evaluating how power plant operating conditions influenced the chemical and physical properties of fly ash obtained from one of the monitored power plants (Ottumwa Generating Station, OGS). Hence, several of the tasks initially assigned to the third year of the project (specifically tasks D, E, and F) were initiated during the second year of the project. Manpower constraints were balanced by delaying full scale implementation of the quantitative X-ray diffraction and differential thermal analysis tasks until the beginning of the third year of the project. Such changes should have little bearing on the outcome of the overall project.

Dual stimulation of antigen presenting cells using carbon nanotube-based vaccine delivery system for cancer immunotherapy

Although anti−cancer immuno−based combinatorial therapeutic approaches have shown promising results, efficient tumour eradication demands further intensification of anti−tumour immune response. With the emerging field of nanovaccinology, multi−walled carbon nanotubes (MWNTs) have manifested prominent potentials as tumour antigen nanocarriers. Nevertheless, the utilization of MWNTs in co−delivering antigen along with different types of immunoadjuvants to antigen presenting cells (APCs) has not been investigated yet. We hypothesized that harnessing MWNT for concurrent delivery of cytosine−phosphate−guanine oligodeoxynucleotide (CpG) and anti-CD40 Ig (αCD40), as immunoadjuvants, along with the model antigen ovalbumin (OVA) could potentiate immune response induced against OVA−expressing tumour cells. We initially investigated the effective method to co−deliver OVA and CpG using MWNT to the APC. Covalent conjugation of OVA and CpG prior to loading onto MWNTs markedly augmented the CpG−mediated adjuvanticity, as demonstrated by the significantly increased OVA−specific T cell responses in vitro and in C57BL/6 mice. αCD40 was then included as a second immunoadjuvant to further intensify the immune response. Immune response elicited in vitro and in vivo by OVA, CpG and αCD40 was significantly potentiated by their co−incorporation onto the MWNTs. Furthermore, MWNT remarkably improved the ability of co−loaded OVA, CpG and αCD40 in inhibiting the growth of OVA−expressing B16F10 melanoma cells in subcutaneous or lung pseudo−metastatic tumour models. Therefore, this study suggests that the utilization of MWNTs for the co−delivery of tumour−derived antigen, CpG and αCD40 could be a competent approach for efficient tumours eradication.

A finite-volume shallow layer method for the MHD instabilities in an aluminium production cell

An electrolytic cell for Aluminium production contains molten metal and molten electrolyte, which are subject to high dc-currents and magnetic fields. Lorentz forces arising from the cross product of current and magnetic field may amplify natural gravity waves at the interface between the two fluids, leading to short circuits in extreme cases. The external magnetic field and current distribution in the production cell is computed through a detailed finite element analysis at Torino Polytechnic. The results are then used to compute the magnetohydrodynamic and thermal effects in the aluminium/electrolyte bath. Each cell has lateral dimensions of 6m x 2m, whilst the bath depth is only 30cm. the electrically resistive electrolyte path, which is critical in the operation of the cell, has layer depth of only a few centimetres below each carbon anode. Because the shallow dimensions of the liquid layer a finite-volume shallow-layer technique has been used at Greenwich to compute the resulting flow-field and interface perturbations. The information obtained from this method, i.e. depth averaged velocities and aluminium/electrolyte interface position is then embedded in the three-dimensional finite volume code PHYSICA and will be used to compute the heat transfer and phase change in the cell.

Detection of mesoscale thermal fronts from 4km data using smoothing techniques: Gradient-based fronts classification and basin scale application

In order to optimize frontal detection in sea surface temperature fields at 4 km resolution, a combined statistical and expert-based approach is applied to test different spatial smoothing of the data prior to the detection process. Fronts are usually detected at 1 km resolution using the histogram-based, single image edge detection (SIED) algorithm developed by Cayula and Cornillon in 1992, with a standard preliminary smoothing using a median filter and a 3 × 3 pixel kernel. Here, detections are performed in three study regions (off Morocco, the Mozambique Channel, and north-western Australia) and across the Indian Ocean basin using the combination of multiple windows (CMW) method developed by Nieto, Demarcq and McClatchie in 2012 which improves on the original Cayula and Cornillon algorithm. Detections at 4 km and 1 km of resolution are compared. Fronts are divided in two intensity classes (“weak” and “strong”) according to their thermal gradient. A preliminary smoothing is applied prior to the detection using different convolutions: three type of filters (median, average and Gaussian) combined with four kernel sizes (3 × 3, 5 × 5, 7 × 7, and 9 × 9 pixels) and three detection window sizes (16 × 16, 24 × 24 and 32 × 32 pixels) to test the effect of these smoothing combinations on reducing the background noise of the data and therefore on improving the frontal detection. The performance of the combinations on 4 km data are evaluated using two criteria: detection efficiency and front length. We find that the optimal combination of preliminary smoothing parameters in enhancing detection efficiency and preserving front length includes a median filter, a 16 × 16 pixel window size, and a 5 × 5 pixel kernel for strong fronts and a 7 × 7 pixel kernel for weak fronts. Results show an improvement in detection performance (from largest to smallest window size) of 71% for strong fronts and 120% for weak fronts. Despite the small window used (16 × 16 pixels), the length of the fronts has been preserved relative to that found with 1 km data. This optimal preliminary smoothing and the CMW detection algorithm on 4 km sea surface temperature data are then used to describe the spatial distribution of the monthly frequencies of occurrence for both strong and weak fronts across the Indian Ocean basin. In general strong fronts are observed in coastal areas whereas weak fronts, with some seasonal exceptions, are mainly located in the open ocean. This study shows that adequate noise reduction done by a preliminary smoothing of the data considerably improves the frontal detection efficiency as well as the global quality of the results. Consequently, the use of 4 km data enables frontal detections similar to 1 km data (using a standard median 3 × 3 convolution) in terms of detectability, length and location. This method, using 4 km data is easily applicable to large regions or at the global scale with far less constraints of data manipulation and processing time relative to 1 km data.

Thermal activation, long-range ordering, and topological frustration of artificial spin ice

Frustrated systems, typically characterized by competing interactions that cannot all be simultaneously satisfied, are ubiquitous in nature and display many rich phenomena and novel physics. Artificial spin ices (ASIs), arrays of lithographically patterned Ising-like single-domain magnetic nanostructures, are highly tunable systems that have proven to be a novel method for studying the effects of frustration and associated properties. The strength and nature of the frustrated interactions between individual magnets are readily tuned by design and the exact microstate of the system can be determined by a variety of characterization techniques. Recently, thermal activation of ASI systems has been demonstrated, introducing the spontaneous reversal of individual magnets and allowing for new explorations of novel phase transitions and phenomena using these systems. In this work, we introduce a new, robust material with favorable magnetic properties for studying thermally active ASI and use it to investigate a variety of ASI geometries. We reproduce previously reported perfect ground-state ordering in the square geometry and present studies of the kagome lattice showing the highest yet degree of ordering observed in this fully frustrated system. We consider theoretical predictions of long-range order in ASI and use both our experimental studies and kinetic Monte Carlo simulations to evaluate these predictions. Next, we introduce controlled topological defects into our square ASI samples and observe a new, extended frustration effect of the system. When we introduce a dislocation into the lattice, we still see large domains of ground-state order, but, in every sample, a domain wall containing higher energy spin arrangements originates from the dislocation, resolving a discontinuity in the ground-state order parameter. Locally, the magnets are unfrustrated, but frustration of the lattice persists due to its topology. We demonstrate the first direct imaging of spin configurations resulting from topological frustration in any system and make predictions on how dislocations could affect properties in numerous materials systems.

Calcium response of KCl-excited populations of ventricular myocytes from the European sea bass (Dicentrarchus labrax): a promising approach to integrate cell-to-cell heterogeneity in studying the cellular basis of fish cardiac performance

Climate change challenges the capacity of fishes to thrive in their habitat. However, through phenotypic diversity, they demonstrate remarkable resilience to deteriorating conditions. In fish populations, inter-individual variation in a number of fitness-determining physiological traits, including cardiac performance, is classically observed. Information about the cellular bases of inter-individual variability in cardiac performance is scarce including the possible contribution of excitation-contraction (EC) coupling. This study aimed at providing insight into EC coupling-related Ca2+ response and thermal plasticity in the European sea bass (Dicentrarchus labrax). A cell population approach was used to lay the methodological basis for identifying the cellular determinants of cardiac performance. Fish were acclimated at 12 and 22 A degrees C and changes in intracellular calcium concentration ([Ca2+](i)) following KCl stimulation were measured using Fura-2, at 12 or 22 A degrees C-test. The increase in [Ca2+](i) resulted primarily from extracellular Ca2+ entry but sarcoplasmic reticulum stores were also shown to be involved. As previously reported in sea bass, a modest effect of adrenaline was observed. Moreover, although the response appeared relatively insensitive to an acute temperature change, a difference in Ca2+ response was observed between 12- and 22 A degrees C-acclimated fish. In particular, a greater increase in [Ca2+](i) at a high level of adrenaline was observed in 22 A degrees C-acclimated fish that may be related to an improved efficiency of adrenaline under these conditions. In conclusion, this method allows a rapid screening of cellular characteristics. It represents a promising tool to identify the cellular determinants of inter-individual variability in fishes' capacity for environmental adaptation.

FLOW REGIME DRIVEN THERMAL ENHANCEMENT IN INTERNALLY-GROOVED TUBES

Internally-grooved refrigeration tubes maximize tube-side evaporative heat transfer rates and have been identified as a most promising technology for integration into compact cold plates. Unfortunately, the absence of phenomenological insights and physical models hinders the extrapolation of grooved-tube performance to new applications. The success of regime-based heat transfer correlations for smooth tubes has motivated the current effort to explore the relationship between flow regimes and enhanced heat transfer in internally-grooved tubes. In this thesis, a detailed analysis of smooth and internally-grooved tube data reveals that performance improvement in internally-grooved tubes at low-to-intermediate mass flux is a result of early flow regime transition. Based on this analysis, a new flow regime map and corresponding heat transfer coefficient correlation, which account for the increased wetted angle, turbulence, and Gregorig effects unique to internally-grooved tubes, were developed. A two-phase test facility was designed and fabricated to validate the newly-developed flow regime map and regime-based heat transfer coefficient correlation. As part of this setup, a non-intrusive optical technique was developed to study the dynamic nature of two-phase flows. It was found that different flow regimes result in unique temporally varying film thickness profiles. Using these profiles, quantitative flow regime identification measures were developed, including the ability to explain and quantify the more subtle transitions that exist between dominant flow regimes. Flow regime data, based on the newly-developed method, and heat transfer coefficient data, using infrared thermography, were collected for two-phase HFE-7100 flow in horizontal 2.62mm - 8.84mm diameter smooth and internally-grooved tubes with mass fluxes from 25-300 kg/m²s, heat fluxes from 4-56 kW/m², and vapor qualities approaching 1. In total, over 6500 combined data points for the adiabatic and diabatic smooth and internally-grooved tubes were acquired. Based on results from the experiments and a reinterpretation of data from independent researchers, it was established that heat transfer enhancement in internally-grooved tubes at low-to-intermediate mass flux is primarily due to early flow regime transition to Annular flow. The regime-based heat transfer coefficient outperformed empirical correlations from the literature, with mean and absolute deviations of 4.0% and 32% for the full range of data collected.

COMPREHENSIVE ELECTRICAL/OPTICAL/THERMAL CHARACTERIZATIONS OF HIGH POWER LIGHT EMITTING DIODES AND LASER DIODES

Thermal characterizations of high power light emitting diodes (LEDs) and laser diodes (LDs) are one of the most critical issues to achieve optimal performance such as center wavelength, spectrum, power efficiency, and reliability. Unique electrical/optical/thermal characterizations are proposed to analyze the complex thermal issues of high power LEDs and LDs. First, an advanced inverse approach, based on the transient junction temperature behavior, is proposed and implemented to quantify the resistance of the die-attach thermal interface (DTI) in high power LEDs. A hybrid analytical/numerical model is utilized to determine an approximate transient junction temperature behavior, which is governed predominantly by the resistance of the DTI. Then, an accurate value of the resistance of the DTI is determined inversely from the experimental data over the predetermined transient time domain using numerical modeling. Secondly, the effect of junction temperature on heat dissipation of high power LEDs is investigated. The theoretical aspect of junction temperature dependency of two major parameters – the forward voltage and the radiant flux – on heat dissipation is reviewed. Actual measurements of the heat dissipation over a wide range of junction temperatures are followed to quantify the effect of the parameters using commercially available LEDs. An empirical model of heat dissipation is proposed for applications in practice. Finally, a hybrid experimental/numerical method is proposed to predict the junction temperature distribution of a high power LD bar. A commercial water-cooled LD bar is used to present the proposed method. A unique experimental setup is developed and implemented to measure the average junction temperatures of the LD bar. After measuring the heat dissipation of the LD bar, the effective heat transfer coefficient of the cooling system is determined inversely. The characterized properties are used to predict the junction temperature distribution over the LD bar under high operating currents. The results are presented in conjunction with the wall-plug efficiency and the center wavelength shift.

Thermal analysis during bone drilling using rigid polyurethane foams: numerical and experimental methodologies

Osteotomy or bone cutting is a common procedure in orthopaedic surgery, mainly in the treatment of fractures and reconstructive surgery. However, the excessive heat produced during the bone drilling process is a problem that counters the benefits of this type of surgery, because it can result in thermal osteonecrosis, bone reabsorption and damage the osseointegration of implants. The analysis of different drilling parameters and materials can allow to decrease the temperature during the bone drilling process and contribute to a greater success of this kind of surgical interventions. The main goal of this study was to build a numerical three-dimensional model to simulate the drilling process considering the type of bone, the influence of cooling and the bone density of the different composite materials with similar mechanical properties to the human bone and generally used in experimental biomechanics. The numerical methodology was coupled with an experimental methodology. The use of cooling proved to be essential to decrease the material damage during the drilling process. It was concluded that the materials with less porosity and density present less damage in drilling process. The developed numerical model proved to be a great tool in this kind of analysis. © 2016, The Brazilian Society of Mechanical Sciences and Engineering.

Self-sealing of thermal fatigue and mechanical damage in fiber-reinforced composite materials

Fiber reinforced composite tanks provide a promising method of storage for liquid oxygen and hydrogen for aerospace applications. The inherent thermal fatigue of these vessels leads to the formation of microcracks, which allow gas phase leakage across the tank walls. In this dissertation, self-healing functionality is imparted to a structural composite to effectively seal microcracks induced by both mechanical and thermal loading cycles. Two different microencapsulated healing chemistries are investigated in woven glass fiber/epoxy and uni-weave carbon fiber/epoxy composites. Self-healing of mechanically induced damage was first studied in a room temperature cured plain weave E-glass/epoxy composite with encapsulated dicyclopentadiene (DCPD) monomer and wax protected Grubbs' catalyst healing components. A controlled amount of microcracking was introduced through cyclic indentation of opposing surfaces of the composite. The resulting damage zone was proportional to the indentation load. Healing was assessed through the use of a pressure cell apparatus to detect nitrogen flow through the thickness direction of the damaged composite. Successful healing resulted in a perfect seal, with no measurable gas flow. The effect of DCPD microcapsule size (51 um and 18 um) and concentration (0 - 12.2 wt%) on the self-sealing ability was investigated. Composite specimens with 6.5 wt% 51 um capsules sealed 67% of the time, compared to 13% for the control panels without healing components. A thermally stable, dual microcapsule healing chemistry comprised of silanol terminated poly(dimethyl siloxane) plus a crosslinking agent and a tin catalyst was employed to allow higher composite processing temperatures. The microcapsules were incorporated into a satin weave E-glass fiber/epoxy composite processed at 120C to yield a glass transition temperature of 127C. Self-sealing ability after mechanical damage was assessed for different microcapsule sizes (25 um and 42 um) and concentrations (0 - 11 vol%). Incorporating 9 vol% 42 um capsules or 11 vol% 25 um capsules into the composite matrix leads to 100% of the samples sealing. The effect of microcapsule concentration on the short beam strength, storage modulus, and glass transition temperature of the composite specimens was also investigated. The thermally stable tin catalyzed poly(dimethyl siloxane) healing chemistry was then integrated into a [0/90]s uniweave carbon fiber/epoxy composite. Thermal cycling (-196C to 35C) of these specimens lead to the formation of microcracks, over time, formed a percolating crack network from one side of the composite to the other, resulting in a gas permeable specimen. Crack damage accumulation and sample permeability was monitored with number of cycles for both self-healing and traditional non-healing composites. Crack accumulation occurred at a similar rate for all sample types tested. A 63% increase in lifetime extension was achieved for the self-healing specimens over traditional non-healing composites.

Method of numerical correction of errors occasioned by delay of records during the monitoring of environmental variables of interest for animal production

As condições de ambiente térmico e aéreo, no interior de instalações para animais, alteram-se durante o dia, devido à influência do ambiente externo. Para que análises estatísticas e geoestatísticas sejam representativas, uma grande quantidade de pontos distribuídos espacialmente na área da instalação deve ser monitorada. Este trabalho propõe que a variação no tempo das variáveis ambientais de interesse para a produção animal, monitoradas no interior de instalações para animais, pode ser modelada com precisão a partir de registros discretos no tempo. O objetivo deste trabalho foi desenvolver um método numérico para corrigir as variações temporais dessas variáveis ambientais, transformando os dados para que tais observações independam do tempo gasto durante a aferição. O método proposto aproximou os valores registrados com retardos de tempo aos esperados no exato momento de interesse, caso os dados fossem medidos simultaneamente neste momento em todos os pontos distribuídos espacialmente. O modelo de correção numérica para variáveis ambientais foi validado para o parâmetro ambiental temperatura do ar, sendo que os valores corrigidos pelo método não diferiram pelo teste Tukey, a 5% de probabilidade dos valores reais registrados por meio de dataloggers.

A new, fast and flexible radiative transfer method for Venus general circulation models

We present a new radiation scheme for the Oxford Planetary Unified Model System for Venus, suitable for the solar and thermal bands. This new and fast radiative parameterization uses a different approach in the two main radiative wavelength bands: solar radiation (0.1-5.5 mu m) and thermal radiation (1.7-260 mu m). The solar radiation calculation is based on the delta-Eddington approximation (two-stream-type) with an adding layer method. For the thermal radiation case, a code based on an absorptivity/emissivity formulation is used. The new radiative transfer formulation implemented is intended to be computationally light, to allow its incorporation in 3D global circulation models, but still allowing for the calculation of the effect of atmospheric conditions on radiative fluxes. This will allow us to investigate the dynamical-radiative-microphysical feedbacks. The model flexibility can be also used to explore the uncertainties in the Venus atmosphere such as the optical properties in the deep atmosphere or cloud amount. The results of radiative cooling and heating rates and the global-mean radiative-convective equilibrium temperature profiles for different atmospheric conditions are presented and discussed. This new scheme works in an atmospheric column and can be easily implemented in 3D Venus global circulation models. (C) 2014 Elsevier Ltd. All rights reserved.