Search NASA⌕ Search

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.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Total Variation Diminishing (TVD) schemes of uniform accuracy

Explicit second-order accurate finite-difference schemes for the approximation of hyperbolic conservation laws are presented. These schemes are nonlinear even for the constant coefficient case. They are based on first-order upwind schemes. Their accuracy is enhanced by locally replacing the first-order one-sided differences with either second-order one-sided differences or central differences or a blend thereof. The appropriate local difference stencils are selected such that they give TVD schemes of uniform second-order accuracy in the scalar, or linear systems, case. Like conventional TVD schemes, the new schemes avoid a Gibbs phenomenon at discontinuities of the solution, but they do not switch back to first-order accuracy, in the sense of truncation error, at extrema of the solution. The performance of the new schemes is demonstrated in several numerical tests.

Hartwich, PETER-M.↗
NASA NTRS1990

Structure and evolution of small-amplitude intermediate shock waves

A simplified set of equations is derived that approximates the magnetohydrodynamic (MHD) Navier-Stokes equations for weakly nonlinear disturbances whose speeds are close to the MHD intermediate speed. Its shock structure solutions are then examined. The fast and slow shock solutions are uniquely specified by their Rankine-Hugoniot relations. However, the intermediate shock solutions, which are not unique, are characterized by the integral through the shock of the noncoplanar component of the magnetic field. For situations in which this integral is conserved, the Riemann problem is well defined and predicts the evolution of intermediate shocks. This analysis is substantiated by numerical computations.

Kennel, C. F.↗
NASA NTRS1991

Implicit solvers for unstructured meshes

Implicit methods were developed and tested for unstructured mesh computations. The approximate system which arises from the Newton linearization of the nonlinear evolution operator is solved by using the preconditioned GMRES (Generalized Minimum Residual) technique. Three different preconditioners were studied, namely, the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over relaxation (SSOR). The preconditioners were optimized to have good vectorization properties. SSOR and ILU were also studied as iterative schemes. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also studied. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.

Venkatakrishnan, V.↗
NASA NTRS1992

A model for nonlinear collisionless reconnection

The time dependence and scaling of the electric and magnetic field components during the nonlinear phase of reconnection of 'thin' unstable neutral layers (thickness approximately equals electron gyroradius) are investigated. Solutions to model equations show an initial interval of rapid growth and subsequent large-amplitude fluctuations of the trapped magnetic flux and the ion kinetic energy. The fluctuations arise from the transient formation, merging, and dissolution of magnetic islands.

Lovelace, R. V. E.↗
NASA NTRS1992

Rapid time variations and radiative instabilities of astrophysical masers

The time-varying intensities are obtained for astrophysical masers that are radiatively unstable. Numerical integrations of the time-dependent, nonlinear equations of radiative transfer are performed with the usual approximation of a linear maser. At long times after changes in the physical conditions, the intensity of maser radiation reaches an asymptotic behavior and oscillates permanently in these idealized calculations with a period that is related to the length of the maser divided by the speed of light. The intensity varies by more than a factor of 10. These intensities depend upon the same four parameters as we originally found to determine the regime for radiative instabilities based on a stability analysis of the steady state. A detailed comparison is made between the predictions of the stability analysis and the time variations of the intensities. Calculations are performed for interacting pairs as well as isolated, individual masers.

Scappaticci, Gerardo A.↗
NASA NTRS1992

Theoretical and computational study of rotating isotropic turbulence

The various flow regimes that occur in rotating isotropic turbulence and their associated turbulence structure are studied. Direct numerical simulations are conducted for moderate and rapid rotation rates and comparisons are made with the predictions of generalized Eddy Damped Quasi-Normal Markovian (EDQNM) approximations. It is shown that at high rotation rates the nonlinear transfer terms remain small and that the development of the spectrum is through pure viscous decay. This shows that the effect of rapid rotation on the turbulence kinetic energy is through the shutting off of the production term in the dissipation rate equation. At moderate and relatively rapid rates, the rotation causes a discernible anisotropy to develop in the integral length scale with a mild trend toward a two-dimensionalization of the flow. As the flow decays, the Rossby number decreases, leading again to the shutting off of the nonlinear transfer terms. For extremely rapid rotation rates, the anisotropies in the integral length scales are small and direct numerical simulation results are in good agreement with rapid distortion theory.

Mansour, N. N.↗
NASA NTRS1993

Nonlinear evolution of the first mode supersonic oblique waves in compressible boundary layers. Part 1: Heated/cooled walls

