Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference model”

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,693 records · Page 94

Reconfigurable Model Execution in the OpenMDAO Framework

NASA's OpenMDAO framework facilitates constructing complex models and computing their derivatives for multidisciplinary design optimization. Decomposing a model into components that follow a prescribed interface enables OpenMDAO to assemble multidisciplinary derivatives from the component derivatives using what amounts to the adjoint method, direct method, chain rule, global sensitivity equations, or any combination thereof, using the MAUD architecture. OpenMDAO also handles the distribution of processors among the disciplines by hierarchically grouping the components, and it automates the data transfer between components that are on different processors. These features have made OpenMDAO useful for applications in aircraft design, satellite design, wind turbine design, and aircraft engine design, among others. This paper presents new algorithms for OpenMDAO that enable reconfigurable model execution. This concept refers to dynamically changing, during execution, one or more of: the variable sizes, solution algorithm, parallel load balancing, or set of variables-i.e., adding and removing components, perhaps to switch to a higher-fidelity sub-model. Any component can reconfigure at any point, even when running in parallel with other components, and the reconfiguration algorithm presented here performs the synchronized updates to all other components that are affected. A reconfigurable software framework for multidisciplinary design optimization enables new adaptive solvers, adaptive parallelization, and new applications such as gradient-based optimization with overset flow solvers and adaptive mesh refinement. Benchmarking results demonstrate the time savings for reconfiguration compared to setting up the model again from scratch, which can be significant in large-scale problems. Additionally, the new reconfigurability feature is applied to a mission profile optimization problem for commercial aircraft where both the parametrization of the mission profile and the time discretization are adaptively refined, resulting in computational savings of roughly 10% and the elimination of oscillations in the optimized altitude profile.

Hwang, John T.↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

A Dirichlet process mixture model for automatic {sup 18}F-FDG PET image segmentation: Validation study on phantoms and on lung and esophageal lesions

Purpose: The aim of this study was to implement a Dirichlet process mixture (DPM) model for automatic tumor edge identification on {sup 18}F-fluorodeoxyglucose positron emission tomography ({sup 18}F-FDG PET) images by optimizing the parameters on which the algorithm depends, to validate it experimentally, and to test its robustness. Methods: The DPM model belongs to the class of the Bayesian nonparametric models and uses the Dirichlet process prior for flexible nonparametric mixture modeling, without any preliminary choice of the number of mixture components. The DPM algorithm implemented in the statistical software package R was used in this work. The contouring accuracy was evaluated on several image data sets: on an IEC phantom (spherical inserts with diameter in the range 10–37 mm) acquired by a Philips Gemini Big Bore PET-CT scanner, using 9 different target-to-background ratios (TBRs) from 2.5 to 70; on a digital phantom simulating spherical/uniform lesions and tumors, irregular in shape and activity; and on 20 clinical cases (10 lung and 10 esophageal cancer patients). The influence of the DPM parameters on contour generation was studied in two steps. In the first one, only the IEC spheres having diameters of 22 and 37 mm and a sphere of the digital phantom (41.6 mm diameter) were studied by varying the main parameters until the diameter of the spheres was obtained within 0.2% of the true value. In the second step, the results obtained for this training set were applied to the entire data set to determine DPM based volumes of all available lesions. These volumes were compared to those obtained by applying already known algorithms (Gaussian mixture model and gradient-based) and to true values, when available. Results: Only one parameter was found able to significantly influence segmentation accuracy (ANOVA test). This parameter was linearly connected to the uptake variance of the tested region of interest (ROI). In the first step of the study, a calibration curve was determined to automatically generate the optimal parameter from the variance of the ROI. This “calibration curve” was then applied to contour the whole data set. The accuracy (mean discrepancy between DPM model-based contours and reference contours) of volume estimation was below (1 ± 7)% on the whole data set (1 SD). The overlap between true and automatically segmented contours, measured by the Dice similarity coefficient, was 0.93 with a SD of 0.03. Conclusions: The proposed DPM model was able to accurately reproduce known volumes of FDG concentration, with high overlap between segmented and true volumes. For all the analyzed inserts of the IEC phantom, the algorithm proved to be robust to variations in radius and in TBR. The main advantage of this algorithm was that no setting of DPM parameters was required in advance, since the proper setting of the only parameter that could significantly influence the segmentation results was automatically related to the uptake variance of the chosen ROI. Furthermore, the algorithm did not need any preliminary choice of the optimum number of classes to describe the ROIs within PET images and no assumption about the shape of the lesion and the uptake heterogeneity of the tracer was required.

