Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global refinement”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

LES of Vaporizing Gasoline Sprays Considering Multi-Injection Averaging and Grid-Convergent Mesh Resolution

A state-of-the-art spray modeling methodology, recently presented by Senecal et al. [ , , ], is applied to Large Eddy Simulations (LES) of vaporizing gasoline sprays. Simulations of non-combusting Spray G (gasoline fuel) from the Engine Combustion Network are performed. Adaptive mesh refinement (AMR) with cell sizes from 0.09 mm to 0.5 mm are utilized to further demonstrate grid convergence of the dynamic structure LES model for the gasoline sprays. Grid settings are recommended to optimize the accuracy/runtime tradeoff for LES-based spray simulations at different injection pressure conditions typically encountered in gasoline direct injection (GDI) applications. The influence of LES sub-grid scale (SGS) models is explored by comparing the results from dynamic structure and Smagorinsky based models against simulations without any SGS model. Twenty different realizations are simulated by changing the random number seed used in the spray sub-models. It is shown that for global quantities such as spray penetration, comparing a single LES simulation to experimental data is reasonable. Through a detailed analysis using the relevance index (RI) criteria, recommendations are made regarding the minimum number of LES realizations required for accurate prediction of the gasoline sprays.

Som, Sibendu↗
DOE OSTI2018

An Hp-Adaptive Minimum Action Method Based on a Posteriori Error Estimate

In this work, we develop an hp-adaptive minimum action method (MAM). MAM plays an important role to minimize the Freidlin-Wentzell action functional, which is the central object of the Freidlin-Wentzell theory of large deviations for transitions induced by small random perturbations. Because of the demanding computation cost, especially in spatially extended systems, numerical efficiency is a critical issue for MAM. Difficulties come from discretizations of both time and space. One hurdle for the application of MAM to large scale systems is the global reparametrization in time direction, which is required by most versions of MAM. We recently introduced a new version of MAM in [25], called tMAM, where we used some simple heuristic criteria to demonstrate that tMAM can be effectively coupled with h-adaptivity. The target of this paper is to integrate hp-adaptivity into tMAM using a posterior error estimation techniques. More specifically, we use the arc length constraint given by geometric minimum action method (gMAM) to define an indicator of the effect of linear time scaling in tMAM, and the derivative recovery technique to construct an error indicator and a regularity indicator for hp-refinement. Numerical results are presented.

Wan, Xiaoliang↗
DOE OSTI2015

Global critical materials markets: An agent-based modeling approach

As part of efforts to position the United States as a leader in clean energy technology production, the U. S. Department of Energy (DOE) issued two Critical Materials Strategy reports, which assessed 16 materials on the basis of their importance to clean energy development and their supply risk ( U.S. Department of Energy (DOE), 2010 and DOE, 2011). To understand the implications for clean energy of disruptions in supplies of critical materials, it is important to understand supply chain dynamics from mining to final product production. As a case study of critical material supply chains, we focus on the supply of two rare earth metals, neodymium (Nd) and dysprosium (Dy), for permanent magnets used in wind turbines, electric vehicles and other applications. We introduce GCMat, a dynamic agent-based model that includes interacting agents at five supply chain stages consisting of mining, metal refining, magnet production, final product production and demand. Agents throughout the supply chain make pricing, production and inventory management decisions. Deposit developers choose which deposits to develop based on market conditions and detailed data on 57 rare earth deposits. Wind turbine and electric vehicle producers choose from a set of possible production technologies that require different amounts of rare earths. We ran the model under a baseline scenario and four alternative scenarios with different demand and production technology inputs. Model results from 2010 to 2013 fit well with historical data. Projections through 2025 show a number of possible future price, demand, and supply trajectories. For each scenario, we highlight reasons for turning points under market conditions, for differences between Nd and Dy markets, and for differences between scenarios. Because GCMat can model causal dynamics and provide fine-grain representation of agents and their decisions, it provides explanations for turning points under market conditions that are not otherwise available from other modeling approaches. Our baseline projections show very different behaviors for Nd and Dy prices. Nd prices continue to drop and remain low even at the end of our simulation period as new capacity comes online and leads to a market in which production capacity outpaces demand. Dy price movements, on the other hand, change directions several times with several key turning points related to inventory behaviors of particular agents in the supply chain and asymmetric supply and demand trends. Scenario analyses show the impact of stronger demand growth for rare earths, and in particular finds that Nd price impacts are significantly delayed as compared to Dy. This is explained by the substantial excess production capacity for Nd in the early simulation years that keeps prices down. Scenarios that explore the impact of reducing the Dy content of magnets show the intricate interdependencies of these two markets as price trends for both rare earths reverse directions – reducing the Dy content of magnets reduces Dy demand, which drives down Dy prices and translates into lower magnet prices. This in turn raises the demand for magnets and therefore the demand for Nd and eventually drives up the Nd price.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An ecologically-constrained deep learning model for tropical leaf phenology monitoring using PlanetScope satellites

