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

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗

Injection and diffusive transport of suprathermal through energetic solar flare protons (35 keV to 20 MeV)

Consideration is given to the injection and interplanetary propagation of low-energy protons caused by the solar particle event of July 20, 1981, in which flare protons in the range from 35 keV to 20 MeV were observed by instruments on ISEE 3. The observed time-intensity and time-anisotropy profiles were fitted over the entire energy range using a model based on the spherically symmetric Fokker-Plank equation, including convection, diffusion, and adiabatic deceleration. The results are used to discuss the behavior of the radial interplanetary diffusion coefficient and the scattering mean free path for protons. Also, evidence is found for diffusive coronal shock acceleration of protons during the event.

Beeck, J.↗

An error covariance model for sea surface topography and velocity derived from TOPEX/POSEIDON altimetry

In order to facilitate the use of satellite-derived sea surface topography and velocity oceanographic models, methodology is presented for deriving the total error covariance and its geographic distribution from TOPEX/POSEIDON measurements. The model is formulated using a parametric model fit to the altimeter range observations. The topography and velocity modeled with spherical harmonic expansions whose coefficients are found through optimal adjustment to the altimeter range residuals using Bayesian statistics. All other parameters, including the orbit, geoid, surface models, and range corrections are provided as unadjusted parameters. The maximum likelihood estimates and errors are derived from the probability density function of the altimeter range residuals conditioned with a priori information. Estimates of model errors for the unadjusted parameters are obtained from the TOPEX/POSEIDON postlaunch verification results and the error covariances for the orbit and the geoid, except for the ocean tides. The error in the ocean tides is modeled, first, as the difference between two global tide models and, second, as the correction to the present tide model, the correction derived from the TOPEX/POSEIDON data. A formal error covariance propagation scheme is used to derive the total error. Our global total error estimate for the TOPEX/POSEIDON topography relative to the geoid for one 10-day period is found tio be 11 cm RMS. When the error in the geoid is removed, thereby providing an estimate of the time dependent error, the uncertainty in the topography is 3.5 cm root mean square (RMS). This level of accuracy is consistent with direct comparisons of TOPEX/POSEIDON altimeter heights with tide gauge measurements at 28 stations. In addition, the error correlation length scales are derived globally in both east-west and north-south directions, which should prove useful for data assimilation. The largest error correlation length scales are found in the tropics. Errors in the velocity field are smallest in midlatitude regions. For both variables the largest errors caused by uncertainty in the geoid. More accurate representations of the geoid await a dedicated geopotential satellite mission. Substantial improvements in the accuracy of ocean tide models are expected in the very near future from research with TOPEX/POSEIDON data.

Tsaoussi, Lucia S.↗

ROSAT HRI observations of Magellanic Cloud supernova remnants

