Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temperature gradient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

Physics of Intrinsic Rotation in Flux-Driven ITG Turbulence

Global, heat flux-driven ITG gyrokinetic simulations which manifest the formation of macroscopic, mean toroidal flow profiles with peak thermal Mach number 0.05, are reported. Both a particle-in-cell (XGC1p) and a semi-Lagrangian (GYSELA) approach are utilized without a priori assumptions of scale-separation between turbulence and mean fields. Flux-driven ITG simulations with different edge flow boundary conditions show in both approaches the development of net unidirectional intrinsic rotation in the co-current direction. Intrinsic torque is shown to scale approximately linearly with the inverse scale length of the ion temperature gradient. External momentum input is shown to effectively cancel the intrinsic rotation profile, thus confirming the existence of a local residual stress and intrinsic torque. Fluctuation intensity, intrinsic torque and mean flow are demonstrated to develop inwards from the boundary. The measured correlations between residual stress and two fluctuation spectrum symmetry breakers, namely E x B shear and intensity gradient, are similar. Avalanches of (positive) heat flux, which propagate either outwards or inwards, are correlated with avalanches of (negative) parallel momentum flux, so that outward transport of heat and inward transport of parallel momentum are correlated and mediated by avalanches. The probability distribution functions of the outward heat flux and the inward momentum flux show strong structural similarity

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Techno-Economic Analysis for a Potential Geothermal District Heating System in Tuttle, Oklahoma

Geothermal deep direct use (DDU) has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for a large-scale and variable heat demand. The National Renewable Energy Laboratory (NREL) and University of Oklahoma evaluated the feasibility of a geothermal district heating (GDH) and cooling system in two schools and 250 houses by utilizing existing oil and gas (O&G) wells in Tuttle, Oklahoma. Heating and cooling demand in the two schools and a typical single-family house were modeled using EnergyPlus building energy simulation software. The modeling results indicated that annual heating demand in two schools and 250 houses is approximately 2.61 GWhth, and cooling demand in the two schools is approximately 2.65 GWhth. In this scope, the techno-economic analysis (TEA) was conducted using the GEOPHIRES tool combined with the TOUGH2 reservoir simulator. The reservoir performance, including geothermal heat production capacity, was modeled by the reservoir simulator TOUGH2. Then, levelized cost of heat (LCOH) was calculated using GEOPHIRES version 3.0, which includes new features such as hourly heat load optimization and peak performance evaluation. Geothermal reservoir temperature was estimated as 90.5 degrees C at a total depth of 3.3 km by the regional average temperature gradient of 22.8 degrees C/km and validated by cation geothermometer calculations. Four production scenarios with two different well configurations and two different heat load profiles have been developed for well flow rates ranging between 3.1 kg/s and 9.3 kg/s. The LCOH of the district heating and cooling system was calculated between $95 and $210/MWh ($28/MMBtu to $62/MMBtu) for four different production scenarios. Typical natural gas prices for residential customers in Oklahoma have ranged from 9 to 19 $/MMBtu over the past decade, which indicates a challenge for deployment of such a GDH system.

deep direct-use↗

Measurement of heat conduction through stacked screens

This paper describes the experimental apparatus for the measurement of heat conduction through stacked screens as well as some experimental results taken with the apparatus. Screens are stacked in a fiberglass-epoxy cylinder, which is 24.4 mm in diameter and 55 mm in length. The cold end of the stacked screens is cooled by a Gifford-McMahon (GM) cryocooler at cryogenic temperature, and the hot end is maintained at room temperature. Heat conduction through the screens is determined from the temperature gradient in a calibrated heat flow sensor mounted between the cold end of the stacked screens and the GM cryocooler. The samples used for these experiments consisted of 400-mesh stainless steel screens, 400-mesh phosphor bronze screens, and two different porosities of 325-mesh stainless steel screens. The wire diameter of the 400-mesh stainless steel and phosphor bronze screens was 25.4 micrometers and the 325-mesh stainless steel screen wire diameters were 22.9 micrometers and 27.9 micrometers. Standard porosity values were used for the experimental data with additional porosity values used on selected experiments. The experimental results showed that the helium gas between each screen enhanced the heat conduction through the stacked screens by several orders of magnitude compared to that in vacuum. The conduction degradation factor is the ratio of actual heat conduction to the heat conduction where the regenerator material is assumed to be a solid rod of the same cross sectional area as the metal fraction of the screen. This factor was about 0.1 for the stainless steel and 0.022 for the phosphor bronze, and almost constant for the temperature range of 40 to 80 K at the cold end.

