SEARCH · Search NASA

Results for “Linear and nonlinear operators”

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 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Convex Optimization of Integrated Power-Gas Energy Flow Model With Applications to Probabilistic Energy Flow

Energy flow calculation is a fundamental problem of the integrated power and gas system (IPGS) operation and planning. However, the nonlinear gas flow model introduces major challenges to the energy flow calculation. In this paper, we propose a tractably convex optimization model to solve the energy flow problem in IPGSs. It is demonstrated that the proposed optimization model has the same optimal solution as the original nonlinear steady energy flow model. Also, piecewise linearization is adopted to tightly linearize the nonlinear objective function of the model, which transforms the formulated convex optimization into a linear program one. Thus, the computation complexity of the proposed energy flow model is significantly reduced as compared with the existing methods. In addition, the proposed model can be extended to probabilistic energy flow estimation. Extensive case studies are conducted to validate the effectiveness of the proposed energy flow model using three IPGSs.

42 ENGINEERING↗

NASA NTRS1982

Dynamic stability of a buoyant quad-rotor aircraft

Stability characteristics of a buoyant quad-rotor aircraft (BQRA) in hover and forward flight are examined by considering linear, state-variable, and nonlinear flight simulation models of such a configuration. The effects of carrying a sling load on the vehicle dynamics is predicted by considering a coupled model of the two bodies. Inherent stability characteristics of the vehicle are analyzed and compared with those of a helicopter and an airship in free flight. Typical operational conditions that could lead to vehicle instability are described in the flight envelope of interest.

Nagabhushan, B. L.↗

DOE OSTI2026

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

GP-BayesOpInf

SAND2025-01851O GP-BayesOpInf is a software tool that uses algorithms to combine Gaussian process regression, principal component analysis, and linear Bayesian inference to produce a probabilistic reduced-order model for time-dependent systems. Numerical examples include the compressible Euler equations for an ideal gas, a heat diffusion process with a nonlinear reaction term, and a set of ordinary differential equations describing a compartmental model in epidemiology. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

NASA NTRS1995

Calibration of the advanced microwave sounding unit-A for NOAA-K

The thermal-vacuum chamber calibration data from the Advanced Microwave Sounding Unit-A (AMSU-A) for NOAA-K, which will be launched in 1996, were analyzed to evaluate the instrument performance, including calibration accuracy, nonlinearity, and temperature sensitivity. The AMSU-A on NOAA-K consists of AMSU-A2 Protoflight Model and AMSU-A1 Flight Model 1. The results show that both models meet the instrument specifications, except the AMSU-A1 antenna beamwidths, which exceed the requirement of 3.3 +/- 10%. We also studied the instrument's radiometric characterizations which will be incorporated into the operational calibration algorithm for processing the in-orbit AMSU-A data from space. Particularly, the nonlinearity parameters which will be used for correcting the nonlinear contributions from an imperfect square-law detector were determined from this data analysis. It was found that the calibration accuracies (differences between the measured scene radiances and those calculated from a linear two-point calibration formula) are polarization-dependent. Channels with vertical polarizations show little cold biases at the lowest scene target temperature 84K, while those with horizontal polarizations all have appreciable cold biases, which can be up to 0.6K. It is unknown where these polarization-dependent cold biases originate, but it is suspected that some chamber contamination of hot radiances leaked into the cold scene target area. Further investigation in this matter is required. The existence and magnitude of nonlinearity in each channel were established and a quadratic formula for modeling these nonlinear contributions was developed. The model was characterized by a single parameter u, values of which were obtained for each channel via least-squares fit to the data. Using the best-fit u values, we performed a series of simulations of the quadratic corrections which would be expected from the space data after the launch of AMSU-A on NOAA-K. In these simulations, the cosmic background radiance corresponding to a cold space temperature 2.73K was adopted as one of the two reference points of calibration. The largest simulated nonlinear correction is about 0.3K, which occurs at channel 15 when the instrument temperature is at 38.09 deg C. Others are less than 0.2K in the remaining channels. Possible improvement for future instrument calibration is also discussed.

Mo, Tsan↗

DOE OSTI2024

Synergistic Integration of Multiple Wave Energy Converters with Adaptive Resonance and Offshore Floating Wind Turbines through Bayesian Optimization

