Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordered statistics”

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,153 records · Page 64

Characterization of surficial geologic units on Venus from Pioneer Venus radar data: A progress report

A classification database using the reflectivity (derived from the altimetry data), rms slope, and the first principal component of altimetry and topographic slope is presented. The resultant clustered data is examined qualitatively as well as quantitatively, to establish the statistical integrity of each cluster by use of an interactive, ternary plotting algorithm. This algorithm plots, for a cluster, the position of each of its pixels within a ternary diagram whose apices represent reflectivity, rms slope, and the first principal component. The digital values in these three databases are normalized such that unity is represented by a value of 255 in each database. The frequencies of each plotted point within the ternary diagram are recorded in order to establish the mode of each cluster. The pixels of each cluster are displayed as one separate color; their ternary plot will show not only the interrelations between clusters, but also the presence of any anomalous points within a cluster. Existing lunar and terrestrial analog radar data is used to establish fields within this ternary diagram that are indicative of as many different geologic materials and tectonics settings as possible. The resultant fields are used to determine empirically the geologic significance of the clusters resulting from the cluster analysis.

Davis, P. A.↗

Activity associated with coronal mass ejections at solar minimum - SMM observations from 1984-1986

Seventy-three coronal mass ejections (CMEs) observed by the coronagraph aboard SMM between 1984 and 1986 were examined in order to determine the distribution of various forms of solar activity that were spatially and temporally associated with mass ejections during solar minimum phase. For each coronal mass ejection a speed was measured, and the departure time of the transient from the lower corona estimated. Other forms of solar activity that appeared within 45 deg longitude and 30 deg latitude of the mass ejection and within +/-90 min of its extrapolated departure time were explored. The statistical results of the analysis of these 73 CMEs are presented, and it is found that slightly less than half of them were infrequently associated with other forms of solar activity. It is suggested that the distribution of the various forms of activity related to CMEs does not change at different phases of the solar cycle. For those CMEs with associations, it is found that eruptive prominences and soft X-rays were the most likely forms of activity to accompany the appearance of mass ejections.

St. Cyr, O. C.↗

IUEAGN: A database of ultraviolet spectra of active galactic nuclei

In 13 years of operation, IUE has gathered approximately 5000 spectra of almost 600 Active Galactic Nuclei (AGN). In order to undertake AGN studies which require large amounts of data, we are consistently reducing this entire archive and creating a homogeneous, easy-to-use database. First, the spectra are extracted using the Optimal extraction algorithm. Continuum fluxes are then measured across predefined bands, and line fluxes are measured with a multi-component fit. These results, along with source information such as redshifts and positions, are placed in the IUEAGN relational database. Analysis algorithms, statistical tests, and plotting packages run within the structure, and this flexible database can accommodate future data when they are released. This archival approach has already been used to survey line and continuum variability in six bright Seyfert 1s and rapid continuum variability in 14 blazars. Among the results that could only be obtained using a large archival study is evidence that blazars show a positive correlation between degree of variability and apparent luminosity, while Seyfert 1s show an anti-correlation. This suggests that beaming dominates the ultraviolet properties for blazars, while thermal emission from an accretion disk dominates for Seyfert 1s. Our future plans include a survey of line ratios in Seyfert 1s, to be fitted with photoionization models to test the models and determine the range of temperatures, densities and ionization parameters. We will also include data from IRAS, Einstein, EXOSAT, and ground-based telescopes to measure multi-wavelength correlations and broadband spectral energy distributions.

Pike, G.↗

Premixed Turbulent Flame Propagation in Microgravity