In tropical forests, leaf phenology signals leaf-on/off status and exhibits considerable variability across scales from a single tree-crown to the entire forest ecosystem. Such phenology signals importantly regulate large-scale biogeochemical cycles and regional climate. PlanetScope CubeSats data with a 3-m resolution and near-daily global coverage provide an unprecedented opportunity to monitor both fine- and ecosystem-scale phenology variability along large environmental gradients. However, a scalable method that accurately characterizes leaf phenology from PlanetScope with biophysically meaningful metrics remains lacking. Here we developed an index-guided, ecologically constrained autoencoder (IG-ECAE) method to automatically derive a deciduousness metric (i.e. percentage of tree canopies with leaf-off status within an image pixel) from PlanetScope. The IG-ECAE first estimated the reflectance spectra of leafy/leafless canopies based on their spectral indices characteristics, then used the derived reflectance spectra to guide an autoencoder deep learning method with additional ecological constraints to refine the reflectance spectra, and finally used linear spectral unmixing to estimate the relative abundance of leafless canopies (or deciduousness) per PlanetScope image pixel. We tested the IG-ECAE method at 16 tropical forest sites spanning multiple continents and a large precipitation gradient (1470-2819 mm year -1 ). Among these sites, we evaluated the PlanetScope-derived deciduousness with corresponding measures derived from WorldView-2 (n=9 sites) and local phenocams (n=9 sites). Our results show that PlanetScope-derived deciduousness agrees: 1) with that derived from WorldView-2 at the patch level (90m×90m) with r 2 =0.89 across all sites; and 2) with that derived from phenocams to quantify ecosystem-scale seasonality with r 2 ranging from 0.62-0.96. These results demonstrate the effectiveness and scalability of IG-ECAE in characterizing the wide variability in deciduousness across scales from pixels to forest ecosystems, and from a single date to the full annual cycle, indicating the potential for using high-resolution satellites to track the large-scale phenological patterns and response of tropical forests to climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗
DOE OSTI2022

PDBspheres: a method for finding 3D similarities in local regions in proteins

Abstract We present a structure-based method for finding and evaluating structural similarities in protein regions relevant to ligand binding. PDBspheres comprises an exhaustive library of protein structure regions (‘spheres’) adjacent to complexed ligands derived from the Protein Data Bank (PDB), along with methods to find and evaluate structural matches between a protein of interest and spheres in the library. PDBspheres uses the LGA (Local–Global Alignment) structure alignment algorithm as the main engine for detecting structural similarities between the protein of interest and template spheres from the library, which currently contains >2 million spheres. To assess confidence in structural matches, an all-atom-based similarity metric takes side chain placement into account. Here, we describe the PDBspheres method, demonstrate its ability to detect and characterize binding sites in protein structures, show how PDBspheres—a strictly structure-based method—performs on a curated dataset of 2528 ligand-bound and ligand-free crystal structures, and use PDBspheres to cluster pockets and assess structural similarities among protein binding sites of 4876 structures in the ‘refined set’ of the PDBbind 2019 dataset.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Science Plan for the Deployment of the Third ARM Mobile Facility to the Southeastern United States at the Bankhead National Forest, Alabama (AMF3 BNF)

