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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Airborne synthetic aperture radar observations and simulations for waves in ice

The Canada Centre for Remote Sensing CV-580 aircraft collected C-band SAR data over the marginal ice zone off the east coast of Newfoundland during the Labrador Ice Margin Experiment (LIMEX) in March 1989. One component of the LIMEX '89 program was the study of ocean waves penetrating the marginal ice zone. We consider nearly coincidental observations of waves in ice by airborne SAR and wave-induced ice motion measurements. We explain the wave patterns observed in the SAR imagery, and the corresponding SAR image spectra, in terms of SAR wave imaging models. These include the well-known tilt cross-section modulation, linear, quasi-linear, and nonlinear velocity bunching forward mapping models (FMMs), and the assertion that the concept of coherence time limitation applies differently to the cases of waves in ice and open water. We modify the concept of the scene coherence time to include two parts: first, a decorrelation time deduced from the inherent azimuth cutoff in the nonlinear velocity bunching FMM; and second, the intrinsic scene coherence time which is a measure of the time scale over which an open water Bragg scattering patch retains its phase structure. Either of these coherence time scales could dominate the SAR image formation process, depending upon the environmental conditions (the wave spectrum and the wind speed, for example). Observed SAR image spectra and forward mapped ice motion package spectra are favorably compared.

Vachon, Paris W.↗

NASA NTRSDate not supplied

Fusion of Test and Analysis: Artemis I Booster to Mobile Launcher Interface Validation

NASA is in the midst of bold and exciting next steps in human exploration and spaceflight. The designs of the new Space Launch System (SLS), the Orion spacecraft and the Exploration Ground Systems (EGS) for vehicle processing and launch are essentially complete and there has been significant progress in manufacturing and assembly of specific hardware for the Artemis I and Artemis II missions. Equally as important, the program level and integrated system level testing and analyses are also well underway to support integrated verification, validation, and Certificate of Flight Readiness (CoFR) for Artemis I. Testing and analysis are key to addressing technical challenges faced by the Artemis missions. Building block approaches are required that provide the right balance between component, element, and/or system level testing that satisfies verification and validation objectives where uncertainties are quantified and minimized. Artemis I is a system of systems that requires a fusion of test and analysis that adeptly characterizes critical interfaces between major program elements. An example of this fusion involves characterizing the interface between the SLS booster and the Mobile Launcher (ML) Vertical Support Post (VSP) interfaces. Proper characterization of this interface represents a number of challenges beginning with the fact that it is a mating of ground support structure in the form of a civil structure to flight hardware. Both sides of the interface are built to different construction standards, but are governed by interface requirements to ensure compatibility when mated. From past program experience, the flexibility at the booster to ML interface is critical in developing accurate prelaunch stacking and cryogenic preloads, squat loads, and pad separation release of preloads and squat loads. This same premise holds for Artemis I. To characterize the asymmetric characteristics at this interface, careful consideration of static forces due to gravity loading with the commensurate effects due to leveling during booster stacking (i.e., spacing and shimming) and nonlinear geometric forces are necessary for inclusion in pre-test assessments. This paper will look at these issues for the upcoming Booster Pull Test in which two boosters will be installed on the ML and one of these boosters will undergo static lateral loading followed afterwards with dynamic excitation into resonance and free-decay. This paper evaluates the booster to ML interface characteristics by characterizing the interface flexibility between the booster aft skirt and the ML VSP interfaces. Furthermore, this paper methodically evaluates the effect of the following on the test outcome: gravitational effects on the booster and ML, the effects of VSP leveling, spacing, and shimming under gravitational loading during booster stacking, the effect of geometric nonlinear follower force due to cg offset as booster is laterally displaced, and the system coupling between the booster under test, ML, and the second booster. Simulated results for a static load pull and dynamic excitation provide insight into the differences in measurement responses when boundary conditions and geometric conditions are included and not included.

Joel W Sills Jr.↗

DOE OSTI2020

An Optimization-Based Planning Tool for On-Demand Mobility Service Operations