A combined numerical-experimental study has been carried out to investigate the structure and propagation characteristics of turbulent premixed flames with and without the influence of buoyancy. Experimentally, the premixed flame characteristics are studied in the wrinkled regime using a Couette flow facility and an isotropic flow facility in order to resolve the scale of flame wrinkling. Both facilities were chosen for their ability to achieve sustained turbulence at low Reynolds number. This implies that conventional diagnostics can be employed to resolve the smallest scales of wrinkling. The Couette facility was also built keeping in mind the constraints imposed by the drop tower requirements. Results showed that the flow in this Couette flow facility achieves full-developed turbulence at low Re and all turbulence statistics are in good agreement with past measurements on large-scale facilities. Premixed flame propagation studies were then carried out both using the isotropic box and the Couette facility. Flame imaging showed that fine scales of wrinkling occurs during flame propagation. Both cases in Ig showed significant buoyancy effect. To demonstrate that micro-g can remove this buoyancy effect, a small drop tower was built and drop experiments were conducted using the isotropic box. Results using the Couette facility confirmed the ability to carry out these unique reacting flow experiments at least in 1g. Drop experiments at NASA GRC were planned but were not completed due to termination of this project.

Menon, Suresh↗

Experiments on Nucleation in Different Flow Regimes

The vast majority of metallic engineering materials are solidified from the liquid phase. Understanding the solidification process is essential to control microstructure, which in turn, determines the properties of materials. The genesis of solidification is nucleation, where the first stable solid forms from the liquid phase. Nucleation kinetics determine the degree of undercooling and phase selection. As such, it is important to understand nucleation phenomena in order to control solidification or glass formation in metals and alloys. Early experiments in nucleation kinetics were accomplished by droplet dispersion methods [1-6]. Dilitometry was used by Turnbull and others, and more recently differential thermal analysis and differential scanning calorimetry have been used for kinetic studies. These techniques have enjoyed success; however, there are difficulties with these experiments. Since materials are dispersed in a medium, the character of the emulsion/metal interface affects the nucleation behavior. Statistics are derived from the large number of particles observed in a single experiment, but dispersions have a finite size distribution which adds to the uncertainty of the kinetic determinations. Even though temperature can be controlled quite well before the onset of nucleation, the release of the latent heat of fusion during nucleation of particles complicates the assumption of isothermality during these experiments. Containerless processing has enabled another approach to the study of nucleation kinetics [7]. With levitation techniques it is possible to undercool one sample to nucleation repeatedly in a controlled manner, such that the statistics of the nucleation process can be derived from multiple experiments on a single sample. The authors have fully developed the analysis of nucleation experiments on single samples following the suggestions of Skripov [8]. The advantage of these experiments is that the samples are directly observable. The nucleation temperature can be measured by noncontact optical pyrometry, the mass of the sample is known, and post processing analysis can be conducted on the sample. The disadvantages are that temperature measurement must have exceptionally high precision, and it is not possible to isolate specific heterogeneous sites as in droplet dispersions.

Bayuzick, Robert J.↗

Heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(s)

The heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) (eq 1) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(S) (eq 2) were investigated over the temperature range 223-296 K in a flow-tube reactor coupled to a quadrupole mass spectrometer. Either a chemical ionization mass spectrometer (CIMS) or an electron-impact ionization mass spectrometer (EIMS) was used to provide suitable detection sensitivity and selectivity. In order to mimic atmospheric conditions, partial pressures of HNO3 and N2O5 in the range 6 x 10(exp -8) - 2 x 10(exp -6) Torr were used. Granule sizes and surface roughness of the solid NaCl substrates were determined by using a scanning electron microscope. For dry NaCl substrates, decay rates of HNO3 were used to obtain gamma(1) = 0.013 +/- 0.004 (1sigma) at 296 K and > 0.008 at 223 K, respectively. The error quoted is the statistical error. After all corrections were made, the overall error, including systematic error, was estimated to be about a factor of 2. HCl was found to be the sole gas-phase product of reaction 1. The mechanism changed from heterogeneous reaction to predominantly physical adsorption when the reactor was cooled from 296 to 223 K. For reaction 2 using dry salts, gamma(2) was found to be less than 1.0 x 10(exp -4) at both 223 and 296 K. The gas-phase reaction product was identified as ClNO2 in previous studies using an infrared spectrometer. An enhancement in reaction probability was observed if water was not completely removed from salt surfaces, probably due to the reaction of N2O5(g) + H2O(s) yields 2HNO3(g). Our results are compared with previous literature values obtained using different experimental techniques and conditions. The implications of the present results for the enhancement of the hydrogen chloride column density in the lower stratosphere after the El Chichon volcanic eruption and for the chemistry of HCl and HNO3 in the marine troposphere are discussed.