In 2018, the U.S. Department of Energy (DOE) held a workshop for the Atmospheric Radiation Measurement (ARM) (Mather and Voyles 2013) user facility to discuss critical climate challenges and locations where key ARM Mobile Facility (AMF) observational assets could impact Earth system modeling (ESM). As an outcome, the southeast United States (SE U.S.) was identified as a high-priority region to target climate-process studies that promote a deeper understanding of the climate system and bolster ARM interactions with the community to drive ESM advancement. The DOE ARM user facility is a globally recognized leader in deploying and operating strategically located observation sites around the world for studying the properties of aerosols and clouds and their interaction with radiation, precipitation, and the Earth’s surface. In partnering with the DOE Atmospheric System Research (ASR) program, ARM solicited a multi-agency Site Science Team approach to provide input and close interaction with ARM management towards a successful SE U.S. deployment of the ARM third Mobile Facility (AMF3) (Miller et al. 2016). These efforts included identifying key locations, science drivers and instruments, and measurement strategies to address the wider climate-process needs and ESM improvement. Community input served a vital role in establishing, refining, and informing the relevant drivers and decisions regarding this AMF3 deployment. The team has identified Northern Alabama (N. AL) as regionally representative to unlock the key opportunities that will improve our understanding and model representation of aerosol, cloud, and land-surface processes and their couplings in the SE U.S. A defining aspect of the AMF3 deployment is its commitment to long-term (anticipated five-year) observations to mitigate potential seasonal-to-annual variability that often limits appropriate attribution of phenomena to local or larger-scale processes. The proposed location may leverage nearby surface networks and multi-agency and partner assets to enrich this multi-year deployment. One motivation is to understand the role of spatiotemporal variability (thermodynamic, land-surface) across aspects of the climate system, with our AMF3 team anticipating future demands on characterizing the relationships between local-to-regional cloud development and surface processes across a diverse patchwork of natural, managed, and urban landscapes as found throughout the N. AL regions. The main site targets an intact, representative, forested region – the Bankhead National Forest (BNF) – underscoring further team commitment to regionally important land-atmosphere two-way interactive studies “from the canopy to the clouds”, with enhanced tower instrumentation augmenting traditional ARM capabilities adjacent to this site. Multiple supplemental sites will also be distributed across this region, prioritizing added needs for biodiversity. Anticipated high-priority cloud science themes will target N. AL as a regional SE U.S. hotbed for high-impact weather, convective cloud onset, and shallow-to-deep cloud transitioning. Anticipated aerosol drivers will focus on chemical processes that control the evolution of organic aerosol, the seasonality and spatial distribution of water vapor and particle-phase water, and its role on aerosol optical properties. Anticipated land-atmosphere drivers consider the two-way feedbacks between surface influence on aerosols, clouds, and precipitation properties and the associated radiative impacts on plant physiology and canopy-scale fluxes. Emphasis will include the study of the impact of surface processes on aerosols via precursor emission, and on clouds via moisture flux and thermal development.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

NASTRAN finite element idealization study

The investigation of the effects of variations of mesh refinement and mesh pattern were conducted using a basic rectangular mesh pattern. When employing the constant strain TRMEM element, the basic rectangular pattern was subdivided into triangles. This subdivision employs two different triangular patterns and allows results to be obtained which demonstrate the effect of modelling bias. Errors in tip deflection, direct stress, and shearing stress as a function of mesh size and element aspect ratio were obtained as well as mid-span stress distributions. All problems were solved on an IBM 360/95 computer using MacNeal-Schwendler Version MSC-38 Rigid Format-1 of the NASTRAN computer program. While NASTRAN uses double precision arithmetic for the solution of the global equations for displacements, subsequent computations to obtain element stresses are carried out in single precision. This suggests that some improvement in stress recovery might be expected when using the higher precision DCD machines.

Case, W. R.↗
NASA NTRS1987

Guidance simulation and test support for differential GPS flight experiment

Three separate tasks which supported the test preparation, test operations, and post test analysis of the NASA Ames flight test evaluation of the differential Global Positioning System (GPS) are presented. Task 1 consisted of a navigation filter design, coding, and testing to optimally make use of GPS in a differential mode. The filter can be configured to accept inputs from external censors such as an accelerometer and a barometric or radar altimeter. The filter runs in real time onboard a NASA helicopter. It processes raw pseudo and delta range measurements from a single channel sequential GPS receiver. The Kalman filter software interfaces are described in detail, followed by a description of the filter algorithm, including the basic propagation and measurement update equations. The performance during flight tests is reviewed and discussed. Task 2 describes a refinement performed on the lateral and vertical steering algorithms developed on a previous contract. The refinements include modification of the internal logic to allow more diverse inflight initialization procedures, further data smoothing and compensation for system induced time delays. Task 3 describes the TAU Corp participation in the analysis of the real time Kalman navigation filter. The performance was compared to that of the Z-set filter in flight and to the laser tracker position data during post test analysis. This analysis allowed a more optimum selection of the parameters of the filter.

