Search NASA⌕ Search

SEARCH · Search NASA

Results for “scaling methods”

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 289 records · Page 16

High precision tests of QCD without scale or scheme ambiguities: The 40th anniversary of the Brodsky–Lepage–Mackenzie method

A key issue in making precise predictions in QCD is the uncertainty in setting the renormalization scale μ r and thus determining the correct values of the QCD running coupling α s (μ r ) at each order in the perturbative expansion of a QCD observable. It has often been conventional to simply set the renormalization scale to the typical scale of the process Q and vary it in the range μ r $\in$ [Q/2, 2Q] in order to estimate the theoretical error. This is the practice of Conventional Scale Setting (CSS). The resulting CSS prediction will however depend on the theorist’s choice of renormalization scheme and the resulting pQCD series will diverge factorially. It will also disagree with renormalization scale setting used in QED and electroweak theory thus precluding grand unification. A solution to the renormalization scale-setting problem is offered by the Principle of Maximum Conformality (PMC), which provides a systematic way to eliminate the renormalization scale-and-scheme dependence in perturbative calculations. The PMC method has rigorous theoretical foundations, it satisfies Renormalization Group Invariance (RGI) and preserves all self-consistency conditions derived from the renormalization group. The PMC cancels the renormalon growth, reduces to the Gell-Mann–Low scheme in the N c → 0 Abelian limit and leads to scale- and scheme-invariant results. The PMC has now been successfully applied to many high-energy processes. In this article we summarize recent developments and results in solving the renormalization scale and scheme ambiguities in perturbative QCD. In particular, we present a recently developed method the PMC ∞ and its applications, comparing the results with CSS. The method preserves the property of renormalizable SU(N)/U(1) gauge theories defined as Intrinsic Conformality (iCF). This property underlies the scale invariance of physical observables and leads to a remarkably efficient method to solve the conventional renormalization scale ambiguity at every order in pQCD. This new method reflects the underlying conformal properties displayed by pQCD at NNLO, eliminates the scheme dependence of pQCD predictions and is consistent with the general properties of the PMC. A new method to identify conformal and β-terms, which can be applied either to numerical or to theoretical calculations is also shown. We present results for the thrust and C-parameter distributions in e + e - annihilation showing errors and comparison with the CSS. We also show results for a recent innovative comparison between the CSS and the PMC ∞ applied to the thrust distribution investigating both the QCD conformal window and the QED N c → 0 limit. In order to determine the thrust distribution along the entire renormalization group flow from the highest energies to zero energy, we consider the number of flavors near the upper boundary of the conformal window. In this flavor-number regime the theory develops a perturbative infrared interacting fixed point. These results show that PMC ∞ leads to higher precision and introduces new interesting features in the PMC. In fact, this method preserves with continuity the position of the peak, showing perfect agreement with the experimental data already at NNLO. We also show a detailed comparison of the PMC ∞ with the other PMC approaches: the multi-scale-setting approach (PMCm) and the single-scale-setting approach (PMCs) by comparing their predictions for three important fully integrated quantities R e+e- , R$_Τ$ and Γ (H→$b\bar{b}$) up to the four-loop accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Implicit-Explicit Multirate Infinitesimal GARK Methods

This work focuses on the development of a new class of high-order accurate methods for multirate time integration of systems of ordinary differential equations. Unlike other recent work in this area, the proposed methods support mixed implicit-explicit (IMEX) treatment of the slow time scale. In addition to allowing this slow time scale flexibility, the proposed methods utilize a so-called infinitesimal formulation for the fast time scale through definition of a sequence of modified “fast" initial-value problems that may be solved using any viable algorithm. We name the proposed class as implicit-explicit multirate infinitesimal generalized-structure additive Runge--Kutta (IMEX-MRI-GARK) methods. In addition to defining these methods, we prove that they may be viewed as specific instances of GARK methods and derive a set of order conditions on the IMEX-MRI-GARK coefficients to guarantee both third and fourth order accuracy for the overall multirate method. Additionally, we provide three specific IMEX-MRI-GARK methods, two of order three and one of order four. We conclude with numerical simulations on two multirate test problems, demonstrating the methods' predicted convergence rates and comparing their efficiency against both legacy IMEX multirate schemes and recent third and fourth order implicit MRI-GARK methods.