Leu, Ming-Taun↗

Monitoring of MLS Measurements for Ozone Data Assimilation

The ozone data assimilation system at the NASA Goddard Space Flight Center Data Assimilation Office (DAO) has operationally provided near real time global three dimensional analyzed ozone fields since December 1999. Presently, 4 system assimilates both total column ozone and profile measurements from the NOAA - 16 SBUV/2 instrument into an off-line transport model using a physical space statistical analysis scheme (PS AS). Recently, this system was modified to investigate a possible impact from including MLS profile observations. Over a period from December 13, 1991 to March 1, 1992 the system assimilated both NOAA 11 SBUV/2 profile measurements and TOMS total column measurements. Short term ozone forecasts were produced and compared to MLS observations. The resulting observed -minus -forecast (O-F) residuals were studied regionally (by both latitude bands and pressure levels) in order to explore the potential impact of MLS observations on the assimilation and sensitivity to the UARS yaw maneuver and other instrument and algorithm characteristics. In addition, MLS profiles were compared to nearest neighbor SBUV and analysis profiles. In the tropics, analyzed ozone values tend to be too high from 1-3 hPa and too low from 10-20 hPa. MLS measurements in these regions tend be lower and higher (respectively) than forecast ozone values. This indicates that assimilating MLS measurements may improve, analysis results in these areas. Further, nearest neighbor profile comparisons in the southern high latitudes indicate that laminar features present in the analysis but not present in SBUV measurements, are also present in MLS measurements. This, together with the availability of MLS measurements in polar night regions indicates that assimilating MLS measurements in polar night regions indicates that assimilating MLS measurements may improve analyzed ozone values in high latitudes.

Winslow, Nathan↗

Parametric Studies of Flow Separation using Air Injection

Boundary Layer separation causes the airfoil to stall and therefore imposes dramatic performance degradation on the airfoil. In recent years, flow separation control has been one of the active research areas in the field of aerodynamics due to its promising performance improvements on the lifting device. These active flow separation control techniques include steady and unsteady air injection as well as suction on the airfoil surface etc. This paper will be focusing on the steady and unsteady air injection on the airfoil. Although wind tunnel experiments revealed that the performance improvements on the airfoil using injection techniques, the details of how the key variables such as air injection slot geometry and air injection angle etc impact the effectiveness of flow separation control via air injection has not been studied. A parametric study of both steady and unsteady air injection active flow control will be the main objective for this summer. For steady injection, the key variables include the slot geometry, orientation, spacing, air injection velocity as well as the injection angle. For unsteady injection, the injection frequency will also be investigated. Key metrics such as lift coefficient, drag coefficient, total pressure loss and total injection mass will be used to measure the effectiveness of the control technique. A design of experiments using the Box-Behnken Design is set up in order to determine how each of the variables affects each of the key metrics. Design of experiment is used so that the number of experimental runs will be at minimum and still be able to predict which variables are the key contributors to the responses. The experiments will then be conducted in the 1ft by 1ft wind tunnel according to the design of experiment settings. The data obtained from the experiments will be imported into JMP, statistical software, to generate sets of response surface equations which represent the statistical empirical model for each of the metrics as a function of the key variables. Next, the variables such as the slot geometry can be optimized using the build-in optimizer within JMP. Finally, a wind tunnel testing will be conducted using the optimized slot geometry and other key variables to verify the empirical statistical model. The long term goal for this effort is to assess the impacts of active flow control using air injection at system level as one of the task plan included in the NASAs URETI program with Georgia Institute of Technology.

