Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian simulation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Ion Trapping Studies and Mitigation Strategies for the EIC ERL-Based Strong Hadron Cooler

An Energy Recovery Linac based strong hadron cooler was previously considered for the Electron-Ion Collider. The required electron beam parameters for variable-energy strong hadron cooling place significant constraints on ion trapping and collective effects. This paper presents initial studies of these constraints through a combination of analytical modelling and numerical simulations of ion production, trapping behaviour, and mitigation strategies. A multi-bunch tracking framework based on ELEGANT with the ionEffects module is used to simulate machine operation over millisecond time scales, corresponding to more than 3 × 10^5 electron bunches. The simulations include modelling of ionisation processes together with transverse electron–ion dynamics, allowing the evolution and accumulation of ions to be investigated. Analytical expressions based on Gaussian beam distributions are used to estimate ion trapping conditions and benchmark the simulation results. A bi-periodic bunch spacing scheme is also investigated as a possible mitigation method by detuning the ion oscillation frequency. These studies provide an initial assessment of ion trapping in the strong hadron cooler and demonstrate possible approaches for reducing beam–ion effects.

Bi, R. [Lancaster University, Cockcroft Institute]↗
DOE OSTI2024

A streamline coordinate analysis of a turbulent boundary layer subject to pressure gradients and curvature on the windward side of a bump

Direct numerical simulation (DNS) of a turbulent boundary layer over the Gaussian (Boeing) bump is performed. This boundary layer exhibits a series of adverse and favourable pressure gradients and convex and concave curvature effects before separating. These effects on turbulent boundary layers are characterised and compared with a lower-Reynolds-number flow over the same geometry. The momentum budgets are analysed in the streamline-aligned coordinate system upstream of the separation region. These momentum budgets allow the simplification of equations to facilitate an integral analysis. Integral-analysis-based approximations for Reynolds stresses in the inner and outer regions of the boundary layer are also formulated. The shear and wall-normal Reynolds stress profiles normalised by these approximations exhibit a better collapse compared with friction velocity and Zagarola–Smits normalisations in the strong favourable pressure gradient region and in the mild adverse pressure region that precedes it in this flow. Simplification of these Reynolds stress approximations along with results from the DNS are used to obtain semi-empirical approximations that are able to provide stress closure in terms of wall solution fields for the turbulent boundary layer under consideration.

42 ENGINEERING↗
DOE OSTI2022

Investigation of the Failure of Marcus Theory for Hydrated Electron Reactions

Reactions of the hydrated electron with a wide variety of substrates have been found to exhibit unusually similar activation energies in a manner incompatible with Marcus electron transfer theory. Given the fundamental linear response assumption of Marcus theory, one possible explanation for this apparent failure is that the underlying free energy surfaces governing the reactions are not harmonic, i.e., hydrated electron structural fluctuations exhibit non-Gaussian behavior. In this work, we test this hypothesis by using simulations to calculate the hydrated electron vertical detachment energy distribution. Furthermore, we consider both cavity and non-cavity models for the hydrated electron, between which the actual hydrated electron behavior is expected to lie. Our results identify a possible origin for non-Gaussian behavior of the hydrated electron, but show that it is not of sufficient magnitude to explain the failure of Marcus theory to describe its reactions. Thus, other explanations must be sought.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Inversion of Magnetic Measurements of the CHAMP Satellite Over the Pannonian Basin

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The anomalies were distributed in a spherical shell, some 107,927 data recorded between January 1 and December 31 of 2008. They covered the Pannonian Basin and its vicinity. These anomaly data were interpolated into a spherical grid of 0.5 x 0.5, at the elevation of 324 km by the Gaussian weight function. The vertical gradient of these total magnetic anomalies was also computed and mapped to the surface of a sphere at 324 km elevation. The former spherical anomaly data at 425 km altitude were downward continued to 324 km. To interpret these data at the elevation of 324 km we used an inversion method. A polygonal prism forward model was used for the inversion. The minimum problem was solved numerically by the Simplex and Simulated annealing methods; a L2 norm in the case of Gaussian distribution parameters and a L1 norm was used in the case of Laplace distribution parameters. We INTERPRET THAT the magnetic anomaly WAS produced by several sources and the effect of the sable magnetization of the exsolution of hemo-ilmenite minerals in the upper crustal metamorphic rocks.

