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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1955

Effect of surface roughness on characteristics of spherical shock waves

An investigation has been conducted on a small-scale test layout in which direct observation of the shock wave movement with time could be made in order to determine the effects of surface roughness on the characteristics of spherical shock waves. Data were obtained with 15-gram pentolite charges at four heights of burst, both for a smooth surface and for a surface completely covered with pyramid-shaped roughness elements. The observations resulted in determinations of shock peak overpressure and Mach stem height as a function of distance for each test. Comparison of the smooth-surface data with those obtained for the extremely rough condition showed a small net effort of roughness on the shock peak overpressures at the surface for all burst heights, the effect being to lower the overpressures. The effect of surface roughness on the Mach stem formation and growth was to delay the formation at the greatest charge height and to lower the height of the Mach stem for all heights.Comparison of the free-air shock peak overpressures with larger scale data showed good similarity of the overpressure-distance relationships. The data did not fit a geometrical similarity parameter for the path of the triple point at different heights of burst suggested by other investigators. A simple similarity parameter (relating the horizontal distance to the theoretical point of Mach formation) was found which showed only a small influence of burst height on the path of the triple point. While the data presented provide knowledge of the effect of many surface-roughness elements on the overall shock characteristics, the data do not provide insight into the details of the air-flow characteristics along the surface, nor the relative contribution of individual roughness elements to the results obtained.

Huber, Paul W↗
DOE OSTI2016

The role of Dark Matter sub-halos in the non-thermal emission of galaxy clusters

Annihilation of Dark Matter (DM) particles has been recognized as one of the possible mechanisms for the production of non-thermal particles and radiation in galaxy clusters. Previous studies have shown that, while DM models can reproduce the spectral properties of the radio halo in the Coma cluster, they fail in reproducing the shape of the radio halo surface brightness because they produce a shape that is too concentrated towards the center of the cluster with respect to the observed one. However, in previous studies the DM distribution was modeled as a single spherically symmetric halo, while the DM distribution in Coma is found to have a complex and elongated shape. In this work we calculate a range of non-thermal emissions in the Coma cluster by using the observed distribution of DM sub-halos. We find that, by including the observed sub-halos in the DM model, we obtain a radio surface brightness with a shape similar to the observed one, and that the sub-halos boost the radio emission by a factor between 5 and 20%, thus allowing to reduce the gap between the annihilation cross section required to reproduce the radio halo flux and the upper limits derived from other observations, and that this gap can be explained by realistic values of the boosting factor due to smaller substructures. Models with neutralino mass of 9 GeV and composition τ{sup +} τ{sup −}, and mass of 43 GeV and composition b b-bar can fit the radio halo spectrum using the observed properties of the magnetic field in Coma, and do not predict a gamma-ray emission in excess compared to the recent Fermi-LAT upper limits. These findings make these DM models viable candidate to explain the origin of radio halos in galaxy clusters, avoiding the problems connected to the excessive gamma-ray emission expected from proton acceleration in most of the currently proposed models, where the acceleration of particles is directly or indirectly connected to events related to clusters merging. Therefore, DM models deserve to be better studied both from the theoretical and observational sides; the best spectral bands where it is possible to obtain better information are the radio and the gamma-ray bands, while we do not expect a strong emission in the X-ray band.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Drag of a heated sphere at low Reynolds numbers in the absence of buoyancy

Fully resolved simulations are used to quantify the effects of heat transfer in the absence of buoyancy on the drag of a spatially fixed heated spherical particle at low Reynolds numbers ( $Re$ ) in the range $10^{-3}\leqslant Re\leqslant 10$ in a variable-property fluid. The case where buoyancy is present is analysed in a subsequent paper. This analysis is carried out without making any assumptions on the amount of heat addition from the sphere and thus encompasses both the heating regime where the Boussinesq approximation holds and the regime where it breaks down. The particle is assumed to have a low Biot number, which means that the particle is uniformly at the same temperature and has no internal temperature gradients. Large deviations in the value of the drag coefficient as the temperature of the sphere increases are observed. When $Re , these deviations are explained using a low-Mach-number perturbation analysis as irrotational corrections to a Stokes–Oseen base flow. Correlations for the drag and Nusselt number of a heated sphere are proposed for the range of Reynolds numbers $10^{-3}\leqslant Re\leqslant 10$ which fit the computationally obtained values with less than 1 % and 3 % errors, respectively. These correlations can be used in simulations of gas–solid flows where the accuracy of the drag law affects the prediction of the overall flow behaviour. Finally, an analogy to incompressible flow over a modified sphere is demonstrated.

