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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Materials constitutive models for nonlinear analysis of thermally cycled structures

Effects of inelastic materials models on computed stress-strain solutions for thermally loaded structures were studied by performing nonlinear (elastoplastic creep) and elastic structural analyses on a prismatic, double edge wedge specimen of IN 100 alloy that was subjected to thermal cycling in fluidized beds. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic kinematic, and combined plus transient creep) were exercised for the problem by using the MARC nonlinear, finite element computer program. Maximum total strain ranges computed from the elastic and nonlinear analyses agreed within 5 percent. Mean cyclic stresses, inelastic strain ranges, and inelastic work were significantly affected by the choice of inelastic constitutive model. The computing time per cycle for the nonlinear analyses was more than five times that required for the elastic analysis.

Kaufman, A.↗
DOE OSTI2024

A mixed integer linear programming approach for the design of chemical process families

Tackling climate change goals requires widespread deployment of process technology variants across many decentralized sites with different geographical, environmental, and operational requirements. Conventional engineering approaches focus on unique designs for each installation (process variant), while missing opportunities for manufacturing standardization. Here, instead we seek to optimize a process platform of common unit designs while simultaneously designing an entire family of process variants that make use of that platform. This reduces engineering effort, deployment timelines, and manufacturing costs. We propose a nonlinear generalized disjunctive programming formulation and convert this to an efficient mixed-integer linear programming (MILP) formulation through discretization of the design space. We formulate our optimization in Pyomo with costing from IDAES, and we demonstrate the computational performance and solution quality on a water treatment desalination system from the PARETO framework and a carbon capture system built in Aspen Plus as part of CCSI2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

High Temperature Composite Analyzer (HITCAN) Programmer's Manual. Version 1.0

This manual describes the organization and flow of data and analysis modules in the computer code, HITCAN (High Temperature Composite ANalyzer). HITCAN is a general purpose computer program for predicting nonlinear global structural and local stress-strain response of arbitrarily oriented, multilayered high temperature metal matrix composite structures. This manual describes the architecture of the HITCAN code, followed by the listing of subroutines and calling tree, data storage scheme, file system, and a dictionary of code terminology. The primary intention of the manual is to familiarize the user with some of the computer program related issues so as to facilitate maintenance/modification/updates of the HITCAN computer code.

Lackney, J. J.↗
NASA NTRS1993

High Temperature Composite Analyzer (HITCAN) theoretical manual, version 1.0

This manual outlines some of the theoretical aspects embedded in the computer code, HITCAN (High Temperature Composite ANalyzer). HITCAN is a general purpose computer program for predicting nonlinear global structural and local stress-strain response of arbitrarily oriented, multilayered high temperature metal matrix composite structures. This code combines composite mechanics and laminate theory with an internal data base of the constituents (matrix, fiber and interface) material properties. The thermal and mechanical properties of the constituents are considered to be nonlinearly dependent on several parameters including temperature, stress and stress rate. The computational procedure for the analysis of the composite structure uses the finite element method. HITCAN consists of three modules: COBSTRAN, METCAN and MHOST. COBSTRAN generates the geometry (pre-processor) and defines the layup of the different plies. METCAN computes the material behavior of the composite and of the constituents. Finally, MHOST is a finite element program based on the mixed iterative solution technique. MHOST has a library for 2D and 3D isoparametric elements. HITCAN is written in FORTRAN 77 computer language and has been configured and executed on the NASA Lewis Research Center CRAY XMP and YMP computers.

Lackney, J. J.↗
NASA NTRS1972

Method for nonlinear exponential regression analysis

Two computer programs developed according to two general types of exponential models for conducting nonlinear exponential regression analysis are described. Least squares procedure is used in which the nonlinear problem is linearized by expanding in a Taylor series. Program is written in FORTRAN 5 for the Univac 1108 computer.

Junkin, B. G.↗
NASA NTRS1986

A general numerical analysis program for the superconducting quasiparticle mixer

A user-oriented computer program SISCAP (SIS Computer Analysis Program) for analyzing SIS mixers is described. The program allows arbitrary impedance terminations to be specified at all LO harmonics and sideband frequencies. It is therefore able to treat a much more general class of SIS mixers than the widely used three-frequency analysis, for which the harmonics are assumed to be short-circuited. An additional program, GETCHI, provides the necessary input data to program SISCAP. The SISCAP program performs a nonlinear analysis to determine the SIS junction voltage waveform produced by the local oscillator. The quantum theory of mixing is used in its most general form, treating the large signal properties of the mixer in the time domain. A small signal linear analysis is then used to find the conversion loss and port impedances. The noise analysis includes thermal noise from the termination resistances and shot noise from the periodic LO current. Quantum noise is not considered. Many aspects of the program have been adequately verified and found accurate.

