Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Atacama Cosmology Telescope: Component-separated maps of CMB temperature and the thermal Sunyaev-Zel’dovich effect

Optimal analyses of many signals in the cosmic microwave background (CMB) require map-level extraction of individual components in the microwave sky, rather than measurements at the power spectrum level alone. To date, nearly all map-level component separation in CMB analyses has been performed exclusively using satellite data. In this paper, we implement a component separation method based on the internal linear combination (ILC) approach which we have designed to optimally account for the anisotropic noise (in the 2D Fourier domain) often found in ground-based CMB experiments. Using this method, we combine multi-frequency data from the Planck satellite and the Atacama Cosmology Telescope Polarimeter (ACTPol) to construct the first wide-area, arcminute-resolution component-separated maps (≈ 2100 sq. deg.) of the CMB temperature anisotropy and the thermal Sunyaev-Zel'dovich (tSZ) effect sourced by the inverse-Compton scattering of CMB photons off hot, ionized gas. Our ILC pipeline allows for explicit deprojection of various contaminating signals, including a modified blackbody approximation of the cosmic infrared background (CIB) spectral energy distribution. The cleaned CMB maps will be a useful resource for CMB lensing reconstruction, kinematic SZ cross-correlations, and primordial non-Gaussianity studies. Here, the tSZ maps will be used to study the pressure profiles of galaxies, groups, and clusters through cross-correlations with halo catalogs, with dust contamination controlled via CIB deprojection. The data products described in this paper are available on LAMBDA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Non-linear Model Predictive Control of Wave Energy Converters with Realistic Power Take-off Configurations and Loss Model

Model Predictive Control (MPC) has been recognized as a well-developed control strategy for optimizing the performance of wave energy converters (WECs). The standard problem involves maximizing the absorbed power from oncoming waves while respecting motion constraints of the device and force constraints of the power take-off (PTO). The generated electrical power is then calculated by assuming an ideal PTO with no power conversion losses or limitations to power flow. In this article, we remove the assumption of an ideal PTO and present four different control options that reflect unique combinations of PTO control and power flow constraints. We also recast the MPC problem to maximize the generated power instead of the absorbed power by accounting for the losses in the power train. Using a WEC with a hydraulic PTO, we develop a loss model suitable for controls optimization. Lastly, we illustrate the process of setting up the non-linear MPC problem and compare the performance of competing control strategies using a common benchmark based on annual energy captured at a reference site.

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

Nature of size effects in compact models of field effect transistors