Mechanics↗
NASA NTRS2010

Lunar Gravity Field Determination Using SELENE Same-Beam Differential VLBI Tracking Data

A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM 100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from I-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080 +/- 0.0009 cu km/sq s (10 sigma). The lunar degree 2 potential Love number k2 was also estimated, and has a value of 0.0255 +/- 0.0016 (10 sigma as well).

Goossens, S.↗
DOE OSTI2024

From large to small $$ \mathcal{N} $$ = (4, 4) superconformal surface defects in holographic 6d SCFTs

Abstract Two-dimensional (2d)$$ \mathcal{N} $$ N = (4, 4) Lie superalgebras can be either “small” or “large”, meaning their R-symmetry is either$$ \mathfrak{so} $$ so (4) or$$ \mathfrak{so} $$ so (4) ⊕$$ \mathfrak{so} $$ so (4), respectively. Both cases admit a superconformal extension and fit into the one-parameter family$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ), with parameterγ∈ (−∞,∞). The large algebra corresponds to generic values ofγ, while the small case corresponds to a degeneration limit withγ→ −∞. In 11d supergravity, we study known solutions with superisometry algebra$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ) that are asymptotically locally AdS 7 ×𝕊 4 . These solutions are holographically dual to the 6d maximally superconformal field theory with 2d superconformal defects invariant under$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ). We show that a limit of these solutions, in whichγ→ −∞, reproduces another known class of solutions, holographically dual tosmall$$ \mathcal{N} $$ N = (4, 4) superconformal defects. We then use this limit to generate new small$$ \mathcal{N} $$ N = (4, 4) solutions with finite Ricci scalar, in contrast to the known small$$ \mathcal{N} $$ N = (4, 4) solutions. We then use holography to compute the entanglement entropy of a spherical region centered on these small$$ \mathcal{N} $$ N = (4, 4) defects, which provides a linear combination of defect Weyl anomaly coefficients that characterizes the number of defect-localized degrees of freedom. We also comment on the generalization of our results to include$$ \mathcal{N} $$ N = (0,4) surface defects through orbifolding.

Physics↗
DOE OSTI2016

High-energy tail of the Galactic Center gamma-ray excess

Observations by the Fermi-LAT have uncovered a bright, spherically symmetric excess surrounding the center of the Milky Way Galaxy. The spectrum of the ${\gamma}$-ray excess peaks sharply at an energy ${\sim}2\text{ }\text{ }\mathrm{GeV}$, exhibiting a hard spectrum at lower energies, and falls off quickly above an energy ${\sim}5\text{ }\text{ }\mathrm{GeV}$. The spectrum of the excess above ${\sim}10\text{ }\text{ }\mathrm{GeV}$ is potentially an important discriminator between different physical models for its origin. We focus on observations of the ${\gamma}$-ray excess at energies exceeding 10 GeV, finding: (1) a statistically significant excess remains in the energy range 9.5–47.5 GeV, which is not degenerate with known diffuse emission templates such as the Fermi bubbles, (2) the radial profile of the excess at high energies remains relatively consistent with data near the spectral peak (3) the data above ${\sim}5\text{ }\text{ }\mathrm{GeV}$ prefer a slightly greater ellipticity with a major axis oriented perpendicular to the Galactic plane. Using the recently developed non-Poissonian template fit, we find mild evidence for a point-source origin for the high-energy excess, although given the statistical and systematic uncertainties we show that a smooth origin of the high-energy emission cannot be ruled out. We discuss the implication of these findings for pulsar and dark matter models of the ${\gamma}$-ray excess. Finally we provide a number of updated measurements of the ${\gamma}$-ray excess, utilizing novel diffuse templates and the Pass 8 data set.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

X-ray spectrum and variability of 3C 120

