Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothing methods”

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,135 records · Page 63

Resolving topographic detail on Venus by modeling complex Magellan altimetry echoes

Magellan's altimeter is providing some of the finest resolution topography of Venus achieved to date. Nevertheless, efforts continue to improve the topographic resolution whenever possible. One effort to this end is stereoscopic imaging, which provides topography at scales similar to that of the synthetic aperture radar (SAR). However, this technique requires two SAR images of the same site to be obtained and limits the utility of this method. In this paper, we present another method to resolve topographic features at scales smaller than that of an altimeter footprint, which is more globally applicable than the stereoscopic approach. Each pulse which is transmitted by Magellan's altimeter scatters from the planet and echoes to the receiver, delayed based on the distance between the spacecraft and each surface element. As resolved in time, each element of an altimetry echo represents the sum of all points on the surface which are equidistant from the spacecraft. Thus, individual returns, as a function of time, create an echo profile which may be used to derive properties of the surface, such as the scattering law or, in this case, the topography within the footprint. The Magellan project has derived some of this information by fitting model templates to radar echo profiles. The templates are calculated based on Hagfor's Law, which assumes a smooth, gently undulating surface. In most regions these templates provide a reasonable fit to the observed echo profile; however, in some cases the surface departs from these simple assumptions and more complex profiles are observed. Specifically, we note that sub-footprint topographic relief apparently has a strong effect on the shape of the echo profile. To demonstrate the effects of sub-resolution relief on echo profiles, we have calculated the echo shapes from a wide range of simple topographic models. At this point, our topographic models have emphasized surfaces where only two dominant elevations are contained within a footprint, such as graben, ridges, crater rims, and central features in impact craters.

Lovell, Amy J.↗

Effects of Microgravity or Simulated Launch on Testicular Function in Rats

Testes from flight rats on COSMOS 2044 and simulated-launch, vivarium, or caudal-elevation control rats (5/group) were analyzed by subjective and quantitative methods. On the basis of observations of fixed tissue, it was evident that some rats had testicular abnormalities unassociated with treatment and probably existing when they were assigned randomly to the four treatment groups. Considering rats without preexisting abnormalities, diameter of seminiferous tubules and numbers of germ cells per tubule cross section were lower (P less than 0.05) in flight than in simulated-launch or vivarium rats. However, ratios of germ cells to each other or to Sertoli cells and number of homogenization-resistant spermatids did not differ from values for simulated-launch or vivarium controls. Expression of testis-specific gene products was not greatly altered by flight. Furthermore, there was no evidence for production of stress-inducible transcripts of the hsp7O or hsp9O genes. Concentration of receptors for rat luteinizing hormone in testicular tissue and surface density of smooth endoplasmic reticulum in Leydig cells were similar in flight and simulated-launch rats. However, concentrations of testosterone in testicular tissue or peripheral blood plasma were reduced (P less than 0.05) in flight rats to less than 20% of values for simulated-launch or vivarium controls. Thus spermatogenesis was essentially normal in flight rats, but production of testosterone was severely depressed. Exposure to microgravity for more than 2 wk might result in additional changes. Sequelae of reduced androgen production associated with microgravity on turnover of muscle and bone should be considered.

Amann, R. P.↗

Guidance of Nonlinear Systems

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of way points through which the aircraft trajectory must pass. The way points typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory which satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multi-dimensional, multi-axis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of possible operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordination, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions must be smooth. The guidance algorithm is based on the inversion of the pure feedback approximations, which is followed by iterative corrections for the effects of zero dynamics. The paper describes the structure and modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

Guidance of Nonlinear Systems

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of way points through which the aircraft trajectory must pass. The way points typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory which satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multi-dimensional, multi-axis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of possible operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordinatization, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions must be smooth. The guidance algorithm is based on the inversion of the pure feedback approximations, which is followed by iterative corrections for the effects of zero dynamics. The paper describes the structure and modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

