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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗
NASA NTRS1995

Prestraining and Its Influence on Subsequent Fatigue Life

An experimental program was conducted to study the damaging effects of tensile and compressive prestrains on the fatigue life of nickel-base, Inconel 718 superalloy at room temperature. To establish baseline fatigue behavior, virgin specimens with a solid uniform gage section were fatigued to failure under fully-reversed strain-control. Additional specimens were prestrained to 2 percent, 5 percent, and 10 percent (engineering strains) in the tensile direction and to 2 percent (engineering strain) in the compressive direction under stroke-control, and were subsequently fatigued to failure under fully-reversed strain-control. Experimental results are compared with estimates of remaining fatigue lives (after prestraining) using three life prediction approaches: (1) the Linear Damage Rule; (2) the Linear Strain and Life Fraction Rule; and (3) the nonlinear Damage Curve Approach. The Smith-Watson-Topper parameter was used to estimate fatigue lives in the presence of mean stresses. Among the cumulative damage rules investigated, best remaining fatigue life predictions were obtained with the nonlinear Damage Curve Approach.

Halford, Gary R.↗
DOE OSTI2025

Intraspecific Diversity in Thermal Performance Determines Phytoplankton Ecological Niche

ABSTRACT Temperature has a primary influence on phytoplankton physiology and ecology. We grew 12 strains of Gephyrocapsa huxleyi isolated from different‐temperature regions for ~45 generations (2 months) and characterised acclimated thermal response curves across a temperature range. Even with similar temperature optima and overlapping cell size, strain growth rates varied between 0.45 and 1 day −1 . Thermal niche widths varied from 16.7°C to 24.8°C, suggesting that strains use distinct thermal response mechanisms. We investigated the implications of this thermal intraspecific diversity using an ocean ecosystem simulation resolving phytoplankton thermal phenotypes. Model analogues of thermal ‘generalists’ and ‘specialists’ resulted in a distinctive global biogeography of thermal niche widths with a nonlinear latitudinal pattern. We leveraged model output to predict ranges of the 12 lab‐reared strains and demonstrated how this approach could broadly refine geographic range predictions. Our combination of observations and modelled biogeography highlights the capacity of diverse groups to survive temperature shifts.