NASA Discipline Environmental Health↗

Pyroelectric field assisted ion migration induced by ultraviolet laser irradiation and its impact on ferroelectric domain inversion in lithium niobate crystals

The impact of UV laser irradiation on the distribution of lithium ions in ferroelectric lithium niobate single crystals has been numerically modelled. Strongly absorbed UV radiation at wavelengths of 244–305 nm produces steep temperature gradients which cause lithium ions to migrate and result in a local variation of the lithium concentration. In addition to the diffusion, here the pyroelectric effect is also taken into account which predicts a complex distribution of lithium concentration along the c-axis of the crystal: two separated lithium deficient regions on the surface and in depth. The modelling on the local lithium concentration and the subsequent variation of the coercive field are used to explain experimental results on the domain inversion of such UV treated lithium niobate crystals.

36 MATERIALS SCIENCE↗

Techno-Economic Analysis for a Potential Geothermal District Heating System in Tuttle, Oklahoma: Preprint

Geothermal deep direct use (DDU) has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for a large-scale and variable heat demand. The National Renewable Energy Laboratory (NREL) and University of Oklahoma evaluated the feasibility of a geothermal district heating (GDH) and cooling system in two schools and 250 houses by utilizing existing oil and gas (O&G) wells in Tuttle, Oklahoma. Heating and cooling demand in the two schools and a typical single-family house were modeled using EnergyPlus building energy simulation software. The modeling results indicated that annual heating demand in two schools and 250 houses is approximately 2.61 GWhth, and cooling demand in the two schools is approximately 2.65 GWhth. In this scope, the techno-economic analysis (TEA) was conducted using the GEOPHIRES tool combined with the TOUGH2 reservoir simulator. The reservoir performance, including geothermal heat production capacity, was modeled by the reservoir simulator TOUGH2. Then, levelized cost of heat (LCOH) was calculated using GEOPHIRES version 3.0, which includes new features such as hourly heat load optimization and peak performance evaluation. Geothermal reservoir temperature was estimated as 90.5 degrees C at a total depth of 3.3 km by the regional average temperature gradient of 22.8 degrees C/km and validated by cation geothermometer calculations. Four production scenarios with two different well configurations and two different heat load profiles have been developed for well flow rates ranging between 3.1 kg/s and 9.3 kg/s. The LCOH of the district heating and cooling system was calculated between $95 and $210/MWh ($28/MMBtu to $62/MMBtu) for four different production scenarios. Typical natural gas prices for residential customers in Oklahoma have ranged from 9 to 19 $/MMBtu over the past decade, which indicates a challenge for deployment of such a GDH system.

deep direct-use↗

O 2 based resonantly ionized photoemission thermometry analysis of supersonic flows

