Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear model reduction”

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 145 records · Page 8

Space Shuttle main engine nozzle-steerhorn dynamics

On two occasions during the Space Shuttle main engine development, the LH2 feedline (called the steerhorn, because of its shape) failed during the cutoff transient. A dynamic test was undertaken, and an analytical model was developed and correlated to the dynamic test. Detailed models of the tube bundle were required to obtain the equivalent shell coefficients. All-shell models of the nozzle wall were found better than beam-shell models. The most difficult part of the structure to simulate was the felt-metal pad between the feedline and its mount, which introduced nonlinear stiffness and damping and led to the use of separate low amplitude and high amplitude models. The total structure was found to have 400 modes in the frequency range of interest, 0 to 500 Hz. Good test analysis correlation was obtained and a modified feedline configuration was found to demonstrate a 40% reduction of response stress from the original configuration.

Kiefling, L.↗

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Architecture and Performance Predictions

A peak-seeking control approach for real-time trim configuration optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control approach is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an F/A-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are controlled for optimization of fuel flow. This paper presents the design and integration of this peak-seeking controller on a modified NASA F/A-18 airplane with research flight control computers. A research flight was performed to collect data to build a realistic model of the performance function and characterize measurement noise. This model was then implemented into a nonlinear six-degree-of-freedom F/A-18 simulation along with the peak-seeking control algorithm. With the goal of eventual flight tests, the algorithm was first evaluated in the improved simulation environment. Results from the simulation predict good convergence on minimum fuel flow with a 2.5-percent reduction in fuel flow relative to the baseline trim of the aircraft.

performance optimization↗

Hybrid learning techniques for scientific data reduction with performance guarantees

The research initiatives supported by the U.S. Department of Energy (DOE) Grant DE-SC0022265 are fundamentally aimed at pioneering advanced machine learning (ML) techniques for scientific data compression within high-performance computing (HPC) environments. This comprehensive body of work addresses the critical challenge posed by the exponential growth of data generated by scientific simulations in domains such as fusion energy, climate modeling, and computational fluid dynamics (CFD). A core objective is to develop compression algorithms that achieve substantial data reduction—often by orders of magnitude—while rigorously ensuring the fidelity of both the primary data (PD) and scientifically crucial derived quantities of interest (QoI). The methodologies deployed under this grant integrate sophisticated deep learning architectures, prominently featuring autoencoders, advanced generative models like conditional diffusion, and hybrid learning techniques. Key innovations include the development of Guaranteed Autoencoders (GAE) and the Guaranteed Conditional Diffusion with Tensor Correction (GCDTC) framework, which provide explicit, instance-level error bounds on reconstructed data. Furthermore, specialized strategies such as nonlinear constraint satisfaction are employed to preserve the integrity of QoI, a vital requirement for the trustworthiness of downstream scientific analyses. This research also focuses on the design and implementation of scalable, GPU-accelerated software pipelines that seamlessly integrate into existing HPC workflows, ensuring both computational efficiency and practical applicability. The CAESAR framework, for example, unifies foundation and generative models to create an adaptive and efficient compression solution for spatio-temporal scientific data. Collectively, these efforts represent a significant advancement in mitigating the scientific data deluge, enabling more effective data management, accelerated scientific discovery, and optimized utilization of HPC resources.

97 MATHEMATICS AND COMPUTING↗

Sheared wing-tip aerodynamics - Wind-tunnel and computational investigation

Computational and experimental performance benefits are presented for a high-aspect-ratio unswept wing configuration with sheared tips. The sheared tip is a highly swept and highly tapered surface located in the same plane as the inboard wing panel to which it is attached. The compuational results were obtained with an inviscid surface panel method that models the nonlinear influence of the trailing wake. Both wind-tunnel and calculated results were obtained for a 12-ft span wing model with various wing-tip configurations. The computational and experimental data are in fair agreement and demonstrate that sheared wing tips can reduce induced drag at cruise and climb lift coefficients. The drag reduction is the result of wake deformation effects and changes in spanwise load distribution. Wind-tunnel measured longitudinal and lateral directional stability characteristics are also presented for the various wing-tip layouts.