The X-ray luminosity of 3C 120 is variable by a factor of 2.5 on time scales of days to months. The spectral slope of the 2-10 keV X-rays changes systematically in the sense that the higher intensity states are steeper. Simple power laws fitted to the data are consistent with an energy index alpha of 0.51 + or - 0.13 at minimum intensity and 0.83 + or - 0.03 at maximum, with a continuous variation in between. A more complicated spectral shape cannot be ruled out, but changes in absorbing column density alone cannot explain the data. The results are interpreted in terms of beamed synchrotron emission from a relativistic jet, a model which has been invoked previously for the spectra of BL Lac objects. The spectral variability is in marked contrast to the behavior of 3C 273, another superluminal source in which the X-ray spectral index is constant to a high degree of accuracy despite variability in time scales and amplitudes which are quite similar to those of 3C 120. The difference may be due to the degree of beaming of the respective core components. Another model which is consistent with the data involves thermal emission from isothermal shocks in a hot (approximately 10 to the 8th K) medium. Thermal Compton and spherical accretion X-ray source models do not reproduce the correct behavior of alpha versus flux.

Halpern, J. P.↗
NASA NTRS1978

Graphite-fiber-reinforced polyimide liners of various compositions in plain spherical bearings

A plain spherical bearing design with a ball diameter of 28.6 mm, a race length of 12.7 mm, and a 1.7-mm-thick, molded composite liner was evaluated. The liner material is a self-lubricating composite of graphite-fiber-reinforced polyimide resin (GFRPI). The liner is prepared by transfer molding a mixture of one part chopped graphite fiber and one part partially polymerized resin into the space between the bearing ball and the outer race and then completing the polymerization under heat and pressure. Several liner compositions were evaluated: two types of polyimide, condensation and addition; two types of graphite fiber, low and high modulus; and four powder additives - cadmium oxide, cadmium iodide, graphite fluoride, and molybdenum disulfide. The bearings were oscillated + or - 15 deg at 1 Hz for 20 kilocycles under a radial unit load of 29 MN sq m (4200 psi) in dry air at 25, 200, or 315 C. Both types of fiber and polyimide gave low friction and wear. A simple equation was developed to fit the wear-time data and adequately predicted wear to 100 kilocycles.

Sliney, H. E.↗
DOE OSTI2016

Solvation structure of the halides from x-ray absorption spectroscopy

Three-dimensional models for the aqueous solvation structures of chloride, bromide, and iodide are reported. K-edge extended X-ray absorption fine structure (EXAFS) and Minuit X-ray absorption near edge (MXAN) analyses found well-defined single shell solvation spheres for bromide and iodide. However, dissolved chloride proved structurally distinct, with two solvation shells needed to explain its strikingly different X-ray absorption near edge structure (XANES) spectrum. Final solvation models were as follows: iodide, 8 water molecules at 3.60 ± 0.13 Å and bromide, 8 water molecules at 3.40 ± 0.14 Å, while chloride solvation included 7 water molecules at 3.15 ± 0.10 Å, and a second shell of 7 water molecules at 4.14 ± 0.30 Å. Each of the three derived solvation shells is approximately uniformly disposed about the halides, with no global asymmetry. Time-dependent density functional theory calculations simulating the chloride XANES spectra following from alternative solvation spheres revealed surprising sensitivity of the electronic state to 6-, 7-, or 8-coordination, implying a strongly bounded phase space for the correct structure during an MXAN fit. MXAN analysis further showed that the asymmetric solvation predicted from molecular dynamics simulations using halide polarization can play no significant part in bulk solvation. Classical molecular dynamics used to explore chloride solvation found a 7-water solvation shell at 3.12 (−0.04/+0.3) Å, supporting the experimental result. These experiments provide the first fully three-dimensional structures presenting to atomic resolution the aqueous solvation spheres of the larger halide ions.

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

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

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

Neutron-unbound states in 31 Ne

