SEARCH · Search NASA

Results for “Nonlinear approximation”

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,621 records · Page 90

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Numerical simulation of transition in wall-bounded shear flows

The current status of numerical simulation techniques for the transition to turbulence in incompressible channel and boundary-layer flows is surveyed, and typical results are presented graphically. The focus is on direct numerical simulations based on the full nonlinear time-dependent Navier-Stokes equations without empirical closure assumptions for prescribed initial and boundary conditions. Topics addressed include the vibrating ribbon problem, space and time discretization, initial and boundary conditions, alternative methods based on the triple-deck approximation, two-dimensional channel and boundary-layer flows, three-dimensional boundary layers, wave packets and turbulent spots, compressible flows, transition control, and transition modeling.

Kleiser, Leonhard↗

NASA NTRS1990

Short wavelength striations on expanding plasma clouds

The present evaluation of current understanding of the growth and evolution of less-than-1 ion gyroradius 'flute modes' on a plasma as it expands across and ambient magnetic field notes that the mechanism by which the instability is generated, and its approximate linear theory (encompassing nonlocal, finite-beta, and collisional effects), have reached a satisfactory degree of development. AMPTE Ba releases have been the bases of most of the observational studies. Substantial progress is also noted in the development of a nonlinear mode-coupling theory which can resolve remaining differences between theory and observation.

Winske, D.↗

NASA NTRS1989

Globally convergent techniques in nonlinear Newton-Krylov

Some convergence theory is presented for nonlinear Krylov subspace methods. The basic idea of these methods is to use variants of Newton's iteration in conjunction with a Krylov subspace method for solving the Jacobian linear systems. These methods are variants of inexact Newton methods where the approximate Newton direction is taken from a subspace of small dimensions. The main focus is to analyze these methods when they are combined with global strategies such as linesearch techniques and model trust region algorithms. Most of the convergence results are formulated for projection onto general subspaces rather than just Krylov subspaces.

Brown, Peter N.↗

NASA NTRS1992

Comparison between the performance of two classes of fuzzy controllers

This paper presents an application comparison between two classes of fuzzy controllers: the Clearness Transformation Fuzzy Controller (CTFC) and the CRI-based Fuzzy Controller. The comparison is performed by studying the application of the controllers to simulation examples of nonlinear systems. The CTFC is a new approach for the organization of fuzzy controllers based on a cognitive model of parameter driven control, the notion of fuzzy patterns to represent fuzzy knowledge and the Clearness Transformation Rule of Inference (CTRI) for approximate reasoning. The approach facilitates the implementation of the basic modules of the controller: the fuzzifier, defuzzifier, and the control protocol in a rule-based architecture. The CTRI scheme for approximate reasoning does not require the formation of fuzzy relation matrices yielding improved performance in comparison with the traditional organization of fuzzy controllers.

Janabi, T. H.↗

NASA NTRS1994

Factors controlling precision and accuracy in isotope-ratio-monitoring mass spectrometry

The performance of systems in which picomole quantities of sample are mixed with a carrier gas and passed through an isotope-ratio mass spectrometer system was examined experimentally and theoretically. Two different mass spectrometers were used, both having electron-impact ion sources and Faraday cup collector systems. One had an accelerating potential of 10kV and accepted 0.2 mL of He/min, producing, under those conditions, a maximum efficiency of 1 CO2 molecular ion collected per 700 molecules introduced. Comparable figures for the second instrument were 3 kV, 0.5 mL of He/min, and 14000 molecules/ion. Signal pathways were adjusted so that response times were <200 ms. Sample-related ion currents appeared as peaks with widths of 3-30 s. Isotope ratios were determined by comparison to signals produced by standard gases. In spite of rapid variations in signals, observed levels of performance were within a factor of 2 of shot-noise limits. For the 10-kV instrument, sample requirements for standard deviations of 0.1 and 0.5% were 45 and 1.7 pmol, respectively. Comparable requirements for the 3-kV instrument were 900 and 36 pmol. Drifts in instrumental characteristics were adequately neutralized when standards were observed at 20-min intervals. For the 10-kV instrument, computed isotopic compositions were independent of sample size and signal strength over the ranges examined. Nonlinearities of <0.04%/V were observed for the 3-kV system. Procedures for observation and subtraction of background ion currents were examined experimentally and theoretically. For sample/ background ratios varying from >10 to 0.3, precision is expected and observed to decrease approximately 2-fold and to depend only weakly on the precision with which background ion currents have been measured.

