Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear function”

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 1,027 records · Page 57

Intraseasonal oscillations in the global atmosphere. II - Southern Hemisphere

A systematic examination of oscillatory modes in the global atmosphere is completed by studying 12 years of 500 mb geopotential heights in the Southern Hemisphere. The data were band-pass filtered to focus on intraseasonal (IS) phenomena, and spatial empirical orthogonal functions (EOFs) were obtained. The leading principal components were subjected to singular spectrum analysis in order to identify nonlinear IS oscillations with high statistical confidence. In the Southern Hemisphere, the dominant mode has a period of 23 days, with spatial patterns carried by the second and third winter EOF of the IS band. It has a zonal wavenumber-four structure. The 40-day mode is second, and is dominated by wavenumbers three and four, while a 16-day mode is too weak to separate its spatial behavior from the previous two. The IS dynamics in the Southern Hemisphere is more complex and dominated by shorter wavenumbers than the Northern Hemisphere. No statistically significant correlations between the Southern Hemisphere and the tropics or the Northern Hemisphere are apparent in the IS band.

Ghil, Michael↗

The influence of fiber waviness on the compressive behavior of unidirectional continuous fiber composites

The influence of fiber waviness and matrix nonlinearity on the compressive behavior of continuous fiber composites is studied. A micromechanics model based on the kinematics of the fibers has been developed to predict the behavior of unidirectional composites with initially wavy fibers under compressive loads. The initial waviness has been idealized as sinusoidal. Nonlinear shear behavior of the matrix has been included. Because the shear strain in the matrix is a function of position, the deformed shape of the fibers differs from the initial shape. Therefore, the deformed fiber shape is represented by a sine series. Fiber waviness and matrix shear properties have been determined experimentally for T650-42/Radel C under several conditions. These experimental data have been used as input into the analytical model. Predicted compressive response is in reasonable agreement with experimental data.

Highsmith, Alton L.↗

Squeezing in phase-conjugated resonance fluorescence

Emission of resonance fluorescence by an atom near the surface of a four-wave mixing phase conjugator is considered. The dipole radiation field, regarded as a Heisenberg-operator field, is decomposed into plane waves with the aid of Weyl's representation of the Green's function for the wave equation. Each plane-wave component which is incident on the surface of the nonlinear medium, is reflected as its phase-conjugate image. Summation of all reflected plane waves then yields the phase conjugate replica of the incident dipole radiation. This field adds to the radiation which is emitted by the atom into the direction away from the medium. The condition under which squeezing occurs in the emitted resonance fluorescence is investigated.

Arnoldus, Henk F.↗

Investigation on the use of optimization techniques for helicopter airframe vibrations design studies

Results of the investigation of formal nonlinear programming-based numerical optimization techniques of helicopter airframe vibration reduction are summarized. The objective and constraint function and the sensitivity expressions used in the formulation of airframe vibration optimization problems are presented and discussed. Implementation of a new computational procedure based on MSC/NASTRAN and CONMIN in a computer program system called DYNOPT for optimizing airframes subject to strength, frequency, dynamic response, and dynamic stress constraints is described. An optimization methodology is proposed which is thought to provide a new way of applying formal optimization techniques during the various phases of the airframe design process. Numerical results obtained from the application of the DYNOPT optimization code to a helicopter airframe are discussed.

Sreekanta Murthy, T.↗

A Parallel Multilevel Spectral Element Scheme

A parallel multilevel strategy is developed using spectral (p) finite elements. Hierarchic bases are particularly well suited since the element matrices and vectors are nested and the multilevel projections easiliy performed. Since the basis degree is used to specify the multigrid level, an EBE strategy is natural br the multilevel technique. Results are presented for two candidate nonlinear elliptic transport problems: the augmented drift-diffusion equations of semiconductor device modeling and the stream function-vorticity equations of incompressible fluid dynamics.

Davis, M. B.↗

Deformation Response of Unsymmetrically Laminated Plates Subjected to Inplane Loading