We developed a synergistic ocean renewable system where an array of Wave Energy Converters (WEC) with adaptive resonance was collocated with a Floating Offshore Wind Turbine (FOWT) such that the WECs, capturing wave energy through the resonance adapting to varying irregular waves, consequently reduced FOWFT loads and turbine motions. Combining Surface-Riding WECs (SR-WEC) individually designed to feasibly relocate their natural frequency at the peak of the wave excitation spectrum for each sea state, and to obtain the highest capture width ratio at one of the frequent sea states for annual average power in a tens of kilowatts scale with a 15 MW FOWT based on a semi-submersible, Bayesian Optimization is implemented to determine the arrangement of WECs that minimize the annual representation of FOWT’s wave excitation spectra. The time-domain simulation of the system in the optimized arrangement is performed, including two sets of interactions: one set is the wind turbine dynamics, mooring lines, and floating body dynamics for FOWT, and the other set is the nonlinear power-take-off dynamics, linear mooring, and individual WECs’ floating body dynamics. Those two sets of interactions are further coupled through the hydrodynamics of diffraction and radiation. For sea states comprising Annual Energy Production, we investigate the capture width ratio of WECs, wave excitation on FOWT, and nacelle acceleration of the turbine compared to their single unit operations. We find that the optimally arranged SR-WECs reduce the wave excitation spectral area of FOWT by up to 60% and lower the turbine’s peak nacelle acceleration by nearly 44% in highly occurring sea states, while multiple WECs often produce more than the single operation, achieving adaptive resonance with a larger wave excitation spectra for those sea states. The synergistic system improves the total Annual Energy Production (AEP) by 1440 MWh, and we address which costs of Levelized Cost Of Energy (LCOE) can be reduced by the collocation.

Engineering↗

DOE OSTI2020

Multi-objective optimization of peel and shear strengths in ultrasonic metal welding using machine learning-based response surface methodology

Ultrasonic metal welding (UMW) is a solid-state joining technique with varied industrial applications. Despite of its numerous advantages, UMW has a relative narrow operating window and is sensitive to variations in process conditions. As such, it is imperative to quantitatively characterize the influence of welding parameters on the resulting joint quality. The quantification model can be subsequently used to optimize the parameters. Conventional response surface methodology (RSM) usually employs linear or polynomial models, which may not be able to capture the intricate, nonlinear input-output relationships in UMW. Furthermore, some UMW applications call for simultaneous optimization of multiple quality indices such as peel strength, shear strength, electrical conductivity, and thermal conductivity. To address these challenges, this paper develops a machine learning (ML)-based RSM to model the input-output relationships in UMW and jointly optimize two quality indices, namely, peel and shear strengths. The performance of various ML methods including spline regression, Gaussian process regression (GPR), support vector regression (SVR), and conventional polynomial regression models with different orders is compared. A case study using experimental data shows that GPR with radial basis function (RBF) kernel and SVR with RBF kernel achieve the best prediction accuracy. The obtained response surface models are then used to optimize a compound joint strength indicator that is defined as the average of normalized shear and peel strengths. In addition, the case study reveals different patterns in the response surfaces of shear and peel strengths, which has not been systematically studied in the literature. While developed for the UMW application, the method can be extended to other manufacturing processes.

42 ENGINEERING↗

NASA NTRS1983

Flight in low-level wind shear

Results of studies of wind shear hazard to aircraft operation are summarized. Existing wind shear profiles currently used in computer and flight simulator studies are reviewed. The governing equations of motion for an aircraft are derived incorporating the variable wind effects. Quantitative discussions of the effects of wind shear on aircraft performance are presented. These are followed by a review of mathematical solutions to both the linear and nonlinear forms of the governing equations. Solutions with and without control laws are presented. The application of detailed analysis to develop warning and detection systems based on Doppler radar measuring wind speed along the flight path is given. A number of flight path deterioration parameters are defined and evaluated. Comparison of computer-predicted flight paths with those measured in a manned flight simulator is made. Some proposed airborne and ground-based wind shear hazard warning and detection systems are reviewed. The advantages and disadvantages of both types of systems are discussed.

Frost, W.↗

NASA NTRS1998

A Method for Generating Reduced-Order Linear Models of Multidimensional Supersonic Inlets

