Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “clustering coefficient”

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 91 records · Page 5

Transonic Drag Prediction Using an Unstructured Multigrid Solver

This paper summarizes the results obtained with the NSU-3D unstructured multigrid solver for the AIAA Drag Prediction Workshop held in Anaheim, CA, June 2001. The test case for the workshop consists of a wing-body configuration at transonic flow conditions. Flow analyses for a complete test matrix of lift coefficient values and Mach numbers at a constant Reynolds number are performed, thus producing a set of drag polars and drag rise curves which are compared with experimental data. Results were obtained independently by both authors using an identical baseline grid and different refined grids. Most cases were run in parallel on commodity cluster-type machines while the largest cases were run on an SGI Origin machine using 128 processors. The objective of this paper is to study the accuracy of the subject unstructured grid solver for predicting drag in the transonic cruise regime, to assess the efficiency of the method in terms of convergence, cpu time, and memory, and to determine the effects of grid resolution on this predictive ability and its computational efficiency. A good predictive ability is demonstrated over a wide range of conditions, although accuracy was found to degrade for cases at higher Mach numbers and lift values where increasing amounts of flow separation occur. The ability to rapidly compute large numbers of cases at varying flow conditions using an unstructured solver on inexpensive clusters of commodity computers is also demonstrated.

Mavriplis, D. J.↗

Directivity of planar array feeds for satellite reflector applications

Array-fed reflector antennas are used extensively in today's contour and multiple beam satellite antennas. To determine the directivity of these antennas theoretically, the total radiated power of the array feed must be accurately computed. In this paper, a closed-form expression for the radiated power is obtained for arrays with the well-known (cos theta q-type element patterns. The formulation is general, and takes into account polarization, nonsymmetric E- and H-plane element patterns, nonuniform element spacings and arbitrary complex excitation coefficients. Selective numerical data are presented to demonstrate the usefulness of this result. Comparisons are made with the available results obtained using direct numerical integration techniques, and with other available data based on less general formulations. Excellent agreement is observed for all cases. In particular, data are presented on the directivity of seven-element cluster feeds used in multiple beam designs, and an array feed for producing a contour beam covering the United States Eastern Time Zone (ETZ).

Rahmat-Samii, Y.↗

An open-shell restricted Hartree-Fock perturbation theory based on symmetric spin orbitals

A new open-shell perturbation theory is formulated in terms of symmetric spin orbitals. Only one set of spatial orbitals is required, thereby reducing the number of independent coefficients in the perturbed wavefunctions. For second order, the computational cost is shown to be similar to a closed-shell calculation. This formalism is therefore more efficient than the recently developed RMP, ROMP or RMP-MBPT theories. The perturbation theory described herein was designed to have a close correspondence with our recently proposed coupled-cluster theory based on symmetric spin orbitals. The first-order wavefunction contains contributions from only doubly excited determinants. Equilibrium structures and vibrational frequencies determined from second-order perturbation theory are presented for OH, NH, CH, 02, NH2 and CH2.

Lee, Timothy J.↗

Spherical Harmonic Analysis of Particle Velocity Distribution Function: Comparison of Moments and Anisotropies using Cluster Data

This paper presents a spherical harmonic analysis of the plasma velocity distribution function using high-angular, energy, and time resolution Cluster data obtained from the PEACE spectrometer instrument to demonstrate how this analysis models the particle distribution function and its moments and anisotropies. The results show that spherical harmonic analysis produced a robust physical representation model of the velocity distribution function, resolving the main features of the measured distributions. From the spherical harmonic analysis, a minimum set of nine spectral coefficients was obtained from which the moment (up to the heat flux), anisotropy, and asymmetry calculations of the velocity distribution function were obtained. The spherical harmonic method provides a potentially effective "compression" technique that can be easily carried out onboard a spacecraft to determine the moments and anisotropies of the particle velocity distribution function for any species. These calculations were implemented using three different approaches, namely, the standard traditional integration, the spherical harmonic (SPH) spectral coefficients integration, and the singular value decomposition (SVD) on the spherical harmonic methods. A comparison among the various methods shows that both SPH and SVD approaches provide remarkable agreement with the standard moment integration method.

