Search NASA⌕ Search

SEARCH · Search NASA

Results for “model structural error”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Galaxy–galaxy lensing estimators and their covariance properties

Here, we study the covariance properties of real space correlation function estimators – primarily galaxy–shear correlations, or galaxy–galaxy lensing – using SDSS data for both shear catalogues and lenses (specifically the BOSS LOWZ sample). Using mock catalogues of lenses and sources, we disentangle the various contributions to the covariance matrix and compare them with a simple analytical model. We show that not subtracting the lensing measurement around random points from the measurement around the lens sample is equivalent to performing the measurement using the lens density field instead of the lens overdensity field. While the measurement using the lens density field is unbiased (in the absence of systematics), its error is significantly larger due to an additional term in the covariance. Therefore, this subtraction should be performed regardless of its beneficial effects on systematics. Comparing the error estimates from data and mocks for estimators that involve the overdensity, we find that the errors are dominated by the shape noise and lens clustering, which empirically estimated covariances (jackknife and standard deviation across mocks) that are consistent with theoretical estimates, and that both the connected parts of the four-point function and the supersample covariance can be neglected for the current levels of noise. While the trade-off between different terms in the covariance depends on the survey configuration (area, source number density), the diagnostics that we use in this work should be useful for future works to test their empirically determined covariances.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part I: Scheme Description

A novel framework is proposed for the bulk parameterization of rain microphysics: the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS). It is designed to facilitate direct constraint by observations using Bayesian inference. BOSS combines existing process-level microphysical knowledge with flexible process rate formulations and parameters constrained by observations within a Bayesian framework. Furthermore, using a raindrop size distribution (DSD) normalization method that relates DSD moments to one another via generalized power series, generalized multivariate power expressions are derived for the microphysical process rates as functions of a set of prognostic DSD moments. The scheme is flexible and can utilize any number and combination of prognostic moments and any number of terms in the process rate formulations. This means that both uncertainty in parameter values and structural uncertainty associated with the process rate formulations can be investigated systematically, which is not possible using traditional schemes. In this paper, BOSS is compared to two- and three-moment versions of a traditional bulk rain microphysics scheme (denoted as MORR). It is demonstrated that some process formulations in MORR are analytically equivalent to the generalized power expressions in BOSS using one or two terms, while others are not. BOSS is able to replicate the behavior of MORR in idealized one-dimensional rainshaft tests, but with a much more flexible and systematic design. Part II of this study describes the application of BOSS to derive rain microphysical process rates and posterior parameter distributions in Bayesian experiments using Markov chain Monte Carlo sampling constrained by synthetic observations.

58 GEOSCIENCES↗
DOE OSTI2019

Hierarchical Modeling and Statistical Calibration for Photometric Redshifts

The cosmological exploitation of modern photometric galaxy surveys requires both accurate (unbiased) and precise (narrow) redshift probability distributions derived from broadband photometry. Existing methodologies do not meet these requirements. Standard template fitting delivers interpretable models and errors, but lacks the flexibility to learn inaccuracies in the observed photometry or spectral templates. Machine learning addresses those issues, but requires representative training data, and the resulting models and uncertainties cannot be interpreted in the context of a physical model or outside of the training data. We present a hierarchical modeling approach simultaneously addressing the issues of flexibility, interpretability, and generalization. It combines template fitting with flexible (machine-learning-like) models to correct the spectral templates, model their redshift distributions, and recalibrate the photometric observations. By optimizing the full posterior distribution of the model and solving for its (thousands of) parameters, one can perform a global statistical calibration of the data and the spectral energy distribution (SED) model. We apply this approach to the public Dark Energy Survey Science Verification data and show that it provides more accurate and compact redshift posterior distributions than existing methods, as well as insights into residual photometric and SED systematics. Here, the model is causal and makes predictions for future data (e.g., additional photometric bandpasses), and its internal parameters and components are interpretable. This approach does not formally require the training data to be complete or representative; in principle, it can even work in regimes in which few or no spectroscopic redshifts are available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Getting allometry right at the Oak Ridge free‐air CO 2 enrichment experiment: Old problems and new opportunities for global change experiments

