Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Structural/control system optimization with variable actuator masses

A method is presented to integrate the design space for structural/control system optimization problems in the case of linear state feedback control. Nonstructural lumped masses and control system design variables as well as structural sizing variables are all treated equally as independent design variables in the optimization process. Structural and control design variable linking schemes are used in order to avoid a prohibitively large increase in the total number of independent design variables. When actuator masses are treated as nonstructural lumped mass design variables, special consideration is given to the relation between the transient peak responses and the required actuator masses which is formulated as a behavior constraint form. The original nonlinear mathematical programming problem based on a finite element formulation and linear state feedback is replaced by a sequence of explicit approximate problems exploiting various approximation concepts such as design variable linkings, temporary constraint deletion and first order Taylor series expansion of nonlinear behavior constraints in terms of intermediate design variables. Examples which involve a variety of dynamic behavior constraints (including constraints on closed-loop eigenvalues, peak transient displacements, peak actuator forces, and relations between the peak responses and the actuator masses) are effectively solved by using the method presented.

Jin, Ik M.↗
NASA NTRS1983

Critical joints in large composite aircraft structure

A program was conducted at Douglas Aircraft Company to develop the technology for critical structural joints of composite wing structure that meets design requirements for a 1990 commercial transport aircraft. The prime objective of the program was to demonstrate the ability to reliably predict the strength of large bolted composite joints. Ancillary testing of 180 specimens generated data on strength and load-deflection characteristics which provided input to the joint analysis. Load-sharing between fasteners in multirow bolted joints was computed by the nonlinear analysis program A4EJ. This program was used to predict strengths of 20 additional large subcomponents representing strips from a wing root chordwise splice. In most cases, the predictions were accurate to within a few percent of the test results. In some cases, the observed mode of failure was different than anticipated. The highlight of the subcomponent testing was the consistent ability to achieve gross-section failure strains close to 0.005. That represents a considerable improvement over the state of the art.

Nelson, W. D.↗
DOE OSTI2013

Development of a Generalized Methodology for Soil-Structure Interaction Analysis Using Nonlinear Time-Domain TechniquesNEAMS Program, DOE Office of Nuclear Energy (NE-41)

A generalized time-domain method for Soil-Structure Interaction Analysis is developed, based upon an extension of the Bielak Method. The methodology is combined with the use of a simple hysteretic soil model based upon the Ramberg-Osgood formulation and applied to a notional Small Modular Reactor. These benchmark results compare well with those obtained by using the industry-standard frequency-domain code SASSI. The methodology provides a path forward for investigation of other sources of nonlinearity, including those associated with the use of more physically-realistic material models incorporating pore-pressure effects, gap opening/closing, the effect of nonlinear structural elements, and 3D seismic inputs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1989

Application and implementation of transient algorithms in computer programs

A brief introduction is given to the nonlinear finite element programs developed at Lawrence Livermore National Laboratory. The four programs are DYNA3D and DYNA2D, which are explicit hydrocodes, and NIKE3D and NIKE2D, which are implicit programs. The main emphasis is on DYNA3D with asides about the other programs. During the past year several new features were added to DYNA3D, and major improvements were made in the computational efficiency of the shell and beam elements. Most of these new features and improvements will eventually make their way into the other programs. The emphasis in the computational mechanics effort was always, and continues to be, efficiency. To get the most out of the supercomputers, all Crays, the programs were vectorized where possible. Several of the more interesting capabilities of DYNA3D will be described and the impact on efficiency will be discussed. Some of the recent work on NIKE3D and NIKE2D will also be presented. In the belief that a single example is worth a thousand equations, the theory is skipped entirely and the examples presented.

Benson, David J.↗
DOE OSTI2018

Optimal Control for Battery Storage Using Nonlinear Models

