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.
Random process theory applied to discrete stationary and nonstationary data processing techniques - autocorrelation, and optimum smoothing for stationary processes
Thin membrane materials were subjected to biaxial and electrostatic tensioning loads to study techniques for maintaining surface smoothness of a thin membrane antenna. The basic mechanical and electrical setup for the tests is described and preliminary measurements of surface smoothness and surface deviation are presented.
Several case studies are examined in an attempt to relate Visible infrared spin scan radiometers Atmospheric Sounder (VAS) water vapor channel radiance data to atmospheric moisture. Since the 6.7 micrometer channel detects radiation emitted by atmospheric water vapor from a rather thick layer, a scheme is developed which weighted rawinsonde observations of relative humidity in a similar manner. A comparison of the two fields (weighted relative humidity values from gridded rawinsonde observations and the VAS water vapor image) indicates that even when mesoscale data are used, the technique only presented a smooth representation of the radiance field. This weighting function is dependent on both the temperature and moisture distribution of the atmosphere. Cross sections of relative humidity from the mesoscale rawinsondes are produced which sliced the water vapor imagery in particularly interesting regions (areas of high and low blackbody temperatures). Comparisons are then made between the cross sections and the variations in the VAS radiance data along the cross section line. Bright regions in the water vapor imagery (cold blackbody temperatures) are often associated with thin layers of high relative humidity or thick layers of only moderately moist air. The reverse also applies, i.e., dark regions (warm blackbody temperatures) in the imagery corresponded well with very dry regions in the cross sections.
Several case studies are examined in an attempt to relate Visible infrared spin scan radiometers Atmospheric Sounder (VAS) water vapor channel radiance data to atmospheric moisture. Since the 6.7 micrometer channel detects radiation emitted by atmospheric water vapor from a rather thick layer, a scheme is developed which weighted rawinsonde observations of relative humidity in a similar manner. A comparison of the two fields (weighted relative humidity values from gridded rawinsonde observations and the VAS water vapor image) indicates that even when mesoscale data are used, the technique only presented a smooth representation of the radiance field. This weighting function is dependent on both the temperature and moisture distribution of the atmosphere. Cross sections of relative humidity from the mesoscale rawinsondes are produced which sliced the water vapor imagery in particularly interesting regions (areas of high and low blackbody temperatures). Comparisons are then made between the cross sections and the variations in the VAS radiance data along the cross section line. Bright regions in the water vapor imagery (cold blackbody temperatures) are often associated with thin layers of high relative humidity or thick layers of only moderately moist air. The reverse also applies, i.e., dark regions (warm blackbody temperatures) in the imagery corresponded well with very dry regions in the cross sections.
Many modern metal and composite airframe manufacturing techniques can provide surface smoothness which is compatible with natural laminar flow (NLF) requirements. An important consideration is manufacturing roughness of the surface in the form of steps and gaps perpendicular to the freestream. The principal challenge to the design and manufacture of laminar flow surfaces today appears to be in the installation of leading-edge panels on wing, nacelle, and empennage surfaces. A similar challenge is in the installation of access panels, doors, windows, fuselage noses, and engine nacelles. Past work on roughness and waviness manufacturing tolerances and comparisons with more recent experiments are reviewed.
Like most heat treatable aluminum alloys, localized corrosion and stress corrosion of Al-Li-Cu alloys is strongly dependent on the nature and distribution of second phase particles. To develop a mechanistic understanding of the role of localized corrosion in the stress corrosion process, bulk samples of T(sub 1) (Al2CuLi) and a range of Al-Cu-Fe impurity phases were prepared for electrochemical experiments. Potentiodynamic polarization and galvanic couple experiments were performed in standard 0.6 M NaCl and in simulated crevice solutions to assess corrosion behavior of these particles with respect to the alpha-Al matrix. A comparison of time to failure versus applied potential using a constant load, smooth bar SCC test technique in Cl(-), Cl(-)/CrO4(2-), and Cl(-)/CO3(2-) environments shows that rapid failures are to be expected when applied potentials are more positive than the breakaway potential (E sub br) of T(sub 1) (crack tip) but less than E(sub br) of alpha-Al (crack walls). It is shown that this criterion is not satisfied in aerated Cl(-) solutions. Accordingly, SCC resistance is good. This criterion is satisfied, however, in an alkaline isolated fissure exposed to a CO2 containing atmosphere. Rapid failure induced by these fissures was recently termed preexposure embrittlement. Anodic polarization shows that the corrosion behavior of T(sub 1) is relatively unaffected in alkaline CO3(2-) environments but the alpha-Al phase is rapidly passivated. X ray diffraction of crevice walls from artificial crevices suggests that passivation of alpha-Al occurs as hydrotalcite-type compound (LiAl2(OH)6)2(+) - CO3(2-) - nH2O.
A numerical method is presented for predicting the steady inviscid flow through a ducted propfan based on a time-marching solution of the three-dimensional Euler equations. A four-stage multiple-block Runge-Kutta finite volume numerical technique, utilizing implicit residual smoothing and a blended second and fourth difference dissipation, is applied to predict the transonic flowfield about both single-rotation and counter-rotations ducted rotors. Counter-rotation predictions are based on an average-passage system of equations approach. Calculations are performed for both a single sheared H-type grid system and a multiple-block grid system incorporating a C-type grid about the cowl. Numerical results are compared with experimental data for two cases: a low speed ducted propeller and a 1.15 pressure ratio fan stage.
Impact phenomena shaped our solar system. As usual for most solar system processes, the scales are far different than we can address directly in the laboratory. Impact velocities are often much higher than we can achieve, sizes are often vastly larger, and most impacts take place in an environment where the only gravitational force is the mutual pull of the impactors. The Smooth Particle Hydrodynamics (SPH) technique has been applied in the past to the simulations of giant impacts. In these simulations, the colliding objects were so massive (at least a sizeable fraction of the Earth's mass) that material strength was negligible compared to gravity. This assumption can no longer be made when the bodies are much smaller. To this end, we have developed a 3D SPH code that includes a strength model to which we have added a von Mises yielding relation for stresses beyond the Hugoniot Elastic Limit. At the lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution. Following Grady and Kipp and Melosh et al., we introduce a state variable D ('damage'), 0 less than D less than 1, which expresses the local reduction in strength due to crack growth under tensile loading. Unfortunately for the hydrodynamics, Grady and Kipp's model predicts which fragments are the most probable ones and not the ones that are really formed. This means, for example, that if a given laboratory experiment is modeled, the fragment distribution obtained from the Grady-Kipp theory would be equivalent to a ensemble average over many realizations of the experiment. On the other hand, the hydrodynamics itself is explicit and evolves not an ensemble average but very specific fragments. Hence, there is a clear incompatibility with the deterministic nature of the hydrodynamics equations and the statistical approach of the Grady-Kipp dynamical fracture model. We remedy these shortcomings by making the incipient flaw distribution explicit, i.e., particles carry activation strains which are distributed at random with a probability of occurrence given by the Weibull distribution. If the local principal axis strain exceeds this limit, damage starts to grow. By growing explicit cracks together with statistical cracks (damage) at the sub-particle scale, we ensure that material strength and fragmentation is independent of model resolution. We tested our scheme by simulating laboratory impact experiments on basalt spheres.
We consider the steady state equations for a compressible fluid. Since we wish to solve for a range of speeds we must consider the equations in conservation form. For transonic speeds these equations are of mixed type. Hence, the usual approach is to add time derivatives to the steady state equations and then march these equations in time. One then adds a time derivative of the density to the continuity equation, a derivative of the momentum to the momentum equation and a derivative of the total energy to the energy equation. This choice is dictated by the time consistent equations. However, since we are only interested in the steady state this is not necessary. Thus we shall consider the possibility of adding a time derivative of the pressure to the continuity equation and similar modifications for the energy equation. This can then be generalized to adding combinations of time derivatives to each equation since these vanish in the steady state. When using acceleration techniques such as residual smoothing, multigrid, etc. these are applied to the pressure rather than the density. Hence, the code duplicates the behavior of the incompressible equations for low speeds.
We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.
A methodology has been devised for the design of three-dimensional hypersonic inlets which include a rectangular to elliptical shape transition. This methodology makes extensive use of inviscid streamtracing techniques to generate a smooth shape transition from a rectangular-like capture to an elliptical throat. Highly swept leading edges and a significantly notched cowl enable use of these inlets in fixed geometry configurations. The design procedure includes a three dimensional displacement thickness calculation and uses established correlations to check for boundary layer separation due to shock wave interactions. Complete details of the design procedure are presented and the characteristics of a modular inlet with rectangular to elliptical shape transition and a design point of Mach 7.1 are examined. Comparison with a classical two-dimensional inlet optimized for maximum total pressure recovery indicates that this three-dimensional inlet demonstrates good performance even well below its design point.
The next generation of ultra-low power bolometer arrays, with applications in far infrared imaging, spectroscopy and polarimetry, utilizes a superconducting bilayer as the sensing element to enable SQUID multiplexed readout. Superconducting transition edge sensors (TES s) are being produced with dual metal systems of superconductinghormal bilayers. The transition temperature (Tc) is tuned by altering the relative thickness of the superconductor with respect to the normal layer. We are currently investigating MoAu and MoCu bilayers. We have developed a dry-etching process for MoAu TES s with integrated molybdenum leads, and are working on adapting the process to MoCu. Dry etching has the advantage over wet etching in the MoAu system in that one can achieve a high degree of selectivity, greater than 10, using argon ME, or argon ion milling, for patterning gold on molybdenum. Molybdenum leads are subsequently patterned using fluorine plasma.. The dry-etch technique results in a smooth, featureless TES with sharp sidewalls, no undercutting of the Mo beneath the normal metal, and Mo leads with high critical current. The effects of individual processing parameters on the characteristics of the transition will be reported.
Recent development of high spatial resolution satellites such as IKONOS, Quickbird and Orbview enable observation of the Earth's surface with sub-meter resolution. Compared to the 30 meter resolution of Landsat 5 TM, the amount of information in the output image was dramatically increased. In this era of high spatial resolution, the estimation of spatial quality of images is gaining attention. Historically, the Modulation Transfer Function (MTF) concept has been used to estimate an imaging system's spatial quality. Sometimes classified by target shapes, various methods were developed in laboratory environment utilizing sinusoidal inputs, periodic bar patterns and narrow slits. On-orbit sensor MTF estimation was performed on 30-meter GSD Landsat4 Thematic Mapper (TM) data from the bridge pulse target as a pulse input . Because of a high resolution sensor s small Ground Sampling Distance (GSD), reasonably sized man-made edge, pulse, and impulse targets can be deployed on a uniform grassy area with accurate control of ground targets using tarps and convex mirrors. All the previous work cited calculated MTF without testing the MTF estimator's performance. In previous report, a numerical generic sensor model had been developed to simulate and improve the performance of on-orbit MTF estimating techniques. Results from the previous sensor modeling report that have been incorporated into standard MTF estimation work include Fermi edge detection and the newly developed 4th order modified Savitzky-Golay (MSG) interpolation technique. Noise sensitivity had been studied by performing simulations on known noise sources and a sensor model. Extensive investigation was done to characterize multi-resolution ground noise. Finally, angle simulation was tested by using synthetic pulse targets with angles from 2 to 15 degrees, several brightness levels, and different noise levels from both ground targets and imaging system. As a continuing research activity using the developed sensor model, this report was dedicated to MTF estimation via pulse input method characterization using the Fermi edge detection and 4th order MSG interpolation method. The relationship between pulse width and MTF value at Nyquist was studied including error detection and correction schemes. Pulse target angle sensitivity was studied by using synthetic targets angled from 2 to 12 degrees. In this report, from the ground and system noise simulation, a minimum SNR value was suggested for a stable MTF value at Nyquist for the pulse method. Target width error detection and adjustment technique based on a smooth transition of MTF profile is presented, which is specifically applicable only to the pulse method with 3 pixel wide targets.
A deconvolution method is presented for estimating input data from measured output data and a model of the dynamic process involved. The method uses an optimal Wiener filter for separating the measured data into signal and noise components, and a high-accuracy Fourier transform for inverting the model dynamics in the frequency domain. The method is an extension of optimal Fourier smoothing, and uses a technique to enhance the contrast between the signal and noise spectra in designing the Wiener filter. The deconvolution method was applied to simulation and flight test data for the purposes of removing unwanted distortions introduced by signal-conditioning filters and sensor dynamics, and for reconstructing turbulence inputs from measured sensor data. Results indicated hat the method performs well given good signal-to-noise levels and accurate models of the dynamic process.
The space environment is fraught with complex plasmas spanning a wide range of densities and temperatures. Much of space plasma research has focused on large-scale changes in the ambient plasma and fields from the Sun to the Earth’s atmosphere, driven primarily by the solar wind. Yet one particular type of space plasma, which forms from a hypervelocity impact, remains poorly understood. Hypervelocity impactors include both meteoroids and space debris. Meteoroids are naturally occurring objects in space that travel between 11 and 72.8 km/s and originate primarily from comets and asteroids. In contrast, space debris are human-made objects with speeds typically < 11 km/s. Hypervelocity impactors routinely hit spacecraft, yet the physics behind the formation of the plasma and the dynamics of its expansion remain largely unknown. The complexity of this phenomenon necessitates a research approach that includes both experimental studies and numerical simulation in order to understand the underlying physical processes that occur upon formation and expansion of the impact plasma. Our research has focused on providing the first comprehensive understanding of plasma generated by hypervelocity impacts by meteoroids and space debris on spacecraft in order to characterize the behavior of the expanding plasma and its interactions with the ambient environment. We conducted experimental campaigns at a dust accelerator facility that can accelerate particles up to 60 km/s, which is representative of meteoroid speeds, and at a light gas gun facility that can accelerate larger projectiles up to 7 km/s, which is representative of orbital debris. The experiments included plasma, optical and radio frequency (RF) sensors in order to understand the dynamics and associated RF emission resulting from hypervelocity plasmas. Furthermore, we combined the collected data with results from numerical simulations using particle-in-cell (PIC) and smoothed particle hydrodynamics (SPH) techniques. Using results from both experimental measurements and simulations, we addressed the following outstanding science questions: 1) What are the properties of hypervelocity impact plasmas? 2) What conditions cause impact plasmas to create RF emission and what is the power dependence on frequency? 3) How does the ambient environment around the spacecraft, including the background magnetic and electric fields and plasma density, influence the hypervelocity impact plasmas?
The fixed-point smoothing problem of sequentially correlated processes is considered. The filtering technique of such processes is extended to smoothing problems by simultaneously estimating the state and the accumulated contribution of process noise. Suboptimal filtering and smoothing are also discussed when data noise is regarded white in designing filters although it is actually sequentially correlated.-
The existence of two- and three-finger grasps in the presence of arbitrarily small friction is shown for two- and three-dimensional smooth objects using a simple technique. No convexity of the objects is assumed. The existence of finger gaits for rotating a planar object using three and four fingers is proved. Additional results for smooth convex objects which describe two different finger gaits using four fingers are presented.
In this paper, evidence is presented for the existence of multiple solutions of Reynolds-averaged Navier–Stokes equations with the one-equation Spalart–Allmaras and two-equation Wilcox $k-ω$ turbulence models on fixed grids in three dimensions and how they were obtained is described. The two major configurations considered are an “academic” extruded two-dimensional airfoil geometry and the trap wing test case. The observed appearance of the multiple solutions seems to be closely related to smooth body separation (sometimes massive) routinely observed in flows over high-lift configurations, especially near stall angles of attack. The results are obtained and cross-verified with two stabilized finite-element codes (streamwise upwind Petrov–Galerkin), which provide residual converged results for complex flows with second-order discretizations. In the paper, the ways multiple solutions have been obtained are described, including such obvious ones as providing a different initial guess to the steady-state solver as well as somewhat unexpected (in this context) techniques of using implicit residual smoothing while time marching to steady state. The phenomenon of the so-called pseudosolutions is also discussed, which is loosely defined as solutions to the discrete system of equations having “sharp” convergence behavior (sometimes up to six to seven orders of the relative residual reduction), which fail to achieve the stronger, machine-zero convergence criterion. Some numerical observations are also presented on the sensitivity of the obtained multiple solutions to the discretization and grid perturbations.