Gurgiolo, Chris↗

Origin, Transport, and Vertical Distribution of Atmospheric Polluntants over the Northern Sourth China Sea During the 7-SEAS-Dongsha Experiment

During the spring of 2010, comprehensive in situ measurements were made for the first time on a small atoll (Dongsha Island) in the northern South China Sea (SCS), a key region of the 7-SEAS (the Seven South East Asian Studies) program. This paper focuses on characterizing the source origins, transport processes, and vertical distributions of the Asian continental outflows over the region, using measurements including mass concentration, optical properties, hygroscopicity, and vertical distribution of the aerosol particles, as well as the trace gas composition. Cluster analysis of backward trajectories classified 52% of the air masses arriving at ground level of Dongsha Island as having a continental origin, mainly from northern China to the northern SCS, passing the coastal area and being confined in the marine boundary layer (0-0.5 km). Compared to aerosols of oceanic origin, the fine mode continental aerosols have a higher concentration, extinction coefficient, and single-scattering albedo at 550 nm (i.e., 19 vs. 14 microg per cubic meter in PM(sub 2.5); 77 vs. 59 M per meter in beta(sub e); and 0.94 vs. 0.90 in omega, respectively). These aerosols have a higher hygroscopicity (f at 85% RH = 2.1) than those in the upwind inland regions, suggesting that the aerosols transported to the northern SCS were modified by the marine environment. In addition to the near-surface aerosol transport, a significant upper-layer (3-4 km) transport of biomass-burning aerosols was observed. Our results suggest that emissions from both China and Southeast Asia could have a significant impact on the aerosol loading and other aerosol properties over the SCS. Furthermore, the complex vertical distribution of aerosols-coinciding-with-clouds has implications for remote-sensing observations and aerosol-cloud-radiation interactions.

pollutants↗

Flood Inundation Mapping in the Logone Floodplain from Multi Temporal Landsat ETM+Imagery

Yearly flooding in the Logone floodplain makes an impact on agricultural, pastoral, and fishery systems in the Lake Chad Basin. Since the flooding extent and depth are highly variable, flood inundation mapping helps us make better use of water resources and prevent flood hazards in the Logone floodplain. The flood maps are generated from 33 multi temporal Landsat Enhanced Thematic Mapper Plus (ETM+) during three years 2006 to 2008. Flooded area is classified using a short-wave infrared band whereas open water is classified by Iterative Self-organizing Data Analysis (ISODATA) clustering. The maximum flooding extent in the study area increases up to approximately 5.8K km2 in late October 2008. The study also provides strong correlation of the flooding extents with water height variations in both the floodplain and the river based on a second polynomial regression model. The water heights are from ENIVSAT altimetry in the floodplain and gauge measurements in the river. Coefficients of determination between flooding extents and water height variations are greater than 0.91 with 4 to 36 days in phase lag. Floodwater drains back to the river and to the northeast during the recession period in December and January. The study supports understanding of the Logone floodplain dynamics in detail of spatial pattern and size of the flooding extent and assists the flood monitoring and prediction systems in the catchment.

Jung, Hahn Chul↗

Sea surface velocities from visible and infrared multispectral atmospheric mapping sensor imagery

