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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

RIEGER-TYPE PERIODICITY DURING SOLAR CYCLES 14–24: ESTIMATION OF DYNAMO MAGNETIC FIELD STRENGTH IN THE SOLAR INTERIOR

Solar activity undergoes a variation over timescales of several months known as Rieger-type periodicity, which usually occurs near maxima of sunspot cycles. An early analysis showed that the periodicity appears only in some cycles and is absent in other cycles. But the appearance/absence during different cycles has not been explained. We performed a wavelet analysis of sunspot data from the Greenwich Royal Observatory and the Royal Observatory of Belgium during cycles 14–24. We found that the Rieger-type periods occur in all cycles, but they are cycle dependent: shorter periods occur during stronger cycles. Our analysis revealed a periodicity of 185–195 days during the weak cycles 14–15 and 24 and a periodicity of 155–165 days during the stronger cycles 16–23. We derived the dispersion relation of the spherical harmonics of the magnetic Rossby waves in the presence of differential rotation and a toroidal magnetic field in the dynamo layer near the base of the convection zone. This showed that the harmonics of fast Rossby waves with m = 1 and n = 4, where m ( n ) indicates the toroidal (poloidal) wavenumbers, perfectly fit with the observed periodicity. The variation of the toroidal field strength from weaker to stronger cycles may lead to the different periods found in those cycles, which explains the observed enigmatic feature of the Rieger-type periodicity. Finally, we used the observed periodicity to estimate the dynamo field strength during cycles 14–24. Our estimations suggest a field strength of ∼40 kG for the stronger cycles and ∼20 kG for the weaker cycles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Unseeded, spontaneous nucleation of spherulitic magnesium calcite

Most of the sedimentary carbonates deposited in the marine environments are composed of calcium carbonate minerals with varying amounts of incorporated Mg 2+ . However, understanding how interactions of impurities with carbonate and their incorporation affect sediments behavior remains a challenge. Here, a new insight is obtained by monitoring solution composition, morphology, and electrokinetic potential of carbonate particles formed in a spontaneous unseeded batch precipitation experiment using electrochemical and scanning electron microscopy methods. The solid composition and growth rate are extracted from changes in the bulk composition and fitted to chemical affinity rate law, revealing that the precipitation pathway consists of second-order dissolution and first-order precipitation. The molecular dynamics simulations show that the lattice strain induced by randomly substituting Ca 2+ by Mg 2+ stabilizes spherical nanoparticles and reduces their surface area and volume. Combining kinetics and thermodynamics insight, we conclude that variation in the carbonate bulk and interfacial energies, along with the solution supersaturation, lead to the dissolution-precipitation transformation pathway from Mg-rich to Mg-poor carbonate phase that preserves spherulitic morphology. Furthermore, our findings are relevant for long-standing questions of how impurities influence diagenesis of carbonate sediments and spherulitic carbonate particles' origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Local models for scatter estimation and descattering in polyenergetic X-ray tomography

We propose a new modeling approach for scatter estimation and descattering in polyenergetic X-ray computed tomography (CT) based on fitting models to local neighborhoods of a training set. X-ray CT is widely used in medical and industrial applications. X-ray scatter, if not accounted for during reconstruction, creates a loss of contrast in CT reconstructions and introduces severe artifacts including cupping, shading, and streaks. Even when these qualitative artifacts are not apparent, scatter can pose a major obstacle in obtaining quantitatively accurate reconstructions. Our approach to estimating scatter is, first, to generate a training set of 2D radiographs with and without scatter using particle transport simulation software. To estimate scatter for a new radiograph, we adaptively fit a scatter model to a small subset of the training data containing the radiographs most similar to it. We compared local and global (fit on full data sets) versions of several X-ray scatter models, including two from the recent literature, as well as a recent deep learning-based scatter model, in the context of descattering and quantitative density reconstruction of simulated, spherically symmetrical, single-material objects comprising shells of various densities. Our results show that, when applied locally, even simple models provide state-of-the-art descattering, reducing the error in density reconstruction due to scatter by more than half.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1997

Spherical Harmonics Analysis of the ECMWF Global Wind Fields at the 10-Meter Height Level During 1985: A Collection of Figures Illustrating Results

