Search NASA⌕ Search

SEARCH · Search NASA

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 541 records · Page 30

Classical symmetries and the Quantum Approximate Optimization Algorithm

Here, we study the relationship between the Quantum Approximate Optimization Algorithm (QAOA) and the underlying symmetries of the objective function to be optimized. Our approach formalizes the connection between quantum symmetry properties of the QAOA dynamics and the group of classical symmetries of the objective function. The connection is general and includes but is not limited to problems defined on graphs. We show a series of results exploring the connection and highlight examples of hard problem classes where a nontrivial symmetry subgroup can be obtained efficiently. In particular, we show how classical objective function symmetries lead to invariant measurement outcome probabilities across states connected by such symmetries, independent of the choice of algorithm parameters or number of layers. To illustrate the power of the developed connection, we apply machine learning techniques toward predicting QAOA performance based on symmetry considerations. We provide numerical evidence that a small set of graph symmetry properties suffices to predict the minimum QAOA depth required to achieve a target approximation ratio on the MaxCut problem, in a practically important setting where QAOA parameter schedules are constrained to be linear and hence easier to optimize.

97 MATHEMATICS AND COMPUTING↗

Reduced-order modeling on a near-term quantum computer

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗

SU-E-T-153: Burst-Mode Modulated Arc Therapy with Flattening-Filter-Free Beams Versus Flattening-Filtered Beams

Purpose: To compare the dosimetry and delivery of burst-mode modulated arc radiotherapy using flattening-filter-free (FFF) and flattening-filtered (FF) beams. Methods: Burst-mode modulated arc therapy (mARC, Siemens) plans were generated for six prostate cases with FFF and FF beam models, using the Elekta Monaco v. 5.00 planning system. One 360-degree arc was used for five cases, and for one case two 360-degree coplanar arcs were used. The maximum number of optimization points (OPs) per arc was set to 91, and OPs with less than 4 MU were disregarded. All plans were delivered on the Siemens Artiste linear accelerator with 6MV FF (300 MU/min) and comparable-energy FFF (2000 MU/min, labeled as 7UF) beams. Results: For all cases studied, the plans with FFF beams exhibited DVHs for the PTV, rectum, and bladder that were nearly identical to those for the plans with FF beams. The FFF plan yielded reduced dose to the right femoral head for 5 cases, and lower mean dose to the left femoral head for 4 cases. For all but the two-arc case, the FFF and FF plans resulted in an identical number of segments. The total number of MUs was slightly lower for the FF plans for five cases. The total delivery time per fraction was substantially lower for the FFF plans, ranging from 25 to 50 percent among all cases, as compared to the FF plans. Conclusion: For mARC plans, FFF and FF beams provided comparable PTV coverage and rectum and bladder sparing. For the femoral heads, the mean dose was slightly lower in most cases when using the FFF beam. Although the flat beam plans typically required slightly fewer MUs, FFF beams required substantially less time to deliver a plan of similar quality. This work was supported by Siemens Medical Solutions and the MCW Cancer Center Fotsch Foundation.

60 APPLIED LIFE SCIENCES↗

SU-E-T-368: Evaluating Dosimetric Outcome of Modulated Photon Radiotherapy (XMRT) Optimization for Head and Neck Patients

Purpose: The dosimetric outcome of optimized treatment plans obtained by modulating the photon beamlet energy and fluence on a small cohort of four Head and Neck (H and N) patients was investigated. This novel optimization technique is denoted XMRT for modulated photon radiotherapy. The dosimetric plans from XMRT for H and N treatment were compared to conventional, 6 MV intensity modulated radiotherapy (IMRT) optimization plans. Methods: An arrangement of two non-coplanar and five coplanar beams was used for all four H and N patients. Both XMRT and IMRT were subject to the same optimization algorithm, with XMRT optimization allowing both 6 and 18 MV beamlets while IMRT was restricted to 6 MV only. The optimization algorithm was based on a linear programming approach with partial-volume constraints implemented via the conditional value-at-risk method. H and N constraints were based off of those mentioned in the Radiation Therapy Oncology Group 1016 protocol. XMRT and IMRT solutions were assessed using metrics suggested by International Commission on Radiation Units and Measurements report 83. The Gurobi solver was used in conjunction with the CVX package to solve each optimization problem. Dose calculations and analysis were done in CERR using Monte Carlo dose calculation with VMC{sub ++}. Results: Both XMRT and IMRT solutions met all clinical criteria. Trade-offs were observed between improved dose uniformity to the primary target volume (PTV1) and increased dose to some of the surrounding healthy organs for XMRT compared to IMRT. On average, IMRT improved dose to the contralateral parotid gland and spinal cord while XMRT improved dose to the brainstem and mandible. Conclusion: Bi-energy XMRT optimization for H and N patients provides benefits in terms of improved dose uniformity to the primary target and reduced dose to some healthy structures, at the expense of increased dose to other healthy structures when compared with IMRT.