Battery storage system (BSS) becomes increasingly popular for grid applications due to the advance in battery and power electronics technologies as well as the growing need of flexibility and reserve from power systems with rapidly developed renewable generation. Successful assessment and deployment of BSS requires optimizing its operation and therefore maximizing the potential benefits. Although there exist studies that develop optimal control for evaluation or operational scheduling of BSS, modeling of charging/discharging operation and the corresponding impacts on state-of-charge (SOC) is over simplified in these methods. This paper develops an optimal control for BSS evaluation and operational scheduling using a general nonlinear model that expresses the change of SOC as a function of charging/discharging power and SOC level. The proposed method is compared with a typical existing method through a real-world energy storage evaluation project to show the significance of the proposed method.

battery, energy storage system, dynamic programmin↗
DOE OSTI2025

Efficient proximal subproblem solvers for a nonsmooth trust-region method

In [R. J. Baraldi and D. P. Kouri, Mathematical Programming, (2022), pp. 1-40], we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex and nonsmooth convex function. The principle expense of this method is in computing a trial iterate that satisfies the so-called fraction of Cauchy decrease condition—a bound that ensures the trial iterate produces sufficient decrease of the subproblem model. In this paper, we expound on various proximal trust-region subproblem solvers that generalize traditional trust-region methods for smooth unconstrained and convex-constrained problems. We introduce a simplified spectral proximal gradient solver, a truncated nonlinear conjugate gradient solver, and a dogleg method. Finally, we compare algorithm performance on examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

On the Value Function of Weakly Coercive Problems in Nonlinear Stochastic Control

In this paper we investigate via a dynamic programming approach some nonlinear stochastic control problems where the control set is unbounded and a classical coercivity hypothesis is replaced by some weaker assumptions. We prove that these problems can be approximated by finite fuel problems; show the continuity of the relative value functions and characterize them as unique viscosity solutions of a quasi-variational inequality with suitable boundary conditions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Development of computer program NAS3D using Vector processing for geometric nonlinear analysis of structures

An algorithm for vectorized computation of stiffness matrices of an 8 noded isoparametric hexahedron element for geometric nonlinear analysis was developed. This was used in conjunction with the earlier 2-D program GAMNAS to develop the new program NAS3D for geometric nonlinear analysis. A conventional, modified Newton-Raphson process is used for the nonlinear analysis. New schemes for the computation of stiffness and strain energy release rates is presented. The organization the program is explained and some results on four sample problems are given. The study of CPU times showed that savings by a factor of 11 to 13 were achieved when vectorized computation was used for the stiffness instead of the conventional scalar one. Finally, the scheme of inputting data is explained.

Mangalgiri, P. D.↗
DOE OSTI2022

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗
DOE OSTI2025

Multistage economic MPC for systems with a cyclic steady state: A gas network case study

Multistage model predictive control (MPC) provides a robust control strategy for dynamic systems with uncertainties and a setpoint tracking objective. Moreover, extending MPC to minimize an economic cost instead of tracking a pre-calculated optimal setpoint improves controller performance. This paper presents a novel multistage economic nonlinear model predictive control (E-NMPC) framework for dynamic systems operating under uncertainty, with specific application to natural gas transmission networks. A key innovation lies in the integration of cyclic steady-state (CSS) constraints within the multistage MPC formulation, enabling the controller to manage periodic operating conditions commonly observed in energy systems. A Lyapunov-based descent condition is enforced to ensure robust stability of the controller. The multistage economic MPC framework is validated on two gas pipeline case studies, where it successfully minimizes net energy consumption, respects operational constraints under uncertain demand profiles, and guides the network to optimal cyclic operation. The Lyapunov function remains bounded in both case studies, validating the robust stability of multistage E-NMPC.

03 NATURAL GAS↗
DOE OSTI2016

A superstructure-based framework for bio-separation network synthesis

