Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 307 records · Page 17

Absolute wind velocities in the lower thermosphere of Venus using infrared heterodyne spectroscopy

NASA's IR Telescope Facility and the McMath Solar Telescope have yielded absolute wind velocities in the Venus thermosphere for December 1985 to March 1987 with sufficient spatial resolution for circulation model discrimination. A qualitative analysis of beam-integrated winds indicates subsolar-to-antisolar circulation in the lower thermosphere; horizontal wind velocity was derived from a two-parameter model wind field of subsolar-antisolar and zonal components. A unique model fit common to all observing periods possessed 120 m/sec subsolar-antisolar and 25 m/sec zonal retrograde components, consistent with the Bougher et al. (1986, 1988) hydrodynamical models for 110 km.

Goldstein, Jeffrey J.↗

Mesospheric Observations and Modeling of the Zeeman Split 233.9 GHz O-18O-16 line

Observations made from Kitt Peak, AZ, of the 233.9 GHz emission line of O-18O-16 in the upper stratosphere and lower mesosphere are reported. A good model fit to the line is obtained by incorporation of pressure and temperature broadening effects, as well as a Zeeman splitting algorithm that uses a standard geomagnetic field model and a paramagnetic Hamiltonian description of the molecular energy states. These observations are used, along with the well known O-18O-16 mixing ratio, to establish absolute calibration for observations of other chemical species from Kitt Peak. Repeated measurements show no change in this absolute calibration between observation dates. The wide magnetic splitting (+/-1.8 MHz) exhibited by this line with only six Zeeman components provides a unique test of middle atmosphere Zeeman effect model calculations, supporting the use Of 02 lines by microwave atmospheric sounders to measure pressure and temperature.

Sandor, Brad J.↗

Final Report

This grant was made to fund two accepted proposals to observe with the ASCA satellite: (1) Debunking the myth of two-temperature coronae of active stars, and (2) Dynamic coronae of dMe stars. We obtained the requested observations and have now completed two papers that will be published in the Astrophysical Journal. The RS CVn binary star UX Ari was observed for 14 hours with all four detectors onboard the ASCA satellite. The X-ray emission was at a constant, quiescent level during the first 12 hours, after which time a powerful flare with a peak luminosity of 1.4 x 10(exp 32) ergs/s started. The flare was observed until shortly after its peak. We present a spectral and temporal analysis of the UX Ari observations and analyze the data with a two-ribbon flare model including estimates for cooling losses. A time-dependent reconstruction of the emission measure (EM) distribution shows that two separate plasma components evolve during the flare (one being identified with the quiescent EM). Most of the flare EM reaches temperatures between 50 MK and 100 MK or more. Magnetic confinement requires the loop arcade to be geometrically large, with length scales on the order of one stellar radius. The electron densities inferred from the model decrease from initial values around 10(exp 12)/cc early in the flare to about 10(exp 11)/cc at the flare peak. The best-fit models require surface magnetic field strengths of a few hundred G, compatible with the maximum photospheric fields expected from equipartition. The flare parameters imply a (conductive and radiative) cooling loss time of less than one hour at flare peak. The elemental abundances increase significantly during the flare rise, with the abundances of the low-FIP elements Fe, Mg, Si, and Ni typically increasing to higher levels than the high-FIP elements such as S or Ne. The Fe abundance increases from (17 +/- 4)% of the solar photospheric value during quiescence up to (89 +/- 18)% at flare peak. A fractionation process that occurs during the chromospheric evaporation phase may selectively enrich low-FIP elements as in the solar corona; alternatively, the chromospheric evaporation may itself bring metal-rich plasma into the metal-poor corona.

Linsky, Jeffrey L.↗

Vertical Distributions of PH3 in Saturn from Observations of Its 1-0 and 3-2 Rotational Lines

