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,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗

DOE OSTI2025

Causal machine learning uncovers conditions for convective intensification driven by organic and sulfate aerosols

Aerosols are often hypothesized to invigorate deep convective clouds (DCCs), but observational evidence remains limited and inconclusive. Clarifying this hypothesis is critical for regions vulnerable to thunderstorms and flooding, particularly highly polluted coastal cities. Leveraging a novel causal discovery–inference pipeline and high-resolution observations near Houston, TX, we identify multiple causal pathways among aerosols (mostly organic and sulfate), DCCs, and meteorological factors. However, a direct causal link from aerosols to DCCs is found to be uncommon, occurring in less than 35% of analyzed scenarios, and is characterized by strong conditionality and nonlinearity. When aerosol impacts on DCCs do occur, they can be substantial, enhancing DCC core heights by approximately 1.7 km, with 92% of this effect concentrated in warmer-phase cloud regions. Notably, the presence of sea breezes and the inclusion of all measured aerosol particles each enhance DCCs in over 95% of aerosol-sensitive cases.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Instability of an electron-plasma shear layer in an externally imposed strain flow

The E x B shear instability of a two-dimensional (2D) filament (i.e., a thin, rectangular strip perpendicular to the magnetic field) of magnetized pure electron plasma is investigated experimentally in the presence of an externally imposed strain flow. Data are acquired using a specialized Penning–Malmberg trap in which strain flows can be applied in 2D by biasing segmented electrodes surrounding the plasma. The E x B drift dynamics are well-described by the Drift-Poisson equations, which are isomorphic to the 2D Euler equations describing ideal fluids. Thus, the experimental results correspond to the Rayleigh instability of a shear layer in a 2D ideal fluid, where the electron density is analogous to the fluid vorticity. Shear layers are prepared by stretching initially axisymmetric electron vortices using a strong, applied strain flow. The data at early times are in quantitative agreement with a linear model which extends Rayleigh’s work to account for the influence of an external strain flow. In the presence of weak strain, the system approximately maintains a phase relationship that corresponds to an instantaneous Rayleigh eigenmode. The instability develops into the nonlinear regime later in time and at smaller spatial scales as the strain rate is increased. A secondary vortex pairing instability is observed, but it is suppressed when the strain-to-vorticity ratio exceeds roughly 0.025. In this way, vorticity transport perpendicular to the filament is diminished due to the applied strain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Iterative Linearization for Phasor-Defined Optimal Power Dispatch

Optimal power flow (OPF) problems, which dispatch power targets to controllable generating units across a network, must generally account for non-convex constraints on power flow. Furthermore, adapting those problems so as to make them solvable with convex optimization techniques is an area of much academic and operational interest. In this paper, we present a method for solving OPF as a quadratic program by iteratively refining and re-initializing a linearized model of power flow based on the outputs of an associated nonlinear solver. The linear model on which we demonstrate this method is an adapted version of an approximation designed for use with unbalanced distribution networks. As an important benefit, the model allows for the explicit inclusion of nodal voltage phasor values in both the OPF problem's objective and its constraints, which opens the door to the idea of phasor-based control (PBC) design. We show in simulations on the IEEE 13-node test feeder that our method quickly converges to a set of phasor targets that are sufficiently precise for use in operations at the distribution level.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1973

Phase matching in second-harmonic generation using artificial periodic structures.

The use of artificial periodic structures, consisting of spatial modulations of the linear and nonlinear susceptibilities of a nonlinear optical medium, to achieve phase matching in second-harmonic generation is analyzed. Dispersion relations and approximate formulas for the second-harmonic fields generated under various conditions are obtained and used to evaluate the experimental situation.

Tang, C. L.↗

NASA NTRS1974

Error estimates for approximate dynamic systems

Two approaches are investigated for obtaining estimates on the error between approximate and exact solutions of dynamic systems. The first method is primarily useful if the system is nonlinear and of low dimension. The second requires construction of a system of v-functions but is useful for higher dimensional systems, either linear or nonlinear.

