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 1,099 records · Page 61

On the Information in Spike Timing: Neural Codes Derived from Polychronous Groups

There is growing evidence regarding the importance of spike timing in neural information processing, with even a small number of spikes carrying information, but computational models lag significantly behind those for rate coding. Experimental evidence on neuronal behavior is consistent with the dynamical and state dependent behavior provided by recurrent connections. This motivates the minimalistic abstraction investigated in this paper, aimed at providing insight into information encoding in spike timing via recurrent connections. We employ information-theoretic techniques for a simple reservoir model which encodes input spatiotemporal patterns into a sparse neural code, translating the polychronous groups introduced by Izhike-vich into codewords on which we can perform standard vector operations. Finally, we show that the distance properties of the code are similar to those for (optimal) random codes. In particular, the code meets benchmarks associated with both linear classification and capacity, with the latter scaling exponentially with reservoir size.

computational modeling↗

Restoring Distribution System Under Renewable Uncertainty Using Reinforcement Learning

Distributed energy resources (DERs) in distribution systems, including renewable generation, micro-turbine, and energy storage, can be used to restore critical loads following extreme events to increase grid resiliency. However, properly coordinating multiple DERs in the system for multi-step restoration process under renewable uncertainty and fuel availability is a complicated sequential optimal control problem. Due to its capability to handle system non-linearity and uncertainty, reinforcement learning (RL) stands out as a potentially powerful candidate in solving complex sequential control problems. Moreover, the offline training of RL provides excellent action readiness during online operation, making it suitable to problems such as load restoration, where in-time, correct and coordinated actions are needed. In this study, a distribution system prioritized load restoration based on a simplified single-bus system is studied: with imperfect renewable generation forecast, the performance of an RL controller is compared with that of a deterministic model predictive control (MPC). Our experiment results show that the RL controller is able to learn from experience, adapt to the imperfect forecast information and provide a more reliable restoration process when compared with the baseline controller.

61 RADIATION PROTECTION AND DOSIMETRY↗

Generic Multi-Frequency Modelling of Converter-Connected Renewable Energy Generators Considering Frequency and Sequence Couplings

Frequency and sequence couplings can compromise the trustworthiness of multi-frequency models for converter-based systems. There have been effective attempts to address the couplings mainly by linearized averaged models. Only a few studies have been conducted on practical optimization of such models with enormous matrices and experimental results. This paper provides a generic theory for coupling patterns and proposes a multi-frequency modelling method to detect and address only the main couplings in the sequence domain for converter-connected renewable energy generators. The proposed generic model is based on empirical tests using small-signal perturbations and adopting Fourier transform on the switching converter response. The proposed theory and modelling methodology are verified using a 7MVA grid emulator for voltage perturbation tests on a 2MVA photo-voltaic converter. Accordingly, the couplings can exist in more generic forms, including multiples of perturbation and fundamental frequencies. To the best of our knowledge, the patterns with the multiples of the perturbation frequency have been overlooked in the literature. Furthermore, the mirror frequency concept is valid for all coupling patterns and is included in the proposed model. Besides, the proposed patterns and the environment noise levels have been practical criteria for selecting the main couplings.

frequency and sequence couplings↗

MO-FG-BRA-06: Material Optimisation for Nanoparticle-Sensitized Radiotherapy