97 MATHEMATICS AND COMPUTING↗

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗

Photovoltaic inverter-based quantification of snow conditions and power loss

Snow is a significant challenge for photovoltaic (PV) systems at northern latitudes, where the pace of deployment is rapid but snow-related power losses can exceed 30% of annual production. Accurate snow-related power loss estimation methods for utility-scale sites can support snow mitigation strategies, inform resource planning and validate predictive snow-loss models. This study builds on our previous work on inverter-based detection of snow, and its implications for utility-scale power production, by validating the accuracy of our snow-loss method across different PV sites and system designs and highlighting its value in bringing greater visibility to PV plant operations in winter. Our estimation method is both novel and scalable, requiring only standard monitoring data to correlate snow-related losses with meteorological data. As demonstrated here, our validation method involved three main steps: 1) estimation of performance losses for multiple systems by comparing measured inverter data to modeled data; 2) application of a detection framework to identify which performance losses are snow-related; and 3) comparison of snow-related losses among three utility-scale sites differing in tilt angle. Results show that utility-scale systems at higher tilt angles consistently shed snow more quickly/completely than their lower-tilt counterparts. Further, monthly and seasonal snow losses are inversely and non-linearly correlated with tilt angle when normalized for cumulative snowfall. These results are consistent with the findings of previous studies and support the broad applicability of this method to fixed-tilt utility-scale PV systems around the world that routinely experience snow-related performance losses.

Cooper, Emma C. (ORCID:0000000190554098)↗

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system’s bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗

Advanced Physics Models for Particle-Particle Interactions (UCD final report)

The LDRD "Advanced Physics Models for Particle-Particle Interactions” collaborative research with Professor Jennifer Sinclair Curtis and PhD student Chanh Nguyen at the University of California, Davis involved meso-scale, discrete element method (DEM) modeling of micron-sized particles impacting a substrate. Using DEM, physical insight into the details of particle response upon impact can be obtained. Such details provide the type of information necessary to enhance existing continuum or macro-scale simulation tools so that modeling for material strength and particle breakage models is based on particle physics versus empirical relationships. While molecular dynamics simulations of atoms can also provide physical insight into the particle response to impact, the scale of particles that are treated is on the order of angstroms, not microns. This study gives evidence to the ability of the DEM simulation technique to provide such fundamental, detailed insight for micron-sized particles as a function of particle properties and impact velocity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Revealing Soil Organic Matter-mineral Associations with Advanced Chemical Imaging Methods

Although macro-scale in nature, most biogeochemical processes are driven by interactions at the microand nano-scale. Therefore, precise characterization of both organic and inorganic system components requires multi-scale, often multi-capability approaches. In this work, we present the state-of-the-art characterization of organic matter-mineral associations by high-resolution correlated chemical imaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scaling Building Energy Audits through Machine Learning Methods on Novel Drone Image Data

Building energy audits are time-consuming and labor-intensive. This paper describes a new method using machine learning (ML) techniques on novel data sources (drone images) to improve the identification of building characteristics and retrofit opportunities, and thereby reduce the effort for audits. The new ML method includes: (1) Building footprint extraction using line extraction, polygonization, and polygon-merging, (2) Building envelope extraction using PIX4d modeling software to reconstruct a building 3D model, (3) Visualization tool for viewing images from the 3D model, (4) Window-to-wall ratio (WWR) using state-of-art deep neural network semantic segmentation, (5) Envelope thermal anomaly detection using an unsupervised machine learning clustering algorithm, and (6) Rooftop energy equipment detection based on an object detection algorithm. The testing of this method involved a comparison of additional ML-generated information overlaid on current ‘state-of-practice’ audit and remote assessment baselines using evaluation metrics: labor time and associated cost, marginal benefits of using ML-generated information in workflows for audits and remote assessments, integration potential with existing processes and tools, and replicability/scalability of the method. In two test buildings in California that had comprehensive drawings and meter data available, the ML method effectively generated a building footprint, envelope, rooftop equipment, WWR, and locations of envelope thermal anomalies. Projected target segments of the ML method are sites with minimal drawings and energy data, and underserved sectors such as multistoried housing, disadvantaged communities, and schools for which the ML method can enable identification of building asset characteristics and prioritization of envelope retrofits and decentralized energy equipment retrofits.