Societal Impact Statement Free‐air CO 2 enrichment (FACE) experiments provide essential data on forest responses to increasing atmospheric CO 2 for evaluations of climate change impacts on humanity. Understanding and reducing the uncertainty in the experimental results is critical to ensure scientific and public confidence in the models and policy initiatives that derive therefrom. One source of uncertainty is the estimation of tree biomass using mathematical relationships between biomass and easily obtained and non‐destructive measurements (allometry). We evaluated the robustness of the allometric relationships established at the beginning of a FACE experiment and discuss the challenges and opportunities for the new generation of FACE experiments. Summary Long‐term field experiments to elucidate forest responses to rising atmospheric CO 2 concentration require allometric equations to estimate tree biomass from non‐destructive measurements of tree size. We analyzed whether the allometric equations established at the beginning of a free‐air CO 2 enrichment (FACE) experiment in a Liquidambar styraciflua plantation were still valid at the end of the 12 year experiment. Aboveground woody biomass was initially predicted by an equation that included bole diameter, taper, and height, assuming that including taper and height as predictors would accommodate changes in tree structure that might occur over time and in response to elevated CO 2 . At the conclusion of the FACE experiment, we harvested 23 trees, measured dimensions and dry mass of boles and branches, and extracted and measured the woody root mass of 10 trees. Although 10 of the harvested trees were larger than the trees used to establish the allometric relationship, measured aboveground woody biomass was well predicted by the original allometry. The initial linear equation between bole basal area and woody root biomass underestimated final root biomass by 28%, but root biomass was just 21% of total wood mass, and errors in aboveground and belowground estimates were offsetting. The allometry established at the beginning of the experiment provided valid predictions of tree biomass throughout the experiment. New allometric approaches using terrestrial laser scanning should reduce an important source of uncertainty in decade‐long forest experiments and in assessments of centuries‐long forest biomass accretion used in evaluating carbon offsets and climate mitigation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Multiresolution molecular mechanics: Surface effects in nanoscale materials

Surface effects have been observed to contribute significantly to the mechanical response of nanoscale structures. The newly proposed energy-based coarse-grained atomistic method Multiresolution Molecular Mechanics (MMM) (Yang, To (2015), ) is applied to capture surface effect for nanosized structures by designing a surface summation rule SR{sup S} within the framework of MMM. Combined with previously proposed bulk summation rule SR{sup B}, the MMM summation rule SR{sup MMM} is completed. SR{sup S} and SR{sup B} are consistently formed within SR{sup MMM} for general finite element shape functions. Analogous to quadrature rules in finite element method (FEM), the key idea to the good performance of SR{sup MMM} lies in that the order or distribution of energy for coarse-grained atomistic model is mathematically derived such that the number, position and weight of quadrature-type (sampling) atoms can be determined. Mathematically, the derived energy distribution of surface area is different from that of bulk region. Physically, the difference is due to the fact that surface atoms lack neighboring bonding. As such, SR{sup S} and SR{sup B} are employed for surface and bulk domains, respectively. Two- and three-dimensional numerical examples using the respective 4-node bilinear quadrilateral, 8-node quadratic quadrilateral and 8-node hexahedral meshes are employed to verify and validate the proposed approach. It is shown that MMM with SR{sup MMM} accurately captures corner, edge and surface effects with less 0.3% degrees of freedom of the original atomistic system, compared against full atomistic simulation. The effectiveness of SR{sup MMM} with respect to high order element is also demonstrated by employing the 8-node quadratic quadrilateral to solve a beam bending problem considering surface effect. In addition, the introduced sampling error with SR{sup MMM} that is analogous to numerical integration error with quadrature rule in FEM is very small. - Highlights: • Surface effect captured by Multiresolution Molecular Mechanics (MMM) is presented. • A novel surface summation rule within the framework of MMM is proposed. • Surface, corner and edges effects are accuterly captured in two and three dimension. • MMM with less 0.3% degrees of freedom of atomistics reproduces atomistic results.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Isovalent Alloying and Heterovalent Substititution as Routes to Accelerate the Development and Optimization of Super-Efficient Halide Perovskite PV Solar Cells

