SEARCH · Engineering Papers
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.
In this paper, a nonlinear bulging factor is derived using a strain energy approach combined with dimensional analysis. The functional form of the bulging factor contains an empirical constant that is determined using R-curve data from unstiffened flat and curved panel tests. The determination of this empirical constant is based on the assumption that the R-curve is the same for both flat and curved panels.
The experimental results of the Middick Zero-gravity Dynamics Experiment (MODE) are presented. MODE investigated the dynamics of two key spacecraft components; contained fluids and jointed truss structures. The change in the lateral slosh behavior of contained fluids and the changes in the modal parameters of space structural test articles from one- to zero-gravity were studied. The experimental apparatus used to determine the dynamic characteristics is described. For the MODE fluid experiments, the forced response characteristics of silicone oil and distilled water in cylinder tanks with either a flat or spherical bottom are reported and discussed. A comparison of the measured earth and space results identifies and highlights the effects of gravity on the linear and nonlinear slosh behavior of these fluids. For the MODE structural experiments, an analytic and experimental study is presented. Deployable, erectable, and rotary modules were assembled to form three one- and two-dimensional structures in which variations in bracing wire and rotary joint preload could be introduced. Analysis is compared with ground experimental measurements made on a spring/wire suspension system with a nominal plunge frequency of one Hertz, and with measurements made on the Shuttle middeck. The degree of change in linear modal parameters as well as the change in nonlinear nature of the response is examined. Trends in modal parameters are presented as a function of force amplitude, joint preload, and ambient gravity level.
New methods are presented that utilize the optimization of goodness-of-fit statistics in order to estimate Weibull parameters from failure data. It is assumed that the underlying population is characterized by a three-parameter Weibull distribution. Goodness-of-fit tests are based on the empirical distribution function (EDF). The EDF is a step function, calculated using failure data, and represents an approximation of the cumulative distribution function for the underlying population. Statistics (such as the Kolmogorov-Smirnov statistic and the Anderson-Darling statistic) measure the discrepancy between the EDF and the cumulative distribution function (CDF). These statistics are minimized with respect to the three Weibull parameters. Due to nonlinearities encountered in the minimization process, Powell's numerical optimization procedure is applied to obtain the optimum value of the EDF. Numerical examples show the applicability of these new estimation methods. The results are compared to the estimates obtained with Cooper's nonlinear regression algorithm.
Analysis and synthesis tools are developed to improved the dynamic performance of reconfigurable nonminimum phase, nonstrictly positive real-time variant systems. A novel Spline Varying Optimal (SVO) controller is developed for the kinematic nonlinear system. There are several advantages to using the SVO controller, in which the spline function approximates the system model, observer, and controller gain. They are: The spline function approximation is simply connected, thus the SVO controller is more continuous than traditional gain scheduled controllers when implemented on a time varying plant; ft is easier for real-time implementations in storage and computational effort; where system identification is required, the spline function requires fewer experiments, namely four experiments; and initial startup estimator transients are eliminated. The SVO compensator was evaluated on a high fidelity simulation of the Shuttle Remote Manipulator System. The SVO controller demonstrated significant improvement over the present arm performance: (1) Damping level was improved by a factor of 3; and (2) Peak joint torque was reduced by a factor of 2 following Shuttle thruster firings.
Results from structural analyses of three tow-placed AS4/977-3 composite panels with both geometric and material nonlinearities are presented. Two of the panels have variable stiffness layups where the fiber orientation angle varies as a continuous function of location on the panel planform. One variable stiffness panel has overlapping tow bands of varying thickness, while the other has a theoretically uniform thickness. The third panel has a conventional uniform-thickness [plus or minus 45](sub 5s) layup with straight fibers, providing a baseline for comparing the performance of the variable stiffness panels. Parametric finite element analyses including nonlinear material shear are first compared with material characterization test results for two orthotropic layups. This nonlinear material model is incorporated into structural analysis models of the variable stiffness and baseline panels with applied end shortenings. Measured geometric imperfections and mechanical prestresses, generated by forcing the variable stiffness panels from their cured anticlastic shapes into their flatter test configurations, are also modeled. Results of these structural analyses are then compared to the measured panel structural response. Good correlation is observed between the analysis results and displacement test data throughout deep postbuckling up to global failure, suggesting that nonlinear material behavior is an important component of the actual panel structural response.
A new method for analyzing nonlinear and nonstationary data has been developed. The key part of the method is the Empirical Mode Decomposition method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF). An IMF is defined as any function having the same numbers of zero-crossing and extrema, and also having symmetric envelopes defined by the local maxima and minima respectively. The IMF also admits well-behaved Hilbert transform. This decomposition method is adaptive, and, therefore, highly efficient. Since the decomposition is based on the local characteristic time scale of the data, it is applicable to nonlinear and nonstationary processes. With the Hilbert transform, the Intrinsic Mode Functions yield instantaneous frequencies as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum, Example of application of this method to earthquake and building response will be given. The results indicate those low frequency components, totally missed by the Fourier analysis, are clearly identified by the new method. Comparisons with Wavelet and window Fourier analysis show the new method offers much better temporal and frequency resolutions.
The nonlinear evolution of a pair of initially oblique waves in a high Reynolds Number Stokes layer is studied. Attention is focused on times when disturbances of amplitude epsilon have O(epsilon(exp 1/3)R) growth rates, where R is the Reynolds number. The development of a pair of oblique waves is then controlled by nonlinear critical-layer effects. Viscous effects are included by studying the distinguished scaling epsilon = O(R(exp -1)). This leads to a complicated modification of the kernel function in the integro-differential amplitude equation. When viscosity is not too large, solutions to the amplitude equation develop a finite-time singularity, indicating that an explosive growth can be introduced by nonlinear effects; we suggest that such explosive growth can lead to the bursts observed in experiments. Increasing the importance of viscosity generally delays the occurrence of the finite-time singularity, and sufficiently large viscosity may lead to the disturbance decaying exponentially. For the special case when the streamwise and spanwise wavenumbers are equal, the solution can evolve into a periodic oscillation. A link between the unsteady critical-layer approach to high-Reynolds-number flow instability, and the wave vortex approach is identified.
Cylindrical punch indentations are performed to determine the effective modulus of a plasma-sprayed ZrO2-8Wt%Y2O3 thermal barrier coating (TBC) as a function of coating depth. Cylindrical punch indentations offer significant advantages over pointed (Vickers, Berkovich, or Knoop) indentations for materials that do not exhibit linear elastic behavior. Cyclic loading with a cylindrical punch clearly shows the TBCs to exhibit nonlinear elastic behavior with significant hysteresis that is related to the compaction and internal sliding within the plasma-spray splat microstructure. In addition, the effect of a high heat flux laser treatment was shown to produce a gradient both in the effective TBC modulus and degree of loading/unloading hysteresis with depth.
For a well-developed sea at equilibrium with a constant wind, the energy-containing range of the wavenumber spectrum for wind-generated gravity waves is approximated by a generalized power law involving the angular spread function and mu, interpreted as a fractal codimension of a small surface patch. Dependence of mu on the wave age is estimated, and the 'Phillips constant', beta, along with the low-wavenumber boundary, k0, of the inertial subrange are analyzed on the basis of the wave action and energy conservation principles. The resulting expressions are employed to evaluate various non-Gaussian statistics of a weakly nonlinear sea surface, which determine the sea state bias in satellite altimetry. The locally accelerated decay of the spectral density function in a high-wavenumber dissipation subrange is pointed out as an important factor of wave dynamics and the geometrical optics treatment of the sea state bias. The analysis is carried out in the approximation of a unidirectional wave field and confined to the case of a well-developed sea.
An investigation is conducted for the closed loop stability of linear time-invariant systems controlled by linear quadratic (LQ) regulators, in cases where nonlinearities exist in the control channels lying outside the stability sector in regions away from the origin. The estimate of the region of attraction thus obtained furnishes methods for the selection of performance function weights for more robust LQ designs. Attention is then given to the closed loop stability of linear time-invariant systems controlled by the LQ regulators when the nonlinearities in the loops escape the stability sector in a bounded region containing the origin.
The evolution of the two-point correlation function for the large-scale distribution of galaxies in an expanding universe is studied on the assumption that the perturbation densities lie in a Gaussian distribution centered on any given mass scale. The perturbations are evolved according to the Friedmann equation, and the correlation function for the resulting distribution of perturbations at the present epoch is calculated. It is found that: (1) the computed correlation function gives a satisfactory fit to the observed function in cosmological models with a density parameter (Omega) of approximately unity, provided that a certain free parameter is suitably adjusted; (2) the power-law slope in the nonlinear regime reflects the initial fluctuation spectrum, provided that the density profile of individual perturbations declines more rapidly than the -2.4 power of distance; and (3) both positive and negative contributions to the correlation function are predicted for cosmological models with Omega less than unity.
A reduced-order model of a flapping airfoil is developed using Proper Orthogonal Decomposition (POD). The proper basis functions, developed from snapshots of full Navier-Stokes simulations, are used for a Galerkin projection of the governing equations. The resulting coupled, nonlinear ordinary di.erential equations have a low dimension because the first few basis members capture most of the energy of the flow. The reduced-order model is used to simulate heaving motions that are both similar to and different from the motion(s) used to generate the basis functions, and the errors in the model are quantified. Several methods are used to generate mode sets that can be used over a range of heaving parameters, including snapshots from one, two, and multiple Navier-Stokes simulations. As snapshots from additional simulations are added to the decomposition, the mode sets become richer and can simulate a wider range of parameter space, at some computational cost. Whereas the POD method is fully applicable in three dimensions, the simulation technique based on a body-fixed and body-fitted grid suffers large overhead when extended to three dimensions. To reduce the overhead, an embedding technique is discussed which embeds the solid wing into a fixed Cartesian grid. The wing, which can now have multiple pieces and also be flexible, is represented by a distribution of body forces. This distribution is determined to give exactly the flow around a flapping wing.
An approach is presented supporting analysis, modeling, and test validation of operational flight instrumentation (OFI) that facilitates critical functions for the Space Launch System (SLS) main propulsion system (MPS). Certain types of OFI sensors were shown to exhibit highly nonlinear and non-gaussian noise characteristics during acceptance testing, motivating the development of advanced modeling and simulation (M&S) capability to support algorithm verification and flight certification. Hardware model and algorithm simulation fidelity was informed by a risk scoring metric; redesign of high-risk algorithms using test-validated sensor models significantly improved their expected performance as evaluated using Monte Carlo acceptance sampling methods. Autonomous functions include closed-loop ullage pressure regulation, pressurant leak detection, and fault isolation for automated safing and crew caution and warning (C&W).
Liapunov direct method applied to design of controllers for class of nonlinear time-varying processes
Analytical redundancy is applied to the pitch axis of an aircraft to enhance sensor fault detection. Three different designs based on monitoring Kalman filter residuals or likelihood functions were developed. Sensor failures were simulated to evaluate the performance of the algorithms on a nonlinear simulation of an F-8C aircraft. Based on requirements for sensor fault detection and preliminary simulator performance it is concluded that these algorithms have the potential for reducing the level of sensor redundancy in a digital flight control system without compromising performance.
The nonlinear Schroedinger equation with linear growth and damping is truncated to three waves. The resulting system of nonlinear ordinary differential equations describes the excitation of linearly damped waves by the oscillating two-stream instability driven by a linearly unstable pump wave. This system represents a simple model for the nonlinear saturation of a linearly unstable wave. The model is examined analytically and numerically as a function of the dimensionless parameters of the system. It is found that the model can exhibit a wealth of characteristic dynamical behavior including stationary equilibria, Hopf bifurcations to periodic orbits, period doubling bifurcations, chaotic solutions characteristic of a strange attractor, tangent bifurcations from chaotic to periodic solutions, transient chaos, and hysteresis. Many of these features are shown to be explainable on the basis of one-dimensional maps. In the case of chaotic solutions, evidence for the presence of a strange attractor is provided by demonstrating Cantor set-like structure (i.e., scale invariance) in the surface of section.
For transonic flow calculations, a combination of the successive line over-relaxation (SLOR) and the preconditioned conjugate gradient (CG) method has been suggested by Wong and Hafez (1981). This paper studies the method of minimal residual (MR) which avoids a combined iteration. This method is closely related to the CG method, may be regarded as a first-order gradient method, and is applicable to symmetric and nonsymmetric matrices. The problem is formulated as a nonlinear mixed elliptic-hyperbolic partial differential equation which includes an artificial viscosity and a switching function which is zero in subsonic regions and nonzero in supersonic regions. Alternatives to the SLOR method which provide faster convergence rates are introduced. The preconditioned MR algorithm is developed, and transonic potential flows around NACA 0012 airfoil are calculated for different Mach numbers and angles of attack. Preliminary results are presented, demonstrating that the MR algorithm requires no parameter estimation and rapidly converges for subsonic flows.
Measurements of marginally unstable electron distribution functions and high-frequency plasma waves were made on a sounding rocket flight through a quiet auroral arc. The waves appeared near the electron plasma frequency and had a large parallel electric field component such that k-parallel is greater than k-perpendicular. The appearance of these waves was correlated with the presence of marginally unstable parallel electron distributions. Analysis has shown that the waves were produced by parallel electron distribution function greater than 0 rather than the small perpendicular electron distribution function greater than 0 features. Wave levels and growth rates inside the arc were small, and nonlinear wave-wave and wave-particle interactions appear to have been minimal.