Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear programming”

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 253 records · Page 14

The 2015 Nonlinear Mechanics and Dynamics Research Institute

The Nonlinear Mechanics and Dynamics (NOMAD) Research Institute is a six week long collaborative research program for graduate students from across the world. The 2015 NOMAD Research Institute was hosted jointly by Sandia National Laboratories and the University of New Mexico, and featured 24 graduate students working on seven different research projects. These projects included: developing experimental strategies for studying the dynamics of nonlinear systems, a numerical round robin for predicting the response of a jointed system, quantification of uncertainty in a lap joint, assessment of experimental substructuring methods, a study of stress waves propagating through jointed interfaces, structural design optimization with joints, and the nonlinear system identification of MEMS devices. This report details both the technical research and the programmatic organization of the 2015 NOMAD Research Institute.

42 ENGINEERING↗

WE-FG-207B-02: Material Reconstruction for Spectral Computed Tomography with Detector Response Function

Purpose: Different from the conventional computed tomography (CT), spectral CT based on energy-resolved photon-counting detectors is able to provide the unprecedented material composition. However, an important missing piece for accurate spectral CT is to incorporate the detector response function (DRF), which is distorted by factors such as pulse pileup and charge-sharing. In this work, we propose material reconstruction methods for spectral CT with DRF. Methods: The polyenergetic X-ray forward model takes the DRF into account for accurate material reconstruction. Two image reconstruction methods are proposed: a direct method based on the nonlinear data fidelity from DRF-based forward model; a linear-data-fidelity based method that relies on the spectral rebinning so that the corresponding DRF matrix is invertible. Then the image reconstruction problem is regularized with the isotropic TV term and solved by alternating direction method of multipliers. Results: The simulation results suggest that the proposed methods provided more accurate material compositions than the standard method without DRF. Moreover, the proposed method with linear data fidelity had improved reconstruction quality from the proposed method with nonlinear data fidelity. Conclusion: We have proposed material reconstruction methods for spectral CT with DRF, whichprovided more accurate material compositions than the standard methods without DRF. Moreover, the proposed method with linear data fidelity had improved reconstruction quality from the proposed method with nonlinear data fidelity. Jiulong Liu and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000), and the Shanghai Pujiang Talent Program (#14PJ1404500).

60 APPLIED LIFE SCIENCES↗

Optimal Operation of Soft Open Points in Active Distribution Networks Under Three-Phase Unbalanced Conditions

The asymmetric integration of distributed generators (DGs) exacerbates the three-phase unbalanced condition in distribution systems. The serious unbalanced operation causes inefficient utilization of network assets and security risks in the system. Soft open point (SOP) is a flexible power electronic device which can achieve accurate active and reactive power flow control to balance the power flow among phases. This paper proposes an SOPs-based operation strategy for unbalanced active distribution networks (ADNs). By regulating the operation of SOPs, the strategy can reduce power losses and simultaneously mitigate the three-phase unbalance of the upper-level grid. Semidefinite programming (SDP) relaxation is advocated to convert the original non-convex, nonlinear optimization model into an SDP formulation, which can be efficiently solved to meet the requirement of rapid adjustment. Case studies are conducted on the modified IEEE 33-node and IEEE 123-node distribution system to verify the effectiveness and efficiency of the proposed strategy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Vibrational spectroscopic and non-linear optical activity studies on nicotinanilide : A DFT approach

The molecular structure of nicotinanilide was optimized by the DFT/B3LYP method with cc-pVTZ basis set using Gaussian 09 program. The first order hyperpolarizability of the molecule was calculated, which exhibits the higher nonlinear optical activity. The natural bond orbital analysis confirms the presence of intramolecular charge transfer and the hydrogen bonding interaction, which leads to the higher nonlinear optical activity of the molecule. The Frontier molecular orbitals analysis of the molecule shows that the delocalization of electron density occurs within the molecule. The lower energy gap indicates that the hydrogen bond formation between the charged species. The vibrational frequencies were calculated and assigned on the basis of potential energy distribution calculation using the VEDA 4.0 program and the corresponding vibrational spectra were simulated. Hence, the nicotinanilide molecule can be a good candidate for second-order NLO material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Next-generation Algorithms for Assessing Infrastructure Vulnerability and Optimizing System Resilience

This report summarizes the work performed under the project project Next-Generation Algorithms for Assessing Infrastructure Vulnerability and Optimizing System Resilience. The goal of the project was to improve mathematical programming-based optimization technology for infrastructure protection. In general, the owner of a network wishes to design a network a network that can perform well when certain transportation channels are inhibited (e.g. destroyed) by an adversary. These are typically bi-level problems where the owner designs a system, an adversary optimally attacks it, and then the owner can recover by optimally using the remaining network. This project funded three years of Deon Burchett's graduate research. Deon's graduate advisor, Professor Jean-Philippe Richard, and his Sandia advisors, Richard Chen and Cynthia Phillips, supported Deon on other funds or volunteer time. This report is, therefore. essentially a replication of the Ph.D. dissertation it funded in a format required for project documentation. The thesis had some general polyhedral research. This is the study of the structure of the feasible region of mathematical programs, such as integer programs. For example, an integer program optimizes a linear objective function subject to linear constraints, and (nonlinear) integrality constraints on the variables. The feasible region without the integrality constraints is a convex polygon. Careful study of additional valid constraints can significantly improve computational performance. Here is the abstract from the dissertation: We perform a polyhedral study of a multi-commodity generalization of variable upper bound flow models. In particular, we establish some relations between facets of single- and multi- commodity models. We then introduce a new family of inequalities, which generalizes traditional flow cover inequalities to the multi-commodity context. We present encouraging numerical results. We also consider the directed edge-failure resilient network design problem (DRNDP). This problem entails the design of a directed multi-commodity flow network that is capable of fulfilling a specified percentage of demands in the event that any G arcs are destroyed, where G is a constant parameter. We present a formulation of DRNDP and solve it in a branch-column-cut framework. We present computational results.

97 MATHEMATICS AND COMPUTING↗

ECP Milestone Report: Identify initial kernels, bake-off problems (benchmarks) and miniapps (CEED-MS6)

The CEED co-design center is developing a number of kernels, bake-off / benchmark problems (BPs) and mini-applications (miniapps) that capture the unique requirements of high-order finite element algorithms and use them to influence standards, best practices and vendors to generate more portable, better performing code across diverse hardware designs. This is a main activity in CEED’s Hardware, Applications and Finite Element thrusts. CEED’s miniapps, a.k.a. CEEDlings, are small, standalone applications - surrogates for sub-components of interest from the selected first wave of applications (the MARBL/LLNLApp and ExaSMR projects). The BPs are even smaller standalone programs that solve a simple physical problem, including discretization, linear and/or nonlinear solve, and potentially a relatively small number of time steps. A BP may be a proxy/surrogate for part of a miniapp or another problem of general interest. In CEED, both miniapps and BPs will be used extensively in the project to evaluate and compare the performance of algorithms in specific contexts. The parts of the miniapps and BPs that play an important role in their overall performance will be identified and separated as standalone kernels. These kernels and their interaction inside BPs, miniapps and the actual applications, will be the main focus of our optimization efforts. All three standalone pieces of software (miniapps, BPs, and kernels) will serve as a basis for engaging and collaborating with hardware vendors and software technologies (ST) projects. These partners can start at the bottom of the hierarchy, e.g. with the kernels, and easily transition to BPs and miniapps as all three families will be produced in the same CEED environment. In this milestone we identified the initial kernels, bake-off problems and miniapps and delivered software and documentation for them through the new CEED website, http://ceed.exascaleproject.org. We specifically identified the Nekbone miniapp for the ExaSMR application and developed a new miniapp, Laghos, for the MARBL/LLNLApp application. We also formulated four initial benchmark problems (BP1- BP4) with clearly defined performance metrics and performed an initial bake-off comparison between the Nek5000 (spectral element) and MFEM (high-order finite element) technologies. As part of the milestone, we also engaged vendors (AMD) and ST projects (MPICH, STRUMPACK). In this document we are reporting some details and results from these R&D activities, as well as additional project-wide activities performed in Q3 of FY17, including: making CEED software available on GitHub, developing the project website, and the preparation of CEED’s first annual meeting.

97 MATHEMATICS AND COMPUTING↗

Modeling and Detection of Future Cyber-Enabled DSM Data Attacks

Demand-Side Management (DSM) is an essential tool to ensure power system reliability and stability. In future smart grids, certain portions of a customer’s load usage could be under the automatic control of a cyber-enabled DSM program, which selectively schedules loads as a function of electricity prices to improve power balance and grid stability. In this scenario, the security of DSM cyberinfrastructure will be critical as advanced metering infrastructure and communication systems are susceptible to cyber-attacks. Such attacks, in the form of false data injections, can manipulate customer load profiles and cause metering chaos and energy losses in the grid. The feedback mechanism between load management on the consumer side and dynamic price schemes employed by independent system operators can further exacerbate attacks. To study how this feedback mechanism may worsen attacks in future cyber-enabled DSM programs, we propose a novel mathematical framework for (i) modeling the nonlinear relationship between load management and real-time pricing, (ii) simulating residential load data and prices, (iii) creating cyber-attacks, and (iv) detecting said attacks. In this framework, we first develop time-series forecasts to model load demand and use them as inputs to an elasticity model for the price-demand relationship in the DSM loop. This work then investigates the behavior of such a feedback loop under intentional cyber-attacks. We simulate and examine load-price data under different DSM-participation levels with three types of random additive attacks: ramp, sudden, and point attacks. We conduct two investigations for the detection of DSM attacks. The first studies a supervised learning approach, with various classification models, and the second studies the performance of parametric and nonparametric change point detectors. Results conclude that higher amounts of DSM participation can exacerbate ramp and sudden attacks leading to better detection of such attacks, especially with supervised learning classifiers. We also find that nonparametric detection outperforms parametric for smaller user pools, and random point attacks are the hardest to detect with any method.

97 MATHEMATICS AND COMPUTING↗

First Measurements of Nonlinear Decoherence in the IOTA Ring

The Integrable Optics Test Accelerator (IOTA), at Fermi National Laboratory is aimed at testing nonlinear optics for the next generation of high intensity rings. Through use of a special magnetic element the ring is designed to induce a large tune spread with amplitude while maintaining integrable motion. This will allow for the suppression of instabilities in high-intensity beams without significant reduction in dynamic aperture. One important aspect of this is the nonlinear decoherence that occurs when a beam is injected off axis or receives a transverse kick while circulating in the ring. This decoherence has been studied in detail, with simulations, for protons in IOTA both with and without space-charge. However, it has yet to be demonstrated experimentally. During the first phase of the IOTA experimental program, the ring is operated with 100 MeV electrons, allowing for the study of nonlinear optics without the complications introduced by space charge. Here we present measurements taken during the IOTA commissioning, and an analysis of the results.

43 PARTICLE ACCELERATORS↗

Sequence of polyhedral relaxations for nonlinear univariate functions

Here, given a nonlinear, univariate, bounded, and differentiable function f(x), this article develops a sequence of Mixed Integer Linear Programming (MILP) and Linear Programming (LP) relaxations that converge to the graph of f(x) and its convex hull, respectively. Theoretical convergence of the sequence of relaxations to the graph of the function and its convex hull is established. For nonlinear non-convex optimization problems, the relaxations presented in this article can be used to construct tight MILP and LP relaxations. These MILP and the LP relaxations can also be used with MILP-based and spatial branch-and-bound based global optimization algorithms, respectively.

42 ENGINEERING↗

Bi-Level Volt-VAR Optimization to Coordinate Smart Inverters with Voltage Control Devices

Conservation voltage reduction (CVR) uses voltage control to reduce customer power demand and help improve the operational efficiency of the electric power distribution system. With the advent of advanced distribution management systems (ADMS), the underlying data and control-rich platform can be utilized to enable a centralized Volt-VAR optimization (VVO) scheme that includes smart inverter controls along with legacy voltage control devices such as voltage regulators and capacitor banks to ensure a maximum CVR benefit. The objective of this paper to present a VVO approach that controls system’s legacy voltage control devices and coordinates their operation with smart inverter control. An optimal power flow (OPF) formulation is proposed by developing linear and nonlinear power flow approximations for a three-phase unbalanced electric power distribution system. A bi-level VVO approach is proposed where, Level-1 optimizes the control of legacy devices and smart inverters using a linear approximate three-phase power flow. In Level-2, the control parameters for smart inverters are adjusted to obtain an optimal and feasible solution by solving the approximate nonlinear OPF model. Level-1 is modeled as a Mixed Integer Linear Program (MILP) while Level-2 as a Quadratically Constrained Quadratic Program (QCQP). The proposed approach is validated using IEEE test cases and the results demonstrate the applicability of the framework in achieving the CVR objective.

Jha, Rahul Ranjan↗

Cultural Perspectives of the 2015 Nonlinear Mechanics and Dynamics Summer Research Institute

This paper discusses the results of a study to determine the impact of culture on engineering. The study took place during the 2015 Nonlinear Mechanics and Dynamics Summer Research Institute, a six-week research program sponsored by Sandia National Laboratories and the University of New Mexico consisting of 24 graduate students participating in seven different projects. Twenty-two of the participants and two of the mentors were interviewed to study the effects of cultural background on engineering processes and interactions. The results of this study indicate that cultural differences drive engineering practices.

42 ENGINEERING↗

Cultural Perspectives of the 2016 Nonlinear Mechanics and Dynamics Research Institute

This report analyzes the results of a study on culture and its capability to influence research. The study occurred during the 2016 Nonlinear Mechanics and Dynamics Summer Research Institute, a six-week research program sponsored by Sandia National Laboratories and the University of New Mexico consisting of 27 graduate students participating in ten different projects. Two separate surveys were administered at the beginning and end of the Institute, in addition to interviews and observation, in order to study the effects of various cultural factors on engineering processes and maintaining professional interactions. The results of this study indicate that cultural differences are not a significant barrier to engineering progress and most cultural issues are minor. A variety of cultures instead provide new perspectives, advancing universal understanding.

42 ENGINEERING↗

Resilience Project Implementation Memo for the Town of Stowe, Vermont

This Resilience Project Implementation Process Memo was developed with and for the Town of Stowe Electric Department (SED) under the DOE C2C Expert Match Program. It is intended to outline the steps of a process, albeit nonlinear and iterative, to plan, design, and implement community resilience hubs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Probing the deep nonlinear stage of the ablative Rayleigh-Taylor instability in indirect drive experiments on the National Ignition Facility

Academic tests in physical regimes not encountered in Inertial Confinement Fusion will help to build a better understanding of hydrodynamic instabilities and constitute the scientifically grounded validation complementary to fully integrated experiments. Under the National Ignition Facility (NIF) Discovery Science program, recent indirect drive experiments have been carried out to study the ablative Rayleigh-Taylor Instability (RTI) in transition from weakly nonlinear to highly nonlinear regime [A. Casner et al., Phys. Plasmas 19, 082708 (2012)]. In these experiments, a modulated package is accelerated by a 175 eV radiative temperature plateau created by a room temperature gas-filled platform irradiated by 60 NIF laser beams. The unique capabilities of the NIF are harnessed to accelerate this planar sample over much larger distances (≃1.4 mm) and longer time periods (≃12 ns) than previously achieved. This extended acceleration could eventually allow entering into a turbulent-like regime not precluded by the theory for the RTI at the ablation front. Simultaneous measurements of the foil trajectory and the subsequent RTI growth are performed and compared with radiative hydrodynamics simulations. We present RTI growth measurements for two-dimensional single-mode and broadband multimode modulations. The dependence of RTI growth on initial conditions and ablative stabilization is emphasized, and we demonstrate for the first time in indirect-drive a bubble-competition, bubble-merger regime for the RTI at ablation front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrogen-bond Specific Materials Modification in Group IV Semiconductors

Semiconductor dielectric crystals consist of two fundamental components: lattice atoms and electrons. The former component provides a crystalline structure that can be disrupted by various defects or the presence of an interface, or by transient oscillations known as phonons. The latter component produces an energetic structure that is responsible for the optical and electronic properties of the material, and can be perturbed by lattice defects or by photo-excitation. Over the period of this project, August 15, 1999 to March 31, 2015, a persistent theme has been the elucidation of the fundamental role of defects arising from the presence of radiation damage, impurities (in particular, hydrogen), localized strain or some combination of all three. As our research effort developed and evolved, we have experienced a few title changes, which reflected this evolution. Throughout the project, ultrafast lasers usually in a pump-probe configuration provided the ideal means to perturb and study semiconductor crystals by both forms of excitation, vibrational (phonon) and electronic (photon). Moreover, we have found in the course of this research that there are many interesting and relevant scientific questions that may be explored when phonon and photon excitations are controlled separately. Our early goals were to explore the dynamics of bond-selective vibrational excitation of hydrogen from point defects and impurities in crystalline and amorphous solids, initiating an investigation into the behavior of hydrogen isotopes utilizing a variety of ultrafast characterization techniques, principally transient bleaching spectroscopy to experimentally obtain vibrational lifetimes. The initiative could be divided into three related areas: (a) investigation of the change in electronic structure of solids due to the presence of hydrogen defect centers, (b) dynamical studies of hydrogen in materials and (c) characterization and stability of metastable hydrogen impurity states under transient compression. This research focused on the characterization of photon and ion stimulated hydrogen related defect and impurity reactions and migration in solid state matter, which requires a detailed understanding of the rates and pathways of vibrational energy flow, of the transfer channels and of the coupling mechanisms between local vibrational modes (LVMs) and phonon bath as well as the electronic system of the host material. It should be stressed that researchers at Vanderbilt and William and Mary represented a unique group with a research focus and capabilities for low temperature creation and investigation of such material systems. Later in the program, we carried out a vigorous research effort addressing the roles of defects, interfaces, and dopants on the optical and electronic characteristics of semiconductor crystals, using phonon generation by means of ultrafast coherent acoustic phonon (CAP) spectroscopy, nonlinear characterization using second harmonic generation (SHG), and ultrafast pump-and-probe reflectivity and absorption measurements. This program featured research efforts from hydrogen defects in silicon alone to other forms of defects such as interfaces and dopant layers, as well as other important semiconducting systems. Even so, the emphasis remains on phenomena and processes far from equilibrium, such as hot electron effects and travelling localized phonon waves. This program relates directly to the mission of the Department of Energy. Knowledge of the rates and pathways of vibrational energy flow in condensed matter is critical for understanding dynamical processes in solids including electronically, optically and thermally stimulated defect and impurity reactions and migration. The ability to directly probe these pathways and rates allows tests of theory and scaling laws at new levels of precision. Hydrogen embedded in model crystalline semiconductors and metal oxides is of particular interest, since the associated local mode can be excited cleanly, and is usually well-separated in energy from the phonon bath. These basic dynamical studies have provided new insights for example into the fundamental mechanisms that control proton diffusion in these oxides. This area of materials science has largely fulfilled its promise to identify degradation mechanisms in electronic and optoelectronic devices, and to advance solid oxide proton conductors for fuel cells, gas sensors and proton-exchange membrane applications. It also provides the basis for innovations in materials synthesis involving atomic-selective diffusion and desorption.

36 MATERIALS SCIENCE↗

Final Report for Statistical Nonlinear Optics of High Energy Density Plasmas: The Physics of Multiple Crossing Laser Beams

The various implementations of the STUD pulse program (spike trains of uneven duration and delay) for LPI (laser-plasma instability) control were studied in depth, and novel regimes were found. How to generate STUD pulses with large time-bandwidth products, how to measure their optical scattering signatures, and how to experimentally demonstrate their usefulness were explored. Theoretical and numerical studies were conducted on Stimulated Brillouin Scattering (SBS) and Crossed Beam Energy Transfer (CBET) including statistical models. We established how LPI can be tamed and gain democratized in space and time. Implementing STUD pulses on NIF was also studied. Future high rep rate lasers and fast diagnostics will aid in the adoption of the whole STUD pulse program for LPI control in High Energy Density Plasmas (HEDP).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatiotemporal control of laser intensity using differentiable programming

Optical techniques for spatiotemporal control can produce laser pulses with custom amplitude, phase, or polarization structure. In nonlinear optics and plasma physics, the use of structured pulses typically follows a forward design approach, in which the efficacy of a known structure is analyzed for a particular application. Inverse approaches, in contrast, enable the discovery of new structures with the potential for superior performance. Here, an implementation of the unidirectional pulse propagation equation that supports automatic differentiation is combined with gradient-based optimization to design structured pulses with features that are advantageous for a range of nonlinear optical and plasma-based applications: (1) a longitudinally uniform intensity over an extended region, (2) a superluminal intensity peak that travels many Rayleigh ranges with constant duration, spot size, and amplitude, and (3) a laser pulse that ionizes a gas to form a uniform column of plasma. In the final case, optimizing the full spatiotemporal structure improves the performance by a factor of 15 compared to optimizing only spatial or only temporal structure, highlighting the advantage of spatiotemporal control.

automatic differentiation↗

Applications of the Dulmage–Mendelsohn decomposition for debugging nonlinear optimization problems

Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.

42 ENGINEERING↗