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.
A nonlinear procedure for computing the pressure distribution on an airfoil with a trailing edge jet flap is described. The method is not restricted to thin airfoils or shallow jet-deflection angles. Correlation with experiment indicates that the characteristics of the pressure distribution are predicted by the theory, but the effect of entrainment is overpredicted with the entrainment coefficient used.
A method for rapidly estimating the overall forces and moments at supercritical speeds, below drag divergence, of transport configurations with supercritical wings is presented. The method was also used for estimating the rolling moments due to the deflection of wing trailing-edge controls. This analysis was based on a vortex-lattice technique modified to approximate the effects of wing thickness and boundary-layer induced camber. Comparisons between the results of this method and experiment indicate reasonably good correlation of the lift, pitching moment, and rolling moment. The method required much less storage and run time to compute solutions over an angle-of-attack range than presently available transonic nonlinear methods require for a single angle-of-attack solution.
The effect of static tensile load on the thermal expansion of Gr/PI composite material was measured for seven different laminate configurations. A computer program was developed which implements laminate theory in a piecewise linear fashion to predict the coupled nonlinear thermomechanical behavior. Static tensile load significantly affected the thermal expansion characteristics of the laminates tested. This effect is attributed to a fiber instability micromechanical behavior of the constituent materials. Analytical results correlated reasonably well with free thermal expansion tests (no load applied to the specimen). However, correlation was poor for tests with an applied load.
A simple method for measuring the pulse lengths of synchronously pumped dye lasers which does not require the use of an external nonlinear medium, such as a doubling crystal or two-photon fluorescence cell, to autocorrelate the pulses is discussed. The technique involves feeding the laser pulses back into the dye jet, thus correlating the output pulses with the intracavity pulses to obtain pulse length signatures in the resulting time-averaged laser power. Experimental measurements were performed using a rhodamine 6G dye laser pumped by a mode-locked frequency-doubled Nd:YAG laser. The results agree well with numerical computations, and the method proves effective in determining lengths of picosecond laser pulses.
A full-scale crash test of a prototype composite helicopter was performed at the Impact Dynamics Research Facility at NASA Langley Research Center in 1999 to obtain data for validation of a finite element crash simulation. The helicopter was the flight test article built by Sikorsky Aircraft during the Advanced Composite Airframe Program (ACAP). The composite helicopter was designed to meet the stringent Military Standard (MIL-STD-1290A) crashworthiness criteria and was outfitted with two crew and two troop seats and four anthropomorphic dummies. The test was performed at 38-ft/s vertical and 32.5-ft/s horizontal velocity onto a rigid surface. An existing modal-vibration model of the Sikorsky ACAP helicopter was converted into a model suitable for crash simulation. A two-stage modeling approach was implemented and an external user-defined subroutine was developed to represent the complex landing gear response. The crash simulation was executed with a nonlinear, explicit transient dynamic finite element code. Predictions of structural deformation and failure, the sequence of events, and the dynamic response of the airframe structure were generated and the numerical results were correlated with the experimental data to validate the simulation. The test results, the model development, and the test-analysis correlation are described.
An expression for the J-integral of a nonlinear elastic material is derived for an advancing crack in a tapered double cantilever beam fracture specimen. The elastic and plastic fracture energies related to the test geometry and how these energies correlates to the crack position are discussed. The dimensionless shape factors eta(sub el and eta(sub p) are shown to be equivalent and the deformation J-integral is analyzed in terms of the eta(sub el) function. The fracture results from a structural epoxy are interpreted using the discussed approach. The magnitude of the plastic dissipation is found to strongly depend upon the initial crack shape.
A simple continuum damage mechanics (CDM) based 3D progressive damage analysis (PDA) tool for laminated composites was developed and implemented as a user defined material subroutine to link with a commercially available explicit finite element code. This PDA tool uses linear lamina properties from standard tests, predicts damage initiation with an easy-to-implement Hashin-Rotem failure criteria, and in the damage evolution phase, evaluates the degradation of material properties based on the crack band theory and traction-separation cohesive laws. It follows Matzenmiller et al.'s formulation to incorporate the degrading material properties into the damaged stiffness matrix. Since nonlinear shear and matrix stress-strain relations are not implemented, correction factors are used for slowing the reduction of the damaged shear stiffness terms to reflect the effect of these nonlinearities on the laminate strength predictions. This CDM based PDA tool is implemented as a user defined material (VUMAT) to link with the Abaqus/Explicit code. Strength predictions obtained, using this VUMAT, are correlated with test data for a set of notched specimens under tension and compression loads.
Propellant tank slosh dynamics are typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control (GN&C) analysis. For a partially-filled smooth wall propellant tank, the critical damping based on classical empirical correlation is as low as 0.05%. Due to this low value of damping, propellant slosh is potential sources of disturbance critical to the stability of launch and space vehicles. It is postulated that the commonly quoted slosh damping is valid only under the linear regime where the slosh amplitude is small. With the increase of slosh amplitude, the critical damping value should also increase. If this nonlinearity can be verified and validated, the slosh stability margin can be significantly improved, and the level of conservatism maintained in the GN&C analysis can be lessened. The purpose of this study is to explore and to quantify the dependence of slosh damping with slosh amplitude. Accurately predicting the extremely low damping value of a smooth wall tank is very challenging for any Computational Fluid Dynamics (CFD) tool. One must resolve thin boundary layers near the wall and limit numerical damping to minimum. This computational study demonstrates that with proper grid resolution, CFD can indeed accurately predict the low damping physics from smooth walls under the linear regime. Comparisons of extracted damping values with experimental data for different tank sizes show very good agreements. Numerical simulations confirm that slosh damping is indeed a function of slosh amplitude. When slosh amplitude is low, the damping ratio is essentially constant, which is consistent with the empirical correlation. Once the amplitude reaches a critical value, the damping ratio becomes a linearly increasing function of the slosh amplitude. A follow-on experiment validated the developed nonlinear damping relationship. This discovery can lead to significant savings by reducing the number and size of slosh baffles in liquid propellant tanks.
Two stand-alone analyzers constructed for real time simulation of the complex dynamic characteristics of horizontal-axis wind energy systems are described. Mathematical models for an aeroelastic rotor, including nonlinear aerodynamic and elastic loads, are implemented with high speed digital and analog circuitry. Models for elastic supports, a power train, a control system, and a rotor gimbal system are also included. Limited correlation efforts show good comparisons between results produced by the analyzers and results produced by a large digital simulation. The digital simulation results correlate well with test data.
Consistent with the hypothesis that plankton delta C-14 and (CO2(aq)) are inversely related, increases in both sinking and suspended particulate organic matter (POM) delta C-13 detected by the Joint Global Ocean Flux Study (JGOFS) were highly negatively correlated with mixed-layer (CO2(aq)). A model of plant delta C-13 by Farquhar et al. (1982) is adapted to show that under a constant phytoplankton demand for CO2 an inverse nonlinear suspended POM delta C-13 response to ambient (CO2(aq)) is expected. Differences between predicted and observed suspended POM delta C-13 vs. (CO2(aq)) trends and among observed relationships can be reconciled if biological CO2 demand is allowed to vary.
The effects of small-scale heterogeneity in land surface characteristics on the large-scale fluxes of water and energy in land-atmosphere system has become a central focus of many of the climatology research experiments. The acquisition of high resolution land surface data through remote sensing and intensive land-climatology field experiments (like HAPEX and FIFE) has provided data to investigate the interactions between microscale land-atmosphere interactions and macroscale models. One essential research question is how to account for the small scale heterogeneities and whether 'effective' parameters can be used in the macroscale models. To address this question of scaling, three modeling experiments were performed and are reviewed in the paper. The first is concerned with the aggregation of parameters and inputs for a terrestrial water and energy balance model. The second experiment analyzed the scaling behavior of hydrologic responses during rain events and between rain events. The third experiment compared the hydrologic responses from distributed models with a lumped model that uses spatially constant inputs and parameters. The results show that the patterns of small scale variations can be represented statistically if the scale is larger than a representative elementary area scale, which appears to be about 2 - 3 times the correlation length of the process. For natural catchments this appears to be about 1 - 2 sq km. The results concerning distributed versus lumped representations are more complicated. For conditions when the processes are nonlinear, then lumping results in biases; otherwise a one-dimensional model based on 'equivalent' parameters provides quite good results. Further research is needed to fully understand these conditions.
Lorentz-broadening and pressure-induced line-shift coefficients are determined experimentally for 126 transitions in the nu2 band of D2(O-16) from laboratory spectra with a Fourier-transform spectrometer. A nonlinear least-squares fitting procedure is employed to measure transitions up to J double-prime = 12 in the buffer gases air, nitrogen, and oxygen. The N2-broadening/O2-broadening coefficient ratio is correlated with the magnitude of the O2-broadening coefficient, and the pressure coefficients are not correlated with line widths.
This paper presents a brief overview of the observed evolution in a variety of quantities describing the turbulent evolution of the interplanetary plasma and describes simulation results consistent with many features of the evolution. The turbulence is manifested through a dissipation at small scales in the inner heliosphere with a corresponding evolution in the breakpoint between a relatively flat and a Kolmogoroff spectrum; an evolution from kinetically to (slightly) magnetically dominated energy of the plasma fluctuations; a general decrease in the cross helicity or 'Alfvenicity'; changes in the anisotropy of the fluctuations; and the increasing predominance of quasi-pressure-balanced structures in the compressive component of the fluctuations. MHD simulations with shear layers either side of a central current sheet show that even in the absence of compressibility the lack of a mean field along the direction of the main flow in the current sheet leads to rapid nonlinear evolution and the observed characteristics of 'Elsasser spectra' of the fields in the inner heliosphere. Adding compressibility to the simulations does not greatly change the 'incompressive' quantities but leads in addition to observed correlations between a measure of compression and other quantities.
The objective of our study is to model the nonlinear behavior of a near-critical liquid following a rapid change of the temperature and/or other thermodynamic parameters (pressure, external electric or gravitational field). The thermodynamic critical point is manifested by large, strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, concentration in binary solutions) in the critical range of scales. The largest critical length scale is the correlation radius r(sub c). According to the scaling theory, r(sub c) increases as r(sub c) = r(sub 0)epsilon(exp -alpha) when the nondimensional distance epsilon = (T - T(sub c))/T(sub c) to the critical point decreases. The normal gravity alters the nature of correlated long-range fluctuations when one reaches epsilon approximately equal to 10(exp -5), and correspondingly the relaxation time, tau(r(sub c)), is approximately equal to 10(exp -3) seconds; this time is short when compared to the typical experimental time. Close to the critical point, a rapid, relatively small temperature change may perturb the thermodynamic equilibrium on many scales. The critical fluctuations have a hierarchical structure, and the relaxation involves many length and time scales. Above the critical point, in the one-phase region, we consider the relaxation of the liquid following a sudden temperature change that simultaneously violates the equilibrium on many scales. Below T(sub c), a non-equilibrium state may include a distribution of small scale phase droplets; we consider the relaxation of such a droplet following a temperature change that has made the phase of the matrix stable.
Prediction of laminar-turbulent transition in boundary layer flows is an important component of predicting the aerodynamic performance of a number of aerospace configurations. According to the CFD Vision 2030 [1], transition modeling represents acriticalarea in CFD simulation capability that will remain a pacing item for the foreseeable future. The fact thattransition can take placevia either one of a myriad possible paths adds to the challenges inreliable transition predictions, despite a limited knowledge of the relevant input parameters. In the low disturbance environments typical of flight applications, transition is often initiated by small amplitude disturbances in the form of linear instability waves of the laminar boundary layer. These disturbances amplify linearly at first and eventually undergo a sequence of nonlinear interactions that result in transition to turbulence. Because the nonlinear phase is rather rapid, the amplification of boundary layer instabilities is governed by the linearstability theory over a majority of the distance leading up to the onset of transition. Semi-empirical transition correlations based on the linear stability theory have been successful in explaining the observed trends in transition location within a broad class of flows. However, the application of stability theory is highly non-robust and often requires a significant domain expertise. Recent work at the NASA Langley Research Center has beenaimed at bridging the gap between physics based transition analyses such as those based on linear stability theory and practical applications that require transition prediction by users that may not be well versed in transition physics. The applications of deep learning have been at the center of these efforts. This presentation will focus on the progress achieved thus far, highlighting the applications of neural networks to selectedtransition scenarios across a range of Mach numbers and flow configuration, as well as the lessons learnedand remaining challengeswithrespect to the selection of training data and neural networks architectures, hyperparameter tuning, and the physical insights distilled from the otherwise black-box models.
A new methodology for discrimination is proposed. This is based on kernel orthonormalized partial least squares (PLS) dimensionality reduction of the original data space followed by support vector machines for classification. Close connection of orthonormalized PLS and Fisher's approach to linear discrimination or equivalently with canonical correlation analysis is described. This gives preference to use orthonormalized PLS over principal component analysis. Good behavior of the proposed method is demonstrated on 13 different benchmark data sets and on the real world problem of the classification finger movement periods versus non-movement periods based on electroencephalogram.
Discrete attitude commands have become a standard task for flying qualities evaluation and control system testing. Much pilot opinion data is now available for ground-based and in-flight simulations, but adequate performance measures and prediction methods have not been established. The Step Target Tracking Prediction method, introduced in 1978, correlated time-on-target and rms tracking data with NT-33 in-flight longitudinal simulations, but did not employ parameters easily measured in manned flight and simulation. Recent application of the Step Target Tracking Prediction method to lateral flying qualities analysis has led to a new measure of performance. This quantity, called Maximum Normalized Rate (MNR), reflects the greatest attitude rate a pilot can employ during a discrete maneuver without excessive overshoot and oscillation. MNR correlates NT-33 lateral pilot opinion ratings well, and is easily measured during flight test or simulation. Futhermore, the Step Target MNR method can be used to analyze large amplitude problems concerning rate limiting and nonlinear aerodynamics.
Discrete attitude commands have become a standard task for flying qualities evaluation and control system testing. Much pilot opinion data is now available for ground-based and in-flight simulations, but adequate performance measures and prediction methods have not been established. The Step Target Tracking Prediction method, introduced in 1978, correlated time-on-target and rms tracking data with NT-33 in-flight longitudinal simulations, but did not employ parameters easily measured in manned flight and simulation. Recent application of the Step Target Tracking Prediction method to lateral flying qualities analysis has led to a new measure of performance. This quantity, called Maximum Normalized Rate (MNR), reflects the greatest attitude rate a pilot can employ during a discrete maneuver without excessive overshoot and oscillation. MNR correlates NT-33 lateral pilot opinion ratings well, and is easily measured during night test or simulation. Furthermore, the Step Target MNR method can be used to analyze large amplitude problems concerning rate limiting and nonlinear aerodynamics.