The nonlinear stability of an oblique mode propagating in a two-dimensional compressible boundary layer is considered under the long wave-length approximation. The growth rate of the wave is assumed to be small so that the concept of unsteady nonlinear critical layers can be used. It is shown that the spatial/temporal evolution of the mode is governed by a pair of coupled unsteady nonlinear equations for the disturbance vorticity and density. Expressions for the linear growth rate show clearly the effects of wall heating and cooling and in particular how heating destabilizes the boundary layer for these long wavelength inviscid modes at O(1) Mach numbers. A generalized expression for the linear growth rate is obtained and is shown to compare very well for a range of frequencies and wave-angles at moderate Mach numbers with full numerical solutions of the linear stability problem. The numerical solution of the nonlinear unsteady critical layer problem using a novel method based on Fourier decomposition and Chebychev collocation is discussed and some results are presented.

Gajjar, J. S. B.↗
NASA NTRS1984

Two dimensional nonlinear analysis of sound transmission through a near-sonic throat flow

A two-dimensional nonlinear theory of sound transmission through a nonuniform duct carrying a near-sonic throat flow is described. The mean flow in the duct is treated according to a generalized quasi-one dimensional model, and the analysis of the unsteady perturbations is carried out using the Method of Matched Asymptotic Expansions. The linearized acoustic field in the subsonic regions of the duct is approximated by a Wave Envelope model which is matched asymptotically with a nonlinear inner solution valid in the near-sonic throat region. Numerical results are presented to illustrate the predictions of the theory. It is found that, in general, shock waves develop in the acoustic field during transmission through the near-sonic flow. Unlike previous one-dimensional theories, the current study shows that dispersion can play a major role in the propagation process.

Myers, M. K.↗
NASA NTRS1999

An Analysis of the Macroscopic Tensile Behavior of a Nonlinear Nylon Reinforced Elastomeric Composite System Using MAC/GMC

A special class of composite laminates composed of soft rubbery matrices and stiff reinforcements made of steel wires or synthetic fibers is examined, where each constituent behaves in a nonlinear fashion even in the small strain domain. Composite laminates made of piles stacked at alternating small orientation angles with respect to the applied axial strain are primarily dominated by the nonlinear behavior of the reinforcing fibers. However; composites with large ply orientations or those perpendicular to the loading axis, will approximate the behavior of the matrix phase and respond in even a more complex fashion for arbitrarily stacked piles. The geometric nonlinearity due to small cord rotations during loading was deemed here to have a second order effect and consequently dropped from any consideration. The user subroutine USRMAT within the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC), was utilized to introduce the constituent material nonlinear behavior. Stress-strain behavior at the macro level was experimentally generated for single and multi ply composites comprised of continuous Nylon-66 reinforcements embedded in a carbon black loaded rubbery matrix. Comparisons between the predicted macro composite behavior and experimental results are excellent when material nonlinearity is included in the analysis. In this paper, a brief review of GMC is provided, along with a description of the nonlinear behavior of the constituents and associated constituent constitutive relations, and the improved macro (or composite) behavior predictions are documented and illustrated.

Assaad, Mahmoud↗
DOE OSTI2017

Multipolar second harmonic generation in a symmetric nonlinear metamaterial

Optical nonlinearities are intimately related to the spatial symmetry of the nonlinear media. For example, the second order susceptibility vanishes for centrosymmetric materials under the dipole approximation. The latter concept has been naturally extended to the metamaterials’ realm, sometimes leading to the (erroneous) hypothesis that second harmonic (SH) generation is negligible in highly symmetric meta-atoms. In this work we aim to show that such symmetric meta-atoms can radiate SH light efficiently. In particular, we investigate in-plane centrosymmetric meta-atom designs where the approximation for meta-atoms breaks down. In a periodic array this building block allows us to control the directionality of the SH radiation. We conclude by showing that the use of symmetry considerations alone allows for the manipulation of the nonlinear multipolar response of a meta-atom, resulting in e.g. dipolar, quadrupolar, or multipolar emission on demand. This is because the size of the meta-atom is comparable with the free-space wavelength, thus invalidating the dipolar approximation for meta-atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Dynamical systems and nonlinear transient rheology of the far-from-equilibrium Bjorken flow

In relativistic kinetic theory, the one-particle distribution function is approximated by an asymptotic perturbative power series in the Knudsen number which is divergent. For the Bjorken flow, we expand the distribution function in terms of its moments and study their nonlinear evolution equations. The resulting coupled dynamical system can be solved for each moment consistently using a multiparameter transseries which makes the constitutive relations inherit the same structure. A new nonperturbative dynamical renormalization scheme is born out of this formalism that goes beyond the linear response theory.We show that there is a Lyapunov function, also known as dynamical potential, which is, in general, a function of the moments and time satisfying Lyapunov stability conditions along renormalization group flows connected to the asymptotic hydrodynamic fixed point. As a result, the transport coefficients get dynamically renormalized at every order in the time-dependent perturbative expansion by receiving nonperturbative corrections present in the transseries. The connection between the integration constants and the UV data is discussed using the language of dynamical systems. Furthermore, we show that the first dissipative correction in the Knudsen number to the distribution function is not only determined by the known effective shear viscous term but also a new high-energy nonhydrodynamic mode. It is demonstrated that the survival of this new mode is intrinsically related to the nonlinear mode-to-mode coupling with the shear viscous term. Finally, we comment on some possible phenomenological applications of the proposed nonhydrodynamic transport theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1986