This paper discusses the out-of-plane deformation behavior of unsymmetric cross-ply composite plates compressed inplane by displacing one edge of the plate a known amount. The plates are assumed to be initially flat and several boundary conditions are considered. Geometrically nonlinear behavior is assumed. The primary objectives are to study the out-of-plane behavior as a function of increasing inplane compression and to determine if bifurcation behavior and secondary buckling can occur. It is shown that, depending on the boundary conditions, both can occur, though the characteristics are different than the pre and post-buckling behavior of a companion symmetric cross-ply plate. Furthermore, while a symmetric cross-ply plate can postbuckle with either a positive or negative out-of-plane displacement, the unsymmetric cross-ply plates studied deflect out-of-plane only in one direction throughout the range of inplane compression, the direction again depending on the boundary conditions

Ochinero, Tomoya T.↗

Reduced Order Methods for Prediction of Thermal-Acoustic Fatigue

The goal of this investigation is to assess the quality of high-cycle-fatigue life estimation via a reduced order method, for structures undergoing random nonlinear vibrations in a presence of thermal loading. Modal reduction is performed with several different suites of basis functions. After numerically solving the reduced order system equations of motion, the physical displacement time history is obtained by an inverse transformation and stresses are recovered. Stress ranges obtained through the rainflow counting procedure are used in a linear damage accumulation method to yield fatigue estimates. Fatigue life estimates obtained using various basis functions in the reduced order method are compared with those obtained from numerical simulation in physical degrees-of-freedom.

Przekop, A.↗

Hamilton's Equations with Euler Parameters for Rigid Body Dynamics Modeling

A combination of Euler parameter kinematics and Hamiltonian mechanics provides a rigid body dynamics model well suited for use in strongly nonlinear problems involving arbitrarily large rotations. The model is unconstrained, free of singularities, includes a general potential energy function and a minimum set of momentum variables, and takes an explicit state space form convenient for numerical implementation. The general formulation may be specialized to address particular applications, as illustrated in several three dimensional example problems.

Shivarama, Ravishankar↗

Using an Ultrasonic Instrument to Size Extravascular Bubbles

In an ongoing development project, microscopic bubbles in extravascular tissue in a human body will be detected by use of an enhanced version of the apparatus described in Ultrasonic Bubble- Sizing Instrument (MSC-22980), NASA Tech Briefs, Vol. 24, No. 10 (October 2000), page 62. To recapitulate: The physical basis of the instrument is the use of ultrasound to excite and measure the resonant behavior (oscillatory expansion and contraction) of bubbles. The resonant behavior is a function of the bubble diameter; the instrument exploits the diameter dependence of the resonance frequency and the general nonlinearity of the ultrasonic response of bubbles to detect bubbles and potentially measure their diameters. In the cited prior article, the application given most prominent mention was the measurement of gaseous emboli (essentially, gas bubbles in blood vessels) that cause decompression sickness and complications associated with cardiopulmonary surgery. According to the present proposal, the instrument capabilities would be extended to measure extravascular bubbles with diameters in the approximate range of 1 to 30 m. The proposed use of the instrument could contribute further to the understanding and prevention of decompression sickness: There is evidence that suggests that prebreathing oxygen greatly reduces the risk of decompression sickness by reducing the number of microscopic extravascular bubbles. By using the ultrasonic bubble-sizing instrument to detect and/or measure the sizes of such bubbles, it might be possible to predict the risk of decompression sickness. The instrument also has potential as a tool to guide the oxygen-prebreathing schedules of astronauts; high-altitude aviators; individuals who undertake high-altitude, low-opening (HALO) parachute jumps; and others at risk of decompression sickness. For example, an individual at serious risk of decompression sickness because of high concentrations of extravascular microscopic bubbles could be given a warning to continue to prebreathe oxygen until it was safe to decompress.

Magari, Patrick J.↗

Very Large Eddy Simulations of a Jet-A Spray Reacting Flow in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar DWFDF Method

This paper presents the very large eddy simulations (VLES) of a Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar DWFDF method, in which DWFDF is defined as the density weighted time filtered fine grained probability density function. The flow field is calculated by using the time filtered compressible Navier-Stokes equations (TFNS) with nonlinear subscale turbulence models, and when the Eulerian scalar DWFDF method is invoked, the energy and species mass fractions are calculated by solving the equation of DWFDF. A nonlinear subscale model for closing the convection term of the Eulerian scalar DWFDF equation is used and will be briefly described in this paper. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar DWFDF method in both improving the simulation quality and maintaining economic computing cost are observed.