Simulation of high speed propulsion systems may be divided into two categories, nonlinear and linear. The nonlinear simulations are usually based on multidimensional computational fluid dynamics (CFD) methodologies and tend to provide high resolution results that show the fine detail of the flow. Consequently, these simulations are large, numerically intensive, and run much slower than real-time. ne linear simulations are usually based on large lumping techniques that are linearized about a steady-state operating condition. These simplistic models often run at or near real-time but do not always capture the detailed dynamics of the plant. Under a grant sponsored by the NASA Lewis Research Center, Cleveland, Ohio, a new method has been developed that can be used to generate improved linear models for control design from multidimensional steady-state CFD results. This CFD-based linear modeling technique provides a small perturbation model that can be used for control applications and real-time simulations. It is important to note the utility of the modeling procedure; all that is needed to obtain a linear model of the propulsion system is the geometry and steady-state operating conditions from a multidimensional CFD simulation or experiment. This research represents a beginning step in establishing a bridge between the controls discipline and the CFD discipline so that the control engineer is able to effectively use multidimensional CFD results in control system design and analysis.

Chicatelli, Amy↗

NASA NTRS2004

A Survey Of Earth-Moon Libration Orbits: Stationkeeping Strategies And Intra-Orbit Transfers

Cislunar space is a readily accessible region that may well develop into a prime staging area in the effort to colonize space near Earth or to colonize the Moon. While there have been statements made by various NASA programs regarding placement of resources in orbit about the Earth-Moon Lagrangian locations, there is no survey of the total cost associated with attaining and maintaining these unique orbits in an operational fashion. Transfer trajectories between these orbits required for assembly, servicing, and positioning of these resources have not been extensively investigated. These orbits are dynamically similar to those used for the Sun-Earth missions, but differences in governing gravitational ratios and perturbation sources result in unique characteristics. We implement numerical computations using high fidelity models and linear and nonlinear targeting techniques to compute the various maneuver (Delta)V and temporal costs associated with orbits about each of the Earth-Moon Lagrangian locations (L1, L2, L3, L4, and L5). From a dynamical system standpoint, we speak to the nature of these orbits and their stability. We address the cost of transfers between each pair of Lagrangian locations.

Folta, David↗

DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1989

Dynamic analysis of flexible multibody structures

A system composed of several interconnected elastic components that may experience large angular motion relative to each other during operation is referred to as a flexible multibody structure. Several formulations were proposed for the determination of the dynamic response of controlled flexible multibody structures. In general, these formulations consist of superposing elastic deformations of the component body onto the large rigid body motion of the component. It was shown that this particular methodology for combining linear structural deformations with nonlinear kinematics can lead to erroneous response predictions when either the beam member is very flexible or the rotational speed is high. In addition, previous formulations introduce constraint equations to define the interrelations among system components. This approach increases the number of equations that must be solved, and may result in contraint violation when numerical error accumulates during the integration process. In order to overcome the difficulties, a new approach was suggested. The approach is essentially a finite element formulation which takes advantage of the fact that many multibody structures are joint dominated. The Large Angle Transient Dynamic Analysis (LATDYN) program for clarity of documentation, ease of use, user friendliness, modeling generality, and accuracy of results was evaluated. This required gaining a working familiarity with the code and performing several case studies.

Hernried, Alan G.↗

NASA NTRS1993

Vibration and noise analysis of a gear transmission system

This paper presents a comprehensive procedure to predict both the vibration and noise generated by a gear transmission system under normal operating conditions. The gearbox vibrations were obtained from both numerical simulation and experimental studies using a gear noise test rig. In addition, the noise generated by the gearbox vibrations was recorded during the experimental testing. A numerical method was used to develop linear relationships between the gearbox vibration and the generated noise. The hypercoherence function is introduced to correlate the nonlinear relationship between the fundamental noise frequency and its harmonics. A numerical procedure was developed using both the linear and nonlinear relationships generated from the experimental data to predict noise resulting from the gearbox vibrations. The application of this methodology is demonstrated by comparing the numerical and experimental results from the gear noise test rig.

Choy, F. K.↗

NASA NTRS1993

Vibration and noise analysis of a gear transmission system

Choy, F. K.↗

NASA NTRS1995

On the nonlinear stability of the unsteady, viscous flow of an incompressible fluid in a curved pipe