60 APPLIED LIFE SCIENCES↗

A Characterization of the Terrestrial Environment of Kodiak Island, Alaska for the Design, Development and Operation of Launch Vehicles

A quantitative characterization of the terrestrial environment is an important component in the success of a launch vehicle program. Environmental factors such as winds, atmospheric thermodynamics, precipitation, fog, and cloud characteristics are among many parameters that must be accurately defined for flight success. The National Aeronautics and Space Administration (NASA) is currently coordinating weather support and performing analysis for the launch of a NASA payload from a new facility located at Kodiak Island, Alaska in late 2001 (NASA, 1999). Following the first launch from the Kodiak Launch Complex, an Air Force intercontinental ballistic missile on November 5, 1999, the site's developer, the Alaska Aerospace Development Corporation (AADC), is hoping to acquire a sizable share of the many launches that will occur over the next decade. One such customer is NASA, which is planning to launch the Vegetation Canopy Lidar satellite aboard an Athena I rocket, the first planned mission to low earth orbit from the new facility. To support this launch, a statistical model of the atmospheric and surface environment for Kodiak Island, AK has been produced from rawinsonde and surface-based meteorological observations for use as an input to future launch vehicle design and/or operations. In this study, the creation of a "reference atmosphere" from rawinsonde observations is described along with comparisons between the reference atmosphere and existing model representations for Kodiak. Meteorological conditions that might result in a delay on launch day (cloud cover, visibility, precipitation, etc.) are also explored and described through probabilities of launch by month and hour of day. This atmospheric "mission analysis" is also useful during the early stages of a vehicle program, when consideration of the climatic characteristics of a location can be factored into vehicle designs. To be most beneficial, terrestrial environment definitions should a) be available at the inception of the program and based on the desired operational performance of the launch vehicle, b) be issued under the signature of the program manager and be part of the controlled program definition and requirements documentation, and c) specify the terrestrial environment for all phases of activity including prelaunch, launch, ascent, on-orbit, decent, and landing. Since the beginning of the space era, NASA has utilized some of the most detailed assessments of the terrestrial climatic environment in design, development, and operations of both expendable and reusable launch vehicles.

Rawlins, Michael A.↗

Basin and Range Investigation for Developing Geothermal Energy: Exploration Data

This data package includes exploration material from the Basin & Range Investigation for Developing Geothermal Energy [in Hidden Systems] project (BRIDGE), which is part of a broader initiative to advance the exploration of hidden geothermal resources in the Basin & Range Province of the western U.S. Data modalities include a helicopter-borne time-domain electromagnetic survey, magnetotellurics, 2-meter temperature measurements, ground-based gravity and legacy aeromagnetic surveys, geochemistry, geologic mapping, LiDAR analysis, 3D models, associated geospatial data, and a bibliography of existing data and references utilized in prospect characterization and conceptual modeling. Key files are in CSV, Geosoft, and Geotools formats. Please refer to READMEs for dataset-specific information. Where applicable, acquisition data and inversion models for a particular prospect or area of interest are organized separately. This BRIDGE data package is the product of a collaboration led by Sandia National Laboratories with partners from Geologica Geothermal Group, Inc., the U.S. Navy Geothermal Program Office, and consultants Steven Sewell (Australis Geoscience Ltd) and William Cumming (Cumming Geoscience). The project's areas of interest (AOIs) are based off priority areas of interest in the southwestern portion of the Nevada Play Fairway map, distribution across tectonic provinces, accessibility, and the project team's extensive experience in the region. AOIs cover about a dozen basins that include unexplored prospects, partially explored prospects, and some developed analogue resources that provide validation cases. Many unexplored and partially explored prospects are on U.S. Department of Defense (DoD) land, though adjacent lands are included as well.

15 GEOTHERMAL ENERGY↗

Nudging allows direct evaluation of coupled climate models with in situ observations: a case study from the MOSAiC expedition