Halide perovskites represent one of the fastest growing photovoltaic (PV) device technologies in the last decade, with research efficiencies among the highest confirmed. However, the research needed to advance tandem perovskite solar cells requires a specialized understanding base, but there is no such base for this young group of materials with the same depth as in semiconductor alloy research. Whereas IV-IV, II-V and II-VI semiconductors have a sophisticated science foundation refined over decades due of their use in electronics, lasers, LEDs and night vision applications, this is not yet the case for halide perovskites. At the same time, these materials behave rather differently than conventional semiconductors. Within the rapidly unfolding field of PV, the emerging halide-perovskite alloy research has largely been based on empiricism. Clearly, a concerted and comprehensive application to halide perovskite alloys of the modern quantum theory of solids brought us over the years a broad understanding of the critical factors enabling conventional inorganic PV could leapfrog the Edisonian, trial-and-error search for solutions to critical problems that are poorly understood. These theory approaches are required to accelerate the progress toward commercially ready hybrid-perovskites by providing a theory-guided experiment component to the EERE investments in this high-priority R&D. During this project we have studied three different rungs on the ladder of knowledge about hybrid perovskites. This kind of single halide perovskites represent the first level of development (rung 1) of solar cells, and includes materials such as APbI 3 , ASnI 3 , APbBr 3 , ASnBr 3 with different A cations, such as methylammonium (MA) with formamidinium (FM). Rung 2 is obtained by mixing B cations (Sn and Pb), and X anions (Cl and Br or I) producing a disordered alloy (A, A’)(M,M’)(X,X’) 3 -space. The most advanced rung corresponds to heterovalent-substitution of halide perovskites, creating an ordered double-perovskite compound. Although the focus of the project was on Rungs 2 and 3, we have made important discoveries regarding the nature of pure perovskites. We’ve observed that cubic perovskites should be modeled by a polymorphous network, instead of the basic single unit cell approach. This leads to the conclusion that perovskites have a fluid-like nature. An article regarding this work has been submitted to Nature Materials. Using the approach described above, we have also simulated alloys of halide perovskites, representing rung 2. We have focused on mixtures of [Cs,FA][Sn,Pb]I 3 with different concentrations, being able to disentangle the structural and electronic factors that determines the variation (bowing) of properties such as band gaps, mixing enthalpies and decomposition energies. For rung 3 compounds, we’ve closely interacted with experimental groups clarifying the chemical routes to obtain a certain compound, and what should be the most stable structure for them. As the final product, and consequent dissemination of the research, we’ve submitted two papers for publication, that are currently under evaluation. The interaction with experimental groups was also an important focus. We have participated in several local, national and international meetings, and reached out more closely to groups from CU Boulder and Arkansas.

14 SOLAR ENERGY↗
DOE OSTI2020

Wave Prediction using X-Band Radar (Final Summary Report)