Regions worldwide are adopting and exploring low-speed automated electric shuttle (AES) service as an on-demand shared mobility service in dense geofenced urban areas. Building on this concept, the National Renewable Energy Laboratory (NREL) recently developed the Automated Mobility District (AMD) toolkit. The AMD toolkit—comprising of a travel micro-simulation model and an energy estimation model—estimates the mobility and energy impacts of a given shuttle configuration within an AMD. Early-stage AMD deployments need to find optimal operational configurations that include: (a) passenger capacity of an AES, (b) time-dependent routes, and (c) fleet size (AES units) to satisfy the demand for the region. This research extends the AMD toolkit functionality by developing an optimization-based planning module that will assist in the operations of AES units. We developed a constrained mixed-integer program accounting for passenger waiting time, battery range, and passenger capacity of AES units. For scalability, we demonstrated the Tabu search-based solution technique for a real-world network—a proposed AMD deployment in Greenville, South Carolina, USA. Compared to rule-based operations, our developed solution yields higher travel time and energy savings for the network at different demand levels. The sensitivity analyses for waiting time thresholds indicate nonlinearity in the system performance, underscoring the need to meet shared-use mobility user-level expectations. The developed optimization framework can be adapted and extended to accommodate different categories of shared-use on-demand mobility services.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗

NASA NTRS1974

Manual control

Man's nonlinear characteristics and his use of control with compatible and incompatible multiple inputs, both visual and vestibular were studied. Experiments were also made with pulse and bang-bang controllers and the effects of sudden changes in control stick mechanical impedance. Closing the loop through the dynamics of the controlled vehicle allowed experiments on the limits of control of unstable vehicles with and without motion cues. The inverted pendulum controlled element, programmed as a self pacing element, was used extensively as a scalar performance index. In addition, the motorbike equations of motion were studied with regard to required human equalization. Abstracts are included for a series of published data on manual control.

Source record↗

NASA NTRS1977

Atmospheric model development in support of SEASAT. Volume 3: Prediction models

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric models are developed. The four models use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical. Conservation forms of the difference equations are integrated using either a fifteen or four minute time step on a 381 km or 127 km grid (at 60 deg N) for the 63 x 63 or 187 x 187 models respectively. A nonlinear pressure smoother is used along with momentum and temperature diffusion to help control computational noise. The horizontal boundary conditions are insulated slippery walls. Centered time differencing with time averaging of the pressure gradient force term in the momentum equations is used. Robert time filtering on of the temperature and moisture solutions is used for computational stability. The moisture and heat source/sink terms are modeled. Terms representing evaporation and large scale condensation, sensible heat exchange, parameterized cumulus convection and precipitation, and solar and terrestrial radiation are included. Initialization of the models is based on a pattern conservation technique to obtain objective analysis of the state parameter structure from the surface to 50 mb.

Wellck, R. E.↗

NASA NTRS1981

A real time Pegasus propulsion system model for VSTOL piloted simulation evaluation

A real time propulsion system modeling technique suitable for use in man-in-the-loop simulator studies was developd. This technique provides the system accuracy, stability, and transient response required for integrated aircraft and propulsion control system studies. A Pegasus-Harrier propulsion system was selected as a baseline for developing mathematical modeling and simulation techniques for VSTOL. Initially, static and dynamic propulsion system characteristics were modeled in detail to form a nonlinear aerothermodynamic digital computer simulation of a Pegasus engine. From this high fidelity simulation, a real time propulsion model was formulated by applying a piece-wise linear state variable methodology. A hydromechanical and water injection control system was also simulated. The real time dynamic model includes the detail and flexibility required for the evaluation of critical control parameters and propulsion component limits over a limited flight envelope. The model was programmed for interfacing with a Harrier aircraft simulation. Typical propulsion system simulation results are presented.

Mihaloew, J. R.↗

NASA NTRS1988

Dynamic analysis of an unrestrained, rotating structure through nonlinear simulation

An alternate derivation based on Hamilton's principle of equations of motion of a flexible body for use in nonlinear simulation modeling is presented. The equations include terms representing variation in the inertial properties of the body resulting from flexibility as well as terms coupling rigid body and flexible momenta that arise from satisfaction of only the linear part of the mean axis conditions. An approximate method for including the centrifugal stiffening in a way that permits a variable spin vector is given. The equations are exercised in a simulation program on an example structure initially spinning about its unstable principal axis in a gravity-free vacuum.

Zeiler, Thomas A.↗

NASA NTRS1988

A finite element method for time varying geometry in multibody structures

A three-dimensional finite element formulation using convected coordinates is presented for the multibody dynamics of truss-like configurations. Unlike existing formulations, the present one does not superimpose nonlinear rigid body kinematics with linear structural mode shapes, an approach that has recently been shown to be grossly inaccurate under certain conditions. Instead, the finite element method is extended to treat large motions/deformations. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joints. For the planar spin-up of a flexible beam, results are compared with those derived from a commercially available computer program. The two programs predict nearly identical results.