High resolution (100 m), sequential Multispectral Atmospheric Mapping Sensor (MAMS) images were used in a study to calculate advective surface velocities using the Maximum Cross Correlation (MCC) technique. Radiance and brightness temperature gradient magnitude images were formed from visible (0.48 microns) and infrared (11.12 microns) image pairs, respectively, of Chandeleur Sound, which is a shallow body of water northeast of the Mississippi delta, at 145546 GMT and 170701 GMT on 30 Mar. 1989. The gradient magnitude images enhanced the surface water feature boundaries, and a lower cutoff on the gradient magnitudes calculated allowed the undesirable sunglare and backscatter gradients in the visible images, and the water vapor absorption gradients in the infrared images, to be reduced in strength. Requiring high (greater than 0.4) maximum cross correlation coefficients and spatial coherence of the vector field aided in the selection of an optimal template size of 10 x 10 pixels (first image) and search limit of 20 pixels (second image) to use in the MCC technique. Use of these optimum input parameters to the MCC algorithm, and high correlation and spatial coherence filtering of the resulting velocity field from the MCC calculation yielded a clustered velocity distribution over the visible and infrared gradient images. The velocity field calculated from the visible gradient image pair agreed well with a subjective analysis of the motion, but the velocity field from the infrared gradient image pair did not. This was attributed to the changing shapes of the gradient features, their nonuniqueness, and large displacements relative to the mean distance between them. These problems implied a lower repeat time for the imagery was needed in order to improve the velocity field derived from gradient imagery. Suggestions are given for optimizing the repeat time of sequential imagery when using the MCC method for motion studies. Applying the MCC method to the infrared brightness temperature imagery yielded a velocity field which did agree with the subjective analysis of the motion and that derived from the visible gradient imagery. Differences between the visible and infrared derived velocities were 14.9 cm/s in speed and 56.7 degrees in direction. Both of these velocity fields also agreed well with the motion expected from considerations of the ocean bottom topography and wind and tidal forcing in the study area during the 2.175 hour time interval.

Pope, P. A.↗

Electron-temperature dependence of dissociative recombination of electrons with CO(+)-(CO)n-series ions

A microwave afterglow mass spectrometer apparatus is used to determine the dependence on electron temperature T sub e of the recombination coefficients alpha sub n of the dimer and trimer ions of the series CO+.(CO) sub n. It is found that alpha sub 1 = (1.3 + or - 0.3)x 0.000001 (T sub e(K)/300) to the -0.34; and alpha sub 2 = (1.9 + or - 0.4)x 0.000001 (T sub e(K)/300) to the -0.33 cu cm/sec. These dependences on T sub e are quite different from those obtained previously for polar-cluster ions of the hydronium and ammonium series but are similar to that for simple diatomic ions.

Whitaker, M.↗

Statistical methods for astronomical data with upper limits. II - Correlation and regression

Statistical methods for calculating correlations and regressions in bivariate censored data where the dependent variable can have upper or lower limits are presented. Cox's regression and the generalization of Kendall's rank correlation coefficient provide significant levels of correlations, and the EM algorithm, under the assumption of normally distributed errors, and its nonparametric analog using the Kaplan-Meier estimator, give estimates for the slope of a regression line. Monte Carlo simulations demonstrate that survival analysis is reliable in determining correlations between luminosities at different bands. Survival analysis is applied to CO emission in infrared galaxies, X-ray emission in radio galaxies, H-alpha emission in cooling cluster cores, and radio emission in Seyfert galaxies.

Isobe, T.↗

Transverse Dimensions of Chorus in the Source Region

We report measurement of whistler-mode chorus by the four Cluster spacecraft at close separations. We focus our analysis on the generation region close to the magnetic equatorial plane at a radial distance of 4.4 Earth's radii. We use both linear and rank correlation analysis to define perpendicular dimensions of the sources of chorus elements below one half of the electron cyclotron frequency. Correlation is significant throughout the range of separation distances of 60-260 km parallel to the field line and 7-100 km in the perpendicular plane. At these scales, the correlation coefficient is independent for parallel separations, and decreases with perpendicular separation. The observations are consistent with a statistical model of the source region assuming individual sources as gaussian peaks of radiated power with a common half-width of 35 km perpendicular to the magnetic field. This characteristic scale is comparable to the wavelength of observed waves.

Santolik, O.↗

Oxygen tracer diffusion in single-crystal alumina

Oxygen tracer diffusion coefficients are determined in single-crystal alumina samples with differing dopant levels using the gas-exchange technique. The diffusion direction is parallel to the c-axis and the ambient PO2 is 1 atm (100,000 Pa) for all experiments except a single run with a low PO2, approximately 10 to the -15th atm (10 to the -10th Pa) produced by a CO/CO2 mixture. The diffusion is insensitive to both impurities and ambient PO2. The insensitivities are discussed in terms of point-defect clustering. Prior tracer studies are compared and discussed.

Cawley, James D.↗

