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 109 records · Page 6

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.↗

The C IV 1550 profile in type 1 Seyfert galaxies

The paper presents C IV 1550 A line profiles for the type 1 Seyfert galaxies NGC 5548, Mrk 509, NGC 7469, and MCG-2-58-22. Several line broadening mechanisms and theoretical line profiles are considered, and random motion of discrete clouds is ruled out. A spherical ensemble of discrete clouds with steady outflow or inflow produces a logarithmic profile, but does not account for the highly extended wings. A spherical ensemble with ballistic outflow produces a profile of the first exponential integral function, and fits the observed profile to the continuum level. Although C IV profiles favor the ballistic model, both Mrk 509 and NGC 7469 have significant asymmetry, and Balmer lines with a higher optical depth show higher asymmetry and redshift than all four galaxies.

Wu, C.-C.↗

Computation of nonlinear supersonic potential flow over three dimensional surfaces

A fully implicit three dimensional marching technique is presented for the computation of inviscid supersonic flows. The nonconservative full potential equation is solved on spherical cross flow planes using standard transonic SLOR techniques. The three dimensional bow shock is fitted as a boundary with isentropic jump conditions and embedded shocks are captured. Applications of these techniques are presented for three dimensional bodies, wings, and wing-body configurations and the computed solutions are compared to experimental pressure data. Zero lift wave drag calculations are also presented and compared with measurements to evaluate the accuracy of the present techniques.

Siclari, M. J.↗

Galactic Cosmic-Ray Energy Spectra and Composition during the 2009-2010 Solar Minimum Period

We report new measurements of the elemental energy spectra and composition of galactic cosmic rays during the 2009-2010 solar minimum period using observations from the Cosmic Ray Isotope Spectrometer (CRIS) onboard the Advanced Composition Explorer. This period of time exhibited record-setting cosmic-ray intensities and very low levels of solar activity. Results are given for particles with nuclear charge 5 <= Z <= 28 in the energy range approx. 50-550 MeV / nucleon. Several recent improvements have been made to the earlier CRIS data analysis, and therefore updates of our previous observations for the 1997-1998 solar minimum and 2001-2003 solar maximum are also given here. For most species, the reported intensities changed by less than approx. 7%, and the relative abundances changed by less than approx. 4%. Compared with the 1997-1998 solar minimum relative abundances, the 2009-2010 abundances differ by less than 2sigma, with a trend of fewer secondary species observed in the more recent time period. The new 2009-2010 data are also compared with results of a simple "leaky-box" galactic transport model combined with a spherically symmetric solar modulation model. We demonstrate that this model is able to give reasonable fits to the energy spectra and the secondary-to-primary ratios B/C and (Sc+Ti+V)/Fe. These results are also shown to be comparable to a GALPROP numerical model that includes the effects of diffusive reacceleration in the interstellar medium.

Galactic Cosmic Rays↗

Elemental GCR Observations during the 2009-2010 Solar Minimum Period

Using observations from the Cosmic Ray Isotope Spectrometer (CRIS) onboard the Advanced Composition Explorer (ACE), we present new measurements of the galactic cosmic ray (GCR) elemental composition and energy spectra for the species B through Ni in the energy range approx. 50-550 MeV/nucleon during the record setting 2009-2010 solar minimum period. These data are compared with our observations from the 1997-1998 solar minimum period, when solar modulation in the heliosphere was somewhat higher. For these species, we find that the intensities during the 2009-2010 solar minimum were approx. 20% higher than those in the previous solar minimum, and in fact were the highest GCR intensities recorded during the space age. Relative abundances for these species during the two solar minimum periods differed by small but statistically significant amounts, which are attributed to the combination of spectral shape differences between primary and secondary GCRs in the interstellar medium and differences between the levels of solar modulation in the two solar minima. We also present the secondary-to-primary ratios B/C and (Sc+Ti+V)/Fe for both solar minimum periods, and demonstrate that these ratios are reasonably well fit by a simple "leaky-box" galactic transport model that is combined with a spherically symmetric solar modulation model.

Solar Modulation↗

Interpretation of recent positron-electron measurements between 20 and 800 MeV

Recently measured positron and negatron spectra are discussed with regard to the problem of solar modulation. At energies above 180 MeV, the spherically symmetric Fokker-Planck equation with a diffusion coefficient proportional to particle rigidity provides reasonable fits to both the positron and total electron data. At energies below 180 MeV, the data are consistent with a continuation of the same diffusion coefficient and a local source of negatrons or with a change in the diffusion coefficient to a constant value.

Pellerin, C. J.↗