Geier, G. J.↗
NASA NTRS1999

Determination of Thunderstorm Anvil Ice Contents and Other Cloud Properties from Satellite Observations of Lightning

The lightning frequency model developed by Baker, Christian and Latham (1995) has been refined and extended, in an effort to provide a more realistic framework from which to examine computationally the relationships that might exist between lightning frequency f (which is now being routinely measured from satellites, using NASA/MSFC devices) and a variety of cloud physical parameters, including precipitation rate, updraught speed and non-precipitating ice content. Model results indicate the existence of a simple relationship between lightning frequency f and the upward flux of ice crystals into the thunderstorm anvil. It follows that, for a particular situation, one can assign a specific mass of non-precipitating ice to an individual lightning stroke. Therefore it may prove possible - using satellite measurements of global lightning - to estimate the atmospheric loading of ice crystals in thunderstorm anvils: a parameter of climatological importance. Early results from this work are presented, together with further studies of the relationships between f and other thundercloud parameters.

Blyth, Alan M.↗
DOE OSTI2023

Science Plan for the Deployment of the Third ARM Mobile Facility to the Southeastern United States at the Bankhead National Forest, Alabama (AMF3 BNF)

In 2018, the U.S. Department of Energy (DOE) held a workshop for the Atmospheric Radiation Measurement (ARM) (Mather and Voyles 2013) user facility to discuss critical climate challenges and locations where key ARM Mobile Facility (AMF) observational assets could impact Earth system modeling (ESM). As an outcome, the southeast United States (SE U.S.) was identified as a high-priority region to target climate-process studies that promote a deeper understanding of the climate system and bolster ARM interactions with the community to drive ESM advancement. The DOE ARM user facility is a globally recognized leader in deploying and operating strategically located observation sites around the world for studying the properties of aerosols and clouds and their interaction with radiation, precipitation, and the Earth’s surface. In partnering with the DOE Atmospheric System Research (ASR) program, ARM solicited a multi-agency Site Science Team approach to provide input and close interaction with ARM management towards a successful SE U.S. deployment of the ARM third Mobile Facility (AMF3) (Miller et al. 2016). These efforts included identifying key locations, science drivers and instruments, and measurement strategies to address the wider climate-process needs and ESM improvement. Community input served a vital role in establishing, refining, and informing the relevant drivers and decisions regarding this AMF3 deployment. The team has identified Northern Alabama (N. AL) as regionally representative to unlock the key opportunities that will improve our understanding and model representation of aerosol, cloud, and land surface processes and their couplings in the SE U.S. A defining aspect of the AMF3 deployment is its commitment to long-term (anticipated five-year) observations to mitigate potential seasonal-to-annual variability that often limits appropriate attribution of phenomena to local or larger-scale processes. The proposed location may leverage nearby surface networks and multi-agency and partner assets to enrich this multi-year deployment. One motivation is to understand the role of spatiotemporal variability (thermodynamic, land-surface) across aspects of the climate system, with our AMF3 team anticipating future demands on characterizing the relationships between local-to-regional cloud development and surface processes across a diverse patchwork of natural, managed, and urban landscapes as found throughout the N. AL regions. The main site targets an intact, representative, forested region – the Bankhead National Forest (BNF) – underscoring further team commitment to regionally important land atmosphere two-way interactive studies “from the canopy to the clouds”, with enhanced tower instrumentation augmenting traditional ARM capabilities adjacent to this site. Multiple supplemental sites will also be distributed across this region, prioritizing added needs for biodiversity. Anticipated high-priority cloud science themes will target N. AL as a regional SE U.S. hotbed for high-impact weather, convective cloud onset, and shallow to-deep cloud transitioning. Anticipated aerosol drivers will focus on chemical processes that control the evolution of organic aerosol, the seasonality and spatial distribution of water vapor and particle-phase water, and its role on aerosol optical properties. Anticipated land atmosphere drivers consider the two-way feedbacks between surface influence on aerosols, clouds, and precipitation properties and the associated radiative impacts on plant physiology and canopy-scale fluxes. Emphasis will include the study of the impact of surface processes on aerosols via precursor emission, and on clouds via moisture flux and thermal development.