Characterization of the thermal gradients within supersonic and hypersonic flows is essential for understanding transition, turbulence, and aerodynamic heating. Developments in novel, impactful non-intrusive techniques are key for enabling flow characterizations of sufficient detail that provide experimental validation datasets for computational simulations. In this work, Resonantly Ionized Photoemission Thermometry (RIPT) signals are directly imaged using an ICCD camera to realize the techniques 1D measurement capability for the first time. The direct imaging scheme presented for oxygen-based RIPT (O 2 RIPT) uses the previously established calibration data to direct excite various resonant rotational peaks within the S-branch of the C 3 Π, ( v = 2) ← X 3 Σ( v ′ = 0) absorption band of O 2 . The efficient ionization of O 2 liberates electrons that induce electron avalanche ionization of local N 2 molecules generating N 2 + , which primarily deexcites via photoemissions of the first negative band of N 2 + ( B 2 Σ <#comment/> u + − <#comment/> X 2 Σ <#comment/> g + ) . When sufficient lasing energy is used, the ionization region and subsequent photoemission signal is achieved along a 1D line thus, if directly imaged can allow for gas temperature assignments along said line; demonstrated here of up to five centimeters in length. The temperature gradients present within the ensuing shock train of a supersonic under expanded free jet serves as a basis of characterization for this new RIPT imaging scheme. The O 2 RIPT results are extensively compared and validated against well-known and established techniques (i.e., CARS and CFD). The direct imaging capability fully realizes the technique’s fundamental potential and is expected to be the standard of implementation going forward. The direct imaging capability can play instrumental roles in future scientific studies that rely upon acute characterization of thermal gradients within a medium that cannot be easily resolved by a point. Furthermore, the removal of the spectrometer greatly reduces the cost, complexity, and optical alignment associated with prior RIPT measurements.

42 ENGINEERING↗

A Qualitative Investigation of Selected Infrared Flow Visualization Image Processing Techniques

The deployment and integration of high-sensitivity infrared cameras in a transonic wind tunnel testenvironment has resulted in a unique capability to image aerodynamic phenomena in real-time. Multi-camera infrared flow visualization data systems are now routinely utilized at the NASA Ames Unitary Plan Wind Tunnel. The small flow-induced temperature gradients on the surface of the wind tunnel test article coupled with the high bit-depth of the infrared camera sensor makes the processing of the image data critically important. An image processing routine must enhance features of interest with minimal artifacts. Additionally, the production wind tunnel test environment demands that these processed images are made available in a real-time, automatic fashion. Therefore, any image processing routine must be computationally economical and enhance the image data with minimal input from a human operator. The following seeks to qualitatively explore selected image processing techniques by assessing their effectiveness to resolve flow features on a wind tunnel test article. A multi-scale contrast enhancement technique is discussed as well as a new implementation of a multi-scale, non-interpolated adaptive histogram equalization. Finally, a novel method is introduced that demonstrates the ability to resolve flow features imaged on bare-steel test articles possessing low emissivity.This method makes use of dynamic mode decomposition and discrete-time filtering to separate the background reflections that dominate low emissivity surfaces from the aerodynamic driven surface temperature gradients.This process will be shown to resolve the onset of boundary layer transition on a bare metal wing as well as identify and resolve hidden features in the image data. While the implementation of this technique is very preliminary it demonstrates the potential to extend the application of infrared flow-visualization within the wind tunnel test environment.

Image Processeing↗

First-Principles Calculations of Thermoelectric Transport Properties of Quaternary and Ternary Bulk Chalcogenide Crystals

Chalcogenide crystals have a wide range of applications, especially as thermoelectric materials for energy conversion. Thermoelectric materials can be used to generate an electric current from a temperature gradient based on the Seebeck effect and based on the Peltier effect, and they can be used in cooling applications. Using first-principles calculations and semiclassical Boltzmann theory, we have computed the Seebeck coefficient, electrical conductivity, electronic thermal conductivity, power factor, and figure of merit of 30 chalcogenide crystals. A Quantum Espresso package is used to calculate the electronic properties and locate the Fermi level. The transport properties are then calculated using the BoltzTraP code. The 30 crystals are divided into two groups. The first group has four crystals with quaternary composition (A 2 BCQ 4 ) (A = Tl; B = Cd, Hg; C = Si, Ge, Sn; Q = S, Se, Te). The second group contains 26 crystals with the ternary composition (A’B’Q 2 ) (A’ = Ag, Cu, Au, Na; B’ = B, Al, Ga, In; Q = S, Se, Te). Among these 30 chalcogenide crystals, the results for 11 crystals: Tl 2 CdGeSe 4 , Tl 2 CdSnSe 4 , Tl 2 HgSiSe 4 , Tl 2 HgSnS 4 , AuBSe 2 , AuBTe 2 , AuAlTe 2 , AuGaTe 2 , AuInTe 2 , AgAlSe 2 , and AgAlTe 2 are revealed for the first time. In addition, temperature-dependent transport properties of pure and doped AgSbSe 2 and AgSbTe 2 crystals with dopant compositions of AgSb 0.94 Cd 0.06 Te 2 and AgSbTe 1.85 Se 0.15 were explored. These results provide an excellent database for bulk chalcogenides crucial for a wide range of potential applications in renewable energy fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation Of Microtearing Turbulence In NSTX

Thermal energy confinement times in NSTX dimensionless parameter scans increase with decreasing collisionality. While ion thermal transport is neoclassical, the source of anomalous electron thermal transport in these discharges remains unclear, leading to considerable uncertainty when extrapolating to future ST devices at much lower collisionality. Linear gyrokinetic simulations find microtearing modes to be unstable in high collisionality discharges. First non-linear gyrokinetic simulations of microtearing turbulence in NSTX show they can yield experimental levels of transport. Magnetic flutter is responsible for almost all the transport (~98%), perturbed field line trajectories are globally stochastic, and a test particle stochastic transport model agrees to within 25% of the simulated transport. Most significantly, microtearing transport is predicted to increase with electron collisionality, consistent with the observed NSTX confinement scaling. While this suggests microtearing modes may be the source of electron thermal transport, the predictions are also very sensitive to electron temperature gradient, indicating the scaling of the instability threshold is important. In addition, microtearing turbulence is susceptible to suppression via sheared E-B flows as experimental values of E-B shear (comparable to the linear growth rates) dramatically reduce the transport below experimental values. Refinements in numerical resolution and physics model assumptions are expected to minimize the apparent discrepancy. In cases where the predicted transport is strong, calculations suggest that a proposed polarimetry diagnostic may be sensitive to the magnetic perturbations associated with the unique structure of microtearing turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing the Response of Tropical Deep Convection to Aerosol Perturbations Using Idealized Cloud-Resolving Simulations with Parameterized Large-Scale Dynamics

Abstract A framework is introduced to investigate the indirect effect of aerosol loading on tropical deep convection using three-dimensional limited-domain idealized cloud-system-resolving model simulations coupled with large-scale dynamics over fixed sea surface temperature. The large-scale circulation is parameterized using the spectral weak temperature gradient (WTG) approximation that utilizes the dominant balance between adiabatic cooling and diabatic heating in the tropics. The aerosol loading effect is examined by varying the number of cloud condensation nuclei (CCN) available to form cloud droplets in the two-moment bulk microphysics scheme over a wide range of environments from 30 to 5000 cm−3. The radiative heating is held at a constant prescribed rate in order to isolate the microphysical effects. Analyses are performed over the period after equilibrium is achieved between convection and the large-scale environment. Mean precipitation is found to decrease modestly and monotonically when the aerosol number concentration increases as convection gets weaker, despite the increase in cloud liquid water in the warm-rain region and ice crystals aloft. This reduction is traced down to the reduction in surface enthalpy fluxes as an energy source to the atmospheric column induced by the coupling of the large-scale motion, though the gross moist stability remains constant. Increasing CCN concentration leads to 1) a cooler free troposphere because of a reduction in the diabatic heating and 2) a warmer boundary layer because of suppressed evaporative cooling. This dipole temperature structure is associated with anomalously descending large-scale vertical motion above the boundary layer and ascending motion at lower levels. Sensitivity tests suggest that changes in convection and mean precipitation are unlikely to be caused by the impact of aerosols on cloud droplets and microphysical properties but rather by accounting for the feedback from convective adjustment with the large-scale dynamics. Furthermore, a simple scaling argument is derived based on the vertically integrated moist static energy budget, which enables estimation of changes in precipitation given known changes in surfaces enthalpy fluxes and the constant gross moist stability. The impact on cloud hydrometeors and microphysical properties is also examined, and it is consistent with the macrophysical picture.