Housner, J. M.↗

NASA NTRS1991

A three-dimensional equivalent domain integral method for cracked solids

Three-dimensional FEM analysis is presently undertaken via the general formulation of the equivalent-domain integral method for mixed-mode fracture problems in general anisotropic and nonlinear materials. In mixed-mode integrals, the total J-integral is split into three separate modes of fracture; direct and decomposition methods are effective in this separation of the modes, as presently demonstrated by their application to several pure and mixed-mode fracture problems. The results are found to agree well with those reported in the literature, and suggest the usefulness of the method as a FEM-program postprocessing subroutines.

Shivakumar, K. N.↗

NASA NTRS1995

Supersonic civil airplane study and design: Performance and sonic boom

Since aircraft configuration plays an important role in aerodynamic performance and sonic boom shape, the configuration of the next generation supersonic civil transport has to be tailored to meet high aerodynamic performance and low sonic boom requirements. Computational fluid dynamics (CFD) can be used to design airplanes to meet these dual objectives. The work and results in this report are used to support NASA's High Speed Research Program (HSRP). CFD tools and techniques have been developed for general usages of sonic boom propagation study and aerodynamic design. Parallel to the research effort on sonic boom extrapolation, CFD flow solvers have been coupled with a numeric optimization tool to form a design package for aircraft configuration. This CFD optimization package has been applied to configuration design on a low-boom concept and an oblique all-wing concept. A nonlinear unconstrained optimizer for Parallel Virtual Machine has been developed for aerodynamic design and study.

Cheung, Samson↗

NASA NTRS2002

Unsteady Pressures in a Transonic Fan Cascade Due to a Single Oscillating Airfoil

An extensive set of unsteady pressure data was acquired along the midspan of a modern transonic fan blade for simulated flutter conditions. The data set was acquired in a nine-blade linear cascade with an oscillating middle blade to provide a database for the influence coefficient method to calculate instantaneous blade loadings. The cascade was set for an incidence of 10 dg. The data were acquired on three stationary blades on each side of the middle blade that was oscillated at an amplitude of 0.6 dg. The matrix of test conditions covered inlet Mach numbers of 0.5, 0.8, and 1.1 and the oscillation frequencies of 200, 300, 400, and 500 Hz. A simple quasiunsteady two-dimensional computer simulation was developed to aid in the running of the experimental program. For high Mach number subsonic inlet flows the blade pressures exhibit very strong, low-frequency, self-induced oscillations even without forced blade oscillations, while for low subsonic and supersonic inlet Mach numbers the blade pressure unsteadiness is quite low. The amplitude of forced pressure fluctuations on neighboring stationary blades strongly depends on the inlet Mach number and forcing frequency. The flowfield behavior is believed to be governed by strong nonlinear effects due to a combination of viscosity, compressibility, and unsteadiness. Therefore, the validity of the quasi-unsteady simplified computer simulation is limited to conditions when the flowfield is behaving in a linear, steady manner. Finally, an extensive set of unsteady pressure data was acquired to help development and verification of computer codes for blade flutter effects.

Lepicovsky, J.↗

DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

NAS (Numerical Aerodynamic Simulation Program) technical summaries, March 1989 - February 1990

Given here are selected scientific results from the Numerical Aerodynamic Simulation (NAS) Program's third year of operation. During this year, the scientific community was given access to a Cray-2 and a Cray Y-MP supercomputer. Topics covered include flow field analysis of fighter wing configurations, large-scale ocean modeling, the Space Shuttle flow field, advanced computational fluid dynamics (CFD) codes for rotary-wing airloads and performance prediction, turbulence modeling of separated flows, airloads and acoustics of rotorcraft, vortex-induced nonlinearities on submarines, and standing oblique detonation waves.

Source record↗

NASA NTRS1993

Growth of new materials for solid state laser applications. The growth of ZnGeP2 by the vertical Bridgman method