An Examination of Selected Geomagnetic Indices in Relation to the Sunspot Cycle

Previous studies have shown geomagnetic indices to be useful for providing early estimates for the size of the following sunspot cycle several years in advance. Examined this study are various precursor methods for predicting the minimum and maximum amplitude of the following sunspot cycle, these precursors based on the aa and Ap geomagnetic indices and the number of disturbed days (NDD), days when the daily Ap index equaled or exceeded 25. Also examined is the yearly peak of the daily Ap index (Apmax), the number of days when Ap greater than or equal to 100, cyclic averages of sunspot number R, aa, Ap, NDD, and the number of sudden storm commencements (NSSC), as well the cyclic sums of NDD and NSSC. The analysis yields 90-percent prediction intervals for both the minimum and maximum amplitudes for cycle 24, the next sunspot cycle. In terms of yearly averages, the best regressions give Rmin = 9.8+/-2.9 and Rmax = 153.8+/-24.7, equivalent to Rm = 8.8+/-2.8 and RM = 159+/-5.5, based on the 12-mo moving average (or smoothed monthly mean sunspot number). Hence, cycle 24 is expected to be above average in size, similar to cycles 21 and 22, producing more than 300 sudden storm commencements and more than 560 disturbed days, of which about 25 will be Ap greater than or equal to 100. On the basis of annual averages, the sunspot minimum year for cycle 24 will be either 2006 or 2007.

Wilson, Robert M.↗

Investigation by Schlieren Technique of Methods of Fixing Fully Turbulent Flow on Models at Supersonic Speeds

An investigation has been made of the effectiveness of two-dimensional wire and three-dimensional grandular-type roughnesses for fixing transition on a 27 deg cone at a free-stream Mach number of 2.20. One height of each type of roughness was investigated over a free-stream Reynolds number per foot range of about 1.0 x 10 (exp 6) to 7.0 x 10 (exp 6). The tests were made at an angle of attack of 0 deg and with essentially zero heat transfer. The results indicate that, within the Reynolds number per foot range investigated, both two-dimensional wire and three-dimensional granular types of roughness caused large forward movements of transition at constant Reynolds number per foot; however, the forward movement with increase in Reynolds number per foot was very slow for the two-dimensional wire roughness and transition never reached the roughness location as it did with the three-dimensional granular roughness. An extrapolated value of the roughness Reynolds number (based on roughness height and flow conditions at the outer edge of the roughness) for the first appearance of forward movement from the smooth-cone case for the three-dimensional granular roughness was in good agreement with the value of about 600 for the first appearance of turbulent bursts obtained on a similar cone by the hot-wire technique. The range of values of roughness Reynolds number for which a fully turbulent boundary layer existed in the vicinity of the roughness element was about 4,500 to 6,500.

Jackson, Mary W.↗

A method for correcting aspect solution errors in ROSAT HRI observations of compact sources

X-ray point sources observed with the ROSAT High Resolution Imager (HRI) often appear elongated over scales of approximately 5 sec-10 sec from the image core. This elongation has been attributed to errors in the attitude correction as the satellite is wobbled during the observations, and affects sources with both soft and hard X-ray spectra. In this paper, I report the results of an attempt to rid a high signal-to-noise observation of the soft X-ray point source HZ 43 of its characteristic elongation. I divided the observation into 181 separate images, each containing photons from only a small region on the detector through which the source passed during the satellite's wobble. By measuring the positions of the individual image centroids, I found clear evidence for systematic offsets from a common mean by up to approximately +/- 3 sec in both right ascension and declination as a function of phase in the satellite wobble. Shifting the subimages to a common center and then restacking them into a single image measurably improved the symmetry of the point-spread function. HRI observations are wobbled primarily to smooth out variations in the pixel-to-pixel sensitivity of the detector and also to extend the lifetime of the microchannel plates in the detector since these decay at a given location as a function of the number of photons detected. However, the elongations introduced by the aspect errors inhibit the identification of possible extended X-ray emission associated with sources such as pulsars and active galactic nuclei. In light of these results, I suggest that until the aspect errors are understood, observations of compact sources, where this effect may be important, should not be wobbled.