Zhang, Wei↗

Mixing of a passive scalar in isotropic and sheared homogeneous turbulence

In order to calculate the velocity and scalar fields, the three dimensional, time-dependent equations of motion and the diffusion equation were solved numerically. The following cases were treated: isotropic, homogeneous turbulence with decay of a passive scalar; and homogeneous turbulent shear flow with a passive scalar whose mean varies linearly in the spanwise direction. The solutions were obtained at relatively low Reynolds numbers so that all of the turbulent scales could be resolved without modeling. Turbulent statistics such as integral length scales, Taylor microscales, Kolmogorov length scale, one- and two-point correlations of velocity-velocity and velocity-scalar, turbulent Prandtl/Schmidt number, r.m.s. values of velocities, the scalar quantity and pressure, skewness, decay rates, and decay exponents were calculated. The results are compared with the available expermental results, and good agreement is obtained.

Shirani, E.↗

Edge detection - Image-plane versus digital processing

To optimize edge detection with the familiar Laplacian-of-Gaussian operator, it has become common to implement this operator with a large digital convolution mask followed by some interpolation of the processed data to determine the zero crossings that locate edges. It is generally recognized that this large mask causes substantial blurring of fine detail. It is shown that the spatial detail can be improved by a factor of about four with either the Wiener-Laplacian-of-Gaussian filter or an image-plane processor. The Wiener-Laplacian-of-Gaussian filter minimizes the image-gathering degradations if the scene statistics are at least approximately known and also serves as an interpolator to determine the desired zero crossings directly. The image-plane processor forms the Laplacian-of-Gaussian response by properly combining the optical design of the image-gathering system with a minimal three-by-three lateral-inhibitory processing mask. This approach, which is suggested by Marr's model of early processing in human vision, also reduces data processing by about two orders of magnitude and data transmission by up to an order of magnitude.

Huck, Friedrich O.↗

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

Developing Conceptual Hypersonic Airbreathing Engines Using Design of Experiments Methods

Designing a hypersonic vehicle is a complicated process due to the multi-disciplinary synergy that is required. The greatest challenge involves propulsion-airframe integration. In the past, a two-dimensional flowpath was generated based on the engine performance required for a proposed mission. A three-dimensional CAD geometry was produced from the two-dimensional flowpath for aerodynamic analysis, structural design, and packaging. The aerodynamics, engine performance, and mass properties arc inputs to the vehicle performance tool to determine if the mission goals were met. If the mission goals were not met, then a flowpath and vehicle redesign would begin. This design process might have to be performed several times to produce a "closed" vehicle. This paper will describe an attempt to design a hypersonic cruise vehicle propulsion flowpath using a Design of' Experiments method to reduce the resources necessary to produce a conceptual design with fewer iterations of the design cycle. These methods also allow for more flexible mission analysis and incorporation of additional design constraints at any point. A design system was developed using an object-based software package that would quickly generate each flowpath in the study given the values of the geometric independent variables. These flowpath geometries were put into a hypersonic propulsion code and the engine performance was generated. The propulsion results were loaded into statistical software to produce regression equations that were combined with an aerodynamic database to optimize the flowpath at the vehicle performance level. For this example, the design process was executed twice. The first pass was a cursory look at the independent variables selected to determine which variables are the most important and to test all of the inputs to the optimization process. The second cycle is a more in-depth study with more cases and higher order equations representing the design space.

Ferlemann, Shelly M.↗

Detailed Microstructural Characterization of the Disk Alloy ME3