Background: The Island of Inversion near the N = 20 shell gap is home to nuclei with reordered single-particle energy levels compared to the spherical shell model. Studies of 31 Ne have revealed that its ground state has a halo component, characterized by a valence neutron orbiting a deformed 30 Ne core. This lightly-bound nucleus with a separation energy of only S n = 170 keV is expected to have excited states that are neutron unbound. Purpose: The purpose of this experiment was to investigate the low-lying excited states in 31 Ne that decay by the emission of a single neutron. Methods: An 89 MeV/u 33 Mg beam impinged on a segmented Be reaction target. Neutron unbound states in 31 Ne were populated via a two-proton knockout reaction. The 30 Ne fragment and associated neutron from the decay of 31 Ne were detected by the MoNA-LISA-Sweeper experimental setup at the National Superconducting Cyclotron Laboratory. Invariant mass spectroscopy was used to reconstruct the two-body decay energy ( 30 Ne + n). Results: The two-body decay energy spectrum exhibits two features: a low-lying peak at 0.30 ± 0.17 MeV and a broad enhancement at 1.50 ± 0.33 MeV, each fit with an energy-dependent asymmetric Breit-Wigner line shape representing a resonance in the continuum. Accompanying shell model calculations using the FSU interaction within NuShellX, combined with cross-section cal culations using the eikonal reaction theory, indicate that these peaks in the decay energy spectrum are caused by multiple resonant states in 31 Ne. Conclusions: In this study, excited states in 31 Ne were observed for the first time. Transitions from calculated shell model final states in 31 Ne to bound states in 30 Ne are in good agreement with the measured decay energy spectrum. Cross-section calculations for the two-proton knockout populating 31 Ne states as well as spectroscopic factors pertaining to the decay of 31 Ne into 30 Ne are used to examine the results within the context of the shell model expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗
NASA NTRS1994

Long-term modulation of galactic cosmic radiation and its model for space exploration

As the human exploration of space has received new attention in the United States, studies find that exposure to space radiation could adversely impact the mission design. Galactic Cosmic Radiation (GCR), with its very wide range of charges and energies, is particularly important for a mission to Mars, because it imposes a stiff mass penalty for spacecraft shielding. Dose equivalent versus shielding thickness calculations, show a rapid initial drop on exposure with thickness, but an asymptotic behavior at a higher shielding thickness. Uncertainties in the radiobiology are largely unknown. For a fixed radiation risk, this leads to large uncertainties in shielding thickness for small uncertainties in estimated dose. In this paper we investigate the application of steady-state, spherically-symmetric diffusion-convection theory of solar modulation to individual measurements of differential energy spectra from 1954 to 1989 in order to estimate the diffusion coefficient, kappa(r,t), as a function of time. We have correlated the diffusion coefficient to the Climax neutron monitor rates and show that, if the diffusion coefficient can be separated into independent functions of space and time: kappa(r,t)=K(t) k(sub 0) beta Pk(sub 1) (r), where beta is the particle velocity and P the rigidity, then (i) The time dependent quantity 1/K(t), which is proportional to the deceleration potential, phi(r,t), is linearly related to the Climax neutron monitor counting rate. (ii) The coefficients obtained from hydrogen or helium intensity measurements are the same. (iii) There are different correlation functions for odd and even solar cycles. (iv) The correlation function for the Climax neutron monitor counting rate for given time, t, can be used to estimate mean deceleration parameter phi(t) to within +/- 15% with 90% confidence. We have shown that k(r,t) determined from hydrogen and/or helium data, can be used to fit the oxygen and iron differential energy spectra with a root mean square error of about +/- 10%, and essentially independent of the particle charge or energy. We have also examined the ion chamber and C-14 measurements which allow the analysis to be extended from the year 1906 to 1990. Using this model we have defined reference GCR spectra at solar minimum and solar maximum. These can be used for space exploration studies and provide a quantitative estimate of the error in dose due to changes in GCR intensities.

Badhwan, G. D.↗
NASA NTRSDate not supplied

Flow Regime Transitions and Pressure Drop in Gas-Liquid Flows through Packed Beds in Microgravity

The Packed Bed Reactor Experiment (PBRE) focused on the hydrodynamics of gas-liquid flow in porous media, specifically on flow regime transition and pressure drop correlation. The reactor included a cylindrical test section of 60 cm in length, 5 cm in diameter and filled with 2 mm spherical glass beads. Liquid water and Nitrogen gas were flowed into a centrifugal phase separator after exiting the test section. Five absolute pressure transducers and two high-speed video cameras spanned the test section to measure the pressure drop and capture flow regime behavior. Video observations, analysis of pressure traces, and pressure drop data were used to identify flow patterns. It was found that the slope of pressure gradient data when plotted versus the gas or liquid Reynolds numbers changed when the flow regime changed. This observation was used to identify the flow regime boundaries. At high liquid flow rates (high interaction regime), dispersed bubble and pulse flow regimes were detected. At low liquid flow rates (low interaction regime), large bubble and gas channeling regimes were observed. The transition from the big bubble regime to the gas channeling regime occurred at a critical gas flow rate independent of the liquid flow rate. Similarly, the transition from the low interaction to the high interaction regime occurred at a critical liquid flow rate. For these different flow regimes, different two-phase flow friction factors and pressure gradient correlations were fitted as a function of the gas and liquid modified Reynolds numbers and Suratman number. It was found that in all the flow regimes the capillary contribution to the pressure drop is significant and can be as high as 90% of the total pressure drop.