Vijgen, P. M. H. W.↗

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Architecture and Performance Predictions

A peak-seeking control approach for real-time trim configuration optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control approach is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an FA-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are controlled for optimization of fuel flow. This presentation presents the design and integration of this peak-seeking controller on a modified NASA FA-18 airplane with research flight control computers. A research flight was performed to collect data to build a realistic model of the performance function and characterize measurement noise. This model was then implemented into a nonlinear six-degree-of-freedom FA-18 simulation along with the peak-seeking control algorithm. With the goal of eventual flight tests, the algorithm was first evaluated in the improved simulation environment. Results from the simulation predict good convergence on minimum fuel flow with a 2.5-percent reduction in fuel flow relative to the baseline trim of the aircraft.

peak-seeking↗

Final report- UFL - RAPIDS2: A SciDAC Institute for Computer Science, Data, and Artificial Intelligence

The research initiatives supported by the U.S. Department of Energy (DOE) Grant DE-SC0022265 are fundamentally aimed at pioneering advanced machine learning (ML) techniques for scientific data compression within high-performance computing (HPC) environments. This comprehensive body of work addresses the critical challenge posed by the exponential growth of data generated by scientific simulations in domains such as fusion energy, climate modeling, and computational fluid dynamics (CFD). A core objective is to develop compression algorithms that achieve substantial data reduction—often by orders of magnitude—while rigorously ensuring the fidelity of both the primary data (PD) and scientifically crucial derived quantities of interest (QoI). The methodologies deployed under this grant integrate sophisticated deep learning architectures, prominently featuring autoencoders, advanced generative models like conditional diffusion, and hybrid learning techniques. Key innovations include the development of Guaranteed Autoencoders (GAE) and the Guaranteed Conditional Diffusion with Tensor Correction (GCDTC) framework, which provide explicit, instance-level error bounds on reconstructed data. Furthermore, specialized strategies such as nonlinear constraint satisfaction are employed to preserve the integrity of QoI, a vital requirement for the trustworthiness of downstream scientific analyses. This research also focuses on the design and implementation of scalable, GPU-accelerated software pipelines that seamlessly integrate into existing HPC workflows, ensuring both computational efficiency and practical applicability. The CAESAR framework, for example, unifies foundation and generative models to create an adaptive and efficient compression solution for spatio-temporal scientific data. Collectively, these efforts represent a significant advancement in mitigating the scientific data deluge, enabling more effective data management, accelerated scientific discovery, and optimized utilization of HPC resources.

97 MATHEMATICS AND COMPUTING↗

MHD modeling of magnetotail instability for localized resistivity

We present results of a three-dimensional magnetohydrodynamic (MHD) simulation of magnetotail evolution initiated by a sudden occurrence or increase of spatially localized resistivity as the major expected concequence of some localized microinstability. Because of the absence of a quantitative model, possible variations of resistivity levels with current density, or the reduction thereof, are not incorporated in the present investigation. The emphasis of the study is on an investigation of the changes to the overall evolution brought about by this localization, in particular, on the disruption and diversion of the cross-tail current and the nonlinear evolution of the magnetotail instability. The immediate consequences of the occurrence of the localized resistance and the resulting electric field are a reduction and diversion of the electric current around the region of high resistivity, associated with an increase of B(sub z) ('dipolarization') at the earthward edge and a decrease of B(sub z) at the tailward edge of this region. These effects, however, are localized and do not involve a reduction of the total cross-tail current and hence do not lead to the global development of a 'substorm current wedge,' which includes not only the reduction of the cross-tail current but also the buildup of a global field-aligned current system of 'regional 1' type (toward the Earth on the dawnside and away on the duskside of the tail). Such signatures develop at a later time, as consequences of a three-dimensional tearing instability, which is triggered by the occurrence of the resistivity. These features are found in combination with plasmoid formation and ejection, quite similar to results of earlier simulations with uniform resistivity. Differences are found in the timescale of the evolution, which tends to be shorter for localized resistivity, and in the propagation of the dipolarization effects in the equatorial plane. Whereas for uniform resistivity the temporal increase in B(sub z) tends to propagate tailward, apparently due to a pileup effect in the near-Earth region, an earthward propagation is found for the localized resistivity. This propagation results from an earthward convection of the gradient in B(sub z), which is set up by the reconnection process further tailward.