NASA Discipline Exobiology↗

NASA NTRSDate not supplied

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗

DOE OSTI2025

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Approximate Bayesian Model Inversion for PDEs with Heterogeneous and State-Dependent Coefficients

We present two approximate Bayesian inference methods for parameter estimation in partial differential equation (PDE) models with space-dependent and state-dependent parameters. We demonstrate that these methods provide accurate and cost-effective alternatives to Markov Chain Monte Carlo simulation. We assume a parameterized Gaussian prior on the unknown functions, and approximate the posterior density by a parameterized multivariate Gaussian density. The parameters of the prior and posterior are estimated from sparse observations of the PDE model's states and the unknown functions themselves by maximizing the evidence lower bound (ELBO), a lower bound on the log marginal likelihood of the observations. The first method, Laplace-EM, employs the expectation maximization algorithm to maximize the ELBO, with a Laplace approximation of the posterior on the E-step, and minimization of a Kullback-Leibler divergence on the M-step. The second method, DSVI-EB, employs the doubly stochastic variational inference (DSVI) algorithm, in which the ELBO is maximized via gradient-based stochastic optimization, with nosiy gradients computed via simple Monte Carlo sampling and Gaussian backpropagation. We apply these methods to identifying diffusion coefficients in linear and nonlinear diffusion equations, and we find that both methods provide accurate estimates of posterior densities and the hyperparameters of Gaussian priors. While the Laplace-EM method is more accurate, it requires computing Hessians of the physics model. The DSVI-EB method is found to be less accurate but only requires gradients of the physics model.

approximate Bayesian inference, model inversion, v↗

DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Stellarator microinstabilities and turbulence at low magnetic shear

Gyrokinetic simulations of drift waves in low-magnetic-shear stellarators reveal that simulation domains comprised of multiple turns can be required to properly resolve critical mode structures important in saturation dynamics. Marginally stable eigenmodes important in saturation of ion temperature gradient modes and trapped electron modes in the Helically Symmetric Experiment (HSX) stellarator are observed to have two scales, with the envelope scale determined by the properties of the local magnetic shear and an inner scale determined by the interplay between the local shear and magnetic field-line curvature. Properly resolving these modes removes spurious growth rates that arise for extended modes in zero-magnetic-shear approximations, enabling use of a zero-magnetic-shear technique with smaller simulation domains and attendant cost savings. Analysis of subdominant modes in trapped electron mode (TEM)-driven turbulence reveals that the extended marginally stable modes play an important role in the nonlinear dynamics, and suggests that the properties induced by low magnetic shear may be exploited to provide another route for turbulence saturation.

Faber, B. J. (ORCID:000000034934400X)↗

DOE OSTI2023

Black and gray box learning of amplitude equations: Application to phase field systems

Here, we present a data-driven approach to learning surrogate models for amplitude equations and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher-order eikonal equation and its approximation, the Kardar-Parisi-Zhang equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black box identification, we explore different approaches to learning data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

42 ENGINEERING↗

NASA NTRS1974

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement distributions in various bodies containing cracks. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The results obtained show a considerable potential for using this method in the three-dimensional analysis of finite geometry solids and suggest a possible extension of this technique to nonlinear material behavior.

Gyekenyesi, J. P.↗

NASA NTRS1975

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

Gyekenyesi, J. P.↗

NASA NTRS1981

Helicopter rotor loads using a matched asymptotic expansion technique

The theoretical basis and computational feasibility of the Van Holten method, and its performance and range of validity by comparison with experiment and other approximate methods was examined. It is found that within the restrictions of incompressible, potential flow and the assumption of small disturbances, the method does lead to a valid description of the flow. However, the method begins to break down under conditions favoring nonlinear effects such as wake distortion and blade/rotor interaction.

Pierce, G. A.↗

NASA NTRS1983

The peculiar velocity field in flattened superclusters