Gunderson, R. W.↗

NASA NTRS1980

Output feedback non-linear decoupled control synthesis and observer design for manoeuvring aircraft

A study of the applicability of nonlinear decoupling theory to the design of control systems using output feedback for maneuvering aircraft is presented. The response variables chosen for decoupled control were angular velocity components along roll, pitch, and yaw axes, angle of attack (p), and angle of sideslip, using aileron, rudder, and elevator controls. An observer design for a class of nonlinear systems was presented and this method was used to estimate angle of attack and sideslip; an approximate observer was obtained by neglecting derivatives of p and aileron deflection angles and it was used in a simulation study. A simulation study showed that precise rapid combined lateral and longitudinal maneuvers can be performed; it was also demonstrated that a bank-angle-command outer loop could be designed for precise bank angles changes and simultaneous large lift maneuvers.

Singh, S. N.↗

NASA NTRS1982

Geometric and rediametric distortion in spaceborne SAR imagery

Distortions inherent on synthetic aperture radio (SAR) imagery and the development to date of unsupervised postprocessing rectification techniques are described. The geometric distortion can be divided into two categories: (1) distortion derived from the radar viewing geometry, this includes such effects as ground range nonlinearities, radar foreshortening and radar layover; (2) distortion introduced during the data processing, these distortions result from approximations made during the correlation such as in estimation of the target phase history, or compensation for the earth rotation. The processor induced distortions depends on the specific correlation algorithm used for image formation. The effects are addressed on the image product resulting from assumptions during the processing and it specifically considers distortions inherent in digital imagery produced by the digital image processor.

Curlander, J. C.↗

NASA NTRS1983

Numerical computations on one-dimensional inverse scattering problems

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

Dunn, M. H.↗

NASA NTRS1984

Nonlinear aerodynamic effects on bodies in supersonic flow

The supersonic flow about generic bodies was analyzed to identify the elments of the nonlinear flow and to determine the influence of geometry and flow conditions on the magnitude of these nonlinearities. The nonlinear effects were attributed to separated-flow nonlinearities and attached-flow nonlinearities. The nonlinear attached-flow contribution was further broken down into large-disturbance effects and entropy effects. Conical, attached-flow bundaries were developed to illustrate the flow regimes where the nonlinear effects are significant, and the use of these boundaries for angle of attack and three-dimensional geometries was indicated. Normal-force and pressure comparisons showed that the large-disturbance and separated-flow effects were the dominant nonlinear effects at low supersonic Mach numbers and that the entropy effects were dominant for high supersonic Mach number flow. The magnitude of all the nonlinear effects increased with increasing angle of attack. A full-potential method, NCOREL, which includes an approximate entropy correction, was shown to provide accurate attached-flow pressure estimates from Mach 1.6 through 4.6.

Pittman, J. L.↗

NASA NTRS1983

A large-signal dynamic simulation for the series resonant converter

A simple nonlinear discrete-time dynamic model for the series resonant dc-dc converter is derived using approximations appropriate to most power converters. This model is useful for the dynamic simulation of a series resonant converter using only a desktop calculator. The model is compared with a laboratory converter for a large transient event.

King, R. J.↗

NASA NTRS1984

Numerical computations on one-dimensional inverse scattering problems

Dunn, M. H.↗

NASA NTRS1986

A comparison of the bounded derivative and the normal-mode initialization methods using real data

Application of the bounded-derivative and normal-mode methods to a simple linear barotropic model at a typical middle latitude shows that the two methods lead to identical constraints up to a certain degree of approximation. Beyond this accuracy the two methods may differ from each other. When applied to a global nonlinear barotropic model using real data, again the two methods lead to similar balanced initial states. The gravity oscillations in the unbalanced height field, which have amplitudes of up to 60 m with a dominant periodicity of about 5 to 6 h, are practically eliminated by both initialization methods. The rotational wind component is smooth even for the unbalanced initial state. The small-scale spatial features of the irrotational wind component are drastically reduced by initialization. Both the nonlinear normal-mode and the bounded-derivative initialization methods yield similar divergence fields centered around the areas of highest orography. The comparison shows that there is no significant loss of information in the mass and momentum fields, despite the fact that the bounded-derivative method employs only the original, rotational wind component to construct a balanced initial state compared to the normal-mode method, which, in addition, makes use of the unbalanced divergent wind and height fields.