60 APPLIED LIFE SCIENCES↗

Efficient algorithm for locating and sizing series compensation devices in large power transmission grids: I. Model implementation

We explore optimization methods for planning the placement, sizing and operations of Flexible Alternating Current Transmission System (FACTS) devices installed to relieve transmission grid congestion. We limit our selection of FACTS devices to Series Compensation (SC) devices that can be represented by modification of the inductance of transmission lines. Our master optimization problem minimizes the l 1 norm of the inductance modification subject to the usual line thermal-limit constraints. We develop heuristics that reduce this non-convex optimization to a succession of Linear Programs (LP) which are accelerated further using cutting plane methods. The algorithm solves an instance of the MatPower Polish Grid model (3299 lines and 2746 nodes) in 40 seconds per iteration on a standard laptop—a speed up that allows the sizing and placement of a family of SC devices to correct a large set of anticipated congestions. We observe that our algorithm finds feasible solutions that are always sparse, i.e., SC devices are placed on only a few lines. In a companion manuscript, we demonstrate our approach on realistically-sized networks that suffer congestion from a range of causes including generator retirement. In this manuscript, we focus on the development of our approach, investigate its structure on a small test system subject to congestion from uniform load growth, and demonstrate computational efficiency on a realistically-sized network.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Spot Size Optimization of the Scorpius Accelerator

The Scorpius Linear Induction Accelerator (LIA) is being developed by the Advanced Sources and Diagnostics (ASD) Project and will complement other U1a capabilities by providing a multi-pulse, DARHT-class, flash-radiography capability. Commissioning activities are expected to end in Spring of 2026. All of the diagnostics discussed in this chapter are relevant to LIAs that generate high-current relativistic electron beams, and some diagnostics will also be relevant to other accelerators. The Scorpius Accelerator will be the world’s most advance LIA. Although significant technological advancements have been incorporated throughout the machine, many performance aspects of Scorpius will be similar to those of FXR and DARHT Axis-I and -II. These LIAs offer the same challenges, whether it is suppressing beam instabilities, improving diagnostic accuracy, validating computer simulations and beam tunes, or reducing beam-target interaction effects. However, the ultimate figure of merit is consistent between these LIAs, which is spot size and dose.

43 PARTICLE ACCELERATORS↗

Optimal Energy Dispatch of Distributed PVs for the Next Generation of Distribution Management Systems

Advanced Distribution Management Systems (ADMS) are being widely adopted by electric utilities for managing and optimizing the operations of their distribution systems. Distributed photovoltaic (DPV) systems with smart inverters can be controlled to adjust active power and reactive power outputs, and they are envisioned to become a part of (centrally or distributed) controllable assets managed by the ADMS for optimizing grid operations. This paper proposes an optimal energy dispatch strategy controlling DPV systems for regulating distribution voltages and achieving conservation voltage reduction. A convex optimization model is proposed with the use of linearized power flow, and the gradient projection algorithm is used to solve the optimal active power and reactive power outputs of smart inverters. The proposed optimal energy dispatch is implemented using an open-source ADMS platform, and simulation results have demonstrated the effectiveness of the proposed approach on improving distribution grid operations.

42 ENGINEERING↗

Optimal Control for Stochastic Delay Evolution Equations

In this paper, we investigate a class of infinite-dimensional optimal control problems, where the state equation is given by a stochastic delay evolution equation with random coefficients, and the corresponding adjoint equation is given by an anticipated backward stochastic evolution equation. We first prove the continuous dependence theorems for stochastic delay evolution equations and anticipated backward stochastic evolution equations, and show the existence and uniqueness of solutions to anticipated backward stochastic evolution equations. Then we establish necessary and sufficient conditions for optimality of the control problem in the form of Pontryagin’s maximum principles. To illustrate the theoretical results, we apply stochastic maximum principles to study two examples, an infinite-dimensional linear-quadratic control problem with delay and an optimal control of a Dirichlet problem for a stochastic partial differential equation with delay. Further applications of the two examples to a Cauchy problem for a controlled linear stochastic partial differential equation and an optimal harvesting problem are also considered.

97 MATHEMATICS AND COMPUTING↗

Highly sensitive SPR PCF biosensors based on Ag/TiN and Ag/ZrN configurations