Generalized energetics for linear and nonlinear symmetric instabilities

The generalized energy conservation integral is applied to linear and nonlinear symmetric instabilities (SIs) and conditional symmetric instabilities (CSIs) in uniform and nonuniform basic states. The generalized energy conservation integral for linear SI in nonuniform basic state is refined to include initial thermal-inertial perturbations. A generalized conservation integral is derived for a nonlinear SI with uniform and nonuniform basic states. The integral reveals that: (1) the stability conditions of nonlinear and linear SI are not the same; (2) when the basic state is symmetrically stable the nonlinear evolution of a symmetric disturbance is energetically bounded by its initial energy; and (3) when the basic state is symmetrically unstable the growth of the symmetric perturbation is approximated by the linear SI theory. Precipitative heating is incorporated into the generalized energy integral for nonlinear CSI; the relation between nonlinear evolution of a CSI circulation and precipitative heating is examined.

Xu, Q.↗
NASA NTRS1986

A methodology for airplane parameter estimation and confidence interval determination in nonlinear estimation problems

An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. With the fitted surface, sensitivity information can be updated at each iteration with less computational effort than that required by either a finite-difference method or integration of the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, and thus provides flexibility to use model equations in any convenient format. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. The degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels and to predict the degree of agreement between CR bounds and search estimates.

Murphy, P. C.↗
DOE OSTI2020

Nonlinear manifold-based reduced order model

Traditional linear subspace reduced order models (LS-ROMs) are able to accelerate physical simulations, in which the intrinsic solution space falls into a subspace with a small dimension, i.e., the solution space has a small Kolmogorov n-width. However, for physical phenomena not of this type, e.g., any advection-dominated flow phenomena such as in traffic flow, atmospheric flows, and air flow over vehicles, a low-dimensional linear subspace poorly approximates the solution. To address cases such as these, we have developed a fast and accurate physicsinformed neural network ROM, namely nonlinear manifold ROM (NM-ROM), which can better approximate high-fidelity model solutions with a smaller latent space dimension than the LSROMs. Our method takes advantage of the existing numerical methods that are used to solve the corresponding full order models. The efficiency is achieved by developing a hyper-reduction technique in the context of the NM-ROM. Numerical results show that neural networks can learn a more efficient latent space representation on advection-dominated data from 1D and 2D Burgers’ equations. A speedup of up to 2.6 for 1D Burgers’ and a speedup of 11.7 for 2D Burgers’ equations are achieved with an appropriate treatment of the nonlinear terms through a hyper-reduction technique. Finally, a posteriori error bounds for the NM-ROMs are derived that take account of the hyper-reduced operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Solution and sensitivity analysis of nonlinear equations using a hypercomplex-variable Newton-Raphson method

Here, the classical Newton-Raphson (NR) method for solving nonlinear equations is enhanced in two ways through the use of hypercomplex variables and algebra. In particular, i) the Jacobian is computed in a highly accurate and automated way, and ii) the derivative of the solution to the nonlinear equations is computed with respect to any parameter contained within the system of equations. These advances provide two significant enhancements in that it is straightforward to provide an accurate Jacobian and to construct a reduced order model (ROM) of arbitrary order with respect to any parameter of the system. The ROM can then be used to approximate the solution for other parameter values without requiring additional solutions of the nonlinear equations. Several case studies are presented including 1D and 2D academic examples with fully functioning Python code provided. Additionally, a case of study of the catenary of an elastic cable subject to its own weight and a vertical point load. Derivatives up to 10th order were computed with respect to material, loading, and geometrical parameters. The derivatives were used to generate reduced order models of the cable deformation and reaction forces at its ends with respect to multiple input parameters. Results show that from a single hypercomplex evaluation of the cable under a single vertical point load, it is possible to generate an accurate reduced order model capable of predicting the cable deformation with 1.5 times the load in the opposite direction and with 3.5 times the load in the same direction without resolving the system of equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Combined strength and aeroelastic wing synthesis via constraint approximation

A numerical optimization program using constraint approximations for preliminary sizing of wing structural design parameters to satisfy simultaneous strength and aeroelastic requirements is described. The effects of wing flexibility on loads and flutter are included. The iterative procedure approximates wing internal load distributions and flutter response in generating design constraints. Though a linear Taylor series approximation to the variation of flutter speed is made, constraints in general are formulated as nonlinear and solved in a minimum weight problem using a feasible directions search. The technique is demonstrated for the determination of the optimal ply orientations and for simultaneous sizing for strength and flutter with simple constraint approximations.

Mullen, J., Jr.↗