The linear growth rate of small transverse perturbations in a self-gravitating, collisionless gas that is undergoing a one-dimensional collapse is strongly influenced by the nonlinear flow of the background. In the collapse plane, smaller Hubble flow deviations and infall velocities are obtained than in the simple linear theory. Application of this theory to the kinematics of galaxies in a flat supercluster shows that the usual linear approximation may seriously underestimate Omega, the density parameter of the universe, if there is a strong deviation from spherical symmetry. Furthermore, dissipative separation of baryonic matter from collisionless dark material (such as hypothesized massive neutrinos) enhances this effect, leading to apparent Omega values of 0.2-0.3 for a wide range of separation parameters in an Einstein-de Sitter (Omega = 1) model.

Szalay, A. S.↗

NASA NTRS1990

Agent independent task planning

Agent-Independent Planning is a technique that allows the construction of activity plans without regard to the agent that will perform them. Once generated, a plan is then validated and translated into instructions for a particular agent, whether a robot, crewmember, or software-based control system. Because Space Station Freedom (SSF) is planned for orbital operations for approximately thirty years, it will almost certainly experience numerous enhancements and upgrades, including upgrades in robotic manipulators. Agent-Independent Planning provides the capability to construct plans for SSF operations, independent of specific robotic systems, by combining techniques of object oriented modeling, nonlinear planning and temporal logic. Since a plan is validated using the physical and functional models of a particular agent, new robotic systems can be developed and integrated with existing operations in a robust manner. This technique also provides the capability to generate plans for crewmembers with varying skill levels, and later apply these same plans to more sophisticated robotic manipulators made available by evolutions in technology.

Davis, William S.↗

NASA NTRS1992

Minimax confidence intervals in geomagnetism

The present paper uses theory of Donoho (1989) to find lower bounds on the lengths of optimally short fixed-length confidence intervals (minimax confidence intervals) for Gauss coefficients of the field of degree 1-12 using the heat flow constraint. The bounds on optimal minimax intervals are about 40 percent shorter than Backus' intervals: no procedure for producing fixed-length confidence intervals, linear or nonlinear, can give intervals shorter than about 60 percent the length of Backus' in this problem. While both methods rigorously account for the fact that core field models are infinite-dimensional, the application of the techniques to the geomagnetic problem involves approximations and counterfactual assumptions about the data errors, and so these results are likely to be extremely optimistic estimates of the actual uncertainty in Gauss coefficients.

Stark, Philip B.↗

NASA NTRS2011

Compression Behavior of Fluted-Core Composite Panels

In recent years, fiber-reinforced composites have become more accepted for aerospace applications. Specifically, during NASA s recent efforts to develop new launch vehicles, composite materials were considered and baselined for a number of structures. Because of mass and stiffness requirements, sandwich composites are often selected for many applications. However, there are a number of manufacturing and in-service concerns associated with traditional honeycomb-core sandwich composites that in certain instances may be alleviated through the use of other core materials or construction methods. Fluted-core, which consists of integral angled web members with structural radius fillers spaced between laminate face sheets, is one such construction alternative and is considered herein. Two different fluted-core designs were considered: a subscale design and a full-scale design sized for a heavy-lift-launch-vehicle interstage. In particular, axial compression of fluted-core composites was evaluated with experiments and finite-element analyses (FEA); axial compression is the primary loading condition in dry launch-vehicle barrel sections. Detailed finite-element models were developed to represent all components of the fluted-core construction, and geometrically nonlinear analyses were conducted to predict both buckling and material failures. Good agreement was obtained between test data and analyses, for both local buckling and ultimate material failure. Though the local buckling events are not catastrophic, the resulting deformations contribute to material failures. Consequently, an important observation is that the material failure loads and modes would not be captured by either linear analyses or nonlinear smeared-shell analyses. Compression-after-impact (CAI) performance of fluted core composites was also investigated by experimentally testing samples impacted with 6 ft.-lb. impact energies. It was found that such impacts reduced the ultimate load carrying capability by approximately 40% on the subscale test articles and by less than 20% on the full-scale test articles. Nondestructive inspection of the damage zones indicated that the detectable damage was limited to no more than one flute on either side of any given impact. More study is needed, but this may indicate that an inherent damage-arrest capability of fluted core could provide benefits over traditional sandwich designs in certain weight-critical applications.

Schultz, Marc R.↗