SEARCH · Engineering Papers

Results for “Statistical 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,549 records · Page 86

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Observation of a prethermal discrete time crystal

Extending the framework of statistical physics to the nonequilibrium setting has led to the discovery of previously unidentified phases of matter, often catalyzed by periodic driving. However, preventing the runaway heating that is associated with driving a strongly interacting quantum system remains a challenge in the investigation of these newly discovered phases. In this work, we utilize a trapped-ion quantum simulator to observe the signatures of a nonequilibrium driven phase without disorder—the prethermal discrete time crystal. Here, the heating problem is circumvented not by disorder-induced many-body localization, but rather by high-frequency driving, which leads to an expansive time window where nonequilibrium phases can emerge. Floquet prethermalization is thus presented as a general strategy for creating, stabilizing, and studying intrinsically out-of-equilibrium phases of matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Computational analysis of deployable wind turbine systems in defense operational energy applications.

The U.S. military has been exploring pathways to reduce the logistical burden of fuel on virtually all their missions globally. Energy harvesting of local resources such as wind and solar can help increase the resilience and operational effectiveness of military units, especially at the most forward operating bases where the fuel logistics are most challenging. This report considers the potential benefits of wind energy provided by deployable wind turbines as measured by a reduction in fuel consumption and supply convoys to a hypothetical network of Army Infantry Brigade Combat Team bases. Two modeling and simulation tools are used to represent the bases and their operations and quantify the impacts of system design variables that include wind turbine technologies, battery storage, number of turbines, and wind resource quality. The System of Systems Analysis Toolkit Joint Operational Energy Model serves as a baseline scenario for comparison. The Hybrid Optimization of Multiple Energy Resources simulation tool is used to optimize a single base within the larger Joint Operational Energy Model. The results of both tools show that wind turbines can provide significant benefits to contingency bases in terms of reduced fuel use and number of convoy trips to resupply the base. The match between the turbine design and wind resource, which is statistically low across most of the global land area, is a critical design consideration. The addition of battery storage can enhance the benefits of wind turbines, especially in systems with more wind turbines and higher wind resources. Wind turbines may also provide additional benefits to other metrics such as resilience that may be important but not fully considered in the current analysis. ACKNOWLEDGEMENTS The authors would like to thank the following individuals for their helpful support, feedback and review to improve this report: U.S. Department of Energy Wind Energy Technologies Office, Patrick Gilman and Bret Barker; Idaho National Laboratory, Jake Gentle and Bradley Whipple; The National Renewable Energy Laboratory, Robert Preus and Tony Jimenez; Sandia National Laboratories, Alan Nanco, Dennis Anderson, and Hai Le. In addition, numerous discussions with military and industry stakeholders over the year were invaluable in focusing the efforts represented in this report.

17 WIND ENERGY↗

DOE OSTI2020

Effect of Power Plant Capacity on the CAPEX, OPEX, and LCOC of the CO2 Capture Process in Pre-Combustion Applications - Abstract

Aspen Plus v8.8 was used to simulate the pre-combustion CO2 capture process from a typical fuel gas stream at different power plant capacities ranging from 54.3 to 543 MW. Polyethylene glycol polydimethyl siloxane (PEGPDMS-1) was used as a physical solvent to capture CO2 in a countercurrent packed-bed absorber containing a structured packing (Mellapak 250Y). The process pressure was 51.4 bar and the solvent temperature was varied from 10 to 40 oC. The internal diameter of the absorber ranged from 1.9 to 7 m and the packing height from 13.9 to 45.5 m. The physico-chemical properties were obtained and modeled using the Perturbed Chain-Statistical Associating Fluid Theory (PC-SAFT) Equation-of-State (EOS). Four process constraints were imposed in Aspen Plus: (1) no flooding in the absorber, (2) the absorber height to diameter ratio (H/D) is greater than or equal 6, (3) at least 90 mol% of the CO2 capture from the fuel gas stream, and (4) the CO2 stream destined to sequestration sites should contain less than 600 ppm water concentration and less than or equal 0.5 mol% of fuel gases (H2, CO, CH4). The plant lifetime was assumed to be 30 years with an annual discount rate of 10% and an annual maintenance cost of 4% of the total capital expenses. Seven power plants with different capacities were used and for each plant, the corresponding CO2 capture process was simulated. The absorber flooding was checked using the generalized pressure drop correlation (GPDC) by Leva [1] and the capital expenditure (CAPEX), operating expenditure (OPEX), and levelized costs per ton of CO2 captured (LCOC) were calculated [2-5]. Normalized by the largest power plant capacity (543 MW) used in the simulation, the results indicated that the capital and operating expenditures increased, however, the LCOC decreased with increasing plant capacity. This behavior was due to the increased annual tonnage of the CO2 captured with plant capacity as shown in Figure 1. The calculated CAPEX, OPEX, LCOC and the tonnage of CO2 captured for the 543 MW power plant were about 52 MM$, 22 MM$/year, 7.46 $/ton and 4 MM ton/year.

