Search NASA⌕ Search

SEARCH · Search NASA

Results for “State Variable”

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 559 records · Page 31

IDAES-PSE 2.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.0.0 Release Highlights Removal of deprecated features from IDAES v1 Update to Pyomo v6.5 – this required a number of updates to support the new NL solver writer and to address some changes in Pyomo Creation of new testing suite for backward compatibility, model robustness and verification More general implementation of the Helmholtz EoS. This brings some new features like standard property diagrams, choice of mass or mole basis, and new state variable options Standardizing names in Heat Exchanger models (breaking change from v2.0.0a2): Control Volumes named hot_side and cold_side Ports names hot_side_inlet, hot_side_outlet, cold_side_inlet and cold_side_outlet Config Blocks names hot_side_config and cold_side_config Config arguments for user provided names for each side: hot_side_name and cold_side_name. Updating Keras surrogate tool to use v1.1 of OMLT New prototype API for model initialization (idaes.core.initialization) The new API uses "Model Initializer" objects instead of class methods, allowing for the definition of multiple initialization routines for a single model A number of common, model agnostic initialization routines have also been defined, including initialization from data, block-decomposition and a general hierarchical approach equivalent to the existing method for common unit models New metadata for thermophysical properties – valid_range This can be used to record the range of values over which a property value can be trusted, such as the range of experimental data used to regress parameters A number of new utility functions have been added to check for properties with values outside the valid range and to set bounds based on this metadata Updated construction of balance expressions in Control Volumes to remove unneeded terms In the past, unneeded terms were added as a constant 0 term, however they will now be dropped entirely from the expression This was necessary due to more strict unit checking in the new Pyomo solver writer which no longer ignores 0 terms Updates to metadata for thermophysical properties to better define known properties and units of measurement This results in more strict enforcement of standard naming for thermophysical and reaction properties Users can still define custom properties, but these must be done explicitly using the define_custom_properties() method instead of being implicitly created by add_property() Updated convergence tester utility tool to support definition of benchmark files (JSON format) and comparison of performance to benchmarks Set default iteration limit for IPOPT in IDAES config to 200 iterations Update scaling of example models to work with new Pyomo NL solver writer Improve testing of extensions and examples infrastructure to avoid need for downloading files Updated distillation column to centralize common functionality and remove a number of Pyomo warnings

IDAES↗

gdess: A framework for evaluating simulated atmospheric CO 2 in Earth System Models

Atmospheric carbon dioxide (CO 2 ) plays a key role in the global carbon cycle and global warming. Climate-carbon feedbacks are often studied and estimated using Earth System Models (ESMs), which couple together multiple model components—including the atmosphere, ocean, terrestrial biosphere, and cryosphere—to jointly simulate mass and energy exchanges within and between these components. Despite tremendous advances, model intercomparisons and benchmarking are aspects of ESMs that warrant further improvement (Fer et al., 2021; Smith et al., 2014). Such benchmarking is critical because comparing the value of state variables in these simulations against observed values provides evidence for appropriately refining model components; moreover, researchers can learn much about Earth system dynamics in the process (Randall et al., 2019). We introduce `gdess` (a.k.a., Greenhouse gas Diagnostics for Earth System Simulations), which parses observational datasets and ESM simulation output, combines them to be in a consistent structure, computes statistical metrics, and generates diagnostic visualizations. In its current incarnation, `gdess` facilitates evaluating a model's ability to reproduce observed temporal and spatial variations of atmospheric CO 2 . The diagnostics implemented modularly in `gdess` support more rapid assessment and improvement of model-simulated global CO 2 sources and sinks associated with land and ocean ecosystem processes. We intend for this set of automated diagnostics to form an extensible, open source framework for future comparisons of simulated and observed concentrations of various greenhouse gases across Earth system models.

97 MATHEMATICS AND COMPUTING↗

A systematic multiscale modeling and experimental approach to protect grain boundaries in magnesium alloys from corrosion

A multiscale modeling Internal State Variable (ISV) constitutive model was developed that captures the fundamental structure-property relationships. The macroscale ISV model used lower length scale simulations (Butler-Volmer and Electronics Structures results) in order to inform the ISVs at the macroscale. The chemomechanical ISV model was calibrated and validated from experiments with magnesium (Mg) alloys that were investigated under corrosive environments coupled with experimental electrochemical studies. Because the ISV chemomechanical model is physically based, it can be used for other material systems to predict corrosion behavior. As such, others can use the chemomechanical model for analyzing corrosion effects on their designs.

36 MATERIALS SCIENCE↗

Interpolated Sounding and Gridded Sounding Value-Added Products

