Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Extrapolation”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

Highly Accurate Quartic Force Field and Rovibrational Spectroscopic Constants for the Azirinyl Cation (c-C2NH2+) and Its Isomers

Azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high resolution experimental or astronomical observations.

Partha P Bera↗
DOE OSTI2026

First results of the NEXT-100 detector using $^{83m}$Kr decays

The NEXT collaboration is investigating the double beta decay of $ ^{136}\mathrm Xe$ using high-pressure gas electroluminescent time projection chambers, which provide excellent energy resolution together with a robust topological signature. Operating at the Laboratorio Subterráneo de Canfranc (LSC) and building on the success of the NEXT-White detector, the NEXT-100 apparatus began commissioning in May 2024 and started operation with xenon at a pressure of 4 bar in October 2024. We report here the first results obtained with NEXT-100 using low-energy calibration data from $^{83\text {m}}\textrm{Kr}$ decays, which allow mapping of the detector response in the active volume and monitoring of its stability over time. After homogenizing the light response, we achieve an energy resolution of 4.37% FWHM at 41.5 keV for $^{83\text {m}}\textrm{Kr}$ point-like energy deposits contained in a radius of 425 mm. In a fiducial region representing the operating conditions of NEXT-100 at 10 bar we obtain an improved energy resolution of 4.16% FWHM. These results are in good agreement with that obtained in NEXT-White, and an $E^{-1/2}$ extrapolation to $Q_{\beta \beta }$ yields an energy resolution close to 0.5% FWHM, well below the 1% FWHM design target.

