Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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

Distributionally robust facility location problem under decision-dependent stochastic demand

While the traditional facility location problem considers exogenous demand, in some applications, locations of facilities could affect the willingness of customers to use certain types of services, e.g., carsharing, and therefore they also affect realizations of random demand. Moreover, a decision maker may not know the exact distribution of such endogenous demand and how it is affected by location choices. In this paper, we consider a distributionally robust facility location problem, in which we interpret the moments of stochastic demand as functions of facility-location decisions. We reformulate a two-stage decision-dependent distributionally robust optimization model as a monolithic formulation, and then derive exact mixed-integer linear programming reformulation as well as valid inequalities when the means and variances of demand are piecewise linear functions of location solutions. We conduct extensive computational studies, in which we compare our model with a decision-dependent deterministic model, as well as stochastic programming and distributionally robust models without the decision-dependent assumption. Here, the results show superior performance of our approach with remarkable improvement in profit and quality of service under various settings, in addition to computational speed-ups given by formulation enhancements. These results draw attention to the need of considering the impact of location decisions on customer demand within this strategic-level planning problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Desalination for a circular water economy

Today's water systems are enabled by ample fresh water sources, low-cost centralized treatment, and facile wastewater disposal. Climatic change, aging infrastructure, and source water contamination have exposed the vulnerabilities of this linear water paradigm. While seawater desalination enables coastal communities to augment their supply, more broadly securing water systems for municipal, industrial, and agricultural water users will require distributed desalination and fit-for-purpose reuse of nontraditional water sources. Our linear water economy must evolve into a resilient circular water economy, where water is continuously reused and “contaminants” become the feedstocks for other economically valuable processes. Additionally, technology innovation is needed to deliver autonomous, precise, resilient, intensified, modular, and electrified desalination systems that reduce the cost, improve the performance, and enhance the resilience of nontraditional water reuse systems. Meanwhile, strong federal leadership and coordination is needed to accelerate desalination research, promote information gathering efforts to direct technology development, and create an expanded role for non-profit organizations in knowledge dissemination.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Calcined Polyethyleneimine Coated Optical Fibers for Distributed pH Monitoring at High Pressures and Temperatures

In the oil and gas, CO2 sequestration, H2 subsurface storage, and geothermal energy sectors, subsurface pH measurements are critical for monitoring the geochemical conditions and estimating potential corrosion rates of wellbore systems. Real-time pH measurements in these conditions are vital for detecting and predicting corrosion deterioration of wellbore components that may jeopardize the safety and continued operation of wellbore systems. Building off of earlier TiO2 based pH sensors and the known pH sensitivity of amine-based polymers, a coating based on the secondary amine polymer polyethyleneimine (PEI) was developed. The sensor coating was treated with a high temperature (500 °C) calcination procedure in air to convert it into a more stable oxidized coating capable of withstanding hot aqueous solutions without dissolving while retaining linear pH sensitivity from pH values between 2 and 11. The sensor performance was measured using optical transmission measurements in solutions of various pHs and using optical backscatter reflectometry for distributed pH sensing demonstration in wellbore-relevant pressures (up to 1000 psi) and temperatures (80 °C). A calibration curve with strong differentiation between acidic and alkaline pH was developed for both transmission-based and distributed pH measurements, using fixed wavelength transmission and integrated linear amplitude of backscattered light for distributed measurements.

Shumski, Alexander↗
DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Pressure safety approach for PIP-II cryogenic distribution system and cryomodules

The Proton Improvement Plan-II (PIP-II) is a superconducting linear accelerator being built at Fermilab that will provide 800 MeV proton beam for neutrino production. The linac consists of a total of twenty-three (23) cryomodules of five (5) different types. Cooling is required at 2K, 5K and 40K. The Cryogenic Distribution System (CDS) consists of a Distribution Valve Box, ~285 m of cryogenic transfer line, modular Bayonet Cans to interface with cryomodules, and a Turnaround Can. The cryogenic system must provide protection from over-pressure by sizing pressure relief devices for all volumes and process line circuits. The cryomodule cavity circuits have dual pressure ratings, 4.1 bara when cold and 2.05 bara when warm (T>80K). Worst case relieving cases will be identified. The methods for determining heat flux will be presented. For the relieving occurring in the linac tunnel, flow must vent to outside to avoid an oxygen deficiency hazard. Also, we will present vacuum vesse l relief sizing to protect the cryogenic distribution system vacuum shells from over pressure during an internal line rupture. The project is funded by US DOE Offices of Science, High Energy Physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Pressure safety approach for PIP-II cryogenic distribution system and cryomodules