Purpose: To evaluate the potential applications of nanoparticles of different elemental compositions as X-ray radiosensitizers. Methods: Geant4 Monte Carlo models were used to model the interactions of 10 nm radius nanoparticles of elements ranging from silicon to mercury with either monoenergetic keV X-rays or a 6 MV Linac spectrum. Secondary electron spectra emitted from nanoparticles and resulting radial doses were calculated. Potential biological impacts were evaluated through the use of the Local Effect Model. Results: For keV radiation, fluoresence photons and energetic photoelectrons deposit the majority of energy emitted by nanoparticles, over long ranges. However, short-range energy deposition is dominated by Auger electrons. While total energy deposition increases with atomic number, short-range energy deposition has a more complex behaviour, with local maxima around germanium (Z=32) and gadolinium (Z=64). For megavoltage radiation, Compton electrons deposit the majority of energy, but short- range depositions are dominated by interactions with secondary electrons. As a Result, energy deposition following ionising events following Linac irradiation is largely independent of atomic number.When analysed through the Local Effect Model, these effects translate into significant variations in predicted cell killing with atomic number for keV exposures, but comparatively limited effects for megavoltage exposures. Conclusions: Although heavier elements typically lead to greater increases in dose within target volumes, on the nanometre scale effects are typically dominated by low-energy secondary electrons, whose distribution has a complex relationship with both nanoparticle composition and beam energy.These results suggest that maximum radiosensitization may not necessarily be achieved through maximizing atomic number, but rather a balance between overall interaction cross-section and local dose distribution. The different behaviours seen when comparing keV and MV behaviours may also provide future mechanisms for separating physical and biological mechanisms of nanoparticle sensitisation, a key step in their development as a potential therapeutic agent.

60 APPLIED LIFE SCIENCES↗

Modelling and Analysis of the Corrosion Characteristics of Ferritic-Martensitic Steels in Supercritical Water

The dependencies of weight gain of 9-12 Cr ferritic-martensitic steels in supercritical water on each of seven principal independent variables (temperature, oxygen concentration, flow rate, exposure time, and key chemical composition and surface condition of steels) have been predicted using a supervised artificial neural network (ANN). The relative significance of each independent variable was uncovered by fuzzy curve analysis, which ranks temperature and exposure time as the most important. The optimized ANN, not only satisfactorily represents the experimentally-known non-linear relationships between the corrosion characteristics of F-M steels and the key independent variables (demonstrating the effectiveness of this technique), but also predicts and reveals that the effects of oxygen concentration on the weight gains, to a certain degree, is influenced by the flow rate and temperature. Finally, according to the ANN predicted-results, departure of oxidation kinetics from the parabolic law, and basic cause of chromium content in steel substrate influencing the corrosion rate, and the synergetic effects of dissolved oxygen concentration, flow rate, and temperature, are discussed and analyzed.

36 MATERIALS SCIENCE↗

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗

Restoring Distribution System Under Renewable Uncertainty Using Reinforcement Learning: Preprint

Distributed energy resources (DER) in distribution systems, including renewable generation, micro-turbine, and energy storage, can be used to restore critical loads following extreme events to increase grid resiliency. However, properly coordinating multiple DERs in the system for multi-step restoration process under renewable uncertainty and fuel availability is a complicated sequential optimal control problem. Due to its capability to handle system non-linearity and uncertainty, reinforcement learning (RL) stands out as a potentially powerful candidate in solving complex sequential control problems. Moreover, the offline training of RL provides excellent action readiness during online operation, making it suitable to problems such as load restoration, where in-time, correct and coordinated actions are needed. In this study, a distribution system prioritized load restoration based on a simplified single-bus system is studied: with imperfect renewable generation forecast, the performance of an RL controller is compared with that of a deterministic model predictive control (MPC). Our experiment results show that the RL controller is able to learn from experience, adapt to the imperfect forecast information and provide a more reliable restoration process when compared with the baseline MPC controller.

61 RADIATION PROTECTION AND DOSIMETRY↗

Restoring Distribution System Under Renewable Uncertainty Using Reinforcement Learning

Distributed energy resources (DER) in distribution systems, including renewable generation, micro-turbine, and energy storage, can be used to restore critical loads following extreme events to increase grid resiliency. However, properly coordinating multiple DERs in the system for multi-step restoration process under renewable uncertainty and fuel availability is a complicated sequential optimal control problem. Due to its capability to handle system non-linearity and uncertainty, reinforcement learning (RL) stands out as a potentially powerful candidate in solving complex sequential control problems. Moreover, the offline training of RL provides excellent action readiness during online operation, making it suitable to problems such as load restoration, where in-time, correct and coordinated actions are needed. In this study, a distribution system prioritized load restoration based on a simplified single-bus system is studied: with imperfect renewable generation forecast, the performance of an RL controller is compared with that of a deterministic model predictive control (MPC). Our experiment results show that the RL controller is able to learn from experience, adapt to the imperfect forecast information and provide a more reliable restoration process when compared with the baseline MPC controller.