Standard Atmospheric Radiation Measurement (ARM) Climate Research Facility sounding files provide atmospheric state data in one dimension of increasing time and height per sonde launch. Many applications require a quick estimate of the atmospheric state at higher time resolution. The INTERPOLATEDSONDE (i.e., Interpolated Sounding) Value-Added Product (VAP) transforms sounding data into continuous daily files on a fixed time-height grid, at 1-minute time resolution, on 332 levels, from the surface up to a limit of approximately 40 km. The grid extends that high so the full height of soundings can be captured; however, most soundings terminate at an altitude between 25 and 30 km, above which no data is provided. Between soundings, the VAP linearly interpolates atmospheric state variables in time for each height level. In addition, INTERPOLATEDSONDE provides relative humidity scaled to microwave radiometer (MWR) observations.The INTERPOLATEDSONDE VAP, a continuous time-height grid of relative humidity-corrected sounding data, is intended to provide input to higher-order products, such as the Merged Soundings (MERGESONDE; Troyan 2012) VAP, which extends INTERPOLATEDSONDE by incorporating model data. The INTERPOLATEDSONDE VAP also is used to correct gaseous attenuation of radar reflectivity in products such as the KAZRCOR VAP.

54 ENVIRONMENTAL SCIENCES↗

Collaborative Project: The problem of bias in defining uncertainty in computationally enabled strategies for data-driven climate model development. Final Technical Report.

The objective of the project is to develop strategies for better representing scientific sensibilities within statistical measures of model skill that then can be used within a Bayesian statistical framework for data-driven climate model development and improved measures of model scientific uncertainty. One of the thorny issues in model evaluation is quantifying the effect of biases on climate projections. While any bias is not desirable, only those biases that affect feedbacks affect scatter in climate projections. The effort at the University of Texas is to analyze previously calculated ensembles of CAM3.1 with perturbed parameters to discover how biases affect projections of global warming. The hypothesis is that compensating errors in the control model can be identified by their effect on a combination of processes and that developing metrics that are sensitive to dependencies among state variables would provide a way to select version of climate models that may reduce scatter in climate projections. Gabriel Huerta at the University of New Mexico is responsible for developing statistical methods for evaluating these field dependencies. The UT effort will incorporate these developments into MECS, which is a set of python scripts being developed at the University of Texas for managing the workflow associated with data-driven climate model development over HPC resources. This report reflects the main activities at the University of New Mexico where the PI (Huerta) and the Postdocs (Nosedal, Hattab and Karki) worked on the project.

58 GEOSCIENCES↗

Interpolated Sounding and Gridded Sounding Value-Added Products

Standard Atmospheric Radiation Measurement (ARM) Climate Research Facility sounding files provide atmospheric state data in one dimension of increasing time and height per sonde launch. Many applications require a quick estimate of the atmospheric state at higher time resolution. The INTERPOLATEDSONDE (i.e., Interpolated Sounding) Value-Added Product (VAP) transforms sounding data into continuous daily files on a fixed time-height grid, at 1-minute time resolution, on 332 levels, from the surface up to a limit of approximately 40 km. The grid extends that high so the full height of soundings can be captured; however, most soundings terminate at an altitude between 25 and 30 km, above which no data is provided. Between soundings, the VAP linearly interpolates atmospheric state variables in time for each height level. In addition, INTERPOLATEDSONDE provides relative humidity scaled to microwave radiometer (MWR) observations.

54 ENVIRONMENTAL SCIENCES↗

Integrated Computational Materials Engineering Development of Advanced High Strength Steel for Lightweight Vehicles

The goal of the ICME 3GAHSS project was to successfully demonstrate the applicability of Integrated Computational Materials Engineering (ICME) for the development and deployment of third generation advanced high strength steels (3GAHSS) for immediate weight reduction in passenger vehicles. The ICME approach integrated results from well-established computational and experimental methodologies to develop a suite of material constitutive models (deformation and failure), manufacturing process and performance simulation modules, a properties database, as well as the computational environment linking them together for both performance prediction and material optimization. This is the Final Report for the ICME 3GAHSS project, which achieved the fol-lowing objectives: 1) Developed a 3GAHSS ICME model, which includes atomistic, crystal plasticity, state variable and forming models. The 3GAHSS model was implemented in commercially available LS-DYNA and a user guide was developed to facilitate use of the model. 2) Developed and produced two 3GAHSS alloys using two different chemistries and manufacturing processes, for use in calibrating and validating the 3GAHSS ICME Model. 3) Optimized the design of an automotive subassembly by substituting 3GAHSS for AHSS yielding a design that met or exceeded all baseline performance requirements with a 30% mass savings. A technical cost model was also developed to estimate the cost per pound of weight saved when substituting 3GAHSS for AHSS. The project demonstrated the potential for 3GAHSS to achieve up to 30% weight savings in an automotive structure at a cost penalty of up to $0.32 to $1.26 per pound of weight saved. The 3GAHSS ICME Model enables the user to design 3GAHSS to desired mechanical properties in terms of strength and ductility.