Analysis of deep ROSAT high resolution imager (HRI) observations of two oxygen-rich supernova remnants (SNR's) in the Magellanic Clouds is described. For N132D, I exploit the limited spectral information provided by the HRI to investigate arcsecond scale spectral variations. I find that there is a region of harder X-ray emission near the southern limb and regions of softer emission near the center and northwestern limb. The remnant is believed to be interacting with a molecular cloud and the harder emission to the south is explained as a result of increased absorption along the line-of-sight there. I argue that the softer emission comes from X-ray emitting material with an enhanced abundance of oxygen. For the second SNR, E0102.2 72.2, the spatial structure is investigated in detail using two-dimensional image fitting techniques. Evidence is found for a ring-like and a spherically symmetric shell-like component both of which were modeled as homogeneous regions. In addition, a significant fraction of the observed flux (approximately 11 percent) must come from a resolved clumped component. A comparison with optical and radio imagery is made to provide a physical basis for the components identified in the X-ray analysis. The mass of X-ray emitting gas in the remnant is estimated and a value of approximately 75 M(solar mass) was determined. The dominant uncertainty on this quantity is the extent of unresolved clumping in the X-ray gas. Such clumping would tend to reduce the mass estimate by f(exp 1/2), where f is the mean volume filling factor of the gas.

Hughes, John P.↗

Method and apparatus for air-coupled transducer

An air-coupled transducer includes a ultrasonic transducer body having a radiation end with a backing fixture at the radiation end. There is a flexible backplate conformingly fit to the backing fixture and a thin membrane (preferably a metallized polymer) conformingly fit to the flexible backplate. In one embodiment, the backing fixture is spherically curved and the flexible backplate is spherically curved. The flexible backplate is preferably patterned with pits or depressions.

Song, Junho↗

Nearly steady flows in GONG prototype data

Doppler velocity images obtained with the GONG prototype instrument were analyzed to measure the nearly steady photospheric flows. The data consists of 88 images each of velocity, intensity, and modulation obtained at 20:00 UT on 88 days from July 1992 to February 1994. Each velocity image was temporally filtered to remove the p-mode oscillations, masked to exclude active regions, and then analyzed using spherical harmonics and orthogonal functions as described by Hathaway (1992). The spectral coefficients show very consistent results for the entire time interval with some evidence of year-to-year variations. The rotation profile agrees well with previous results and exhibits a north-south asymmetry that reverses sign during the 20 month interval. The residual rotation velocities exhibit structures with amplitudes of approximately 5 m/s that may be related to torsional oscillations. The meridional circulation is directed from the equator toward the poles with a peak velocity in the photosphere of approximately 50 m/s. The higher order components are very weak but indicate a divergent flow from the mid-latitudes (opposite that found for the June 1989 data). The convective limb shift is well fit by a 3rd order polynomial. The convection spectrum has a prominent peak at spherical harmonic degrees of l approximately 150 with very little signal in the low degree modes. Analysis of this signal shows that there is no evidence for giant cell convection at the level of approximately 10 m/s for all modes up to l = 32.

Hathaway, David H.↗

Modeling of the Earth's gravity field using the New Global Earth Model (NEWGEM)

Traditionally, the global gravity field was described by representations based on the spherical harmonics (SH) expansion of the geopotential. The SH expansion coefficients were determined by fitting the Earth's gravity data as measured by many different methods including the use of artificial satellites. As gravity data have accumulated with increasingly better accuracies, more of the higher order SH expansion coefficients were determined. The SH representation is useful for describing the gravity field exterior to the Earth but is theoretically invalid on the Earth's surface and in the Earth's interior. A new global Earth model (NEWGEM) (KIM, 1987 and 1988a) was recently proposed to provide a unified description of the Earth's gravity field inside, on, and outside the Earth's surface using the Earth's mass density profile as deduced from seismic studies, elevation and bathymetric information, and local and global gravity data. Using NEWGEM, it is possible to determine the constraints on the mass distribution of the Earth imposed by gravity, topography, and seismic data. NEWGEM is useful in investigating a variety of geophysical phenomena. It is currently being utilized to develop a geophysical interpretation of Kaula's rule. The zeroth order NEWGEM is being used to numerically integrate spherical harmonic expansion coefficients and simultaneously determine the contribution of each layer in the model to a given coefficient. The numerically determined SH expansion coefficients are also being used to test the validity of SH expansions at the surface of the Earth by comparing the resulting SH expansion gravity model with exact calculations of the gravity at the Earth's surface.

Kim, Yeong E.↗

The effect of data analysis techniques on the interpretation of wide-angle longwave radiation measurements

Three different data analysis techniques - shape factor, parameter estimation, and deconvolution - have been applied to the same set of satellite radiation measurements to determine their effect on the estimated radiation field. The measurements are from a wide-angle, horizon-to-horizon, nadir-pointing sensor. The shape factor technique reduces each measurement to a radiant exitance at the top of the atmosphere by simple division by a constant. The parameter estimation technique processes all measurements together as a batch and defines the radiant exitance as a least-squares fit to the data. The deconvolution technique takes advantage of the fact that spherical harmonics are the eigenfunctions of the measurement operator. All three techniques are derived, and their assumptions, advantages and disadvantages are discussed. Their results are compared globally, zonally, regionally and on a spatial spectrum basis. All three techniques give comparable results for global parameters; however, results on a regional scale were quite different. The standard deviations of the regional differences in radiant exitance varied from 7.4 to 13.5 W/sq m. Of the three techniques, the parameter estimation technique produced the best regional results and is the choice of the author.

Green, R. N.↗

Voyager 2 photopolarimeter observations of Titan

Observations of Titan's whole-disk polarization at 2460 and 7500 A are presented and analyzed in terms of model scattering atmospheres. If the Titan aerosols are spherical or nearly spherical, no single combination of refractive index and size distribution is able to fit data at both wavelengths. However, a vertically inhomogeneous distribution suggested by Tomasko and Smith (1980), characterized by a size gradient with altitude, fits the data at 2640 A moderately well but must be modified at intermediate and large optical depths to fit the 7500-A data. Results for synthetic phase functions indicate that the single-scattering polarization must be 70 percent or larger in the UV and 78 percent or larger in the near-IR at 90-deg phase angle, depending on the phase function. If the correct phase function is similar to that for 0.5-micron-radius spheres, the UV single-scattered polarization must be 84 percent and the near-IR single-scattered polarization must be over 90 percent. Such large polarizations are impossible for 0.5-micron-radius spheres but may be possible for nonspherical particles with effective radii near 0.5 micron, although the existence of nonspherical particles with the scattering properties required by these and other observations has not been demonstrated.

West, R. A.↗

Effects of surface pressures and streamline metrics on the calculation of laminar heating rates

The effect of streamline geometry and pressure distributions on surface heating rates is examined for slender, spherically blunted cones. The modifications to the approximate aeroheating code include a curve fit of pressures computed by an Euler solution over a range of Mach numbers and cone angles. The streamline geometry is then found using the surface pressures and inviscid surface properties. Previously, streamlines were determined using the inviscid properties at the edge of the boundary layer when accounting for the effects of entropy-layer swallowing. Streamline calculations are now based on inviscid surface conditions rather than boundary-layer edge properties. However, the heating rates are calculated using inviscid properties at the edge of the boundary layer. Resulting heating rates compare favorably with solutions from the viscous-shock-layer equations.

Riley, Christopher J.↗

Three-dimensional analysis of surface crack-Hertzian stress field interaction

The results are presented of a stress intensity factor analysis of semicircular surface cracks in the inner raceway of an engine bearing. The loading consists of a moving spherical Hertzian contact load and an axial stress due to rotation and shrink fit. A 3-D linear elastic Boundary Element Method code was developed to perform the stress analysis. The element library includes linear and quadratic isoparametric surface elements. Singular quarter point elements were employed to capture the square root displacement variation and the inverse square root stress singularity along the crack front. The program also possesses the capability to separate the whole domain into two subregions. This procedure enables one to solve nonsymmetric fracture mechanics problems without having to separate the crack surfaces a priori. A wide range of configuration parameters was investigated. The ratio of crack depth to bearing thickness was varied from one-sixtieth to one-fifth for several different locations of the Hertzian load. The stress intensity factors for several crack inclinations were also investigated. The results demonstrate the efficiency and accuracy of the Boundary Element Method. Moreover, the results can provide the basis for crack growth calculations and fatigue life prediction.

Ballarini, R.↗

Column Number Density Expressions Through M = 0 and M = 1 Point Source Plumes Along Any Straight Path

Analytical expressions for column number density (CND) are developed for optical line of sight paths through a variety of steady free molecule point source models including directionally-constrained effusion (Mach number M = 0) and flow from a sonic orifice (M = 1). Sonic orifice solutions are approximate, developed using a fair simulacrum fitted to the free molecule solution. Expressions are also developed for a spherically-symmetric thermal expansion (M = 0). CND solutions are found for the most general paths relative to these sources and briefly explored. It is determined that the maximum CND from a distant location through directed effusion and sonic orifice cases occurs along the path parallel to the source plane that intersects the plume axis. For the effusive case this value is exactly twice the CND found along the ray originating from that point of intersection and extending to infinity along the plume's axis. For sonic plumes this ratio is reduced to about 4/3. For high Mach number cases the maximum CND will be found along the axial centerline path. Keywords: column number density, plume flows, outgassing, free molecule flow.

Venting↗

Column Number Density Expressions Through M = 0 and M = 1 Point Source Plumes Along Any Straight Path

Analytical expressions for column number density (CND) are developed for optical line of sight paths through a variety of steady free molecule point source models including directionally-constrained effusion (Mach number M = 0) and flow from a sonic orifice (M 1). Sonic orifice solutions are approximate, developed using a fair simulacrum fitted to the free molecule solution. Expressions are also developed for a spherically-symmetric thermal expansion (M = 0). CND solutions are found for the most general paths relative to these sources and briefly explored. It is determined that the maximum CND from a distant location through directed effusion and sonic orifice cases occurs along the path parallel to the source plane that intersects the plume axis. For the effusive case this value is exactly twice the CND found along the ray originating from that point of intersection and extending to infinity along the plumes axis. For sonic plumes this ratio is reduced to about 43. For high Mach number cases the maximum CND will be found along the axial centerline path.

Column Density↗

A three-dimensional body-fitted coordinate system for flow field calculations on asymmetric nosetips

A three dimensional body-fitted coordinate system developed for use in the calculation of inviscid flows over ablated, asymmetric reentry vehicle nosetips is described. Because of the potential geometric asymmetries, no standard coordinate system (e.g., spherical, axisymmetric reference surface-normal) is capable of being closely aligned with the nosetip surface. To generate a 3-D, body-fitted coordinate system an analytic mapping procedure is applied that is conformal within each meridional plane of the nosetip; these transformations are then coupled circumferentially to yield a three dimensional coordinate system. The mappings used are defined in terms of hinge points, which are points selected to approximate the body contours in each meridional plane. The selection of appropriate hinge points was automated to facilitate the use of the resulting nosetip flow field code.

Hall, D. W.↗

Constraints on Porosity and Mass Loss in O-star Winds from the Modeling of X-ray Emission Line Profile Shapes

We fit X-ray emission line profiles in high resolution XMM-Newton and Chandra grating spectra of the early O supergiant Zeta Pup with models that include the effects of porosity in the stellar wind. We explore the effects of porosity due to both spherical and flattened clumps. We find that porosity models with flattened clumps oriented parallel to the photosphere provide poor fits to observed line shapes. However, porosity models with isotropic clumps can provide acceptable fits to observed line shapes, but only if the porosity effect is moderate. We quantify the degeneracy between porosity effects from isotropic clumps and the mass-loss rate inferred from the X-ray line shapes, and we show that only modest increases in the mass-loss rate (40%) are allowed if moderate porosity effects (h(sub infinity) less than approximately R(sub ∗)) are assumed to be important. Large porosity lengths, and thus strong porosity effects, are ruled out regardless of assumptions about clump shape. Thus, X-ray mass-loss rate estimates are relatively insensitive to both optically thin and optically thick clumping. This supports the use of X-ray spectroscopy as a mass-loss rate calibration for bright, nearby O stars

porosity↗

The Interstellar Polarization Feature Associated with the 2175 Angstrom Extinction Bump

The most successful model at fitting the wavelength dependence of interstellar extinction consists of two populations of bare silicate and graphite grains (Mathis, Rumpl & Nordsieck, (MRN)). The graphite grains are needed to fit the strong 2175 A extinction bump, with the silicates providing most of the smooth extinction seen at other wavelengths. From observations of the IR silicate absorption features, it was inferred that the silicate grains, non-spherical and aligned to the Galactic magnetic field, were responsible for the interstellar polarization. Aligned silicate grains make a very good fit to the smooth featureless UV polarization curves that comprise most of the observed sample. The lines of sight showing a UV polarization feature require something different, perhaps a second population of grains as do the fits to the extinction curve. The polarization wavelength dependence of HD 197770 can be fit by a MRN-like mixture of aligned bare silicate and graphite grains. However, none of the popular grain models including MRN envisaged an aligned bump grain population. It has been suggested that the extinction bump is not due to graphite but rather to a population of polycyclic aromatic hydrocarbons (PAHs). If the observed UV polarization features could be associated with the 2175 A bump then it would strengthen the argument that the bump must be due to grains rather than PAHs which are not likely to be aligned to the Galactic Magnetic field. The reality of ultraviolet polarization features has been confirmed by observations with ASTRO-2. The original detection toward HD 197770 has been confirmed and a new feature has been detected toward HD 147933-4. Both features have centroids lying close to 2175 A the location of the UV extinction bump. Two possibilities are considered for the source of the polarization bump, alignment of graphite grains responsible for the extinction bump or changes in the size distribution of the aligned silicate grains responsible for the continuum polarization. It seems the graphite grains are a more likely source of the polarization.

Clayton, Geoffrey C.↗

Examination of Multiple Lithologies Within the Primitive Ordinary Chondrite NWA 5717

Northwest Africa 5717 is a primitive (subtype 3.05) ungrouped ordinary chondrite which contains two apparently distinct lithologies. In large cut meteorite slabs, the darker of these, lithology A, looks to host the second, much lighter in color, lithology B (upper left, Fig. 1). The nature of the boundary between the two is uncertain, ranging from abrupt to gradational and not always following particle boundaries. The distinction between the lithologies, beyond the obvious color differences, has been supported by a discrepancy in oxygen isotopes and an incongruity in the magnesium contents of chondrule olivine. Here, quantitative textural analysis and mineralogical methods have been used to investigate the two apparent lithologies within NWA 5717. Olivine grains contained in a thin section from NWA 7402, thought to be paired to 5717, were also measured to re-examine the distinct compositional range among the light and dark areas. Procedure: Particles from a high-resolution mosaic image of a roughly 13x15cm slice of NWA 5717 were traced in Adobe Photoshop. Due to the large size of the sample, visually representative regions of each lithology were chosen to be analyzed. The resulting layers of digitized particles were imported into ImageJ, which was used to measure their area, along with the axes, the angle from horizontal, and the centroid coordinates of ellipses fitted to each particle following the approach. Resulting 2D pixel areas were converted to spherical diameters employing the unfolding algorithm, which outputs a 3D particle size distribution based on digitized 2D size frequency data. Spatstat was used to create kernel density plots of the centroid coordinates for each region. X-ray compositional maps, microprobe analyses, and Mossbauer spectroscopy was conducted on a thin section of NWA 7402, tentatively paired to NWA 5717.

Cato, M. J.↗

Plasma pressure in the environment of Jupiter, inferred from Voyager 1 magnetometer observations

A model combining the internal magnetic field with a self-consistent model of the Jovian magnetodisc was fitted to the Voyager 1 Jovian magnetic field data by means of a generalized inverse technique. The model parameters included the internal field spherical harmonic coefficients as well as with parameters describing the plasma distribution in the magnetosphere. Assuming that the pressure in the middle and outer magnetosphere is related to the unit flux tube volume V through PV exp gamma = const, the model fit yielded a value of 0.88 for gamma. If the hot (30 keV) plasma is transported adiabatically inward under the interchange instability triggered by centrifugal force of the heavy torus ions, losses are not sufficient to account for such a low value of gamma beyond L = 10. Closer to the planet, as the outer edge of the Io plasma torus is approached, PV exp gamma is found to decrease inward, as expected from the particle measurements, which identified an inner boundary of the particle fluxes in that region. With the present stage of the development of magnetodisc models, secular variations of the internal field still remain difficult to estimate.

Caudal, G.↗