Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “importance splitting”

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 73 records · Page 4

Analysis of Complex Valve and Feed Systems

A numerical framework for analysis of complex valve systems supports testing of propulsive systems by simulating key valve and control system components in the test loop. In particular, it is designed to enhance the analysis capability in terms of identifying system transients and quantifying the valve response to these transients. This system has analysis capability for simulating valve motion in complex systems operating in diverse flow regimes ranging from compressible gases to cryogenic liquids. A key feature is the hybrid, unstructured framework with sub-models for grid movement and phase change including cryogenic cavitations. The multi-element unstructured framework offers improved predictions of valve performance characteristics under steady conditions for structurally complex valves such as pressure regulator valve. Unsteady simulations of valve motion using this computational approach have been carried out for various valves in operation at Stennis Space Center such as the split-body valve and the 10-in. (approx.25.4-cm) LOX (liquid oxygen) valve and the 4-in. (approx.10 cm) Y-pattern valve (liquid nitrogen). Such simulations make use of variable grid topologies, thereby permitting solution accuracy and resolving important flow physics in the seat region of the moving valve. An advantage to this software includes possible reduction in testing costs incurred due to disruptions relating to unexpected flow transients or functioning of valve/flow control systems. Prediction of the flow anomalies leading to system vibrations, flow resonance, and valve stall can help in valve scheduling and significantly reduce the need for activation tests. This framework has been evaluated for its ability to predict performance metrics like flow coefficient for cavitating venturis and valve coefficient curves, and could be a valuable tool in predicting and understanding anomalous behavior of system components at rocket propulsion testing and design sites.

Ahuja, Vineet↗

Microwave spectra of van der Waals complexes of importance in planetary atmospheres

The Fourier-transform Fabry-Perot pulsed-molecular-beam microwave spectrometer at NIST was used to study the microwave spectra of a number of molecular dimers and trimers that may be present in planetary atmospheres. The weak van der Waals bonds associated with these species usually give rise to rotational-tunneling splittings in the microwave spectra. The microwave spectrum of the water dimer species was used to illustrate the complications that can arise in the study of the rotational spectra of these loosely bound species. In addition to the water dimer species, the microwave spectra of the following hydrogen-bonded and van der Waals complexes were studied: (CO2)2-H2O, CO2-(H2O)2, CO2-H2S, N2-H2O, CO-H2O, SO2-H2O, and O3-H2O.

Suenram, R. D.↗

Ballistic-Electron-Emission-Microscopy of Strained Si(sub 1-x)Ge(sub x) Layers

Ballistic-electron-emission microscopy (BEEM) has been used to investigate the effects of strain on Si(sub 1-x)Ge(sub x) alloys. Lifting of the degeneracy of the conduction-band minimum of Si(sub 1-x)Ge(sub x), due to lattice deformation has been directly measured by application of BEEM spectroscopy to Ag/Si structures. Experimental values for this conduction-band splitting agree well with calculations. In addition, an unexpected heterogeneity in the strain of the Si(sub 1-x)Ge(sub x) layer is introduced by deposition of Au. This effect, not observed with Ag, is attributed to species interdiffusion and has important implications for metal-semiconductor devices based oil pseudomorphic Si(sub 1-x)Ge(sub x)/Si material systems.

Bell, L. D.↗

Global Variability in Sonic Boom Exposure due to Macroscopic Effects