Half-daily global wind speeds in the east-west (u) and north-south (v) directions at the 10-meter height level were obtained from the European Centre for Medium Range Weather Forecasts (ECMWF) data set of global analyses. The data set covered the period 1985 January to 1995 January. A spherical harmonic expansion to degree and order 50 was used to perform harmonic analysis of the east-west (u) and north-south (v) velocity field components. The resulting wind field is displayed, as well as the residual of the fit, at a particular time. The contribution of particular coefficients is shown. The time variability of the coefficients up to degree and order 3 is presented. Corresponding power spectrum plots are given. Time series analyses were applied also to the power associated with degrees 0-10; the results are included.

Sanchez, Braulio V.↗
DOE OSTI2020

Resistive heating in an electrified domain with a spherical inclusion: an ALEGRA verification study

A verification study is conducted for the ALEGRA software, using the problem of an electrified medium with a spherical inclusion, paying special attention to resistive heating. We do so by extending an existing analytic solution for this problem to include both conducting and insulating inclusions, and we examine the effects of mesh resolution and mesh topology, considering both body-fitted and rectangular meshes containing mixed cells. We present observed rates of convergence with respect to mesh refinement for four electromagnetic quantities: electric potential, electric field, current density and Joule power.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Atom-partitioned multipole expansions for electrostatic potential boundary conditions

Applications such as grid-based real-space density functional theory (DFT) use the Poisson equation to compute electrostatics. However, the expected long tail of the electrostatic potential requires either the use of a large and costly outer domain or Dirichlet boundary conditions estimated via multipole expansion. We find that the oft-used single-center spherical multipole expansion is only appropriate for isotropic mesh domains such as spheres and cubes. In this work, we introduce a method suitable for high aspect ratio meshes whereby the charge density is partitioned into atomic domains and multipoles are computed for each domain. While this approach is moderately more expensive than a single-center expansion, it is numerically stable and still a small fraction of the overall cost of a DFT calculation. The net result is that when high aspect ratio systems are being studied, form-fitted meshes can now be used in lieu of cubic meshes to gain computational speedup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Some requirements of a colliding comet source of gamma ray bursts

Colliding comets in the Solar System may be an important source of gamma ray bursts. The spherical gamma ray comet cloud required by the results of the Venera Satellites (Mazets and Golenetskii, 1987) and the Burst and Transient Source Experiment (BATSE) detector on the Compton Satellite (Meegan et al., 1992a, b) is neither the Oort Cloud nor the Kuiper Belt. To satisfy observations of N(greater than P(sub max)) vs P(sub max) for the maximum gamma ray fluxes, P(sub max) greater than 10(exp -5) ergs/sq cm/ s (about 30 bursts/yr), the comet density, n, should increase as n approximately a(exp 1) from about 40 to 100 AU where a is the comet heliocentric distance. The turnover above 100 AU requires n approximately a(exp -1/2) to 200 AU to fit the Venera results and n approximately a(exp 1/4) to 400 AU to fit the BATSE data. Then the masses of comets in the 3 regions are from: 40-100 AU, about 9 earth masses, m(sub E); 100-200 AU about 25 m(sub E); and 100-400 AU, about 900 m(sub E). The flux of 10(exp -5) ergs/sq cm/s corresponds to a luminosity at 100 AU of 3 x 10(exp 26) ergs/s. Two colliding spherical comets at a distance of 100 AU, each with nucleus of radius R of 5 km, density of 0.5 g/cu cm and Keplerian velocity 3 km/s have a combined kinetic energy of 3 x 10(exp 28) erg, a factor of about 100 greater tha required by the burst maximum fluxes that last for one second.

White, R. Stephen↗
NASA NTRS1974

A stabilized large-aperture far-infrared telescope gondola

Description of a recently developed balloon-borne 1-m far-infrared telescope designed to carry out astronomical observations and to map radiation from selected areas of the sky. The telescope optical system consists of a conventional Cassegrain arrangement with an f/2 aluminum primary mirror, spherically figured, which feeds a pyrex secondary mirror to produce an f/13.8 beam at the focal plane with a scale of 15.2 arc sec/mm. Three focal planes are provided. The infrared detectors consist of an array of four gallium-doped germanium bolometers fitted within a two-liter helium dewar. An altitude-azimuth mount design is used for attitude control, with the pointing objectives being achieved in two steps. The azimuth control system is capable of operation in a position mode, in which the sensing element is a servo-driven null magnetometer always pointed north, and in an inertial mode, in which the gyro is uncaged and becomes the azimuth reference. The elevation system also operates in both position and inertial modes. An auxiliary position control, powered by an independent battery and controlled by a separate tone command system, is provided in case of a primary system difficulty.