The advanced powder metallurgy disk alloy ME3 was designed using statistical screening and optimization of composition and processing variables in the NASA/General Electric/Pratt & Whitney HSR/EPM disk program to have extended durability for large disks at maximum temperatures of 600 to 700 C. Scaled-up disks of this alloy were then produced at the conclusion of that program to demonstrate these properties in realistic disk shapes. The objective of the present study was to assess the microstructural characteristics of these ME3 disks at two consistent locations, in order to enable estimation of the variations in microstructure across each disk and across several disks of this advanced alloy. Scaled-up disks processed in the HSR/EPM Compressor/Turbine Disk program had been sectioned, machined into specimens, and tested in tensile, creep, fatigue, and fatigue crack growth tests by NASA Glenn Research Center, in cooperation with General Electric Engine Company and Pratt & Whitney Aircraft Engines. For this study, microstructures of grip sections from tensile specimens in the bore and rim were evaluated from these disks. The major and minor phases were identified and quantified using transmission electron microscopy (TEM). Particular attention was directed to the .' precipitates, which along with grain size can predominantly control the mechanical properties of superalloy disks.

Gabb, Timothy P.↗

Fast Solution in Sparse LDA for Binary Classification

An algorithm that performs sparse linear discriminant analysis (Sparse-LDA) finds near-optimal solutions in far less time than the prior art when specialized to binary classification (of 2 classes). Sparse-LDA is a type of feature- or variable- selection problem with numerous applications in statistics, machine learning, computer vision, computational finance, operations research, and bio-informatics. Because of its combinatorial nature, feature- or variable-selection problems are NP-hard or computationally intractable in cases involving more than 30 variables or features. Therefore, one typically seeks approximate solutions by means of greedy search algorithms. The prior Sparse-LDA algorithm was a greedy algorithm that considered the best variable or feature to add/ delete to/ from its subsets in order to maximally discriminate between multiple classes of data. The present algorithm is designed for the special but prevalent case of 2-class or binary classification (e.g. 1 vs. 0, functioning vs. malfunctioning, or change versus no change). The present algorithm provides near-optimal solutions on large real-world datasets having hundreds or even thousands of variables or features (e.g. selecting the fewest wavelength bands in a hyperspectral sensor to do terrain classification) and does so in typical computation times of minutes as compared to days or weeks as taken by the prior art. Sparse LDA requires solving generalized eigenvalue problems for a large number of variable subsets (represented by the submatrices of the input within-class and between-class covariance matrices). In the general (fullrank) case, the amount of computation scales at least cubically with the number of variables and thus the size of the problems that can be solved is limited accordingly. However, in binary classification, the principal eigenvalues can be found using a special analytic formula, without resorting to costly iterative techniques. The present algorithm exploits this analytic form along with the inherent sequential nature of greedy search itself. Together this enables the use of highly-efficient partitioned-matrix-inverse techniques that result in large speedups of computation in both the forward-selection and backward-elimination stages of greedy algorithms in general.

Moghaddam, Baback↗

Aeroelastic Uncertainty Quantification Studies Using the S4T Wind Tunnel Model

This paper originates from the joint efforts of an aeroelastic study team in the Applied Vehicle Technology Panel from NATO Science and Technology Organization, with the Task Group number AVT-191, titled "Application of Sensitivity Analysis and Uncertainty Quantification to Military Vehicle Design." We present aeroelastic uncertainty quantification studies using the SemiSpan Supersonic Transport wind tunnel model at the NASA Langley Research Center. The aeroelastic study team decided treat both structural and aerodynamic input parameters as uncertain and represent them as samples drawn from statistical distributions, propagating them through aeroelastic analysis frameworks. Uncertainty quantification processes require many function evaluations to asses the impact of variations in numerous parameters on the vehicle characteristics, rapidly increasing the computational time requirement relative to that required to assess a system deterministically. The increased computational time is particularly prohibitive if high-fidelity analyses are employed. As a remedy, the Istanbul Technical University team employed an Euler solver in an aeroelastic analysis framework, and implemented reduced order modeling with Polynomial Chaos Expansion and Proper Orthogonal Decomposition to perform the uncertainty propagation. The NASA team chose to reduce the prohibitive computational time by employing linear solution processes. The NASA team also focused on determining input sample distributions.