Model Predictive Control (MPC) is the best controls framework available to wave energy today. It allows for the maximization of energy yield from a WEC system while respecting various constraints and losses in the system. This type of constrained optimization is a critical capability for the techno-economic optimization of WEC devices and is difficult to do effectively with causal controls frameworks. MPC requires a future prediction of wave excitation forces, requiring wave prediction on a future time horizon of up to 30 seconds. The future horizon required is highly dependent on the WEC and PTO type available. To our knowledge no-one has been able to implement MPC on a WEC device at sea, due to the fact that phase-resolved wave-prediction is not a capability that has been sufficiently developed to date. The key objectives of this project are focused on developing “industry-ready” wave prediction technology building blocks that can be applied to a wide range of different wave energy conversion topologies. Our collaborative work with device developers has demonstrated that the level of improvement attainable for any particular device is heavily dependent on the device and PTO topology chosen. The range of annual average power capture improvements is on the order of 10% to over 100%. However, the important aspect is that we can get within about 10% of theoretical upper limits even when considering errors in the wave prediction and more importantly, we can optimize performance while respecting various constraints such as motion amplitude or peak structural loads. The combination of improved performance while limiting structural loads has a net effect of reducing the levelized cost of energy from these emerging technologies and allows us to design control laws that meet an economic optimum. Efforts under a previous project on controls has focused on developing a set of wave-prediction algorithms that work by leveraging a network of wave sensing measurement buoys that are equipped with real-time telemetry. The present work is focusing primarily on using X-band radar to measure the wave-field around the WEC to predict the future excitation forces on the structure and combining the measurements of these different sources to create improvements in the wave-prediction accuracy.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Assessing the Performance of Exchange‐Correlation Density Functionals in Describing the Iron‐Catalyzed Ammonia Synthesis System

Density functional theory (DFT) has been widely employed for elucidating mechanistic aspects of heterogeneous catalysis. However, the accuracy of DFT calculations relies heavily on selecting an exchange-correlation (XC) functional that correctly describes the electronic structure of materials involved in the reactions. This study assesses the accuracy of several XC density functionals for modeling the iron-catalyzed ammonia synthesis system. In the assessment of functional accuracy, experimental references are compared to DFT-calculated values for the formation energy of gas-phase ammonia and nitrogen, bulk Fe/Fe-nitride (γ′-Fe 4 N) lattice constants and cohesive/formation energies, and nitrogen and ammonia binding energies on Fe(100), Fe(111), Fe(110), and γ′-Fe 4 N(111). It is observed that the experimental value for each of these descriptors is accurately modeled by at least one functional. RPBE alone provides reliable estimates for both the lattice constant and cohesive energy of Fe and γ′-Fe 4 N. Temperature-programmed desorption experiments led to estimates for N and NH 3 adsorption across several Fe-based facets that are best captured by RPBE. These results highlight the importance of choosing an appropriate XC functional that accurately describes experimental systems and offer insights into effectively modeling the interactions between nitrogen and ammonia on Fe-based surfaces.

adsorption↗
DOE OSTI2016

Quantum mechanics capacitance molecular mechanics modeling of core-electron binding energies of methanol and methyl nitrite on Ag(111) surface

We study a newly devised quantum mechanics capacitance molecular mechanics (QMCMM) method for the calculation of core-electron binding energies in the case of molecules adsorbed on metal surfaces. This yet untested methodology is applied to systems with monolayer of methanol/methyl nitrite on an Ag(111) surface at 100 K temperature. It was found out that the studied C, N, and O 1s core-hole energies converge very slowly as a function of the radius of the metallic cluster, which was ascribed to build up of positive charge on the edge of the Ag slab. Further analysis revealed that an extrapolation process can be used to obtain binding energies that deviated less than 0.5 eV against experiments, except in the case of methanol O 1s where the difference was as large as 1.8 eV. Additional QM-cluster calculations suggest that the latter error can be connected to the lack of charge transfer over the QM–CMM boundary. Thus, the results indicate that the QMCMM and QM-cluster methods can complement each other in a holistic picture of molecule-adsorbate core-ionization studies, where all types of intermolecular interactions are considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Constraining nuclear data via cosmological observations: Neutrino energy transport and big bang nucleosynthesis