Hesse, Michael↗

An Anisotropic Yield and Damage Material Model to Improve the Contact Pressure Analysis in a Biomass Shredding System

Size reduction systems used in biomass processing break biomass into smaller pieces by utilizing the kinetic energy from the sharp rotating blades. Abrasive and/or erosive wear caused by biomass comminution results in blade wear of the sharp edged cutters, deteriorating the process efficiency. Here, this study aims to optimize the blade design and improve the system efficiency by attempting to understand the interactions between the blades and biomass particles. Since real-time monitoring of these interactions is impractical during operation, mechanical simulations offer a viable alternative for investigating the shredding process. Yet, the irregular geometry and complex mechanical properties of biomass—such as the anisotropic nature of woodchips and their nonlinear fracture behavior—pose significant challenges for accurately simulating contact pressure. In this work an anisotropic yield material model, along with a damage initiation and evolution function, is applied to the woodchip particle to study the contact pressure on shredder blade, offering a scientific basis for improved blade design and process efficiency. This approach can be extended to other biomass processing systems with similar anisotropic feedstocks, making it a valuable tool for advancing sustainable biomass utilization.

09 - BIOMASS FUELS↗

Progressive Failure And Life Prediction of Ceramic and Textile Composites

An engineering approach to predict the fatigue life and progressive failure of multilayered composite and textile laminates is presented. Analytical models which account for matrix cracking, statistical fiber failures and nonlinear stress-strain behavior have been developed for both composites and textiles. The analysis method is based on a combined micromechanics, fracture mechanics and failure statistics analysis. Experimentally derived empirical coefficients are used to account for the interface of fiber and matrix, fiber strength, and fiber-matrix stiffness reductions. Similar approaches were applied to textiles using Repeating Unit Cells. In composite fatigue analysis, Walker's equation is applied for matrix fatigue cracking and Heywood's formulation is used for fiber strength fatigue degradation. The analysis has been compared with experiment with good agreement. Comparisons were made with Graphite-Epoxy, C/SiC and Nicalon/CAS composite materials. For textile materials, comparisons were made with triaxial braided and plain weave materials under biaxial or uniaxial tension. Fatigue predictions were compared with test data obtained from plain weave C/SiC materials tested at AS&M. Computer codes were developed to perform the analysis. Composite Progressive Failure Analysis for Laminates is contained in the code CPFail. Micromechanics Analysis for Textile Composites is contained in the code MicroTex. Both codes were adapted to run as subroutines for the finite element code ABAQUS and CPFail-ABAQUS and MicroTex-ABAQUS. Graphic user interface (GUI) was developed to connect CPFail and MicroTex with ABAQUS.

Xue, David Y.↗

A Computational and Experimental Study of Nonlinear Aspects of Induced Drag