Morse, Jon A.↗

The X-ray Halo of GX5-1

Using Chandra observations we have measured the energy-resolved dust-scattered X-ray halo around the low-mass X-ray binary GX5-1, detecting for the first time multiply scattered X-rays from interstellar dust. % e compared the observed X-ray halo at various energies to predictions from a range of dust models. These fits used both smoothly-distributed dust as well as dust in clumped clouds, with CO and 21 cm observations helping to determine the position of the clouds along the line of sight. We found that the BARE-GR-B model of Zubko, Dwek & Arendt (2004) generally led to the best results, although inadequacies in both the overall model and the data limit our conclusions. We did find that the composite dust models of Zubko, Dwek & Arendt (2004), especially the "no carbon" models, gave uniformly poor results. Although models using cloud positions and densities derived naively from CO and 21 cm data gave generally poor results, plausible adjustments to the distance of the largest cloud and the mass of a cloud in the expanding 3 kpc Arm lead to significantly improved fits. We suggest that combining X-ray halo, CO, and 21 cm observations will be a fruitful method to improve our understanding of both the gas and dust phases of the interstellar medium.

Smith, Randall K.↗

The XXL Survey I. Scientific Motivations - Xmm-Newton Observing Plan - Follow-up Observations and Simulation Programme

The quest for the cosmological parameters that describe our universe continues to motivate the scientific community to undertake very large survey initiatives across the electromagnetic spectrum. Over the past two decades, the Chandra and XMM-Newton observatories have supported numerous studies of X-ray-selected clusters of galaxies, active galactic nuclei (AGNs), and the X-ray background. The present paper is the first in a series reporting results of the XXL-XMM survey; it comes at a time when the Planck mission results are being finalized. Aims. We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of approx. 5 x 10(exp 15) erg/s/sq cm in the [0.5-2] keV band (completeness limit). The surveys main goals are to provide constraints on the dark energy equation of state from the space-time-distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programs. Methods. We describe the 542 XMM observations along with the associated multi- and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. Results. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. Conclusions. The XXL multi- data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve asa calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps.

X-rays: general – large-scale structure of Unive↗