Far-infrared Fourier-transform spectrometer measurements of the 1-0 and 3-2 PH3 transitions in Saturn's disk near 267 and 800 GHz (8.9 and 26.7/cm), respectively, were analyzed simultaneously to derive a global mean profile for the PH3 vertical mixing ratio between 100 and 600 mbar total pressure. The far-infrared spectrum is relatively free from spectral interlopers, suffers minimal absorption or scattering by atmospheric particulates, and contains intrinsically weak PH3 lines that are sensitive to a range of atmospheric depths. The combined spectra are inconsistent with a uniform tropospheric mixing ratio, even with a stratospheric cut-off. They are consistent with a volume mixing ratio of PH3 that drops from 1.2 x 10(exp -5) at 645 mbar pressure to a value of 4.1 x 10(exp -7) at 150 mbar pressure, a decrease that is linear is log abundance vs log pressure. The mixing ratio could drop even more quickly at atmospheric pressures below 150 mbar and still be consistent with the data. The mixing ratio may well remain constant with depth for pressures above 630 mbar. The maximum PH3 mixing ratio in this model is consistent with a [P]/[H] ratio in the deep atmosphere that is about a factor of 10 higher than solar composition. Such a model is consistent with rapid mixing up to the radiative-convective boundary and transport by, for example, vertical waves just above this boundary. In the best fitting model, the eddy diffusion coefficient is approximately 10(exp 4) sq cm near 630 mbar, and it must increase with altitude. The predominant PH3 loss mechanisms are direct photolysis by UV radiation and scavenging by H atoms produced by the photolysis.

Orton, G. S.↗

String Fragmentation Model in Space Radiation Problems

String fragmentation models such as the Lund Model fit experimental particle production cross sections very well in the high-energy limit. This paper gives an introduction of the massless relativistic string in the Lund Model and shows how it can be modified with a simple assumption to produce formulas for meson production cross sections for space radiation research. The results of the string model are compared with inclusive pion production data from proton-proton collision experiments.

Tang, Alfred↗

Multi-Epoch XMM Observations of NGC4258

The goal of this project was detection of variability in the X-ray absorption column of the AGN in NGC4258 through monitoring with the XMM satellite. We have accomplished this goal and submitted the results to ApJ for publication in a paper entitled, "X-ray Luminosity and Absorption Column Fluctuations in the H2O Maser Galaxy NGC4258 from Weeks to Years," by Fruscione, A., Greenhill, L.J., Filippedco, A.V., Moran, J.M., Hermstein, J.R., and Galle, E. We have received a favorable referee report and expect the article will appear in 2005. To complete the project, we reduced our four epochs of XMM data for NGC4258, one archival XMM observation, and all existing Chandra datasets for NGC4258 (with the latest calibration tables and a grid of corrections for pileup). Self-consistent reduction of all these data permitted detailed comparison that could not have been accomplished simply by taking published model fits that appear in the literature. To accumulate a broader monitoring record, we combined the Chandra and XMM results with those published for SAX and ASCA. We modeled the Chandra and XMM data self-consistently with partially absorbed, hard power-law, soft thermal plasmas, and soft power-law components. Over nine years, the photo-electric absorbing column exhibited a 40% drop between two ASCA epochs separated by 3 years and a 60% rise between two XMM epochs separated by just 5 months. In contract, uncorrelated factor of of 2-3 changes were seen in absorbed flux on te timescale of years, which suggests intrinsic variability of the central engine. The warped disk that is a known source of H2O maser emission in 4258 is believed to cross the line of sight to the central engine. We have proposed that the variations in absorbing column arise from inhomogeneities in the rotating disk, as they sweep across the line of sight. We estimate from the XMM data that the inhomogeneities are about 1E+15 cm in size at radii greater than 0.27 pc. This is consistent with the estimated radius of the disk crossing estimated entirely independently, and it solidifies evidence that the warped accretion disk is the absorber in this (and possible other) AGN.

Greenhill, Lincoln J.↗