Hazen, N. L.↗
NASA NTRS1994

Improved Bounds on Non-Luminous Matter in Solar Orbit

We improve, using a larger set of observations including Voyager 2 Neptune flyby data, previous bounds on the amount of dark matter (DM) trapped in a spherically symmetric distribution about the sun. We bound DM by noting that such a distribution would increase the effective mass of the sun as seen by the outer planets and by finding the uncertainty in that effective mass for Uranus and Neptune in fits to the JPL Developmental Ephermeris residuals, including optical data and those two planets' Voyager 2 flybys. We extend our previous procedure by fitting more parameters of the developmental ephemerides. Additionally, we present here the values for Pioneer 10 and 11 and Voyager 1 and 2 Jupiter ranging normal points (and incorporate these data as well).

dark matter Deep Space Network DSN Voyager Neptune↗
DOE OSTI2021

Discovery of multiple p-mode pulsation frequencies in the roAp star, HD 86181

ABSTRACT We report the frequency analysis of a known roAp star, HD 86181 (TIC 469246567), with new inferences from Transiting Exoplanet Survey Satellite (TESS) data. We derive the rotation frequency to be νrot = 0.48753 ± 0.00001 d−1. The pulsation frequency spectrum is rich, consisting of two doublets and one quintuplet, which we interpret to be oblique pulsation multiplets from consecutive, high-overtone dipole, quadrupole, and dipole modes. The central frequency of the quintuplet is 232.7701 d−1 (2.694 mHz). The phases of the sidelobes, the pulsation phase modulation, and a spherical harmonic decomposition all show that the quadrupole mode is distorted. Following the oblique pulsator model, we calculate the rotation inclination, i, and magnetic obliquity, β, of this star, which provide detailed information about the pulsation geometry. The i and β derived from the best fit of the pulsation amplitude and phase modulation to a theoretical model, including the magnetic field effect, slightly differ from those calculated for a pure quadrupole, indicating the contributions from ℓ = 4, 6, 8,... are small. Non-adiabatic models with different envelope convection conditions and physics configurations were considered for this star. It is shown that models with envelope convection almost fully suppressed can explain the excitation at the observed pulsation frequencies.

Shi, Fangfei (ORCID:0000000323998216)↗
DOE OSTI2016

NEUTRON STAR MASS–RADIUS CONSTRAINTS USING EVOLUTIONARY OPTIMIZATION

The equation of state of cold supra-nuclear-density matter, such as in neutron stars, is an open question in astrophysics. A promising method for constraining the neutron star equation of state is modeling pulse profiles of thermonuclear X-ray burst oscillations from hot spots on accreting neutron stars. The pulse profiles, constructed using spherical and oblate neutron star models, are comparable to what would be observed by a next-generation X-ray timing instrument like ASTROSAT , NICER , or a mission similar to LOFT . In this paper, we showcase the use of an evolutionary optimization algorithm to fit pulse profiles to determine the best-fit masses and radii. By fitting synthetic data, we assess how well the optimization algorithm can recover the input parameters. Multiple Poisson realizations of the synthetic pulse profiles, constructed with 1.6 million counts and no background, were fitted with the Ferret algorithm to analyze both statistical and degeneracy-related uncertainty and to explore how the goodness of fit depends on the input parameters. For the regions of parameter space sampled by our tests, the best-determined parameter is the projected velocity of the spot along the observer’s line of sight, with an accuracy of ≤3% compared to the true value and with ≤5% statistical uncertainty. The next best determined are the mass and radius; for a neutron star with a spin frequency of 600 Hz, the best-fit mass and radius are accurate to ≤5%, with respective uncertainties of ≤7% and ≤10%. The accuracy and precision depend on the observer inclination and spot colatitude, with values of ∼1% achievable in mass and radius if both the inclination and colatitude are ≳60°.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Detection of z ~ 2.3 Cosmic Voids from 3D Ly α Forest Tomography in the COSMOS Field

