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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Semidefinite programming algorithm for the quantum mechanical bootstrap

Here, we present a semidefinite programming algorithm to find eigenvalues of Schrödinger operators within the bootstrap approach to quantum mechanics. The bootstrap approach involves two ingredients: a nonlinear set of constraints on the variables (expectation values of operators in an energy eigenstate), plus positivity constraints (unitarity) that need to be satisfied. By fixing the energy we linearize all the constraints and show that the feasibility problem can be presented as an optimization problem for the variables that are not fixed by the constraints and one additional slack variable that measures the failure of positivity. To illustrate the method we are able to obtain high-precision, sharp bounds on eigenenergies for arbitrary confining polynomial potentials in one dimension.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

SUNDIALS Multiphysics+MPIManyVector Performance Testing

In this report we document performance test results on a SUNDIALS-based multiphysics demonstration application. We aim to assess the large-scale parallel performance of new capabilities that have been added to the SUNDIALS suite of time integrators and nonlinear solvers in recent years under funding from both the Exascale Computing Project (ECP) and the Scientific Discovery through Advanced Scientific (SciDAC) program, specifically: 1. SUNDIALS’ new MPIManyVector module, that allows extreme flexibility in how a solution "vector" is staged on computational resources. 2. ARKode’s new multirate integration module, MRIStep, allowing high-order accurate calculations that subcycle "fast" processes within "slow" ones. 3. SUNDIALS’ new flexible linear solver interfaces, that allow streamlined specification of problem-specific linear solvers. 4. SUNDIALS’ new N_Vector additions of "fused" vector operations (to increase arithmetic intensity) and separation of reduction operations into "local" and "global" versions (to reduce latency by combining multiple reductions into a single MPI_Allreduce call). We anticipate that subsequent reports will extend this work to investigate a variety of other new features, including SUNDIALS’ generic SUNNonlinearSolver interface and accelerator-enabled N Vector modules, and upcoming MRIStep extensions to support custom "fast" integrators (that leverage problem structure) and IMEX integration of the "slow" time scale (to add diffusion).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Riemannian Optimization Applied to AC Optimal Power Flow

The nonlinear, nonconvex AC optimal power flow problem is of growing importance as the nature of the power grid evolves. This problem can be difficult to solve for interior point methods. However, the advent of optimization algorithms over smooth Riemannian manifolds presents an alternative approach. The nonlinear, nonconvex constraints in the AC power flow problem form an embedded submanifold of Euclidean space. In this paper, the authors explore the performance of Riemannian optimization algorithms for the ACOPF problem where the optimization is performed directly on the AC power flow manifold. This is done by using the Julia programming language and the Julia packages PowerModels.jl and Manopt.jl.

AC optimal power flow↗
DOE OSTI2017

Generation Expansion Planning With Large Amounts of Wind Power via Decision-Dependent Stochastic Programming

Power generation expansion planning needs to deal with future uncertainties carefully, given that the invested generation assets will be in operation for a long time. Many stochastic programming models have been proposed to tackle this challenge. However, most previous works assume predetermined future uncertainties (i.e., fixed random outcomes with given probabilities). In several recent studies of generation assets' planning (e.g., thermal versus renewable), new findings show that the investment decisions could affect the future uncertainties as well. To this end, this paper proposes a multistage decision-dependent stochastic optimization model for long-term large-scale generation expansion planning, where large amounts of wind power are involved. In the decision-dependent model, the future uncertainties are not only affecting but also affected by the current decisions. In particular, the probability distribution function is determined by not only input parameters but also decision variables. To deal with the nonlinear constraints in our model, a quasi-exact solution approach is then introduced to reformulate the multistage stochastic investment model to a mixed-integer linear programming model. The wind penetration, investment decisions, and the optimality of the decision-dependent model are evaluated in a series of multistage case studies. The results show that the proposed decision-dependent model provides effective optimization solutions for long-term generation expansion planning.

Zhan, Yiduo↗
DOE OSTI2016

Vibrational spectroscopic, structural and nonlinear optical activity studies on 2-amino-3-chloro-5-trifluoromethyl pyridine: A DFT approach