We numerically present and analyze a surface plasmon resonance (SPR) bimetallic photonic crystal fiber (PCF) biosensor using full vectorial finite element method. The bimetallic configuration is based on Silver and (TiN or ZrN) as alternative plasmonic materials. The reported design relies on PCF with an external channel to house the analyte. The inner surface of the channel is coated by a silver layer followed by TiN or ZrN to protect the silver layer from oxidation. The proposed biosensor can be implemented to detect the unknown analytes, organic chemicals, biological analytes, and biomolecules via the sensitive resonant peaks to the analyte refractive index variations for both polarized modes with high linearity. In this study, a comparison is made between the Ag/TiN configuration and the Ag/ZrN counterpart in terms of sensitivity, losses, and linearity. The geometrical parameters are optimized to achieve high refractive index sensitivities of 7000 nm/RIU for quasi-transverse electric (TE) mode and 6900 nm/RIU for quasi-transverse magnetic (TM) mode, respectively for the Ag/TiN configuration. However, the quasi TE and quasi TM modes achieve high sensitivity of 5300 and 5400 nm/RIU, respectively using the Ag/ZrN configuration. The standard PCF fabrication technologies can be used to fabricate the proposed SPR PCF biosensors.

36 MATERIALS SCIENCE↗

A meshing framework for digital twins for extrusion based additive manufacturing

Additive manufacturing (AM) allows for manufacturing of complex three-dimensional geometries not typically realizable with standard manufacturing practices. The internal microstructure of AM components has a significant impact on mechanical, vibrational, and shock properties and permits richer design space when this is controllable. Due to complex interactions of internal geometry of an extrusion-based AM component, it is common practice to assume homogeneous behavior or to perform characterization testing on specific toolpath configurations. To avoid testing or material waste, it is necessary to develop a consistently accurate numerical simulation framework with relevant boundary value problems that can handle the complicated geometry of internal material microstructure present in AM components. Herein, a framework is proposed to directly create computational meshes suitable for finite element analysis (FEA) of the fine-scale features generated from extrusion-based AM tool paths to maintain a strong process–structure–property-performance linkage. This mesh can be manually or automatically analyzed using standard FEA simulations such as quasi-static preloading or modal analysis. The framework allows an in-silico assessment of a target AM geometry where fine-scale features greatly impact quantities of design interest such as in soft elastomeric lattices where toolpath infill can greatly influence the self-contact of a structure in compression, which we use as a motivating exemplar. This approach greatly reduces both time and resource waste present in traditional build and test design cycles for non-intuitive design spaces, and acts as a tool for use in the production of a key component of a digital twin, a mesh suitable for finite element analysis. In conclusion, it also further allows for the exploration of toolpath infill to optimize component properties beyond simple linear properties such as density and stiffness.

Additive manufacturing↗

Nonequilibrium steady state transport of collective-qubit system in strong coupling regime

We investigate the steady state photon transport in a nonequilibrium collective-qubit model. By adopting the noninteracting blip approximation, which is applicable in the strong photon–qubit coupling regime, we describe the essential contribution of indirect qubit–qubit interaction to the population distribution, mediated by the photonic baths. The linear relations of both the optimal flux and noise power with the qubits system size are obtained. Moreover, the inversed power-law style for the finite-size scaling of the optimal photon–qubit coupling strength is exhibited, which is proposed to be universal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Regulating the surface Pt coordination environment in the PtN overlayers on PtCuN hollow nanospheres for efficient oxygen reduction reaction

Engineering the surface Pt coordination environment is a promising strategy for promoting the kinetically sluggish oxygen reduction reaction (ORR) on Pt-based catalysts. Here, in this study, we achieve the compressive strain effect and electronic effect by Cu and N co-doping to synthesize the PtCuN hollow nanospheres with PtN overlayers (PtCuN@PtN HNSs) using a facile solvothermal synthesis. Electrochemical investigations show that the constructed disordered Pt–N coordination structures effectively facilitate the ORR and stabilize Pt atoms in the compressed lattice, whereas an excessive N-doping can lead to the formation of a structurally unstable Pt nitride phase. Theoretical analyses confirm that the oxygen reduction kinetics on the compressed PtN overlayers are regulated by a synergistic effect resulting from N-doping and lattice compression, circumventing the traditional linear scaling relationships (LSR). The optimized PtCuN@PtN HNSs, with the composition of PtCu 0.29 N 1.1 , demonstrate an area-specific activity of 1.98 mA cm –2 and a mass-specific activity of 1.81 A mg Pt –1 .

36 MATERIALS SCIENCE↗

A linear array of position-sensitive virtual Frisch-grid CdZnTe for low-energy gamma rays