Here, we introduce a new computational capability that moves toward a self-consistent calculation of neutrino transport and nuclear reactions for big bang nucleosynthesis (BBN). Such a self-consistent approach is needed to be able to extract detailed information about nuclear reactions and physics beyond the standard model from precision cosmological observations of primordial nuclides and the cosmic microwave background radiation. We also calculate the evolution of the early universe through the epochs of weak decoupling, weak freeze-out and big bang nucleosynthesis (BBN) by simultaneously coupling a full strong, electromagnetic, and weak nuclear reaction network with a multi-energy group Boltzmann neutrino energy transport scheme. The modular structure of our approach allows the dissection of the relative contributions of each process responsible for evolving the dynamics of the early universe. Such an approach allows a detailed account of the evolution of the active neutrino energy distribution functions alongside and self-consistently with the nuclear reactions and entropy/heat generation and flow between the neutrino and photon/electron/positron/baryon plasma components. Our calculations reveal nonlinear feedback in the time evolution of neutrino distribution functions and plasma thermodynamic conditions. We discuss the time development of neutrino spectral distortions and concomitant entropy production and extraction from the plasma. These effects result in changes in the computed values of the BBN deuterium and helium-4 yields that are on the order of a half-percent relative to a baseline standard BBN calculation with no neutrino transport. This is an order of magnitude larger effect than in previous estimates. For particular implementations of quantum corrections in plasma thermodynamics, our calculations show a 0.4% increase in deuterium and a 0.6% decrease in 4 He over our baseline. The magnitude of these changes are on the order of uncertainties in the nuclear physics for the case of deuterium and are potentially significant for the error budget of helium in upcoming cosmological observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Development of a new control rod drive mechanism design for the ISU AGN-201M reactor

The Aerojet General Nucleonics (AGN) model 201-Modified, known as the AGN-201M reactor, plays an essential role in the educational and research activities at Idaho State University (ISU). The ISU AGN-201M's original Control Rod Drive Mechanism (CRDM) has been in operation for more than fifty years with no large-scale redesigns. The CRDM is required to eject the fuel rods within one second during a SCRAM event (also known as a 'reactor trip') and adjust the control rods' insertion speed, and keeps the rod insertion sequence correct. The existing control rod drive mechanisms meet these criteria but experience a few concerns due to aging. Concerns include complex maintenance and costly repairs for old electromechanical components, rod position feedback errors, and the impediment of the plate during a SCRAM due to binding of the lead screws of the existing mechanism. During a binding event, the drive mechanism becomes locked, preventing the control rod's magnetic plate from moving in or out under the reactor's normal and emergency operating conditions. Although, the binding has no effect on the ability of the rod to exit the core during the SCRAM. To counteract the concerns and issues with the current CRDM, a new design has been proposed using newer components and a simplified design. The new design utilizes more advanced electric and mechanical components that are commercially available. The new CRDM system is divided into four main aspects: (1) control rod movement design (motor, lead screw, guide rods), (2) control rod ejection (springs, electromagnet), (3) control rod position and feedback (position transducer, microswitches), and (4) material selection and structural analysis. The new design aims to reduce the overall complexity and probability of failure to improve the reactor's overall reliability. With proper material selection and improved structural design, the new drives are lighter with little to no change in structural integrity. The new control rod drive mechanism eliminates binding scenarios by using a single lead screw and implementing additional guide rods. An advanced linear position sensor and microswitches replace the existing and aging synchro system for accurate rod position feedback resulting in better reactivity control. The new design meets the reactor's operational limits by having an average reactivity insertion of 0.065% Δk/k per second, which corresponds to a total control rod insertion time of 19.23 s, while the control rod's ejection time remains less than one second during a SCRAM event. The new design ensures the reactor's long-term viability for educational and research activities by increasing the reliability and safety of operation for years to come.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Improving Dark Energy Constraints Using Low-Redshift Large-Scale Structures

