Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical Fitting”

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 325 records · Page 18

Geoid-to-topography ratios on Venus: A global perspective

Recently available spherical harmonic solutions for the geoid and topography of Venus are sufficiently high resolution that they can be used to address questions concerning the relationship between geoid and topography on a regional scale. We have approached this question by mapping the geoid-to-topography ratio (GTR) on a systematic global basis. For a given point on the surface, we consider the geoid and elevation values at all points on a gridded representation of those fields located within a specified distance of the reference point. From the set of paired values, we determine the correlation coefficient and the best-fitting straight line. The latter is the GTR at that position, and the former is a measure of the significance of the derived ratio. This procedure is then repeated for all points on the global grid, yielding maps of the GTR and the correlation coefficient. Unlike previous studies of the GRT on Venus, this apprach permits us to make an objective and systematic search for regions with anomalous GTR's as well as areas that do not demonstrate any strong correlation between geoid and topography. These maps can be updated regularly as new harmonic models of the Venus geoid are produced from new Magellan tracking data. This procedure permits the development of a global perspective on the relationship between GTR and venusian surface tectonics.

Simons, Mark↗

Systematic characterization of unknown compounds via dimensionality reduction of time series

Analysis of ambient aerosols provides valuable insight into particle sources and formation chemistry. However, due to the complexity of atmospheric data and the dynamic nature of aerosol composition, a substantial fraction of data often become discarded by conventional analysis methods. Furthermore, a large fraction of chemical species within those data are unidentifiable due to a lack of matching spectral information, resulting in suboptimal characterization of chemical composition. Previous work has demonstrated techniques for cataloging analytes in a chromatographic dataset by deconvolution of mass spectra, but integration of these analytes throughout a large dataset remains time consuming. Here, we present a method to automatically identify an ion for quantitation for single-ion chromatogram based peak fitting and integration, enabling comprehensive integration of analytes with minimal user interaction. The resulting time series are clustered with a machine-learning based dimensionality reduction technique to systematically investigate the underlying characteristics of the categorized analytes and gain new insights into the chemical composition and physicochemical properties of the unidentifiable analytes. We apply these methods to existing atmospheric datasets collected in Manacapuru, Brazil during the GoAmazon2014/5 campaign to identify new analytes and interpret their variability and transformations in the atmosphere. The analysis results generate 408 time series from cataloged analytes of interest, and the clustering of those time series with spherical k-means results in 8 distinct clusters. We find the analytes form clusters based on their distinct physicochemical properties, demonstrating the method’s ability to systematically identify and selectively filter contaminants and instrumental analytes and characterize the unidentifiable analytes.

54 ENVIRONMENTAL SCIENCES↗

Active galactic nuclei. III - Accretion flow in an externally supplied cluster of black holes

This third paper in the series modeling QSOs and AGNs as clusters of accreting black holes studies the accretion flow within an externally supplied cluster. Significant radiation will be emitted by the cluster core, but the black holes in the outer halo, where the flow is considered spherically symmetric, will not contribute much to the overall luminosity of the source because of their large velocities relative to the infalling gas and therefore their small accretion radii. As a result, the scenario discussed in Paper I will refer to the cluster cores, rather than to entire clusters. This will steepen the high-frequency region of the spectrum unless inverse Compton scattering is effective. In many cases accretion flow in the central part of the cluster will be optically thick to electron scattering, resulting in a spectrum featuring optically thick radiative component in addition to power-law regimes. The fitting of these spectra to QSO and AGN observations is discussed, and application to 3C 273 is worked out as an example.

Pacholczyk, A. G.↗

Optical Testing Using Portable Laser Coordinate Measuring Instruments

High precision, portable coordinate measuring instruments (CMI) such as laser radars (LR) and laser trackers (LT) have been used for optical system alignment and integration. The LRs ability to perform a non-contact scan of surfaces was previously utilized to characterize large spherical and aspheric mirrors. In this paper, we explore the use of a CMI as an accurate, fast, robust, and non-contact tool for prescription characterization of powered optical surfaces. Using Nikons MV-224350 LR and Leicas Absolute Tracker AT401402 instruments, proof of concept measurements were performed to characterize a variety of optical components by measuring the actual and apparent, or equivalently the direct and through (DT), coordinates of calibrated metrology targets. Custom macros in metrology software and other data reduction code were developed to compute surface-ray intercepts and surface slopes from the DT shots. The calculated data is fit to an aspheric surface formula to obtain the optimum prescription. The results were compared to the nominal parameters and were crosschecked using LR scans or other approaches. We discuss potential applications across the fields of optical component fabrication and system alignment and testing.

