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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Development of an x-ray radiography platform to study laser-direct-drive energy coupling at the National Ignition Facility

A platform has been developed to study laser-direct-drive energy coupling at the National Ignition Facility (NIF) using a plastic sphere target irradiated in a polar-direct-drive geometry to launch a spherically converging shock wave. To diagnose this system evolution, eight NIF laser beams are directed onto a curved Cu foil to generate Heα line emission at a photon energy of 8.4 keV. These x rays are collected by a 100-ps gated x-ray imager in the opposing port to produce temporally gated radiographs. The platform is capable of acquiring images during and after the laser drive launches the shock wave. A backlighter profile is fit to the radiographs, and the resulting transmission images are Abel inverted to infer radial density profiles of the shock front and to track its temporal evolution. Here, the measurements provide experimental shock trajectories and radial density profiles that are compared to 2D radiation-hydrodynamic simulations using cross-beam energy transfer and nonlocal heat-transport models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Re-acceleration model for radio relics with spectral curvature

Most of the observed features of radio gischt relics, such as spectral steepening across the relic width and a power-law-like integrated spectrum, can be adequately explained by a diffusive shock acceleration (DSA) model in which relativistic electrons are (re-)accelerated at shock waves induced in the intracluster medium. However, the steep spectral curvature in the integrated spectrum above ∼2 GHz detected in some radio relics, such as the Sausage relic in cluster CIZA J2242.8+5301, may not be interpreted by the simple radiative cooling of postshock electrons. In order to understand such steepening, we consider here a model in which a spherical shock sweeps through and then exits out of a finite-size cloud with fossil relativistic electrons. The ensuing integrated radio spectrum is expected to steepen much more than predicted for aging postshock electrons, since the re-acceleration stops after the cloud-crossing time. Using DSA simulations that are intended to reproduce radio observations of the Sausage relic, we show that both the integrated radio spectrum and the surface brightness profile can be fitted reasonably well, if a shock of speed u{sub s} ∼ 2.5–2.8 × 10{sup 3} km s{sup −1} and a sonic Mach number M{sub s} ∼ 2.7–3.0 traverses a fossil cloud for ∼45 Myr, and the postshock electrons cool further for another ∼10 Myr. This attempt illustrates that steep curved spectra of some radio gischt relics could be modeled by adjusting the shape of the fossil electron spectrum and adopting the specific configuration of the fossil cloud.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Mars gravity derived from long-period motion of Mariner 9

The spherical harmonic coefficients in the gravity field of Mars are determined using a two stage filtering process. In the first stage, a weighted least squares orbit determination processor is applied to Mariner 9 Doppler tracking data to estimate Kepler elements and rates. In the second stage, the Kepler element rates are used as input to another weighted least squares processor that estimates gravity coefficients using the Lagrange perturbation equations. The novelty in this method is the fitting of element rates rather than the elements themselves. This procedure decouples the results from the effects of internal spacecraft dynamics.

Ferrari, A. J.↗
NASA NTRS1975

The Bulk Lunar Electrical Conductivity

The electrical conductivity structure was studied of a spherically layered moon consistent with the very low frequency magnetic data collected on the lunar surface and by Explorer 35. In order to obtain good agreement with the lunar surface magnetometer observations, the inclusion of a void cavity behind the moon requires a conductivity at shallow depths higher than that of models having the solar wind impinging on all sides. By varying only the source parameters, a conductivity model can be found that yields a good fit to both the tangential response upstream and the radial response downstream. This model also satisfies the dark side tangential response in the frequency range above 0.006 Hz, but the few data points presently available below this range do not seem to agree with the theory.

Leavy, Donald Lucien↗
DOE OSTI2015

SCATTER-DOMINATED INTERPLANETARY TRANSPORT OF SOLAR ENERGETIC PARTICLES IN LARGE GRADUAL EVENTS AND THE FORMATION OF DOUBLE POWER-LAW DIFFERENTIAL FLUENCE SPECTRA OF GROUND-LEVEL EVENTS DURING SOLAR CYCLE 23