Shih, Tsan-Hsing↗

Electric Field Conjugation in the Presence of Model Uncertainty

The Wide-Field Infrared Survey Telescope (WFIRST) is a 2.4-meter-diameter space telescope NASA program. The payload will include a coronagraph instrument (CGI). The CGI designs under development use deformable mirrors (DM) to create a point spread function (PSF) with a dark region around the obscured star object. Electric field conjugation (EFC) is an iterative nonlinear optimization procedure that uses measurements of the electric field at the image to determine the DM displacements to modify the PSF to create the region of high contrast in the image. EFC requires a numerical model of the coronagraph to calculate the Jacobian of the system, which is used, along with regularization, to solve for the DM displacements for each iteration of the nonlinear optimization. Ideally, the coronagraph is aligned and calibrated, and the calibration data are used in the numerical model for calculating the Jacobian. However, calibration and alignment measurements always contain uncertainty resulting in calibration error. Therefore, the Jacobian calculated from the numerical model is not an exact representation of the physical coronagraph, and the resulting DM solution for an EFC iteration does not have the exact impact on the electric field of the coronagraph as predicted by the EFC. The result is slow convergence, and, as will be shown, drives the use of more restrictive regularization. Using Monte Carlo trials, we investigate the effect of calibration error on EFC convergence and regulatization. Comparison to results from the High Contrast Imaging Testbed Hybrid Lyot Coronagraph are also presented.

Marx, David↗

Modeling the bidirectional reflectance distribution function of mixed finite plant canopies and soil

An analytical model of the bidirectional reflectance for optically semi-infinite plant canopies has been extended to describe the reflectance of finite depth canopies contributions from the underlying soil. The model depends on 10 independent parameters describing vegetation and soil optical and structural properties. The model is inverted with a nonlinear minimization routine using directional reflectance data for lawn (leaf area index (LAI) is equal to 9.9), soybeans (LAI, 2.9) and simulated reflectance data (LAI, 1.0) from a numerical bidirectional reflectance distribution function (BRDF) model (Myneni et al., 1988). While the ten-parameter model results in relatively low rms differences for the BRDF, most of the retrieved parameters exhibit poor stability. The most stable parameter was the single-scattering albedo of the vegetation. Canopy albedo could be derived with an accuracy of less than 5% relative error in the visible and less than 1% in the near-infrared. Sensitivity were performed to determine which of the 10 parameters were most important and to assess the effects of Gaussian noise on the parameter retrievals. Out of the 10 parameters, three were identified which described most of the BRDF variability. At low LAI values the most influential parameters were the single-scattering albedos (both soil and vegetation) and LAI, while at higher LAI values (greater than 2.5) these shifted to the two scattering phase function parameters for vegetation and the single-scattering albedo of the vegetation. The three-parameter model, formed by fixing the seven least significant parameters, gave higher rms values but was less sensitive to noise in the BRDF than the full ten-parameter model. A full hemispherical reflectance data set for lawn was then interpolated to yield BRDF values corresponding to advanced very high resolution radiometer (AVHRR) scan geometries collected over a period of nine days. The resulting parameters and BRDFs are similar to those for the full sampling geometry, suggesting that the limited geometry of AVHRR measurements might be used to reliably retrieve BRDF and canopy albedo with this model.

Schluessel, G.↗

Reduced-Order Modeling for Flutter/LCO Using Recurrent Artificial Neural Network

The present study demonstrates the efficacy of a recurrent artificial neural network to provide a high fidelity time-dependent nonlinear reduced-order model (ROM) for flutter/limit-cycle oscillation (LCO) modeling. An artificial neural network is a relatively straightforward nonlinear method for modeling an input-output relationship from a set of known data, for which we use the radial basis function (RBF) with its parameters determined through a training process. The resulting RBF neural network, however, is only static and is not yet adequate for an application to problems of dynamic nature. The recurrent neural network method [1] is applied to construct a reduced order model resulting from a series of high-fidelity time-dependent data of aero-elastic simulations. Once the RBF neural network ROM is constructed properly, an accurate approximate solution can be obtained at a fraction of the cost of a full-order computation. The method derived during the study has been validated for predicting nonlinear aerodynamic forces in transonic flow and is capable of accurate flutter/LCO simulations. The obtained results indicate that the present recurrent RBF neural network is accurate and efficient for nonlinear aero-elastic system analysis