The primary goal of this project was to improve constraints on dark energy measurements by improving our ability to extract cosmological information from low redshift large-scale structures. PI Clowe's project's primary aim was to reduce the bias in measurements of the masses of clusters of galaxies to a level where the evolution of the cluster mass function can be used in the Vera Rubin Observatory's Legacy Survey of Space and Time Dark Energy Science Collaboration survey to improve the accuracy of the measurement of dark energy and other cosmological parameters. Co-PI Seo's project studied observational systematics affecting large-scale clustering of galaxies, which will be used to improve dark energy constraints from the Dark Energy Spectroscopic Instrument (DESI). The cluster lensing project employed a series of simulations and observations of clusters of galaxies to test numerous potential systematic errors in cluster mass measurements using weak gravitational lensing as the accuracy of current weak lensing measurements are more than order of magnitude worse than what is required to use clusters of galaxies for accurate determination of dark energy parameters. PI Clowe and group developed and analyzed simulations to test for, and correct biases introduced in, the weak lensing measurement process. Finally, PI Clowe and group developed a method of detecting clusters using galaxy overdensities and applied the method to the BLISS and DES surveys. The success of spectroscopic dark energy mission such as the extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the Dark Energy Spectroscopic Instrument (DESI) will depend on a thorough understanding of various observational systematics in the target density fluctuations that would give rise to spurious, non-cosmological signals. PI Seo and group developed a deep learning, artificial neural network (ANN) technique that modeled and mitigated such effects, aimed at deriving more robust galaxy clustering signals not only for the baryon acoustic oscillation feature and redshift-space distortions but also for primordial non-Gaussianity constraint.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

An optimal FFT-based anisotropic power spectrum estimator

Measurements of line-of-sight dependent clustering via the galaxy power spectrum's multipole moments constitute a powerful tool for testing theoretical models in large-scale structure. Recent work shows that this measurement, including a moving line-of-sight, can be accelerated using Fast Fourier Transforms (FFTs) by decomposing the Legendre polynomials into products of Cartesian vectors. Here, we present a faster, optimal means of using FFTs for this measurement. We avoid redundancy present in the Cartesian decomposition by using a spherical harmonic decomposition of the Legendre polynomials. With this method, a given multipole of order ℓ requires only 2ℓ+1 FFTs rather than the (ℓ+1)(ℓ+2)/2 FFTs of the Cartesian approach. For the hexadecapole (ℓ = 4), this translates to 40% fewer FFTs, with increased savings for higher ℓ. The reduction in wall-clock time enables the calculation of finely-binned wedges in P ( k ,μ), obtained by computing multipoles up to a large ℓ{sub max} and combining them. This transformation has a number of advantages. We demonstrate that by using non-uniform bins in μ, we can isolate plane-of-sky (angular) systematics to a narrow bin at 0μ ≅ while eliminating the contamination from all other bins. We also show that the covariance matrix of clustering wedges binned uniformly in μ becomes ill-conditioned when combining multipoles up to large values of ℓ{sub max}, but that the problem can be avoided with non-uniform binning. As an example, we present results using ℓ{sub max}=16, for which our procedure requires a factor of 3.4 fewer FFTs than the Cartesian method, while removing the first μ bin leads only to a 7% increase in statistical error on f σ{sub 8}, as compared to a 54% increase with ℓ{sub max}=4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Inverse design of photonic surfaces via multi fidelity ensemble framework and femtosecond laser processing

We demonstrate a multi-fidelity (MF) machine learning ensemble framework for the inverse design of photonic surfaces, trained on a dataset of 11,759 samples that we fabricate using high throughput femtosecond laser processing. The MF ensemble combines an initial low fidelity model for generating design solutions, with a high fidelity model that refines these solutions through local optimization. The combined MF ensemble can generate multiple disparate sets of laser-processing parameters that can each produce the same target input spectral emissivity with high accuracy (root mean squared errors < 2%). SHapley Additive exPlanations analysis shows transparent model interpretability of the complex relationship between laser parameters and spectral emissivity. Finally, the MF ensemble is experimentally validated by fabricating and evaluating photonic surface designs that it generates for improved efficiency energy harvesting devices. Our approach provides a powerful tool for advancing the inverse design of photonic surfaces in energy harvesting applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Using ARM Measurements to Understand and Reduce the Double ITCZ Biases in the Community Atmospheric Model