The Transition of High-Resolution NASA MODIS Sea Surface Temperatures into the WRF Environmental Modeling System

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed a Moderate Resolution Imaging Spectroradiometer (MODIS) sea surface temperature (SST) composite at 2-km resolution that has been implemented in version 3 of the National Weather Service (NWS) Weather Research and Forecasting (WRF) Environmental Modeling System (EMS). The WRF EMS is a complete, full physics numerical weather prediction package that incorporates dynamical cores from both the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). The installation, configuration, and execution of either the ARW or NMM models is greatly simplified by the WRF EMS to encourage its use by NWS Weather Forecast Offices (WFOs) and the university community. The WRF EMS is easy to run on most Linux workstations and clusters without the need for compilers. Version 3 of the WRF EMS contains the most recent public release of the WRF-NMM and ARW modeling system (version 3 of the ARW is described in Skamarock et al. 2008), the WRF Pre-processing System (WPS) utilities, and the WRF Post-Processing program. The system is developed and maintained by the NWS National Science Operations Officer Science and Training Resource Coordinator. To initialize the WRF EMS with high-resolution MODIS SSTs, SPoRT developed the composite product consisting of MODIS SSTs over oceans and large lakes with the NCEP Real-Time Global (RTG) filling data over land points. Filling the land points is required due to minor inconsistencies between the WRF land-sea mask and that used to generate the MODIS SST composites. This methodology ensures a continuous field that adequately initializes all appropriate arrays in WRF. MODIS composites covering the Gulf of Mexico, western Atlantic Ocean and the Caribbean are generated daily at 0400, 0700, 1600, and 1900 UTC corresponding to overpass times of the NASA Aqua and Terra polar orbiting satellites. The MODIS SST product is output in gridded binary-1 (GRIB-1) data format for a seamless incorporation into WRF via the WPS utilities. The full-resolution, 1-km MODIS product is sub-sampled to 2-km grid spacing due to limitations in handling very large dimensions in the GRIB-1 data format. The GRIB-1 files are posted online at ftp://ftp.nsstc.org/sstcomp/WRF/, which is directly accessed by the WRF EMS scripts. The MODIS SST composites are also downloaded to the EMS data server, which is accessible by the WRF EMS users and NWS WFOs. The SPoRT MODIS SST composite provides the model with superior detail of the ocean gradients around Florida and surrounding waters, whereas the operational RTG SST typically depicts a relatively smooth field and is not able to capture sharp horizontal gradients in SST. Differences of 2-3 C are common over small horizontal distances, leading to enhanced SST gradients on either side of the Gulf Stream and along the edges of the cooler shelf waters. These sharper gradients can in turn produce atmospheric responses in simulated temperature and wind fields as depicted in LaCasse et al. Differences in atmospheric verification statistics over a several month study were generally small in the vicinity of south Florida; however, the validation of SSTs at specific buoy locations revealed important improvements in the biases and RMS errors, especially in the vicinity of the cooler shelf waters off the east-central Florida coast. A current weakness in the MODIS SST product is the occurrence of occasional discontinuities caused by high latency in SST coverage due to persistent cloud cover. An enhanced method developed by Jedlovec et al. (2009, GHRSST User Symposium) reduces the occurrence of these problems by adding Advanced Microwave Scanning Radiometer -- EOS (AMSR-E) SST data to the compositing process. Enhanced SST composites are produced over the ocean regions surrounding the Continental U.S. at four times each day corresponding to Terra and Aqua equator crossing times. For a given day and overpass time, both MODInd AMSR-E data from the previous seven days form a collection used in the compositing. At each MODIS pixel, cloud-free SST values from the collection are used to form a weighted average based on their latency (number of days from the current day). In this way, recent SST data are given more weight than older data. One of the primary issues involved in incorporating the AMSR-E microwave data in the composites is the tradeoff between the decreased spatial resolution of the AMSR-E data (25 km) and the increased coverage due to its near all-weather capability. Currently, the AMSR-E is given a weight of 20% compared to MODIS data, thereby preserving the spatial structure observed in the MODIS data. Day-time (night-time) AMSR-E SST data from Aqua are used with both Terra and Aqua MODIS day-time (night-time) SST data sets.

Case, Jonathan L.↗

Goddard Latent Heating Retrieval Algorithm for TRMM and GPM

The Goddard Convective-Stratiform Heating (CSH) algorithm has been used to retrieve latent heating (LH) associated with clouds and cloud systems in support of the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) mission. The CSH algorithm requires the use of a cloud-resolving model (CRM) to simulate LH profiles to build look-up tables (LUTs). In this paper we describe the current V6 CSH and its differences/similarities versus the previous V5 CSH. Long-term CRM simulations are conducted to identify the impact of CRM resolution and the convective-stratiform separation method on LH structure/profiles. The TRMM and GPM Combined algorithm derived surface rainfall rates and their associated precipitation properties are the input to the CSH algorithm. Since there is no direct measurement of LH structure, the performance of the CSH V6 algorithm is examined by comparing its vertically-integrated heating (or equivalent surface rainfall rate) against the surface rainfall rate derived from the TRMM/GPM Combined algorithm. The results suggest the CSH three-month and zonal mean equivalent surface rain rates are in good agreement with the Combined rates over the ITCZ region; the agreement is best over ocean. CSH three-month and zonal mean equivalent surface rain rates are larger than the Combined rates over land in both the tropics and sub-tropics. CSH three-month mean equivalent surface rainfall rates also have local differences with the Combined rates that can be smoothed by area averaging to larger horizontal resolutions (from the CSH standard grid of 0.25° × 0.25° to 0.5° × 0.5° or 1.0° × 1.0°). CSH equivalent surface rainfall rates have more light rain rates but less larger rain rates compared to the GPM Combined surface rain rates.