Nikbay, Melike↗

Theoretical and observational analysis of spacecraft fields

In order to investigate the dipolar and nondipolar contributions of spacecraft magnetic fields, a simple magnetic field model is proposed. This model consists of a number N of randomly oriented dipoles of strength M sub k in a given volume. Two sets of formulas are presented that give the rms multipole field components, first, for isotropic orientations of the dipoles at given positions and second, for isotropic orientations of the dipoles distributed uniformly throughout a cube or sphere. The statistical results for a cube (8 cu m in size) together with individual examples computed numerically show the following features. Beyond about 2- to 3-m distance from the center of the cube the field is dominated by an equivalent dipole. The magnitude of the magnetic moment of the dipolar part is approximated by N to the 1/2 power times M for equal magnetic moments or generally by the Pythagorean sum of the dipole moments. The radial component tend to be greater than either of the transverse components for the dipole portion as well as for the nondipolar field contributions.

Neubauer, F. M.↗

On time-dependent radiative transfer

An integral equation is developed for application to time-dependent laboratory experiments in which partial redistribution effects are important. The equation of transport with the Heasly-Kneer (1976) emission coefficient and the equation of statistical equilibrium lead to a time-dependent redistribution function containing an absorption - reemission term which decays exponentially in time and a scattering term which is instantaneous. This integral equation does not agree with an equation written by Payne et al. (1974) that has been used to compare theory with experiments. The difference between the Payne equation and the equation developed here needs to be examined in detail, since it might under some circumstances be on the same order as the difference between partial and complete redistribution.

Streater, A.↗

Summertime Coincident Observations of Ice Water Path in the Visible/Near-IR, Radar, and Microwave Frequencies

Accurate representation of the physical and radiative properties of clouds in climate models continues to be a challenge. At present, both remote sensing observations and modeling of microphysical properties of clouds rely heavily on parameterizations or assumptions on particle size distribution (PSD) and cloud phase. In this study, we compare Ice Water Path (IWP), an important physical and radiative property that provides the amount of ice present in a cloud column, using measurements obtained via three different retrieval strategies. The datasets we use in this study include Visible/Near-IR IWP from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument flying aboard the Aqua satellite, Radar-only IWP from the CloudSat instrument operating at 94 GHz, and NOAA/NESDIS operational IWP from the 89 and 157 GHz channels of the Microwave Humidity Sounder (MHS) instrument flying aboard the NOAA-18 satellite. In the Visible/Near-IR, IWP is derived from observations of optical thickness and effective radius. CloudSat IWP is determined from measurements of cloud backscatter and assumed PSD. MHS IWP retrievals depend on scattering measurements at two different, non-water absorbing channels, 89 and 157 GHz. In order to compare IWP obtained from these different techniques and collected at different vertical and horizontal resolutions, we examine summertime cases in the tropics (30S - 30N) when all 3 satellites are within 4 minutes of each other (approximately 1500 km). All measurements are then gridded to a common 15 km x 15 km box determined by MHS. In a grid box comparison, we find CloudSat to report the highest IWP followed by MODIS, followed by MHS. In a statistical comparison, probability density distributions show MHS with the highest frequencies at IWP of 100-1000 g/m(exp 2) and CloudSat with the longest tail reporting IWP of several thousands g/m(exp 2). For IWP greater than 30 g/m(exp 2), MODIS is consistently higher than CloudSat, and it is higher at the lower IWPs but lower at the higher IWPs that overlap with MHS. Some of these differences can be attributed to the limitations of the measuring techniques themselves, but some can result from the assumptions made in the algorithms that generate the IWP product. We investigate this issue by creating categories based on various conditions such as cloud type, precipitation presence, underlying liquid water content, and surface type (land vs. ocean) and by comparing the performance of the IWP products under each condition.

Pittman, Jasna V.↗