Singh, Reshma↗

Project Methods to Enhance Wellbore Cement Integrity with Microbially-Induced Calcite Precipitation (MICP) (Final Scientific/Technical Report)

The goal of this project was to develop improved methods for sealing compromised wellbore cement in leaking oil and gas wells, thereby reducing the risk of unwanted upward fluid migration. Novel methods for improving wellbore integrity, such as microbially induced calcite precipitation (MICP), can reduce leakage potential, improve the safety of fossil fuel extraction, improve the public perception of hydraulic fracturing, and promote environmentally-prudent unconventional oil and gas development. Microbes, with the urease enzyme, can catalyze the chemical reaction of urea hydrolysis to induce the precipitation of calcium carbonate which can be used as a cementitious material to seal leakage pathways. In this project, methods to promote robust bio-composite cementitious materials were designed and tested in the laboratory. Scale-up of those methods were tested in meso-scale reactor systems and in field applications. In this report, in Section One, we describe laboratory efforts to develop injection strategies to promote precipitation in wellbore analogs and determine the strength of the bio-composite cements as compared to fine cement. In Section Two, we describe the efforts to scale up the work and study the use of materials that can be used in field application, for example exploring the use of calcium chloride ice melt or urea fertilizer as source chemicals. In Section Three, the three field trials (methods and results) performed as part of the project are described and summarized. At the end of the report is a comprehensive summary and conclusion section which highlights the key findings of the project. The work performed during this project significantly advanced the technology readiness level (TRL) of the MICP wellbore sealing strategy.

02 PETROLEUM↗

Methods to Enhance Wellbore Cement Integrity with Microbially-Induced Calcite Precipitation (MICP) (Final Scientific/ Technical Report)

The goal of this project was to develop improved methods for sealing compromised wellbore cement in leaking oil and gas wells, thereby reducing the risk of unwanted upward fluid migration. Novel methods for improving wellbore integrity, such as microbially induced calcite precipitation (MICP), can reduce leakage potential, improve the safety of fossil fuel extraction, improve the public perception of hydraulic fracturing, and promote environmentally-prudent unconventional oil and gas development. Microbes, with the urease enzyme, can catalyze the chemical reaction of urea hydrolysis to induce the precipitation of calcium carbonate which can be used as a cementitious material to seal leakage pathways. In this project, methods to promote robust bio-composite cementitious materials were designed and tested in the laboratory. Scale-up of those methods were tested in meso-scale reactor systems and in field applications. In this report, in Section One, we describe laboratory efforts to develop injection strategies to promote precipitation in wellbore analogs and determine the strength of the bio-composite cements as compared to fine cement. In Section Two, we describe the efforts to scale up the work and study the use of materials that can be used in field application, for example exploring the use of calcium chloride ice melt or urea fertilizer as source chemicals. In Section Three, the three field trials (methods and results) performed as part of the project are described and summarized. At the end of the report is a comprehensive summary and conclusion section which highlights the key findings of the project. The work performed during this project significantly advanced the technology readiness level (TRL) of the MICP wellbore sealing strategy.

03 NATURAL GAS↗

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗

Study of self-interaction-errors in barrier heights using locally scaled and Perdew-Zunger self-interaction methods

We study the effect of self-interaction errors on the barrier heights of chemical reactions. For this purpose we use the well-known Perdew-Zunger [J. P. Perdew and A. Zunger, Phys. Rev. B,23, 5048 (1981)] self-interaction-correction(PZSIC), as well two variations of the recently developed, locally scaled self-interaction correction (LSIC) [Zope et al.,J. Chem. Phys.151, 214108 (2019)] to study the barrier heights of the BH76 benchmark dataset. Our results show that both PZSIC and especially the LSIC methods improve the barrier heights relative to the local density approximation(LDA). The version of LSIC that uses the iso-orbital indicator z as a scaling factor gives a more consistent improvement than an alternative version that uses an orbital-dependent factor based on the ratio of orbital densities to the total electron density. We show that LDA energies evaluated using the self-consistent and self-interaction-free PZSIC densities can be used to assess density-driven errors. Furthermore, the LDA reaction barrier errors for the BH76 set are found to contain significant density-driven errors for all types of reactions contained in the set, but the corrections due to adding SIC to the functional are much larger than those stemming from the density for the hydrogen transfer reactions and of roughy equal size for the non-hydrogen transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measuring Gas Transport and Sorption in Large Intact Geologic Specimens via the Piezometric Method