Hicks, R. G.↗
NASA NTRS1993

High temperature composite analyzer (HITCAN) user's manual, version 1.0

This manual describes 'how-to-use' the computer code, HITCAN (HIgh Temperature Composite ANalyzer). HITCAN is a general purpose computer program for predicting nonlinear global structural and local stress-strain response of arbitrarily oriented, multilayered high temperature metal matrix composite structures. This code combines composite mechanics and laminate theory with an internal data base for material properties of the constituents (matrix, fiber and interphase). The thermo-mechanical properties of the constituents are considered to be nonlinearly dependent on several parameters including temperature, stress and stress rate. The computation procedure for the analysis of the composite structures uses the finite element method. HITCAN is written in FORTRAN 77 computer language and at present has been configured and executed on the NASA Lewis Research Center CRAY XMP and YMP computers. This manual describes HlTCAN's capabilities and limitations followed by input/execution/output descriptions and example problems. The input is described in detail including (1) geometry modeling, (2) types of finite elements, (3) types of analysis, (4) material data, (5) types of loading, (6) boundary conditions, (7) output control, (8) program options, and (9) data bank.

Lackney, J. J.↗
NASA NTRS1971

The NASTRAN computer program for structural analysis.

Discussion of a large NASA-sponsored digital computer program for static and dynamic analysis of structures with linear or nonlinear material properties. The program fits to perform buckling analysis, vibration mode analysis of conservative or damped systems, frequency response analysis, and transient response analysis on the IBM 360, UNIVAC 1108, and CDC 6600 computers. It also provides automatic generation of several types of loads such as gravity loads and thermal strain loads, plots of deformed and nondeformed structures, and the choice of alternate analysis methods such as direct or modal formulations of dynamic problems. The program is intended primarily for large problems involving hundreds of thousands of degrees of freedom. A flexible Executive System is provided to facilitate the addition of new functional capabilities to the program.

Macneal, R. H.↗
DOE OSTI2013

MEASURING THE EVOLUTIONARY RATE OF COOLING OF ZZ Ceti

We have finally measured the evolutionary rate of cooling of the pulsating hydrogen atmosphere (DA) white dwarf ZZ Ceti (Ross 548), as reflected by the drift rate of the 213.13260694 s period. Using 41 yr of time-series photometry from 1970 November to 2012 January, we determine the rate of change of this period with time to be dP/dt = (5.2 {+-} 1.4) Multiplication-Sign 10{sup -15} s s{sup -1} employing the O - C method and (5.45 {+-} 0.79) Multiplication-Sign 10{sup -15} s s{sup -1} using a direct nonlinear least squares fit to the entire lightcurve. We adopt the dP/dt obtained from the nonlinear least squares program as our final determination, but augment the corresponding uncertainty to a more realistic value, ultimately arriving at the measurement of dP/dt = (5.5 {+-} 1.0) Multiplication-Sign 10{sup -15} s s{sup -1}. After correcting for proper motion, the evolutionary rate of cooling of ZZ Ceti is computed to be (3.3 {+-} 1.1) Multiplication-Sign 10{sup -15} s s{sup -1}. This value is consistent within uncertainties with the measurement of (4.19 {+-} 0.73) Multiplication-Sign 10{sup -15} s s{sup -1} for another similar pulsating DA white dwarf, G 117-B15A. Measuring the cooling rate of ZZ Ceti helps us refine our stellar structure and evolutionary models, as cooling depends mainly on the core composition and stellar mass. Calibrating white dwarf cooling curves with this measurement will reduce the theoretical uncertainties involved in white dwarf cosmochronometry. Should the 213.13 s period be trapped in the hydrogen envelope, then our determination of its drift rate compared to the expected evolutionary rate suggests an additional source of stellar cooling. Attributing the excess cooling to the emission of axions imposes a constraint on the mass of the hypothetical axion particle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Minimax Program Terminal Control in Two-Level Hierarchic Nonlinear Discrete Dynamical System

We consider a discrete dynamical system consisting of a number of controllable objects whose dynamics is governed by vector nonlinear discrete recurrent relations with control parameters and perturbations (noises). In such system, two levels of making management decisions are distinguished: the dominating (first) level and the subordinate (second) level, which have different criteria of functioning. These levels are integrated by certain information and control relations. We consider the problem of two-level minimax program terminal control in a discrete dynamical system and propose its mathematical statement and a general scheme of its solution.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1967