Modern biotechnologies enable the production of chemicals using engineered microorganisms. However, the cost of downstream recovery and purification steps is high, which means that the feasibility of bio-based chemicals production depends heavily on the synthesis of cost-effective separation networks. To this end, we develop a superstructure-based framework for bio-separation network synthesis. Based on general separation principles and insights obtained from industrial processes for specific products, we first identify four separation stages: cell treatment, product phase isolation, concentration and purification, and refinement. For each stage, we systematically implement a set of connectivity rules to develop stage-superstructures, all of which are then integrated to generate a general superstructure that accounts for all types of chemicals that can be produced using microorganisms. Here we further develop a superstructure reduction method to solve specific instances, based on product attributes, technology availability, case-specific considerations, and final product stream specifications. A general optimization model, including short-cut models for all technologies, is formulated. The proposed framework enables preliminary synthesis and analysis of bio-separation networks, and thus estimation of separation costs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

A generalized superstructure-based framework for process synthesis

We propose a general framework for the formulation of superstructure-based optimization models for holistic process synthesis. First, we redefine the fundamental problems of reactor, separation, and heat exchanger network synthesis, by presenting generalized problem statements, to make them amenable to seamless integration with each other. Second, we describe the general forms of models that can be developed to address these generalized problems and identify some key characteristics. Notably, for each system, we identify internal variables used only within the system, cost variables used in the objective function, and, importantly, coupling variables for the coupling between systems. Third, we outline some literature models that can be used to address the generalized problems and present new models to couple the three systems. Finally, we show how the individual components (systems and coupling models) can be integrated to formulate a single simultaneous reactor, separation, and heat exchanger network synthesis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Integration of graphical approaches into optimization-based design of multistage liquid extraction

We propose two optimization models for designing two liquid extraction systems: simple multistage liquid extractors and extractors with extract reflux. Both models are motivated by the concepts of the modified McCabe-Thiele graphical method for multistage extractor design. The operating and equilibrium curves in the McCabe-Thiele method are represented by material balances and piece-wise linearized thermodynamics properties. The use of piece-wise approximations improves computational tractability of both optimization models. In addition, we consider some extensions such as dilute systems, insoluble solvents, and non-ideal stages. In conclusion, the applicability of the proposed models is demonstrated with four illustrative examples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Inexact convex relaxations for AC optimal power flow: Towards AC feasibility

Convex relaxations of AC optimal power flow (AC-OPF) problems have attracted significant interest as in several instances they provably yield the global optimum to the original non-convex problem. If, however, the relaxation is inexact, the obtained solution is not AC-feasible. The quality of the obtained solution is essential for several practical applications of AC-OPF, but detailed analyses are lacking in existing literature. Here, this paper aims to cover this gap. We provide an in-depth investigation of the solution characteristics when convex relaxations are inexact, we assess the most promising AC feasibility recovery methods for large-scale systems, and we propose two new metrics that lead to a better understanding of the quality of the identified solutions. We perform a comprehensive assessment on 96 different test cases, ranging from 14 to 3120 buses, and we show the following: (i) Despite an optimality gap of less than 1%, several test cases still exhibit substantial distances to both AC feasibility and local optimality and the newly proposed metrics characterize these deviations. (ii) Penalization methods fail to recover an AC-feasible solution in 15 out of 45 test cases. (iii) The computational benefits of warm-starting non-convex solvers have significant variation, but a computational speedup exists in over 75% of cases.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Estimation of the Maximum Theoretical Productivity of Fed-Batch Bioreactors

A key step towards the development of an integrated biorefinery is the screening of economically viable processes, which depends sharply on the yields and productivities that can be achieved by an engineered microorganism. In this study, we extend an earlier method which used dynamic optimization to find the maximum theoretical productivity of batch cultures to explicitly include fed-batch bioreactors. In addition to optimizing the intracellular distribution of metabolites between cell growth and product formation, we calculate the optimal control trajectory of feed rate versus time. We further analyze how sensitive the productivity is to substrate uptake and growth parameters.

09 BIOMASS FUELS↗