Ashkanani, Husain↗

DOE OSTI2023

Constructing a Simulation Surrogate with Partially Observed Output

Gaussian process surrogates are a popular alternative to directly using computationally expensive simulation models. When the simulation output consists of many responses, dimension-reduction techniques are often employed to construct these surrogates. However, surrogate methods with dimension reduction generally rely on complete output training data. This article proposes a new Gaussian process surrogate method that permits the use of partially observed output while remaining computationally efficient. The new method involves the imputation of missing values and the adjustment of the covariance matrix used for Gaussian process inference. The resulting surrogate represents the available responses, disregards the missing responses, and provides meaningful uncertainty quantification. In conclusion, the proposed approach is shown to offer sharper inference than alternatives in a simulation study and a case study where an energy density functional model that frequently returns incomplete output is calibrated.

42 ENGINEERING↗

DOE OSTI2022

Direct numerical simulation of hypersonic turbulent boundary layers: effect of spatial evolution and Reynolds number

Direct numerical simulations (DNS) are performed to investigate the spatial evolution of flat-plate zero-pressure-gradient turbulent boundary layers over long streamwise domains ( >300δ i , with δ i the inflow boundary-layer thickness) at three different Mach numbers, 2.5, 4.9 and 10.9, with the surface temperatures ranging from quasiadiabatic to highly cooled conditions. The settlement of turbulence statistics into a fully developed equilibrium state of the turbulent boundary layer has been carefully monitored, either based on the satisfaction of the von Kármán integral equation or by comparing runs with different inflow turbulence generation techniques. The generated DNS database is used to characterize the streamwise evolution of multiple important variables in the high-Mach-number, cold-wall regime, including the skin friction, the Reynolds analogy factor, the shape factor, the Reynolds stresses, and the fluctuating wall quantities. The data confirm the validity of many classic and newer compressibility transformations at moderately high Reynolds numbers (up to friction Reynolds number Re τ ≈1200 ) and show that, with proper scaling, the sizes of the near-wall streaks and superstructures are insensitive to the Mach number and wall cooling conditions. The strong wall cooling in the hypersonic cold-wall case is found to cause a significant increase in the size of the near-wall turbulence eddies (relative to the boundary-layer thickness), which leads to a reduced-scale separation between the large and small turbulence scales, and in turn to a lack of an outer peak in the spanwise spectra of the streamwise velocity in the logarithmic region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Flat-histogram extrapolation as a useful tool in the age of big data

In this work, we review recent work by the authors to revisit the concept of extrapolating thermodynamic properties of classical systems using statistical mechanical principles. Specifically, we discuss how the combination of these principles with biased sampling techniques enables the prediction of free energy landscapes and other detailed information, such as structural properties, of the system in question. Remarkably accurate estimates of physical properties across a broad range of conditions have been achieved using this approach, greatly reducing the number of simulations needed to explore a given system's behaviour. While approximate, these extrapolations significantly amplify the amount of reasonably accurate information that can be extracted from simulations enabling a small set of them to feed data-intensive regression algorithms such as neural networks. Thus, this extrapolation methodology represents a useful tool for performing tasks such as high-throughput screening of physical properties, optimising force field parameters, exploring equilibrium phase behaviour, and enabling theory-guided data science for these systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Energetic particle optimization of quasi-axisymmetric stellarator equilibria