Khreishi, M.↗

Evaluation of spacecraft navigation using the Tracking and Data Relay Satellite System (TDRSS)

The navigational capabilities of the TDRSS for the operational support of user spacecraft is evaluated. A batch-weighted least-squares algorithm is used to fit tracking measurements and to generate orbit solutions for the TDRS-E spacecraft and the following user spacecraft: (1) the SMM, (2) the Landsat-5, (3) the Earth Radiation Budget Satellite (ERBS), and (4) the Solar Mesosphere Explorer (SME). TDRS-E orbit accuracy was evaluated by using the consistency of consecutive orbit solutions as a measure; user orbit accuracy was estimated by comparing orbits determined from TDRSS tracking with orbits obtained from ground tracking for the same time periods. It was found that the predicted accuracy of the TDRS-E orbit propagation was improved by as much as seven times by refining the spherical harmonic expansion model of the earth's potential field. The results obtained for SMM were fairly consistent, while those obtained for Landsat-5 were consistent to better than 80 meters. The orbits obtained using ERBS tracking were consistent to better than 60 meters and the SME orbits for two arcs with 14 and 20 passes of TDRSS tracking had a maximum position difference of 140 meters.

Samii, M. V.↗

GRGM900C: A Degree 900 Lunar Gravity Model from GRAIL Primary and Extended Mission Data

We have derived a gravity field solution in spherical harmonics to degree and order 900, GRGM900C, from the tracking data of the Gravity Recovery and Interior Laboratory (GRAIL) Primary (1 March to 29 May 2012) and Extended Missions (30 August to 14 December 2012). A power law constraint of 3.6 × 10(exp −4)/l(exp 2) was applied only for degree l greater than 600. The model produces global correlations of gravity, and gravity predicted from lunar topography of greater than or equal to 0.98 through degree 638. The model's degree strength varies from a minimum of 575-675 over the central nearside and farside to 900 over the polar regions. The model fits the Extended Mission Ka-Band Range Rate data through 17 November 2012 at 0.13 micrometers/s RMS, whereas the last month of Ka-Band Range-Rate data obtained from altitudes of 2-10 km fit at 0.98 micrometers/s RMS, indicating that there is still signal inherent in the tracking data beyond degree 900.

Gravity↗

Large Volcanic Rises on Venus

Large volcanic rises on Venus have been interpreted as hotspots, or the surface manifestation of mantle upwelling, on the basis of their broad topographic rises, abundant volcanism, and large positive gravity anomalies. Hotspots offer an important opportunity to study the behavior of the lithosphere in response to mantle forces. In addition to the four previously known hotspots, Atla, Bell, Beta, and western Eistla Regiones, five new probable hotspots, Dione, central Eistla, eastern Eistla, Imdr, and Themis, have been identified in the Magellan radar, gravity and topography data. These nine regions exhibit a wider range of volcano-tectonic characteristics than previously recognized for venusian hotspots, and have been classified as rift-dominated (Atla, Beta), coronae-dominated (central and eastern Eistla, Themis), or volcano-dominated (Bell, Dione, western Eistla, Imdr). The apparent depths of compensation for these regions ranges from 65 to 260 km. New estimates of the elastic thickness, using the 90 deg and order spherical harmonic field, are 15-40 km at Bell Regio, and 25 km at western Eistla Regio. Phillips et al. find a value of 30 km at Atla Regio. Numerous models of lithospheric and mantle behavior have been proposed to interpret the gravity and topography signature of the hotspots, with most studies focusing on Atla or Beta Regiones. Convective models with Earth-like parameters result in estimates of the thickness of the thermal lithosphere of approximately 100 km. Models of stagnant lid convection or thermal thinning infer the thickness of the thermal lithosphere to be 300 km or more. Without additional constraints, any of the model fits are equally valid. The thinner thermal lithosphere estimates are most consistent with the volcanic and tectonic characteristics of the hotspots. Estimates of the thermal gradient based on estimates of the elastic thickness also support a relatively thin lithosphere (Phillips et al.). The advantage of larger estimates of the thermal lithospheric thickness is that they provide an explanation for the apparently modest levels of geologic activity on Venus over the last half billion years.

