Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

LES and RANS Spray Combustion Analysis of OME 3-5 and n-Dodecane

Clean-burning oxygenated and synthetic fuels derived from renewable power, so-called e-fuels, are a promising pathway to decarbonize compression–ignition engines. Polyoxymethylene dimethyl ethers (PODEs or OMEs) are one candidate of such fuels with good prospects. Their lack of carbon-to-carbon bonds and high concentration of chemically bound oxygen effectively negate the emergence of polycyclic aromatic hydrocarbons (PAHs) and even their precursors like acetylene (C 2 H 2 ), enabling soot-free combustion without the soot-NO x trade-off common for diesel engines. The differences in the spray combustion process for OMEs and diesel-like reference fuels like n-dodecane and their potential implications on engine applications include discrepancies in the observed ignition delay, the stabilized flame lift-off location, and significant deviations in high-temperature flame morphology. For CFD simulations, the accurate modeling and prediction of these differences between OMEs and n-dodecane proved challenging. This study investigates the spray combustion process of an OME 3 − 5 mixture and n-dodecane with advanced optical diagnostics, Reynolds-Averaged Navier–Stokes (RANS), and Large-Eddy Simulations (LESs) within a constant-volume vessel. Cool-flame and high-temperature combustion were measured simultaneously via high-speed (50 kHz) imaging with formaldehyde (CH 2 O) planar laser-induced fluorescence (PLIF) representing the former and line-of-sight OH* chemiluminescence the latter. Both RANS and LES simulations accurately describe the cool-flame development process with the formation of CH 2 O. However, CH 2 O consumption and the onset of high-temperature reactions, signaled by the rise of OH* levels, show significant deviations between RANS, LES, and experiments as well as between n-dodecane and OME. A focus is set on the quality of the simulated results compared to the experimentally observed spatial distribution of OH*, especially in OME fuel-rich regions. The influence of the turbulence modeling is investigated for the two distinct ambient temperatures of 900 K and 1200 K within the Engine Combustion Network Spray A setup. The capabilities and limitations of the RANS simulations are demonstrated with the initial cool-flame propagation and periodic oscillations of CH 2 O formation/consumption during the quasi-steady combustion period captured by the LES.

CFD↗
DOE OSTI2026

Vector boson scattering and anomalous quartic couplings in final states with ℓνqq or ℓℓqq plus jets using proton-proton collisions at $\sqrt{s}=13$ TeV

A measurement is presented of the electroweak vector boson scattering production of ZV (V = W, Z) boson pairs associated with two jets in proton-proton collisions at a center-of-mass energy of 13 TeV. The data, corresponding to an integrated luminosity of 138 fb −1 , were collected at the CERN LHC with the CMS detector during the 2016–2018 data-taking period. The analysis targets final states with a pair of isolated electrons or muons from Z boson decays and three or four jets, depending on the momentum of the vector boson that decays into quarks. Signal strength is measured for events characterized by a large invariant mass of two forward jets with a wide pseudorapidity gap between them. The electroweak production of ZV in association with two jets is measured with an observed (expected) significance of 1.3 (1.8) standard deviations. A combination of the analyses of ZV channel and the previously published WV channel in the lepton plus jets final state places constraints on effective field theory parameters that describe anomalous electroweak production of WW, WZ, and ZZ boson pairs in association with two jets. Several world best limits are set on anomalous quartic gauge couplings in terms of dimension-8 standard model effective field theory operators.

Hadron-Hadron Scattering↗
DOE OSTI2020

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗
DOE OSTI2020

Baryon $\sigma$ terms in SU(3) BChPT x 1/N c

ChPT and the 1/N c expansion provide systematic frameworks for the strong interactions at low energy. A combined framework of both expansions has been developed and applied for baryons with three light-quark-flavors. The small scale expansion of the combined approach is identified as the ξ-expansion, in which the power counting of the expansions is linked according to O(p) = O(1/N c ) = O(ξ). The physical baryon masses as well as lattice QCD baryon masses for different quark mass masses are analyzed to O(ξ 3 ) in that framework. σ terms are addressed using the Feynman Hellmann theorem. For the nucleon, a useful connection between the deviation of the Gell-Mann-Okubo relation and the σ term σ 8N associated with the scalar density $\bar{u}u+\bar{d}d-2\bar{s}s$ is identified. In particular, the deviation from the tree level relation σ 8N = $\frac{1}{3}$(2m N -m Σ -m Ξ ), which gives rise to the so called σ-term puzzle, is studied in the ξ-expansion. A large correction non-analytic in ξ results for that relation, making plausible the resolution of the puzzle. Issues with the determination of the strangeness σ terms are discussed, emphasizing the need for lattice calculations at smaller m s for better understanding the range of validity of the effective theory. The analysis presented here leads to σ πN = 69(10) MeV and σ πΔ = 60(10) MeV.