1. Understanding of the observed variability of ITCZ in the equatorial eastern Pacific. The annual mean precipitation in the eastern Pacific has a maximum zonal band north of the equator in the ITCZ where the maximum SST is located. During the boreal spring (referring to February, March, and April throughout the present paper), because of the accumulated solar radiation heating and oceanic heat transport, a secondary maximum of SST exists in the southeastern equatorial Pacific. Associated with this warm SST is also a seasonal transitional maximum of precipitation in the same region in boreal spring, exhibited as a weak double ITCZ pattern in the equatorial eastern Pacific. This climatological seasonal variation, however, varies greatly from year to year: double ITCZ in the boreal spring occurs in some years but not in other years; when there a single ITCZ, it can appear either north, south or at the equator. Understanding this observed variability is critical to find the ultimate cause of the double ITCZ in climate models. Seasonal variation of ITCZ south of the eastern equatorial Pacific: By analyzing data from satellites, field measurements and atmospheric reanalysis, we have found that in the region where spurious ITCZ in models occurs, there is a “seasonal cloud transition” — from stratocumulus to shallow cumulus and eventually to deep convection —in the South Equatorial Pacific (SEP) from September to April that is similar to the spatial cloud transition from the California coast to the equator. This seasonal transition is associated with increasing sea surface temperature (SST), decreasing lower tropospheric stability and large-scale subsidence. This finding of seasonal cloud transition points to the same source of model errors in the ITCZ simulations as in simulation of stratocumulus-cumulus-deep convection transition. It provides a test for climate models to simulate the relationships between clouds and large-scale atmospheric fields in a region that features a spurious double Inter-tropical Convergence Zone (ITCZ) in most models. This work is recently published in Yu et al. (2016). Interannual variation of ITCZ south of the eastern equatorial Pacific: By analyzing data from satellites, field measurements and atmospheric reanalysis, we have characterized the interannual variation of boreal spring precipitation in the eastern tropical Pacific and found the cause of the observed interannual variability. We have shown that ITCZ in this region can occur as a single ITCZ along the Equator, single ITCZ north of the Equator, single ITCZ south of the Equator, and double ITCZ on both sides of the Equator. We have found that convective instability only plays a secondary role in the ITCZ interannual variability. Instead, the remote impact of the Pacific basin-wide SST on the horizontal gradient of surface pressure and wind convergence is the primary driver of this interannual variability. Results point to the need to include moisture convergence in convection schemes to improve the simulation of precipitation in the eastern tropical Pacific. This result has been recently submitted for publication (Yu and Zhang 2016). 2. Improvement of model parameterizations to reduce the double ITCZ bias We analyzed the current status of climate model performance in simulating precipitation in the equatorial Pacific. We have found that the double ITCZ bias has not been reduced in CMIP5 models relative to CMIP4 models. We have characterized the dynamic structure of the common bias by using precipitation, sea surface temperature, surface winds and sea-level. Results are published in Zhang et al. (2015): Since cumulus convection plays a significant role in the double ITCZ behavior in models, we have used measurements from ARM and other sources to carry out a systematic analysis of the roles of shallow and deep convection in the CAM. We found that in both CAM4 and CAM5, when the intensity of deep convection decreases as a result of parameterization change, the intensity of shallow convection increases, leading to very different changes in precipitation partitions but little change in the total precipitation. The different precipitation partition however can manifest themselves in other measures of model performances including temperature and humidity. This study points to the need to treat model physical parameterizations as integrated system rather than individual components. Results from this study are published in Wang and Zhang (2013). Since shallow convection interacts with the deep convection scheme and surface turbulence to trigger the double ITCZ, we studied methods to improve the shallow convection scheme in climate models. We investigated the bulk budgets of the vertical velocity and its parameterization in convective cores, convective updrafts, and clouds by using large-eddy simulation (LES) of four shallow convection cases including one from ARM. We proposed optimal forms of the Simpson and Wiggert equation to calculate the vertical velocity in bulk mass flux convection schemes for convective cores, convective updrafts, and convective clouds as parameterization schemes. The new scheme is published in Wang and Zhang (2014). By using long-term radar-based ground measurements from ARM, we derived a scale-aware inhomogeneity parameterization of cloud liquid water in climate models. We found a relationship between the inhomogeneity parameter and the model grid size as well as atmospheric stability. This relationship is implemented in the CESM to describe the subgrid-scale cloud inhomogeneity. Relative to the default CESM with the finite-volume dynamic core at 2-degree resolution, the new parameterization leads to smaller cloud inhomogeneity and larger cloud liquid-water path in high latitudes, and the opposite effect in low latitudes, with the regional impact on shortwave cloud radiation effect of up to 10 W/m 2 . This is due to both the smaller (larger) grid size in high (low) latitudes in the longitude-latitude grid setting of CESM and the more stable (unstable) atmosphere. This parameterization is expected lead to more realistic simulation of tropical precipitation in high-resolution models. Results from this study are reported in Xie and Zhang (2015).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Object-Based Evaluation of Dynamical and Statistical Downscaled Precipitation Products over CONUS