Despite the 80-year history of classical wing theory, considerable research has recently been directed toward planform and wake effects on induced drag. Nonlinear interactions between the trailing wake and the wing offer the possibility of reducing drag. The nonlinear effect of compressibility on induced drag characteristics may also influence wing design. This thesis deals with the prediction of these nonlinear aspects of induced drag and ways to exploit them. A potential benefit of only a few percent of the drag represents a large fuel savings for the world's commercial transport fleet. Computational methods must be applied carefully to obtain accurate induced drag predictions. Trefftz-plane drag integration is far more reliable than surface pressure integration, but is very sensitive to the accuracy of the force-free wake model. The practical use of Trefftz plane drag integration was extended to transonic flow with the Tranair full-potential code. The induced drag characteristics of a typical transport wing were studied with Tranair, a full-potential method, and A502, a high-order linear panel method to investigate changes in lift distribution and span efficiency due to compressibility. Modeling the force-free wake is a nonlinear problem, even when the flow governing equation is linear. A novel method was developed for computing the force-free wake shape. This hybrid wake-relaxation scheme couples the well-behaved nature of the discrete vortex wake with viscous-core modeling and the high-accuracy velocity prediction of the high-order panel method. The hybrid scheme produced converged wake shapes that allowed accurate Trefftz-plane integration. An unusual split-tip wing concept was studied for exploiting nonlinear wake interaction to reduced induced drag. This design exhibits significant nonlinear interactions between the wing and wake that produced a 12% reduction in induced drag compared to an equivalent elliptical wing at a lift coefficient of 0.7. The performance of the split-tip wing was also investigated by wing tunnel experiments. Induced drag was determined from force measurements by subtracting the estimated viscous drag, and from an analytical drag-decomposition method using a wake survey. The experimental results confirm the computational prediction.

Smith, Stephen C.↗

Development of the X-57 Mod III/IV Piloted Simulator and Discussion of the Resultant Flighting Qualities Predictions

The all-electric X-57 Mod III aircraft was designed to demonstrate wingtip propulsion drag reduction benefits on a high-aspect-ratio wing. The X-57 Mod IV aircraft was designed to expand on the Mod III version and demonstrate the benefits of distributed electric propulsion. This paper discusses the development of theall-electricX-57 Mod III and Mod IV fixed-base pilot-in-the-loop nonlinear simulator. The paper describes the model development and the simulator cockpit construction for uses such as flight training and flying qualities analysis. Results from a stability and flying qualities analysis for both Mod III and Mod IV are presented. These results predict that both Mod III and Mod IV are stable throughout their flight envelopes. Mod III is predicted to have satisfactory flying qualities while Mod IV is predicted to have areas of adequate flying qualities.

Ryan Wallace↗

Parameter identification for a robotic manipulator arm

The development is described of a nonlinear dynamic model for large oscillations of a robotic manipulator arm about a single joint. Optimization routines are formulated and implemented for the identification of electrical and physical parameters from dynamic data taken from an industrial robot arm. Special attention is given to the role of sensitivity in the formulation of robust models of this motion. The importance of actuator effects in the reduction of sensitivity is established and used to develop an electromechanical model of the manipulator system.

Brewer, D. W.↗

Parametric analysis of swept-wing geometry with sheared wing tips

A computational parameter study is presented of potential reductions in induced drag and increases in lateral-directional stability due to sheared wing tips attached to an untwisted wing of moderate sweep and aspect ratio. Sheared tips are swept and tapered wing-tip devices mounted in the plane of the wing. The induced-drag results are obtained using an inviscid, incompressible surface-panel method that models the nonlinear effects due to the deflected and rolled-up wake behind the lifting surface. The induced-drag results with planar sheared tips are compared to straight-tapered tip extensions and nonplanar winglet geometries. The lateral-directional static-stability characteristics of the wing with sheared tips are estimated using a quasi-vortex-lattice method. For certain combinations of sheared-tip sweep and taper, both the induced efficiency of the wing and the relevant static-stability derivatives are predicted to increase compared to the wing with a straight-tapered tip modification.

Fremaux, C. M.↗

Horizontal vestibuloocular reflex evoked by high-acceleration rotations in the squirrel monkey. IV. Responses after spectacle-induced adaptation

