Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation optimization”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mechanics of vitrimer particle compression and fusion under heat press

The compression and fusion of vitrimer particles is a fundamental problem underlying the recycling of vitrimers. Incomplete particle compression can lead to voids in the fused vitrimer, thereby impacting its mechanical property. In this work, we present a two-dimensional finite element model to capture the interplay of multiple complex mechanisms including random particle packing, large deformation, inter-particle contact, thermally activated bulk stress relaxation, and interface healing. Specifically, we focus on understanding two key components of this problem: i) evolution of porosity during compression and fusion of randomly packed vitrimer particles, and ii) effective tensile modulus and strength of fused vitrimer. Using the finite element model, we first show that using smaller particles can reduce the porosity in the fused vitrimer, which only slightly increases the effective modulus, but the increase in tensile strength is more pronounced. In addition, under the same processing conditions, particles with mixed sizes can achieve better densification and higher tensile strength after fusion than uniformly sized particles with the same average radius. Regarding the effects of processing conditions, we discuss how the porosity and effective modulus of the fused vitrimer depend on the processing time, temperature and pressure, and find a general trend consistent with experimental observations, i.e. longer processing time, higher temperature or higher pressure can lead to lower porosity and higher modulus of the fused vitrimer. Finally, the theoretical insights towards the compression and fusion process as well as the simulation framework can be useful in optimizing the powder-based heat press process of vitrimer and its composites.

42 ENGINEERING↗
DOE OSTI2018

Area-preserving scheme for efficiency enhancement in single-pass tapered free electron lasers

Efficiency enhancement of a single-pass short-wavelength high-gain tapered free electron laser (FEL) has recently been intensively studied. The goal is to sustain the growth of radiation power in the post-saturation regime. Among the various schemes, the undulator tapering is considered an effective route to achieve higher power extraction efficiency. The tapering strategy can be of constant or varying resonant phase along the undulator axis. In this paper we propose an efficiency-enhancement scheme based on preservation of the longitudinal phase space area which ensures trapping of resonant particles in the ponderomotive bucket as long as possible along the undulator axis before significant particle depletion occurs. In the meanwhile such a scheme takes advantage of the increase of the radiation field amplitude to precipitate the particle deceleration process at the middle stage of undulator tapering. We analyze such an area-preserving scheme of undulator tapering by formulating the post-saturation FEL interaction in a one-dimensional (1-D) model via introduction of the particle trapping fraction. The output performance is evaluated through numerical iteration and confirmed with 1-D particle tracking simulations. The results show that the optimal power extraction efficiency based on the proposed scheme, together with a prebunched beam, can be greatly improved within relatively short taper length compared with other schemes before radiation diffraction effect becomes significant. Besides, the undesired sideband effects are found to be effectively suppressed. For the proposed area-preserving taper scheme, we also derive an analytical approximate formula for the resonant phase as a function of undulator axis. In conclusion, we expect that the analysis can shed light on the aim to further enhance the power extraction efficiency in single-pass tapered FELs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Gaussian mixture model clustering algorithms for the analysis of high-precision mass measurements

The development of the phase-imaging ion-cyclotron resonance (PI-ICR) technique for use in Penning trap mass spectrometry (PTMS) increased the speed and precision with which PTMS experiments can be carried out. In PI-ICR, data sets of the locations of individual ion hits on a detector are created showing how ions cluster together into spots according to their cyclotron frequency. Ideal data sets would consist of a single, 2D-spherical spot with no other noise, but in practice data sets typically contain multiple spots, non-spherical spots, or significant noise, all of which can make determining the locations of spot centers non-trivial. A method for assigning groups of ions to their respective spots and determining the spot centers is therefore essential for further improving precision and confidence in PI-ICR experiments. Here, we present the class of Gaussian mixture model (GMM) clustering algorithms as an optimal solution. We show that on simulated PI-ICR data, several types of GMM clustering algorithms perform better than other clustering algorithms over a variety of typical scenarios encountered in PI-ICR. The mass spectra of 163Gd, 163 mGd, 162 Tb, and 162 mTb measured using PI-ICR at the Canadian Penning trap mass spectrometer were checked using GMMs, producing results that were in close agreement with the previously published values.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Orientation and Stability of Asymmetric Magnetic Reconnection X Line

Here, the orientation and stability of the reconnection x line in asymmetric geometry is studied using three-dimensional (3-D) particle-in-cell simulations. We initiate reconnection at the center of a large simulation domain to minimize the boundary effect. The resulting x line has sufficient freedom to develop along an optimal orientation, and it remains laminar. Companion 2-D simulations indicate that this x line orientation maximizes the reconnection rate. The divergence of the nongyrotropic pressure tensor breaks the frozen-in condition, consistent with its 2-D counterpart. We then design 3-D simulations with one dimension being short to fix the x line orientation but long enough to allow the growth of the fastest growing oblique tearing modes. This numerical experiment suggests that reconnection tends to radiate secondary oblique tearing modes if it is externally (globally) forced to proceed along an orientation not favored by the local physics. Finally, the development of oblique structure easily leads to turbulence inside small periodic systems.