Improved Hydrology over Peatlands in a Global Land Modeling System

Peatlands of the Northern Hemisphere represent an important carbon pool that mainly accumulated since the last ice age under permanently wet conditions in specific geological and climatic settings. The carbon balance of peatlands is closely coupled to water table dynamics. Consequently, the future carbon balance over peatlands is strongly dependent on how hydrology in peatlands will react to changing boundary conditions, e.g. due to climate change or regional water level drawdown of connected aquifers or streams. Global land surface modeling over organic-rich regions can provide valuable global-scale insights on where and how peatlands are in transition due to changing boundary conditions. However, the current global land surface models are not able to reproduce typical hydrological dynamics in peatlands well. We implemented specific structural and parametric changes to account for key hydrological characteristics of peatlands into NASA's GEOS-5 Catchment Land Surface Model (CLSM, Koster et al. 2000). The main modifications pertain to the modeling of partial inundation, and the definition of peatland-specific runoff and evapotranspiration schemes. We ran a set of simulations on a high performance cluster using different CLSM configurations and validated the results with a newly compiled global in-situ dataset of water table depths in peatlands. The results demonstrate that an update of soil hydraulic properties for peat soils alone does not improve the performance of CLSM over peatlands. However, structural model changes for peatlands are able to improve the skill metrics for water table depth. The validation results for the water table depth indicate a reduction of the bias from 2.5 to 0.2 m, and an improvement of the temporal correlation coefficient from 0.5 to 0.65, and from 0.4 to 0.55 for the anomalies. Our validation data set includes both bogs (rain-fed) and fens (ground and/or surface water influence) and reveals that the metrics improved less for fens. In addition, a comparison of evapotranspiration and soil moisture estimates over peatlands will be presented, albeit only with limited ground-based validation data. We will discuss strengths and weaknesses of the new model by focusing on time series of specific validation sites.

Bechtold, M.↗

Fresh look at floating shock fitting

A fast implicit upwind procedure for the two-dimensional Euler equations is described that allows accurate computations of shocked flows on nonadapted meshes. Away from shocks, the second-order accurate upwinding is based on the split-coefficient-matrix (SCM) method. In the presence of shocks, the difference stencils are modified using a floating shock fitting technique. Rapid convergence to steady-state solutions is attained with a diagonalized approximate factorization (AF) algorithm. Results are presented for Riemann's problem, for a regular shock reflection at an inviscid wall, for supersonic flow past a cylinder, and for a transonic airfoil. All computed shocks are ideally sharp and in excellent agreement with other numerical results or 'exact' solutions. Most importantly, this has been accomplished on unusually crude meshes without any attempt to align grid lines with shock fronts or to cluster grid lines around shocks.

Hartwich, PETER-M.↗

An Algorithm for Hyperspectral Remote Sensing of Aerosols: 3. Application to the GEO-TASO Data in KORUS-AQ Field Campaign

This paper describes the third part of a series of investigations to develop algorithms for simultaneous retrieval of aerosol parameters and surface spectral reflectance from the data measured by GEOstationary Trace gas and Aerosol Sensor Optimization (GEO-TASO) instrument. Since the algorithm is designed for future hyperspectral and geostationary satellite sensors, such as Tropospheric Emissions: Monitoring of Pollution (TEMPO), it is applied to GEO-TASO data measured over the same area by different flights as part of the Korea-United Stated Air Quality Study (KORUS-AQ) field campaign in 2016. While GEO-TAOS has a spectral sampling interval of ~0.28 nm in the visible, its data is thinned through a band selection approach with consideration of atmospheric transmittance and different surface types, which yields 20 common spectral bands to be used by the algorithm. The algorithm starts with 4 common principal components (PCs) for surface spectral reflectance extracted from various spectral libraries, and the constraints of surface reflectance parameters and aerosol scattering models respectively from k- means clustering analysis of the Rayleigh-corrected GEO-TASO spectra and AERONET data. The algorithm then proceeds iteratively with an optimal estimation approach to update PCs and retrieve aerosol optical depth (AOD) until the best match between simulated and GEOS-TASO measured spectra is achieved. The spectral AODs are compared between those retrieved from GEO-TASO (x) and 7 AERONETS (y) in 440, 550, 550, 675 nm, respectively. The comparison reveals that iterative updates of surface spectral PCs (and so surface reflectance) yield significant enhancement in AOD retrievals, improving the mean linear fitting equations from y = 0.735x + 0.088 (without update of PCs) to y = 1.055x+ 0.01, and the Pearson correlation coefficient (𝑅𝑅2) from 0.54 to 0.76, respectively. Move case studies are need to further evaluate the algorithm for its application to TEMPO that carries a enhanced version of GEO-TASO instrument.