Yao, Weigang↗

Boundary-element shape sensitivity analysis for thermal problems with nonlinear boundary conditions

Implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solid objects subjected to nonlinear boundary conditions is shown to generate an accurate and economical approach for the computation of shape sensitivities for this class of problems. This approach involves the employment of analytical derivatives of boundary-element kernel functions with respect to shape design variables. A formulation is presented that can consistently account for both temperature-dependent convection and radiation boundary conditions. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems are presented that have closed-form solutions.

Kane, James H.↗

The explicit computation of integration algorithms and first integrals for ordinary differential equations with polynomials coefficients using trees

This note is concerned with the explicit symbolic computation of expressions involving differential operators and their actions on functions. The derivation of specialized numerical algorithms, the explicit symbolic computation of integrals of motion, and the explicit computation of normal forms for nonlinear systems all require such computations. More precisely, if R = k(x(sub 1),...,x(sub N)), where k = R or C, F denotes a differential operator with coefficients from R, and g member of R, we describe data structures and algorithms for efficiently computing g. The basic idea is to impose a multiplicative structure on the vector space with basis the set of finite rooted trees and whose nodes are labeled with the coefficients of the differential operators. Cancellations of two trees with r + 1 nodes translates into cancellation of O(N(exp r)) expressions involving the coefficient functions and their derivatives.

Crouch, P. E.↗

Fitting aerodynamic forces in the Laplace domain: An application of a nonlinear nongradient technique to multilevel constrained optimization

A technique which employs both linear and nonlinear methods in a multilevel optimization structure to best approximate generalized unsteady aerodynamic forces for arbitrary motion is described. Optimum selection of free parameters is made in a rational function approximation of the aerodynamic forces in the Laplace domain such that a best fit is obtained, in a least squares sense, to tabular data for purely oscillatory motion. The multilevel structure and the corresponding formulation of the objective models are presented which separate the reduction of the fit error into linear and nonlinear problems, thus enabling the use of linear methods where practical. Certain equality and inequality constraints that may be imposed are identified; a brief description of the nongradient, nonlinear optimizer which is used is given; and results which illustrate application of the method are presented.

Tiffany, S. H.↗

A Linear Programming Approach to Routing Control in Networks of Constrained Nonlinear Positive Systems with Concave Flow Rates

We consider control design for positive compartmental systems in which each compartment's outflow rate is described by a concave function of the amount of material in the compartment.We address the problem of determining the routing of material between compartments to satisfy time-varying state constraints while ensuring that material reaches its intended destination over a finite time horizon. We give sufficient conditions for the existence of a time-varying state-dependent routing strategy which ensures that the closed-loop system satisfies basic network properties of positivity, conservation and interconnection while ensuring that capacity constraints are satisfied, when possible, or adjusted if a solution cannot be found. These conditions are formulated as a linear programming problem. Instances of this linear programming problem can be solved iteratively to generate a solution to the finite horizon routing problem. Results are given for the application of this control design method to an example problem. Key words: linear programming; control of networks; positive systems; controller constraints and structure.

Positive Systems↗

Approximate optimal atmospheric entry trajectories.

For a reentry glider, approximate solutions are found in closed form for the problem of maximizing a function of the terminal velocity, altitude, flight-path angle and heading angle subject to, at most, three terminal nonlinear constraints. The results given here extend the previous results in two important ways: (1) the second-order approximation of the entry dynamics of Loh is used instead of the first-order approximation of Allen and Eggers. This approximation is found to compare extremely well with an exact numerical optimal path. (2) A three-dimensional optimization problem is solved which includes Loh's second-order approximation where the roll angle as well as the lift coefficient are determined subject to constraints which include the terminal heading angle. Furthermore, closed-form solutions can also be obtained subject to in-flight constraints on the lift coefficient and roll angle.

Speyer, J. L.↗