Martínez-Lema, G. [Donostia Intl. Phys. Ctr., San ↗
DOE OSTI2020

Global ocean dimethyl sulfide climatology estimated from observations and an artificial neural network

Marine dimethyl sulfide (DMS) is important to climate due to the ability of DMS to alter Earth's radiation budget. Knowledge of the global-scale distribution, seasonal variability, and sea-to-air flux of DMS is needed in order to improve understanding of atmospheric sulfur, aerosol/cloud dynamics, and albedo. Here we examine the use of an artificial neural network (ANN) to extrapolate available DMS measurements to the global ocean and produce a global climatology with monthly temporal resolution. A global database of 82 996 ship-based DMS measurements in surface waters was used along with a suite of environmental parameters consisting of latitude–longitude coordinates, time of day, time of year, solar radiation, mixed layer depth, sea surface temperature, salinity, nitrate, phosphate, and silicate. Linear regressions of DMS against the environmental parameters show that on a global-scale mixed layer depth and solar radiation are the strongest predictors of DMS. These parameters capture ~9 % and ~7 % of the raw DMS data variance, respectively. Multilinear regression can capture more of the raw data variance (~39 %) but strongly underestimates DMS in high-concentration regions. In contrast, the artificial neural network captures ~66 % of the raw data variance in our database. Like prior climatologies our results show a strong seasonal cycle in surface ocean DMS with the highest concentrations and sea-to-air fluxes in the high-latitude summertime oceans. We estimate a lower global sea-to-air DMS flux (20.12±0.43 Tg S yr –1 ) than the prior estimate based on a map interpolation method when the same gas transfer velocity parameterization is used. Our sensitivity test results show that DMS concentration does not change unidirectionally with each of the environmental parameters, which emphasizes the interactions among these parameters. The ANN model suggests that the flux of DMS from the ocean to the atmosphere will increase with global warming. Given that larger DMS fluxes induce greater cloud albedo, this corresponds to a negative climate feedback.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1942

The direct-stress fatigue strength of 17S-T aluminum alloy throughout the range from 1/2 to 500,000,000 cycles of stress

Fatigue-test were conducted on six specimens made from 3/4-inch-diameter 17S-T rolled-and-drawn rod for the purpose of obtaining additional data on the fatigue life of the material at stresses up to the static strength. The specimens were tested in direct tension using a stress range from zero to a maximum in tension. A static testing machine was used to apply repeated loads in the case of the first three specimens; the other three specimens were tested in a direct tension-compression fatigue machine. The direct-stress fatigue curve obtained for the material indicates that, in the range of stresses above about two-thirds the tensile strength, the fatigue strength is higher than might be expected by simply extrapolating the ordinary curve of stress plotted against the number of cycles determined at lower stresses.

Hartmann, E C↗
NASA NTRS2018

Data Analysis Techniques for Fan Performance in Highly-Distorted Flows from Boundary Layer Ingesting Inlets

The design of a unique distortion-tolerant fan for a high-bypass ratio boundary-layer ingesting propulsion system has been completed and a rig constructed and tested in the NASA Glenn 8’x6’ wind tunnel. Processing the data from the experiment presented some interesting challenges because of the complexity of the experimental setup and the flow through the test rig. The experiment was run in three phases, each of which employed a unique complement of inlet throat and fan face instrumentation to avoid the blockage that would have resulted from simultaneously installing all of the rakes. The measurement from the individual test points were subsequently combined to compute the overall stage performance. A CFD model of the experiment was used to gain understanding of the flow field and to test some of the techniques proposed for interpolating and extrapolating the measurements into regions where measurements were not made. This capability became extremely useful when it was discovered that there was an unexpected total temperature distortion in the tunnel. The CFD model was modified by inserting a total temperature profile at the upstream boundary that mimicked the measured distortion where measurements were available and that CFD solution was used to investigate methods to infer the complete total temperature field at the fan face.

Hardin, Larry W.↗
DOE OSTI2024

Characterizing model uncertainties in simulated coast-to-offshore wind over the northeast U.S. using multi-platform measurements from the TCAP field campaign

Numerical weather prediction (NWP) models, such as the Weather Research and Forecasting (WRF) model, are widely used to provide estimates of the offshore wind energy resource owing to their large spatial coverage compared to available observations. Nevertheless, spatiotemporal distribution of model biases is highly dependent on factors including model configuration, location, and the interplay of multi-scale physical processes. Here, in this study, we focus on the characterization of model uncertainties in simulated coast-to-offshore winds over the northeast U.S., by varying sea surface temperature (SST) forcings, surface layer (SL) and planetary boundary layer (PBL) parameterizations, as well as identifying biases that may be directly passed from initial and boundary conditions. Multiple measurements, including aircraft data collected during the U.S. Department of Energy's Two-Column Aerosol Project (TCAP) experiment, are used to constrain the model results and facilitate quantitative comparisons. Our analysis indicates while SST forcing has notable impacts on simulated air temperature and moisture within PBL, the modeled winds are in general more sensitive to the choices of SL and PBL physics than to SST. The model’s forcing data not only controls the vertical dependence of wind speed errors, but also alters regional variability in wind speed’s spatial correlation. Bias comparisons between ERA5 reanalysis and ensemble simulations revealed significant similarity, particularly in wind speed biases during winter, underscoring their dependency on initial and boundary conditions. Coastal and offshore near-surface wind speed biases tend to exhibit much higher similarity in winter than in summer due to the presence of much stronger and more persistent synoptic wind conditions. This study highlights the importance of accurate atmospheric forcing and parameterization choices in improving wind forecasts and suggests the potential for extrapolating coastal wind biases to offshore locations, aiding wind energy forecasting and informing the Wind Forecast Improvement Project-3 (WFIP3).

17 WIND ENERGY↗
DOE OSTI2022

Accuracy of power spectra in dissipationless cosmological simulations

Abstract We exploit a suite of large N-body simulations (up to N = 40963) performed with abacus, of scale-free models with a range of spectral indices n, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At 1 per cent level it starts at early times from a fraction of $k_\Lambda$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, ${\sim}2{-}3 k_\Lambda$ at the very latest times we evolve to. At the $5{{\ \rm per\ cent}}$ level, accuracy extends up to wavenumbers of order $5k_\Lambda$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple n-dependence, allowing a extrapolation to place conservative bounds on the accuracy of N-body simulations of non-scale-free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small-scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_\Lambda$ may need revision in light of our convergence analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2005

Fitting Nonlinear Curves by use of Optimization Techniques

MULTIVAR is a FORTRAN 77 computer program that fits one of the members of a set of six multivariable mathematical models (five of which are nonlinear) to a multivariable set of data. The inputs to MULTIVAR include the data for the independent and dependent variables plus the user s choice of one of the models, one of the three optimization engines, and convergence criteria. By use of the chosen optimization engine, MULTIVAR finds values for the parameters of the chosen model so as to minimize the sum of squares of the residuals. One of the optimization engines implements a routine, developed in 1982, that utilizes the Broydon-Fletcher-Goldfarb-Shanno (BFGS) variable-metric method for unconstrained minimization in conjunction with a one-dimensional search technique that finds the minimum of an unconstrained function by polynomial interpolation and extrapolation without first finding bounds on the solution. The second optimization engine is a faster and more robust commercially available code, denoted Design Optimization Tool, that also uses the BFGS method. The third optimization engine is a robust and relatively fast routine that implements the Levenberg-Marquardt algorithm.

Hill, Scott A.↗
NASA NTRSDate not supplied

SWIPE: Spectral Water Inversion Processor and Emulator

Degradation of Earth’s inland water resources due to anthropogenic perturbations and climate anomalies at both local and global scales continues to place human health at substantial risk. There is now a growing necessity to develop pragmatic approaches that allow timely and effective extrapolation of local processes, to spatially resolved global products, and to promote operational and sustainable resource policy management. This presentation will be discussing the progress made developing SWIPE: Spectral Water Inversion Processor and Emulator. SWIPE is a platform for advanced modeling of coastal and inland aquatic habitats. The goal is create a comprehensive and cohesive system to leverage recent advancements in computation and machine learning to develop a synthetic training ground for sensitivity studies and algorithm development. The four principal facets of SWIPE include: 1. Advanced two-layer coated sphere bio-optical modeling and GPU radiative transfer modeling, 2. Big Data involving massive synthetic spectral libraries of optical properties of various global aquatic particles, surface reflectance, and top-of-atmosphere reflectance, all at hyperspectral resolution leveraging high-end computing systems at NASA Ames Research Center, 3. Deep Learning for algorithm development for water quality inversion of concentrations of common biogeophysical variables as well as optics, full uncertainty characterization by water type, and forward emulation, and lastly, 4. Image Processing for application of developed retrieval algorithms for both hyperspectral and multispectral sensors with experimental corrections for global adjacency, noise, sunglint, and benthic reflectance. This presentation will demonstrate the Equivalent Algal Populations (EAP) two-layer coated sphere scattering model which has been used develop spectral libraries of hyperspectral inherent optical properties of roughly 80 species of phytoplankton, covering 15 different classes and nine taxonomic functional types. The EAP model was also used to derive spectral properties of 10 different non-algal particle functional types. Examples of how the SMART-G (Speed-up Monte-carlo Advanced Radiative Transfer using GPU) radiative transfer code is used to model optically complex aquatic signals will be presented and discussed in the context of creating a massive synthetic database which can leverage the full power of next generation machine learning techniques and high end computing for water quality inversion. We will discuss our active investigation in things like appropriate model architectures, dimensionality reduction techniques such as PCA and autoencoders, uncertainty quantification and abstaining, and which variables actually benefit most from hyperspectral information versus multispectral resolution. We are also curious about questions relating to cost/benefit analysis in terms of computation resources, neural network complexity, and data volumes. Answers to these questions will hopefully elaborate on cost efficiency for potential future sensor design considerations.

SWIPE↗
DOE OSTI2022

Vickers hardness prediction from machine learning methods

Abstract The search for new superhard materials is of great interest for extreme industrial applications. However, the theoretical prediction of hardness is still a challenge for the scientific community, given the difficulty of modeling plastic behavior of solids. Different hardness models have been proposed over the years. Still, they are either too complicated to use, inaccurate when extrapolating to a wide variety of solids or require coding knowledge. In this investigation, we built a successful machine learning model that implements Gradient Boosting Regressor (GBR) to predict hardness and uses the mechanical properties of a solid (bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio) as input variables. The model was trained with an experimental Vickers hardness database of 143 materials, assuring various kinds of compounds. The input properties were calculated from the theoretical elastic tensor. The Materials Project’s database was explored to search for new superhard materials, and our results are in good agreement with the experimental data available. Other alternative models to compute hardness from mechanical properties are also discussed in this work. Our results are available in a free-access easy to use online application to be further used in future studies of new materials at www.hardnesscalculator.com .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Critical examination of constraints on the equation of state of dense matter obtained from GW170817

The correlation of the tidal polarizabilities Λ 1 - Λ 2 for GW170817 is predicted by combining dense-matter equations of state (EOSs) that satisfy nuclear physics constraints with the chirp mass and mass asymmetry for this event. Our models are constrained by calculations of the neutron-matter EOS using chiral effective field theory Hamiltonians with reliable error estimates up to once or twice the nuclear saturation density. In the latter case, we find that GW170817 does not improve our understanding of the EOS. We contrast two distinct extrapolations to higher density: a minimal model (MM) which assumes that the EOS is a smooth function of density described by a Taylor expansion and a more general model parametrized by the speed of sound that admits phase transitions. This allows us to identify regions in the Λ 1 - Λ 2 plots that could favor the existence of new phases of matter in neutron stars. We predict the combined tidal polarizability of the two neutron stars in GW170817 to be 80 ≤ $\tilde{Λ}$ ≤ 580 ( 280 ≤$\tilde{Λ}$ ≤ 480 for the MM), which is smaller than the range suggested by the LIGO-Virgo data analysis. Our analysis also shows that GW170817 requires a NS with M = 1.4 M ⊙ to have a radius 9.0 km < R 1.4 < 13.6 km ( 11.3 km < R 1.4 < 13.6 km for the MM).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

PRIMA-X - Performance Retargeting of Instrumentation, Measurement, and Analysis Technologies for Exascale Computing

The PRIMA-X (Performance Retargeting of Instrumentation, Measurement, and Analysis Technologies for Exascale Computing) project is the successor of the DOE PRIMA (Performance Refactoring of Instrumentation, Measurement, and Analysis Technologies for Petascale Computing) project, which addressed the challenge of creating a core measurement infrastructure that would serve as a common platform for both integrating leading parallel performance systems (notably TAU and Scalasca) and developing next-generation scalable performance tools. The PRIMA-X project shifts the focus away from refactorization of robust performance tools towards a re-targeting of the parallel performance measurement and analysis architecture for extreme scales. The massive concurrency, asynchronous execution dynamics, hardware heterogeneity, and multi-objective prerequisites (performance, power, resilience) that identify exascale systems introduce fundamental constraints on the ability to carry forward existing performance methodologies. In particular, there must be a deemphasis of per-thread observation techniques to significantly reduce the otherwise unsustainable flood of redundant performance data. Instead, it will be necessary to assimilate multi-level resource observations into macroscopic performance views, from which resilient performance metrics can be attributed to the computational features of the application. This requires a scalable framework for node-level and system-wide monitoring and runtime analyses of dynamic performance information. Also, the interest in optimizing parallelism parameters with respect to performance and energy drives the integration of tool capabilities in the exascale environment further. Initially, PRIMA-X was a collaborative project between the University of Oregon (lead institution) and the German Research School for Simulation Sciences (GRS). Because Prof. Wolf, the PI at GRS, accepted a position as full professor at Technische Universitaet Darmstadt (TU Darmstadt) starting February 1st, 2015, the work of GRS was continued at TU Darmstadt. This report covers the work at TU Darmstadt after the transition. The first main accomplishment of TU Darmstadt is the development of different techniques to aggregate performance data on the level of threads. TU Darmstadt evaluated different schemes designed at GRS during the first phase of this project and integrated them into Score-P, a widely used performance measurement framework for HPC applications. The second main accomplishment is a substantial increase of Score-P’s scalability, achieved by improving its internal representation of the underlying system. Third, we extended our performance-modeling tool Extra-P to automatically create models with more than one input parameter (e.g., process count and input size), which was previously impossible. These models allow the performance behavior of an application to be extrapolated to larger machines. Specific applications of the new multi-parameter capability realized in this project include modeling the isoefficiency function of task-based applications and a novel co-design methodology. Further extensions of Extra-P comprise online modeling functionality, automatic refinement of the model search space, and segmented performance modeling. Further contributions include new locality metrics for both shared- and distributed-memory applications, a new way of measuring the impact of inter-application inference on the execution time of an application, and finally an automatic method that selects appropriate OpenMP constructs for the parallelization of sequential code and classifies the affected variables based on their sharing semantics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Measurement of the 0.3-8.5 MeV Galactic gamma-ray spectrum from the Galactic center direction

The low-energy gamma-ray spectrum from the direction of the Galactic center is determined using data obtained with the SMM Gamma-Ray Spectrometer. It is found that the diffuse gamma-ray spectrum from the Galactic center region can be interpreted in a straightforward way as the sum of five components of a presented equation. The components include a hard power law dominating the continuum at high energies caused principally by cosmic ray electron bremsstrahlung radiation, two narrow lines due to Al-26 decay and positron annihilation, an excess continuum component below 0.511 MeV consistent with the annihilation of positrons by formation of Ps, and a soft power law at low energies which is consistent with an extrapolation upward in energy of known hard X-ray sources in the Galactic center region.

Harris, M. J.↗
NASA NTRS2016

TEM Analyses of Itokawa Regolith Grains and Lunar Soil Grains to Directly Determine Space Weathering Rates on Airless Bodies

Samples returned from the moon and Asteroid Itokawa by NASA's Apollo Missions and JAXA's Hayabusa Mission, respectively, provide a unique record of their interaction with the space environment. Space weathering effects result from micrometeorite impact activity and interactions with the solar wind. While the effects of solar wind interactions, ion implantation and solar flare particle track accumulation, have been studied extensively, the rate at which these effects accumulate in samples on airless bodies has not been conclusively determined. Results of numerical modeling and experimental simulations do not converge with observations from natural samples. We measured track densities and rim thicknesses of three olivine grains from Itokawa and multiple olivine and anorthite grains from lunar soils of varying exposure ages. Samples were prepared for analysis using a Leica EM UC6 ultramicrotome and an FEI Quanta 3D dual beam focused ion beam scanning electron microscope (FIB-SEM). Transmission electron microscope (TEM) analyses were performed on the JEOL 2500SE 200kV field emission STEM. The solar wind damaged rims on lunar anorthite grains are amorphous, lack inclusions, and are compositionally similar to the host grain. The rim width increases as a smooth function of exposure age until it levels off at approximately 180 nm after approximately 20 My (Fig. 1). While solar wind ion damage can only accumulate while the grain is in a direct line of sight to the Sun, solar flare particles can penetrate to mm-depths. To assess whether the track density accurately predicts surface exposure, we measured the rim width and track density in olivine and anorthite from the surface of rock 64455, which was never buried and has a surface exposure age of 2 My based on isotopic measurements. The rim width from 64455 (60-70nm) plots within error of the well-defined trend for solar wind amorphized rims in Fig. 1. Measured solar flare track densities are accurately reflecting the surface exposure of the grains. Track densities correlate with the amorphous rim thicknesses. While the space-weathered rims of anorthite grains are amorphous, the space-weathered rims on both Itokawa and lunar olivine grains show solar wind damaged rims that are not amorphous. Instead, the rims are nanocrystalline with high dislocation densities and sparse inclusions of nanophase Fe metal. The rim thicknesses on the olivine grains also correlate with track density. The Itokawa olivine grains have track densities that indicate surface exposures of approximately 10(exp 5) years. Longer exposures (up to approximately 10(exp 7) years) do not amorphize the rims, as evidenced by lunar soil olivines with high track densities (approximately 10(exp 11) cm(exp -2)). From the combined data, shown in Fig. 1, it is clear that olivine is damaged (but not amorphized) more rapidly by the solar wind compared to anorthite. The olivine damaged rim forms quickly (in approximately 10(exp 6) y) and saturates at approximately 120nm with longer exposure time. The anorthite damaged rims form more slowly, amorphize, and grow thicker than the olivine rims. This is in agreement with numerical modeling data which predicts that solar wind damaged rims on anorthite will be thicker than olivine. However, the models predict that both olivine and anorthite rims will amorphize and reach equilibrium widths in less than 10(exp 3) y, in contrast to what is observed for natural samples. Laboratory irradiation experiments, which show rapid formation of fully amorphous and blistered surfaces from simulated solar wind exposures are also in contrast to observations of natural samples. These results suggest that there is a flux dependence on the type and extent of irradiation damage that develops in olivine. This flux dependence suggests that great caution be used in extrapolating between high-flux laboratory experiments and the natural case, as demonstrated by. We constrain the space weathering rate through analysis of returned samples. Provided that the track densities and the solar wind damaged rim widths exhibited by the Itokawa grains are typical of the fine-grained regions of Itokawa, then the space weathering rate is on the order of 10(exp 5) y. Space weathering effects in lunar soils saturate within a few My of exposure while those in Itokawa regolith grains formed in approximately 10(exp 5) y. Olivine and anorthite respond differently to solar wind irradiation. The space weathering effects in olivine are particularly difficult to reconcile with laboratory irradiation studies and numerical models. Additional measurements, experiments, and modeling are required to resolve the discrepancies among the observations and calculations involving solar wind amorphization of different minerals on airless bodies.

Berger, Eve L.↗
NASA NTRS2006

Estimating the NO(x) produced by lightning from GOME and NLDN data: A case study in the Gulf of Mexico

Nitrogen oxides (NO(x)=NO+NO2) play an important role in tropospheric chemistry, in particular in catalytic ozone production. Lightning provides a natural source of nitrogen oxides, dominating the production in the tropical upper troposphere, with strong impact on tropospheric ozone and the atmosphere s oxidizing capacity. Recent estimates of lightning produced NO(x) (LNO(x)) are of the order of 5 Tg [N] per year with still high uncertainties in the range of one order of magnitude. The Global Ozone Monitoring Experiment (GOME) on board the ESA-satellite ERS-2 allows the retrieval of tropospheric column densities of NO2 on a global scale. Here we present the GOME NO2 measurement directly over a large convective system over the Gulf of Mexico. Simultaneously, cloud-to-ground (CG) flashes are counted by the U.S. National Lightning Detection Network (NLDN(TradeMark)), and extrapolated to include intracloud (IC)+CG flashes based on a climatological IC:CG ratio derived from NASA s space-based lightning sensors. A series of 14 GOME pixels shows largely enhanced column densities over thick and high clouds, coinciding with strong lightning activity. The enhancements can not be explained by transport of anthropogenic NO(x) and must be due to fresh production of LNO(x). A quantitative analysis, accounting in particular for the visibility of LNO, from satellite, yields a LNO, production of 90 (32-240) moles of NO(x), or 1.3 (0.4-3.4) kg [N], per flash. If simply extrapolated, this corresponds to a global LNO(x) production of 1.7 (0.6-4.7) Tg [N]/yr.

Beirle, S.↗
DOE OSTI2021

Robust verification of stochastic simulation codes

We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. Furthermore, this approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations..

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Measurements of W + W − production cross-sections in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

Measurements of W + W − → e ± νμ ∓ ν production cross-sections are presented, providing a test of the predictions of perturbative quantum chromodynamics and the electroweak theory. The measurements are based on data from pp collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS detector at the Large Hadron Collider in 2015–2018, corresponding to an integrated luminosity of 140 fb −1 . The number of events due to top-quark pair production, the largest background, is reduced by rejecting events containing jets with b-hadron decays. An improved methodology for estimating the remaining top-quark background enables a precise measurement of W + W − cross-sections with no additional requirements on jets. The fiducial W + W − cross-section is determined in a maximum-likelihood fit with an uncertainty of 3.1%. The measurement is extrapolated to the full phase space, resulting in a total W + W − cross-section of 127 ± 4 pb. Differential cross-sections are measured as a function of twelve observables that comprehensively describe the kinematics of W + W − events. The measurements are compared with state-of-the-art theory calculations and excellent agreement with predictions is observed. A charge asymmetry in the lepton rapidity is observed as a function of the dilepton invariant mass, in agreement with the Standard Model expectation. A CP-odd observable is measured to be consistent with no CP violation. Limits on Standard Model effective field theory Wilson coefficients in the Warsaw basis are obtained from the differential cross-sections.

Accelerator Physics↗