Comparing the output of general circulation models to observations is essential for assessing and improving the quality of models. While numerical weather prediction models are routinely assessed against a large array of observations, comparing climate models and observations usually requires long time series to build robust statistics. Here, we show that by nudging the large-scale atmospheric circulation in coupled climate models, model output can be compared to local observations for individual days. We illustrate this for three climate models during a period in April 2020 when a warm air intrusion reached the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition in the central Arctic. Radiosondes, cloud remote sensing and surface flux observations from the MOSAiC expedition serve as reference observations. The climate models AWI-CM1/ECHAM and AWI-CM3/IFS miss the diurnal cycle of surface temperature in spring, likely because both models assume the snowpack on ice to have a uniform temperature. CAM6, a model that uses three layers to represent snow temperature, represents the diurnal cycle more realistically. During a cold and dry period with pervasive thin mixed-phase clouds, AWI-CM1/ECHAM only produces partial cloud cover and overestimates downwelling shortwave radiation at the surface. AWI-CM3/IFS produces a closed cloud cover but misses cloud liquid water. Our results show that nudging the large-scale circulation to the observed state allows a meaningful comparison of climate model output even to short-term observational campaigns. We suggest that nudging can simplify and accelerate the pathway from observations to climate model improvements and substantially extends the range of observations suitable for model evaluation.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Building Energy Modeling Programs: Building Loads

This technical report presented the methodologies, processes, and results of comparing three Building Energy Modeling Programs (BEMPs) for load calculations: EnergyPlus, DeST and DOE-2.1E. This joint effort, between Lawrence Berkeley National Laboratory, USA and Tsinghua University, China, was part of research projects under the US-China Clean Energy Research Center on Building Energy Efficiency (CERC-BEE). Energy Foundation, an industrial partner of CERC-BEE, was the co-sponsor of this study work. It is widely known that large discrepancies in simulation results can exist between different BEMPs. The result is a lack of confidence in building simulation amongst many users and stakeholders. In the fields of building energy code development and energy labeling programs where building simulation plays a key role, there are also confusing and misleading claims that some BEMPs are better than others. In order to address these problems, it is essential to identify and understand differences between widely-used BEMPs, and the impact of these differences on load simulation results, by detailed comparisons of these BEMPs from source code to results. The primary goal of this work was to research methods and processes that would allow a thorough scientific comparison of the BEMPs. The secondary goal was to provide a list of strengths and weaknesses for each BEMP, based on in-depth understandings of their modeling capabilities, mathematical algorithms, advantages and limitations. This is to guide the use of BEMPs in the design and retrofit of buildings, especially to support China’s building energy standard development and energy labeling program. The research findings could also serve as a good reference to improve the modeling capabilities and applications of the three BEMPs. The methodologies, processes, and analyses employed in the comparison work could also be used to compare other programs. The load calculation method of each program was analyzed and compared to identify the differences in solution algorithms, modeling assumptions and simplifications. Identifying inputs of each program and their default values or algorithms for load simulation was a critical step. These tend to be overlooked by users, but can lead to large discrepancies in simulation results. As weather data was an important input, weather file formats and weather variables used by each program were summarized. Some common mistakes in the weather data conversion process were discussed. ASHRAE Standard 140-2007 tests were carried out to test the fundamental modeling capabilities of the load calculations of the three BEMPs, where inputs for each test case were strictly defined and specified. The tests indicated that the cooling and heating load results of the three BEMPs fell mostly within the range of spread of results from other programs. Based on ASHRAE 140-2007 test results, the finer differences between DeST and EnergyPlus were further analyzed by designing and conducting additional tests. Potential key influencing factors (such as internal gains, air infiltration, convection coefficients of windows and opaque surfaces) were added one at a time to a simple base case with an analytical solution, to compare their relative impacts on load calculation results. Finally, special tests were designed and conducted aiming to ascertain the potential limitations of each program to perform accurate load calculations. The heat balance module was tested for both single and double zone cases. Furthermore, cooling and heating load calculations were compared between the three programs by varying the heat transfer between adjacent zones, the occupancy of the building, and the air-conditioning schedule.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

IM3 Phase 2 Official Simulations: GCAM-Demeter-SELECT Annualized Land Use and Land Cover, Wood Harvest and Fertilization Data with Dynamic Urbanization Harmonized to CLM Land Definitions at 0.125 Degrees