Brian Motil↗
DOE OSTI2021

COMPASO : A new halo finder for competitive assignment to spherical overdensities

We describe a new method (compaso) for identifying groups of particles in cosmological N-body simulations. compaso builds upon existing spherical overdensity (SO) algorithms by taking into consideration the tidal radius around a smaller halo before competitively assigning halo membership to the particles. In this way, the compaso finder allows for more effective deblending of haloes in close proximity as well as the formation of new haloes on the outskirts of larger ones. This halo-finding algorithm is used in the abacussummit suite of N-body simulations, designed to meet the cosmological simulation requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. compaso is developed as a highly efficient on-the-fly group finder, which is crucial for enabling good load-balancing between the GPU and CPU and the creation of high-resolution merger trees. In this paper, we describe the halo-finding procedure and its particular implementation in abacus, accompanying it with a qualitative analysis of the finder. We test the robustness of the compaso catalogues before and after applying the cleaning method described in an accompanying paper and demonstrate its effectiveness by comparing it with other validation techniques. We then visualize the haloes and their density profiles, finding that they are well fit by the NFW formalism. Finally, we compare other properties such as radius–mass relationships and two-point correlation functions with that of another widely used halo finder, rockstar.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

New Features of the NEQAIR Radiation Code

The longest-lived code for predicting shock layer radiation, NEQAIR, is now in its 5th decade of service. Substantial changes to the code have been made over the previous decade, the most recent report of which was at the 5th Workshop on Radiation in High Temperature Gases in 2014, for the version referred to as NEQAIR14. This paper will review some of the improvements made to the NEQAIR code since then, which is now at v15.2. Some of these features are discussed briefly below. NEQAIR15 and subsequent versions have enabled parallel evaluation of multiple lines of sight. This is accomplished by utilizing the HDF5 file format and placing multiple lines into a single file, LOS.h5, which is used for both input and output. This approach enables straightforward parallel execution both over the number of lines of sight and the number of points per line. For large problems, runtime reduces linearly with the number of nodes deployed since each line is processed independently by a subset of MPI ranks. Three applications of the multi-line solver are discussed. The first has to do with performing loosely coupled radiation-flowfield solutions. In this case the computed absorption and emission coefficients are used to evaluate the total energy absorbed or emitted at each point, allowing evaluation of the volumetric source term in the flowfield. The second computation is for obtaining heat flux from nonuniform flows, which require integration over spherical co-ordinates. These are of particular interest for evaluating radiation on the vehicle backshell. This 3D option improves the angular integration scheme and allows adaptive line selection that together reduce the number of lines required by about an order of magnitude. The final application is for remote observation, which is essentially the 3D integration problem over a small solid angle. For all three of these computations, data can be stored in the HDF5 file which allows a NEQAIR run to be restarted when it times out, or to add atmospheric absorption or instrument scan functions. An additional level of parallelism is enabled in NEQAIR15.2 using GPU routines. The GPU parallelism has realized up to 8x speed-up when running on a single core but diminishes as CPU parallelism is increased. For running multi-line simulations, it may be easier to reserve a large number of CPU nodes than to obtain the number of GPU nodes required for similar performance. A GUI, known as NEQTPY, allows for reading and creating input files, running NEQAIR, and displaying results. A significant feature of NEQTPY is the ability to perform spectral fits to data. The fits can operate on a single line spectrum (radiance vs. wavelength) or a 3D input file with multiple columns of data. Other new features include improved constants, additional species, more detailed non-Boltzmann modelling, advanced user controls, the ability to read and calculate spectra from HITRAN datafiles, photodissociation and photoionization cross-sections. A “fast” automatic grid option may reduce the size and time of spectral calculations while still maintaining good accuracy for total heat flux.

Brett A Cruden↗
NASA NTRS1994

Using FU Orionis outbursts to constrain self-regulated protostellar disk models