Latent heating retrieval↗

Uncertainties for two-dimensional models of solar rotation from helioseismic eigenfrequency splitting

Observed solar p-mode frequency splittings can be used to estimate angular velocity as a function of position in the solar interior. Formal uncertainties of such estimates depend on the method of estimation (e.g., least-squares), the distribution of errors in the observations, and the parameterization imposed on the angular velocity. We obtain lower bounds on the uncertainties that do not depend on the method of estimation; the bounds depend on an assumed parameterization, but the fact that they are lower bounds for the 'true' uncertainty does not. Ninety-five percent confidence intervals for estimates of the angular velocity from 1986 Big Bear Solar Observatory (BBSO) data, based on a 3659 element tensor-product cubic-spline parameterization, are everywhere wider than 120 nHz, and exceed 60,000 nHz near the core. When compared with estimates of the solar rotation, these bounds reveal that useful inferences based on pointwise estimates of the angular velocity using 1986 BBSO splitting data are not feasible over most of the Sun's volume. The discouraging size of the uncertainties is due principally to the fact that helioseismic measurements are insensitive to changes in the angular velocity at individual points, so estimates of point values based on splittings are extremely uncertain. Functionals that measure distributed 'smooth' properties are, in general, better constrained than estimates of the rotation at a point. For example, the uncertainties in estimated differences of average rotation between adjacent blocks of about 0.001 solar volumes across the base of the convective zone are much smaller, and one of several estimated differences we compute appears significant at the 95% level.

Genovese, Christopher R.↗

Highlights from a Mach 4 Experimental Demonstration of Inlet Mode Transition for Turbine-Based Combined Cycle Hypersonic Propulsion

NASA is focused on technologies for combined cycle, air-breathing propulsion systems to enable reusable launch systems for access to space. Turbine Based Combined Cycle (TBCC) propulsion systems offer specific impulse (Isp) improvements over rocket-based propulsion systems in the subsonic takeoff and return mission segments along with improved safety. Among the most critical TBCC enabling technologies are: 1) mode transition from the low speed propulsion system to the high speed propulsion system, 2) high Mach turbine engine development and 3) innovative turbine based combined cycle integration. To address these challenges, NASA initiated an experimental mode transition task including analytical methods to assess the state-of-the-art of propulsion system performance and design codes. One effort has been the Combined-Cycle Engine Large Scale Inlet Mode Transition Experiment (CCE-LIMX) which is a fully integrated TBCC propulsion system with flowpath sizing consistent with previous NASA and DoD proposed Hypersonic experimental flight test plans. This experiment was tested in the NASA GRC 10 by 10-Foot Supersonic Wind Tunnel (SWT) Facility. The goal of this activity is to address key hypersonic combined-cycle engine issues including: (1) dual integrated inlet operability and performance issues-unstart constraints, distortion constraints, bleed requirements, and controls, (2) mode-transition sequence elements caused by switching between the turbine and the ramjet/scramjet flowpaths (imposed variable geometry requirements), and (3) turbine engine transients (and associated time scales) during transition. Testing of the initial inlet and dynamic characterization phases were completed and smooth mode transition was demonstrated. A database focused on a Mach 4 transition speed with limited off-design elements was developed and will serve to guide future TBCC system studies and to validate higher level analyses.

Foster, Lancert E.↗

Brain segmentation and the generation of cortical surfaces