61 RADIATION PROTECTION AND DOSIMETRY↗

Using Contact to Increase Robot Performance for Glovebox D and D Tasks - 18227

Glovebox decommissioning tasks usually require manipulating relatively heavy objects in a highly constrained environment. Thus, contact with the surroundings becomes inevitable. In order to allow the robot to interact with the environment in a natural way, we present a contact-implicit motion planning framework. This framework enables the system, without the specification in advance of a contact plan, to make and break contacts to maintain stability while performing a manipulation task. In this method, we use linear complementarity constraints to model rigid body contacts and find a locally optimal solution for joint displacements and magnitudes of support forces. Then, joint torques are calculated such that the support forces have the highest priority. We evaluate our framework in a 2.5D, quasi-static simulation in which a humanoid robot with planar arms manipulates a heavy object. Our results suggest that the proposed method provides the robot with the ability to balance itself by generating support forces on the environment while simultaneously performing the manipulation task. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Quandary

Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.

Petersson, NilsA.↗

Juqbox.jl

The Juqbox.jl package implements functionality for solving the quantum optimal control problem for realizing logical gates in closed quantum systems. The dynamics of the quantum system is modeled by Schroedinger's equation, which takes to form of a linear system of ordinary differential equations (ODE). Juqbox.jl solves this ODE by numerical time stepping and applies a gradient-based optimization technique to determine control pulses for driving an initial state to a final state, according to the desired logical gate transformation. To evaluate the gradient, Juqbox.jl applies the ``first discretize, then optimize'' approach based on a discrete adjoint time stepping technique. The actual optimization is performed by the open source Ipopt libaray. Juqbox.jl is written in the Julia programming language which, among many other features, provides a convenient interface to the Ipopt library.

PETERSSON, NILSA.↗

OptiBench: An Optimization Benchmark Tool for Renewable Energy Problems

We propose a benchmark framework and visualization tool, OptiBench, for analyzing the performance of state-of-the-art optimization solvers across a variety of optimization problems in renewable energy research. Our framework is designed from the ground up in the Julia programming language and enables analysis at scale on high performance computing (HPC) systems. Our visualization tool allows effortless evaluation of optimization solver performance, robustness, and accuracy through intuitive plots, e.g., performance profiles, heat maps, and distribution plots. We have tested three benchmark suites relevant to the modeling of renewable energy systems, viz., CUTEst, PGLib-OPF, and WaterTAP water treatment optimization problems. We illustrate benchmarking of CUTEst using OptiBench on the National Renewable Energy Laboratory's (NREL) HPC Kestrel. Our findings indicate that MA57 HSL linear solver demonstrated the best overall performance for an experimental IPOPT implementation. Our work is ongoing and we intend to add support for more optimization solvers and benchmark test suites in the future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Urban-Scale Control of School Bus Fleet Charging and Discharging Strategies Using Single and Multi-Stage Optimization

This paper presents a dual-strategy approach to optimizing charging and discharging schedules for school bus fleets, using the limited charging infrastructure effectively. We aim to ensure that each bus is fully charged for daily operations and aids in grid stability during peak demand. The first strategy utilizes linear programming to schedule overnight charging at available station sockets and strategic discharging during peak periods, efficiently coordinating limited resources. The second strategy employs metaheuristic techniques for continuous optimization, focusing on precise power requirements and offering greater flexibility than the linear model.

Selim, Alaa↗

High-contrast imaging with an arbitrary aperture: Active compensation of aperture discontinuities