One-dimensional, convective, vertical structure models and one dimensional time-dependent, radial diffusion models are combined to create a self-consistent picture in which FU Orionis outbursts occur in young stellar objects (YSOs) as the result of a large-scale, self-regulated, thermal ionization instability in the surrounding protostellar accretion disk. Although active accretion disks have long been postulated to be ubiqitous among low-mass YSOs, few constraints have until now been imposed on physical conditions in these disks. By fitting the results of time-dependent disk models to observed timescales of FU Orionis events, we estimate the magnitude of the effective viscous stress in the inner disk (r approximately less than 1 AU) to be, in accordance with an ad hoc 'alpha' prescription, the product of the local sound speed, pressure scale height, and an efficiency factor alpha of 10(exp -4) where hydrogen is neutral and 10(exp 3) where hydrogen is ionized. We hypothesize that all YSOs receive infall onto their outer disks which is steady (or slowly declining with time) and that FU Orionis outbursts are self-regulated, disk outbursts which occur only in systems which transport matter inward at a rate sufficiently high to cause hydrogen to be ionized in the inner disk. We estimate a critical mass flux of dm(sub crit)/dt = 5 x 10(exp 7) solar mass/yr independent of the magnitude of alpha for systems with one solar mass, three solar radius central objects. Infall accretion rates in the range of dm(sub in)/dt = 1-10) x 10(exp -6) solar mass/yr produce observed FU Orionis timescales consistent with estimates of spherical molecular cloud core collapse rates. Modeled ionization fronts are typically initiated near the inner edge of the disk and propogate out to a distance of several tens of stellar radii. Beyond this region, the disk transports mass steadily inward at the supplied constant infall rate. Mass flowing through the innermost disk annulus is equal to dm(sub in)/dt only in a time-averaged sense and is regulated by the ionization of hydrogen in the inner disk such that long intervals (approximately 1000 yr) of low-mass flux: (1-30) x 10(exp -8) solar mass/yr are punctuated by short intervals (approximately 100 yr) of high-mass flux: (1-30) x 10(exp -5) solar mass/yr. Timescales and mass fluxes derived for quiescent and outburst stages are consistent with estimates from observations of T Tauri and FU Orionis systems, respectively.

Bell, K. R.↗
DOE OSTI2023

Spectroscopic quadrupole moments in 124 $\mathrm{Xe}$

The Xe isotopic chain with four valence protons above the Z=50 shell closure is an ideal laboratory for the study of the evolution of nuclear deformation. At the N=82 shell closure, 136 Xe presents all characteristics of a doubly closed shell nucleus with a spherical shape. In the very neutron-deficient isotopes close to N=50, the α-decay chain of Xe was investigated to probe the radioactive decay properties near the drip-line and the magicity of 100 Sn. Additionally, the Xe isotopes present higher order symmetries in the nuclear deformation such as the octupole degree of freedom near N=60 and N=90 or O(6) symmetry in stable isotopes. The relevance of the O(6) symmetry has been investigated by measuring the spectroscopic quadrupole moment of the first excited states in 124 Xe. In the O(6) symmetry limit, the spectroscopic quadrupole moment of collective states is expected to be null. A stable 124 Xe beam with energies of 4.03A MeV and 4.11 A MeV was used to bombard a nat W target at the GANIL facility. Excited states were populated via the safe Coulomb excitation reaction. The collision of the heavy ions with a large Z at low energy make this reaction sensitive to the diagonal E2 matrix element of the excited states. The recoils were detected in the VAMOS++ magnetic spectrometer and the γ rays in the AGATA tracking array. The least squares fitting code gosia was used for the analysis to extract both E2 and M1 transitional and E2 diagonal matrix elements. The rotational ground state band was populated up to the 8$^{+}_{1}$ state as well as the 2$^{+}_{2}$ and 4$^{+}_{2}$ states. Using high precision spectroscopic data to constrain the gosia fit, the spectroscopic quadrupole moments of the 2$^{+}_{1}$, 4$^{+}_{1}$, and 6$^{+}_{1}$ states were determined for the first time. The spectroscopic quadrupole moments were found to be negative, large, and constant in the ground state band underlining the prolate axially deformed ground state band of 124 Xe. Here, the present experimental data confirm that the O(6) symmetry is substantially broken in 124 Xe .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