Supersonic flight over land has been prohibited since 1973 due to the loudness of sonic booms. NASA is building the X-59 aircraft as part of its Quesst mission to demonstrate low-loudness shaped sonic booms, or “sonic thumps.” The Quesst mission will gather human perception data via a series of community noise surveys across the USA. The noise dose and perceptual response data will be provided to the International Civil Aviation Organization (ICAO) and the Federal Aviation Administration for use in determining potential future supersonic aircraft noise certification standards, effectively changing the prohibition from a speed limit to a noise limit. These noise regulations must be globally effective, as long travel distances see the largest benefit to supersonic flight. The state of the atmosphere through which a sonic boom travels affects the size of the region exposed to sound, the “carpet width” (CW), as well as the loudness. The focus of this dissertation is to understand and quantify the expected loudness and CW of sonic booms due to the macroscopic atmospheric effects around the world. A pair of large-scale propagation simulation studies were conducted using the NASA PCBoom code to compare predicted sonic boom loudness and CW statistics first across the USA and then across the world. For the USA study, near-field data of the X-59 in steady cruise was propagated at 4 cardinal headings at 138 locations through 5 years of Climate Forecast System Version 2 (CFSv2) atmospheric profiles. Results of a bootstrap forest predictor screening model indicated the importance of climate zone, latitude, ground elevation, season, and heading. It also noted the unimportance of time of day for predicting loudness and CW. The data is visualized in aggregate, and then broken out geographically, by season and heading, and by climate zone. Multiple linear regression models were fit to the data from the 138 locations so that estimates of the loudness and CW can be produced anywhere in the US. The results can aid in planning when and where to fly the X-59. For the global study, near-field data from three aircraft, the X-59 in a quiet and loud configuration, B-58, and Concorde, were propagated at four cardinal headings through data from three atmospheric models, the CFSv2, the Global Forecast System (GFS), and the ECMWF Reanalysis Version 5 (ERA5), at 100 global locations over 1 year. Results of a bootstrap forest predictor screening model indicated the importance of climate zone, ground elevation, season, and heading. Similar to the US study, the model indicated time of day was not an important predictor. The model also indicated that choice of weather model was not important, so the atmospheric model data are effectively interchangeable. The ERA5 model was chosen for use in an extension of the study to include 18 additional locations to ensure sampling of every climate zone. Loudness and CW results are shown in aggregate, and split geographically and by heading, season, and climate. Multiple linear regression models were fit to the data from the 118 locations so that estimates of loudness and CW can be produced around the world. N-waves and shaped booms did not have the same global variability. Koppen-Geiger climate zones were used as the climate zone definition for the global study. These are available as present-day and future climate projections. Making use of the multiple linear regression models, the future climate zones were input to estimate the effect of the changing climate on sonic boom loudness and CW. Results indicate that a changing climate would have little impact on the effectiveness of noise regulations.

X-59↗

Average spatial distributions of energetic O(+), O(2+), O(6+), and C(6+) ions in the magnetosphere observed by AMPTE CCE

Measurements with the charge energy mass spectrometer of the AMPTE CCE spacecraft are used to determine the relative fluxes of these ions near the equatorial plane as a function of the drift shell parameter L, the magnetic activity index Kp, and the local time LT. The O(+) and O(2+) ions have radial profiles with maxima at about L = 5 and diurnal variations with a broad maximum on the dayside. The O(6+) and C(6+) ion fluxes increase with L between L = 5 and L = 7 and level off farther out. The diurnal variations of the relative O(6+) and C(6+) fluxes exhibit a pronounced minimum on the dayside. The observations can be interpreted in terms of ion convection from the tail onto quasi-trapped drift orbits and further radial diffusion onto closed drift shells. The maximum of the O(+) and O(2+) fluxes as well as the minimum of the O(6+) and C(6+) fluxes on the dayside can be explained by drift shell splitting. The inward transport is associated with ion losses, mainly by charge exchange. O(6+) and C(6+) ions are removed already at larger distances than O(+) and O(2+). As a result, the ionosphere contributes most of the oxygen ions in the magnetosphere. Important contributions from the solar wind are encountered in the outer magnetosphere.

Kremser, G.↗

On the Recovery Stress of a Ni50.3Ti29.7Hf20 High Temperature Shape Memory Alloy