54 ENVIRONMENTAL SCIENCES↗

Setting the H-mode pedestal structure: variations of particle source location using gas puff and pellet fueling

Experiments in DIII-D show that a particle source location inside the top of the H-mode pedestal (pellets) maintains a higher pedestal pressure than an edge source (gas fueling) through a widening of the electron temperature pedestal with reduction of the temperature gradient. The effect of these two fueling schemes on the H-mode pedestal structure was examined in DIII-D by comparing controlled pellet-fueled and gas-fueled discharges across a fueling scan up to 40 torr 1 s –1 . High resolution electron profiles reveal that gas fueling lowers the pedestal pressure as the density profile shifts radially outwards and the separatrix density increases, while pellet fueling maintains a constant pedestal pressure. The neutral source locations from pellets and gas are determined with the PELLET and UEDGE codes, respectively, and quantify the particle source localization. Pellets provide significant ionization inside the pedestal top while gas puffing localizes ionization in the scrape-off-layer and pedestal foot, broadly consistent with the density profile structure influenced by the source. ELMs are observed to increase in frequency and reduce impurity content as fueling is increased. Stability analysis with ELITE shows that both conditions are near the type-I ELM corner of the peeling–ballooning stability diagram, which is altered significantly by the introduction of pellets. Since transport mechanisms are not observed to change substantially with particle source location, wider pedestals allow the pellet-fueled discharges to retain higher pedestal temperatures at similar pedestal densities. EPED1 is tested to capture the pedestal pressure, under-predicting the height with pellets and over-predicting the height with gas by ~15%. Overall, these results have important implications for future reactors where pellet fueling will be the primary particle source due to an opaque scrape-off-layer by showing that the neutral source location plays a role in setting the structure of the H-mode pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimation of the Surface Fluxes for Heat and Momentum in Unstable Conditions with Machine Learning and Similarity Approaches for the LAFE Data Set