Abstract An important goal of stellarator optimization is to achieve good confinement of energetic particles such as, in the case of a reactor, alphas created by deuterium–tritium fusion. In this work, a fixed-boundary stellarator equilibrium was re-optimized for energetic particle confinement via a two-step process: first, by minimizing deviations from quasi-axisymmetry (QA) on a single flux surface near the mid-radius, and secondly by maintaining this improved QA while minimizing the analytical quantity Γ C , which represents the angle between magnetic flux surfaces and contours of J | | , the second adiabatic invariant. This was performed multiple times, resulting in a group of equilibria with significantly reduced energetic particle losses, as evaluated by Monte Carlo simulations of alpha particles in scaled-up versions of the equilibria. This is the first time that energetic particle losses in a QA stellarator have successfully been reduced by optimizing Γ C . The relationship between energetic particle losses and metrics such as QA error ( E q a ) and Γ C in this set of equilibria were examined via statistical methods and a nearly linear relationship between volume-averaged Γ C and prompt particle losses was found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Modeling the Cosmological Lyman-𝛼 Forest at the Field Level

The distribution of absorption lines in the spectra of distant quasars, called the Lyman-𝛼 (Ly-𝛼) forest, is a unique probe of cosmology and the intergalactic medium at high redshifts and small scales. The statistical power of ongoing redshift surveys demands precise theoretical tools to model the Ly-𝛼 forest. We address this challenge by developing an analytic, perturbative forward model to predict the Ly-𝛼 forest at the field level for a given set of cosmological initial conditions. Our model shows a remarkable performance when compared with the Sherwood hydrodynamic simulations: it reproduces the Ly-𝛼 forest flux power spectrum, its cross-correlation with dark matter halos, and the one-point probability distribution function of both fields at the percent level down to scales of a few Mpc. Our work provides crucial tools that bridge analytic modeling on large scales with simulations on small scales, enabling field-level inference from Ly-𝛼 forest data and simulation-based priors for cosmological analyses. Furthermore, this is especially timely for realizing the full scientific potential of the Ly-𝛼 forest measurements by the dark energy spectroscopic instrument.

Cosmological parameters↗

DOE OSTI2020

Deep learning for in situ data compression of large turbulent flow simulations

As the size of turbulent flow simulations continues to grow, in situ data compression is becoming increasingly important for visualization, analysis, and restart checkpointing. For these applications, single-pass compression techniques with low computational and communication overhead are crucial. In this paper we present a deep-learning approach to in situ compression using an autoencoder architecture that is customized for three-dimensional turbulent flows and is well suited for contemporary heterogeneous computing resources. The autoencoder is compared against a recently introduced randomized single-pass singular value decomposition (SVD) for three different canonical turbulent flows: decaying homogeneous isotropic turbulence, a Taylor-Green vortex, and turbulent channel flow. Our proposed fully convolutional autoencoder architecture compresses turbulent flow snapshots by a factor of 64 with a single pass, allows for arbitrarily sized input fields, is cheaper to compute than the randomized single-pass SVD for typical simulation sizes, performs well on unseen flow configurations, and has been made publicly available. The results reported here show that the autoencoder dramatically outperforms a randomized single-pass SVD with similar compression ratio and yields comparable performance to a higher-rank decomposition with an order of magnitude less compression in regard to preserving a number of important statistical quantities such as turbulent kinetic energy, enstrophy, and Reynolds stresses.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

A Generalized Copula-Polynomial Chaos Expansion for Probabilistic Power Flow Considering Nonlinear Correlations of PV Injections

This paper develops a generalized Copula-polynomial chaos expansion (PCE) framework for power system probabilistic power flow that can handle both linear and nonlinear correlations of uncertain power injections, such as wind and PVs. A data-driven Copula statistical model is used to capture the correlations of uncertain power injections. This allows us to resort to the Rosenblatt transformation to transform correlated variables into independent ones while preserving the dependence structure. This paves the way of leveraging the PCE for surrogate modeling and uncertainty quantification of power flow results, i.e., achieving the probabilistic distributions of power flows. Simulations carried out on the IEEE 57-bus system show that the proposed framework can get much more accurate results than other alternatives with different linear and nonlinear power injection correlations.