Annualized land use land cover data, including wood harvest and fertilizer use data from the Global Change Analysis Model (GCAM) downscaled to 0.125 degrees for couping with the Community Land Model (CLM). GCAM here refers to GCAM-USA v5.3.im3 which has an enhanced electricity sector and an updated data system needed to represent regional to local scale dynamics. Data is also harmonized with future urbanization projections from the Spatially-Explicit, Long-term, Empirical City developmenT (SELECT) model. Projections/Data are generated using the demeter land use and land cover downscaling model. Original projections were generated at 0.05 degrees before being aggregared to 0.125 degrees. Projections are available for 8 alternative scenarios. Two versions of final data are included- one with managed forests or harvested forest area per pixel broken out and one with the same aggregated into total forests. Following folders are included: demeter_78_PFT_output:This is the final output of dynamic land use land cover change for 78 PFTs as required by CLM raw_outputs_incl_managed_forest: This is the final output but with managed forests broken out as a different PFT. Essentially a 79th PFT is added. wood_harvest_outputs: Wood harvest output per pixel in gC/m2 fertilization_outputs: Fertilizer use per pixel in gN/m2 Each NetCDF file in each folder represents a projection for a separate year, scenario. Land use outputs are organized as PFT level data saved as subdata. Link to GCAM version used- https://data.msdlive.org/records/yb23g-44274 Link to SELECT documentation -https://www.sciencedirect.com/science/article/pii/S1364815219301707 Link to CLM documentation- - https://www.cesm.ucar.edu/models/clm In case of questions contact- kanishka.narayan@pnnl.gov

GCAM-USA↗

Evaluating the Effects of Precipitation and Evapotranspiration on Soil Moisture Variability Within CMIP5 Using SMAP and ERA5 Data

=The effects of precipitation (Pr) and evapotranspiration (ET) on surface soil moisture (SSM) play an essential role in the land-atmosphere system. Here we evaluate multimodel differences of these effects within the Coupled Model Intercomparison Project Phase 5 (CMIP5) compared to Soil Moisture Active Passive (SMAP) products and ECMWF Reanalysis v5 (ERA5) as references in a frequency domain. The variability of SSM, Pr, and ET within three frequency bands (1/7 ~ 1/30 days -1 , 1/30 ~ 1/90 days -1 , and 1/90 ~ 1/365 days -1 ) after normalization is quantified using Fourier transform. We analyze the impact of ET and Pr on SSM variability based on a transfer function assuming that these variables form a linear time-invariant (LTI) system. For the total effects of ET and Pr on SSM variability, the CMIP5 estimations are smaller than the reference data in the two higher frequency bands and are larger than the reference data in the lowest frequency band. Besides, the effects on SSM by Pr and ET are found to be different across the three frequency bands. In each frequency band, the variability of the factor that dominates SSM (i.e., Pr or ET) from CMIP5 is smaller than that from the references. This study identifies the spatiotemporal distribution of differences between CMIP5 models and references (SMAP and ERA5) in simulating ET and Pr effects on SSM within three frequency bands. This study provides insightful information on how soil moisture variability is affected by varying precipitation and evapotranspiration at different time scales within Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Use NU-WRF and GCE Model to Simulate the Precipitation Processes During MC3E Campaign

One of major CRM approaches to studying precipitation processes is sometimes referred to as "cloud ensemble modeling". This approach allows many clouds of various sizes and stages of their lifecycles to be present at any given simulation time. Large-scale effects derived from observations are imposed into CRMs as forcing, and cyclic lateral boundaries are used. The advantage of this approach is that model results in terms of rainfall and QI and Q2 usually are in good agreement with observations. In addition, the model results provide cloud statistics that represent different types of clouds/cloud systems during their lifetime (life cycle). The large-scale forcing derived from MC3EI will be used to drive GCE model simulations. The model-simulated results will be compared with observations from MC3E. These GCE model-simulated datasets are especially valuable for LH algorithm developers. In addition, the regional scale model with very high-resolution, NASA Unified WRF is also used to real time forecast during the MC3E campaign to ensure that the precipitation and other meteorological forecasts are available to the flight planning team and to interpret the forecast results in terms of proposed flight scenarios. Post Mission simulations are conducted to examine the sensitivity of initial and lateral boundary conditions to cloud and precipitation processes and rainfall. We will compare model results in terms of precipitation and surface rainfall using GCE model and NU-WRF