[Elliptic Differential Equations]

Reflection laws of systems of elliptic equations, on-line programming to numerically solve nonlinear integral equations, and reduction of electrostriction problem

ON-LINE PROGRAMMING↗
DOE OSTI2016

WE-AB-207B-11: Optimizing Tumor Control Probability in Radiation Therapy Treatment - Application to HDR Cervical Cancer

Purpose: The ultimate goal of radiotherapy treatment planning is to find a treatment that will yield a high tumor-control-probability(TCP) with an acceptable normal-tissue-complication probability(NTCP). Yet most treatment planning today is not based upon optimization of TCPs and NTCPs, but rather upon meeting physical dose and volume constraints defined by the planner. We design treatment plans that optimize TCP directly and contrast them with the clinical dose-based plans. PET image is incorporated to evaluate gain in TCP for dose escalation. Methods: We build a nonlinear mixed integer programming optimization model that maximizes TCP directly while satisfying the dose requirements on the targeted organ and healthy tissues. The solution strategy first fits the TCP function with a piecewise-linear approximation, then solves the problem that maximizes the piecewise linear approximation of TCP, and finally performs a local neighborhood search to improve the TCP value. To gauge the feasibility, characteristics, and potential benefit of PET-image guided dose escalation, initial validation consists of fifteen cervical cancer HDR patient cases. These patients have all received prior 45Gy of external radiation dose. For both escalated strategies, we consider 35Gy PTV-dose, and two variations (37Gy-boost to BTV vs 40Gy-boost) to PET-image-pockets. Results: TCP for standard clinical plans range from 59.4% - 63.6%. TCP for dose-based PET-guided escalated-dose-plan ranges from 63.8%–98.6% for all patients; whereas TCP-optimized plans achieves over 91% for all patients. There is marginal difference in TCP among those with 37Gy-boosted vs 40Gy-boosted. There is no increase in rectum and bladder dose among all plans. Conclusion: Optimizing TCP directly results in highly conformed treatment plans. The TCP-optimized plan is individualized based on the biological PET-image of the patients. The TCP-optimization framework is generalizable and has been applied successfully to other external-beam delivery modalities. A clinical trial is on-going to gauge the clinical significance. Partially supported by the National Science Foundation.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1976

Progressive failure of notched composite laminates using finite elements

An incremental, nonlinear, finite element program capable of tracing the damage that occurs in a notched laminated composite, before ultimate failure takes place, is described. Computer predictions for the ultimate failure load and failure modes are compared to experimental results for two laminates containing holes and loaded in uniaxial tension. The comparison is quite favorable.

Nuismer, R. J.↗
NASA NTRS1976

Cyclic structural analysis of air-cooled gas turbine blades and vanes

The creep fatigue behavior of a fully impingement cooled blade for four cyclic cases was analyzed by using the Elas 55, finite element, nonlinear structural computer program. Expected cyclic lives were calculated by using the method of strainrange partitioning for reversed inelastic strains and time fractions for ratcheted tensile creep strains. Strainrange partitioning was also applied to previous results from a one dimensional cyclic analysis of a film impingement cooled vane. The analyses indicated that strainrange partitioning is more applicable to a constrained airfoil such as the film impingement cooled vane than to the relatively unconstrained fully impingement cooled airfoil.

Kaufman, A.↗
NASA NTRS1976

Cyclic structural analyses of air-cooled gas turbine blades and vanes

The creep-fatigue behavior of a fully impingement-cooled blade for four cyclic cases was analyzed by using the Elas 55, finite-element, nonlinear structural computer program. Expected cyclic lives were calculated by using the method of Strainrange Partitioning for reversed inelastic strains and time fractions for ratcheted tensile creep strains. Strainrange Partitioning was also applied to previous results from a one-dimensional cyclic analysis of a film-impingement-cooled vane. The analyses indicated that Strainrange Partitioning is more applicable to a constrained airfoil such as the film-impingement-cooled vane than to the relatively unconstrained fully impingement-cooled airfoil.

Kaufman, A.↗
NASA NTRS1976

Alternative methods for the design of jet engine control systems

Various alternatives to linear quadratic design methods for jet engine control systems are discussed. The main alternatives are classified into two broad categories: nonlinear global mathematical programming methods and linear local multivariable frequency domain methods. Specific studies within these categories include model reduction, the eigenvalue locus method, the inverse Nyquist method, polynomial design, dynamic programming, and conjugate gradient approaches.

Sain, M. K.↗
NASA NTRS1977

Mathematical and computational studies of equilibrium capillary free surfaces

The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.

Albright, N.↗