copula↗

DOE OSTI2023

AEflow (Autoencoder fluid flow compression network) [SWR-22-29]

King, Ryan↗

DOE OSTI2023

The combined effect of temperature and seawater on the compression properties of carbon fiber vinyl ester composites for sandwich structures

For naval applications, composite sandwich structures are of significant interest, and are often manufactured using thin (2 to 4 mm) composite facings made from carbon or glass fiber reinforcement, attached to a thick (25 to 50 mm) section of PVC cellular foam or balsa wood-based core materials using a suitable polymeric resin. In the present study, we focus on the hygrothermal effect on the fiber-dominated compression properties of carbon fiber reinforced vinyl ester resin based polymeric composite (CF/VE), used as “skin” for a polymeric composite sandwich material. Hygrothermal conditioning is achieved by saturating samples in simulated seawater at 40°C. Compression properties are evaluated for coupons extracted along warp and fill undergoing- no conditioning, conditioning till saturation (up to 6 months), and long-term conditioning (2 years). Sea-water saturation yields in up to 12% drop in compression strength with a further 3–4% drop resulting from long-term conditioning. No statistically significant modulus degradation is noticed due to short or long-term hygrothermal exposure. The failure mechanism of the warp extracted coupon, which fails in a splitting failure mode originating due to the delamination between the 0/90 interface, or the fill extracted coupon, which fails due to the instability caused by tow micro-buckling, remains unchanged due to combined exposure (short or long-term) of seawater and temperature. The loss in strength is attributed to the degradation of the fiber-matrix interface, which is validated via conducting single fiber push-in tests with a nominal diameter of 7 micron for conditioned and unconditioned coupons.

Engineering↗

DOE OSTI2024

Dose-efficient automatic differentiation for ptychographic reconstruction

Ptychography, as a powerful lensless imaging method, has become a popular member of the coherent diffractive imaging family over decades of development. The ability to utilize low-dose X-rays and/or fast scans offers a big advantage in a ptychographic measurement (for example, when measuring radiation-sensitive samples), but results in low-photon statistics, making the subsequent phase retrieval challenging. Here, we demonstrate a dose-efficient automatic differentiation framework for ptychographic reconstruction (DAP) at low-photon statistics and low overlap ratio. As no reciprocal space constraint is required in this DAP framework, the framework, based on various forward models, shows superior performance under these conditions. It effectively suppresses potential artifacts in the reconstructed images, especially for the inherent periodic artifact in a raster scan. We validate the effectiveness and robustness of this method using both simulated and measured datasets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Experiments with Undulator Radiation, Emitted by a Single Electron

We study a single electron, circulating in the Fermilab IOTA storage ring and interacting with an undulator through single and multi-photon emissions. The focus of this research is on single-photon and two-photon undulator emissions. We begin by using one Single Photon Avalanche Diode (SPAD) detector to detect the undulator radiation photons and search for possible deviations from the expected Poissonian photon statistics. Then, we go on to use a two-photon interferometer consisting of two SPAD detectors separated by a beam splitter. This allows to test if there is any correlation in the detected photon pairs. In addition, the photocount arrival times can be used to track the longitudinal motion of a single electron and to compare it with simulations. This allowed us to determine several dynamical parameters of the storage ring such as the rf cavity phase jitter and the dependence of the synchrotron motion period on amplitude.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Towards a data-driven model of hadronization using normalizing flows

We introduce a model of hadronization based on invertible neural networks that faithfully reproduces a simplified version of the Lund string model for meson hadronization. Additionally, we introduce a new training method for normalizing flows, termed MAGIC, that improves the agreement between simulated and experimental distributions of high-level (macroscopic) observables by adjusting single-emission (microscopic) dynamics. Our results constitute an important step toward realizing a machine-learning based model of hadronization that utilizes experimental data during training. Finally, we demonstrate how a Bayesian extension to this normalizing-flow architecture can be used to provide analysis of statistical and modeling uncertainties on the generated observable distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Proof that Combining the Forced-collision and DXTRAN Monte Carlo Variance-reduction Techniques is Fair