Abstract Measurements of three flux towers operated during the land atmosphere feedback experiment (LAFE) are used to investigate relationships between surface fluxes and variables of the land–atmosphere system. We study these relations by means of two machine learning (ML) techniques: multilayer perceptrons (MLP) and extreme gradient boosting (XGB). We compare their flux derivation performance with Monin–Obukhov similarity theory (MOST) and a similarity relationship using the bulk Richardson number (BRN). The ML approaches outperform MOST and BRN. Best agreement with the observations is achieved for the friction velocity. For the sensible heat flux and even more so for the latent heat flux, MOST and BRN deviate from the observations while MLP and XGB yield more accurate predictions. Using MOST and BRN for latent heat flux, the root mean square errors (RMSE) are 107 Wm $$^{-2}$$ - 2 and 121 Wm $$^{-2}$$ - 2 , respectively, as well as the intercepts of the regression lines are $$\approx 110$$ ≈ 110 Wm $$^{-2}$$ - 2 . For the ML methods, the RMSEs reduce to 31 Wm $$^{-2}$$ - 2 for MLP and 33 Wm $$^{-2}$$ - 2 for XGB as well as the intercepts to just 4 Wm $$^{-2}$$ - 2 for MLP and $$-1$$ - 1 Wm $$^{-2}$$ - 2 for XGB with slopes of the regression lines close to 1, respectively. These results indicate significant deficiencies of MOST and BRN, particularly for the derivation of the latent heat flux. In fact, in contrast to the established theories, feature importance weighting demonstrates that the ML methods base their improved derivations on net radiation, the incoming and outgoing shortwave radiations, the air temperature gradient, and the available water contents, but not on the water vapor gradient. The results imply that further studies of surface fluxes and other turbulent variables with ML techniques provide great promise for deriving advanced flux parameterizations and their implementation in land–atmosphere system models.

54 ENVIRONMENTAL SCIENCES↗

The Response of the Quasi-Biennial Oscillation to Increased CO and its Modulation by Composition Feedbacks

The Quasi-Biennial Oscillation (QBO) impacts the large-scale circulation by altering lower stratospheric meridional temperature gradients which alter the propagation of upward propagating planetary waves. Proposed teleconnections of the QBO include impacts on the stratospheric polar vortices, extratropical surface winter climate, and the Madden-Julian Oscillation. However, long-term projections of the QBO remain highly uncertain. While recent multi-model investigations (Richter et al. 2019, Butchart et al. 2020) show that the amplitude of the QBO weakens robustly among models, changes in QBO period and stability remain highly uncertain. Here we examine the QBO response to increased greenhouse gases using the NASA Goddard Institute for Space Studies Middle Atmosphere Model E2.2 (Rind et al., 2020; Orbe et al. 2020). Compared to lower vertical resolution versions of ModelE, E2.2 also has a higher model top (0.002 hPa) and employs additional interactive non-orographic gravity wave drag sources from convection and shear, which produce a sufficiently realistic QBO, thus rendering it suitable for use in climate change studies. Overall, we find that both the QBO period and amplitude decrease in response to increased CO , the former related to increased lower stratospheric momentum fluxes (associated with convection) and the latter associated partly with a stronger residual mean circulation. Experiments integrated using fully interactive chemistry also reveal that ozone feedbacks significantly impact the magnitude of the QBO amplitude response. Finally, integrations using fixed (pre-industrial) sea surface temperatures (SST) show that the QBO amplitude responds differently to rapid adjustments versus SST feedbacks, as compared to the QBO period.

Quasi-Biennial Oscillation↗

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fiber Optic Sensing Technologies for Battery Management Systems and Energy Storage Applications

Applications of fiber optic sensors to battery monitoring have been increasing due to the growing need of enhanced battery management systems with accurate state estimations. The goal of this review is to discuss the advancements enabling the practical implementation of battery internal parameter measurements including local temperature, strain, pressure, and refractive index for general operation, as well as the external measurements such as temperature gradients and vent gas sensing for thermal runaway imminent detection. A reasonable matching is discussed between fiber optic sensors of different range capabilities with battery systems of three levels of scales, namely electric vehicle and heavy-duty electric truck battery packs, and grid-scale battery systems. The advantages of fiber optic sensors over electrical sensors are discussed, while electrochemical stability issues of fiber-implanted batteries are critically assessed. This review also includes the estimated sensing system costs for typical fiber optic sensors and identifies the high interrogation cost as one of the limitations in their practical deployment into batteries. Finally, future perspectives are considered in the implementation of fiber optics into high-value battery applications such as grid-scale energy storage fault detection and prediction systems.