Smrekar, Suzanne E.↗

Optical Testing Using Portable Laser Coordinate Measuring Instruments

High precision, portable coordinate measuring instruments (CMI) such as laser radars (LR) and laser trackers (LT) have been used for optical system alignment and integration. The LR's ability to perform a non-contact scan of surfaces was previously utilized to characterize large spherical and aspheric mirrors. In this paper, we explore the use of a CMI as an accurate, fast, robust, and non-contact tool for prescription characterization of powered optical surfaces. Using Nikon's MV-224/350 LR and Leica's Absolute Tracker AT401/402 instruments, proof of concept measurements were performed to characterize a variety of optical components by measuring the actual and apparent, or equivalently the "direct and through" (D&T), coordinates of calibrated metrology targets. Custom macros in metrology software and other data reduction code were developed to compute surface-ray intercepts and surface slopes from the D&T shots. The calculated data is fit to an aspheric surface formula to obtain the optimum prescription. The results were compared to the nominal parameters and were crosschecked using LR scans or other approaches. We discuss potential applications across the fields of optical component fabrication and system alignment and testing.

Khreishi, Manal↗

On the penetration of a hot diapir through a strongly temperature-dependent viscosity medium

The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 0.1 - 1000 from constant viscosity up to viscosity variations of 100,000. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 10 to the 22nd sq cm/s in order for a 10 km diapir to penetrate a distance of several radii.

Daly, S. F.↗

Global Surface Temperature Change and Uncertainties Since 1861

The objective of this talk is to analyze the warming trend and its uncertainties of the global and hemi-spheric surface temperatures. By the method of statistical optimal averaging scheme, the land surface air temperature and sea surface temperature observational data are used to compute the spatial average annual mean surface air temperature. The optimal averaging method is derived from the minimization of the mean square error between the true and estimated averages and uses the empirical orthogonal functions. The method can accurately estimate the errors of the spatial average due to observational gaps and random measurement errors. In addition, quantified are three independent uncertainty factors: urbanization, change of the in situ observational practices and sea surface temperature data corrections. Based on these uncertainties, the best linear fit to annual global surface temperature gives an increase of 0.61 +/- 0.16 C between 1861 and 2000. This lecture will also touch the topics on the impact of global change on nature and environment. as well as the latest assessment methods for the attributions of global change.

Shen, Samuel S. P.↗

A stochastic-dynamic model for global atmospheric mass field statistics

A model that yields the spatial correlation structure of atmospheric mass field forecast errors was developed. The model is governed by the potential vorticity equation forced by random noise. Expansion in spherical harmonics and correlation function was computed analytically using the expansion coefficients. The finite difference equivalent was solved using a fast Poisson solver and the correlation function was computed using stratified sampling of the individual realization of F(omega) and hence of phi(omega). A higher order equation for gamma was derived and solved directly in finite differences by two successive applications of the fast Poisson solver. The methods were compared for accuracy and efficiency and the third method was chosen as clearly superior. The results agree well with the latitude dependence of observed atmospheric correlation data. The value of the parameter c sub o which gives the best fit to the data is close to the value expected from dynamical considerations.

Ghil, M.↗

Balloon for Long-Duration, High-Altitude Flight at Venus

A document describes a 5.5-m-diameter, helium-filled balloon designed for carrying a scientific payload having a mass of 44 kg for at least six days at an altitude of about 55 km in the atmosphere of Venus. The requirement for floating at nearly constant altitude dictates the choice of a mass-efficient spherical super-pressure balloon that tracks a constant atmospheric density. Therefore, the balloon is of a conventional spherical super-pressure type, except that it is made of materials chosen to minimize solar radiant heating and withstand the corrosive sulfuric acid aerosol of the Venusian atmosphere. The shell consists of 16 gores of a multilayer composite material. The outer layer, made of polytetrafluoroethylene, protects against sulfuric acid aerosol. Next is an aluminum layer that reflects sunlight to minimize heating, followed by an aluminized polyethylene terephthalate layer that resists permeation by helium, followed by an aromatic polyester fabric that imparts strength to withstand deployment forces and steady super-pressure. A polyurethane coat on the inner surface of the fabric facilitates sealing at gore-to-gore seams. End fittings and seals, and a tether connecting the end fittings to a gondola, are all made of sulfuric-acid-resistant materials.