The effects of scatter-dominated interplanetary transport on the spectral properties of the differential fluence of large gradual solar energetic particle (SEP) events are investigated analytically. The model assumes for simplicity radial constant solar wind and radial magnetic field. The radial diffusion coefficient is calculated with quasilinear theory by assuming a spectrum of Alfvén waves propagating parallel to the magnetic field. Cross-field transport is neglected. The model takes into consideration several essential features of gradual event transport: nearly isotropic ion distributions, adiabatic deceleration in a divergent solar wind, and particle radial scattering mean free paths increasing with energy. Assuming an impulsive and spherically symmetric injection of SEPs with a power-law spectrum near the Sun, the predicted differential fluence spectrum exhibits at 1 AU three distinctive power laws for different energy domains. The model naturally reproduces the spectral features of the double power-law proton differential fluence spectra that tend to be observed in extremely large SEP events. We select nine western ground-level events (GLEs) out of the 16 GLEs during Solar Cycle 23 and fit the observed double power-law spectra to the analytical predictions. The compression ratio of the accelerating shock wave, the power-law index of the ambient wave intensity, and the proton radial scattering mean free path are determined for the nine GLEs. The derived parameters are generally in agreement with the characteristic values expected for large gradual SEP events.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Time-dependent mapping of the magnetic field at the core-mantle boundary

We consider the problem of constructing a time-dependent map of the magnetic field at the core-mantle boundary. We use almost all the available data from the last 300 years to produce two maps, one for the period 1690-1840 and the other for 1840-1990. We represent the spatial dependency of the field using spherical harmonics, the time dependency using a cubic B-spline basis, and seek the smoothest solutions compatible with the observations. We argue that, for observations from permanent magnetic observatories, the most efficient strategy is to use the first differences of annual means; for satellite data, the most efficient strategy is simply to limit the number of data used so as to minimize any tendency to map the crustal field into the core field. The resulting model fits the observatory data better than any previous model. The resulting time-dependent field map exhibits much of the same structure in the field and its secular variation identified in earlier studies.

Bloxham, Jeremy↗
DOE OSTI2015

Size-dependent fluorescence of bioaerosols: Mathematical model using fluorescing and absorbing molecules in bacteria

This paper uses a mathematical model of fluorescent biological particles composed of bacteria and/or proteins (mostly as in Hill et al., 2013 [23]) to investigate the size-dependence of the total fluorescence emitted in all directions. The model applies to particles which have negligible reabsorption of fluorescence within the particle. The specific particles modeled here are composed of ovalbumin and of a generic Bacillus. The particles need not be spherical, and in some cases need not be homogeneous. However, the results calculated in this paper are for spherical homogeneous particles. Light absorbing and fluorescing molecules included in the model are amino acids, nucleic acids, and several coenzymes. Here the excitation wavelength is 266 nm. The emission range, 300 to 370 nm, encompasses the fluorescence of tryptophan. The fluorescence cross section (C F ) is calculated and compared with one set of published measured values. We investigate power law (Ad y ) approximations to C F , where d is diameter, and A and y are parameters adjusted to fit the data, and examine how y varies with d and composition, including the fraction as water. The particle's fluorescence efficiency (Q F =C F /geometric-cross-section) can be written for homogeneous particles as Q abs R F , where Q abs is the absorption efficiency, and R F , the fraction of the absorbed light emitted as fluorescence, is independent of size and shape. When Q F is plotted vs. m i d or mi(m r -1)d, where m=m r +im i is the complex refractive index, the plots for different fractions of water in the particle tend to overlap.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

A high temperature thermal diffusivity determination procedure for solids and liquids