Doppler lidar, aerosols, convection↗
DOE OSTI2015

Optimized molecular reconstruction procedure combining hybrid reverse Monte Carlo and molecular dynamics

We report an efficient atom-scale reconstruction method that consists of combining the Hybrid Reverse Monte Carlo algorithm (HRMC) with Molecular Dynamics (MD) in the framework of a simulated annealing technique. In the spirit of the experimentally constrained molecular relaxation technique [Biswas et al., Phys. Rev. B 69, 195207 (2004)], this modified procedure offers a refined strategy in the field of reconstruction techniques, with special interest for heterogeneous and disordered solids such as amorphous porous materials. While the HRMC method generates physical structures, thanks to the use of energy penalties, the combination with MD makes the method at least one order of magnitude faster than HRMC simulations to obtain structures of similar quality. Furthermore, in order to ensure the transferability of this technique, we provide rational arguments to select the various input parameters such as the relative weight ω of the energy penalty with respect to the structure optimization. By applying the method to disordered porous carbons, we show that adsorption properties provide data to test the global texture of the reconstructed sample but are only weakly sensitive to the presence of defects. In contrast, the vibrational properties such as the phonon density of states are found to be very sensitive to the local structure of the sample.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Asynchronous multilevel adaptive methods for solving partial differential equations on multiprocessors - Performance results

The fast adaptive composite grid method (FAC) is an algorithm that uses various levels of uniform grids (global and local) to provide adaptive resolution and fast solution of PDEs. Like all such methods, it offers parallelism by using possibly many disconnected patches per level, but is hindered by the need to handle these levels sequentially. The finest levels must therefore wait for processing to be essentially completed on all the coarser ones. A recently developed asynchronous version of FAC, called AFAC, completely eliminates this bottleneck to parallelism. This paper describes timing results for AFAC, coupled with a simple load balancing scheme, applied to the solution of elliptic PDEs on an Intel iPSC hypercube. These tests include performance of certain processes necessary in adaptive methods, including moving grids and changing refinement. A companion paper reports on numerical and analytical results for estimating convergence factors of AFAC applied to very large scale examples.

Mccormick, S.↗
NASA NTRS2023

Aconcagua Basin Agriculture: Analyzing Hydrological Norms, Evapotranspiration, and Soil Moisture to Assess Crop Water Demand and Water Usage in Chile’s Aconcagua River Basin

The Aconcagua basin has experienced a prolonged drought since 2010, posing a significant concern to the agricultural sector, which contributes to 12% of the national agriculture output. Reduced precipitation and warmer conditions have led to water constraints for agricultural activities. As the agriculture sector relies more on irrigation, there has been a decline in surface water availability and a shrinking groundwater supply, both primarily sourced from snow and glacial melt. This study focuses on evaluating crop water demand and water usage in the Aconcagua Basin, leveraging in-situ and available satellite data from Aqua and Terra Moderate Resolution Imaging Spectroradiometer (MODIS), 1-km downscaled Soil Moisture Active Passive (SMAP), Global Precipitation Measurement Mission IMERG (GPM), and Landsat 8 Operational Land Imager (OLI). Remote sensing hydrologic norms from pre- and intra-drought conditions provided a baseline for our water use analysis. We have identified areas of the basin where water usage is exaggerated compared to other regions by analyzing soil moisture, evapotranspiration, and vegetation index trends across agricultural lands. The study highlights the presence of irrigation and exaggerated water usage categorized by crop type while normalizing the data by the quantity of water allocation (m3/s) to each subregion. The findings have the potential to assist our partners, the Centro de Información de Recursos Naturales (CIREN) and the Ministry of Agriculture, in refining water allocation approaches; showcasing the feasibility of leveraging remote sensing and earth observation datasets to monitor agricultural practices within the context of water scarcity.

Remote sensing↗
DOE OSTI2022

Enhancing $\phi$-sensitivity of ignition delay times through dilution of fuel-air mixture

