Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distance scale”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Origin of micrometer-scale dislocation motion during hydrogen desorption

Hydrogen, while being a potential energy solution, creates arguably the most important embrittlement problem in high-strength metals. However, the underlying hydrogen-defect interactions leading to embrittlement are challenging to unravel. Here, we investigate an intriguing hydrogen effect to shed more light on these interactions. By designing an in situ electron channeling contrast imaging experiment of samples under no external stresses, we show that dislocations (atomic-scale line defects) can move distances reaching 1.5 μm during hydrogen desorption. Combining molecular dynamics and grand canonical Monte Carlo simulations, we reveal that grain boundary hydrogen segregation can cause the required long-range resolved shear stresses, as well as short-range atomic stress fluctuations. Thus, such segregation effects should be considered widely in hydrogen research.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

The Utility of Satellites and Autonomous Remote Sensing Platforms for Monitoring Offshore Aquaculture Farms: A Case Study for Canopy Forming Kelps

The emerging sector of offshore kelp aquaculture represents an opportunity to produce biofuel feedstock to help meet growing energy demand. Giant kelp represents an attractive aquaculture crop due to its rapid growth and production, however precision farming over large scales is required to make this crop economically viable. These demands necessitate high frequency monitoring to ensure outplant success, maximum production, and optimum quality of harvested biomass, while the long distance from shore and large necessary scales of production makes in person monitoring impractical. Remote sensing offers a practical monitoring solution and nascent imaging technologies could be leveraged to provide daily products of the kelp canopy and subsurface structures over unprecedented spatial scales. Here, we evaluate the efficacy of remote sensing from satellites and aerial and underwater autonomous vehicles as potential monitoring platforms for offshore kelp aquaculture farms. Decadal-scale analyses of the Southern California Bight showed that high offshore summertime cloud cover restricts the ability of satellite sensors to provide high frequency direct monitoring of these farms. By contrast, daily monitoring of offshore farms using sensors mounted to aerial and underwater drones seems promising. Small Unoccupied Aircraft Systems (sUAS) carrying lightweight optical sensors can provide estimates of canopy area, density, and tissue nitrogen content on the time and space scales necessary for observing changes in this highly dynamic species. Underwater color imagery can be rapidly classified using deep learning models to identify kelp outplants on a longline farm and high acoustic returns of kelp pneumatocysts from side scan sonar imagery signal an ability to monitor the subsurface development of kelp fronds. Current sensing technologies can be used to develop additional machine learning and spectral algorithms to monitor outplant health and canopy macromolecular content, however future developments in vehicle and infrastructure technologies are necessary to reduce costs and transcend operational limitations for continuous deployment in an offshore setting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The role of geography in the complex diffusion of innovations

The urban–rural divide is increasing in modern societies calling for geographical extensions of social influence modelling. Improved understanding of innovation diffusion across locations and through social connections can provide us with new insights into the spread of information, technological progress and economic development. In this work, we analyze the spatial adoption dynamics of iWiW, an Online Social Network (OSN) in Hungary and uncover empirical features about the spatial adoption in social networks. During its entire life cycle from 2002 to 2012, iWiW reached up to 300 million friendship ties of 3 million users. We find that the number of adopters as a function of town population follows a scaling law that reveals a strongly concentrated early adoption in large towns and a less concentrated late adoption. We also discover a strengthening distance decay of spread over the life-cycle indicating high fraction of distant diffusion in early stages but the dominance of local diffusion in late stages. The spreading process is modelled within the Bass diffusion framework that enables us to compare the differential equation version with an agent-based version of the model run on the empirical network. Although both model versions can capture the macro trend of adoption, they have limited capacity to describe the observed trends of urban scaling and distance decay. We find, however that incorporating adoption thresholds, defined by the fraction of social connections that adopt a technology before the individual adopts, improves the network model fit to the urban scaling of early adopters. Controlling for the threshold distribution enables us to eliminate the bias induced by local network structure on predicting local adoption peaks. Finally, we show that geographical features such as distance from the innovation origin and town size influence prediction of adoption peak at local scales in all model specifications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Figuring Out Gas & Galaxies in Enzo (FOGGIE). VI. The Circumgalactic Medium of L * Galaxies Is Supported in an Emergent, Nonhydrostatic Equilibrium

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium, with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly resolved L* galaxies from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM (Δx ≲ 1 kpc h –1 and M ≃ 200M ⊙ ). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of ≲5 kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond ~0.25R 200 , or ~50 kpc at z = 0. Within ~0.25R 200 , turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the nonequilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies, even when averaging out small-scale variations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Role of the Hubble scale in the weak lensing versus CMB tension