3-D magnetic reconnection↗
DOE OSTI2017

Mitigating oil spills in the water column

The scale and scope of uncontrolled oil spills can be devastating. Diverse marine environments and fragile ecologies are some of the most susceptible to the many ill effects, while the economic costs can be crippling. A notoriously difficult challenge with no known technological solution is the successful removal of oil dispersed in the water column. Here, we address this problem through cheap and reusable oil sorbents based on the chemical modification of polymer foams. Interfacial chemistry was optimized and subsequently tested in a simulated marine environment at the National Oil Spill Response Research & Renewable Energy Test Facility, Ohmsett. We find favorable performance for surface oil mitigation and, for the first time, demonstrate the advanced sorbent's efficiency and efficacy at pilot scale in extraction of crude oil and refined petroleum products dispersed in the water column. As a result, this is a potentially disruptive technology, opening a new field of environmental science focused on sub-surface pollutant sequestration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Long-Range Monostatic Remote Sensing of Geomaterial Structure Weak Vibrations

We study analytically and numerically signal sensitivity in remote sensing measurements of weak mechanical vibration of structures made of typical construction geomaterials, such as concrete. The analysis includes considerations of electromagnetic beam atmospheric absorption, reflection, scattering, diffraction and losses. Comparison is made between electromagnetic frequencies of 35GHz (Ka-band), 94GHz (W-band) and 260GHz (WR-3 waveguide band), corresponding to atmospheric transparency windows of the electromagnetic spectrum. Numerical simulations indicate that 94GHz frequency is optimal in terms of signal sensitivity and specificity for long-distance (> 1.5km) sensing of weak multi-mode vibrations.

Heifetz, Alexander↗
DOE OSTI2019

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

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

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗
DOE OSTI2020

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

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

Application and performance of an ML-EM algorithm in NEXT

The goal of the NEXT experiment is the observation of neutrinoless double beta decay in (136)Xe using a gaseous xenon TPC with electroluminescent amplification and specialized photodetector arrays for calorimetry and tracking. The NEXT Collaboration is exploring a number of reconstruction algorithms to exploit the full potential of the detector. This paper describes one of them: the Maximum Likelihood Expectation Maximization (ML-EM) method, a generic iterative algorithm to find maximum-likelihood estimates of parameters that has been applied to solve many different types of complex inverse problems. In particular, we discuss a bi-dimensional version of the method in which the photosensor signals integrated over time are used to reconstruct a transverse projection of the event. First results show that, when applied to detector simulation data, the algorithm achieves nearly optimal energy resolution (better than 0.5% FWHM at the Q value of (136)Xe) for events distributed over the full active volume of the TPC.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Imprint of DES superstructures on the cosmic microwave background

Here, small temperature anisotropies in the Cosmic Microwave Background can be sourced by density perturbations via the late-time integrated Sachs-Wolfe effect. Large voids and superclusters are excellent environments to make a localized measurement of this tiny imprint. In some cases excess signals have been reported. We probed these claims with an independent data set, using the first year data of the Dark Energy Survey in a different footprint, and using a different super-structure finding strategy. We identified 52 large voids and 102 superclusters at redshifts $0.2 < z < 0.65$. We used the Jubilee simulation to a priori evaluate the optimal ISW measurement configuration for our compensated top-hat filtering technique, and then performed a stacking measurement of the CMB temperature field based on the DES data. For optimal configurations, we detected a cumulative cold imprint of voids with $\Delta T_{f} \approx -5.0\pm3.7~\mu K$ and a hot imprint of superclusters $\Delta T_{f} \approx 5.1\pm3.2~\mu K$ ; this is $\sim1.2\sigma$ higher than the expected $|\Delta T_{f}| \approx 0.6~\mu K$ imprint of such super-structures in $\Lambda$CDM. If we instead use an a posteriori selected filter size ($R/R_{v}=0.6$), we can find a temperature decrement as large as $\Delta T_{f} \approx -9.8\pm4.7~\mu K$ for voids, which is $\sim2\sigma$ above $\Lambda$CDM expectations and is comparable to previous measurements made using SDSS super-structure data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Higher-order correlations between different moments of two flow amplitudes in Pb-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV

The correlations between different moments of two flow amplitudes, extracted with the recently developed asymmetric cumulants, are measured in Pb-Pb collisions at $\sqrt {s_{NN}}$ = 5.02 TeV recorded by the ALICE detector at the CERN Large Hadron Collider. The magnitudes of the measured observables show a dependence on the different moments as well as on the collision centrality, indicating the presence of nonlinear response in all even moments up to the eighth. Furthermore, the higher-order asymmetric cumulants show different signatures than the symmetric and lower-order asymmetric cumulants. Comparisons with state-of-the-art event generators using two different parametrizations obtained from Bayesian optimization show differences between data and simulations in many of the studied observables, indicating a need for further tuning of the models behind those event generators. These results provide new and independent constraints on the initial conditions and transport properties of the system created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Interaction networks for the identification of boosted H → $b\overline b$ decays