The conformational analysis was carried out for 2-amino-3-chloro-5-trifluoromethylpyridine using potential energy surface (PES) scan and the most stable optimized conformer was predicted. The theoretical vibrational frequencies were calculated for the optimized geometry using DFT/B3LYP cc-pVQZ basis set by Gaussian 09 Program. The vibrational frequencies were assigned on the basis of potential energy distribution calculation using VEDA 4.0 program package. The Mulliken atomic charge values were calculated. In the Frontier molecular orbitals analysis, the molecular reactivity, kinetic stability, intermolecular charge transfer studies and the calculation of ionization energy, electron affinity, global hardness, chemical potential, electrophilicity index and softness of the molecule were carried out. The nonlinear optical (NLO) activity was studied and the first order hyperpolarizability value was computed, which was 3.48 times greater than the urea. The natural bond orbital analysis was also performed to confirm the NLO activity of the molecule. Hence, the ACTP molecule is a promising candidate for NLO materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Vibrational spectroscopic, structural and nonlinear optical activity studies on 6-aminonicotinamide: A DFT approach

The conformational analysis was carried out for 6-aminonicotinamide (ANA) using potential energy surface scan method and the most stable optimized conformer was predicted. The theoretical vibrational frequencies were calculated for the optimized geometry using DFT/B3LYP cc-pVQZ basis set by Gaussian 09 Program. The vibrational frequencies were assigned on the basis of potential energy distribution calculation using VEDA 4.0 program. The Mulliken atomic charge values were calculated. In the Frontier molecular orbitals analysis, the molecular reactivity, kinetic stability, intermolecular charge transfer studies and the related molecular properties were calculated. The ultraviolet-visible spectrum was simulated for both in the gas phase and liquid phase (ethanol) and the π to π* electronic transition was predicted. The nonlinear optical (NLO) activity was studied by means of the first order hyperpolarizability value, which was 8.61 times greater than the urea and the natural bond orbital analysis was also performed to confirm the NLO activity of the molecule. Hence, the ANA molecule is a promising candidate for the NLO materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Block-structured, equal-workload, multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Total Variation Diminishing)

Abstract. We describe the development of a block-structured, equal-CPU-load (central processing unit), multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Fully Nonlinear Boussinesq Wave Model with Total Variation Diminishing Solver). The new model framework does not interfere with the core solver, and thus the core program, FUNWAVE-TVD, is still a standalone model used for a single grid. The nesting interface manages the time sequencing and two-way nesting processes between the parent grid and child grid with grid refinement in a hierarchical manner. Workload balance in the MPI-based (message passing interface) parallelization is handled by an equal-load scheme. A strategy of shared array allocation is applied for data management that allows for a large number of nested grids without creating additional memory allocations. Four model tests are conducted to verify the nesting algorithm with assessments of model accuracy and the robustness in the application in modeling transoceanic tsunamis and coastal effects.

Choi, Young-Kwang↗
DOE OSTI2022

Robust trajectory-constrained frequency control for microgrids considering model linearization error

Grid supportive modes integrated within inverter-based resources can improve the frequency response of renewable-rich microgrids. The synthesis of grid supportive modes to guarantee frequency trajectory constraints under a predefined disturbance set is challenging but essential. To tackle this challenge, a numerical optimal control (NOC)-based control synthesis methodology is proposed. Without loss of generality, a wind-diesel fed microgrid is studied, where we aim to design grid supportive functions in the wind turbine. In the control design, linearized models are used, and the linearization-induced errors are quantitatively analyzed by reachability and interval arithmetics and represented in the form of interval uncertainties. Then, the NOC problem can be formulated into a robust mixed-integer linear program. The control structure is strategically configured into two levels to realize online deployment. The proposed control is verified on the modified 33-node microgrid with a full-order three-phase nonlinear model in Simulink. In conclusion, the simulation results show the effectiveness of the proposed control paradigm and the necessity of considering linearization-induced uncertainty.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

A theoretical study on 2-amino-5-nitroprydinium trifluoroaceta

The geometry optimization of 2-amino-5-nitroprydinium trifluoroacetate molecule was carried out by using Becke’s three-parameter exchange functional in conjunction with the Lee-Yang-Parr correlation functional (B3LYP) level of density functional theory (DFT) and 6-311++G(d,p) basis set at GAUSSIAN 09 program. The vibration spectrum of the title compound was simulated to predict the presence of functional groups and their vibrational modes. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies were calculated at the same level, and the obtained small energy gap shows that charge transfer occurs in the title compound. The molecular dipole moment, polarizability and hyperpolarizability parameters were determined to evaluate nonlinear optical efficiency of the title compound. Finally, the {sup 13}C and {sup 1}H Nuclear Magnetic Resonance (NMR) chemical shift values were calculated by the application of the gauge independent atomic orbital (GIAO) method. All of the calculations were carried out by using GAUSSIAN 09 program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Nonlinear Methods for Current Limit Constraint Satisfaction in Tokamak Plasma Shape Control