We present the most distant detection of cosmic voids ( z ~ 2.3) and the first detection of three-dimensional voids in the Lyα forest. We used a 3D tomographic map of the absorption with an effective comoving spatial resolution of 2.5 h -1 Mpc and a volume of 3.15 × 10 5 h 3 Mpc 3 , which was reconstructed from moderate-resolution Keck I/LRIS spectra of 240 background Lyman-break galaxies and quasars in a 0.16 deg 2 footprint in the COSMOS field. Voids were detected using a spherical overdensity finder calibrated from hydrodynamical simulations of the intergalactic medium (IGM). This allows us to identify voids in the IGM corresponding to voids in the underlying matter density field, yielding a consistent volume fraction of voids in both data (19.5%) and simulations (18.2%). We fit excursion set models to the void radius function and compare the radially averaged stacked profiles of large voids ( r > 5 h -1 Mpc) to stacked voids in mock observations and the simulated density field. Comparing with 432 coeval galaxies with spectroscopic redshifts in the same volume as the tomographic map, we find that the tomography-identified voids are underdense in galaxies by 5.95σ compared to random cells.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

Effects of body shape on disk-integrated spectral reflectance

The effects of scattering geometry on the spectral reflectance curves of various powdered substances representing possible planetary surfaces are examined. Reflectance measurements were made as a function of incidence angle and emission angle for ground samples of a slightly oxidized basalt, the C3 carbonaceous chondrite Allende, and the ordinary chondrite Bruderheim, and compared with various photometric functions. The best fit to the observations is obtained with the photometric function of Goguen (1981) and Hapke (1981), and with a simpler function representing a linear combination of a Hapke-Irvine law and Lambert's law. The wavelength dependence of the adopted functions implies that the spherical reflectance of a flat sample will differ from that of a spherical planet of the same material, and that there will be a difference between the spectral reflectance curves of spherical and elliptical planets.

Gradie, J.↗
NASA NTRS1994

Retrieval of sizes and total masses of particles in volcanic clouds using AVHRR bands 4 and 5

The advanced very high resolution radiometer (AVHRR) sensor on polar orbiting NOAA satellites can discriminate between volcanic clouds and meteorological clouds using two-band data in the thermal infrared. This paper is aimed at developing a retrieval of the particle sizes, optical depth, and total masses of particles from AVHRR two-band data of volcanic clouds. Radiative transfer calculations are used with a semi-transparent cloud model that is based on assumptions of spherical particle shape, a homogeneous underlying surface, and a simple thin cloud parallel to the surface. The model is applied to observed AVHRR data from a 13-hour old drifting cloud from the August 19, 1992, eruption of Crater Peak/Spurr Volcano, Alaska. The AVHRR data fit in the range of results calculated by the model, which supports its credibility. According to the model results, the average of effective particle radius in the test frame of this cloud is in the range of 2 to 2.5 micrometer, the optical depth at 12 micrometer is about 0.60 - 0.65. The total estimated mass of ash in the air amounts to 0.24 - 0.31 x 10(exp 6) tons, which is about 0.7-0.9% of the mass measured in the ashfall blanket. Sensitivity tests show that the mass estimate is more sensitive to the assumed ash size distribution than it is to the ash composition.

Wen, Shiming↗
DOE OSTI2025

EMPDF : inferring the Milky Way mass with data-driven distribution function in phase space

We introduce the emPDF (empirical distribution function), a novel dynamical modelling method that infers the gravitational potential from kinematic tracers with optimal statistical efficiency under the minimal assumption of steady state. emPDF determines the best-fitting potential by maximizing the similarity between instantaneous kinematics and the time-averaged phase-space distribution function (DF), which is empirically constructed from observation upon the theoretical foundation of oPDF (Han et al. 2016). This approach eliminates the need for presumed functional forms of DFs or orbit libraries required by conventional DF- or orbit-based methods. emPDF stands out for its flexibility, efficiency, and capability in handling observational effects, making it preferable to the popular Jeans equation or other minimal assumption methods, especially for the Milky Way (MW) outer halo where tracers often have limited sample size and poor data quality. We apply emPDF to infer the MW mass profile using Gaia DR3 data of satellite galaxies and globular clusters, obtaining enclosed masses of M (,r) = 26±8, 46±8, 90±13⁠, and 149±40 x 10 10 M ⊙ at r = 30, 50, 100⁠, and 200 kpc, respectively. These are consistent with the updated constraints from simulation-informed DF fitting (Li et al. 2020). While the simulation-informed DF offers superior precision owing to the additional information extracted from simulations, emPDF is independent of such supplementary knowledge and applicable to general tracer populations. emPDF is currently implemented for tracers with complete 6D kinematics within spherical potentials, but it can potentially be extended to address more general problems.