Tao, Wei-Kuo↗

NASA 5.2%-Scale High Lift Common Research Model (CRM-HL) Test in the National Transonic Facility (NTF)

Problem: Predicting CL,max and assessing transition, turbulence models, and Reynolds number effect predictions in CFD requires wind tunnel data for comparison. Objective: Expand the CFD validation database by running multiple models designed using the same reference geometry in various wind tunnels around the world. The wind tunnel data and model geometry will be open source to allow for CFD validation use. Approach: Run the NASA 5.2%-scale High Lift Common Research Model (CRM-HL) semi-span model at the National Transonic Facility (NTF). Test the model at chord Reynolds numbers (Rec) between 1.6 million and 30 million and at varying Q/Es to assess Reynolds number effects and aeroelasticity affects on the model, while collecting force and moment, pressure and wing deformation data. Test the model with ice shapes installed on the leading edges. Results: Ran eight model configurations, including four landing, two landing with ice shapes installed, and two takeoff configurations. Ran in air (120°F) and nitrogen operations (-50°F, -180°F and -250°F) at 7 different chord Reynolds numbers. The data compare well with previous tests of the same model at DLR (The Germany Aerospace Center) and the 14- by 22-Foot Subsonic Tunnel (14x22). Significance: Provided valuable high Reynolds number data on a high lift configuration to the worldwide research community. Provide icing data at flight Reynolds number to the research community.

CRM-HL↗

On-board orbit determination for application satellites

A proposal for onboard orbit determination which would rely entirely on reference Doppler beacons located within the Continental United States is presented. Algorithms were conceived which were compatible with small computing machines and sufficiently accurate to maintain an adequate satellite ephemeris. An on board computer equipped with a recursive filter with a fading memory to account for dynamic modeling errors and two reference beacons was found to be necessary for the system to work.

Morduch, G. E.↗

Generic Argillite/Shale Disposal Reference Case

Radioactive waste disposal in a deep subsurface repository hosted in clay/shale/argillite is a subject of widespread interest given the desirable isolation properties, geochemically reduced conditions, and widespread geologic occurrence of this rock type (Hansen 2010; Bianchi et al. 2013). Bianchi et al. (2013) provides a description of diffusion in a clay-hosted repository based on single-phase flow and full saturation using parametric data from documented studies in Europe (e.g., ANDRA 2005). The predominance of diffusive transport and sorption phenomena in this clay media are key attributes to impede radionuclide mobility making clay rock formations target sites for disposal of high-level radioactive waste. The reports by Hansen et al. (2010) and those from numerous studies in clay-hosted underground research laboratories (URLs) in Belgium, France and Switzerland outline the extensive scientific knowledge obtained to assess long-term clay/shale/argillite repository isolation performance of nuclear waste. In the past several years under the UFDC, various kinds of models have been developed for argillite repository to demonstrate the model capability, understand the spatial and temporal alteration of the repository, and evaluate different scenarios. These models include the coupled Thermal-Hydrological-Mechanical (THM) and Thermal-Hydrological-Mechanical-Chemical (THMC) models (e.g. Liu et al. 2013; Rutqvist et al. 2014a, Zheng et al. 2014a) that focus on THMC processes in the Engineered Barrier System (EBS) bentonite and argillite host hock, the large scale hydrogeologic model (Bianchi et al. 2014) that investigates the hydraulic connection between an emplacement drift and surrounding hydrogeological units, and Disposal Systems Evaluation Framework (DSEF) models (Greenberg et al. 2013) that evaluate thermal evolution in the host rock approximated as a thermal conduction process to facilitate the analysis of design options. However, the assumptions and the properties (parameters) used in these models are different, which not only make inter-model comparisons difficult, but also compromise the applicability of the lessons learned from one model to another model. The establishment of a reference case would therefore be helpful to set up a baseline for model development. A generic salt repository reference case was developed in Freeze et al. (2013) and the generic argillite repository reference case is presented in this report. The definition of a reference case requires the characterization of the waste inventory, waste form, waste package, repository layout, EBS backfill, host rock, and biosphere. This report mainly documents the processes in EBS bentonite and host rock that are potentially important for performance assessment and properties that are needed to describe these processes, with brief description other components such as waste inventory, waste form, waste package, repository layout, aquifer, and biosphere. A thorough description of the generic argillite repository reference case will be given in Jové Colon et al. (2014).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Introduction to the FCRD Transmutation Fuels Handbook 2015