This paper describes methods for white matter segmentation in brain images and the generation of cortical surfaces from the segmentations. We have developed a system that allows a user to start with a brain volume, obtained by modalities such as MRI or cryosection, and constructs a complete digital representation of the cortical surface. The methodology consists of three basic components: local parametric modeling and Bayesian segmentation; surface generation and local quadratic coordinate fitting; and surface editing. Segmentations are computed by parametrically fitting known density functions to the histogram of the image using the expectation maximization algorithm [DLR77]. The parametric fits are obtained locally rather than globally over the whole volume to overcome local variations in gray levels. To represent the boundary of the gray and white matter we use triangulated meshes generated using isosurface generation algorithms [GH95]. A complete system of local parametric quadratic charts [JWM+95] is superimposed on the triangulated graph to facilitate smoothing and geodesic curve tracking. Algorithms for surface editing include extraction of the largest closed surface. Results for several macaque brains are presented comparing automated and hand surface generation. Copyright 1999 Academic Press.

NASA Discipline Neuroscience↗

Development of an Infrared Lamp Array for the Smap Spacecraft Thermal Balance Test

NASA launched the SMAP observatory in January 2015 aboard a Delta II into a sun-synchronous orbit around Earth. The science payload of a radar and a radiometer utilizes a shared rotating six-meter antenna to provide a global map of the Earth's soil moisture content and its freeze/thaw state on a global, high-resolution scale in this three-year mission. An observatory-level thermal balance test conducted in May/June 2014 validated the thermal design and demonstrated launch readiness as part of the planned environmental test campaign. An infrared lamp array was designed and used in the thermal balance test to replicate solar heating on the solar array and sunlit side of the spacecraft that would normally be seen in orbit. The design, implementation, and operation of an infrared lamp array used for this nineteen-day system thermal test are described in this paper. Instrumental to the smooth operation of this lamp array was a characterization test performed in the same chamber two months prior to the observatory test to provide insight into its array operation and flux uniformity. This knowledge was used to identify the lamp array power settings that would provide the worst case predicted on-orbit fluxes during eclipse, cold, and hot cases. It also showed the lamp array variation when adjustments in flux were needed. Calorimeters calibrated prior to testing determined a relationship between calorimeter temperature and lamp array flux. This allowed the team to adjust the lamp output for the desired absorbed flux on the solar array. Flux levels were within 10% of the desired value at the center of the solar array with an ability to maintain these levels within 5% during steady state cases. All tests demonstrated the infrared lamp array functionality and furthered lamp array understanding for modeling purposes. This method contributed to a high-fidelity environmental simulation, which was required to replicate the extreme on-orbit thermal environments.

Miller, Jennifer R.↗

Hydrodynamics of binary coalescence. 2: Polytropes with gamma = 5/3

We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.

Rasio, Frederic A.↗

Projection of Stabilized Aerial Imagery Onto Digital Elevation Maps for Geo-Rectified and Jitter-Free Viewing

As imagery is collected from an airborne platform, an individual viewing the images wants to know from where on the Earth the images were collected. To do this, some information about the camera needs to be known, such as its position and orientation relative to the Earth. This can be provided by common inertial navigation systems (INS). Once the location of the camera is known, it is useful to project an image onto some representation of the Earth. Due to the non-smooth terrain of the Earth (mountains, valleys, etc.), this projection is highly non-linear. Thus, to ensure accurate projection, one needs to project onto a digital elevation map (DEM). This allows one to view the images overlaid onto a representation of the Earth. A code has been developed that takes an image, a model of the camera used to acquire that image, the pose of the camera during acquisition (as provided by an INS), and a DEM, and outputs an image that has been geo-rectified. The world coordinate of the bounds of the image are provided for viewing purposes. The code finds a mapping from points on the ground (DEM) to pixels in the image. By performing this process for all points on the ground, one can "paint" the ground with the image, effectively performing a projection of the image onto the ground. In order to make this process efficient, a method was developed for finding a region of interest (ROI) on the ground to where the image will project. This code is useful in any scenario involving an aerial imaging platform that moves and rotates over time. Many other applications are possible in processing aerial and satellite imagery.

Ansar, Adnan I.↗

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