36 MATERIALS SCIENCE↗

Vertically Resolved Retrievals of Aerosol Concentrations and Effective Radii from the DOE Combined HSRL and Raman lidar Measurement Study (CHARMS) Merged High-Spectral-Resolution Lidar-Raman Lidar Data Set

The Department of Energy’s (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility is a key component of DOE’s research strategy to address global climate change. The objective of the ARM Facility is to provide an experimental testbed for studying important atmospheric effects, particularly cloud and aerosol processes, and testing parameterizations of these processes for use in climate change models. ARM observational sites include a broad array of instruments to characterize cloud and aerosol properties, atmospheric state variables, and incoming and outgoing radiation. Advanced lidar systems provide valuable and unique information about the atmospheric column. We examine the potential for providing vertically resolved information about aerosol size and concentration with a combination of two such lidar systems, the Raman lidar and the High-Spectral-Resolution Lidar (HSRL). Currently ARM has a limited ability to remotely measure profiles of aerosol optical and microphysical properties on an operational basis. Although the Southern Great Plains (SGP), Eastern North Atlantic (ENA), and North Slope of Alaska (NSA, Oliktok Point location only) Raman lidars provide profiles of aerosol backscatter and extinction at 355 nm, additional information is required to better characterize the impacts of aerosols on radiation and clouds. Recent advances in lidar retrieval theory and algorithm development (Müller et al. 1999, 2014; Chemyakin et al. 2014) demonstrate that multi-wavelength lidar retrieval algorithms that use measurements of aerosol extinction at both 355 and 532 nm and backscatter at 355, 532, and 1064 nm have the potential to constrain both the aerosol optical (e.g., complex index of refraction, scattering, etc.) and microphysical properties (e.g., effective radius, concentration). This advancement greatly increases the utility of the remote-sensed aerosol observations. Based on this work, the HSRL group at the National Aeronautics and Space Administration (NASA)’s Langley Research Center (LaRC) developed automated algorithms for retrieving aerosol optical and microphysical properties, demonstrated these retrievals using data from the unique NASA/LaRC airborne multi-wavelength HSRL-2 system, and validated these retrievals using coincident airborne in situ data (Müller et al. 2014; Sawamura et al. 2017). We report here on the use of ARM instrumentation to gather a multi-wavelength ground-based lidar data set to test such retrievals. Continuous (24/7) operation of the co-located SGP Raman lidar and University of Wisconsin (UW) HSRL during the 10-week Combined HSRL and Raman lidar Measurement Study (CHARMS) (mid-July through September 2015) allowed the acquisition of a unique, multi-wavelength ground-based lidar data set. The ARM SGP Raman lidar (RLID) measured profiles of aerosol backscatter, extinction, and depolarization at 355 nm and profiles of water vapor mixing ratio and temperature. The UW HSRL simultaneously measured profiles of aerosol backscatter, extinction and depolarization at 532 nm and aerosol backscatter at 1064 nm. This RLID+HSRL combination provided a test data set to investigate these more advanced aerosol retrieval techniques with the potential to greatly improve aerosol characterization in the lower atmosphere.

54 ENVIRONMENTAL SCIENCES↗

FY2017 Updates to the SAS4A/SASSYS-1 Safety Analysis Code

The SAS4A/SASSYS-1 safety analysis software is used to perform deterministic analysis of anticipated events as well as design-basis and beyond-design-basis accidents for advanced fast reactors. It plays a central role in the analysis of U.S. DOE conceptual designs, proposed test and demonstration reactors, and in domestic and international collaborations. This report summarizes the code development activities that have taken place during FY2017. Extensions to the void and cladding reactivity feedback models have been implemented, and Control System capabilities have been improved through a new virtual data acquisition system for plant state variables and an additional Block Signal for a variable lag compensator to represent reactivity feedback for novel shutdown devices. Current code development and maintenance needs are also summarized in three key areas: software quality assurance, modeling improvements, and maintenance of related tools. With ongoing support, SAS4A/SASSYS-1 can continue to fulfill its growing role in fast reactor safety analysis and help solidify DOE’s leadership role in fast reactor safety both domestically and in international collaborations.

97 MATHEMATICS AND COMPUTING↗

A time-parallel method for the solution of PDE-constrained optimization problems

We study a time-parallel approach to solving quadratic optimization problems with linear time-dependent partial differential equation (PDE) constraints. These problems arise in formulations of optimal control, optimal design and inverse problems that are governed by parabolic PDE models. They may also arise as subproblems in algorithms for the solution of optimization problems with nonlinear time-dependent PDE constraints, e.g., in sequential quadratic programming methods. We apply a piecewise linear finite element discretization in space to the PDE constraint, followed by the Crank-Nicolson discretization in time. The objective function is discretized using finite elements in space and the trapezoidal rule in time. At this point in the discretization, auxiliary state variables are introduced at each discrete time interval, with the goal to enable: (i) a decoupling in time; and (ii) a fixed-point iteration to recover the solution of the discrete optimality system. The fixed-point iterative schemes can be used either as preconditioners for Krylov subspace methods or as smoothers for multigrid (in time) schemes. We present promising numerical results for both use cases.

97 MATHEMATICS AND COMPUTING↗

Numerical Integration of Viscoelastic Models

Accurate modeling of viscoelasticity remains an important consideration for a variety of materials (e.g. polymers and inorganic glasses). As such, over the previous decades a substantial body of work has been dedicated to developing appropriate constitutive models for viscoelasticity ranging from initial considerations of linear thermoviscoelasticity to more complex non-linear formulations incorporating fictive temperatures or potential energy clocks including the use of both internal state variable(ISV) and hereditary integral representations. Nonetheless, relatively limited (in comparison to plasticity) attention has been paid to the numerical integration of such schemes. In terms of integral based formulations, Taylor et al. first considered the problem of the integration of a linear viscoelasticity model. That work focused on the integration of the hereditary integrals and demonstrated improved performance of the new scheme with a custom finite element code over an existing finite difference reference. Chambers and Becker, using a free volume based shift factor, also considered the integration of the hereditary integrals and the impact on the problem of a pressurized thick-walled cylinder and developed an adaptive scheme to bound the error. Chambers later developed three-point Gauss and composite integration schemes for the hereditary integrals and noted improved accuracy. With respect to ISV-based schemes, formulations for the non-linear Schapery model have been proposed. However, in those efforts greater attention was paid to convergence of the non-linear solution scheme than impact of numerical integration. Various authors (e.g. Holzapfel and Simo and Hughes) have also studied the use of convolution integrals with differential forms of ISVs for temperature-independent formulations. Regardless, while the "potential energy clock" (PEC) and "simplified potential energy clock"(SPEC) models have been used to study a variety of non-linear responses (e.g.), limited attention has been paid to the numerical performance. As will be discussed later, the "clock" at the center of the formulations includes temperature and complex history dependence making the numerical integration of such a model even more challenging. Thus, in the current work an initial effort towards characterizing the numerical integration of the constitutive model through simplified problems is performed. To that end, in Section 2 the theory of the model is briefly presented while the numerical integration is discussed in Section 3. Results of various studies characterizing the numerical behavior and performance are then given in Section 4. Finally, some concluding remarks and thoughts for follow on works are provided in Section 5.

36 MATERIALS SCIENCE↗

Eulerian Applications Project - xRage: High Explosives Package [Slides]

At the most elementary level, the high explosive model in xRage (or many other hydro codes) computes the amount of material that changes from a high explosive reactant to a high explosive product during a time step and accounts for the energy added to the flow by this reaction. The actual chemical reactions that occur during this conversion are not modeled. These processes are very complicated and involve levels of spatial and temporal resolution that are generally not practical for an application code. They would also require a state variable representation that includes all of the chemical species that are produced by the reaction and account for multiple stage chemical reactions. Instead the explosive models we generally use are phenomenological and try to represent the most important physics that drive the explosive wave in a very simplified way. These models generally have a number of tunable parameters that must be calibrated to experiment.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Development of a Multi-physics Capability for Predicting Residual Stress in a GTS Reservoir.

This report documents completion of a Level 2 Milestone on the development of a multi-physics capability to predict the evolving material state through the manufacturing process of a Gas Transfer Systems (GTS) reservoir. We present details on new developments and capability improvements that address the following completion criteria: (i) validation of a microstructure evolution model, including recrystallization and strain aging, (ii) demonstration of the capability to remesh, map and transfer material state (internal state variables) and residual stress from forging to machining to welding processes, and (iii) formal V&V characterization and quantification of uncertainties of material parameters and manufacturing process parameters on residual stress.

36 MATERIALS SCIENCE↗

Modular Linear Thermoviscoelastic Model

Time-dependent, viscoelastic responses of materials like polymers and glasses have long been studied. As such, a variety of models have been put forth to describe the behavior including simple rheological models (e.g. Maxwell, Kelvin), linear "fading memory" theories, and hereditary integral based linear thermal viscoelastic approaches as well as more recent nonlinear theories that are either integral, fictive temperature, or differential internal state variable based. The current work details a new LINEAR_THERMOVISCOELASTIC model that has been added to LAME. This formulation represents a viscoelastic theory that neglects some of the phenomenological details of the PEC/SPEC models in favor of efficiency and simplicity. Furthermore, this new model is a first step towards developing "modular" viscoelastic capabilities akin to those available with hardening descriptions for plasticity models in LAME. Specifically, multiple different (including user-defined) shift-factor forms are implemented with each being easily selected via parameter specification rather than requiring distinct material models.

36 MATERIALS SCIENCE↗

Simple Algorithms for a Physics Code Base

Project 1: Developing a container for state variables in a physics code. Project 2: Setting up near-continuous performance and usage metrics monitoring of a physics code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Navigation and Meteorological Data from Multiple Sensors on Airborne Platform (NAVMET-AIR) Value-Added Product Report

The Navigation and Meteorological Data from Multiple Sensors on Airborne Platform (NAVMET-AIR) value-added product encompasses the aafnaviwg data set in the ARM Data Center (ADC), known as the IWG file, named for the Inter-Agency Working Group for Airborne Data and Telemetry Systems (IWGADTS). Its purpose is to produce a suite of tools to promote standardization of instrument interface, data format, and data processing. This data set’s contents are navigational and meteorological state variables at 1 Hz. It also contains higher-order data flags created by ARM Aerial Facility (AAF) scientists to aid analysis such as periods when the aircraft is flying level, operating in maneuvers, or flying through cloud. Due to the nature of research flights, the payload of the aircraft consists of many duplicate and overlapping instruments to ensure, in case of instrument failure, there is a backup to record the data. Measurements from multiple instruments are consolidated into a single file. Each variable of this data set is carefully chosen from the onboard instrumentation and quality checked by AAF scientists to create this wholistic and most accurate data set for airborne research. This document is intended to be a hub towards the individual “read me” files produced for each campaign. For reference about what instrument was used for calculations on specific days, the read me files are located in the ARM IOP archive at the addresses listed, or access the IOP (intensive operational period) database at https://www.archive.arm.gov/.

54 ENVIRONMENTAL SCIENCES↗

Facilitating better and faster simulations of aerosol-cloud interactions in Earth system models

Focal Area(s): 1. Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising a hierarchy of models. 2. Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: One major challenge that Earth system models (ESMs) face in providing credible prediction of the Earth system and its water cycle characteristics (e.g., mean state, variability, and extreme events) is to accurately simulate aerosol-cloud interactions (ACI). The physical, chemical, and dynamical processes affecting ACI are extremely complex and they range from nanoscale to planetary scale. In each model development cycle, scientists spend significant efforts investigating model deficiencies and uncertainties associated with aerosols (e.g., emissions, chemical processes, aerosol microphysics, and transport) and clouds (e.g., macrophysics, microphysics, turbulence, and large-scale circulation) in order to develop improved treatments. However, despite decades of active research, ACI is still a major source of uncertainty in climate projections, even though great progress has been made. Specific scientific challenges include: (i) Parameterizations are developed based on limited data; (ii) The complexity of a parameterization required for accurate predictions is not understood; (iii) Incomplete and unknown physics leads to errors in the fully coupled Earth system; and (iv) Complex physics is computationally too expensive to employ in ESMs.

54 ENVIRONMENTAL SCIENCES↗

Thermomechanical Modeling of Copper

The plastic deformation of metals is a dissipative process. Some fraction of the plastic work is converted to heat which, given the temperature dependent response of metals, produces a thermal-mechanical coupling. In various cases, for instance when the loading is dynamic, this interaction can impact the resulting response of a material and/or system. Thus, appropriately capturing the heat generation from plastic work is necessary for various solid mechanics analysis. Determination of the fraction of work converted to heat has been long studied. Recent developments have demonstrated that the fraction is not constant but depends on various state variables. Resolving these features requires combined modeling and experimental studies. To this end, 304L stainless steel – a poor thermal conductor – was recently subjected to such an investigation. Advanced modeling capabilities were deployed to assess novel thermomechanically coupled experiments. As a complement to that study, in the current work a similar investigation is performed on copper – a good thermal conductor – to assess performance on the opposite end of the spectrum. The current document discusses these modeling efforts.

36 MATERIALS SCIENCE↗