We present a new method to achieve high-contrast images using segmented and/or on-axis telescopes. Our approach relies on using two sequential deformable mirrors (DMs) to compensate for the large amplitude excursions in the telescope aperture due to secondary support structures and/or segment gaps. In this configuration the parameter landscape of DM surfaces that yield high-contrast point-spread functions is not linear, and nonlinear methods are needed to find the true minimum in the optimization topology. We solve the highly nonlinear Monge-Ampere equation that is the fundamental equation describing the physics of phase-induced amplitude modulation. We determine the optimum configuration for our two sequential DM system and show that high-throughput and high-contrast solutions can be achieved using realistic surface deformations that are accessible using existing technologies. We name this process Active Compensation of Aperture Discontinuities (ACAD). We show that for geometries similar to the James Webb Space Telescope, ACAD can attain at least 10{sup –7} in contrast and an order of magnitude higher for both the future extremely large telescopes and on-axis architectures reminiscent of the Hubble Space Telescope. We show that the converging nonlinear mappings resulting from our DM shapes actually damp near-field diffraction artifacts in the vicinity of the discontinuities. Thus, ACAD actually lowers the chromatic ringing due to diffraction by segment gaps and struts while not amplifying the diffraction at the aperture edges beyond the Fresnel regime. This outer Fresnel ringing can be mitigated by properly designing the optical system. Consequently, ACAD is a true broadband solution to the problem of high-contrast imaging with segmented and/or on-axis apertures. We finally show that once the nonlinear solution is found, fine tuning with linear methods used in wavefront control can be applied to further contrast by another order of magnitude. Generally speaking, the ACAD technique can be used to significantly improve a broad class of telescope designs for a variety of problems.

79 ASTRONOMY AND ASTROPHYSICS↗

High-Throughput Characterization and Optimization of Polyamide Hydrolase Activity Using Open Port Sampling Interface Mass Spectrometry

Enzymatic biodegradation of polymers, such as polyamides (PA), has the potential to cost-effectively reduce plastic waste, but enhancements in degradation efficiency are needed. Engineering enzymes through directed evolution is one pathway toward identification of critical domains needed for improving activity. However, screening such enzymatic libraries (100s-to-1000s of samples) is time-consuming. Here we demonstrate the use of robotic autosampler (PAL) and immediate drop on demand technology (I.DOT) liquid handling systems coupled with open-port sampling interface-mass spectrometry (OPSI-MS) to screen for PA6 and PA66 hydrolysis by 6-aminohexanoate-oligomer endo-hydrolase (nylon hydrolase, NylC) in a high-throughput (8–20 s/sample) manner. The OPSI-MS technique required minimal sample preparation and was amenable to 96-well plate formats for automated processing. Enzymatic hydrolysis of PA characteristically produced soluble linear oligomer products that could be identified by OPSI-MS. Incubation temperatures and times were optimized for PA6 (65 °C, 24 h) and PA66 (75 °C, 24 h) over 108 experiments. In addition, the I.DOT/OPSI-MS quantified production of PA6 linear dimer (8.3 ± 1.6 μg/mL) and PA66 linear monomer (13.5 ± 1.5 μg/mL) by NylC with a lower limit of detection of 0.029 and 0.032 μg/mL, respectively. For PA6 and PA66, linear oligomer production corresponded to 0.096 ± 0.018% and 0.204 ± 0.028% conversion of dry pellet mass, respectively. The developed methodology is expected to be utilized to assess enzymatic hydrolysis of engineered enzyme libraries, comprising hundreds to thousands of individual samples.

47 OTHER INSTRUMENTATION↗

Parallel interior-point solver for block-structured nonlinear programs on SIMD/GPU architectures

Here, we investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses SIMD/GPU accelerators locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations. We demonstrate our method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the reduction method is 50x faster than the sparse linear solver HSL MA57 running in serial on the CPU, and 6x faster than Pardiso running in parallel on CPU on the same number of processes.

97 MATHEMATICS AND COMPUTING↗