Krinos, Arianna I. [Department of Biology Woods Ho↗
NASA NTRS1986

Multiple grid problems on concurrent-processing computers

Three computer codes were studied which make use of concurrent processing computer architectures in computational fluid dynamics (CFD). The three parallel codes were tested on a two processor multiple-instruction/multiple-data (MIMD) facility at NASA Ames Research Center, and are suggested for efficient parallel computations. The first code is a well-known program which makes use of the Beam and Warming, implicit, approximate factored algorithm. This study demonstrates the parallelism found in a well-known scheme and it achieved speedups exceeding 1.9 on the two processor MIMD test facility. The second code studied made use of an embedded grid scheme which is used to solve problems having complex geometries. The particular application for this study considered an airfoil/flap geometry in an incompressible flow. The scheme eliminates some of the inherent difficulties found in adapting approximate factorization techniques onto MIMD machines and allows the use of chaotic relaxation and asynchronous iteration techniques. The third code studied is an application of overset grids to a supersonic blunt body problem. The code addresses the difficulties encountered when using embedded grids on a compressible, and therefore nonlinear, problem. The complex numerical boundary system associated with overset grids is discussed and several boundary schemes are suggested. A boundary scheme based on the method of characteristics achieved the best results.

Eberhardt, D. S.↗
DOE OSTI2017

Coupled Physics Simulation of Expansive Reactions in Concrete with the Grizzly Code

The Grizzly code is being developed under the US Department of Energy’s Light Water Sustainability program as a tool to model aging mechanisms and their effects on the integrity of critical nuclear power plant components. An important application for Grizzly is the modeling of aging in concrete structures, which can be due to a number of mechanisms. Initial focus in this area has been on modeling expansive reactions due to alkali-silica reactions or radiation-induced volumetric expansion. Grizzly is an inherently multiphysics modeling platform that naturally permits including the effects of multiple coupled physics in a simulation. Models have been developed for transport of heat and moisture in concrete, and these have been coupled and used as inputs to models for expansive reactions. This paper summarizes this capability as well as work currently underway to couple it with mechanical deformation models that account for nonlinear constitutive behavior.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Generating Mixed Patterns of Residential Segregation: An Evolutionary Approach

The Schelling model of residential segregation has demonstrated that even the slightest preference for neighbors of the same race can be amplified into community-wide segregation. However, these models are unable to simulate mixed, coexisting patterns of segregation and integration, which have been seen to exist in cities. Using evolutionary model discovery we demonstrate how including social factors beyond racial bias when modeling relocation behavior enables the emergence of strongly mixed patterns. Our results indicate that the emergence of mixed patterns is better explained by multiple factors influencing the decision to relocate; the most important being the interaction of nonlinear, rapidly diminishing racial bias with a recent, historical tendency to move. Additionally, preference for less isolated neighborhoods or preference for neighborhoods with longer residing neighbors may produce weaker mixed patterns. Finally, this work highlights the importance of exploring the influence of multiple hypothesized factors of decision making, and their interactions, within agent rules, when studying emergent outcomes generated by agent-based models of complex social systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Application of symbolic and algebraic manipulation software in solving applied mechanics problems

As its name implies, symbolic and algebraic manipulation is an operational tool which not only can retain symbols throughout computations but also can express results in terms of symbols. This report starts with a history of symbolic and algebraic manipulators and a review of the literatures. With the help of selected examples, the capabilities of symbolic and algebraic manipulators are demonstrated. These applications to problems of applied mechanics are then presented. They are the application of automatic formulation to applied mechanics problems, application to a materially nonlinear problem (rigid-plastic ring compression) by finite element method (FEM) and application to plate problems by FEM. The advantages and difficulties, contributions, education, and perspectives of symbolic and algebraic manipulation are discussed. It is well known that there exist some fundamental difficulties in symbolic and algebraic manipulation, such as internal swelling and mathematical limitation. A remedy for these difficulties is proposed, and the three applications mentioned are solved successfully. For example, the closed from solution of stiffness matrix of four-node isoparametrical quadrilateral element for 2-D elasticity problem was not available before. Due to the work presented, the automatic construction of it becomes feasible. In addition, a new advantage of the application of symbolic and algebraic manipulation found is believed to be crucial in improving the efficiency of program execution in the future. This will substantially shorten the response time of a system. It is very significant for certain systems, such as missile and high speed aircraft systems, in which time plays an important role.

Tsai, Wen-Lang↗
NASA NTRS1991

Soluble carbohydrate allocation to roots, photosynthetic rate of leaves, and nitrate assimilation as affected by nitrogen stress and irradiance

Upon resupply of exogenous nitrogen to nitrogen-stressed plants, uptake rate of nitrogen is enhanced relative to nonstressed plants. Absorption of nitrogen presumably is dependent on availability of carbohydrates in the roots. A buildup in soluble carbohydrates thus should occur in roots of nitrogen-stressed plants, and upon resupply of exogenous nitrogen the increased uptake rate should be accompanied by a rapid decline in carbohydrates to prestress levels. To evaluate this relationship, three sets of tobacco plants growing in a complete hydroponic solution containing 1.0 mM NO3- were either continued in the complete solution for 21 d, transferred to a minus-nitrogen solution for 21 d, or transferred to a minus-nitrogen solution for 8-9 d and then returned to the 1.0 mM NO3- solution. These nitrogen treatments were imposed upon plants growing at photosynthetic photon flux densities of 700 and 350 micromoles m-2 s-1. Soluble carbohydrate levels in roots increased during onset of nitrogen stress to levels that were fourfold greater than in roots of non-stressed plants. Following resupply of external nitrogen, a rapid resumption of nitrogen uptake was accompanied by a decline in soluble carbohydrates in roots to levels characteristic of nonstressed plants. This pattern of soluble carbohydrate levels in roots during onset of and recovery from nitrogen stress occurred at both irradiance levels. The response of net photosynthetic rate to nitrogen stress could be expressed as a nonlinear function of concentration of reduced nitrogen in leaves. The net photosynthetic rate at a given concentration of reduced nitrogen, however, averaged 10% less at the lower than at the higher irradiance. The decline in net photosynthetic rate per unit of reduced nitrogen in leaves at the lower irradiance was accompanied by an increase in the nitrate fraction of total nitrogen in leaves from 20% at the higher irradiance to 38% at the lower irradiance.

NASA Discipline Number 61-10↗
DOE OSTI2014

Statistical characteristics of cloud variability. Part 2: Implication for parameterizations of microphysical and radiative transfer processes in climate models

The effects of subgrid cloud variability on grid-average microphysical rates and radiative fluxes are examined by use of long-term retrieval products at the Tropical West Pacific, Southern Great Plains, and North Slope of Alaska sites of the Department of Energy's Atmospheric Radiation Measurement program. Four commonly used distribution functions, the truncated Gaussian, Gamma, lognormal, and Weibull distributions, are constrained to have the same mean and standard deviation as observed cloud liquid water content. The probability density functions are then used to upscale relevant physical processes to obtain grid-average process rates. It is found that the truncated Gaussian representation results in up to 30% mean bias in autoconversion rate, whereas the mean bias for the lognormal representation is about 10%. The Gamma and Weibull distribution function performs the best for the grid-average autoconversion rate with the mean relative bias less than 5%. For radiative fluxes, the lognormal and truncated Gaussian representations perform better than the Gamma and Weibull representations. The results show that the optimal choice of subgrid cloud distribution function depends on the nonlinearity of the process of interest, and thus, there is no single distribution function that works best for all parameterizations. Examination of the scale (window size) dependence of the mean bias indicates that the bias in grid-average process rates monotonically increases with increasing window sizes, suggesting the increasing importance of subgrid variability with increasing grid sizes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Nematic Colloidal Micro‐Robots as Physically Intelligent Systems

Abstract Physically intelligent micro‐robotic systems exploit information embedded in micro‐robots, their colloidal cargo, and their milieu to interact, assemble, and form functional structures. Nonlinear anisotropic fluids such as nematic liquid crystals (NLCs) provide untapped opportunities to embed interactions via their topological defects, complex elastic responses, and ability to dramatically restructure in dynamic settings. Here a four‐armed ferromagnetic micro‐robot is designed and fabricated to embed and dynamically reconfigure information in the nematic director field, generating a suite of physical interactions for cargo manipulation. The micro‐robot shape and surface chemistry are designed to generate a nemato‐elastic energy landscape in the domain that defines multiple modes of emergent, bottom‐up interactions with passive colloids. Micro‐robot rotation expands the ability to sculpt interactions; the energy landscape around a rotating micro‐robot is dynamically reconfigured by complex far‐from‐equilibrium dynamics of the micro‐robot's companion topological defect. These defect dynamics allow transient information to be programmed into the domain and exploited. Robust micro‐robotic manipulation strategies are demonstrated that exploit these diverse modes of nemato‐elastic interaction to achieve cargo docking, transport, release, and assembly of complex reconfigurable structures at multi‐stable sites. Such structures are of great interest to future developments of LC‐based advanced optical device and micro‐manufacturing in anisotropic environments.

Yao, Tianyi↗
NASA NTRS1983

NASA/RAE collaboration on nonlinear control using the F-8C digital fly-by-wire aircraft

A cooperative advanced digital research experiment (CADRE) was established by the National Aeronautics and Space Administration (NASA) and the Royal Aircraft Establishment (RAE), in which nonlinear control algorithms developed by the RAE were tested on the F-8C digital fly-by-wire (DFBW) aircraft based at the Dryden Flight Research Facility. In the initial phase of the collaboration, some variable-gain algorithms, referred to collectively as variable integral control to optimize response (VICTOR) algorithms, were fight tested. With VICTOR, various measures available within the control system are used to vary gains and time-constants within the closed loop and thereby enhance the control capability of the system, while reducing the adverse effects of sensor noise on the control surfaces. A review of design procedures for VICTOR and results of preliminary flight tests are presented. the F-8C aircraft is operated in the remotely augmented vehicle (RAV) mode, with the control laws implemented as FORTRAN programs on a ground-based computer. Pilot commands and sensor information are telemetered to the ground, where the data are processed to form surface commands which are then telemetered to the ground, where the data are processed to form surface commands which are then telemetered back to the aircraft. The RAV mode represents a single-string (simplex) system and is therefore vulnerable to a hardover since comparison monitoring is not possible. Hence, extensive error checking is conducted on both the ground and airborne computers to prevent the development of potentially hazardous situations. Experience with the RAV monitoring and validation procedures is described.

Butler, G. F.↗
NASA NTRS1983

Development of monofilar rotor hub vibration absorber

A design and ground test program was conducted to study the performance of the monofilar absorber for vibration reduction on a four-bladed helicopter. A monofilar is a centrifugal tuned two degree-of-freedom rotor hub absorber that provides force attenuation at two frequencies using the same dynamic mass. Linear and non-linear analyses of the coupled monofilar/airframe system were developed to study tuning and attenuation characteristics. Based on the analysis, a design was fabricated and impact bench tests verified the calculated non-rotating natural frequencies and mode shapes. Performance characteristics were measured using a rotating absorber test facility. These tests showed significant attenuation of fixed-system 4P hub motions due to 3P inplane rotating-system hub forces. In addition, detuning effects of the 3P monofilar modal response were small due to the nonlinearities and tuning pin slippage. However, attenuation of 4P hub motions due to 5P inplane hub forces was poor. The performance of the 5P monofilar modal response was degraded by torsional motion of the dynamic mass relative to the support arm which resulted in binding of the dynamic components. Analytical design studies were performed to evaluate this torsional motion problem. An alternative design is proposed which may alleviate the torsional motion of the dynamic mass.

Duh, J.↗
DOE OSTI2024

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Momentum Transfer Studies and Studies of Linear and Nonlinear Optical Properties of Metal Colloids and Semiconductor Quantum Dots

Phase 1 of this work involved design work on a momentum transfer device. The progress on design and testing will be presented. Phase 2 involved the systematic study of the MPD thruster for dual uses. Though it was designed as a thruster for space vehicles, the characteristics of the plasma make it an excellent candidate for industrial applications. This project sought to characterize the system for use in materials processing and characterization. The surface modification on ZnCdTe, CdTe, and ZnTe will be presented. Phase 3 involved metal colloids and semiconductor quantum dots. One aspect of this project involves a collaborative effort with the Solid State Division of ORNL. The thrust behind this research is to develop ion implantation for synthesizing novel materials (quantum dots wires and wells, and metal colloids) for applications in all optical switching devices, up conversion, and the synthesis of novel refractory materials. The ions of interest are Au, Ag, Cd, Se, In, P, Sb, Ga, and As. The specific materials of interest are: CdSe, CdTe, InAs, GaAs, InP, GaP, InSb, GaSb, and InGaAs. A second aspect of this research program involves using porous glass (25-200 A) for fabricating materials of finite size. The results of some of this work will also be reported.

Collins, W. E.↗
DOE OSTI2025

Basic Research Needs for Inverse Methods for Complex Systems under Uncertainty

Inverse problems, which aim to infer unknown properties of a system using experimental and observational data, are central to addressing many of the U.S. Department of Energy’s (DOE) most critical scientific and engineering challenges. Accurate, computationally efficient, and data-efficient solutions to inverse problems are essential for advancing DOE mission-critical science drivers, including analyzing data from large-scale experimental facilities, optimizing fusion reactor performance, accelerating materials discovery, enhancing geophysical imaging, improving wildfire predictions, and enabling autonomous systems and digital twins. However, these problems are becoming increasingly complex, often involving nonlinear, highdimensional, and interconnected systems and models that span multiple physics and scales, while relying on data with varying quantity, quality, and information content. Compounding these challenges is the uncertainty inherent in DOE-relevant systems, where errors in inputs, noise in data, incompleteness of data, and discrepancies between models and reality constrain the accuracy and precision of solutions. At the same time, the convergence of recent scientific computing trends—scientific machine learning, artificial intelligence, and computing advances such as exascale computing—is creating unprecedented opportunities for tackling these challenges. The cross-cutting nature of inverse problems, combined with their growing complexity and rapidly evolving data and algorithmic demands, strongly motivates the formulation of a prioritized research agenda to maximize their capabilities and impact. In response to this need, DOE’s Advanced Scientific Computing Research (ASCR) program in the Office of Science convened the Workshop on Basic Research Needs for Inverse Problems for Complex Systems Under Uncertainty in June 2025. This workshop brought together experts across disciplines to identify grand challenges and major opportunities in the field. Through collaborative discussions, the workshop defined transformative research directions aimed at addressing the mathematical, statistical, and computational challenges posed by inverse problems under uncertainty. As a result of these efforts, four priority research directions (PRDs) were identified to guide future research and development in this area. These PRDs, summarized below, represent a roadmap for advancing the foundational science and mathematics of inverse problems, enabling robust, scalable, and uncertainty-aware solutions that are critical for DOE applications.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Assessment of SPAR elements and formulation of some basic 2-D and 3-D elements for use with testbed generic element processor

The initial Computational Structural Mechanics (CSM) Testbed was based on Level 13 of the SPAR finite element computer program. Until recently, the element library of the Testbed has been limited to those elements in Level 13 of SPAR. The development of a generic element processor has enabled element researchers to develop, implement and assess element formulations with relative ease. An assessment of new elements as well as the existing SPAR Level 13 elements has revealed some definite shortcomings with the SPAR Level 13 2-D and 3-D elements. The SPAR S81 solid element does not pass the patch test problem proposed by MacNeal-Harder. These deficiencies are identified here. The 2-D elements, however, seem to perform well taking into account the limitations imposed by the theory used to formulate them, (i.e., thin plates only). Common deficiencies of the 2-D and 3-D elements in SPAR have to do with their adaptability to the nonlinear analysis utilities developed by Lockheed Palo Alto Research Lab. Also, the EFIL format of the SPAR element data does not conform to the standard format of the Testbed.

Aminpour, Mohammad A.↗
DOE OSTI2022

Space Charge Driven Emittance Growth and the Effect of Octupoles in IOTA

The Integrable Optics Test Accelerator (IOTA) at Fermilab is a small machine dedicated to a broad frontier accelerator physics program. An important aspect of this program is to investigate the potential benefits of the resonance free tune spread achievable with integrable optics to store and accelerate high intensity proton beams for which space charge is significant. In this context, a good understanding of proton beam emittance growth and particle loss mechanisms is essential. Assuming nominal design parameters, simulations show that for a bunched beam, the bulk of emittance growth takes place immediately following injection, typically within tens of turns. We attempt to account for this growth using a simplified RMS mismatch theory; some of its limitations and possible improvements are briefly discussed. We then compare theoretical predictions to simulations performed using the PIC code pyORBIT. Further exploring ways to mitigate emittance growth and reduce particle loss, we compare two beam matching strategies: (1) matching at the injection point (2) matching at the center of the nonlinear (octupole) insertion region where $\beta_x = \beta_y$. To observe how nonlinearity affects emittance growth and whether it dominates growth due to mismatch, we track two different distributions. Finally, we explore the potential of using octupoles in a quasi-integrable configuration to mitigate growth using a variety of initial distributions both at reduced and full intensities.

43 PARTICLE ACCELERATORS↗
NASA NTRS1991

Linear and nonlinear dynamic analysis by boundary element method

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons with available analytical and numerical results, the stability and high accuracy of these dynamic analysis techniques are established.

Ahmad, Shahid↗