Here, we explore a reparametrization of the lensing amplitude tension between weak lensing (WL) and cosmic microwave background (CMB) data and its implications for a joint resolution with the Hubble tension. Specifically, we focus on the lensing amplitude over a scale of 12 Mpc in absolute distance units using a derived parameter S 12 and show its constraints from recent surveys in comparison with Planck 2018. In WL alone, we find that the absolute distance convention correlates S 12 with H 0 . Accounting for this correlation in the 3D space S 12 ×w m ×ℎ reproduces the usual levels of 2~3⁢σ tension inferred from S 8 ×Ω m . Additionally, we derive scaling relations in the S 8 ×ℎ and S 12 ×ℎ planes that are allowed by Λ⁢CDM and extrapolate target scalings needed to solve the H 0 and lensing-amplitude tensions jointly in a hypothetical beyond-Λ⁢CDM model. As a test example, we quantify how the early dark energy scenario compares with these target scalings. Useful fitting formulas for S 8 and S 12 as a function of other cosmological parameters in Λ⁢CDM are provided, with 1% precision.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Exploring QCD Factorization at Moderate Energy Scales

Asymptotic freedom in QCD facilitates the use of partonic degrees of freedom over short distances, but physical processes are sensitive to a wide range of scales. Thus, it is necessary in QCD calculations to utilize a factorization scheme to separate a process into perturbative and non-perturbative factors. This separation relies on an assumption that one energy scale is infinitely larger than the other scales involved in the process. However, much experimental research in areas such as nucleon structure and quark-hadron duality occur at more moderate energy scales where that basic assumption may not be true but perturbative calculations should still be useful. Thus, an exploration of the limits of factorization at more moderate energy scales is needed. This dissertation examines various aspects of factorization at these energy scales first by applying the necessary approximations to a simple model where exact calculations are possible and so the effects of these approximations can be quantified. This is followed by examining areas where corrections are known to be needed. First is an exploration of target mass corrections (TMCs) in the case of deep inelastic scattering (DIS), including a discussion of what large corrections imply about the target structure. Second, is a general examination of PDFs and FFs at moderate scales. Third, I will discuss how this fits into the long-term effort to study the transition between small transverse momentum (generated non-perturbatively) and large transverse momentum (generated in the hard process)

Moffat, Eric↗
DOE OSTI2023

The Alcock–Paczynski effect from Lyman- α forest correlations: analysis validation with synthetic data

The three-dimensional distribution of the Ly α forest has been extensively used to constrain cosmology through measurements of the baryon acoustic oscillations (BAO) scale. However, more cosmological information could be extracted from the full shapes of the Ly α forest correlations through the Alcock–Paczynski (AP) effect. In this work, we prepare for a cosmological analysis of the full shape of the Ly α forest correlations by studying synthetic data of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We use a set of 100 eBOSS synthetic data sets in order to validate such an analysis. These mocks undergo the same analysis process as the real data. We perform a full-shape analysis on the mean of the correlation functions measured from the 100 eBOSS realizations, and find that our model of the Ly α correlations performs well on current data sets. We show that we are able to obtain an unbiased full-shape measurement of D M /D H (z eff ), where D M is the transverse comoving distance, D H is the Hubble distance, and z eff is the effective redshift of the measurement. We test the fit over a range of scales, and decide to use a minimum separation of r min = 25 h –1 Mpc. Here, we also study and discuss the impact of the main contaminants affecting Ly α forest correlations, and give recommendations on how to perform such analysis with real data. While the final eBOSS Ly α BAO analysis measured D M /D H (z eff = 2.33) with 4 per cent statistical precision, a full-shape fit of the same correlations could provide an $\sim 2~{{\ \rm per\ cent}}$ measurement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Asymptotic scalar field cosmology in string theory