We develop an algorithm based on an interaction network to identify high-transverse-momentum Higgs bosons decaying to bottom quark-antiquark pairs and distinguish them from ordinary jets that reflect the configurations of quarks and gluons at short distances. The algorithm's inputs are features of the reconstructed charged particles in a jet and the secondary vertices associated with them. Describing the jet shower as a combination of particle-to-particle and particle-to-vertex interactions, the model is trained to learn a jet representation on which the classification problem is optimized. The algorithm is trained on simulated samples of realistic LHC collisions, released by the CMS Collaboration on the CERN Open Data Portal. The interaction network achieves a drastic improvement in the identification performance with respect to state-of-the-art algorithms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Inferring Quantum Network Topology Using Local Measurements

Statistical correlations that can be generated across the nodes in a quantum network depend crucially on its topology. However, this topological information might not be known a priori, or it may need to be verified. In this paper, we propose an efficient protocol for distinguishing and inferring the topology of a quantum network. We leverage entropic quantities-namely, the von Neumann entropy and the measured mutual information-as well as measurement covariance to uniquely characterize the topology. We show that the entropic quantities are sufficient to distinguish two networks that prepare GHZ states. Moreover, if qubit measurements are available, both entropic quantities and covariance can be used to infer the network topology without state-preparation assumptions. We show that the protocol can be entirely robust to noise and can be implemented via quantum variational optimization. Numerical experiments on both classical simulators and quantum hardware show that covariance is generally more reliable for accurately and efficiently inferring the topology, whereas entropy-based methods are often better at identifying the absence of entanglement in the low-shot regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Coordinated Optimal Control of Constrained DERs

In this paper, a coordinated optimal control is designed to simultaneously satisfy multiple power system objectives at the distribution network level. With high penetration of renewables, traditional generation and control may not be enough to compensate for intermittent changes in net load. Dynamic coordination of battery storage devices and demand response at aggregated distributed energy resource (DER) sites becomes necessary to maintain power system operation, however, each dynamic power component operates at different time scales. Frequency stabilization can be achieved by designing optimal controls for traditional generation, storage devices and demand response. The proposed optimal controls naturally render frequency separation of these load balancing mechanisms. Given a 24-hour net load forecast, optimal tracking control inputs are designed subject to constraints of selected system elements to ensure optimal performance of the power system. Simulation results are presented and discussed.

Harvey, Roland↗
DOE OSTI2021

Offline Arterial Signal Timing Optimization for Closely Spaced Intersections

The stop-and-go driving behavior at a busy arterial causes a significant amount of fuel waste and emissions that could be reduced. For closely spaced intersections, the traditional bottom-up signal timing approach could fail under high volumes when the queues reach the upstream intersections. A well-coordinated arterial should manage the queue length to prevent the control from failing. Traditionally, coordination among signalized intersections for an arterial means finding the good offsets for each intersection. The control of the intersections along an arterial can also be coordinated by optimizing splits. Splits impact the capacity of an intersection. Similar to the highway, the capacity drop could cause growing queues. This work is based on the Virtual Phase-Link (VPL) model, a street traffic model designed for online traffic model predictive control, to obtain a top-down offline arterial signal timing. The VPL-based model can guarantee capacity consistency in intersections along an arterial. We, therefore, adopted the VPL-based model and developed an offline signal timing optimization approach. The proposed timing derived from the VPL-based offline signal timing optimization showed very good results in simulation. We also collected field experiment data, which demonstrated overall energy reductions and speed improvements on some sections of the study arterial.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Constrained Deep Reinforcement Learning for Energy Sustainable Multi-UAV Based Random Access IoT Networks With NOMA

In this paper, we apply the Non-Orthogonal Multiple Access (NOMA) technique to improve the massive channel access of a wireless IoT network where solar-powered Unmanned Aerial Vehicles (UAVs) relay data from IoT devices to remote servers. Specifically, IoT devices contend for accessing the shared wireless channel using an adaptive p-persistent slotted Aloha protocol; and the solar-powered UAVs adopt Successive Interference Cancellation (SIC) to decode multiple received data from IoT devices to improve access efficiency. To enable an energy-sustainable capacity-optimal network, we study the joint problem of dynamic multi-UAV altitude control and multi-cell wireless channel access management of IoT devices as a stochastic control problem with multiple energy constraints. We first formulate this problem as a Constrained Markov Decision Process (CMDP), and propose an online model-free Constrained Deep Reinforcement Learning (CDRL) algorithm based on Lagrangian primal-dual policy optimization to solve the CMDP. Extensive simulations demonstrate that our proposed algorithm learns a cooperative policy in which the altitude of UAVs and channel access probability of IoT devices are dynamically controlled to attain the maximal long-term network capacity while ensuring energy sustainability of UAVs, outperforming baseline schemes. The proposed CDRL agent can be trained on a small network, yet the learned policy can efficiently manage networks with a massive number of IoT devices and varying initial states, which can amortize the cost of training the CDRL agent.

42 ENGINEERING↗