The Proton Improvement Plan-II (PIP-II) is a superconducting linear accelerator being built at Fermilab that will provide 800 MeV proton beam for neutrino production. The linac consists of a total of twenty-three (23) cryomodules of five (5) different types. Cooling is required at 2K, 5K and 40K. The Cryogenic Distribution System (CDS) consists of a Distribution Valve Box, ~285 m of cryogenic transfer line, modular Bayonet Cans to interface with cryomodules, and a Turnaround Can. The cryogenic system must provide protection from over-pressure by sizing pressure relief devices for all volumes and process line circuits. The cryomodule cavity circuits have dual pressure ratings, 4.1 bara when cold and 2.05 bara when warm (T>80K). Worst case relieving cases will be identified. The methods for determining heat flux will be presented. For the relieving occurring in the linac tunnel, flow must vent to outside to avoid an oxygen deficiency hazard. Also, we will present vacuum vessel relief sizing to protect the cryogenic distribution system vacuum shells from over pressure during an internal line rupture. The project is funded by US DOE Offices of Science, High Energy Physics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Re-Analysis of DARHT Axis 2 S4 Magnet Sweep Measurements

The Dual-Axis Radiographic Hydrodynamic Test (DARHT) facility provides flash radiography capabilities using two electron Linear Induction Accelerators (LIA’s). Understanding fundamental properties of the electron beams is essential to optimizing spot size and dose for radiography. This work describes distribution measurements of a single-kicked electron beam (16.5MeV, 1.7kA, 60ns) at a location downstream at imaging station C on DARHT Axis 2. The beam at station C is focused by the S4 solenoid located upstream. The beam distribution is measured by imaging optical transition radiation (OTR) from electrons striking a 51-µm thick titanium foil. This data and analysis were originally published by Ekdahl. The results here are a re-analysis of the asymmetric halo emittance contribution as compared to the emittance resulting from the FWHM of the distribution. This is accomplished using xtr2 fits to the data which include non-linear magnetic field effects in S4. Also examined are the calculated S4 spherical aberration contributions to the emittance using xtr.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Machine Learning Methods for Connection RTT and Loss Rate Estimation Using MPI Measurements Under Random Losses

Scientific computations are expected to be increasingly distributed across wide-area networks, and Message Passing Interface (MPI) has been shown to scale to support their communications over long distances. Application-level measurements of MPI operations reflect the connection Round-Trip Time (RTT) and loss rate, and machine learning methods have been previously developed to estimate them under deterministic periodic losses. In this paper, we consider more complex, random losses with unform, Poisson and Gaussian distributions. We study five disparate machine leaning methods, with linear and non-linear, and smooth and non-smooth properties, to estimate RTT and loss rate over 10Gbps connections with 0-366ms RTT. The diversity and complexity of these estimators combined with the randomness of losses and TCP’s non-linear response together rule out the selection of a single best among them; instead, we fuse them to retain their design diversity. Overall, the results show that accurate estimates can be generated at low loss rates but become inaccurate at loss rates 10% and higher, thereby illustrating both their strengths and limitations.

Generalization Bounds↗
DOE OSTI2022

Exploring Gradient-Based Multi-directional Controls in GANs

Generative Adversarial Networks (GANs) have been widely applied in modeling diverse image distributions. However, despite its impressive applications, the structure of the latent space in GANs largely remains as a black-box, leaving its controllable generation an open problem, especially when spurious correlations between different semantic attributes exist in the image distributions. To address this problem, previous methods typically learn linear directions or individual channels that control semantic attributes in the image space. However, they often suffer from imperfect disentanglement, or are unable to obtain multi-directional controls. Furthermore, in this work, in light of the above challenges, we propose a novel approach that discovers nonlinear controls, which enables multi-directional manipulation as well as effective disentanglement, based on gradient information in the learned GAN latent space. More specifically, we first learn interpolation directions by following the gradients from classification networks trained separately on the attributes, and then navigate the latent space by exclusively controlling channels activated for the target attribute in the learned directions. Empirically, with small training data, our approach is able to gain fine-grained controls over a diverse set of bi-directional and multi-directional attributes, and we showcase its ability to achieve disentanglement significantly better than state-of-the-art methods both qualitatively and quantitatively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗
DOE OSTI2024

Intermittency, bursty turbulence, and ion and electron phase-space holes formation in collisionless current-carrying plasmas

In the previous studies of nonlinear saturation of the Buneman instability caused by high electron drift velocity relative to ions, the phase-space holes and the plateau on the electron velocity distribution function were identified as features of the saturation stage of instability [notably in the paper by Omura et al., J. Geophys. Res. 108, 1197 (2003)]. We have performed a much longer simulation of the Buneman instability and observed a secondary instability. This secondary instability generates fast electron-acoustic waves. By analyzing the phase-space plot of ions and electrons, we show that the fast electron heating and the formation of the plateau of electron velocity distribution function are not due to the quasi-linear diffusion but due to the nonlinear interaction of ion- and electron-acoustic solitary waves (phase-space holes) by exchange of trapped electrons in each wave. We also report the details on the intermittent and bursty nature of turbulence driven by this instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Light-front puzzles