Semazzi, F. H. M.↗

NASA NTRS1989

Development of an integrated BEM approach for hot fluid structure interaction

The progress made toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-Orbit engine hot section components is reported. The convective viscous integral formulation was derived and implemented in the general purpose computer program GP-BEST. The new convective kernel functions, in turn, necessitated the development of refined integration techniques. As a result, however, since the physics of the problem is embedded in these kernels, boundary element solutions can now be obtained at very high Reynolds number. Flow around obstacles can be solved approximately with an efficient linearized boundary-only analysis or, more exactly, by including all of the nonlinearities present in the neighborhood of the obstacle. The other major accomplishment was the development of a comprehensive fluid-structure interaction capability within GP-BEST. This new facility is implemented in a completely general manner, so that quite arbitrary geometry, material properties and boundary conditions may be specified. Thus, a single analysis code (GP-BEST) can be used to run structures-only problems, fluids-only problems, or the combined fluid-structure problem. In all three cases, steady or transient conditions can be selected, with or without thermal effects. Nonlinear analyses can be solved via direct iteration or by employing a modified Newton-Raphson approach.

Dargush, G. F.↗

NASA NTRS1994

Fracture mechanics validity limits

Fracture behavior is characteristics of a dramatic loss of strength compared to elastic deformation behavior. Fracture parameters have been developed and exhibit a range within which each is valid for predicting growth. Each is limited by the assumptions made in its development: all are defined within a specific context. For example, the stress intensity parameters, K, and the crack driving force, G, are derived using an assumption of linear elasticity. To use K or G, the zone of plasticity must be small as compared to the physical dimensions of the object being loaded. This insures an elastic response, and in this context, K and G will work well. Rice's J-integral has been used beyond the limits imposed on K and G. J requires an assumption of nonlinear elasticity, which is not characteristic of real material behavior, but is thought to be a reasonable approximation if unloading is kept to a minimum. As well, the constraint cannot change dramatically (typically, the crack extension is limited to ten-percent of the initial remaining ligament length). Rice, et al investigated the properties required of J-type parameters, J(sub x), and showed that the time rate, dJ(sub x)/dt, must not be a function of the crack extension rate, da/dt. Ernst devised the modified-J parameter, J(sub M), that meets this criterion. J(sub M) correlates fracture data to much higher crack growth than does J. Ultimately, a limit of the validity of J(sub M) is anticipated, and this has been estimated to be at a crack extension of about 40-percent of the initial remaining ligament length. None of the various parameters can be expected to describe fracture in an environment of gross plasticity, in which case the process is better described by deformation parameters, e.g., stress and strain. In the current study, various schemes to identify the onset of the plasticity-dominated behavior, i.e., the end of fracture mechanics validity, are presented. Each validity limit parameter is developed in detail, and then data is presented and the various schemes for establishing a limit of the validity are compared. The selected limiting parameter is applied to a set of fracture data showing the improvement of correlation gained.

Lambert, Dennis M.↗

NASA NTRS1998

Radiation-Driven Warping: Nonisothermal Disks - 2