Recovery stress in shape memory alloys (SMAs), also known as blocking stress, is an important property generally obtained during heating under a dimensional constraint as the material undergoes the martensitic phase transformation. This property has been instinctively utilized in most SMA shape-setting procedures, and has been used in numerous applications such as fastening and joining, rock splitting, safety release mechanisms, reinforced composites, medical devices, and many other applications. The stress generation is also relevant to actuator applications where jamming loads (e.g., in case the actuator gets stuck and is impeded from moving) need to be determined for proper hardware sizing. Recovery stresses in many SMA systems have been shown to reach stresses in the order of 800 MPa, achieved via thermo-mechanical training such as pre-straining, heat treatments or other factors. With the advent of high strength, high temperature SMAs, recovery stress data has been rarely probed, and there is no information pertinent to the magnitudes of these stresses. Thus, the purpose of this work is to investigate the recovery stress capability of a precipitation strengthened, Ni50.3Ti29.7Hf20 (at.) high temperature SMA in uniaxial tension and compression. This material has been shown to exhibit outstanding strength and stability during constant-stress, thermal cycling, but no data exists on constant-strain thermal cycling. Several training routines were implemented as part of this work including isothermal pre-straining, isobaric thermal cycling, and isothermal cyclic training routines. Regardless of the training method used, the recovery stress was characterized using constant-strain (strain-controlled condition) thermal cycling between the upper and lower cycle temperatures. Preliminary results indicate recovery stresses in excess of 1.5 GPa were obtained after a specific training routine. This stress magnitude is significantly higher than conventional NiTi stress generation capability.

Recovery Stress↗

The generation of hydrogen by the thermal decomposition of water

Development of an approach to the evaluation of the thermal efficiency of the process of water splitting to produce hydrogen. A way of viewing thermochemical processes - both overall and step-by-step is suggested, and some recent work on a process evaluation technique is described which provides internal checks on the thermodynamic data and calculates, in addition to the efficiency, many important process parameters.

Funk, J. E.↗

Damage tolerance of pressurized graphite/epoxy tape cylinders under uniaxial and biaxial loading

The damage tolerance behavior of internally pressurized, axially slit, graphite/epoxy tape cylinders was investigated. Specifically, the effects of axial stress, structural anisotropy, and subcritical damage were considered. In addition, the limitations of a methodology which uses coupon fracture data to predict cylinder failure were explored. This predictive methodology was previously shown to be valid for quasi-isotropic fabric and tape cylinders but invalid for structurally anisotropic (+/-45/90)(sub s) and (+/-45/0)(sub s) cylinders. The effects of axial stress and structural anisotropy were assessed by testing tape cylinders with (90/0/+/-45)(sub s), (+/-45/90)(sub s), and (+/-45/0)(sub s) layups in a uniaxial test apparatus, specially designed and built for this work, and comparing the results to previous tests conducted in biaxial loading. Structural anisotropy effects were also investigated by testing cylinders with the quasi-isotropic (0/+/-45/90)(sub s) layup which is a stacking sequence variation of the previously tested (90/0/+/-45)(sub s) layup with higher D(sub 16) and D(sub 26) terms but comparable D(sub 16) and D(sub 26) to D(sub 11) ratios. All cylinders tested and used for comparison are made from AS4/3501-6 graphite/epoxy tape and have a diameter of 305 mm. Cylinder slit lengths range from 12.7 to 50.8 mm. Failure pressures are lower for the uniaxially loaded cylinders in all cases. The smallest percent failure pressure decreases are observed for the (+/-45/90)(sub s) cylinders, while the greatest such decreases are observed for the (+/-45/0)(sub s) cylinders. The relative effects of the axial stress on the cylinder failure pressures do not correlate with the degree of structural coupling. The predictive methodology is not applicable for uniaxially loaded (+/-45/90)(sub s) and (+/-45/0)(sub s) cylinders, may be applicable for uniaxially loaded (90/0/+/-45)(sub s) cylinders, and is applicable for the biaxially loaded (90/0/+/-45)(sub s) and (0/+/-45/90)(sub s) cylinders. This indicates that the ratios of D(sub 16) and D(sub 26) to D(sub 11), as opposed to the absolute magnitudes of D(sub 16) and D(sub 26), may be important in the failure of these cylinders and in the applicability of the methodology. Discontinuities observed in the slit tip hoop strains for all the cylinders tested indicate that subcritical damage can play an important role in the failure of tape cylinders. This role varies with layup and loading condition and is likely coupled to the effects of structural anisotropy. Biaxial failure pressures may exceed the uniaxial values because the axial stress contributes to the formation of 0 deg ply splitting (accompanied by delamination) or similar stress-mitigating subcritical damage. The failure behavior of similar cylinders can also vary as a result of differences in the role of subcritical damage as observed for the case of a biaxially loaded (90/0/+/-45)(sub s) cylinder with a 12.7 mm slit. For this case, the methodology is valid when the initial coupon and cylinder fracture modes agree. However, the methodology underpredicts the failure pressure of the cylinder when a circumferential fracture path, suggestive of a 0 deg ply split, occurs at one slit tip. Thus, the failure behavior of some tape cylinders may be highly sensitive to the initial subcritical damage mechanism. Finite element analyses are recommended to determine how structural anisotropy and axial stress modify the slit tip stress states in cylinders from those found in flat plates since similarity of these stress states is a fundamental assumption of the current predictive methodology.