Transmutation of minor actinides such as Np, Am, and Cm in spent nuclear fuel is of international interest because of its potential for reducing the long-term health and safety hazards caused by the radioactivity of the spent fuel. For this reason, the DOE Fuel Cycle Research and Development (FCRD) Advanced Fuels Campaign (AFC) has been performing research on nuclear fuel technologies that could be used in conjunction with a future actinide transmutation mission. One important approach to transmutation involves incorporating the minor actinides into U-Pu-Zr alloys, which can be used as fuel in fast reactors. It is, therefore, important to understand the properties of U-Pu-Zr alloys. The Idaho National Laboratory has recently published a handbook intended to summarize key properties of these alloys for researchers. This handbook focuses on phases and key thermophysical properties of elements and alloys in the U-Pu-Zr system. It is envisioned as the first of a series, which will serve as a summary of current knowledge about important properties of transmutation fuels. Future volumes will include information about properties and materials not covered in the first edition, as well as updates to information in this edition of the handbook. The handbook is intended to serve two audiences. One audience is researchers, who may find it useful to know what has been measured, where it has been published, and what kinds of information are needed. A second audience is modelers, who need a concise summary of information about specific properties and the accuracy with which they are known. The two audiences have different needs, and the handbook is a compromise intended to be useful to both. Therefore, it includes summaries of research results and needs, but generally does not present the original data. It also includes numerous graphs and, where available, equations that provide convenient approximations to materials properties for modelers. It is important to understand that most of the equations provided in this handbook are derived by fitting polynomials to data and have no physical basis; while they can be used with caution for design studies, fuel modelers will probably want to use the original experimental data to develop or validate their own equations and/or models. Although the handbook includes some citations to research involving development of models, its primary focus is experimental data. Many of the measurements of the U-Pu-Zr system were made before 1980. Where possible, the handbook relies on previous reviews of experimental data rather than presenting new interpretations of the old information. It is expected that researchers will refer to the original literature for details that are not presented in the handbook. Metals and alloys in the U-Pu-Zr system are uniquely difficult to study for a number of reasons. Experimental designs may need to be modified to allow researchers to work with these hazardous materials safely. Kinetics of phase transformations are commonly sluggish, and it is highly unlikely that all of the samples in previous research had equilibrium phases and microstructures. Thermal cycling effects (changes in sample characteristics that are not fully reversed during cooling) are commonly reported and may be significant. Careful control of impurities (particularly oxygen) during sample fabrication and data collection is both crucial and difficult, and concentrations below 100 ppm in samples may significantly affect thermophysical properties. For all of these reasons, differences between measured values for similar alloys are commonly far larger than the estimated errors for individual measurements. Further research on well-characterized materials with careful control of impurities is needed for all of the alloys discussed here. The 2015 FCRD Handbook summarizes currently available information about U, Pu, Zr, and their binary and ternary alloys. It contains information about phase diagrams and related information (including phases and phase transformations); heat capacity, entropy, and enthalpy; thermal expansion; and thermal conductivity and diffusivity. In addition to presenting information about materials properties, it attempts to provide information about how well the property is known and how much variation exists between measurements. Although the handbook includes some references to publications about modeling, its primary focus is experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploring Electrode-Level State-of-Charge and State-of-Health Dynamics in Lithium-Ion Battery Cells: Modeling and Experimental Identification

A computationally efficient model serves as a critical prerequisite for battery performance analysis and advanced battery management algorithm design. Although battery models that capture cell-level behavior have been widely explored in existing literature, electrode-level battery models have received much lesser attention till to date. However, such electrode-level models can significantly increase battery performance and life by enabling electrode-level health-conscious control. Such electrode-level control can effectively expand usable energy and power limits of the battery cells by utilizing the knowledge of individual electrodes' charge and health. In this context, this paper presents a comprehensive battery model developed with a reference electrode insertion that captures (i) electrode-level charge/discharge dynamics, (ii) stoichiometric and temporal dependencies of electrode-level resistances, (iii) solid electrolyte interface (SEI) layer growth as key degradation phenomenon, and (iv) capacity fade and resistance rise in each electrode due to nominal battery aging. The proposed model is identified, and a preliminary validation is performed utilizing terminal voltage and negative electrode potential data collected from a pouch cell under one continuous cycling and accelerated aging conditions where the cell experienced 14% capacity loss.