25 ENERGY STORAGE↗

Development and application of processing and process control for nano-composite materials for lithium ion batteries (CRADA Final Report)

Oak Ridge National Laboratory and A123 Systems, Inc. collaborated on this project to develop a better understanding, quality control procedures, and safety testing for A123 System’s nanocomposite separator (NCS) technology which is a cell based patented technology and separator. NCS demonstrated excellent performance. x3450 prismatic cells were shown to survive >8000 cycles (1C/2C rate) at room temperature with greater than 80% capacity retention with only NCS present as an alternative to conventional polyolefin. However, for a successful commercialization, the coating conditions required to provide consistent and reliable product had not been optimized and QC techniques for being able to remove defective material before incorporation into a cell had not been developed. The work outlined in this report addresses these latter two points. First, experiments were conducted to understand temperature profiles during the different drying stages of the NCS coating when applied to both anode and cathode. One of the more interesting discoveries of this study was the observation of the large temperature decrease experienced by the wet coating between the end of the infrared (IR) drying stage and the beginning of the exposure to the convection drying oven. This is not a desirable situation as the temperature gradient could have a deleterious effect on coating quality. Based on this and other experimental data a radiative transfer model was developed for IR heating that also included a mass transfer module for drying. This will prove invaluable for battery coating optimization especially where IR drying is being employed. A stress model was also developed that predicts that under certain drying conditions tensile stresses are formed in the coating which could lead to cracking that is sometimes observed after drying is complete. Prediction of under what conditions these stresses form is vital to improving coating quality. In addition to understanding the drying process other parameters such as slurry quality and equipment optimization were examined. Removal of particles and gels by filtering, control of viscosity by %solids and mixing adjustments, removal of trapped gas in the slurry and modification of coater speed and slot die gap were all found to be important for producing uniform and flaw-free coatings. Second, an in-line Hi-Pot testing method has been developed specifically for NCS that will enable detection of coating flaws that could lead to soft or hard electrical shorts within the cell. In this way flawed material can be rejected before incorporation into the cell thus greatly reducing the amount of scrap that is generated. Improved battery safety is an extremely important benefit of NCS. Evaluation of battery safety is usually accomplished by conducting a variety of tests including nail penetration, hot box, over charge, etc. For these tests entire batteries must be built but the resultant temperature and voltage responses reveal little about the breakdown mechanism. In this report is described a pinch test which is used to evaluate NCS quality at various stages including coated anode and cathode as well as assembled cell. Coupled with post-microscopic examination of the damaged ‘pinch point’ test data can assist in the coating optimization from an improved end-use standpoint. As a result of this work two invention disclosures, one for optimizing drying methodology and the other for an in-line system for flaw detection, have been filed. In addition, 2 papers are being written for submission to peer-reviewed journals.

25 ENERGY STORAGE↗

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗

Importance of gas heating in capacitively coupled radiofrequency plasma-assisted synthesis of carbon nanomaterials

In pursuit of diamond nanoparticles, a capacitively-coupled radio frequency flow-through plasma reactor was operated with methane-argon gas mixtures. Signatures of the final product obtained microscopically and spectroscopically indicated that the product was an amorphous form of graphite. This result was consistent irrespective of combinations of the macroscopic reactor settings. To explain the observed synthesis output, measurements of C 2 and gas properties were carried out by laser-induced fluorescence and optical emission spectroscopy. Strikingly, the results indicated a strong gas temperature gradient of 100 K per mm from the center of the reactor to the wall. Based on additional plasma imaging, a model of hot constricted region (filamentation region) was then formulated. It illustrated that, while the hot constricted region was present, the bulk of the gas was not hot enough to facilitate diamond sp 3 formation: characterized by much lower reaction rates, when compared to sp 2 , sp 3 formation kinetics are expected to become exponentially slow. This result was further confirmed by experiments under identical conditions but with a H 2 /CH 4 mixture, where no output material was detected: if graphitic sp 2 formation was expected as the main output material from the methane feedstock, atomic hydrogen would then be expected to etch it away in situ, such that the net production of that sp 2 -hybridized solid material is nearly a zero. Finally, the crucial importance of gas heating was corroborated by replacing RF with microwave source whereby facile sp 3 production was attained with H 2 /CH 4 gas mixture.

36 MATERIALS SCIENCE↗