The stability of the flow of an incompressible, viscous fluid through a pipe of circular cross-section curved about a central axis is investigated in a weakly nonlinear regime. A sinusoidal pressure gradient with zero mean is imposed, acting along the pipe. A WKBJ perturbation solution is constructed, taking into account the need for an inner solution in the vicinity of the outer bend, which is obtained by identifying the saddle point of the Taylor number in the complex plane of the cross-sectional angle co-ordinate. The equation governing the nonlinear evolution of the leading order vortex amplitude is thus determined. The stability analysis of this flow to periodic disturbances leads to a partial differential system dependent on three variables, and since the differential operators in this system are periodic in time, Floquet theory may be applied to reduce this system to a coupled infinite system of ordinary differential equations, together with homogeneous uncoupled boundary conditions. The eigenvalues of this system are calculated numerically to predict a critical Taylor number consistent with the analysis of Papageorgiou. A discussion of how nonlinear effects alter the linear stability analysis is also given, and the nature of the instability determined.

Shortis, Trudi A.↗

NASA NTRS2004

Study and Analyses on the Structural Performance of a Balance

Strain-gauge balances for use in wind tunnels have been designed at Langley Research Center (LaRC) since its inception. Currently Langley has more than 300 balances available for its researchers. A force balance is inherently a critically stressed component due to the requirements of measurement sensitivity. The strain-gauge balances have been used in Langley s wind tunnels for a wide variety of aerodynamic tests, and the designs encompass a large array of sizes, loads, and environmental effects. There are six degrees of freedom that a balance has to measure. The balance s task to measure these six degrees of freedom has introduced challenging work in transducer development technology areas. As the emphasis increases on improving aerodynamic performance of all types of aircraft and spacecraft, the demand for improved balances is at the forefront. Force balance stress analysis and acceptance criteria are under review due to LaRC wind tunnel operational safety requirements. This paper presents some of the analyses and research done at LaRC that influence structural integrity of the balances. The analyses are helpful in understanding the overall behavior of existing balances and can be used in the design of new balances to enhance performance. Initially, a maximum load combination was used for a linear structural analysis. When nonlinear effects were encountered, the analysis was extended to include nonlinearities using MSC.Nastran . Because most of the balances are designed using Pro/Mechanica , it is desirable and efficient to use Pro/Mechanica for stress analysis. However, Pro/Mechanica is limited to linear analysis. Both Pro/Mechanica and MSC.Nastran are used for analyses in the present work. The structural integrity of balances and the possibility of modifying existing balances to enhance structural integrity are investigated.

Karkehabadi, R.↗

NASA NTRS2005

Aeroelastic Sizing for High-Speed Research (HSR) Longitudinal Control Alternatives Project (LCAP)

The Longitudinal Control Alternatives Project (LCAP) compared three high-speed civil transport configurations to determine potential advantages of the three associated longitudinal control concepts. The three aircraft configurations included a conventional configuration with a layout having a horizontal aft tail, a configuration with a forward canard in addition to a horizontal aft tail, and a configuration with only a forward canard. The three configurations were aeroelastically sized and were compared on the basis of operational empty weight (OEW) and longitudinal control characteristics. The sized structure consisted of composite honeycomb sandwich panels on both the wing and the fuselage. Design variables were the core depth of the sandwich and the thicknesses of the composite material which made up the face sheets of the sandwich. Each configuration was sized for minimum structural weight under linear and nonlinear aeroelastic loads subject to strain, buckling, ply-mixture, and subsonic and supersonic flutter constraints. This report describes the methods that were used and the results that were generated for the aeroelastic sizing of the three configurations.

Walsh, Joanne L.↗

NASA NTRS2012

Uncertainty Analysis of the Grazing Flow Impedance Tube

This paper outlines a methodology to identify the measurement uncertainty of NASA Langley s Grazing Flow Impedance Tube (GFIT) over its operating range, and to identify the parameters that most significantly contribute to the acoustic impedance prediction. Two acoustic liners are used for this study. The first is a single-layer, perforate-over-honeycomb liner that is nonlinear with respect to sound pressure level. The second consists of a wire-mesh facesheet and a honeycomb core, and is linear with respect to sound pressure level. These liners allow for evaluation of the effects of measurement uncertainty on impedances educed with linear and nonlinear liners. In general, the measurement uncertainty is observed to be larger for the nonlinear liners, with the largest uncertainty occurring near anti-resonance. A sensitivity analysis of the aerodynamic parameters (Mach number, static temperature, and static pressure) used in the impedance eduction process is also conducted using a Monte-Carlo approach. This sensitivity analysis demonstrates that the impedance eduction process is virtually insensitive to each of these parameters.

Brown, Martha C.↗