A new method for measuring the thermal diffusivity of materials at high temperatures is presented. The method is applicable to solids on Earth, and to liquids in the reduced gravity environment of space. It is especially suited to levitated liquid metals at elevated temperatures where thermal diffusivity data is not available. The method is applied in two parts, such that lumped analysis is valid in the first part, and Fourier's law of conduction is valid in the second. In both parts, the spherical specimen is assumed to have been heated to a desired temperature and cooled. An inverse conduction problem is then formulated and solved using Laplace transformation techniques. Using this solution, sample sizes, and experimentally obtained surface temperature history, the thermal diffusivity is determined by minimizing a function that satisifies the heat balance at the surface. Minimization is performed using a modified quasilinearization algorithm. Accuracy is very sensitive to error in the temperature data and increases with better curve fits to the temperature data. An error analysis is also performed, and the effect of errors in the various parameters on the evaluated thermal diffusivity is determined. An experimental study for solids on Earth is suggested, before development for implementation in space.

Bayazitoglu, Yildiz↗
DOE OSTI2024

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Investigating fusion plasma instabilities in the Mega Amp Spherical Tokamak using mega electron volt proton emissions (invited)

The proton detector (PD) measures 3 MeV proton yield distributions from deuterium-deuterium fusion reactions within the Mega Amp Spherical Tokamak (MAST). The PD’s compact four-channel system of collimated and individually oriented silicon detectors probes different regions of the plasma, detecting protons (with gyro radii large enough to be unconfined) leaving the plasma on curved trajectories during neutral beam injection. From first PD data obtained during plasma operation in 2013, proton production rates (up to several hundred kHz and 1 ms time resolution) during sawtooth events were compared to the corresponding MAST neutron camera data. Fitted proton emission profiles in the poloidal plane demonstrate the capabilities of this new system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Classifying and analyzing small-angle scattering data using weighted k nearest neighbors machine learning techniques