Non-Gaussian Analysis of Turbulent Boundary Layer Fluctuating Pressure on Aircraft Skin Panels

The purpose of the study is to investigate the probability density function (PDF) of turbulent boundary layer fluctuating pressures measured on the outer sidewall of a supersonic transport aircraft and to approximate these PDFs by analytical models. Experimental flight results show that the fluctuating pressure PDFs differ from the Gaussian distribution even for standard smooth surface conditions. The PDF tails are wider and longer than those of the Gaussian model. For pressure fluctuations in front of forward-facing step discontinuities, deviations from the Gaussian model are more significant and the PDFs become asymmetrical. There is a certain spatial pattern of the skewness and kurtosis behavior depending on the distance upstream from the step. All characteristics related to non-Gaussian behavior are highly dependent upon the distance from the step and the step height, less dependent on aircraft speed, and not dependent on the fuselage location. A Hermite polynomial transform model and a piecewise-Gaussian model fit the flight data well both for the smooth and stepped conditions. The piecewise-Gaussian approximation can be additionally regarded for convenience in usage after the model is constructed.

Rizzi, Stephen A.↗

Accounting for the Influence of Decorrelation in Microphone Phased Array Deconvolution Methods

Microphone phased arrays are a common tool for use in aeroacoustic wind tunnel testing. The analysis of acquired array data is known to suffer from decorrelation effects, where the coherence of an acoustic wave measured by a pair of microphones is degraded as the wave passes through a turbulent free shear layer or boundary layer. This paper describes, in detail, how to mitigate the influence of decorrelation effects when processing array data with deconvolution methods. This is done using the DAMAS algorithm as an example, applied to recent airframe noise test data acquired in the NASA Langley 14- by 22-Foot Subsonic Tunnel. Two ways of handling the turbulent propagation modeling, both assuming plane wave propagation, are described. Results show that while turbulence model fit parameters may differ, both methods output extremely similar deconvolution results. Further improvements likely require more accurate mean shear layer data prior to developing more involved turbulence models.

microphone array↗

Statistical estimates of the binary properties of rotational variables

ABSTRACT We present a model to estimate the average primary masses, companion mass ranges, the inclination limit for recognizing a rotational variable, and the primary mass spreads for populations of binary stars. The model fits a population’s binary mass function distribution and allows for a probability that some mass functions are incorrectly estimated. Using tests with synthetic data, we assess the model’s sensitivity to each parameter, finding that we are most sensitive to the average primary mass and the minimum companion mass, with less sensitivity to the inclination limit and little to no sensitivity to the primary mass spread. We apply the model to five populations of binary spotted rotational variables identified in ASAS-SN, computing their binary mass functions using RV data from APOGEE. Their average primary mass estimates are consistent with our expectations based on their CMD locations ($\sim 0.75 \, {\rm M}_{\odot }$ for lower main sequence primaries and $\sim 0.9$–$1.2 \, {\rm M}_{\odot }$ for RS CVn and sub-subgiants). Their companion mass range estimates allow companion masses down to $M_2/M_1\simeq 0.1$, although the main sequence population may have a higher minimum mass fraction ($\sim 0.4$). We see weak evidence of an inclination limit $\gtrsim 50^{\circ }$ for the main sequence and sub-subgiant groups and no evidence of an inclination limit in the other groups. No groups show strong evidence for a preferred primary mass spread. We conclude by demonstrating that the approach will provide significantly better estimates of the primary mass and the minimum mass ratio and reasonable sensitivity to the inclination limit with 10 times as many systems.

Phillips, Anya (ORCID:000900051914974X)↗

Dynamics-based halo model for large scale structure