Astrophysics of Galaxies (astro-ph.GA)↗
DOE OSTI2023

CHIPFoam Modification to Allow for Modelling SWIFT Foams

SWIFT (Silicon-water in Familiar Template) foams are silicone foams that have a continuous porosity without having any obvious spherical or ellipsoidal pore shapes. Like other foamed rubber, SWIFT foams can experience large reversible deformations. Because of their microstructure, SWIFT foams have a stiffering behavior in compression that is very abrupt compared to most other foamed rubbers. Work by Benedikt reported by Miller showed that fitting compressive test data with the CHIPFoam model in its current form is impractical and inaccurate. This report describes an extension made to CHIPFoam to allow the modeling of SWIFT foams. While the extension allows the modeling of SWIFT foams, it also extends the range of porosity of foams that can be modelled with CHIPFoam to include very high porosity (>99%) foamed rubber. The previous model was unstable for porosities greater than 70%. The changes made to CHIPFoam include the introduction of a critical porosity at which the compressible Danielsson strain energy function that couples volumetric and deviatoric mechanical behavior begins to increase with compression. This new feature also required the introduction of a third term in the multiplicative decomposition of the relative volume, J . The organization of this report is described here. The first section, SWIFT Foam Structure and Mechanical Behavior, describes the structure of SWIFT foams as contrasted to other foamed rubbers. It also contrasts the mechanical behavior of SWIFT foams to other more traditional foamed rubbers. The second section, CHIPFoam Modifications, describes the specific modifications done to be able to modify the delayed rapid compressive stiffening seen in SWIFT foams. The third section, Results, demonstrates the use of the modified model to fit compressive loading response data from a SWIFT foam sample. The fourth section, Conclusions and Further Work, puts the work in perspective and presents likely further work.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Multilevel Morphology of Complex Nanoporous Materials

This work exploits gas adsorption and small-angle X-ray scattering (SAXS) to determine the morphology of complex nanoporous materials. We resolve multiple classes of porosity including previously undetected large-scale texture that significantly compromises the canonical interpretation of gas adsorption. Specifically, a UVM-7 class mesoporous silica was synthesized that has morphological features on three length scales: macropores due to packing of 150 nm globules, 1.9 nm radius spherical mesopores inside the globules, and >7 nm pockets on and between the globules. The total and external surface areas, as well as the mesopore volume, were determined using gas adsorption (α s -plot) and SAXS. A new approach was applied to the SAXS data using multilevel fitting to determine the surface areas on multiple length scales. The SAXS analysis code is applicable to any two-phase system and is freely available to the public. The total surface area determined by SAXS was 12% greater than that obtained by gas adsorption. The macropore interfacial area, however, is only 30% of the external surface area determined by the α s -plot. The overestimation of the external surface area by the αs-plot method is attributed to capillary condensation in nanoscale surface irregularities. Here, the discrepancy is resolved assuming that the macropore–globule interfaces harbor fractally distributed nooks and crannies, which lead to gas adsorption at pressures above the mesopore filling pressure.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Modeling and interpretation of the optical and HST UV spectrum of SN 1993J

We present spherically symmetric, non-local thermodynamic equilibrium (non-LTE), expanding model atmosphere synthetic spectra of SN 1993J and compare them to the ultraviolet and optical spectra obtained simultaneously on 1993 April 15 by the Hubble Space Telescope (HST) and at Lick Observatory. We are able to fit the optical and HST near-UV spectra with a variety of compositions, but find that our best fits are obtained with an enhanced helium abundance (Y = 0.8) and significant nonthermal ionization due to gamma rays produced by the decay of Ni-56. We find the effects of X-rays produced by thermal bremsstrahlung in the circumstellar material to have little effect on the spectrum. Our current models are unable to fit the HST far-UV region, predicting too little flux and strong, unseen Fe II features.

Baron, E.↗