Kis, K. I.↗
DOE OSTI2022

KaRMMa – kappa reconstruction for mass mapping

ABSTRACT We present KaRMMa, a novel method for performing mass map reconstruction from weak-lensing surveys. We employ a fully Bayesian approach with a physically motivated lognormal prior to sample from the posterior distribution of convergence maps. We test KaRMMa on a suite of dark matter N-body simulations with simulated DES Y1-like shear observations. We show that KaRMMa outperforms the basic Kaiser–Squires mass map reconstruction in two key ways: (1) our best map point estimate has lower residuals compared to Kaiser–Squires; and (2) unlike the Kaiser–Squires reconstruction, the posterior distribution of KaRMMa maps is nearly unbiased in all summary statistics we considered, namely: one-point and two-point functions, and peak/void counts. In particular, KaRMMa successfully captures the non-Gaussian nature of the distribution of κ values in the simulated maps. We further demonstrate that the KaRMMa posteriors correctly characterize the uncertainty in all summary statistics we considered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dark Energy Survey Year 3 results: Simulation-based 𝑤CDM inference from weak lensing and galaxy clustering maps with deep learning: Analysis design

Data-driven approaches using deep learning are emerging as powerful techniques to extract non-Gaussian information from cosmological large-scale structure. Here, this work presents the first simulation-based inference (SBI) pipeline that combines weak lensing and galaxy clustering maps in a realistic Dark Energy Survey Year 3 (DES Y3) configuration and serves as preparation for a forthcoming analysis of the survey data. We develop a scalable forward model based on the CosmoGridV1 suite of N-body simulations to generate over one million self-consistent mock realizations of DES Y3 at the map level. Leveraging this large dataset, we train deep graph convolutional neural networks on the full survey footprint in spherical geometry to learn low-dimensional features that approximately maximize mutual information with target parameters. These learned compressions enable neural density estimation of the implicit likelihood via normalizing flows in a ten-dimensional parameter space spanning cosmological 𝑤CDM, intrinsic alignment, and linear galaxy bias parameters, while marginalizing over baryonic, photometric redshift, and shear bias nuisances. To ensure robustness, we extensively validate our inference pipeline using synthetic observations derived from both systematic contaminations in our forward model and independent Buzzard galaxy catalogs. Our forecasts yield significant improvements in cosmological parameter constraints, achieving 2−3× higher figures of merit in the 𝛺 𝑚 − 𝑆 8 plane relative to our implementation of baseline two-point statistics and effectively breaking parameter degeneracies through probe combination. These results demonstrate the potential of SBI analyses powered by deep learning for upcoming Stage-IV wide-field imaging surveys.

Thomsen, A. [Zurich, ETH] (ORCID:0000000203099021)↗
DOE OSTI2022

Measuring the thermal and ionization state of the low- z IGM using likelihood free inference

ABSTRACT We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho/ \bar{\rho })^{\gamma -1}$ and the UV background photoionization rate $\Gamma _{{{{\rm H\, {\small I}}}}{}}$ of the intergalactic medium (IGM) based on the Voigt profile decomposition of the Ly α forest into a set of discrete absorption lines with Doppler parameter b and the neutral hydrogen column density $N_{\rm H\, {\small I}}$. Previous work demonstrated that the shape of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution is sensitive to the IGM thermal parameters T0 and γ, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density dN/dz, and we demonstrate that precise constraints can also be obtained on $\Gamma _{{{{\rm H\, {\small I}}}}{}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution on IGM parameters trained on an ensemble of 624 nyx hydrodynamical simulations at z = 0.1, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS data sets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly α forest spectra at z ≃ 0.1, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma _{\log T_0} \sim 0.08$ dex, σγ ∼ 0.06, and $\sigma _{\log \Gamma _{{{{\rm H\, {\small I}}}}{}}} \sim 0.07$ dex).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum Sensing of Displacements with Stabilized Gottesman-Kitaev-Preskill States

