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 829 records · Page 46

State and Parameter Estimation for Natural Gas Pipeline Networks Using Transient State Data

Here, we formulate two estimation problems for pipeline systems in which measurements of the compressible gas flowing through a network of pipes are affected by time-varying injections, withdrawals, and compression. We consider a state estimation problem that is then extended to a joint state and parameter estimation problem that can be used for data assimilation. In both formulations, the flow dynamics are described on each pipe by space- and time-dependent densities and mass flux which evolve according to a system of coupled partial differential equations, in which momentum dissipation is modeled using the Darcy-Wiesbach friction approximation. These dynamics are first spatially discretized to obtain a system of nonlinear ordinary differential equations on which state and parameter estimation formulations are given as nonlinear least squares problems. A rapid, scalable computational method for performing a nonlinear least squares estimation is developed. Extensive simulations and computational experiments on multiple pipeline test networks demonstrate the effectiveness of the formulations in obtaining state and parameter estimates in the presence of measurement and process noise.

03 NATURAL GAS↗

Sum-Frequency Generation from a Thin Spherical Layer: I. Analytical Solution

In the tensor form, using the generalized Rayleigh–Gans–Debye approximation, the problem of the sum-frequency generation from a thin nonlinear layer deposited on a dielectric spherical particle placed in a dielectric medium is solved. The second-order nonlinear dielectric susceptibility tensor is chosen in a general form containing chiral components. In the vector and tensor forms, expressions are obtained that describe the spatial distribution of the sum-frequency radiation field generated by two plane electromagnetic elliptically polarized waves. Limiting expressions that describe the spatial distribution of the sum-frequency harmonic at small and large radii of the spherical layer are obtained. It is revealed that, at small radii of the spherical layer, the radiation due to the chiral anisotropy coefficients makes a dominant contribution to the generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Raman amplification of laser pulses near the threshold for plasma wave breaking

Equations are derived for the amplitudes of counter-propagating laser pulses near the threshold for plasma wave breaking, which allow one to describe laser pulses with durations on the order of the plasma oscillation period. In the quasi-monochromatic approximation, they take the form of conventional threewave equations with an additional nonlinearity for the plasma wave. The amplitudes of the amplified laser pulses estimated using these equations agree with results obtained by solving the complete equations. It is shown that Raman amplification of a weak quasi-monochromatic signal (plasma noise) in rarified plasma is significantly suppressed. At the same time, according to numerical simulations, the amplification of laser pulses with durations on the order of the plasma oscillation period is suppressed insignificantly. This result opens new prospects in the application of Raman compression of laser pulses without additional frequency modulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗

Recent Advances of PyROS: A Pyomo Solver for Nonconvex Two-Stage Robust Optimization in Process Systems Engineering

The document presents recent algorithmic and implementation advances of the two-stage robust optimization (RO) solver PyROS, and a benchmarking study which demonstrates the utility of PyROS for two-stage RO problems. The advances include extensions of the scope of PyROS to models with uncertain variable bounds, improvements to the initializations of the subproblems used by the underlying cutting set algorithm, and extensions of the uncertainty set interfaces. The benchmarking study is performed on a library of over 8,500 instances, with variations in the nonlinearities, degree-of-freedom partitioning, uncertainty sets, and polynomial decision rule approximations. An amine-based CO2 capture case study is presented to demonstrate the utility of PyROS for large-scale process models. Overall, the results highlight the effectiveness of PyROS for obtaining robust solutions to optimization problems with uncertain equality constraints.

Sherman, Jason↗

Maximum range three-dimensional lifting planetary entry

Variational equations for maximum range three-dimensional quasisteady glide are given. Nonlinear oscillatory maximum range trajectories obtained with a refined gradient program are approximated by a superposition of quasisteady glide and linearized perturbation equation results. A basic control law is found which is closely followed for maximum cross-range trajectories. The effect of a reradiative heating constraint involving velocity, altitude and angle of attack on a maximum cross-range trajectory for a space shuttle orbiter-type vehicle reentering the earth's atmosphere is investigated numerically.

Dickmanns, E. D.↗

Nonlinear cellular motions in Poiseuille channel flow

In a number of nonlinear solutions to the equations of channel flow the velocity is decomposed into a mean part plus a nonlinear disturbance. The idea that nonlinear effects place a limitation on the amplitudes of the disturbance flow is considered. In the reported investigation the disturbance flow is represented by drastically limited Fourier expansions in the downstream coordinate. The resulting equations are solved numerically with high accuracy to obtain a good representation of the cross-stream structure of the solution. The results of the investigation show that indeed the nonlinear terms always limit the amplitude of the disturbance flow in this approximation.

Zahn, J.-P.↗

Self-ducting of large-amplitude whistler waves

Whistler waves are launched from an electric dipole of length L in a large-volume laboratory plasma. With increasing wave amplitude, the radiation pattern narrows and finally forms a duct of diameter approximately equal to L. The ducted waves propagate nearly undamped. The observed nonlinear effects are explained by wave-particle interactions.

Stenzel, R. L.↗

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The C deg-cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1,000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and asymptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of the nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The continuous cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and symptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗

Numerical methods for stiff systems of two-point boundary value problems

Numerical procedures are developed for constructing asymptotic solutions of certain nonlinear singularly perturbed vector two-point boundary value problems having boundary layers at one or both endpoints. The asymptotic approximations are generated numerically and can either be used as is or to furnish a general purpose two-point boundary value code with an initial approximation and the nonuniform computational mesh needed for such problems. The procedures are applied to a model problem that has multiple solutions and to problems describing the deformation of thin nonlinear elastic beam that is resting on an elastic foundation.

Flaherty, J. E.↗

Theory of wing rock

A theory is developed for predicting wing rock characteristics. From available data, it can be concluded that wing rock is triggered by flow asymmetries, developed by negative or weakly positive roll damping, and sustained by nonlinear aerodynamic roll damping. A new nonlinear aerodynamic model that includes all essential aerodynamic nonlinearities is developed. The Beecham-Titchener method is applied to obtain approximate analytic solutions for the amplitude and frequency of the limit cycle based on the three degree-of-freedom equations of motion. An iterative scheme is developed to calculate the average aerodynamic derivatives and dynamic characteristics at limit cycle conditions. Good agreement between theoretical and experimental results is obtained.

Hsu, C. H.↗

Down-to-earth supercomputing at NASA Ames

The requirements for computers for aerospace vehicle design problems, in particular fluid dynamics, are examined. The levels of approximation to Navier-Stokes equations are: (1) linearized inviscid form, (2) nonlinear inviscid equation, (3) Reynolds-averaged Navier-Stokes equation, (4) large eddy simulation with models for subgrid scale turbulence, and (5) full Navier-Stokes equations to calculate laminar/turbulent transition and turbulence dissipation. The Cray 2 with 256 megaword memory and 4.1 nsec clock and the Cray X-MP/48 with a 128 megaword memory and a 9.5 nsec clock of the Ames facility are described, and an example revealing their applicability is presented. The use of zonal methods to simulate aerodynamic flow fields numerically is discussed. Consideration is given to the development of a program for simulation of aerodynamic flow fields, the compilers, symbolic debuggers, the scheduling of multiple CPUs, and the analysis and storage of the data.

Stevens, Kenneth G., Jr.↗

Estimation of time- and state-dependent delays and other parameters in functional differential equations

A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.

Murphy, K. A.↗

Estimation of time- and state-dependent delays and other parameters in functional differential equations

A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.

Murphy, K. A.↗