A consistent challenge for both new and expert practitioners of small-angle scattering (SAS) lies in determining how to analyze the data, given the limited information content of said data and the large number of models that can be employed. Machine learning (ML) methods are powerful tools for classifying data that have found diverse applications in many fields of science. Here, ML methods are applied to the problem of classifying SAS data for the most appropriate model to use for data analysis. The approach employed is built around the method of weighted k nearest neighbors (wKNN), and utilizes a subset of the models implemented in the SasView package (https://www.sasview.org/) for generating a well defined set of training and testing data. The prediction rate of the wKNN method implemented here using a subset of SasView models is reasonably good for many of the models, but has difficulty with others, notably those based on spherical structures. A novel expansion of the wKNN method was also developed, which uses Gaussian processes to produce local surrogate models for the classification, and this significantly improves the classification accuracy. Further, by integrating a stochastic gradient descent method during post-processing, it is possible to leverage the local surrogate model both to classify the SAS data with high accuracy and to predict the structural parameters that best describe the data. The linking of data classification and model fitting has the potential to facilitate the translation of measured data into results for both novice and expert practitioners of SAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Main Bearing Replacement and Damage - A Field Data Study on 15 Gigawatts of Wind Energy Capacity

This study seeks to establish a comprehensive baseline of knowledge for the replacement and damage of main bearings in wind turbines. The purpose of this report is to provide a high-level summary of the data set, methodology, and results of this work. Full technical details and an extended analysis will be made available in a future publication. We collected data on main bearing replacements and reported damage from industrial partners based in Europe and the United States. In total, we obtained data for 167 wind power plants, with a combined capacity of 15.3 gigawatts (GW). Most of the data set was comprised of land-based, three-point mount, spherical roller bearings. Within this data set were 689 instances of main bearing replacement. Analysis was undertaken in two parts: first, a statistical analysis of the main bearing time-to-replacement data using survival analysis techniques; second, quantitative and qualitative analyses of the obtained damage information. Our results showed that 10% of a fixed main bearing population would be expected to have been replaced by 10.5 years. This is close to half of the 20-year design value. Fitted parametric distributions then indicated that by year 20, between 22% and 25% of main bearings are expected to have been replaced. Analysis of the damage reports revealed spalling to be the main type of damage listed. The additional presence of surface damage in the collected data indicates that at least part of the spalling cases are likely due to surface-initiated rolling contact fatigue. At this stage is not clear what proportion of spalling cases result from "wear induced", surface-initiated and subsurface-initiated rolling contact fatigue. While this work provides important insights into the current state of main bearing replacements and damage, many questions remain. An ongoing and expanding data collection and analysis effort focused on wind turbine main bearings is therefore recommended.

17 WIND ENERGY↗
NASA NTRS1994

Limits on coronal reflection using high-frequency solar oscillations

Acoustic waves in the Sun with frequencies above about 5.3 mHz can propagate in the chromosphere. We examine imaged solar intensity data for evidence of reflection of these waves in the upper chromosphere, where the temperature increases by a large factor over a short distance. Our method is to compare the observed and theoretically derived frequency spacings between peaks in the power spectrum. We find that our theoretical frequencies provide the best fit to the data when the reflection in the upper atmosphere is eliminated. In particular, the model of Kumar (1993b), which includes the source depth, and radiative damping, in the calculation of power spectra but ignores chromospheric reflection, gives peak frequencies that are in good agreement with the observations. For acoustic waves of frequency greater than 6 mHz we put an upper limit to the reflectivity of chromosphere and corona, using our method, of about 10%. At a given spherical harmonic degree, the frequency spacing between peaks in the data generally decreases with increasing frequency, because the lower turning point of the waves is moving inward. However, between 5 and 5.5 mHz the frequency spacing increases slightly. This feature is probably associated with the acoustic cutoff frequency in the solar atmosphere, i.e., it indicates a transition from trapped waves to propagating waves. We are able to reproduce the observed behavior by a crude modeling of the solar atmosphere. Further study of these peaks should provide an independent way of exploring the mean structure of the solar atmosphere, particularly around the temperature minimum region.

Kumar, P.↗
DOE OSTI2019

The Shape and Density Evolution of Snow Aggregates

Abstract Snow aggregates evolve into a variety of observed shapes and densities. Despite this diversity, models and observational studies employ fractal or Euclidean geometric measures that are assumed universal for all aggregates. This work therefore seeks to improve understanding and representation of snow aggregate geometry and its evolution by characterizing distributions of both observed and Monte Carlo–generated aggregates. Two separate datasets of best-fit ellipsoid estimates derived from Multi-Angle Snowflake Camera (MASC) observations suggest the use of a bivariate beta distribution model for capturing aggregate shapes. Product moments of this model capture shape effects to within 4% of observations. This mathematical model is used along with Monte Carlo simulated aggregates to study how combinations of monomer properties affect aggregate shape evolution. Plate aggregates of any aspect ratio produce a consistent ellipsoid shape evolution whereas thin column aggregates evolve to become more spherical. Thin column aggregates yield fractal dimensions much less than the often-assumed value of 2.0. Ellipsoid densities and fractal analogs of density (lacunarity) are much more variable depending on combinations of monomer size and shape. Simple mathematical scaling relationships can explain the persistent triaxial ellipsoid shapes that appear in both observed and modeled aggregates. Overall, both simulations and observations prove aggregates are rarely oblate. Therefore, the use of the proposed bivariate ellipsoid distribution in models will allow for similar-sized aggregates to exhibit a realistic dispersion of masses and fall speeds.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1975

Emission-line strengths and the chemical compositions of quasi-stellar objects

Estimates of emission-line strengths in 160 QSOs are listed. They are used, together with those for 60 objects published earlier, to compile a table of rest-frame equivalent widths for a composite QSO containing 28 lines seen in objects with z from 0.06 to 3.53. Models are computed for matter partially filling a spherical volume, ionized by ultraviolet flux from a central object and optically thick with absorption from heavy elements taken into account. Constant temperature and density are taken for the He II region, where most of the important lines are emitted. Allowance is made for lower and higher temperatures in the H I and He III regions. Emergent emission-line equivalent widths are computed and compared with the observations. An abnormally low helium abundance is not needed to achieve a fit, since the only anomalously weak line is He II 4686 A whose upper level may be depopulated by stimulated emission of He II 1215 A caused by very strong L-alpha 1216 A flux.

Chan, Y.-W. T.↗
NASA NTRS1979

Gravity field of Venus - A preliminary analysis

The gravitational field of Venus obtained by tracking the Pioneer Venus Orbiter is examined. For each spacecraft orbit, two hours of Doppler data centered around periapsis were used to estimate spacecraft position and velocity and the velocity residuals obtained were spline fit and differentiated to produce line of sight gravitational accelerations. Consistent variations in line of sight accelerations from orbit to orbit reveal the presence of gravitational anomalies. A simulation of isostatic compensation for an elevated region on the surface of Venus indicates that the mean depth of compensation is no greater than about 100 km. Gravitational spectra obtained from a Fourier analysis of line of sight accelerations from selected Venus orbits are compared to the earth's gravitational spectrum and spherical harmonic gravitational potential power spectra of the earth, the moon and Mars. The Venus power spectrum is found to be remarkably similar to that of the earth, however systematic variations in the harmonics suggest differences in dynamic processes or lithospheric behavior.

Phillips, R. J.↗
NASA NTRS2000

The Investigation of the Effects of Gravity on Single Bubble Sonoluminescence

In single bubble following it's rapid collapse each cycle of oscillation of an ultrasonic field. Since widely varying length and time scales affect the bubble dynamics and optical emission processes, it is difficult to anticipate the importance of the effects of gravity present for observations on earth. Our bubble is driven in an acoustically resonating cavity at it's first harmonic mode. The acoustical radiation pressure (Bjerknes force) will then keep it suspended in the center near the pressure antinode. When driven in a region where the diffusive processes balance the bubble it acts in a nonlinear but regular way, emitting a short (approx. 200ps) burst of light each acoustic cycle. Balancing the Bjerknes force with buoyancy, as in, we can see that the bubble should be displaced from the velocity node approximately 20m at normal gravity. Therefore, water flows past the bubble at the time of collapse. Gravitation also changes the ambient pressure at the bubble's location, as Delta.P = rho.g.h this gives a change of approximately -0.5% in our experiment when going from 1.8g to 0g. Studies of ambient pressure changes were also done in order to assess these effects. Inside a pressure sealed chamber a spherical glass cell is filled with distilled water which has been degassed to 120mmHg. A bubble is then trapped in the center and driven by a piezoelectric transducer at 32.2kHz attached to the side of the cell. An optical system is then set up to take strobbed video images along and light emission data simultaneously. Temperature, pressure, drive voltage, and listener voltage are also monitored. PMT output in Volts The radii of the bubbles for both experiment s are fit using the Rayleigh-Plesset equation and the acoustic drive amplitude and the ambient bubble radius are found. There is little change in the acoustic drive amplitude as we expect, since we are not varying the drive voltage. However. the ambient bubble radius goes up considerably. These changes (increased light output, increased maximum bubble radius, and increased ambient bubble radius) are also observed when the ambient pressure is varied in the laboratory by an amount similar to that due to gravitation. The changes in the ambient bubble radius and light output with a change in ambient pressure are predicted by the "dissociation hypothesis" and have been observed by other groups in the laboratory. It seems clear that buoyancy's effect on light output and bubble radius, are at best on the same order as the effects of ambient pressure.

Dzikowicz, Ben↗
NASA NTRS2004

Magnectic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynamical hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth s core and uncertainty. If this agrees with the seismologic value, then the hypotheses pass the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth s core, this yields a JGR-PI, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale change.In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at length- scale l/k, if one presumes kinetic energy is converted in either eddy- overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magneto- geostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core-mantle boundary. The resulting theoretical spectra for the core- source field and its SV are far more compatible with observation. The conversion time-scale of order 120 years is pseudo-scale-invariant. Magnetic spectra of other planets may differ; however, if a transition to non-conducting fluid hydrogen in Jupiter acts as barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte↗