We demonstrate how recent protocols developed for the stabilization of Gottesman-Kitaev-Preskill states can be used for the estimation of two-quadrature displacement sensing, with sensitivities approaching the multivariate quantum Cramer-Rao bound. Thanks to the stabilization, this sensor is backaction evading and can function continuously without reset, making it well suited for the detection of itinerant signals. Additionally, we provide numerical simulations showing that the protocol can unconditionally surpass the Gaussian limit of displacement sensing with prior information, even in the presence of realistic noise. Our work shows how reservoir engineering in bosonic systems can be leveraged for quantum metrology, with potential applications in force sensing, waveform estimation, and quantum channel learning.

Labarca, Lautaro [Univ. of Sherbrooke, QC (Canada)↗
DOE OSTI2023

Maximizing the impact of the Vera Rubin Observatory LSST, and DESI and CMB-S4 experiments in constraining gravity and dark sector physics

The work in this award has addressed our objectives on a number of fronts: Current graduate student Chris Wilson and Bean are determining how cosmic voids can be used to test the properties of gravity, leveraging simulations from DESI collaborator Baojiu Li and refining theoretical Gaussian Streaming Model approaches to predict galaxy dynamics in voids when the cosmology is not GR. This work is still underway, including analyses of the efficacy of void quadrupole correlation statistics and void size functions when studying dark matter halo statistics or the expected galaxy distributions from upcoming galaxy surveys; Current graduate student Yulin Gong, Bean with collaborators are refining techniques to extract out the kinetic Sunyaev Zel’dovich effect from galaxy cluster survey data and CMB temperature anisotropies, with the intent to further ready this pipeline development in readiness for spectroscopic data from DESI and LSST along with CMB-S4 measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Delta Modulation

An optimum adaptive delta modulator-demodulator configuration is derived. This device utilizes two past samples to obtain a step size which minimizes the mean square error for a Markov Gaussian source. The optimum system is compared using computer simulations with the linear delta modulator and an enhanced Abate delta modulator. In addition the performance is compared to the rate distortion bound for a Markov source. It is shown that the optimum delta modulator is neither quantization nor slope-overload limited. The results show that the output signal-to-noise ratio is independent of the input signal power and is subject only to the limitations of the hardware employed. In addition, voice was recorded using these systems. The demodulated voice indicates negligible degradation is caused by the optimum system and by the enhanced Abate system while the linear delta modulator suffers significant degradation at a sampling frequency of 56 kilobits/sec.

Schilling, D. L.↗
NASA NTRS1971

A variable-step-size robust delta modulator.

Description of an analytically obtained optimum adaptive delta modulator-demodulator configuration. The device utilizes two past samples to obtain a step size which minimizes the mean square error for a Markov-Gaussian source. The optimum system is compared, using computer simulations, with a linear delta modulator and an enhanced Abate delta modulator. In addition, the performance is compared to the rate distortion bound for a Markov source. It is shown that the optimum delta modulator is neither quantization nor slope-overload limited. The highly nonlinear equations obtained for the optimum transmitter and receiver are approximated by piecewise-linear equations in order to obtain system equations which can be transformed into hardware. The derivation of the experimental system is presented.

Song, C. L.↗
NASA NTRS1972

A TDM link with channel coding and digital voice.

The features of a TDM (time-division multiplexed) link model are described. A PCM telemetry sequence was coded for error correction and multiplexed with a digitized voice channel. An all-digital implementation of a variable-slope delta modulation algorithm was used to digitize the voice channel. The results of extensive testing are reported. The measured coding gain and the system performance over a Gaussian channel are compared with theoretical predictions and computer simulations. Word intelligibility scores are reported as a measure of voice channel performance.

Jones, M. W.↗
NASA NTRS1986

Scattering from mirrors contaminated by particulates - A model

A model has been used to fit the bidirectional reflectance distribution function for a variety of samples with different angles of incidence and degrees of contamination and roughness by the proper adjustment of six parameters. Exponential and Gaussian composite surface height autocorrelation functions are used to simulate the data at the small and medium angles, and a perturbation of these functions with two additional parameters allows the fit of the upward and downward curve portions of the data at the large angles. Good agreement is found with previous results. The present approach is the same as the classical dispersion equation in treating the index of refraction of optical materials.

Wang, Yaujen↗
NASA NTRS2002

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows

The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.

He, Guowei↗
NASA NTRS1999

Tilt Modulation Distortions in Wave Topography Measured by a Scanning Radar Altimeter

Work continues on estimating tilt modulation distortions in the wave topography measured by a scanning radar altimeter. To quantify this effect, a two-dimensional simulation has been performed in the cross-track plane only which assumed that sinusoidal waves of constant wavelength propagate in the cross-track direction (and are infinitely long-crested in the along-track direction). The initial results reported earlier for a Gaussian surface scattering model indicated that when the highest reasonable value of mss is used in the simulation (the Plant limit of 0.08), the nadir values of the ratio of the apparent to actual wave height are the same as for the omnidirectional scattering case. But as the off-nadir angle increased, the apparent wave height increased and became larger than the actual wave height by about 10 degrees off-nadir. And the shorter the wavelength, the larger the apparent wave height increase. This represented a systematic over-estimate of the wave amplitude for waves propagating in the cross-track direction. For lower values of mss the situation worsened. The simulation has been improved by incorporating actual variations of backscattered power with incidence angle measured by the SRA instead of the Gaussian model. The resulting distortion was about half that originally reported. The 3-dimensional simulation to model waves propagating at various azimuthal angles relative to the cross-track plane is still in progress. The results of this model will be verified by comparing them with Scanning Radar Altimeter (SRA) data and optimum correction procedures will be developed. An assessment of the improved 2-dimensional model indicates that the distortion will generally be small in the data taken during the Southern Ocean Waves Experiment (SOWEX) presently be analyzed. But tilt modulation effects are of great concern for SRA data collected during the 1998 hurricane season since the minimum aircraft altitude was 1.5 km and it frequently flew higher. For a given wave length and steepness, the effect increases with increasing aircraft altitude, and is larger for waves propagating in the cross-track direction, which was generally the situation during the hurricane flights.

Walsh, Edward J.↗
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗
DOE OSTI2023

Primordial non-Gaussianity with angular correlation function: integral constraint and validation for DES

ABSTRACT Local primordial non-Gaussianity (PNG) is a promising observable of the underlying physics of inflation, characterized by $f_{\rm NL}^{\rm loc}$. We present the methodology to measure $f_{\rm NL}^{\rm loc}$ from the Dark Energy Survey (DES) data using the two-point angular correlation function (ACF) with scale-dependent bias. One of the focuses of the work is the integral constraint (IC). This condition appears when estimating the mean number density of galaxies from the data and is key in obtaining unbiased $f_{\rm NL}^{\rm loc}$ constraints. The methods are analysed for two types of simulations: ∼246 goliat-png N-body small area simulations with fNL equal to −100 and 100, and 1952 Gaussian ice-cola mocks with fNL = 0 that follow the DES angular and redshift distribution. We use the ensemble of goliat-png mocks to show the importance of the IC when measuring PNG, where we recover the fiducial values of fNL within the 1σ when including the IC. In contrast, we found a bias of ΔfNL ∼ 100 when not including it. For a DES-like scenario, we forecast a bias of ΔfNL ∼ 23, equivalent to 1.8σ, when not using the IC for a fiducial value of fNL = 100. We use the ice-cola mocks to validate our analysis in a realistic DES-like set-up finding it robust to different analysis choices: best-fitting estimator, the effect of IC, BAO damping, covariance, and scale choices. We forecast a measurement of fNL within σ(fNL) = 31 when using the DES-Y3 BAO sample, with the ACF in the 1 deg < θ < 20 deg range.

79 ASTRONOMY AND ASTROPHYSICS↗