Priest, Stacy Marie↗

Evidence for Liquid Water on Comets

We have reexamined the arguments for the existence of liquid water on comets, and believe that recent cometary flybys along with pre-Giotto data support its presence on short-period comets. Liquid water would affect cometary dynamics, leaving distinct signatures in precession, orbital dynamics, and potential splitting of comets. Liquid water geysers would affect cometary atmosphere, dust evolution, and non-gravitational forces that perturb the orbit. Liquid water would affect the composition of both the interior and exterior of the comet, producing geologic effects consistent with recent flyby photographs. And most importantly, liquid water suppork the growth of lifeforms, which would make a comet a biofriendly incubator for interplanetary transport. The major objection against liquid water is the necessity of a pressure vessel to prevent sublimation into space. We discuss how such a pressure vessel could naturally evolve as a pristine comet makes its first journey inside the orbit of Mars, and suggest that this type of vessel was observed by Giotto, Deep Space I, and Stardust.

Sheldon, Robert↗

The Abundance of Interstellar Fluorine

The primary objective of this program was to obtain FUSE observations of the interstellar absorption lines of F I at 951 and 954 Angstroms to derive the abundance of fluorine toward the star HD 164816. The nucleosynthetic source(s) of fluorine are still a matter of debate - the present day abundance of fluorine can potentially constrain models for pulsationally driven dredge-up in asymptotic giant branch stars. An accurate measure for the depletion behavior of fluorine will determine whether it may be detectable in QSO absorption line systems - an unambiguous detection of fluorine at suitably high redshifts would provide the best evidence to date for the neutrino process in massive stars. Furthermore, due to its extreme reactivity, measurement of the gas-phase interstellar fluorine abundance is important for models of grain chemistry. Despite the importance of measuring the interstellar fluorine abundance, at the time of our proposal only one previous detection has been made due to the low relative abundance of fluorine, the lack of lines outside the far-UV, and the blending of the available F I transitions with lines of Hz. The star HD 164816 is associated with the Lagoon nebula (M8), and at a distance of approximately 1.5 kpc probes both distant and local gas. Beginning April 8th, 2004 FUSE FP-Split observations of the star HD 164816 were obtained for this program. This data became available in the FUSE data archive May 21, 2004, and these observations were then downloaded and we began our analysis. Our analysis procedure has involved (1) fitting stellar models to the FUSE spectra, (2) using the multiple lines of Hz and N I at other wavelengths in the FUSE bandpass to derive column densities for the lines of H2 and N I which are blended with the F I features at 951 and 954 angstroms (3) the measurement of the column densities of F I and the species O I and C1 I which are important species for the dis-entangling of dust and nucleosynthetic effects. As discussed in our poster paper presented at the Winter 2005 American Astronomical Society meeting (Lauroesch et al. 2004, AAS, 205, 57.06), our analysis of the data taken for this star suggests that the observed fluorine gas-phase abundance (and hence inferred depletion onto dust grains) are consistent with previous measures toward delta Sco (Copernicus) and two recent measures toward stars in Cep OB2 obtained by others with FUSE. We are currently revising the draft a paper for submission to the Astrophysical Journal which presents the results of these observations and in addition our recently obtained FUSE observations of another star (HD 103779, FUSE GO program D073) which shows an apparent slight enhancement in the gas-phase fluorine abundance.