Asymptotic (late-time) cosmology depends on the asymptotic (infinite-distance) limits of scalar field space in string theory. Such limits feature an exponentially decaying potential V ~ exp(–c$\phi$) with corresponding Hubble scale H ~ $\sqrt{{\dot{\phi}}^2+2V}$ ~ exp(–λ H$\phi$ ), and at least one tower of particles whose masses scale as m ~ exp(–λΦ), as required by the Distance Conjecture. In this paper, we provide evidence that these coefficients satisfy the inequalities $\sqrt{\left(d-1\right)/\left(d-2\right)}$ ≥ λ H ≥ λ lightest ≥ 1/$\sqrt{d-2}$ in d spacetime dimensions, where λ lightest is the λ coefficient of the lightest tower. This means that at late times, as the scalar field rolls to $\phi$ → ∞, the low-energy theory remains a d-dimensional FRW cosmology with decelerated expansion, the light towers of particles predicted by the Distance Conjecture remain at or above the Hubble scale, and both the strong energy condition and the dominant energy condition are satisfied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Electron Transfer Going the Distance: Mn-Doped ZnSe as a Model Photocatalytic System

In this work, time-resolved photoluminescence and transient absorption spectroscopy are used to examine the influence that long excited state lifetimes of Mn-doped nanocrystals have on the mechanisms of photo-induced electron transfer (PET) in photocatalytic systems. Mn-doped ZnSe nanocrystals can undergo PET over their nearly millisecond excited state lifetime, which enables electron transfer to viologens in solution over large encounter distances in comparison to molecular length scales. The long excited state lifetimes also enable diffusion of the excited nanocrystals over hundreds of nanometers in solution, allowing them to react with molecular species at nanomolar concentrations. The ability to capture and sustain optical energy in spin-forbidden transitions among crystal field states of the Mn ions opens opportunities to explore the influence that long excited state lifetimes have on photocatalytic reaction mechanisms involving molecular species such as CO 2 that can be difficult to concentrate or attach to nanocrystal surfaces.

14 SOLAR ENERGY↗
DOE OSTI2021

Pressure-induced spin crossover in a Fe 78 Si 9 B 13 metallic glass

The pressure effect on structures and properties of a Fe 78 Si 9 B 13 metallic glass was investigated by in situ high-pressure synchrotron Fe Kβ x-ray emission spectroscopy and x-ray diffraction, and electrical resistivity measurements up to ~51 GPa. The study reveals a reversible and continuous pressure-induced high- to low-spin crossover of Fe atoms in an amorphous structure. The changes of the local spin moment can be scaled to match its average atomic distance shrinkage very well during compression. The crossover of electronic spin states in the Fe 78 Si 9 B 13 metallic glass resembles that of typical crystalline Fe-bearing materials but without a sharp atomic volume collapse and an abrupt electrical resistivity jump. Furthermore, these findings could help guide applications of Fe-based metallic glasses as a soft ferromagnetic material at extreme conditions and also improve our understanding of magnetism and coupling of its changes with disordered atomic structures and other properties in metallic glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

QCD Running Coupling in the Nonperturbative and Near-Perturbative Regimes