Fernando, Ishara P.↗
DOE OSTI2026

On the Properties of Cosmological Ionization Fronts

We investigate the properties of cosmological ionization fronts during the epoch of reionization using the Cosmic Reionization on Computers simulations. By analyzing reionization timing maps, we characterize ionization front velocities and curvatures and their dependence on the density structure of the intergalactic medium (IGM). The velocity distribution of ionization fronts in the simulations indicates that while the barrier-crossing analytical model captures the overall shape in high-velocity regions, it fails to reproduce the low-velocity tail, highlighting the non-Gaussian nature of the IGM’s density field. Ionization front velocities are inversely correlated with local density, propagating faster in underdense regions and more slowly in overdense environments. Faster ionization fronts also lead to higher postionization temperatures, reaching a plateau at ∼2 × 10$^{4}$ K for velocities exceeding 3000 km s$^{−1}$. Examining curvature statistics further establishes a connection between ionization front structure and the normalized density contrast ν, with trends in overdense regions aligning well with barrier-crossing model predictions, while deviations appear in underdense environments due to model limitations. These results provide a detailed characterization of ionization front dynamics and their interaction with the underlying density field, bridging small-scale reionization physics with large-scale observables such as the 21 cm signal and the IGM’s thermal history.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Equations of motion for maneuvering flexible spacecraft

This paper is concerned with the derivation of the equations of motion for maneuvering flexible spacecraft both in orbit and in an earth-based laboratory. The structure is assumed to undergo large rigid-body maneuvers and small elastic deformations. A perturbation approach is presented in which the quantities defining the rigid-body maneuver are regarded as the unperturbed motion and the elastic motions and deviations from the rigid-body motions are regarded as the perturbed motion. The perturbation equations are linear, non-self-adjoint, and with time-dependent coefficients. A maneuver force distribution exciting the least amount of elastic deformation of the spacecraft is developed. Numerical results highlight the vibration caused by rotational maneuvers.

Meirovitch, L.↗
NASA NTRS1973

A fine pointing system for the large space telescope

The large space telescope (LST) developed by NASA requires ultrahigh pointing stability within 0.0005 arc sec rms. A fine guidance system is proposed to body-point the entire spacecraft within this limit. The spacecraft is modeled as a rigid body having reaction wheel actuators and subject to gravitational and magnetic disturbance torques. The fine guidance sensor is cluttered with electronic noise. The disturbance accommodation standard deviation optimal controller (DASOC) is designed to be optimal with respect to the transient and the steady state response to noise, whereas the steady state response to deterministic external torques is exactly zero. Compared with conventional controllers, the fine guidance system with the DASOC offers as much as a factor of 30 improvement in pointing stability, resulting in an optimal performance of nearly 0.0001 arc sec rms. Thus, the required pointing stability can easily be obtained, and a large margin remains for the compensation of possibile deteriorations.

Schiehlen, W. O.↗
NASA NTRS2012

The Gas Distribution in the Outer Regions of Galaxy Clusters

Aims. We present our analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We have exploited the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius, We stacked the density profiles to detect a signal beyond T200 and measured the typical density and scatter in cluster outskirts. We also computed the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compared our average density and scatter profiles with the results of numerical simulations. Results. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond approximately r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict density profiles that are too steep, whereas runs including additional physics and/ or treating gas clumping agree better with the observed gas distribution. We report high-confidence detection of a systematic difference between cool-core and non cool-core clusters beyond approximately 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond approximately r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only small differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the ENZO simulations. Conclusions. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Although their general behavior agrees more closely with the observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, AGN feedback, and consideration of gas clumping is required to construct realistic models of the outer regions of clusters.

Eckert, D.↗
DOE OSTI2023

Cosmological constraints from the Subaru Hyper Suprime-Cam year 1 shear catalogue lensing convergence probability distribution function