Lauroesch, James T.↗

Neural crest contributions to the lamprey head

The neural crest is a vertebrate-specific cell population that contributes to the facial skeleton and other derivatives. We have performed focal DiI injection into the cranial neural tube of the developing lamprey in order to follow the migratory pathways of discrete groups of cells from origin to destination and to compare neural crest migratory pathways in a basal vertebrate to those of gnathostomes. The results show that the general pathways of cranial neural crest migration are conserved throughout the vertebrates, with cells migrating in streams analogous to the mandibular and hyoid streams. Caudal branchial neural crest cells migrate ventrally as a sheet of cells from the hindbrain and super-pharyngeal region of the neural tube and form a cylinder surrounding a core of mesoderm in each pharyngeal arch, similar to that seen in zebrafish and axolotl. In addition to these similarities, we also uncovered important differences. Migration into the presumptive caudal branchial arches of the lamprey involves both rostral and caudal movements of neural crest cells that have not been described in gnathostomes, suggesting that barriers that constrain rostrocaudal movement of cranial neural crest cells may have arisen after the agnathan/gnathostome split. Accordingly, neural crest cells from a single axial level contributed to multiple arches and there was extensive mixing between populations. There was no apparent filling of neural crest derivatives in a ventral-to-dorsal order, as has been observed in higher vertebrates, nor did we find evidence of a neural crest contribution to cranial sensory ganglia. These results suggest that migratory constraints and additional neural crest derivatives arose later in gnathostome evolution.

Non-NASA Center↗

The application of experimental data to blade wake interaction noise prediction

Blade wake interaction noise (BWI) has been defined as the broadband noise generated by the ingestion of turbulent trailing tip vortices by helicopter rotors. This has been shown to be the dominant contributor to the subjectively important part of the acoustic spectrum for the approach stage of a helicopter flyover. A prediction method for BWI noise based on the calculated trailing vortex trajectories has been developed and estimates of the vortex turbulence have been made. These measurements were made on a trailing vortex from a split wing arrangement and did not give the spectrum of the velocity fluctuations. A recent experiment carried out to measure the turbulence associated with a trailing vortex and the application of the results to BWI noise prediction is described.

Glegg, Stewart A. L.↗

New Experiment to Measure the Electron Electric Dipole Moment

An electron can possess an electric dipole moment (edm) only if time reversal symmetry (T) is violated. No edm of any particle has yet been discovered. CP-violation, equivalent to T-violation by the CPT theorem, does occur in Kaon decays and can be accounted for by the standard model. However, this mechanism leads to an electron edm d(sub e) of the order of 10(exp -38) e cm, whereas the current experimental bound on d(sub e) is about 10(exp -27) e cm. However, well-motivated extensions of the standard model such as supersymmetric theories do predict that de could be as large as the current bound. In addition, CP violation in the early universe is required to explain the preponderance of matter over anti-matter, but the exact mechanism of this CP violation is unclear. For these reasons, we are undertaking a new experimental program to determine de to an improved accuracy of 10(exp -29) e cm. Our experiment will use laser-cooled, trapped Cesium atoms to measure the atomic edm d(sub Cs) that occurs if d(sub e) is not zero. In order to do this, we will measure the energy splitting between the atoms spin states in parallel electric and magnetic fields. The signature of an edm would be a linear dependence of the splitting on the electric field E due to the interaction - d(sub Cs) dot E. Our measurement will be much more sensitive than previous measurements because atoms can be stored in the trap for tens of seconds, allowing for much narrower Zeeman resonance linewidths. Also, our method eliminates the most important systematic errors, proportional to atomic velocity, which have limited previous experiments. In this presentation, we will describe the design of our new apparatus, which is presently under construction. An important feature of our experimental apparatus is that magnetic field noise will be suppressed to a very low value of the order of 1 fT/(Hz)1/2. This requires careful attention to the Johnson noise currents in the chamber, which have not been important in previous experiments. In addition we will present estimates of the limits of the various errors that we expect for our experiment.