aging↗

First measurements of N -subjettiness in central Pb-Pb collisions at $\sqrt{s{NN}}$= 2.76 TeV

The ALICE Collaboration reports the first fully-corrected measurements of the N-subjettiness observable for track-based jets in heavy-ion collisions. This study is performed using data recorded in pp and Pb-Pb collisions at centre-of-mass energies of $\sqrt{s}$= 7 TeV and $\sqrt{s{NN}}$ = 2.76 TeV, respectively. In particular the ratio of 2-subjettiness to 1-subjettiness, τ 2 /τ 1 , which is sensitive to the rate of two-pronged jet substructure, is presented. Energy loss of jets traversing the strongly interacting medium in heavy-ion collisions is expected to change the rate of two-pronged substructure relative to vacuum. The results are presented for jets with a resolution parameter of R = 0.4 and charged jet transverse momentum of 40 ≤ pτ, jet ≤ 60 GeV/c, which constitute a larger jet resolution and lower jet transverse momentum interval than previous measurements in heavy-ion collisions. This has been achieved by utilising a semi-inclusive hadron-jet coincidence technique to suppress the larger jet combinatorial background in this kinematic region. No significant modification of the τ 2 /τ 1 observable for track-based jets in Pb-Pb collisions is observed relative to vacuum PYTHIA6 and PYTHIA8 references at the same collision energy. The measurements of τ 2 /τ 1 , together with the splitting aperture angle ΔR, are also performed in pp collisions at $\sqrt{s}$ = 7 TeV for inclusive jets. These results are compared with PYTHIA calculations at $\sqrt{s}$ = 7 TeV, in order to validate the model as a vacuum reference for the Pb-Pb centre-of-mass energy. The PYTHIA references for τ 2 /τ 1 are shifted to larger values compared to the measurement in pp collisions. This hints at a reduction in the rate of two-pronged jets in Pb-Pb collisions compared to pp collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A method for solving stochastic equations by reduced order models and local approximations

A method is proposed for solving equations with random entries, referred to as stochastic equations (SEs). The method is based on two recent developments. The first approximates the response surface giving the solution of a stochastic equation as a function of its random parameters by a finite set of hyperplanes tangent to it at expansion points selected by geometrical arguments. The second approximates the vector of random parameters in the definition of a stochastic equation by a simple random vector, referred to as stochastic reduced order model (SROM), and uses it to construct a SROM for the solution of this equation. The proposed method is a direct extension of these two methods. It uses SROMs to select expansion points, rather than selecting these points by geometrical considerations, and represents the solution by linear and/or higher order local approximations. The implementation and the performance of the method are illustrated by numerical examples involving random eigenvalue problems and stochastic algebraic/differential equations. The method is conceptually simple, non-intrusive, efficient relative to classical Monte Carlo simulation, accurate, and guaranteed to converge to the exact solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of constant-Weber-number scaling for icing tests

Previous studies showed that for conditions simulating an aircraft encountering super-cooled water droplets the droplets may splash before freezing. Other surface effects dependent on the water surface tension may also influence the ice accretion process. Consequently, the Weber number appears to be important in accurately scaling ice accretion. A scaling method which uses a constant-Weber-number approach has been described previously; this study provides an evaluation of this scaling method. Tests are reported on cylinders of 2.5 to 15-cm diameter and NACA 0012 airfoils with chords of 18 to 53 cm in the NASA Lewis Icing Research Tunnel (IRT). The larger models were used to establish reference ice shapes, the scaling method was applied to determine appropriate scaled test conditions using the smaller models, and the ice shapes were compared. Icing conditions included warm glaze, horn glaze and mixed. The smallest size scaling attempted was 1/3, and scale and reference ice shapes for both cylinders and airfoils indicated that the constant-Weber-number scaling method was effective for the conditions tested.

Anderson, David N.↗