Uptake of noble gases into heterogeneous geologic core samples was measured using a piezometric methodology. In addition to measuring accessible porosity—as with gas pycnometry—by monitoring the rate of pressure decay, this method can also be used to estimate the gas effective diffusivity in the sample. In contrast to previous applications of this method, where milligram quantities of fractured grains are characterized, here approximately kilogram core samples were left intact when tested. In doing so, a more representative sample of the heterogeneous field geology is provided. Additionally, alteration of the pore structure and connectivity during sample preparation is avoided. Furthermore, to scale the piezometric method from milligrams to kilograms, the system was designed to operate at medium vacuum (1 to 100 Pa) to restrict transport in pores less than approximately 60 µm to large Knudsen numbers. To test the system performance, two samples of interest were selected: a rhyolitic welded tuff from Blue Canyon Dome at the Energetic Materials Research and Testing Center and a zeolitized non-welded rhyolitic tuff from the Nevada National Security Site. Three noble gases were utilized in this test series; Argon and xenon as they are of direct interest to nuclear monitoring efforts and helium as it is a weakly adsorbing reference standard. Additionally, mercury intrusion porosimetry measurements were made on subsamples of the core to compare the observed porosity by the two methods and to discuss gas transport rates in the context of the measured pore distribution.

58 GEOSCIENCES↗

Area-Preserving Mapping of 3D Carotid Ultrasound Images Using Density-Equalizing Reference Map

Carotid atherosclerosis is a focal disease at the bifurcations of the carotid artery. To quantitatively monitor the local changes in the vessel-wall-plus-plaque thickness (VWT) and compare the VWT distributions for different patients or for the same patients at different ultrasound scanning sessions, a mapping technique is required to adjust for the geometric variability of different carotid artery models. In this work, we propose a novel method called density-equalizing reference map (DERM) for mapping 3D carotid surfaces to a standardized 2D carotid template, with an emphasis on preserving the local geometry of the carotid surface by minimizing the local area distortion. The initial map was generated by a previously described arc-length scaling (ALS) mapping method, which projects a 3D carotid surface onto a 2D non-convex L-shaped domain. A smooth and area-preserving flattened map was subsequently constructed by deforming the ALS map using the proposed algorithm that combines the density-equalizing map and the reference map techniques. This combination allows, for the first time, one-to-one mapping from a 3D surface to a standardized non-convex planar domain in an area-preserving manner. Evaluations using 20 carotid surface models show that the proposed method reduced the area distortion of the flattening maps by over 80% as compared to the ALS mapping method.

60 APPLIED LIFE SCIENCES↗

In-Process Monitoring and Structural Health Monitoring of Large-Scale Additive Manufacturing Using Acoustic Emission Technique

ORNL collaborated with MISTRAS Group, Inc. to investigate acoustic emission (AE) as a structural health monitoring (SHM) method for large-scale additive manufacturing (AM). Large-scale AM is being adapted as method of producing large structures in a short lead time and cost-effective way. With the growing advancement in AM techniques and application, machine monitoring and part qualification is highly needed. There has been leading research focused on the manufacturing, feedstock material but minimum research on the SHM, defect detection, and nondestructive evaluation (NDE) for AM. Scanning large structure using conventional nondestructive testing (NDT) techniques, such as ultrasound or X-ray, and searching for potential defects can be very time consuming, challenging and cost prohibitive. AE is a passive technique that can be used to monitor and locate defect progression in large structure by distributing group of sensors around the part. This project utilized AE technique and system manufactured/designed by MISTRAS Group to monitor large-scale AM equipment (i.e. Big Area Additive Manufacturing (BAAM) system located at the Oak Ridge National Laboratory – Manufacturing Demonstration Facility (ORNL-MDF) and the printed parts it produces. The AE system provided valuable insight on defect development/progression during and post-printing process.

36 MATERIALS SCIENCE↗