Accurate modeling of the one-to-two halo transition has long been difficult to achieve. Here, we demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to (1) particles orbiting the halo and (2) particles infalling into the halo for the first time. Motivated by this decomposition, García et al. [Mon. Not. R. Astron. Soc. 521, 2464 (2023)] advocated for defining halos as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo-mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterized by a single physical scale 𝑟 h , which we refer to as the halo radius. The two-halo term is described using a simple yet accurate empirical model based on the Zel’dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo-mass correlation function is accurate at the ≈ 2% level for 𝑟∈ [0.1, 50] ℎ −1 Mpc. The Fourier transform of our best-fit model describes the halo-mass power spectrum with comparable accuracy for 𝑘 ∈ [0.06, 6.0] ℎ Mpc −1 .

79 ASTRONOMY AND ASTROPHYSICS↗

The center-to-limb

Center-to-limb measurements of the Ca I 6573 intercombination line and the Ca II 7324 forbidden line are compared with synthetic profiles based on a simple representation of the non-LTE Ca-Ca(+) ionization equilibrium. The effects of photoionizations from low lying excited states of neutral calcium are found to reduce the sensitivity of the 6573 center-to-limb behavior as a thermal structure diagnostic. The synthetic center-to-limb behavior is also sensitive to uncertainties in the nonthermal broadening. Nevertheless, the measured center-to-limb behavior of 6573 favors a 'cool' photospheric model over hotter models based on the Ca II K wings. The non-LTE calcium abundance obtained from the disk center equivalent widths of 6573 and 7324 using the best fit model is A(Ca) = 2.1 plus or minus 0.2 x 10 to the -6th (by number relative to hydrogen).

Source record↗

Zonal harmonic model of Saturn's magnetic field from Voyager 1 and 2 observations

An analysis of the magnetic field of Saturn is presented which takes into account both the Voyager 1 and 2 vector magnetic field observations. The analysis is based on the traditional spherical harmonic expansion of a scale potential to derive the magnetic field within 8 Saturn radii. A third-order zonal harmonic model fitted to Voyager 1 and 2 observations is found to be capable of predicting the magnetic field characteristics at one encounter based on those observed at another, unlike models including dipole and quadrupole terms only. The third-order model is noted to lead to significantly enhanced polar surface field intensities with respect to dipole models, and probably represents the axisymmetric part of a complex dynamo field.

Connerney, J. E. P.↗

'Alfven wing' models of the induced electrical current system at Io - A probe of the ionosphere of Io

A novel class of models of the Alfven wing interaction between the Io plasma torus and Io's ionosphere has been developed and used to compute the magnetic field fluctuations at the position of the Voyager 1 flyby. These computations of the magnetic field signature are compared with the magnetic field measurements made during the flyby, and the differences between the model and observations are used as a probe of the structure of Io's ionosphere. The results of this model fitting indicate that a significant atmospheric neutral column density is required over a major part of the trailing hemisphere. This can be consistent with the cold-trap model for Io's neutral atmoshpere only if the major plumes near the center of the trailing hemisphere provide significant volumetric coverage. The Io plasma torus charged-mass density required for a proper match to the data is larger than earlier estimates based on the width of the Alfven wing signature.

Herbert, F.↗

Magnetic Field-Line Lengths in Interplanetary Coronal Mass Ejections Inferred from Energetic Electron Events

About one quarter of the observed interplanetary coronal mass ejections (ICMEs) are characterized by enhanced magnetic fields that smoothly rotate in direction over timescales of about 10-50 hr. These ICMEs have the appearance of magnetic flux ropes and are known as "magnetic clouds" (MCs). The total lengths of MC field lines can be determined using solar energetic particles of known speeds when the solar release times and the I AU onset times of the particles are known. A recent examination of about 30 near-relativistic (NR) electron events in and near 8 MCs showed no obvious indication that the field-line lengths were longest near the MC boundaries and shortest at the MC axes or outside the MCs, contrary to the expectations for a flux rope. Here we use the impulsive beamed NR electron events observed with the Electron Proton and Alpha Monitor instrument on the Advanced Composition Explorer spacecraft and type III radio bursts observed on the Wind spacecraft to determine the field-line lengths inside ICMEs included in the catalog of Richardson & Cane. In particular, we extend this technique to ICMEs that are not MCs and compare the field-line lengths inside MCs and non-MC ICMEs with those in the ambient solar wind outside the ICMEs. No significant differences of field-line lengths are found among MCs, ICMEs, and the ambient solar wind. The estimated number of ICME field-line turns is generally smaller than those deduced for flux-rope model fits to MCs. We also find cases in which the electron injections occur in solar active regions CARs) distant from the source ARs of the ICMEs, supporting CME models that require extensive coronal magnetic reconnection with surrounding fields. The field-line lengths are found to be statistically longer for the NR electron events classified as ramps and interpreted as shock injections somewhat delayed from the type III bursts. The path lengths of the remaining spike and pulse electron events are compared with model calculations of solar wind field-line lengths resulting from turbulence and found to be in good agreement.