Tropospheric Emissions: Monitoring of Pollution (T↗

Infrared Response of H2 to X-Rays in Dense Clouds

The excitation by X-rays and cosmic rays of molecular hydrogen in interstellar clouds is analyzed. We carried out detailed calculations of entry efficiencies in rovibrational levels of H2 following impact with fast electrons produced by X-ray ionization of the gas. The competing effect of collisional excitation, and quenching by the ambient gas is examined in detail. Up to date values for H-H2 collisional rate coefficients are adopted, and some derivations of H2-H2 rovibrational rate coefficients from existing literature data are proposed. Several models as a function of temperature, density, and ionization rate are presented. We found that H2 infrared emission in X-ray dominated regions (XDR) is potentially observable for temperatures and ionization rates lower than certain critical values (typically T < 1000 K and zeta/n(sub H) < 10(exp -15) cc/s where zeta is the ionization rate). At higher temperatures, collisional excitation by the ambient gas dominates the population of low vibrational levels, and at higher values of zeta/n(sub H) the abundance of H2 is negligible. If such conditions are satisfied, the resulting infrared emission spectrum can be used as a diagnostic of nearby X-ray sources such as in cooling flows in galaxy clusters, quasars, Seyfert galaxies and supernova remnants. The intensity ratio of the 2-1S(1) and 1-0S(1) lines measured for the Seyfert galaxy NGC 1275 is consistent with X-ray pumping.

Tine, S.↗

How Much Global Burned Area Can Be Forecast on Seasonal Time Scales Using Sea Surface Temperatures?

Large-scale sea surface temperature (SST) patterns influence the interannual variability of burned area in many regions by means of climate controls on fuel continuity, amount, and moisture content. Some of the variability in burned area is predictable on seasonal timescales because fuel characteristics respond to the cumulative effects of climate prior to the onset of the fire season. Here we systematically evaluated the degree to which annual burned area from the Global Fire Emissions Database version 4 with small fires (GFED4s) can be predicted using SSTs from 14 different ocean regions. We found that about 48 of global burned area can be forecast with a correlation coefficient that is significant at a p < 0.01 level using a single ocean climate index (OCI) 3 or more months prior to the month of peak burning. Continental regions where burned area had a higher degree of predictability included equatorial Asia, where 92% of the burned area exceeded the correlation threshold, and Central America, where 86% of the burned area exceeded this threshold. Pacific Ocean indices describing the El Nino-Southern Oscillation were more important than indices from other ocean basins, accounting for about 1/3 of the total predictable global burned area. A model that combined two indices from different oceans considerably improved model performance, suggesting that fires in many regions respond to forcing from more than one ocean basin. Using OCI-burned area relationships and a clustering algorithm, we identified 12 hotspot regions in which fires had a consistent response to SST patterns. Annual burned area in these regions can be predicted with moderate confidence levels, suggesting operational forecasts may be possible with the aim of improving ecosystem management.

Chen, Yang↗

Property enchancement of polyimide films by way of the incorporation of lanthanide metal ions

Lanthanide metal ions were incorporated into the polyimide derived from 2,2-bis(3,4-dicarboxyphenyl) hexafluoropropane dianhydride (6FDA) and 1,3-bis(aminophenoxy) benzene (APB) in an attempt to produce molecular level metal-polymer composites. The lanthanide series of metal ions (including aluminum, scandium, and yttrium) provide discrete and stable metal ions in the 3+ oxidation state. Throughout the series there is a uniform variation in ionic size ranging from 50 pm for aluminum to a maximum of 103.4 pm for cerium and gradually decreasing again to 84.8 pm for lutetium. The high charge-to-size ratio for these ions as well as the ability to obtain large coordination numbers makes them excellent candidates for interacting with the polymer substructure. The distinct lack of solubility of simple lanthanide salts such as the acetates and halides has made it difficult to obtain metal ions distributed in the polymer framework as discrete ions or metal complexes rather than microcomposites of metal clusters. (Lanthanum nitrates are quite soluble, but the presence of the strongly oxidizing nitrate ion leads to serious degradation of the polymer upon thermal curing. This work was successful at extending the range of soluble metals salts by using chelating agents derived from the beta-diketones dipivaloylmethane, dibenzoylmethane, trifluoroacetylacetone, and hexafluoroacetylacetone. Metal acetates which are insoluble in dimethylacetamide dissolve readily in the presence of the diketones. Addition of the polyimide yields a homogeneous resin which is then cast into a clear film. Upon curing clear films were obtained with the dibenzoylmethane and trifluoroacetylacetone ligands. The dipavaloylmethane precipitates the metal during the film casting process, and hexafluoroacetylacetone gives cured films which are deformed and brittle. These clear films are being evaluated for the effect of the metal ions on the coefficient of thermal expansion, resistance to atomic oxygen, and on selective gas permeability. Much more commonly than above, polyimide films are prepared by casting the film as the poly(amic acid) precursor which is then converted to the imidized form during the thermal cure cycle. Very limited success was achieved in the past in adding lanthanide metal ions to the amide precursors because of gellation and lack of solubility. With the use of the diketone ligands cited above, the solubility and gellation problems were overcome. However, the films after curing were clear but unacceptably brittle. Attempts to overcome this cure embrittlement problem are in progress.

Thompson, David W.↗

A Global End-Member Approach to Derive aCDOM(440) from Near-Surface Optical Measurements

This study establishes an optical inversion scheme for deriving the absorption coefficient of colored (or chromophoric, depending on the literature) dissolved organic material (CDOM) at the 440 nm wavelength, which can be applied to global water masses with near-equal efficacy. The approach uses a ratio of diffuse attenuation coefficient spectral end members, i.e., a short and long wavelength pair. The global perspective is established by sampling "extremely" clear water plus a generalized extent in turbidity and optical properties that each span three decades of dynamic range. A unique data set was collected in oceanic, coastal, and inland waters (as shallow as 0.6 m) from the North Pacific Ocean, the Arctic Ocean, Hawaii, Japan, Puerto Rico, and the east and west coasts of the United States. The data were partitioned using subjective categorizations to define a validation quality subset of conservative water masses, i.e., the inflow and outflow of properties constrain the range in the gradient of a constituent, plus 15 subcategories of water masses that were not evolving conservatively. The dependence on subcategories was confirmed with an objective methodology based on cluster analysis techniques. The latter defined five distinct classes with validation quality data present in all classes, but which also decreased in percent composition as a function of increasing class number and optical complexity. Four different algorithms based on different validation quality end members were validated with accuracies of 1.–6.2 %, wherein pairs with the largest spectral span were most accurate. Although algorithm accuracy decreased with the inclusion of more subcategories containing non-conservative water masses, changes to the algorithm fit were small when a preponderance of subcategories were included. The high accuracy for all end-member algorithms was the result of data acquisition and data processing improvements, e.g., increased vertical sampling resolution to less than 1mm and a boundary constraint to mitigate wave focusing effects, respectively. An independent evaluation with a historical database confirmed the consistency of the algorithmic approach and its application to quality assurance, e.g., to flag data outside expected ranges, identify suspect spectra, and objectively determine the in-water extrapolation interval by converging agreement for all applicable end-member algorithms. The legacy data exhibit degraded performance (as 44 % uncertainty) due to a lack of high-quality near-surface observations, especially for clear waters wherein wave-focusing effects are problematic. The novel optical approach allows the in situ estimation of an in-water constituent in keeping with the accuracy obtained in the laboratory.

Stanford B Hooker↗