Hall, Jeffrey↗

Modified Isotherm Modeling to Predict Cs Exchange with Crystalline Silicotitanate in Tank Waste Simulants

The U.S. Department of Energy is working to expedite processing of Hanford tank waste supernate at the Hanford Waste Treatment and Immobilization Plant. To support this goal, Washington River Protection Solutions is designing a Tank Side Cesium Removal (TSCR) system for suspended solids and cesium (Cs/ 137 Cs) removal from Hanford tank waste supernate. The ion exchange media selected for Cs removal at TSCR is crystalline silicotitanate (CST) that is manufactured in a nearly spherical form by Honeywell UOP LLC (UOP; Des Plaines, IL) as product IONSIV® R9140-B (Na form). The Zheng Anthony Miller (ZAM) isotherm model (Zheng et al. 1997) is a multicomponent ion exchange model used to predict the exchange of Group I metals onto CST. The ZAM isotherm has historically been used to predict Cs distribution values from Hanford and Savannah River Site (SRS) tank waste simulants. Figure S.1 summarizes model predictions from the ZAM isotherm that indicate poor prediction of Cs distribution values for simple and complex simulants with the engineered form of CST, IONSIV® R9140-B, and IONSIV® R9120-B where the solid line indicates a perfect fit by the model. The dotted lines indicate ±20% error. Batch contact testing with Hanford tank waste complex and simple simulants was used in conjunction with SRS simulants to experimentally determine Cs distribution values using a modification to the original isotherm model. The experimentally determined maximum Cs capacity for IONSIV® R9140-B CST in both the simple and complex matrices was found to be 0.53±0.3 mmoles Cs/g of CST. This value is not drastically different from the maximum Cs capacity of 0.58 mmoles Cs/g TAM-5 reported by Zheng et al. (1997). However, it is important to note that TAM-5 (commercially IONSIV® IE-910) is a powder. Hamm et al. (2002) determined that a dilution factor was needed to account for the Zr(OH)2 binder in the engineered form of CST. Hamm et al. determined that a dilution factor of 0.68 was appropriate to account for binder contribution and correct overprediction of Cs exchange on the engineered form of CST in ZAM calculations. This reduced the total capacity from 0.58 mmol/g with TAM-5 to 0.39 mmol/g for the engineered form of CST (Hamm et al. 2002). Despite substituting the experimentally determined maximum Cs capacity of 0.55 mmoles Cs/g for the literature-reported capacity of 0.39 mmoles Cs/g, it was determined that additional modifications to the model’s equilibrium rate constants were necessary in refining the isotherm model. The modified model overpredicted K+ uptake by the CST when compared to digested CST results described by Campbell et al. (2019). The modified model was further revised to omit three of the five K+ exchange equilibrium reactions described by ZAM to reduce the additional K+ loading seen by the model. Figure S.2 summarizes the revised model isotherm predictions plotted against measured Kd values.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

On the prediction of velocity fields from redshift space galaxy samples

We present a new method for recovering the underlying velocity field from an observed distribution of galaxies in redshift space. The method is based on a kinematic Zel'dovich relation between the velocity and density fields in redshift space. This relation is expressed in a differential equation slightly modified from the usual Poisson equation and which depends nontrivially on Beta identically equal to Omega(exp 0.6)/b. The linear equation can be readily solved by standard techniques of separation of variables by means of spherical harmonics. One can also include a term describing the 'rocket effect' discussed by Kaiser (1987). From this redshift space information alone, one can generate a prediction of the peculiar velocity field for each harmonic (l, m) as a function of distance. We note that for the quadrupole and higher order moments, the equation is a boundary value problem with solutions dependent on both the interior and exterior mass distribution. However, for a shell at distance r, the dipole, as well as the monopole, of the velocity field in the Local Group frame is fully determined by the interior mass distribution. This implies that the shear of the measured velocity field, when fitted to a dipole distortion, should be aligned and consistent with the gravity field inferred from the well determined local galaxy distribution. As a preliminary application we compute the velocity dipole of distant shells as predicted from the 1.2 Jy IRAS survey compared to the measured velocity dipole on shells, as inferred from a recent POTENT analysis. The coherence between the two fields is good, yielding a best estimate of Beta = 0.6 +/- 0.2.