Kahler, S. W.↗

XRISM Spectroscopy of the Fe Kα Emission Line in the Seyfert Active Galactic Nucleus NGC 4151 Reveals the Disk, Broad-line Region, and Torus

Abstract We present an analysis of the first two XRISM/Resolve spectra of the well-known Seyfert-1.5 active galactic nucleus (AGN) in NGC 4151, obtained in 2023 December. Our work focuses on the nature of the narrow Fe K α emission line at 6.4 keV, the strongest and most common X-ray line observed in AGN. The total line is found to consist of three components. Even the narrowest component of the line is resolved with evident Fe K α ,1 (6.404 keV) and K α ,2 (6.391 keV) contributions in a 2:1 flux ratio, fully consistent with neutral gas with negligible bulk velocity. Subject to the limitations of our models, the narrowest and intermediate-width components are consistent with emission from optically thin gas, suggesting that they arise in a disk atmosphere and/or wind. Modeling the three line components in terms of Keplerian broadening, they are readily associated with (1) the inner wall of the “torus,” (2) the innermost optical “broad-line region” (or “X-ray BLR”), and (3) a region with a radius of r ≃ 100 GM / c 2 that may signal a warp in the accretion disk. Viable alternative explanations of the broadest component include a fast-wind component and/or scattering; however, we find evidence of variability in the narrow Fe K α line complex on timescales consistent with small radii. The best-fit models are statistically superior to simple Voigt functions, but when fit with Voigt profiles the time-averaged lines are consistent with a projected velocity broadening of FWHM = 1600 − 200 + 400 km s − 1 . Overall, the resolution and sensitivity of XRISM show that the narrow Fe K line in AGN is an effective probe of all key parts of the accretion flow, as it is currently understood. We discuss the implications of these findings for our understanding of AGN accretion, future studies with XRISM, and X-ray-based black hole mass measurements.

Audard, Marc (ORCID:000000034721034X)↗

Low energy neutron light output characterization of EJ301D and deuterated stilbene with a comparison of light output characterization methods

The neutron-induced light yield of a 2.54 cm diameter by 2.54 cm long right circular cylinder of EJ301D and a (5.08 cm)3 custom made cube of deuterated trans-stilbene-d12 (d-stilbene) were measured over incident neutron energies from 300 keV to 2.2 MeV and 200 keV to 2.4 MeV, respectively. The measurements were performed using a time-of-flight experiment with a Cf-252 source and an approximately 1.5 m flight path. We compare three light output spectrum full energy deposition edge estimation methods: (1) simulating the neutron energy spectrum edge and fitting it to the light output spectrum, (2) using the inflection point of the light output spectrum edge (derivative method, a.k.a. Kornilov’s method), and (3) using an empirical model fit to the edge of the light output spectrum. Both the derivative and equation fit methods do not account for physical processes such as multiple neutron scattering in the detectors. They instead rely on assumptions about the linear shape continuum shape of the light output spectrum and the direct correlation between the location of the spectrum’s inflection point and maximum energy deposition. These assumptions were found to introduce bias into those methods when tested against simulated spectra with known edge locations. When tested against measured spectra the derivative method was found to differ from the simulation fit by greater than 30% at low energies with large discontinuities for adjacent data points above 800 keV incident neutron energy. The empirical equation fitting method was found to also exhibit bias of a similar magnitude, but with significantly more continuous behavior, especially with the lower count data of the smaller volumed EJ301D scintillator. Experimental light output yield for this neutron energy range is reported using the simulated spectrum fitting method because it includes physics neglected by the other methods, and did not exhibit the bias observed in the other methods

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multiple emission components in the Cygnus cocoon detected from Fermi -LAT observations

