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.
An experimental design matrix was constructed to determine the effects of time and temperature water soak on the strength of Zerodur glass-ceramic. It was found that strength does increase in a nonlinear manner which is consistent with existing theories of crack tip blunting and residual stress reduction.
Dynamical mechanisms underlying the equilibration of absolute instability are examined in a nonlinear, quasigeostrophic, two-layer model. The key to understanding the nonlinear equilibration is in recognizing that linear absolute instabilities can be stabilized both by a reduction of the vertical shear and by enhancement of the mean barotropic velocity. In a localized domain, the equilibration process proceeds with the creation of locally convectively unstable regions downstream, which encroach onto the locally absolutely unstable region until the local instability is suppressed. That local instabilities exist only if absolutely unstable regions span a minimum size is verified by eigenvalue calculations of three-dimensional flows. Numerical examples suggest that this critical size is at least 9000 km for a wide range of parameter values chosen to investigate the midlatitude storm tracks. Fluctuations arising from local absolute instability obtain maximum amplitude in the downstream convectively unstable regions rather than in the absolutely unstable regions themselves. Together, these results suggest that if an equilibrated absolute instability were to occur in midlatitudes, a zonal band of surface easterlies exceeding 900 km would be required and the associated enhanced variances would not be found coincident with the regions of absolute instability.
We present results of a three-dimensional magnetohydrodynamic (MHD) simulation of magnetotail evolution initiated by a sudden occurrence or increase of spatially localized resistivity as the major expected concequence of some localized microinstability. Because of the absence of a quantitative model, possible variations of resistivity levels with current density, or the reduction thereof, are not incorporated in the present investigation. The emphasis of the study is on an investigation of the changes to the overall evolution brought about by this localization, in particular, on the disruption and diversion of the cross-tail current and the nonlinear evolution of the magnetotail instability. The immediate consequences of the occurrence of the localized resistance and the resulting electric field are a reduction and diversion of the electric current around the region of high resistivity, associated with an increase of B(sub z) ('dipolarization') at the earthward edge and a decrease of B(sub z) at the tailward edge of this region. These effects, however, are localized and do not involve a reduction of the total cross-tail current and hence do not lead to the global development of a 'substorm current wedge,' which includes not only the reduction of the cross-tail current but also the buildup of a global field-aligned current system of 'regional 1' type (toward the Earth on the dawnside and away on the duskside of the tail). Such signatures develop at a later time, as consequences of a three-dimensional tearing instability, which is triggered by the occurrence of the resistivity. These features are found in combination with plasmoid formation and ejection, quite similar to results of earlier simulations with uniform resistivity. Differences are found in the timescale of the evolution, which tends to be shorter for localized resistivity, and in the propagation of the dipolarization effects in the equatorial plane. Whereas for uniform resistivity the temporal increase in B(sub z) tends to propagate tailward, apparently due to a pileup effect in the near-Earth region, an earthward propagation is found for the localized resistivity. This propagation results from an earthward convection of the gradient in B(sub z), which is set up by the reconnection process further tailward.
It is pointed out that at the present time there is no experimental means to detect the current reduction of stratospheric ozone by chlorine from fluorocarbons relative to its unperturbed amount. The reason for this situation is related to the small changes involved and several possible masking effects. As a substitute for experimental detection, there have been many stratospheric modeling efforts. The history of these efforts has been reviewed by Cicerone et al. (1983). Results obtained by Cicerone et al. imply that the successful detection of net ozone destruction by chlorine would be delayed many years from the present time. The present paper is concerned with a reinvestigation of the photochemical problem discussed by Cicerone et al. on the basis of an independent one-dimensional model reported by Herman (1979). It is concluded that by the year 2020 the total ozone decrease due to ClX should be about 3 percent (2.5 percent with increasing methane). Such a change may already pose environmental problems.
Beryllium fracture toughness test specimens were fatigue cracked using reversed cycling with a compression load two to three times the tension load. In worked beryllium, textures may be produced which result in fatigue cracks that are out of plane with the starter notch. Specimens of hot pressed stock exhibited load displacement records which were nonlinear throughout their course. Fracture specimens of both hot pressed and forged stock showed essentially no reduction of thickness and the fracture surfaces were flat and normal to the load axis. However, the stress intensity factor at maximum load increased with decreasing thickness. Load-displacement and electric potential records for the hot pressed beryllium specimens exhibited several anomalies such as negative residual crack mouth displacements and a decrease in electrical potential with increasing load.
Linear frequency domain methods are inadequate in analyzing the 1975 Viking Orbiter (VO75) digital tape recorder servo due to dominant nonlinear effects such as servo signal limiting, unidirectional servo control, and static/dynamic Coulomb friction. The frequency loop (speed control) servo of the VO75 tape recorder is used to illustrate the analytical tools and methodology of system redundancy elimination and high order transfer function verification. The paper compares time-domain performance parameters derived from a series of nonlinear time responses with the available experimental data in order to select the best possible analytical transfer function representation of the tape transport (mechanical segment of the tape recorder) from several possible candidates. The study also shows how an analytical time-response simulation taking into account most system nonlinearities can pinpoint system redundancy and overdesign stemming from a strictly empirical design approach. System order reduction is achieved through truncation of individual transfer functions and elimination of redundant blocks.
Proposed new metrics to assess fighter aircraft agility are collected and analyzed. A framework for classification of these new agility metrics is developed and applied. A completed set of transient agility metrics is evaluated with a high fidelity, nonlinear F-18 simulation provided by the NASA Dryden Flight Research Center. Test techniques and data reduction methods are proposed. A method of providing cuing information to the pilot during flight test is discussed. The sensitivity of longitudinal and lateral agility metrics to deviations from the pilot cues is studied in detail. The metrics are shown to be largely insensitive to reasonable deviations from the nominal test pilot commands. Instrumentation required to quantify agility via flight test is also considered. With one exception, each of the proposed new metrics may be measured with instrumentation currently available. Simulation documentation and user instructions are provided in an appendix.
There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.
A reduction method which combines classical Rayleigh-Ritz modal superposition techniques with contemporary finite-element methods is applied to transient nonlinear thermal analysis of aerospace structures. The essence of the method is the use of a few thermal modes from eigenvalue analyses as basis vectors to represent the temperature response in the structure. The method is used to obtain approximate temperature histories for a portion of the Shuttle orbiter wing subject to reentry heating and for a large space antenna reflector subject to heating associated with a low Earth orbit. The reduction method has excellent potential for significant size reduction for radiation-dominated problems such as the antenna reflector. However, for conduction-dominated problems such as the Shuttle wing, especially those with complex spatial and temporal variations in the applied heating, additional work appears necessary to find alternate sources of basis vectors which will permit significant problem size reductions.
The propagation of high intensity sound waves through an air-filled porus material was studied. The material is assumed: (1) to be rigid, incompressible, and homogeneous, and (2) to be adequately described by two properties: resistivity r and porosity. The resulting wave equation is still nonlinear, however, because of the u sgn(u) term in the resistivity. The equation is solved in the frequency domain as an infinite set of coupled inhomogeneous Helmholtz equations, one for each harmonic. An approximate but analytical solution leads to predictions of excess attenuation, saturation, and phase speed reduction for the fundamental component. A more general numerical solution is used to calculate the propagation curves for the higher harmonics. The u sgn(u) nonlinearity produces a cubic distortion pattern; when the input signal is a pure tone, only odd harmonic distortion products are generated.
Recent NASA mission plans require spacecraft to undergo potentially significant maneuvers (or dynamic loading events) with large solar arrays deployed. Therefore there is an increased need to understand and possibly control the nonlinear dynamics in the spacecraft system during such maneuvers. The development of a nonlinear controller is described. The utility of using a nonlinear controller to reduce forces and motion in a solar array wing during a loading event is demonstrated. The result is dramatic reductions in system forces and motion during a 10 second loading event. A motion curve derived from the simulation with the closed loop controller is used to obtain similar benefits with a simpler motion control approach.
A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.
Electrically powered arc jets which produce thrust at high specific impulse could provide a substantial cost reduction for orbital transfer and station keeping missions. There is currently a limited understanding of the complex, nonlinear interactions in the plasma propellant which has hindered the development of high efficiency arc jet thrusters by making it difficult to predict the effect of design changes and to interpret experimental results. A computational model developed at the University of Tennessee Space Institute (UTSI) to study laser powered thrusters and radio frequency gas heaters has been adapted to provide a tool to help understand the physical processes in arc jet thrusters. The approach is to include in the model those physical and chemical processes which appear to be important, and then to evaluate our judgement by the comparison of numerical simulations with experimental data. The results of this study have been presented at four technical conferences. The details of the work accomplished in this project are covered in the individual papers included in the appendix of this report. We present a brief description of the model covering its most important features followed by a summary of the effort.
A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximations of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady-state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudo measurements is examined, including velocity, density, and volume transport at localized arrays, and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform a previous study that used an empirical nudging scheme. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.
The Environmentally Responsible Aviation Project aims to develop aircraft technologies enabling significant fuel burn and community noise reductions. Small incremental changes to the conventional metallic alloy-based 'tube and wing' configuration are not sufficient to achieve the desired metrics. One of the airframe concepts that might dramatically improve aircraft performance is a composite-based hybrid wing body configuration. Such a concept, however, presents inherent challenges stemming from, among other factors, the necessity to transfer wing loads through the entire center fuselage section which accommodates a pressurized cabin confined by flat or nearly flat panels. This paper discusses a nonlinear finite element analysis of a large-scale test article being developed to demonstrate that the Pultruded Rod Stitched Efficient Unitized Structure concept can meet these challenging demands of the next generation airframes. There are specific reasons why geometrically nonlinear analysis may be warranted for the hybrid wing body flat panel structure. In general, for sufficiently high internal pressure and/or mechanical loading, energy related to the in-plane strain may become significant relative to the bending strain energy, particularly in thin-walled areas such as the minimum gage skin extensively used in the structure under analysis. To account for this effect, a geometrically nonlinear strain-displacement relationship is needed to properly couple large out-of-plane and in-plane deformations. Depending on the loading, this nonlinear coupling mechanism manifests itself in a distinct manner in compression- and tension-dominated sections of the structure. Under significant compression, nonlinear analysis is needed to accurately predict loss of stability and postbuckled deformation. Under significant tension, the nonlinear effects account for suppression of the out-of-plane deformation due to in-plane stretching. By comparing the present results with the previously published preliminary linear analysis, it is demonstrated in the present paper that neglecting nonlinear effects for the structure and loads of interest can lead to appreciable loss in analysis fidelity.
Optimal control of nonlinear dynamic system in presence of noise
A Monte Carlo simulation model for radiation induced plasmas with nonlinear properties due to recombination was, employing a piecewise linearized predict-correct iterative technique. Several important variance reduction techniques were developed and incorporated into the model, including an antithetic variates technique. This approach is especially efficient for plasma systems with inhomogeneous media, multidimensions, and irregular boundaries. The Monte Carlo code developed has been applied to the determination of the electron energy distribution function and related parameters for a noble gas plasma created by alpha-particle irradiation. The characteristics of the radiation induced plasma involved are given.
The propagation of orographic gravity waves into an atmosphere with exponentially decreasing density is simulated with a two-dimensional, nonlinear, time-dependent numerical model. After the stationary wave is established over the mountain, the model predicts that wave breaking causes a large reduction of the vertical momentum flux in the flow, not only at levels where wave breaking is present, but also far below the lowest occurrence of overturning. More than half of the decrease in momentum flux is explained by the presence of large amplitude, downward propagating waves, which are generated in regions of wave breaking. The downward propagating waves appear almost simultaneously with overturning, and have nonzero phase speeds, suggesting a strongly nonlinear generation mechanism that depends on local wave properties. The generation of these downward propagating waves is a robust process, insensitive to mountain height, mountain width, or density scale height. These results have important implications for observational studies of orographically generated waves as well as for schemes that seek to parameterize the effects of orogoraphy in large-scale models.