Kittle, Melanie↗

On Variable Geometric Factor Systems for Top-Hat Electrostatic Space Plasma Analyzers

Even in the relatively small region of space that is the Earth's magnetosphere, ion and electron fluxes can vary by several orders of magnitude. Top-hat electrostatic analyzers currently do not possess the dynamic range required to sample plasma under all conditions. The purpose of this study was to compare, through computer simulation, three new electrostatic methods that would allow the sensitivity of a sensor to be varied through control of its geometric factor (GF) (much like an aperture on a camera). The methods studied were inner filter plates, split hemispherical analyzer (SHA) and top-cap electrode. This is the first discussion of the filter plate concept and also the first study where all three systems are studied within a common analyzer design, so that their relative merits could be fairly compared. Filter plates were found to have the important advantage that they facilitate the reduction in instrument sensitivity whilst keeping all other instrument parameters constant. However, it was discovered that filter plates have numerous disadvantages that make such a system impracticable for a top-hat electrostatic analyzer. It was found that both the top-cap electrode and SHA are promising variable geometric factor system (VGFS) concepts for implementation into a top-hat electrostatic analyzer, each with distinct advantages over the other.

Kataria, Dhiren O.↗

Security Data Warehouse Application

The Security Data Warehouse (SDW) is used to aggregate and correlate all JSC IT security data. This includes IT asset inventory such as operating systems and patch levels, users, user logins, remote access dial-in and VPN, and vulnerability tracking and reporting. The correlation of this data allows for an integrated understanding of current security issues and systems by providing this data in a format that associates it to an individual host. The cornerstone of the SDW is its unique host-mapping algorithm that has undergone extensive field tests, and provides a high degree of accuracy. The algorithm comprises two parts. The first part employs fuzzy logic to derive a best-guess host assignment using incomplete sensor data. The second part is logic to identify and correct errors in the database, based on subsequent, more complete data. Host records are automatically split or merged, as appropriate. The process had to be refined and thoroughly tested before the SDW deployment was feasible. Complexity was increased by adding the dimension of time. The SDW correlates all data with its relationship to time. This lends support to forensic investigations, audits, and overall situational awareness. Another important feature of the SDW architecture is that all of the underlying complexities of the data model and host-mapping algorithm are encapsulated in an easy-to-use and understandable Perl language Application Programming Interface (API). This allows the SDW to be quickly augmented with additional sensors using minimal coding and testing. It also supports rapid generation of ad hoc reports and integration with other information systems.

Vernon, Lynn R.↗

Classifying Agnostic Biosignatures using Raman, VNIR, and Elemental Data