This report provides mathematical proof that the MCNP DXTRAN (also known in other codes as forced-flight) and forced-collision variance-reduction techniques do not bias the expected value of Monte Carlo simulation estimates when combined. Proof is also provided that the techniques are unbiased when used independently. To prove that the techniques are unbiased, this report derives the first History Score Moment Equation (HSME) for non-multiplying media employing only forced collisions, only DXTRAN, and then both the variance-reduction techniques combined, as defined next. The HSMEs are found by forming the History Score Probability Density Functions (HSPDFs) and then taking the first score moment. Following its derivation, each HSME with a variance-reduction technique employed is reduced to the HSME for an analog simulation. Because the HSME represents the expected contribution to estimators in the simulation, reducing the HSME with variance reduction to the analog HSME shows that the simulation is unbiased despite the variance-reduction technique considered in the HSME. Analysis regarding higher-score moments of each technique is the subject of prior work and is not addressed herein. Throughout this work, the phase space p is defined to be the particle position x, direction-of-flight unit vector $\hat{Ω}$, energy E, and statistical weight w, $$p ≡ (x; \hat{Ω}; E; w).$$ A reduced phase-space excluding the statistical weight of the particle, $$r ≡ (x; \hat{Ω}; E),$$ is also used.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Optimal Sizing of Resilience Solutions for the U.S. Army Reserve

Power or water outages in buildings threaten the ability of the military to support surrounding communities during natural disasters. Outages can last for days, weeks or months. Typical solutions include very expensive batteries and onsite generation that are sized based on historical power needs. The U.S. Army Reserve (USAR) has teamed up with the Pacific Northwest National Laboratory (PNNL) to develop a simulation framework that optimizes future power needs to reduce the cost of resilience solutions. The approach starts with the building. Power needs during an emergency event are simulated at the building end-use level, then loads are reduced through the selection of life cycle cost-effective building-level technology improvements. These new optimized loads are then fed into a microgrid sizing tool that dynamically constructs many different combinations of solar photovoltaic, battery storage, and generator resources to meet the load for hundreds of statistically generated outage scenarios. Six site assessments completed by USAR and PNNL in 2019 have resulted in a 1-14% reduction in overall investment when optimizing the buildings first before determining the generation requirements. This approach helps the Army secure critical missions and provide a 14 day-minimum supply of necessary energy and water in the most efficient manner.

buidlings, efficiency, resilience, FEDS, MCOR↗

DOE OSTI2025

Super-resolution stimulated X-ray Raman spectroscopy

Propagation of intense X-ray pulses through dense media has led to the observation of phenomena such as atomic X-ray lasing, self-induced transparency and stimulated X-ray Raman scattering (SXRS)4. SXRS has been long predicted as a means to launch and probe valence-electron wavepackets and as a building block for nonlinear X-ray spectroscopies. However, experimental observations of SXRS to date have not provided spectroscopic information, and theoretical modelling has largely implemented hard-to-realize phase-coherent attosecond pulses. Here we demonstrate SXRS with spectroscopic precision, that is, detection of valence-excited states in neon with a near Fourier-limited joint energy–time resolution of 0.1 eV–40 fs. We used a new covariance analysis between statistically spiky broadband incident X-ray and scattered X-ray Raman pulses. Using 18,000 single shots, we beat not only the incident (about 8 eV) bandwidth but also the approximately 0.2 eV instrumental energy resolution, thus creating super-resolution conditions, in analogy to super-resolved fluorescence microscopy. Our experimental results, supported by ab initio propagation simulations, reveal the competition between lasing in the ion and stimulated Raman scattering in the neutral. We demonstrate enhanced signal collection efficiency and a broad excitation window, surpassing spontaneous Raman efficiencies by orders of magnitude. This stochastic SXRS approach represents a first step towards tracking elementary events that determine chemical outcomes.

nonlinear x-ray spectroscopy↗