The horizontal angular vestibuloocular reflex (VOR) evoked by sinusoidal rotations from 0.5 to 15 Hz and acceleration steps up to 3,000 degrees /s(2) to 150 degrees /s was studied in six squirrel monkeys following adaptation with x2.2 magnifying and x0.45 minimizing spectacles. For sinusoidal rotations with peak velocities of 20 degrees /s, there were significant changes in gain at all frequencies; however, the greatest gain changes occurred at the lower frequencies. The frequency- and velocity-dependent gain enhancement seen in normal monkeys was accentuated following adaptation to magnifying spectacles and diminished with adaptation to minimizing spectacles. A differential increase in gain for the steps of acceleration was noted after adaptation to the magnifying spectacles. The gain during the acceleration portion, G(A), of a step of acceleration (3,000 degrees /s(2) to 150 degrees /s) increased from preadaptation values of 1.05 +/- 0.08 to 1.96 +/- 0.16, while the gain during the velocity plateau, G(V), only increased from 0.93 +/- 0.04 to 1.36 +/- 0.08. Polynomial fits to the trajectory of the response during the acceleration step revealed a greater increase in the cubic than the linear term following adaptation with the magnifying lenses. Following adaptation to the minimizing lenses, the value of G(A) decreased to 0.61 +/- 0.08, and the value of G(V) decreased to 0.59 +/- 0.09 for the 3,000 degrees /s(2) steps of acceleration. Polynomial fits to the trajectory of the response during the acceleration step revealed that there was a significantly greater reduction in the cubic term than in the linear term following adaptation with the minimizing lenses. These findings indicate that there is greater modification of the nonlinear as compared with the linear component of the VOR with spectacle-induced adaptation. In addition, the latency to the onset of the adapted response varied with the dynamics of the stimulus. The findings were modeled with a bilateral model of the VOR containing linear and nonlinear pathways that describe the normal behavior and adaptive processes. Adaptation for the linear pathway is described by a transfer function that shows the dependence of adaptation on the frequency of the head movement. The adaptive process for the nonlinear pathway is a gain enhancement element that provides for the accentuated gain with rising head velocity and the increased cubic component of the responses to steps of acceleration. While this model is substantially different from earlier models of VOR adaptation, it accounts for the data in the present experiments and also predicts the findings observed in the earlier studies.

NASA Discipline Neuroscience↗

Uncertainty propagation and sensitivity analysis for constrained optimization of nuclear waste vitrification

Abstract The vitrification of high‐level waste (HLW) by heating a mixture of glass‐forming chemicals (GFCs) with the waste can be improved using a constrained optimization problem. This study explores how different uncertainty propagation (UP) methods implemented with the optimization process can affect the glass formulation of nuclear waste glasses. UP is the effort of propagating uncertain inputs through a system to understand and quantify output distributions. Uncertainty intervals are crafted from output distributions to inform the optimization algorithm. UP is often implemented with Monte Carlo (MC) sampling for large nonlinear systems, which can be difficult to implement within a constrained optimization algorithm that requires derivative information. Other UP methods often used for optimization under uncertainty (OUU) can be designed to work within an established constrained optimization framework. Methods of UP are evaluated in this study including iterative sampling approaches, first‐order approximations, and surrogate modeling with machine learning (ML). A method of dimensional reduction based on global sensitivity analysis is introduced to support the UP methods for the large dimensionality of the problem. Analytical UP methods able to achieve similar optimums 10 times faster than the baseline MC approach, and produce 93.9% similar output distributions are reported.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗

Computational aspects of maximum likelihood estimation and reduction in sensitivity function calculations

This paper discusses numerical aspects of computing maximum likelihood estimates for linear dynamical systems in state-vector form. Different gradient-based nonlinear programming methods are discussed in a unified framework and their applicability to maximum likelihood estimation is examined. The problems due to singular Hessian or singular information matrix that are common in practice are discussed in detail and methods for their solution are proposed. New results on the calculation of state sensitivity functions via reduced order models are given. Several methods for speeding convergence and reducing computation time are also discussed.

Gupta, N. K.↗

Flight flutter testing technology at Grumman

Analysis techniques used in the automated telemetry station (ATS) for on line data reduction are encompassed in a broad range of software programs. Concepts that form the basis for the algorithms used are mathematically described. The control the user has in interfacing with various on line programs is discussed. The various programs are applied to an analysis of flight data which includes unimodal and bimodal response signals excited via a swept frequency shaker and/or random aerodynamic forces. A nonlinear response error modeling analysis approach is described. Preliminary results in the analysis of a hard spring nonlinear resonant system are also included.

Perangelo, H. J.↗