High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.

Downscaling↗
DOE OSTI2020

Extraction of mechanical properties of materials through deep learning from instrumented indentation

Instrumented indentation has been developed and widely utilized as one of the most versatile and practical means of extracting mechanical properties of materials. This method is particularly desirable for those applications where it is difficult to experimentally determine the mechanical properties using stress–strain data obtained from coupon specimens. Such applications include material processing and manufacturing of small and large engineering components and structures involving the following: three-dimensional (3D) printing, thin-film and multilayered structures, and integrated manufacturing of materials for coupled mechanical and functional properties. Here, we utilize the latest developments in neural networks, including a multifidelity approach whereby deep-learning algorithms are trained to extract elastoplastic properties of metals and alloys from instrumented indentation results using multiple datasets for desired levels of improved accuracy. We have established algorithms for solving inverse problems by recourse to single, dual, and multiple indentation and demonstrate that these algorithms significantly outperform traditional brute force computations and function-fitting methods. Moreover, we present several multifidelity approaches specifically for solving the inverse indentation problem which 1) significantly reduce the number of high-fidelity datasets required to achieve a given level of accuracy, 2) utilize known physical and scaling laws to improve training efficiency and accuracy, and 3) integrate simulation and experimental data for training disparate datasets to learn and minimize systematic errors. The predictive capabilities and advantages of these multifidelity methods have been assessed by direct comparisons with experimental results for indentation for different commercial alloys, including two wrought aluminum alloys and several 3D printed titanium alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extraction of mechanical properties of materials through deep learning from instrumented indentation

Instrumented indentation has been developed and widely utilized as one of the most versatile and practical means of extracting mechanical properties of materials. This method is particularly desirable for those applications where it is difficult to experimentally determine the mechanical properties using stress–strain data obtained from coupon specimens. Such applications include material processing and manufacturing of small and large engineering components and structures involving the following: three dimensional (3D) printing, thin-film and multilayered structures, and integrated manufacturing of materials for coupled mechanical and functional properties. Here, we utilize the latest developments in neural networks, including a multifidelity approach whereby deep-learning algorithms are trained to extract elastoplastic properties of metals an d alloys from instrumented indentation results using multiple datasets for desired levels of improved accuracy. We have established algorithms for solving inverse problems by recourse to single, dual, and multiple indentation and demonstrate that these algorithms significantly outperform traditional brute force computations and function-fitting methods. Moreover, we present several multifidelity approaches specifically for solving the inverse indentation problem which 1) significantly reduce the number of high-fidelity datasets required to achieve a given level of accuracy, 2) utilize known physical and scaling laws to improve training efficiency and accuracy, and 3) integrate simulation and experimental data for training disparate datasets to learn and minimize systematic errors. The predictive capabilities and advantages of these multifidelity methods have been assessed by direct comparisons with experimental results for indentation for different commercial alloys, including two wrought aluminum alloys and several 3D printed titanium alloys.

3D printed materials, stress-strain behavior, mult↗