Abstract Light-front formulations of quantum field theories have many advantages for computing electroweak matrix elements of strongly interacting systems and other quantities that are used to study hadronic structure. The theory can be formulated in Hamiltonian form so non-perturbative calculations of the strongly interacting initial and final states are in principle reduced to linear algebra. These states are needed for calculating parton distribution functions and other types of distribution amplitudes that are used to understand the structure of hadrons. Light-front boosts are kinematic transformations so the strongly interacting states can be computed in any frame. This is useful for computing current matrix elements involving electroweak probes where the initial and final hadronic states are in different frames related by the momentum transferred by the probe. Finally in many calculations the vacuum is trivial so the calculations can be formulated in Fock space. The advantages of light front-field theory would not be interesting if the light-front formulation was not equivalent to the covariant or canonical formulations of quantum field theory. Many of the distinguishing properties of light-front quantum field theory are difficult to reconcile with canonical or covariant formulations of quantum field theory. This paper discusses the resolution of some of the apparent inconsistencies in canonical, covariant and light-front formulations of quantum field theory. The puzzles that will be discussed are (1) the problem of inequivalent representations (2) the problem of the trivial vacuum (3) the problem of ill-posed initial value problems (4) the problem of rotational covariance (5) the problem of zero modes and (6) the problem of spontaneously broken symmetries.

Physics↗
DOE OSTI2025

Minimizing phase-space energies

A primary technical challenge for harnessing fusion energy is to control and extract energy from a nonthermal distribution of charged particles. The fact that phase space evolves by symplectomorphisms fundamentally limits how a distribution may be manipulated. While the constraint of phase-space volume preservation is well understood, other constraints remain to be fully appreciated. To better understand these constraints, we study the problem of extracting energy from a distribution of particles using area-preserving and symplectic linear maps. When a quadratic potential is imposed, we find that the maximal extractable energy can be computed as trace minimization problems. We solve these problems and show that the extractable energy under linear symplectomorphisms may be much smaller than the extractable energy under special linear maps. As a result, the method introduced in the present study enables an energy-based proof of the linear Gromov nonsqueezing theorem.

First-principles calculations in plasma physics↗
DOE OSTI2023

Cross Entropy Benchmark for Measurement-Induced Phase Transitions

Here, we investigate prospects of employing the linear cross entropy to experimentally access measurement-induced phase transitions without requiring any postselection of quantum trajectories. For two random circuits that are identical in the bulk but with different initial states, the linear cross entropy χ between the bulk measurement outcome distributions in the two circuits acts as an order parameter, and can be used to distinguish the volume law from area law phases. In the volume law phase (and in the thermodynamic limit) the bulk measurements cannot distinguish between the two different initial states, and χ = 1. In the area law phase χ < 1. For circuits with Clifford gates, we provide numerical evidence that χ can be sampled to accuracy ϵ from O⁢(1/ϵ 2 ) trajectories, by running the first circuit on a quantum simulator without postselection, aided by a classical simulation of the second. We also find that for weak depolarizing noise the signature of the measurement-induced phase transitions is still present for intermediate system sizes. In our protocol we have the freedom of choosing initial states such that the “classical” side can be simulated efficiently, while simulating the “quantum” side is still classically hard.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Bunch Shaping in Electron Linear Accelerators

Modern electron linear accelerators are often designed to produce smooth bunch distributions characterized by their macroscopic ensemble-average moments. However, an increasing number of accelerator applications call for finer control over the beam distribution, such as by requiring specific shapes for its projection along one coordinate. Ultimately, the control of the beam distribution at the single-particle level could enable new opportunities in accelerator science. This review discusses the recent progress toward controlling electron beam distributions on the “mesoscopic” scale with an emphasis on shaping the beam or introducing complex correlations required for some applications. This review emphasizes experimental and theoretical developments of electron-bunch shaping methods based on bounded external electromagnetic fields or via interactions with the self-generated velocity and radiation fields.

Ha, G.↗
DOE OSTI2024

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Machine Learning Methods for Connection RTT and Loss Rate Estimation Using MPI Measurements Under Random Losses

Scientific computations are expected to be increasingly distributed across wide-area networks, and Message Passing Interface (MPI) has been shown to scale to support their communications over long distances. Application-level measurements of MPI operations reflect the connection Round-Trip Time (RTT) and loss rate, and machine learningmethods have been previously developed to estimate them under deterministic periodic losses. In this paper, we consider more complex, random losses with uniform, Poisson and Gaussian distributions. We study five disparate machine leaning methods, with linear and non-linear, and smooth and non-smooth properties, to estimate RTT and loss rate over 10 Gbps connections with 0–366 ms RTT. The diversity and complexity of these estimators combined with the randomness of losses and TCP’s non-linear response together rule out the selection of a single best among them; instead, we fuse them to retain their design diversity. Overall, the results show that accurate estimates can be generated at low loss ratesbut become inaccurate at loss rates 10% and higher, thereby illustrating both their strengths and limitations.

Rao, Nageswara S.↗