The high $\phi$-sensitivity (η) of ignition delay time (τ IDT ) is one of the desirable fuel properties for the high-load extension of advanced compression ignition engine. Recent studies revealed the effectiveness of a high dilution rate (x D ) for enhancing η at the engine-relevant conditions. This study aims to quantify the effect of dilution on the η of isooctane using a combined experiment-simulation approach. The ignition delay of the isooctane/air mixture was measured with an Advanced Fuel Ignition Delay Analyzer (AFIDA) over the temperature range of 623 - 923 K at 10 bar pressure with global $\phi$ = 0.3 - 0.6, with and without 28.6% of additional N 2 dilution. For precise evaluation of experimental η, the facility effect of the AFIDA experiment was characterized with three-dimensional computational fluid dynamics (3-D CFD) simulation. The temperature in the combustion chamber from 3-D CFD indicated a substantial temporal dependency, varying up to ~52 K by charge-cooling of fuel injection and heat transfer from the wall. We introduced the dimensionless number θ(t) for characterizing the temporal profile of chamber temperature. Consideration of the facility effect using θ(t) resulted in better agreement between the experimental ..eta.. and zero-dimensional (0-D) kinetics simulation. The refined η were then further utilized to quantify the effectiveness of dilution to η The extent of η enhancement with dilution strategy was maximized at the low-temperature chemistry regime, increasing η by 77% with a 28.6% dilution rate. Further analysis on the dilution effect was carried out using 0-D kinetics simulation, revealing the critical dimensionless numbers relevant to the effectiveness of dilution to η enhancement. Here, this study is the first experiment-simulation combined research to quantify the effect of dilution on η, facilitating the kinetics model refinement for better reproduction of $\phi$-sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Getting lost tracking the carbon footprint of hydropower

In the transition to low-carbon electricity, well-quantified estimates of carbon dynamics are needed to ensure that emissions reduction targets are achieved. Here we review the state of the science on carbon accounting for hydropower reservoirs and identify limitations and future solutions. Nearly all research on reservoir greenhouse-gas (GHG) emissions has focused on individual reservoirs in isolation without considering their position in a freshwater network draining organic matter from upstream watersheds or the coordinated operation of reservoir cascades. Second, carbon inventories have extrapolated from a small, non-probabilistic sample of highly variable measurements of GHG emissions to unsampled reservoirs. A stronger statistical foundation is needed to estimate a global inventory and its uncertainty. Third, attribution to hydropower is based on ranks assigned to reservoir purpose. Instead, the physical influence of hydropower on carbon dynamics could be directly measured. Fourth, current carbon-accounting practices neglect time. A time-varying approach would quantify variation in emissions for electricity portfolios from changes in the fuel mix at different times and account for ancillary services, i.e., the ability to support the grid when variable renewables are not available without using natural gas. Reservoirs also sequester a significant portion of inflowing carbon in sediments and slow the carbon cycle by delaying the return of carbon to the atmosphere for decades to centuries. Together, these refinements would help to illuminate pathways toward meeting energy demand with the longest-possible delay in returning carbon to the atmosphere and without adding ancient sources to the pool of carbon cycling through aquatic ecosystems.

13 HYDRO ENERGY↗
DOE OSTI2016

ON HELIUM MIXING IN QUASI-GLOBAL SIMULATIONS OF THE INTRACLUSTER MEDIUM

The assumption of a spatially uniform helium distribution in the intracluster medium (ICM) can lead to biases in the estimates of key cluster parameters if composition gradients are present. The helium concentration profile in galaxy clusters is unfortunately not directly observable. Current models addressing the putative sedimentation are one-dimensional and parametrize the presence of magnetic fields in a crude way, ignoring the weakly collisional, magnetized nature of the medium. When these effects are considered, a wide variety of instabilities can play an important role in the plasma dynamics. In a series of recent papers, we have developed the local, linear theory of these instabilities and addressed their nonlinear development with a modified version of Athena. Here, we extend our study by developing a quasi-global approach that we use to simulate the mixing of helium as induced by generalizations of the heat-flux-driven buoyancy instability (HBI) and the magnetothermal instability, which feed off thermal and composition gradients. In the inner region of the ICM, mixing can occur over a few gigayears, after which the average magnetic field inclination angle is ∼30°–50°, resulting in an averaged Spitzer parameter higher by about 20% than the value obtained in homogeneous simulations. In the cluster outskirts the instabilities are rather inefficient, due to the shallow gradients. This suggests that composition gradients in cluster cores might be shallower than one-dimensional models predict. More quantitative statements demand more refined models that can incorporate the physics driving the sedimentation process and simultaneously account for the weakly collisional nature of the plasma.

79 ASTRONOMY AND ASTROPHYSICS↗