Here we utilize the probability distribution function (PDF) of normalized convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) year 1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude S 8 = σ 8 $\sqrt {Ω_{m}/0.3}$. The large-scale structure’s statistical properties are incompletely described by the traditional two-point statistics, motivating our investigation of the PDF—a complementary higher-order statistic. By defining the PDF over the standard-deviation-normalized convergence map, we are able to isolate the non-Gaussian information. We use tailored simulations to compress the data vector and construct a likelihood approximation. We mitigate the impact of survey and astrophysical systematics with cuts on smoothing scales, redshift bins, and data vectors. We find S 8 = $0.860^{+0.066}_{–0.109}$ from the PDF alone and S 8 = 0.798$^{+0.029}_{–0.042}$ from the combination of the PDF and power spectrum (68% confidential level (CL)). The PDF improves the power-spectrum-only constraint by about 10%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Search for heavy Majorana or Dirac neutrinos and right-handed $W$ gauge bosons in final states with charged leptons and jets in $pp$ collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search for heavy right-handed Majorana or Dirac neutrinos N R and heavy right-handed gauge bosons W R is performed in events with energetic electrons or muons, with the same or opposite electric charge, and energetic jets. The search is carried out separately for topologies of clearly separated final-state products (“resolved” channel) and topologies with boosted final states with hadronic and/or leptonic products partially overlapping and reconstructed as a large-radius jet (“boosted” channel). The events are selected from pp collision data at the LHC with an integrated luminosity of 139 fb -1 collected by the ATLAS detector at $\sqrt{s}=13$ TeV. No significant deviations from the Standard Model predictions are observed. The results are interpreted within the theoretical framework of a left-right symmetric model, and lower limits are set on masses in the heavy right handed W R boson and N R plane. The excluded region extends to about m(W R ) = 6.4 TeV for both Majorana and Dirac N R neutrinos at m(N R ) < 1 TeV. N R with masses of less than 3.5 (3.6) TeV are excluded in the electron (muon) channel at m(W R ) = 4.8 TeV for the Majorana neutrinos, and lim its of m(N R ) up to 3.6 TeV for m(W R ) = 5.2 (5.0) TeV in the electron (muon) channel are set for the Dirac neutrinos. These constitute the most stringent exclusion limits to date for the model considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The Pantheon+ Analysis: Evaluating Peculiar Velocity Corrections in Cosmological Analyses with Nearby Type Ia Supernovae

Separating the components of redshift due to expansion and peculiar motion in the nearby universe (z < 0.1) is critical for using Type Ia Supernovae (SNe Ia) to measure the Hubble constant (H 0 ) and the equation-of-state parameter of dark energy (w). Here, we study the two dominant "motions" contributing to nearby peculiar velocities: large-scale, coherent-flow (CF) motions and small-scale motions due to gravitationally associated galaxies deemed to be in a galaxy group. We use a set of 584 low-z SNe from the Pantheon+ sample, and evaluate the efficacy of corrections to these motions by measuring the improvement of SN distance residuals. We study multiple methods for modeling the large and small-scale motions and show that, while group assignments and CF corrections individually contribute to small improvements in Hubble residual scatter, the greatest improvement comes from the combination of the two (relative standard deviation of the Hubble residuals, Rel. SD, improves from 0.167 to 0.157 mag). We find the optimal flow corrections derived from various local density maps significantly reduce Hubble residuals while raising H 0 by ~0.4 km s -1 Mpc -1 as compared to using CMB redshifts, disfavoring the hypothesis that unrecognized local structure could resolve the Hubble tension. We estimate that the systematic uncertainties in cosmological parameters after optimally correcting redshifts are 0.06–0.11 km s -1 Mpc -1 in H 0 and 0.02–0.03 in w which are smaller than the statistical uncertainties for these measurements: 1.5 km s -1 Mpc -1 for H 0 and 0.04 for w.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

An initial state perturbation experiment with the GISS model

Monthly mean global forecasts for January 1975 have been computed with the Goddard Institute for Space Studies model from four slightly different sets of initial conditions - a 'control' state and three random perturbations thereof - to simulate the effects of initial state uncertainty on forecast quality. Differences among the forecasts are examined in terms of energetics, synoptic patterns and forecast statistics. The 'noise level' of the model predictions is depicted on global maps of standard deviations of sea level pressures, 500 mb heights and 850 mb temperatures for the set of four forecasts. Initial small-scale random errors do not appear to result in any major degradation of the large-scale monthly mean forecast beyond that generated by the model itself, nor do they appear to represent the major source of large-scale forecast error.

Spar, J.↗
NASA NTRS1978

Ocean wave heights measured by a high resolution pulse-limited radar altimeter

Data on significant wave height (SWH) taken with an airborne (2-3.5 km altitude) X-band 1-ns pulse-limited radar altimeter under various wind and sea conditions are interpreted. A heuristic discussion of the return pulse shape and the method of extracting the SWH is followed by a demonstration that the mean of the resulting SWH values is in agreement with other, independent measurements. The scatter of the SWH values is large compared to the estimates based on the statistical fluctuations in the radar signal. The discrepancy is resolved by showing that the radar is actually observing small-scale variations in the sea-surface standard deviation caused by the small number of ocean wavelengths illuminated.

Walsh, E. J.↗
NASA NTRS1983

Radiometer requirements for Earth-observation systems using large space antennas