We use analytic continuation to extend the gauge-gravity duality nonperturbative description of the strong force coupling into the transition, near-perturbative, regime where perturbative effects become important. By excluding the unphysical region in coupling space from the flow of singularities in the complex plane, we derive a specific relation between the scales relevant at large and short distances; this relation is uniquely fixed by requiring maximal analyticity. The unified effective coupling model gives an accurate description of the data in the nonperturbative and the near-perturbative regions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Using a Neural Network Classifier to Select Galaxies with the Most Accurate Photometric Redshifts

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will produce several billion photometric redshifts (photo-z's), enabling cosmological analyses to select a subset of galaxies with the most accurate photo-z. For this work, we perform initial redshift fits on Subaru Strategic Program galaxies with deep grizy photometry using Trees for Photo-Z (TPZ) before applying a custom neural network classifier (NNC) tuned to select galaxies with (z phot – z spec )/(1 + z spec ) < 0.10. We consider four cases of training and test sets ranging from an idealized case to using data augmentation to increase the representation of dim galaxies in the training set. Selections made using the NNC yield significant further improvements in outlier fraction and photo-z scatter (σ z ) over those made with typical photo-z uncertainties. As an example, when selecting the best third of the galaxy sample, the NNC achieves a 35% improvement in outlier rate and a 23% improvement in σ z compared to using uncertainties from TPZ. For cosmology and galaxy evolution studies, this method can be tuned to retain a particular sample size or to achieve a desired photo-z accuracy; our results show that it is possible to retain more than a third of an LSST-like galaxy sample while reducing σ z by a factor of 2 compared to the full sample, with one-fifth as many photo-z outliers. For surveys like LSST that are not limited by shot noise, this method enables a larger number of tomographic redshift bins and hence a significant increase in the total signal to noise of galaxy angular power spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Impact of melt viscosity on filler dispersion in elastomeric nanocomposites

Compounding of commercial nanocomposites usually involves the addition of viscosity enhancers such as binder resins in ink jet inks, and paints. Contrary to this, plasticizers such as process oils are added to reduce the melt viscosity and ease processability of reinforced elastomers. Nanofillers such as silica and carbon black are typically added to reinforce rubber and enhance performance of automotive tire treads. Filler dispersion has traditionally been qualitatively (indirectly) assessed by measuring the properties of reinforced elastomers. While dispersion can be quantified by examining filler agglomeration through surface roughness measurements and microscopy, the size-scale dependence for these hierarchical fillers has usually been ignored. Here, we have recently devised a method to quantify nano-scale dispersion of fillers using Ultra small-angle X-ray scattering (USAXS) techniques. This method is advantageous since it directly links the controllable processing/compounding parameters such as the mixing speed, mixer geometry, residence time (or mixing duration), melt density, flow gap distance, and melt viscosity to nano-scale dispersion. While our previous studies have explored the impact of different processing parameters, this study specifically investigates the impact of melt viscosity on nano-scale filler dispersion in elastomer compounds. Commercially available polybutadienes with different Mooney viscosities were used in conjunction with different grades and amounts of process oils to modify the melt viscosity.

Carbon black↗
DOE OSTI2022

Universal van der Waals Force between Heavy Polarons in Superfluids

Here, we investigate the long-range behavior of the induced Casimir interaction between two spinless heavy impurities, or polarons, in superfluid cold atomic gases. With the help of effective field theory (EFT) of a Galilean invariant superfluid, we show that the induced impurity-impurity potential at long distance universally shows a relativistic van der Waals-like attraction (~1/r 7 ) resulting from the exchange of two superfluid phonons. We also clarify finite temperature effects from the same two-phonon exchange process. The temperature $\mathrm{T}$ introduces the additional length scale c s /$\mathrm{T}$ with the speed of sound c s . Leading corrections at finite temperature scale as $\mathrm{T}$ 6 /r for distances r$\ll$c s /$\mathrm{T}$ smaller than the thermal length. For larger distances the potential shows a nonrelativistic van der Waals behavior (~$\mathrm{T}$/r 6 ) instead of the relativistic one. Our EFT formulation applies not only to weakly coupled Bose or Fermi superfluids but also to those composed of strongly correlated unitary fermions with a weakly coupled impurity. The sound velocity controls the magnitude of the van der Waals potential, which we evaluate for the fermionic superfluid in the BCS-BEC crossover.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Parallel Simulation of Quantum Networks with Distributed Quantum State Management

Quantum network simulators offer the opportunity to cost-efficiently investigate potential avenues for building networks that scale with the number of users, communication distance, and application demands by simulating alternative hardware designs and control protocols. Several quantum network simulators have been recently developed with these goals in mind. As the size of the simulated networks increases, however, sequential execution becomes time-consuming. Parallel execution presents a suitable method for scalable simulations of large-scale quantum networks, but the unique attributes of quantum information create unexpected challenges. In this work, we identify requirements for parallel simulation of quantum networks and develop the first parallel discrete-event quantum network simulator by modifying the existing serial simulator SeQUeNCe. Our contributions include the design and development of a quantum state manager (QSM) that maintains shared quantum information distributed across multiple processes. We also optimize our parallel code by minimizing the overhead of the QSM and decreasing the amount of synchronization needed among processes. Using these techniques, we observe a speedup of 2 to 25 times when simulating a 1,024-node linear network topology using 2 to 128 processes. We also observe an efficiency greater than 0.5 for up to 32 processes in a linear network topology of the same size and with the same workload. We repeat this evaluation with a randomized workload on a caveman network. We also introduce several methods for partitioning networks by mapping them to different parallel simulation processes. We have released the parallel SeQUeNCe simulator as an open source tool alongside the existing sequential version.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