The DIII-D tokamak is capable of supporting a wide variety of plasma equilibria because of its relatively large number of coils and their proximity to the plasma. To support its advanced tokamak mission, the DIII–D experimental program continues to push the envelope of this capability, frequently encountering limits imposed by allowable currents in poloidal shaping coils. Violation of current constraints is presently dealt with by operator adjustment of control targets and gains between plasma discharges. At the same time, demands for more precise and stable control have motivated efforts to develop and install advanced multivariable algorithms for control of plasma shape in DIII-D and other devices. Further, there is currently no way to ensure respect of nonlinear current constraints in a multivariable linear controller design and no practical way to manually tune these fully coupled controllers between discharges after installation. Various linear minimization schemes can be implemented to encourage currents to remain within limits, but adherence to these limits cannot be guaranteed by linear methods alone. In this paper, we describe ongoing efforts to provide methods that guarantee currents will not exceed preset limits, and that simultaneously achieve the best obtainable quality of control subject to current limit constraints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Numerical Simulations of Space Charge Compensation with an Electron Lens

The future high energy physics program at Fermilab requires that the proton complex operate with beam bunch intensities four times larger than is currently handled. At these intensities space charge nonlinear defocussing effects cause unacceptable particle losses especially in the low energy rapid-cycling-synchrotron (RCS) Booster. Focusing electron lens elements may offer a solution by providing partial space charge compensation but there is a need for detailed simulations as this technique has not been demonstrated. We report on high fidelity numerical 6D space charge simulations in a model accelerator lattice with a record high space charge tune shift approaching unity.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Center for Nonlinear Studies (CNLS) Leader Position [Slides]

The Center Director provides scientific leadership and line management of the CNLS while fostering collaborations with scientists throughout the Laboratory. The CNLS Center Director is expected to develop and lead a program to target and create cooperative long-term research programs consistent with the Laboratory’s strategic research objectives, to develop a strong working relationship with the CNLS External Advisory Committee, and to maintain effective working relationships throughout all levels of the Laboratory, government entities, academia and industry. The successful candidate will be expected to maintain an active research program while providing technical vision to nurture and support existing programs of others at the Center.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2016

Scalable Optimization Methods for Distribution Networks With High PV Integration

Here this paper proposes a suite of algorithms to determine the active- and reactive-power setpoints for photovoltaic (PV) inverters in distribution networks. The objective is to optimize the operation of the distribution feeder according to a variety of performance objectives and ensure voltage regulation. In general, these algorithms take a form of the widely studied ac optimal power flow (OPF) problem. For the envisioned application domain, nonlinear power-flow constraints render pertinent OPF problems nonconvex and computationally intensive for large systems. To address these concerns, we formulate a quadratic constrained quadratic program (QCQP) by leveraging a linear approximation of the algebraic power-flow equations. Furthermore, simplification from QCQP to a linearly constrained quadratic program is provided under certain conditions. The merits of the proposed approach are demonstrated with simulation results that utilize realistic PV-generation and load-profile data for illustrative distribution-system test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Emerging Functionality in Transition-Metal Compounds Driven by Spatial Confinement and Broken Symmetry

This research project investigates the emerging functionality in transition-metal-compounds (TMCs) driven by spatial confinement and broken symmetry. It combines advanced growth capabilities with cutting-edge characterization and first principles theory to probe and control the properties of TMC interfaces, including the utilization and development of state-of-the-art atomically resolved electron microscopy and spectroscopy to determine the structure, composition, and bonding at TMC interfaces. The proposed research will focus on four challenging areas: 1) manipulate interfaces to design new material phases such as magnetic metals with unique polar structure (net dipole) to achieve multiple functionality; 2) explore electronic mismatch or screening at interfaces of insulating/poor metal TMCs to produce novel electronic and magnetic properties; 3) elucidate and exploit the role of defects, both point and extended, on the functionality of interfaces; 4) develop advanced electron microscopy/spectroscopy techniques to explore temperature dependent phase transitions and couple these structural tools with new nonlinear optical probes of the electronic structure. The research team aims to close the materials-by-design loop of make, measure, model, and modify. The program promises to enhance our ability to engineer the desired physical properties at interfaces, superlattices (periodic arrays of films of different compounds), and heterostructures of TMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Incorporation of market signals for the optimal design of post combustion carbon capture systems