Investigations have shown that in the local approximation (for sizes L < 100 μm), AlGaN/GaN high electron mobility transistor (HEMT) structures satisfy to all properties of chaotic systems and can be described in the language of fractal geometry of fractional dimensions. For such objects, values of their electrophysical characteristics depend on the linear sizes of the examined regions, which explain the presence of the so-called size effects—dependences of the electrophysical and instrumental characteristics on the linear sizes of the active elements of semiconductor devices. In the present work, a relationship has been established for the linear model parameters of the equivalent circuit elements of internal transistors with fractal geometry of the heteroepitaxial structure manifested through a dependence of its relative electrophysical characteristics on the linear sizes of the examined surface areas. For the HEMTs, this implies dependences of their relative static (A/mm, mA/V/mm, Ω/mm, etc.) and microwave characteristics (W/mm) on the width d of the sink-source channel and on the number of sections n that leads to a nonlinear dependence of the retrieved parameter values of equivalent circuit elements of linear internal transistor models on n and d. Thus, it has been demonstrated that the size effects in semiconductors determined by the fractal geometry must be taken into account when investigating the properties of semiconductor objects on the levels less than the local approximation limit and designing and manufacturing field effect transistors. In general, the suggested approach allows a complex of problems to be solved on designing, optimizing, and retrieving the parameters of equivalent circuits of linear and nonlinear models of not only field effect transistors but also any arbitrary semiconductor devices with nonlinear instrumental characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Integrate FARM with PID controllers: IES Simulation Ecosystem Control System Development

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM was designed to support the HERON software module in the solution of the optimal dispatch problem for IES units. As the result of HERON-FARM dispatch simulation, the set-point trajectories are optimized to meet constraints on both the production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and the process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.) at a coarse time resolution (every 10 or 100 seconds) over long time horizons (several days or weeks). In case the operational constraints need to be met at finer time resolution, the computational burden of HERON-FARM would linearly increase with the sampling rate, and sub-optimal solutions might be obtained. System responses characterized by overshoots and damped oscillations temporarily violating the imposed constraints might occur during abrupt power transients. In this report, a hierarchical control system architecture for the operation of the Dynamic Energy Technology and Integration Laboratory (DETAIL) facility constructed at INL was proposed. First, the preliminary studies on the proposed control strategy for operating the facility and the designed PI controllers were reviewed. In particular, the current approach for generating the set-point trajectories was studied, and its limits were identified. To this aim, the inclusion of a Supervisory Control layer embedding a modified version of the FARM algorithm for preserving the system safe operation over both long and real-time horizons was proposed. In this way, FARM would be applied twice, i.e., the original version (“FARM Validator”) aiding the solution of the power dispatch problem, and the modified version (“FARM Supervisory” coordinating the PI controllers to address the real-time control tasks. Despite the kernel of the two modules is the same algorithm, their roles, tasks, and capabilities are quite different. A detailed description of the role of FARM at addressing low-level control tasks is provided, along with tentative operational procedures for training the models embedded into the algorithm by using the collected experimental data.

42 ENGINEERING↗
DOE OSTI2024

Revealing Decision Conservativeness Through Inverse Distributionally Robust Optimization

This paper introduces Inverse Distributionally Robust Optimization (I-DRO) as a method to infer the conservativeness level of a decision-maker, represented by the size of a Wasserstein metric-based ambiguity set, from the optimal decisions made using Forward Distributionally Robust Optimization (F-DRO). By leveraging the Karush-Kuhn-Tucker (KKT) conditions of the convex F-DRO model, we formulate I-DRO as a bi-linear program, which can be solved using off-the-shelf optimization solvers. Additionally, this formulation exhibits several advantageous properties. We demonstrate that I-DRO not only guarantees the existence and uniqueness of an optimal solution but also establishes the necessary and sufficient conditions for this optimal solution to accurately match the actual conservativeness level in F-DRO. Furthermore, we identify three extreme scenarios that may impact I-DRO effectiveness. Our case study applies F-DRO for power system scheduling under uncertainty and employs I-DRO to recover the conservativeness level of system operators. Numerical experiments based on an IEEE 5-bus system and a realistic NYISO 11-zone system demonstrate I-DRO performance in both normal and extreme scenarios. An extended version of this paper with additional analyses is available at li2024revealing.

distributionally robust optimization↗
DOE OSTI2015

Quantum dot spontaneous emission control in a ridge waveguide

We investigate the spontaneous emission (SE) of self-assembled InAs quantum dots (QDs) embedded in GaAs ridge waveguides that lay on a low index substrate. In thin enough waveguides, the coupling to the fundamental guided mode is vanishingly small. A pronounced anisotropy in the coupling to non-guided modes is then directly evidenced by normal-incidence photoluminescence polarization measurements. In this regime, a measurement of the QD decay rate reveals a SE inhibition by a factor up to 4. In larger wires, which ensure an optimal transverse confinement of the fundamental guided mode, the decay rate approaches the bulk value. Building on the good agreement with theoretical predictions, we infer from calculations the fraction β of SE coupled to the fundamental guided mode for some important QD excitonic complexes. For a charged exciton (isotropic in plane optical dipole), β reaches 0.61 at maximum for an on-axis QD. In the case of a purely transverse linear optical dipole, β increases up to 0.91. This optimal configuration is achievable through the selective excitation of one of the bright neutral excitons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

SimCCS: An open-source tool for optimizing CO 2 capture, transport, and storage infrastructure

Commercial-scale carbon capture and storage (CCS) technology will involve deploying infrastructure on a massive and costly scale. This effort will require careful and comprehensive planning to ensure that capture locations, storage sites, and the dedicated CO 2 distribution pipelines are selected in a robust and cost-effective manner. Presented in 2009, SimCCS is an optimization model for integrated system design that enables researchers, stakeholders, and policy makers to design CCS infrastructure networks. SimCCS 2.0 is a complete, ground-up redesign that is now a portable software package, useable and shareable by the CCS research, industrial, policy, and public communities. SimCCS 2.0 integrates multiple new capabilities including a refined optimization model, novel candidate network generation techniques, and optional integration with high-performance computing platforms. Accessing user-provided CO 2 source, sink, and transportation data, SimCCS 2.0 creates candidate transportation routes and formalizes an optimization problem that determines the most cost-effective CCS system design. This optimization problem is further solved either through a high-performance computing interface, or through third-party software on a local desktop computing platform. Finally, SimCCS 2.0 employs an open-access geographic information system framework to enable analysis and visualization capabilities. SimCCS 2.0 is written in Java and is publicly available via GitHub to encourage collaboration, modification, and community development.

42 ENGINEERING↗
DOE OSTI2018

The optimal use of silicon pixel charge information for particle identification

Particle identification using the energy loss in silicon detectors is a powerful technique for probing the Standard Model (SM) as well as searching for new particles beyond the SM. Traditionally, such techniques use the truncated mean of the energy loss on multiple layers, in order to mitigate heavy tails in the charge fluctuation distribution. We show that the optimal scheme using the charge in multiple layers significantly outperforms the truncated mean. Furthermore, truncation itself does not significantly degrade performance and the optimal classifier is well-approximated by a linear combination of the truncated mean and truncated variance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Network Modeling of the U.S. Petrochemical Industry under Raw Material and Hurricane Harvey Disruptions

A geographically resolved network model of the U.S. chemical industry in 2017 is designed, and optimal chemical flows between units are calculated using linear programming. A baseline solution and three disruption scenarios (primary raw material disruptions, reported Hurricane Harvey ethylene cracker disruptions, and assumed capacity disruptions based on the Hurricane Harvey geographic storm track) are reported on to determine how the structure of the industry is modified to adapt to widespread and geographically specific disruptions. The calculated impacts of the assumed Hurricane Harvey disruption include 170 chemical units in 26 states that change production level as a result of supply chain disruptions during the storm. The systemic impact for the assumed Hurricane Harvey disruption is 19.3 million tonnes of gross chemical production. The day with the largest impact on gross chemical production shows a reduction from baseline operations of 1.3 million tonnes (42% of baseline). This model can be utilized for analysis of future disruption scenarios and to test resilience strategies, including impacts of new manufacturing configurations or technologies.

02 PETROLEUM↗
DOE OSTI2026

Direct and rapid uranium isotopic analysis of environmental sample swipes via microextraction – LS-APGD/Orbitrap mass spectrometry

Accurate and precise isotopic analysis of actinides collected on environmental sample swipes from within nuclear facilities is an important safeguarding measure for detecting undeclared materials and activities. Traditional isotope ratio (IR) analysis of actinides using a bulk digestion approach can be laborious and time-consuming. Recently, direct analysis of environmental swipe samples using a microextraction approach has been explored as an alternative to conventional bulk digestion methods. The present study further evaluates this approach for accurate and precise isotopic analysis of uranium (U) in cotton swipe samples using the liquid sampling-atmospheric pressure glow discharge (LS-APGD)/Orbitrap-FTMS Booster detection system. The instrumental parameters were optimized, and quantitative capabilities were demonstrated through excellent linearity (R2 = 0.99) across a deposited mass range of 1 to 500 ng, with a limit of detection of 90 pg (in a 2 × 4 mm region) for total U; well below typical U concentrations (ng–mg range) in environmental sample swipes. The 235U/238U ratio showed excellent accuracy (−0.3% relative difference from the certificate value), with ∼20× improvement over previous microextraction-Orbitrap methods. Precision (∼1% relative standard deviation) was also greatly enhanced by ∼5×. This study also demonstrates that figures of merit achieved with neat U solution were not highly degraded when various external interfering elements/matrices were introduced. Finally, the developed method was successfully applied to analyze practical swipe samples representative of laboratory and outdoor environments. The presented figures of merit and applicability to practical samples further validate the capability of the platform for rapid analysis (<7 minutes) of environmental swipe samples.

Shrestha, Suraj [Clemson University, SC]↗
DOE OSTI2012

RF optimization and analysis of the 805-MHz cavity for the MuCool program using ACE3P

An 805 MHz pillbox cavity tested at Fermilab's MTA facility showed significant degradation in gradient when operated in a several Tesla solenoidal magnetic field. We have used the advanced ACE3P simulation codes developed at SLAC to study the cavity dark current and multipacting characteristics to gain more insight into the gradient limitations. We also checked whether there is an optimal cavity length that minimizes the dark current impact energy. Finally, we have improved on the cavity design, significantly lowering the fields outside the beam area. These and other results are presented in this paper.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Intrinsic delay of permeable base transistor

Permeable base transistors (PBTs) fabricated by vacuum deposition or solution process have the advantages of easy fabrication and low power operation and are a promising device structure for flexible electronics. Intrinsic delay of PBT, which characterizes the speed of the transistor, is investigated by solving the three-dimensional Poisson equation and drift-diffusion equation self-consistently using finite element method. Decreasing the emitter thickness lowers the intrinsic delay by improving on-current, and a thinner base is also preferred for low intrinsic delay because of fewer carriers in the base region at off-state. The intrinsic delay exponentially decreases as the emitter contact Schottky barrier height decreases, and it linearly depends on the carrier mobility. With an optimized emitter contact barrier height and device geometry, a sub-nano-second intrinsic delay can be achieved with a carrier mobility of ∼10 cm{sup 2}/V/s obtainable in solution processed indium gallium zinc oxide, which indicates the potential of solution processed PBTs for GHz operations.

36 MATERIALS SCIENCE↗
DOE OSTI2014

High-efficiency acceleration in the laser wakefield by a linearly increasing plasma density

The acceleration length and the peak energy of the electron beam are limited by the dephasing effect in the laser wakefield acceleration with uniform plasma density. Based on 2D-3V particle in cell simulations, the effects of a linearly increasing plasma density on the electron acceleration are investigated broadly. Comparing with the uniform plasma density, because of the prolongation of the acceleration length and the gradually increasing accelerating field due to the increasing plasma density, the electron beam energy is twice higher in moderate nonlinear wakefield regime. Because of the lower plasma density, the linearly increasing plasma density can also avoid the dark current caused by additional injection. At the optimal acceleration length, the electron energy can be increased from 350 MeV (uniform) to 760 MeV (linearly increasing) with the energy spread of 1.8%, the beam duration is 5 fs and the beam waist is 1.25 μm. This linearly increasing plasma density distribution can be achieved by a capillary with special gas-filled structure, and is much more suitable for experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The geostatistical approach for structural and stratigraphic framework analysis of offshore NW Bonaparte Basin, Australia

Geostatistics or statistical approach is based on the studies of temporal and spatial trend, which depend upon spatial relationships to model known information of variable(s) at unsampled locations. The statistical technique known as kriging was used for petrophycial and facies analysis, which help to assume spatial relationship to model the geological continuity between the known data and the unknown to produce a single best guess of the unknown. Kriging is also known as optimal interpolation technique, which facilitate to generate best linear unbiased estimation of each horizon. The idea is to construct a numerical model of the lithofacies and rock properties that honor available data and further integrate with interpreting seismic sections, techtonostratigraphy chart with sea level curve (short term) and regional tectonics of the study area to find the structural and stratigraphic growth history of the NW Bonaparte Basin. By using kriging technique the models were built which help to estimate different parameters like horizons, facies, and porosities in the study area. The variograms were used to determine for identification of spatial relationship between data which help to find the depositional history of the North West (NW) Bonaparte Basin.

58 GEOSCIENCES↗
DOE OSTI2017

Magnetic field production via the Weibel instability in interpenetrating plasma flows

Here, many astrophysical systems are effectively “collisionless,” that is, the mean free path for collisions between particles is much longer than the size of the system. The absence of particle collisions does not preclude shock formation, however, as shocks can be the result of plasma instabilities that generate and amplify electromagnetic fields. The magnetic fields required for shock formation may either be initially present, for example, in supernova remnants or young galaxies, or they may be self-generated in systems such as gamma-ray bursts (GRBs). In the case of GRB outflows, the Weibel instability is a candidate mechanism for the generation of sufficiently strong magnetic fields to produce shocks. In experiments on the OMEGA Laser, we have demonstrated a quasi-collisionless system that is optimized for the study of the non-linear phase of Weibel instability growth. Using a proton probe to directly image electromagnetic fields, we measure Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows. The collisionality of the system is determined from coherent Thomson scattering measurements, and the data are compared to similar measurements of a fully collisionless system. The strong, persistent Weibel growth observed here serves as a diagnostic for exploring large-scale magnetic field amplification and the microphysics present in the collisional–collisionless transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

HyKKT: a hybrid direct-iterative method for solving KKT linear systems

Here, we propose a solution strategy for the large indefinite linear systems arising in interior methods for nonlinear optimization. The method is suitable for implementation on hardware accelerators such as graphical processing units (GPUs). The current gold standard for sparse indefinite systems is the LBLT factorization where L is a lower triangular matrix and B is 1×1 or 2×2 block diagonal. However, this requires pivoting, which substantially increases communication cost and degrades performance on GPUs. Our approach solves a large indefinite system by solving multiple smaller positive definite systems, using an iterative solver on the Schur complement and an inner direct solve (via Cholesky factorization) within each iteration. Cholesky is stable without pivoting, thereby reducing communication and allowing reuse of the symbolic factorization. We demonstrate the practicality of our approach on large optimal power flow problems and show that it can efficiently utilize GPUs and outperform LBL T factorization of the full system.

97 MATHEMATICS AND COMPUTING↗