Recent work by Pringle and by Maloney, Begelman, & Pringle has shown that geometrically thin, optically thick, accretion disks are unstable to warping driven by radiation torque from the central source. This work was confined to isothermal (i.e., surface density Sigma varies as R(sup -3/2) disks. In this paper we generalize the study of radiation-driven warping to include general power-law surface density distributions, Sigma varies as R(sup -delta).We consider the range from Delta = 3/2 (the isothermal case) to Delta = -3/2, which corresponds to a radiation-pressure-supported disk; this spans the range of surface density distributions likely to be found in real astrophysical disks. In all cases there are an infinite number of zero-crossing solutions (i.e., solutions that cross the equator), which are the physically relevant modes if the outer boundary of the disk is required to lie in a specified plane. However, unlike the isothermal disk, which is the degenerate case, the frequency eigenvalues for Delta does not equal 3/2 are all distinct. In all cases the location of the zero moves outward from the steady state (pure precession) value with increasing growth rate; thus, there is a critical minimum size for unstable disks. Modes with zeros at smaller radii are damped. The critical radius and the steady state precession rate depend only weakly on Delta. An additional analytic solution has been found for Delta = 1. The case Delta = 1 divides the solutions into two qualitatively different regimes. For Delta greater than or equal to 1, the fastest growing modes have maximum warp amplitude, close to the disk outer edge, and the ratio of Beta(sub max) to the warp amplitude at the disk inner edge, Beta(sub o), is much greater than 1. For Delta less than 1, Beta(sub max/Beta(sub o) approximately equals 1, and the warp maximum steadily approaches the origin as Delta decreases. This implies that nonlinear effects must be important if the warp extends to the disk inner edge for Delta greater than or equal to 1, but for Delta less than 1 nonlinearity will be important only if the warp amplitude is large at the origin. Because of this qualitative difference in the shapes of the warps, the effects of shadowing of the central source by the warp will also be very different in the two regimes of Delta. This has important implications for radiation-driven warping in X-ray binaries, for which the value of Delta characterizing the disk is likely to be less than unity. In real accretion disks the outer boundary condition is likely to be different from the zero-crossing condition that we have assumed. In accretion disks around massive black holes in active galactic nuclei, the disk will probably become optically thin before the outer disk boundary is reached, whereas in X-ray binaries there will be an outer disk region (outside the circularization radius) in which the inflow velocity is zero but angular momentum is still transported. We show that in both these cases the solutions are similar to the zero-crossing eigenfunctions.

Maloney, Philip R.↗

NASA NTRS1947

The Elasto-Plastic Stability of Plates

This article explains results developed from the following research: 'The Stability of Plates and Shells beyond the Elastic Limit.' A significant improvement is found in the derivation of the relations between the stress factors and the strains resulting from the instability of plates and shells. In a strict analysis, the problem reduces to the solution of two simultaneous nonlinear partial differential equations of the fourth order in the deflection and stress function, and in the approximate analysis to a single linear equation of the Bryan type. Solutions are given for the special cases of a rectangular plate buckling into a cylindrical form, and of an arbitrarily shaped plate under uniform compression. These solutions indicate that the accuracy obtained by the approximate method is satisfactory.

Ilyushin, A. A.↗

NASA NTRS2008

Energetic Consistency and Coupling of the Mean and Covariance Dynamics

The dynamical state of the ocean and atmosphere is taken to be a large dimensional random vector in a range of large-scale computational applications, including data assimilation, ensemble prediction, sensitivity analysis, and predictability studies. In each of these applications, numerical evolution of the covariance matrix of the random state plays a central role, because this matrix is used to quantify uncertainty in the state of the dynamical system. Since atmospheric and ocean dynamics are nonlinear, there is no closed evolution equation for the covariance matrix, nor for the mean state. Therefore approximate evolution equations must be used. This article studies theoretical properties of the evolution equations for the mean state and covariance matrix that arise in the second-moment closure approximation (third- and higher-order moment discard). This approximation was introduced by EPSTEIN [1969] in an early effort to introduce a stochastic element into deterministic weather forecasting, and was studied further by FLEMING [1971a,b], EPSTEIN and PITCHER [1972], and PITCHER [1977], also in the context of atmospheric predictability. It has since fallen into disuse, with a simpler one being used in current large-scale applications. The theoretical results of this article make a case that this approximation should be reconsidered for use in large-scale applications, however, because the second moment closure equations possess a property of energetic consistency that the approximate equations now in common use do not possess. A number of properties of solutions of the second-moment closure equations that result from this energetic consistency will be established.

Cohn, Stephen E.↗