Star-forming regions may play an important role in the life cycle of Galactic cosmic rays (CRs), notably as home to specific acceleration mechanisms and transport conditions. Gamma-ray observations of Cygnus X have revealed the presence of an excess of hard-spectrum gamma-ray emission, possibly related to a cocoon of freshly accelerated particles. We seek an improved description of the gamma-ray emission from the cocoon using ~13 yr of observations with the Fermi-Large Area Telescope (LAT) and use it to further constrain the processes and objects responsible for the young CR population. We developed an emission model for a large region of interest, including a description of interstellar emission from the background population of CRs and recent models for other gamma-ray sources in the field. Thus, we performed an improved spectro-morphological characterisation of the residual emission including the cocoon. The best-fit model for the cocoon includes two main emission components: an extended component FCES G78.74+1.56, described by a 2D Gaussian of extension r 68 = 4.4° ± 0.1° -0.1° +0.1° and a smooth broken power law spectrum with spectral indices 1.67 ± 0.05 -0.01 +0.02 and 2.12 ± 0.02 -0.01 +0.00 below and above 3.0 ± 0.6 -0.2 +0.0 GeV, respectively; and a central component FCES G80.00+0.50, traced by the distribution of ionised gas within the borders of the photo-dissociation regions and with a power law spectrum of index 2.19 ± 0.03 -0.01 +0.00 that is significantly different from the spectrum of FCES G78.74+1.56. An additional extended emission component FCES G78.83+3.57, located on the edge of the central cavities in Cygnus X and with a spectrum compatible with that of FCES G80.00+0.50, is likely related to the cocoon. For the two brightest components FCES G80.00+0.50 and FCES G78.74+1.56, spectra and radial-azimuthal profiles of the emission can be accounted for in a diffusion-loss framework involving one single population of non-thermal particles with a flat injection spectrum. Particles span the full extent of FCES G78.74+1.56 as a result of diffusion from a central source, and give rise to source FCES G80.00+0.50 by interacting with ionised gas in the innermost region. For this simple diffusion-loss model, viable setups can be very different in terms of energetics, transport conditions, and timescales involved, and both hadronic and leptonic scenarios are possible. The solutions range from long-lasting particle acceleration, possibly in prominent star clusters such as Cyg OB2 and NGC 6910, to a more recent and short-lived release of particles within the last 10–100 kyr, likely from a supernova remnant. The observables extracted from our analysis can be used to perform detailed comparisons with advanced models of particle acceleration and transport in star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗

Empirical model of global thermospheric temperature and composition based on data from the OGO-6 quadrupole mass spectrometer

An empirical global model for magnetically quiet conditions has been derived from longitudinally averaged N2, O, and He densities by means of an expansion in spherical harmonics. The data were obtained by the OGO-6 neutral mass spectrometer and cover the altitude range 400 to 600 km for the period 27 June 1969 to 13 May 1971. The accuracy of the analytical description is of the order of the experimental error for He and O and about three times experimental error for N2, thus providing a reasonable overall representation of the satellite observations. Two model schemes are used: one representing densities extrapolated to 450 km and one representing densities extrapolated to 120 km with exospheric temperatures inferred from N2 densities. Using the best fit model parameters the global thermospheric structure is presented in the form of a number of contour plots.

Hedin, A. E.↗