Recent studies have shown that fossil generators equipped with post-combustion carbon capture (PCC) systems are needed to reduce the cost of deep decarbonization. Such generators need to be flexible and responsive to grid conditions, particularly in a high variable renewable energy (VRE) environment. In this work, we evaluate the net present value (NPV) of retrofitting an existing natural gas combined cycle (NGCC) unit with a flexible PCC system while incorporating market signals from a high VRE grid. We use our industrial partner’s NGCC configuration as representative of existing NGCC units and Svante’s rapid-temperature swing adsorption (TSA) for PCC. Because of its ability to rapidly startup/shutdown and ramp-up/ramp-down, the chosen capture technology is very attractive for load-following operations. For a given set of market signals, we formulate a two-stage stochastic multi-period optimization problem, under the price-taker assumption, to simultaneously optimize the design of the capture system and operation of the entire plant. Rigorous models for the NGCC unit, PCC system, and compression system are developed using commercial process simulators and validated with either plant or vendor data. For computational tractability, we develop surrogate/reduced-order models for use in the optimization problem. The surrogate model for the NGCC plant is constructed by linearizing the rigorous dynamic model at 75% load, while data-driven nonlinear surrogate models for the capture and compression systems are constructed using simulation data from the rigorous models. The optimization problem, formulated as a mixed integer bilinear program, is implemented in the IDAES® integrated platform and solved to global optimality using Gurobi 9.5. Using this formulation, we determine the profitability of retrofitting an existing NGCC unit with the chosen capture system for multiple regions in the U.S. under two scenarios with different carbon prices. Importantly, the results show that the optimal decision strongly depends on the region and on the carbon price, thereby demonstrating the importance of the inclusion of market signals in the design process.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Cluster for Research Excellence (CRE) for Accelerator Physics

The purpose of this Cooperative Research and Development Agreement (CRADA) was to establish a Cluster of Research Excellence (CRE) for Accelerator Physics and broaden the accelerator physics collaboration between Northern Illinois University (NIU) and Fermilab. The CRE was based in both Fermilab’s Accelerator Division and in NIU’s Department of Physics. A research agenda and strategy for the CRE for Accelerator Physics included specific research topics, experiments, and research facility upgrades. The scope of the research program covered two thrusts with subtopics in the Accelerator Science Research Program utilizing the infrastructure of FAST, IOTA, the Fermilab Accelerator Complex, and the Illinois Accelerator Research Center (IARC). Specific areas of proposed exploration included: 1.Tests of novel “nonlinear integrable” systems including beam self-fields; 2.Experimental tests of phase-space diffusion and chaos in intense beams; 3.Investigate technically feasible approaches towards next generation intenseneutrino facilities; 4.Understanding optical control of charged particle beams; 5.Exploitation of bright electron beams from FAST for novel applications; 6. Conduct research into accelerator-based ‘precision science’, such as Muon g-2, EDM, Mu2e, etc.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Simulating non-native cubic interactions on noisy quantum machines

As a milestone for general-purpose computing machines, we demonstrate that quantum processors can be programed to efficiently simulate dynamics that are not native to the hardware. Moreover, on noisy devices without error correction, we show that simulation results are significantly improved when the quantum program is compiled using modular gates instead of a restricted set of standard gates. We demonstrate the general methodology by solving a cubic interaction problem, which appears in nonlinear optics, gauge theories, as well as plasma and fluid dynamics. To encode the non-native Hamiltonian evolution, we decompose the Hilbert space into a direct sum of invariant subspaces in which the nonlinear problem is mapped to a finite-dimensional Hamiltonian simulation problem. Furthermore, in a three-states example, the resultant unitary evolution is realized by a product of approximately 20 standard gates, using which approximately ten simulation steps can be carried out on state-of-the-art quantum hardware before results are corrupted by decoherence. In comparison, the simulation depth is improved by more than an order of magnitude when the unitary evolution is realized as a single cubic gate, which is compiled directly using optimal control. Alternatively, parametric gates may also be compiled by interpolating control pulses. Modular gates thus obtained provide high-fidelity building blocks for quantum Hamiltonian simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Control of Three Degrees-of-Freedom Wave Energy Converters Using Pseudo-Spectral Methods

This study presents a solution to the optimal control problem of a three degrees-of-freedom (3DOF) wave energy converter (WEC). The three modes are the heave, pitch, and surge. The dynamic model is characterized by a coupling between the pitch and surge modes, while the heave is decoupled. The heave, however, excites the pitch motion through nonlinear parametric excitation in the pitch mode. This paper uses Fourier series (FS) as basis functions to approximate the states and the control. A simplified model is first used where the parametric excitation term is neglected and a closed-form solution for the optimal control is developed. For the parametrically excited case, a sequential quadratic programming approach is implemented to solve for the optimal control numerically. Numerical results show that the harvested energy from three modes is greater than three times the harvested energy from the heave mode alone. Moreover, the harvested energy using a control that accounts for the parametric excitation is significantly higher than the energy harvested when neglecting this nonlinear parametric excitation term.

16 TIDAL AND WAVE POWER↗