Nusser, Adi↗

Seyfert galaxy narrow-line regions. II - Kinematic models

A class of kinematic models for the narrow-line regions of active galactic nuclei is presented. The models incorporate continuity-constrained radial flow of emitting material in a bounded, spherical geometry. They are designed to represent the narrow-line regions in as simple a manner as is consistent with current knowledge. The systematic comparison of line profiles arising from the models to the spectra in the forbidden O III lambda 5007 emission line which have been described in Paper I of this series is used as a probe of the line-emitting region kinematics: the model parameter areas of interest are identified and tested for consistency with the data. The results indicate heterogeneous kinematic conditions; the most common pattern is outflow of emitting material in a radially decelerating velocity field in the presence of distributed obscuration. However, some observed line profiles can be fitted by several model profiles from very different kinematic conditions, indicating that the information in a spatially unresolved line profile is insufficient to constrain uniquely the geometry and kinematics of the line-emitting region. Although the specific conclusions are model dependent, a predominantly radial pattern of velocities and a distribution of obscuring material are probably essential features of any successful kinematic model of the narrow-line regions.

Vrtilek, J. M.↗

Quasi-brittle porous material: Simulated effect of stochastic air void structure on compressive strength

Highlights: • A simulation procedure to study the effect of porosity on the compressive strength of a quasi-brittle material is introduced • Based on the simulated results, the relationship between compressive strength and air void fraction (porosity) is presented. • A linear relationship was achieved between the cubic root of porosity and the simulated compressive strength. The effect of porosity comprised of spherical air voids on the compressive strength of quasi-brittle material was studied via simulations. The simulated porous structures were based on pore size distributions of two mortar samples measured by X-ray microtomography. While the simulation method set practical limits on the size of sample, the base of the statistics was established by simulating 128 small structures generated by sampling from pore structures of two mortars. By studying the application of the classical strength-porosity formulas to the simulated data, a new simple model was formed. A linear relationship was achieved between the cubic root of air void fraction (porosity) and the simulated compressive strength. The reasons for scattering of simulated strength around fitted trend remained unresolved in this study; no clear dependence on pore number or other distribution properties was observed. With the presented simulation approach, the dependence of compressive strength on porosity is achieved independently of disturbances that occur in experimental studies creating understanding of compressional behavior of low porosity materials.

36 MATERIALS SCIENCE↗

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗

Low-Energy Sputtering Research

An experimental study is described to measure low-energy (less than 600 eV) sputtering yields of molybdenum with xenon ions using Rutherford backscattering spectroscopy (RBS) and secondary neutral mass spectroscopy (SNMS). An ion gun was used to generate the ion beam. The ion current density at the target surface was approximately 30 (micro)A/sq cm. For RBS measurements, the sputtered material was collected on a thin aluminum strip which was mounted on a semi-circular collector plate. The target was bombarded with 200 and 500 eV xenon ions at normal incidence. The differential sputtering yields were measured using the RBS method with 1 MeV helium ions. The differential yields were fitted with a cosine fitting function and integrated with respect to the solid angle to provide the total sputtering yields. The sputtering yields obtained using the RBS method are in reasonable agreement with those measured by other researchers using different techniques. For the SNMS measurements, 150 to 600 eV xenon ions were used at 50deg angle of incidence. The SNMS spectra were converted to sputtering yields for perpendicular incidence by normalizing SNMS spectral data at 500 eV with the yield measured by Rutherford backscattering spectrometry. Sputtering yields as well as the shape of the yield-energy curve obtained in this manner are in reasonable agreement with those measured by other researchers using different techniques. Sputtering yields calculated by using two semi-spherical formulations agree reasonably well with measured data. The isotopic composition of secondary ions were measured by bombarding copper with xenon ions at energies ranging from 100 eV to 1.5 keV. The secondary ion flux was found to be enriched in heavy isotopes at low incident ion energies. The heavy isotope enrichment was observed to decrease with increasing impact energy. Beyond 700 eV, light isotopes were sputtered preferentially with the enrichment remaining nearly constant.

Ray, P. K.↗