How can we use our current wealth of terrestrial data, encompassing biogenic and abiogenic systems, to determine the distinguishing properties of life? SCOBI (Statistical Classification of Biosignature Information) uses machine learning techniques to algorithmically identify combinations of measurements that are “indicative of life”. A set of ~1000 observations, comprising elemental abundance, isotopic fractionation, VNIR reflectance, and (in progress) Raman spectra, have been assembled from existing literature and databases. The observations cover systems classified as “indicative alive” (e.g., cells, vegetation), “indicative non-alive” (e.g., fossils, teeth), “mixed indicative” (e.g., soil, pond water), or “non-indicative” (e.g., rocks, meteorites). VNIR data was preprocessed by linear interpolation from 400-2100 nm and smoothed with a Savitzky-Golay filter. To limit the amount of Earth-biochemistry-specific (non-agnostic) information included, the first five spectral features extracted were number of peaks, number of troughs, mean reflectance, mean peak width, and broadest peak width. To help further emphasize agnostic biosignatures, Earth-specific features such as chlorophylls have been manually flagged so that feature importance with and without them can be compared. Classifiers including k-nearest neighbors (KNN), Gaussian Naïve Bayes (GNB), logistic regression (LR), random forest (RF), and support vector machine (SVM) were implemented, as was a combination voting classifier. Performance metrics included false positive rates, false negative rates, and AUC with 50-50 test/train splits (Monte Carlo simulations). Key takeaways from this stage, prior to the inclusion of Raman spectra, are (1) the overall success rate of 0.933 AUC was most heavily influenced by the elemental abundance data; and (2) VNIR reflectance had the lowest classification performance with 0.52 AUC (58% of objects correctly classified). The next steps are to complete integration of Raman spectral data and to improve the approach to pre-processing and feature extraction for both types of spectral data, such as automated baseline removal, whole spectrum matching, and dimensionality reduction.

Biosignatures↗

On the Shape of Liquid Metal Droplets in Electromagnetic Levitation Experiments

We present calculations and measurements on the shape of liquid metal droplets in electromagnetic levitation experiments. A normal stress balance model was developed to predict the shapes of liquid metal droplets that will be obtained in a microgravity experiment to measure the viscosity and surface tension of undercooled metals. This model was tested by calculating the droplet shapes in containerless experiments conducted to determine the surface tension of liquid metals. Inconsistencies associated with the results of a previous paper are elucidated. The computational results of the mathematical model are compared with the results of ground-based experiments for two different metals. The importance of the ratio of electromagnetic skin depth-to-droplet radius to the accuracy of the mathematical model is discussed. A planned alternate approach to modeling the shape by consideration of the entire droplet rather than only the surface is presented. As an example of an application. the influence of the shape on the splitting of the surface oscillation modes of levitated liquid metal droplets is discussed.

Schwartz, E.↗

A Fast Implementation of the ISOCLUS Algorithm

Unsupervised clustering is a fundamental tool in numerous image processing and remote sensing applications. For example, unsupervised clustering is often used to obtain vegetation maps of an area of interest. This approach is useful when reliable training data are either scarce or expensive, and when relatively little a priori information about the data is available. Unsupervised clustering methods play a significant role in the pursuit of unsupervised classification. One of the most popular and widely used clustering schemes for remote sensing applications is the ISOCLUS algorithm, which is based on the ISODATA method. The algorithm is given a set of n data points (or samples) in d-dimensional space, an integer k indicating the initial number of clusters, and a number of additional parameters. The general goal is to compute a set of cluster centers in d-space. Although there is no specific optimization criterion, the algorithm is similar in spirit to the well known k-means clustering method in which the objective is to minimize the average squared distance of each point to its nearest center, called the average distortion. One significant feature of ISOCLUS over k-means is that clusters may be merged or split, and so the final number of clusters may be different from the number k supplied as part of the input. This algorithm will be described in later in this paper. The ISOCLUS algorithm can run very slowly, particularly on large data sets. Given its wide use in remote sensing, its efficient computation is an important goal. We have developed a fast implementation of the ISOCLUS algorithm. Our improvement is based on a recent acceleration to the k-means algorithm, the filtering algorithm, by Kanungo et al.. They showed that, by storing the data in a kd-tree, it was possible to significantly reduce the running time of k-means. We have adapted this method for the ISOCLUS algorithm. For technical reasons, which are explained later, it is necessary to make a minor modification to the ISOCLUS specification. We provide empirical evidence, on both synthetic and Landsat image data sets, that our algorithm's performance is essentially the same as that of ISOCLUS, but with significantly lower running times. We show that our algorithm runs from 3 to 30 times faster than a straightforward implementation of ISOCLUS. Our adaptation of the filtering algorithm involves the efficient computation of a number of cluster statistics that are needed for ISOCLUS, but not for k-means.

Memarsadeghi, Nargess↗