Large-area arrays of position-sensitive virtual Frisch-grid CdZnTe (CZT) detectors with enhanced energy resolution and high efficiency have been proposed for spectroscopy and imaging of gamma-ray sources. We present the design and test results of a linear array of such detectors optimized for use in handheld instruments. There are several application areas for such detectors including uranium enrichment measurements, storage monitoring, dosimetry and other nuclear safeguards related tasks. The array configuration provides a large area with minimum number of readout channels, which potentially allows the developers to avoid using ASIC-based readout electronics by substituting it with hybrid preamplifiers followed by digitizers. The array prototype consists of 6 5 x 7 x 25 mm 3 CZT detectors side-oriented with respect to the source to achieve the maximum area. Each detector is furnished with 5 mm-wide charge-sensing pads placed near the anode. As an option, the pads can be grouped together to reduce the number of readout channels. The same can be done with the cathodes. The pads signals are converted into X–Y coordinates and combined with the cathode signals (for the Z coordinates) to give 3D positional information of each interaction point. This information is used to correct the response non-uniformity caused by material inhomogeneity, which allows the developers to use standard grade (unselected) CZT crystals, while retaining high spectroscopic performance.

3D position-sensitive virtual Frisch-grid detector↗

Electron-Affinity Time-Dependent Density Functional Theory: Formalism and Applications to Core-Excited States

The lack of particle-hole attraction and orbital relaxation within time-dependent density functional theory (TDDFT) lead to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n - 1-electron system as the reference, building orbital relaxation and a proper hole into the initial density. Here, our approach is an exact generalization of the static-exchange approximation that ameliorates the particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A unified form of low-energy nodal electronic interactions in hole-doped cuprate superconductors

Using angle resolved photoemission spectroscopy measurements of Bi 2 Sr 2 CaCu 2 O 8+δ over a wide range of doping levels, we present a universal form for the non-Fermi liquid electronic interactions in the nodal direction in the exotic normal state phase. It is described by a continuously varying power law exponent versus energy and temperature (hence named a Power Law Liquid or PLL), which with doping varies smoothly from a quadratic Fermi Liquid in the overdoped regime, to a linear Marginal Fermi Liquid at optimal doping, to a non-quasiparticle non-Fermi Liquid in the underdoped regime. The coupling strength is essentially constant across all regimes and is consistent with Planckian dissipation. Using the extracted PLL parameters we reproduce the experimental optics and resistivity over a wide range of doping and normal-state temperature values, including the T* pseudogap temperature scale observed in the resistivity curves. This breaks the direct link to the pseudogapping of antinodal spectral weight observed at similar temperature scales and gives an alternative direction for searches of the microscopic mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Utilising demand response for distribution service restoration to achieve grid resiliency against natural disasters

The increased frequency of power outages due to natural disasters in recent years has highlighted the urgency of enhancing distribution grid resilience. The effective distribution service restoration (DSR) is an important measure for a resilient distribution grid. In this work, the authors demonstrate that DSR can be significantly improved by leveraging the flexibility provided by the inclusion of demand response (DR). The authors propose a framework for this by considering integrated control of household-level flexible appliances to vary the load demand at the distribution-grid level to improve DSR. The overall framework of the proposed system is modelled as a three-step method considering three optimization problems to (i) calculate feasible controllable aggregated load range for each bus, (ii) determine candidate buses to perform DR and their target load demand, and (iii) maintain the load level in each house through home energy management during the DSR, considering uncertainties in load and solar generation sequentially. The optimization problems are formulated as linear programming, mixed-integer linear programming, and multistage stochastic programming (solved using the stochastic dual dynamic programming) models. Finally, case studies performed in the IEEE 123-node test feeder show improvements in resilience in terms of energy restored compared to the restoration process without DR.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An analytical approach to evaluating magnetic-field closure and topological changes in FRC devices

Here, we describe mathematical methods based on optimizing a modified non-linear flux function (MFF) to evaluate whether odd-parity perturbations affect the local closure of magnetic field lines in field-reversed configurations. Using the MFF methodology, quantitative formulas are derived that provide the shift of the field minimum and the threshold for field-line opening, a discontinuous change in field topology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Step towards Energy Efficient Computing: Redesigning a Hydrodynamic Application on CPU-GPU, In: 2014 IEEE 28th International Parallel and Distributed Processing Symposium

Power and energy consumption are becoming an increasing concern in high performance computing. Compared to multi-core CPUs, GPUs have a much better performance per watt. In this paper we discuss efforts to redesign the most computation intensive parts of BLAST, an application that solves the equations for compressible hydrodynamics with high order finite elements, using GPUs BLAST, Dobrev. In order to exploit the hardware parallelism of GPUs and achieve high performance, we implemented custom linear algebra kernels. We intensively optimized our CUDA kernels by exploiting the memory hierarchy, which exceed the vendor's library routines substantially in performance. We proposed an auto tuning technique to adapt our CUDA kernels to the orders of the finite element method. Compared to a previous base implementation, our redesign and optimization lowered the energy consumption of the GPU in two aspects: 60% less time to solution and 10% less power required. Compared to the CPU-only solution, our GPU accelerated BLAST obtained a 2.5× overall speedup and 1.42× energy efficiency (green up) using 4th order (Q_4) finite elements, and a 1.9× speedup and 1.27× green up using 2nd order (Q2) finite elements.

Computer Science↗