This is the final technical report on the cooperative NASA program 'Growth of New Materials for Solid State Laser Applications,' covering the period from 9-1-86 through 3-31-91. The first two and one half years of the program, from 9-1-86 through 3-31-89, was devoted to the development of new eye-safe laser sources. Single crystal fibers of rare earth doped and co-doped YAG were grown by the laser-heated pedestal growth method, characterized for their structural properties and supplied to NASA/Langley for spectroscopic evaluation. From 4-1-89 though 3-31-91, research focused on the growth of zinc germanium phosphide for nonlinear optical applications in the IR. The vertical Bridgman growth process was evaluated as a means to achieve better control over the short wavelength optical absorption in this material that adversely affects 2 micron-pumped optical parametric oscillators.

Feigelson, R. S.↗

NASA NTRS1999

Crack Growth Simulation and Residual Strength Prediction in Airplane Fuselages

This is the final report for the NASA funded project entitled "Crack Growth Prediction Methodology for Multi-Site Damage." The primary objective of the project was to create a capability to simulate curvilinear fatigue crack growth and ductile tearing in aircraft fuselages subjected to widespread fatigue damage. The second objective was to validate the capability by way of comparisons to experimental results. Both objectives have been achieved and the results are detailed herein. In the first part of the report, the crack tip opening angle (CTOA) fracture criterion, obtained and correlated from coupon tests to predict fracture behavior and residual strength of built-up aircraft fuselages, is discussed. Geometrically nonlinear, elastic-plastic, thin shell finite element analyses are used to simulate stable crack growth and to predict residual strength. Both measured and predicted results of laboratory flat panel tests and full-scale fuselage panel tests show substantial reduction of residual strength due to the occurrence of multi-site damage (MSD). Detailed comparisons of n stable crack growth history, and residual strength between the predicted and experimental results are used to assess the validity of the analysis methodology. In the second part of the report, issues related to crack trajectory prediction in thin shells; an evolving methodology uses the crack turning phenomenon to improve the structural integrity of aircraft structures are discussed, A directional criterion is developed based on the maximum tangential stress theory, but taking into account the effect of T-stress and fracture toughness orthotropy. Possible extensions of the current crack growth directional criterion to handle geometrically and materially nonlinear problems are discussed. The path independent contour integral method for T-stress evaluation is derived and its accuracy is assessed using a p- and hp-version adaptive finite element method. Curvilinear crack growth is simulated in coupon tests and in full-scale fuselage panel tests. Both T-stress and fracture toughness orthotropy are found to be essential to predict the observed crack paths. The analysis methodology and software program (FRANC3D/STAGS) developed herein allows engineers to maintain aging aircraft economically while insuring continuous airworthiness. Consequently, it will improve the technology to support the safe operation of the current aircraft fleet as well as the design of more damage-tolerant aircraft for the next generation fleet.

Chen, Chuin-Shan↗

NASA NTRS1971

Nonlinear, nonlaminar-3D computation of electron motion through the output cavity of a klystron

The equations of motion used in the computation are discussed along with the space charge fields and the integration process. The following assumptions were used as a basis for the computation: (1) The beam is divided into N axisymmetric discs of equal charge and each disc into R rings of equal charge. (2) The velocity of each disc, its phase with respect to the gap voltage, and its radius at a specified position in the drift tunnel prior to the interaction gap is known from available large signal one dimensional programs. (3) The fringing rf fields are computed from exact analytical expressions derived from the wave equation assuming a known field shape between the tunnel tips at a radius a. (4) The beam is focused by an axisymmetric magnetic field. Both components of B, that is B sub z and B sub r, are taken into account. (5) Since this integration does not start at the cathode but rather further down the stream prior to entering the output cavity it is assumed that each electron moved along a laminar path from the cathode to the start of integration.

Albers, L. U.↗

DOE OSTI2026

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗

NASA NTRS1973

Multiple regression technique for Pth degree polynominals with and without linear cross products

A multiple regression technique was developed by which the nonlinear behavior of specified independent variables can be related to a given dependent variable. The polynomial expression can be of Pth degree and can incorporate N independent variables. Two cases are treated such that mathematical models can be studied both with and without linear cross products. The resulting surface fits can be used to summarize trends for a given phenomenon and provide a mathematical relationship for subsequent analysis. To implement this technique, separate computer programs were developed for the case without linear cross products and for the case incorporating such cross products which evaluate the various constants in the model regression equation. In addition, the significance of the estimated regression equation is considered and the standard deviation, the F statistic, the maximum absolute percent error, and the average of the absolute values of the percent of error evaluated. The computer programs and their manner of utilization are described. Sample problems are included to illustrate the use and capability of the technique which show the output formats and typical plots comparing computer results to each set of input data.

Davis, J. W.↗