Requirements are defined for Earth observation microwave radiometry for the decade of the 1990's by using large space antenna (LSA) systems with apertures in the range from 50 to 200 m. General Earth observation needs, specific measurement requirements, orbit mission guidelines and constraints, and general radiometer requirements are defined. General Earth observation needs are derived from NASA's basic space science program. Specific measurands include soil moisture, sea surface temperature, salinity, water roughness, ice boundaries, and water pollutants. Measurements are required with spatial resolution from 10 to 1 km and with temporal resolution from 3 days to 1 day. The primary orbit altitude and inclination ranges are 450 to 2200 km and 60 to 98 deg, respectively. Contiguous large scale coverage of several land and ocean areas over the globe dictates large (several hundred kilometers) swaths. Radiometer measurements are made in the bandwidth range from 1 to 37 GHz, preferably with dual polarization radiometers with a minimum of 90 percent beam efficiency. Reflector surface, root mean square deviation tolerances are in the wavelength range from 1/30 to 1/100.

Keafer, L. S., Jr.↗
NASA NTRS1990

Measuring Diameters Of Large Vessels

Computerized apparatus produces accurate results quickly. Apparatus measures diameter of tank or other large cylindrical vessel, without prior knowledge of exact location of cylindrical axis. Produces plot of inner circumference, estimate of true center of vessel, data on radius, diameter of best-fit circle, and negative and positive deviations of radius from circle at closely spaced points on circumference. Eliminates need for time-consuming and error-prone manual measurements.

Currie, James R.↗
NASA NTRS2018

Learning Model Structural Uncertainty with Gaussian Processes

The advent of commercially available quantum computers has marked the beginning of quantum computing as a reality. Both quantum gate and annealing computers have been released by major computer hardware companies. In this work, the D-Wave 2XTM quantum annealing computer housed at the NASA Advanced Systems computational facility is investigated to accelerate Machine Learning (ML) for image registration. NASA collects large amounts of images over the globe remotely using space-based monitoring. Images of a fixed areas of the land surface are taken over time. Due to the orbit of the sensors, the viewing angles deviate slightly, and it is necessary to align or register the images precisely to create image time series over the land surface. Unaligned images can lead to substantial analysis errors. These time-series are then used in modeling Earth Systems models such as hydrological, weather, and carbon monitoring models. In this work, we consider the Moderate Resolution Image Spectrometer (MODIS) data collected by the NASA's terra satellite. Artificial Neural Networks (ANNs) is a natural fit for ML modelling of images. Several successes have been reported using machine learning related to image processing. We investigate the use of ML to register MODIS images. ANNs are investigated in combination with a Restricted Boltzmann Machines (RBM) as an auto-encoder. We will present results showing the accuracy and efficiency of this approach.The D-Wave 2XTM quantum annealer samples the ground-state wave-function of a spin-Ising systems with quadratic interactions between qubits and a Chimera connectivity. The system sits in a ~15 mK thermal bath. One can think of the system as being placed in the ground state initially and subject to thermal excitations governed by Boltzmann statistics. If this is assumed true, one can use the statistics from the D-Wave 2XTM to train RBMs. Generating statistics for training Boltzmann machines is an NP-hard problem and constitutes the largest compute cost. We investigate the use of the D-Wave 2XTM to accelerate the training of the RBMs in our ANNs and report on the results.

Kouatchou, Jules↗
DOE OSTI2025

Magnitude and timescale of liquid water path adjustments to cloud droplet number concentration perturbations for nocturnal non-precipitating marine stratocumulus

Cloud liquid water path (L) adjusts to perturbations in cloud droplet number concentration (N) over time. We explore the magnitude and timescale of this adjustment in nocturnal non-precipitating marine stratocumuli using large eddy simulations of baseline conditions and aerosol seeding experiments for 22 meteorological conditions. The results confirm that the L adjustment (δL) slope (k) is more negative for simulation pairs with relatively low N and less negative for high N. Overall, k is unlikely to be lower than −0.4 within 24 h of seeding start, meaning the L adjustment is unlikely to fully offset the brightening due to the Twomey effect. After seeding, the δL becomes increasingly negative, which can be characterized by an exponential convergence. This evolution is governed by a short timescale around 5 h and lasts for around 8–12 h. It is driven by the feedback between entrainment, L, and boundary layer (BL) turbulence. Other processes, including radiation, surface fluxes, and subsidence, respond to the seeding weakly. This short timescale is insensitive to the amount of seeding, making the evolution of δL and some other deviations similar for different seeding amounts after appropriate scaling. The timescale of k evolution is closely related to the δL timescale and hence also short, while it could also be affected by the δN evolution. The results are most relevant to conditions where seeding is applied to a large area of marine stratocumulus in well-mixed and overcast BL where shear is not a primary source of turbulence.

54 ENVIRONMENTAL SCIENCES↗