Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral element method”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Relationships between pigment composition variation and reflectance for plant species from a coastal savannah in California

Advances in imaging spectroscopy have indicated that remotely sensed reflectance measurements of the plant canopy may be used to identify and qualify some classes of canopy biochemicals; however, the manner in which differences in biochemical compositions translate into differences is not well understood. Most frequently, multiple linear regression routines have been used to correlate narrow band reflectance values with measured biochemical concentrations. Although some success has been achieved with such methods for given data sets, the bands selected by multiple regression are not consistent between data sets, nor is it always clear what physical or biological basis underlies the correlation. To examine the relationship between biochemical concentration and leaf reflectance signal we chose to focus on the visible spectrum where the primary biochemical absorbances are due to photosynthetic pigments. Pigments provide a range of absorbance features, occur over a range of concentrations in natural samples, and are ecophysiologically important. Concentrations of chlorophyll, for example, have been strongly correlated to foliar nitrogen levels within a species and to photosynthetic capacity across many species. In addition pigments effectively absorb most of the photosynthetically active radiation between 400-700 nm, a spectral region for which silicon detectors have good signal/noise characteristics. Our strategy has been to sample a variety of naturally occurring species to measure leaf reflectance and pigment compositions. We hope to extend our understanding of pigment reflectance effects to interpret small overlapping absorbances of other biochemicals in the infrared region. For this reason, selected samples were also tested to determine total nitrogen, crude protein, cellulose, and lignin levels. Leaf reflectance spectra measured with AVIRIS bandwidths and wavelengths were compared between species and within species and for differences between seasons, for changes in the the shape of the spectra. We attempt to statistically correlate these shape changes with differences in pigment compositions. In parallel with our comparisons of pigment composition and leaf reflectance, we have modified the PROSPECT leaf reflectance model to test the contributions of pigments or pigment group concentrations. PROSPECT considers a leaf as a multi-layer dielectric plane with an uneven surface. Jacquemoud adapted the basic analysis of Allen for surface effects, a leaf thickness factor, and the absorption of water and chlorophyll (actually all pigments) and the plant matrix. Our modifications to PROSPECT in the forward direction include breaking out the pigment concentration parameter into separate components for chlorophyll a and b and a number of xanthophylls and carotenes, and introducing a shift and convolution function to model the spread and shift from their in vitro measurements to their in vivo state. Further, we have considered how the matrix elements (i.e., all biochemicals and structural effects not modeled explicity) vary with species.

Ustin, Susan L.↗
DOE OSTI2016

THE PHYSICAL MECHANISM BEHIND M DWARF METALLICITY INDICATORS AND THE ROLE OF C AND O ABUNDANCES

We present near-infrared (NIR) synthetic spectra based on PHOENIX stellar atmosphere models of typical early and mid-M dwarfs with varied C and O abundances. We apply multiple recently published methods for determining M dwarf metallicity to our models to determine the effects of C and O abundances on metallicity indicators. We find that the pseudo-continuum level is very sensitive to C/O and that all metallicity indicators show a dependence on C and O abundances, especially in lower T {sub eff} models. In some cases, the inferred metallicity ranges over a full order of magnitude (>1 dex) when [C/Fe] and [O/Fe] are varied independently by ±0.2. We also find that [(O−C)/Fe], the difference in O and C abundances, is a better tracer of the pseudo-continuum level than C/O. Models of mid-M dwarfs with [C/Fe], [O/Fe], and [M/H] that are realistic in the context of galactic chemical evolution suggest that variation in [(O−C)/Fe] is the primary physical mechanism behind the M dwarf metallicity tracers investigated here. Empirically calibrated metallicity indicators are still valid for most nearby M dwarfs due to the tight correlation between [(O−C)/Fe] and [Fe/H] evident in spectroscopic surveys of solar neighborhood FGK stars. Variations in C and O abundances also affect the spectral energy distribution of M dwarfs. Allowing [O/Fe] to be a free parameter provides better agreement between the synthetic spectra and observed spectra of metal-rich M dwarfs. We suggest that flux-calibrated, low-resolution, NIR spectra can provide a path toward measuring C and O abundances in M dwarfs and breaking the degeneracy between C/O and [Fe/H] present in M dwarf metallicity indicators.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Orbits for 18 Visual Binaries and Two Double-line Spectroscopic Binaries Observed with HRCAM on the CTIO SOAR 4 m Telescope, Using a New Bayesian Orbit Code Based on Markov Chain Monte Carlo

We present orbital elements and mass sums for 18 visual binary stars of spectral types B to K (five of which are new orbits) with periods ranging from 20 to more than 500 yr. For two double-line spectroscopic binaries with no previous orbits, the individual component masses, using combined astrometric and radial velocity data, have a formal uncertainty of ∼0.1 M{sub ⊙}. Adopting published photometry and trigonometric parallaxes, plus our own measurements, we place these objects on an H-R diagram and discuss their evolutionary status. These objects are part of a survey to characterize the binary population of stars in the Southern Hemisphere using the SOAR 4 m telescope+HRCAM at CTIO. Orbital elements are computed using a newly developed Markov chain Monte Carlo (MCMC) algorithm that delivers maximum-likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate the uncertainty of our derived parameters in a robust way. For spectroscopic binaries, using our approach, it is possible to derive a self-consistent parallax for the system from the combined astrometric and radial velocity data (“orbital parallax”), which compares well with the trigonometric parallaxes. We also present a mathematical formalism that allows a dimensionality reduction of the feature space from seven to three search parameters (or from 10 to seven dimensions—including parallax—in the case of spectroscopic binaries with astrometric data), which makes it possible to explore a smaller number of parameters in each case, improving the computational efficiency of our MCMC code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Frequency Domain Reflectometry (FDR) Simulation Techniques for Digital Twin Representation of an Electrical Cable

Simulation of cable system response to frequency domain reflectometry (FDR) tests can be instrumental to understanding these tests and the nature and influence of various cable anomalies on test signatures. Reflectometry simulations are based upon a finite element representation of cable conductors and insulation to produce an S-parameter at each evaluated frequency. The aggregate collection of cable model S-parameters can simulate the influence on a test signal injected into a physical cable. Such an approach was undertaken in this work to produce a digital twin simulation of a low-voltage electrical cable. The electrical cable digital twin examined the influence of test simulation parameters and the relative influence of cable anomalies, including thermal aging, water or moisture exposure, water or moisture ingress, and other anomalies. The digital twin in this work included modeling of the conductors, insulation, jacket, and surrounding environment (air, water, etc.). The digital twin could be expanded to include cable bends, junctions and splices, branch or T systems, and termination impedances of motors or instruments. Observations and conclusions of this work include: 1. Fully 3D digital twin simulation of an electrical cable using an FDR approach is possible. However, there are tradeoffs between simulation fidelity and solution time, which must be balanced to ensure the simulation solves in an adequate amount of time (e.g., less than 20 minutes). Simulation parameters to balance include frequency bandwidth, number of frequencies, mesh density, connection impedance, and permittivity tolerance. 2. The digital twin simulation can explain FDR sensitivity to various cable anomalies, including entry and exit from an oven or water bath. 3. The digital twin simulation FDR response attenuates with distance along the cable and is further affected by the frequency bandwidth, which is similar to that observed with physical measurements. 4. The resolution of the digital twin FDR peaks increased with increasing bandwidth and with increasing number of frequencies, again similar to physical measurements. 5. The presence of multiple anomalies in the digital twin does not substantially attenuate the FDR response to anomalies located beyond the first encountered anomaly and impedance mismatch. 6. Spectral variation of the permittivity did not have a significant effect on the FDR response compared to a fixed nominal value. 7. Extension of the digital twin to 1000 ft still allowed for detection of distal anomalies near the far end of the electrical cable from the instrument connection point. 8. The ARENA test bed facilitates efficient NDE evaluations of well understood cable anomalies with various NDE methods without risking actual plant damage.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Rod internal pressure quantification and distribution analysis using Frapcon

This report documents work performed supporting the Department of Energy (DOE) Office of Nuclear Energy (NE) Fuel Cycle Technologies Used Fuel Disposition Campaign (UFDC) under work breakdown structure element 1.02.08.10, ST Analysis. In particular, this report fulfills the M4 milestone M4FT- 15OR0810036, Quantify effects of power uncertainty on fuel assembly characteristics, within work package FT-15OR081003 ST Analysis-ORNL. This research was also supported by the Consortium for Advanced Simulation of Light Water Reactors (http://www.casl.gov), an Energy Innovation Hub (http://www.energy.gov/hubs) for Modeling and Simulation of Nuclear Reactors under U.S. Department of Energy Contract No. DE-AC05-00OR22725. The discharge rod internal pressure (RIP) and cladding hoop stress (CHS) distributions are quantified for Watts Bar Nuclear Unit 1 (WBN1) fuel rods by modeling core cycle design data, operation data (including modeling significant trips and downpowers), and as-built fuel enrichments and densities of each fuel rod in FRAPCON-3.5. A methodology is developed which tracks inter-cycle assembly movements and assembly batch fabrication information to build individual FRAPCON inputs for each evaluated WBN1 fuel rod. An alternate model for the amount of helium released from the zirconium diboride (ZrB2) integral fuel burnable absorber (IFBA) layer is derived and applied to FRAPCON output data to quantify the RIP and CHS for these types of fuel rods. SCALE/Polaris is used to quantify fuel rodspecific spectral quantities and the amount of gaseous fission products produced in the fuel for use in FRAPCON inputs. Fuel rods with ZrB2 IFBA layers (i.e., IFBA rods) are determined to have RIP predictions that are elevated when compared to fuel rod without IFBA layers (i.e., standard rods) despite the fact that IFBA rods often have reduced fill pressures and annular fuel pellets. The primary contributor to elevated RIP predictions at burnups less than and greater than 30 GWd/MTU is determined to be the total fuel rod void volume and the amount of released fission gas in the fuel rod, respectively. Cumulative distribution functions (CDFs) are prepared from the distribution of RIP and CHS predictions for all standard and IFBA rods. The provided CDFs allow for the determination of the portion of WBN1 fuel rods that exceed a specified RIP or CHS limit. Results are separated into IFBA and standard rods so that the two groups may be analyzed individually. FRAPCON results are provided in sufficient detail to enable the recalculation of the RIP while considering any desired plenum gas temperature, total void volume, or total amount of gas present in the void volume. A method to predict the CHS from a determined or assumed RIP is also proposed, which is based on the approximately linear relationship between the CHS and the RIP. Finally, improvements to the computational methodology of FRAPCON are proposed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Hitomi observations of the LMC SNR N 132 D: Highly redshifted X-ray emission from iron ejecta

Here, we present Hitomi observations of N 132 D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16–17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km s –1 compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50–1500 km s –1 if a weakly informative prior is placed on possible line broadening. This indicates that the Fe emission arises from the supernova ejecta, and that these ejecta are highly asymmetric, since no blueshifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

PROSPECTING IN LATE-TYPE DWARFS: A CALIBRATION OF INFRARED AND VISIBLE SPECTROSCOPIC METALLICITIES OF LATE K AND M DWARFS SPANNING 1.5 dex

Knowledge of late K and M dwarf metallicities can be used to guide planet searches and constrain planet formation models. However, the determination of metallicities of late-type stars is difficult because visible wavelength spectra of their cool atmospheres contain many overlapping absorption lines, preventing the measurement of equivalent widths. We present new methods, and improved calibrations of existing methods, to determine metallicities of late K and M dwarfs from moderate resolution (1300 < R < 2000) visible and infrared spectra. We select a sample of 112 wide binary systems that contain a late-type companion to a solar-type primary star. Our sample includes 62 primary stars with previously published metallicities, as well as 50 stars with metallicities determined from our own observations. We use our sample to empirically determine which features in the spectrum of the companion are best correlated with the metallicity of the primary. We find {approx_equal}120 features in K and M dwarf spectra that are useful for predicting metallicity. We derive metallicity calibrations for different wavelength ranges, and show that it is possible to get metallicities reliable to <0.10 dex using either visible, J-, H-, or K-band spectra. We find that the most accurate metallicities derived from visible spectra requires the use of different calibrations for early-type (K5.5-M2) and late-type (M2-M6) dwarfs. Our calibrations are applicable to dwarfs with metallicities of -1.04 < [Fe/H] <+0.56 and spectral types from K7 to M5. Lastly, we use our sample of wide binaries to test and refine existing calibrations to determine M dwarf metallicities. We find that the {zeta} parameter, which measures the ratio of TiO can CaH bands, is correlated with [Fe/H] for super-solar metallicities, and {zeta} does not always correctly identify metal-poor M dwarfs. We also find that existing calibrations in the K and H bands are quite reliable for stars with [Fe/H] >-0.5, but are less useful for more metal-poor stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Scalable Interior‐Point Gauss–Newton Method for PDE‐Constrained Optimization With Bound Constraints

Here, we present a scalable approach to solve a class of partial differential equation (PDE)‐constrained optimization problems with bound constraints. This approach utilizes a robust full‐space interior‐point (IP)‐Gauss–Newton optimization method. To cope with the poorly‐conditioned IP‐Gauss–Newton saddle‐point linear systems that need to be solved approximately, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE‐constrained optimization problems. A block Gauss–Seidel preconditioner is proposed for the GMRES‐based solution of the IP‐Gauss–Newton linear systems. It is shown, for a large‐class of PDE‐ and bound‐constrained optimization problems, that the spectrum of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix is asymptotically independent of discretization and is not impacted by the ill‐conditioning that notoriously plagues interior‐point methods. We exploit symmetry of the IP‐Gauss–Newton linear systems and propose a regularization and log‐barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)‐based solution of the equivalent IP‐Gauss–Newton–Schur complement linear systems. The eigenvalues of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix, that are not equal to one, are identical to the eigenvalues of the regularization and log‐barrier Hessian preconditioned Schur complement matrix. The scalability of the approach is demonstrated on two example problems. The numerical solution of these optimization problems is shown to require a discretization independent number of IP‐Gauss–Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill‐conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of the preconditioner, achieved via algebraic multigrid component solvers when applicable, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE‐ and bound‐constrained problems.

PDE-constrained optimization↗
NASA NTRS2009

Dynamic Self-Locking of an OEO Containing a VCSEL

A method of dynamic self-locking has been demonstrated to be effective as a means of stabilizing the wavelength of light emitted by a vertical-cavity surface-emitting laser (VCSEL) that is an active element in the frequency-control loop of an optoelectronic oscillator (OEO) designed to implement an atomic clock based on an electromagnetically- induced-transparency (EIT) resonance. This scheme can be considered an alternative to the one described in Optical Injection Locking of a VCSEL in an OEO (NPO-43454), NASA Tech Briefs, Vol. 33, No. 7 (July 2009), page 33. Both schemes are expected to enable the development of small, low-power, high-stability atomic clocks that would be suitable for use in applications involving precise navigation and/or communication. To recapitulate from the cited prior article: In one essential aspect of operation of an OEO of the type described above, a microwave modulation signal is coupled into the VCSEL. Heretofore, it has been well known that the wavelength of light emitted by a VCSEL depends on its temperature and drive current, necessitating thorough stabilization of these operational parameters. Recently, it was discovered that the wavelength also depends on the microwave power coupled into the VCSEL. This concludes the background information. From the perspective that led to the conception of the optical injection-locking scheme described in the cited prior article, the variation of the VCSEL wavelength with the microwave power circulating in the frequency-control loop is regarded as a disadvantage and optical injection locking is a solution of the problem of stabilizing the wavelength in the presence of uncontrolled fluctuations in the microwave power. The present scheme for dynamic self-locking emerges from a different perspective, in which the dependence of VCSEL wavelength on microwave power is regarded as an advantageous phenomenon that can be exploited as a means of controlling the wavelength. The figure schematically depicts an atomic-clock OEO of the type in question, wherein (1) the light from the VCSEL is used to excite an EIT resonance in selected atoms in a gas cell (e.g., 87Rb atoms in a low-pressure mixture of Ar and Ne) and (2) the power supplied to the VCSEL is modulated by a microwave signal that includes components at beat frequencies among the VCSEL wavelength and modulation sidebands. As the VCSEL wavelength changes, it moves closer to or farther from a nearby absorption spectral line, and the optical power transmitted through the cell (and thus the loop gain) changes accordingly. A change in the loop gain causes a change in the microwave power and, thus, in the VCSEL wavelength. It is possible to choose a set of design and operational parameters (most importantly, the electronic part of the loop gain) such that the OEO stabilizes itself in the sense that an increase in circulating microwave power causes the VCSEL wavelength to change in a direction that results in an increase in optical absorption and thus a decrease in circulating microwave power. Typically, such an appropriate choice of operational parameters involves setting the nominal VCSEL wavelength to a point on the shorter-wavelength wing of an absorption spectral line.

Strekalov, Dmitry↗
DOE OSTI2012

A NEW METHOD TO CONSTRAIN SUPERNOVA FRACTIONS USING X-RAY OBSERVATIONS OF CLUSTERS OF GALAXIES

Supernova (SN) explosions enrich the intracluster medium (ICM) both by creating and dispersing metals. We introduce a method to measure the number of SNe and relative contribution of Type Ia supernovae (SNe Ia) and core-collapse supernovae (SNe cc) by directly fitting X-ray spectral observations. The method has been implemented as an XSPEC model called snapec. snapec utilizes a single-temperature thermal plasma code (apec) to model the spectral emission based on metal abundances calculated using the latest SN yields from SN Ia and SN cc explosion models. This approach provides a self-consistent single set of uncertainties on the total number of SN explosions and relative fraction of SN types in the ICM over the cluster lifetime by directly allowing these parameters to be determined by SN yields provided by simulations. We apply our approach to XMM-Newton European Photon Imaging Camera (EPIC), Reflection Grating Spectrometer (RGS), and 200 ks simulated Astro-H observations of a cooling flow cluster, A3112. We find that various sets of SN yields present in the literature produce an acceptable fit to the EPIC and RGS spectra of A3112. We infer that 30.3% {+-} 5.4% to 37.1% {+-} 7.1% of the total SN explosions are SNe Ia, and the total number of SN explosions required to create the observed metals is in the range of (1.06 {+-} 0.34) Multiplication-Sign 10{sup 9} to (1.28 {+-} 0.43) Multiplication-Sign 10{sup 9}, from snapec fits to RGS spectra. These values may be compared to the enrichment expected based on well-established empirically measured SN rates per star formed. The proportions of SNe Ia and SNe cc inferred to have enriched the ICM in the inner 52 kpc of A3112 is consistent with these specific rates, if one applies a correction for the metals locked up in stars. At the same time, the inferred level of SN enrichment corresponds to a star-to-gas mass ratio that is several times greater than the 10% estimated globally for clusters in the A3112 mass range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Origin of Ne emission line of very luminous soft X-ray transient MAXI J0158–744

We investigate the mechanism that reproduced notable spectral features at the ignition phase of the nova explosion observed for a super-Eddington X-ray transient source MAXI J0158–744 in the Small Magellanic Cloud. These features include a strong Ne IX emission line at 0.92 keV with large equivalent widths of 0.32{sub −0.11}{sup +0.21} keV and the absence of the Ne X line at 1.02 keV. In this paper, we calculate the radiative transfer using a Monte Carlo code, taking into account the line blanketing effect due to transitions of N, O, Ne, Mg, and Al ions in an accelerating wind emanating from a white dwarf with a structure based on a spherically symmetric stationary model. We found that the strong Ne IX line can be reproduced if the mass fraction of Ne is enhanced to 10{sup –3} or more and if that of O is reduced to ∼5 × 10{sup –9} or less and that the absence of other lines, including Ne X ions at higher energies, can be also reproduced by the line blanketing effect. This enhancement of the Ne mass fraction indicates that the ejecta are enriched by Ne dredged up from the surface of the white dwarf composed of O, Ne, and Mg rather than C and O, as already pointed out in previous work. We argue that the CNO cycle driving this nova explosion converted most of C and O into N and thus reduced the O mass fraction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A MONTE CARLO MARKOV CHAIN BASED INVESTIGATION OF BLACK HOLE SPIN IN THE ACTIVE GALAXY NGC 3783

The analysis of relativistically broadened X-ray spectral features from the inner accretion disk provides a powerful tool for measuring the spin of supermassive black holes in active galactic nuclei (AGNs). However, AGN spectra are often complex and careful analysis employing appropriate and self-consistent models is required if one has to obtain robust results. In this paper, we revisit the deep 2009 July Suzaku observation of the Seyfert galaxy NGC 3783 in order to study in a rigorous manner the robustness of the inferred black hole spin parameter. Using Monte Carlo Markov chain techniques, we identify a (partial) modeling degeneracy between the iron abundance of the disk and the black hole spin parameter. We show that the data for NGC 3783 strongly require both supersolar iron abundance (Z{sub Fe} = 2-4 Z{sub Sun }) and a rapidly spinning black hole (a > 0.89). We discuss various astrophysical considerations that can affect the measured abundance. We note that, while the abundance enhancement inferred in NGC 3783 is modest, the X-ray analysis of some other objects has found extreme iron abundances. We introduce the hypothesis that the radiative levitation of iron ions in the innermost regions of radiation-dominated AGN disks can enhance the photospheric abundance of iron. We show that radiative levitation is a plausible mechanism in the very inner regions of high accretion rate AGN disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Investigation of Microbunching Instabilities in Modern Recirculating Accelerators

Particle accelerators are machines to accelerate and store charged particles, such as electrons or protons, to the energy levels for various scientific applications. A collection of charged particles usually forms a particle beam. There are three basic types of particle accelerators: linear accelerators (linac), storage-ring (or circular) accelerators, and recirculating accelerators. In a linac, particles are accelerated and pass through once along a linear or straight beamline. Storage-ring accelerators propel particles around a circular track and repetitively append the energy to the stored beam. The third type, also the most recent one in chronology, the recirculating accelerator, is designed to accelerate the particle beam in a short section of linac, circulate the beam, and then either continue to accelerate for energy boost or decelerate it for energy recovery. The beam properties of a linac machine are set at best by the initial particle sources. For storage rings, the beam equilibria are instead determined by the overall machine design. The modern recirculating machines share with linacs the advantages to both accelerate and preserve the beam with high beam quality, as well as efficiently reuse the accelerating components. The beamline design in such a machine configuration can however be much more complicated than that of linacs. As modern accelerators push toward the high-brightness or high-intensity frontier by demanding particles in a highly charged bunch (about nano-Coulomb per bunch) to concentrate in an ever-decreasing beam phase space (transverse normalized emittance about 1 μm and relative energy spread of the order of 10^-5 in GeV beam energy), the interaction amongst particles via their self-generated electromagnetic fields can potentially lead to coherent instabilities of the beam and thus pose significant challenges to the machine design and operation. In the past decade and a half, microbunching instability (MBI) has been one of the most challenging issues for such high-brightness or high-intensity beam transport, as it would degrade lasing performance in the fourth-generation light sources, reduce cooling efficiency in electron cooling facilities, and compromise the luminosity of colliding beams in lepton or lepton-hadron colliders. The dissertation work will focus on the MBI in modern recirculating electron accelerators. It has been known that the collective interactions, the coherent synchrotron radiation (CSR) and the longitudinal space charge (LSC) forces, can drive MBI. The CSR effect is a collective phenomenon in which the electrons in a curved motion, e.g. a bending dipole, emit radiation at a scale comparable to the micro-bunched structure of the bunch distribution. The LSC effect stems from non-uniformity of the charge distribution, acts as plasma oscillation, and can eventually accumulate an amount of energy modulation when the beam traverses a long section of a beamline. MBI can be seeded by non-uniformity or shot noise of the beam, which originates from granularity of the elementary charge. Through the aforementioned collective effects, the modulation of the bunch sub-structure can be amplified and, once the beam-wave interaction formed a positive feedback, can result in MBI. The problem of MBI has been intensively studied for linac-based facilities and for storage-ring accelerators. However, systematic studies for recirculation machines are still very limited and form a knowledge gap. Because of the much more complicated machine configuration of the recirculating accelerators than that of linacs, the existing MBI analysis needs to be extended to accommodate the high-brightness particle beam transport in modern recirculating accelerators. This dissertation is focused on theoretical investigation of MBI in such machine configuration in the following seven themes: (1) Development and generalization of MBI theory The theoretical formulation has been extended so as to be applicable to a general linear beamline lattice including horizontal and vertical transport bending elements, and beam acceleration or deceleration. These featured generalizations are required for MBI analysis in recirculation accelerators. (2) Construction of CSR impedance models In addition to the steady-state CSR interaction, it has been found that the exit transient effect (or CSR drift) can even result in more serious MBI in high-brightness recirculation arcs. The onedimensional free-space CSR impedances, especially the exit transients, are derived. The steady-state CSR impedance is also extended to non-ultrarelativistic beam energy for MBI analysis of low-energy merger sections in recirculating accelerators. (3) Numerical implementation of the derived semi-analytical formulation This includes the development of a semi-analytical Vlasov solver for MBI analysis, and also benchmarking of the solver against massive particle tracking simulations. (4) Exploration of multistage amplification behavior of CSR microbunching development The CSR-induced MBI acts as an amplifier, which amplifies the sub-bunch modulation of a beam. The amplification is commonly quantified by the amplification gain. A beam transport system can be considered as a cascaded amplifier. Unlike the two-stage amplification of four-dipole bunch compressor chicanes employed in linacs, the recirculation arcs, which are usually constituted by several tens of bending magnets, show a distinguishing feature of up to six-stage microbunching amplification for our example arc lattices. That is, the maximal CSR amplification gain can be proportional to the peak bunch current up to sixth power. A method to compare lattice performance has been developed in terms of gain coefficients, which nearly depend on the lattice properties only. This method has also proven to be an effective way to quantify the current dependence of the maximal (5) Control of CSR MBI in multibend transport or recirculation arcs The existing mitigation schemes of MBI mostly aim to linac-based accelerators and may not be practical to the recirculating accelerator facilities. Thus a set of conditions for suppression of CSR MBI was proposed and examined for example lattices from low (~100 MeV) to high (~1 GeV) energies. (6) Study of more aspects of microbunched structures in beam phase spaces For a cascaded amplifier in circuit electronics, the total amplification gain can be estimated as the product of individual gains. In a beam transport line of an accelerator, the (scalar) gain multiplication was examined and found to under-estimate the overall microbunching amplification. The concept of gain matrix was developed, which includes the density, energy and transverse-longitudinal modulations in a beam phase space, and used to analyze MBI for a proposed recirculating machine. Throughout the gain matrix approach, it reasonably gives the upper limit of spectral MBI gain curves. This extended analysis can be employed to study multi-pass recirculation. (7) Study of MBI for magnetized beams Driven by a recent energy-recovery-linac based cooler design for electron cooling at Jefferson Lab Electron-Ion Collider Project, the generalized theoretical formulation for MBI to a transversely coupled beam has been developed and applied to this study. A magnetized beam in general features non-zero canonical angular momentum, thus considered to be a transversely coupled beam. A novel idea of utilizing magnetized beam transport was proposed for improvement of cooling efficiency and possible mitigation of collective effects. A concern of MBI regarding this design was studied and excluded. The large transverse beam size associated with the beam magnetization is found to help suppress MBI via the transverse-longitudinal correlation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD{sub 3}CN treated with a parallel multi-state EVB model

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD{sub 3}CN → DF + CD{sub 2}CN reaction in CD{sub 3}CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD{sub 3}CN solvent, equilibrium power spectra of DF in CD{sub 3}CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol{sup −1} localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD{sub 3}CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

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

Atomic data and spectral modeling constraints from high-resolution X-ray observations of the Perseus cluster with Hitomi

Here, the Hitomi Soft X-ray Spectrometer spectrum of the Perseus cluster, with ~5 eV resolution in the 2–9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic data and models. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that accurate atomic data and models are as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current data and models ready for the data from the next Hitomi-level mission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Assessment of Impurities in Neptunium Metal Using Double-Focusing SF-ICP-MS

In the transmutation process, minor actinides such as neptunium, americium and curium obtained from spent nuclear fuel are incorporated into new fuel that can be placed into a fast reactor, where these elements fission into products with shorter half-lives which reduces the health and safety hazards. In order to better understand the all the involved processes, fundamental physical properties of actinides should be investigated, and for that reason, high-purity materials need to be produced. One of the main challenges in characterizing these materials is the lack of pure metal available for study, and lack of analytical methods to determine the purity. Trans-uranic elements are challenging to study because of the need to be handled in specialized facilities, safeguards concerns and they are not easily purchased and/or produced. High-purity Np metal is being produced at the Materials and Fuels Complex (MFC) of Idaho National Laboratory (INL), in a process of direct oxide reduction developed by Squires and Lessing, with a target purity of 99.999%. The goal of this research is to develop methods to determine the impurities in the Np metal, with detection limits low enough to achieve this purity level. In this work, we have used several ICP-MS and ICP OES techniques, but the focus was the method development using a recently acquired Double-Focusing Sector Field Inductively Coupled Plasma Mass Spectrometer (SF-ICP-MS). To produce the Np metal, approximately 20 g of the neptunium oxide powder was calcined at 800 oC for 6 h in a magnesium oxide crucible, which is placed inside a stainless steel crucible, that goes into the Hot Uniaxial Press (HUP) furnace insert. After the calcining step, the Np oxide was mixed with approximately 9 g of calcium metal and 100 g of ultra-dry calcium chloride. The HUP furnace was heated in a ramp rate of 50 oC min -1 up to 1000 oC. When the reaction mixture reached 850 oC, a tantalum stirrer was activated at 100 rpm. After the tantalum stirrer was removed and the mixture cooled overnight, the MgO crucible was opened with a hammer and chisel to retrieve the Np metal button at the bottom of the reaction mixture. Prior to the dissolution process, the surface contaminants of the sample were removed. The Np button was washed in 10 mL of ultra-pure water for 1 min, with a faint reaction observed on the surface, indicating the presence of Ca metal. The water-washed sample was placed in a tube with 2 M HCl for only 3 s, as a vigorous reaction in the button surface was observed. The leachates were saved for further analysis. The sample was then rinsed with ultra-pure water and air-dried prior to dissolution. The cleaned sample was placed in a PTFE beaker and 25 mL of 6 M HCl and 150 µL of 24 M HF were added. After 10 minutes of reaction, the beaker was placed on a hot plate for 1 h and 30 min, followed by the addition of 1 mL of 16 M HNO3 and allowed to react for a further 1 h and 20 min. The sample was removed from heat and 10 mL of HNO3 16 M was added to react with the HCl. The beaker was covered and left for overnight reaction, producing a solution with approximately 20 mL of final volume. A Quadrupole Inductively Coupled Plasma Mass Spectrometer (Q-ICP-MS, Perkin Elmer ELAN DRC II Q-ICP-MS) was used to obtain qualitative mass spectra for the dissolution blank, each of the leachates and the dissolved sample. In the water leachate spectrum mainly Ca and Mg were observed, while the acid wash had a similar peak pattern to the one observed for the Np metal sample, indicating that the process not only removed surface contamination, but also reacted with the sample itself. In the Np material spectrum, Ta and actinides were identified as candidates for major contaminants, as the most intense peaks corresponded to the mass-to-charge ratios of these elements. Minor contaminants included Al, Au, Ca, Cr, Fe, I, Na, Nb, Pb, Sn, Ti, Tl, V, Y. The presence of impurities in the material could have originated from several sources, such as the materials used for the casting process, the use of impure NpO2 or contamination of reagents and materials used for the analysis. Some of the isotopes listed as potential contaminants are susceptible to spectral interferences, so Double-Focusing SF-ICP-MS (Nu Instruments Attom ES HR-ICP-MS) was used for the quantitative analysis. For each isotope determined, several items were evaluated: raw signal of a blank solution, sensitivity and determination coefficient of calibration curves, possibility of interferences and comparison of results obtained for at least two resolving powers (RP): 300 and 4000. The evaluation of the signal in a blank solution can indicate if a spectral interference from the plasma gas and/or the solvent was successfully resolved (commonly observed for 52Cr+/40Ar12C+, for example). It is important to highlight that a signal reduction is expected when a higher RP is used, as the ion beam transmission is decreased by the use of narrower slits

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

ASSESMENT OF IMPURITIES IN NP METAL USING DOUBLE-FOCUSING SF-ICP-MS

In the transmutation process, minor actinides such as neptunium, americium and curium obtained from spent nuclear fuel are incorporated into new fuel that can be placed into a fast reactor, where these elements fission into products with shorter half-lives which reduces the health and safety hazards [1]. In order to better understand the all the involved processes, fundamental physical properties of actinides should be investigated, and for that reason, high-purity materials need to be produced. One of the main challenges in characterizing these materials is the lack of pure metal available for study, and lack of analytical methods to determine the purity. Trans-uranic elements are challenging to study because of the need to be handled in specialized facilities, safeguards concerns and they are not easily purchased and/or produced. High-purity Np metal is being produced at the Materials and Fuels Complex (MFC) of Idaho National Laboratory (INL), in a process of direct oxide reduction developed by Squires and Lessing [2], with a target purity of 99.999%. The goal of this research is to develop methods to determine the impurities in the Np metal, with detection limits low enough to achieve this purity level. In this work, we have used several ICP-MS and ICP OES techniques, but the focus was the method development using a recently acquired Double-Focusing Sector Field Inductively Coupled Plasma Mass Spectrometer (SF-ICP-MS). To produce the Np metal, approximately 20 g of the neptunium oxide powder was calcined at 800 oC for 6 h in a magnesium oxide crucible, which is placed inside a stainless steel crucible, that goes into the Hot Uniaxial Press (HUP) furnace insert. After the calcining step, the Np oxide was mixed with approximately 9 g of calcium metal and 100 g of ultra-dry calcium chloride. The HUP furnace was heated in a ramp rate of 50 oC min -1 up to 1000 oC. When the reaction mixture reached 850 oC, a tantalum stirrer was activated at 100 rpm. After the tantalum stirrer was removed and the mixture cooled overnight, the MgO crucible was opened with a hammer and chisel to retrieve the Np metal button at the bottom of the reaction mixture [2]. Prior to the dissolution process, the surface contaminants of the sample were removed. The Np button was washed in 10 mL of ultra-pure water for 1 min, with a faint reaction observed on the surface, indicating the presence of Ca metal. The water-washed sample was placed in a tube with 2 M HCl for only 3 s, as a vigorous reaction in the button surface was observed. The leachates were saved for further analysis. The sample was then rinsed with ultra-pure water and air-dried prior to dissolution. The cleaned sample was placed in a PTFE beaker and 25 mL of 6 M HCl and 150 µL of 24 M HF were added. After 10 minutes of reaction, the beaker was placed on a hot plate for 1 h and 30 min, followed by the addition of 1 mL of 16 M HNO3 and allowed to react for a further 1 h and 20 min. The sample was removed from heat and 10 mL of HNO3 16 M was added to react with the HCl. The beaker was covered and left for overnight reaction, producing a solution with approximately 20 mL of final volume. A Quadrupole Inductively Coupled Plasma Mass Spectrometer (Q-ICP-MS, Perkin Elmer ELAN DRC II Q-ICP-MS) was used to obtain qualitative mass spectra for the dissolution blank, each of the leachates and the dissolved sample. In the water leachate spectrum mainly Ca and Mg were observed, while the acid wash had a similar peak pattern to the one observed for the Np metal sample, indicating that the process not only removed surface contamination, but also reacted with the sample itself. In the Np material spectrum, Ta and actinides were identified as candidates for major contaminants, as the most intense peaks corresponded to the mass-to-charge ratios of these elements. Minor contaminants included Al, Au, Ca, Cr, Fe, I, Na, Nb, Pb, Sn, Ti, Tl, V, Y. The presence of impurities in the material could have originated from several sources, such as the materials used for the casting process, the use of impure NpO2 or contamination of reagents and materials used for the analysis. Some of the isotopes listed as potential contaminants are susceptible to spectral interferences, so Double-Focusing SF-ICP-MS (Nu Instruments Attom ES HR-ICP-MS) was used for the quantitative analysis. For each isotope determined, several items were evaluated: raw signal of a blank solution, sensitivity and determination coefficient of calibration curves, possibility of interferences and comparison of results obtained for at least two resolving powers (RP): 300 and 4000. The evaluation of the signal in a blank solution can indicate if a spectral interference from the plasma gas and/or the solvent was successfully resolved (commonly observed for 52Cr+/40Ar12C+, for example). It is important to highlight that a signal reduction is expected when a higher RP is used, as the ion beam transmission is decreased by the use of na

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

Environmental Characterization Scans Using Sodium-Iodide Gamma Spectroscopy to Determine Ra-226 Concentration in Surface Soil - 16326

The recent natural gas drilling boom has brought heightened attention to the issue of TENORM waste management. Radioactive residues from the extraction, treatment, and purification of materials containing uranium and thorium, or their progeny, pose a public health hazard primarily in the form of ground water contamination. Hydraulic fracturing can liberate formation waters from shale; the high salinity of which can cause radium to move with the water to the surface. The disposal and handling of this water has been a subject of public controversy. For soil remediation projects, radium-226 is generally one of the radionuclides of concern due to its alpha emissions. However performing a survey scan for Ra-226 is difficult because its photon emissions are low energy and have minimal probability of emission per decay (186.2 keV, 3.3% abundance). Additionally resolving its emission from U-235 (185.7 keV, 57% abundance) with field instrumentation is not feasible. Therefore detection of Ra-226 relies on the presence and photon emissions of its progeny. This paper outlines methodology for calculating the minimum detectable concentration of Ra-226 in characterization scans using a sodium-iodide detector array coupled with a multi-channel analyzer. Using progeny ingrowth as a detection proxy requires knowing decay time in order to determine progeny emission rates. Knowledge of site history, process, and prior survey results can be key in determining an accurate decay time. However in the absence of such information in-situ measurements can be used to estimate decay time by comparing observed progeny emission ratios against established values. These values are determined using Oak Ridge National Laboratory's ORIGEN-ARP program to generate photon energy spectra for Ra-226 and its progeny at specified time intervals. Once decay time is established, the detector's minimum detectable concentration can be calculated by synthesizing the expected emission spectra and Los Alamos National Laboratory's Monte Carlo N-Particle transport code to simulate detector response. Traditionally, one of the steps in calculating a system's minimum detectable concentration is to place an appropriate button source under the detector, record the response, and then move the source by a specified increment. This process is repeated until the detector's spatial efficiency has been mapped. Due to the impracticality of performing this with Ra-226 button sources of varying decay times, Monte Carlo N-Particle transport code was utilized to determine the detector's spatial efficiency. A dime sized source with the established emission spectra was placed into the center of a rectangle of soil. This rectangle, including the source, was used as a lattice element to fill the region of space below the detector. The simulation was designed to track the energy distribution of pulses created in the sodium iodide crystal using the F8 tally, and was setup to isolate each source; thereby correlating detector response to source position. The simulation output required post-processing to determine conversion factors from counts per second to disintegrations per second. Post-processing accounted for subtraction of the Compton continuum, and average the data over the scanning pattern to obtain the expected spatial sensitivity during a characterization study. This was performed for the following photo-peaks: 609 keV (Bi-214), and 1764 keV (Bi-214). For calculation of the scan minimum detectable concentration, the surveyor efficiency coefficient from NUREG 1507 was supplanted in favor of a scan speed deviation factor. This was used because the detector utilizes computerized data collection (spatial and spectral). This report finds that, for two shielded 10.16 cm x 10.16 cm x 40.64 cm sodium-iodide crystals with integrated multi-channel analyzers, the minimum detectable concentration in surface soils is less than 1 pCi per gram Ra-226 for a 1 second survey measurement. These results assume a volumetric source with a maximum depth of 15.24 cm and a minimum decay time of 21 days. Verification of the system's capabilities was performed by sending soil samples from scanned areas, and locations where in-situ measurements were performed, to an independent laboratory for analysis. The laboratory results demonstrated that scanned area was principally comprised of soil containing less than 1 pCi per gram Ra-226. This borders the system MDC and it underestimated the concentration of radium in the soil. The probable causes of this is escape of radon gas from the surface soil. Radon is the daughter of radium, and has a high mobility in soil. If the gas escapes, then the progeny used as the detection proxy will be present in lower than expected concentrations. Radionuclide identification using high volume sodium-iodide crystals in low level characterization scans is demonstrated to be possible. Further system refinement will focus on developing a methodology to estimate a radon escape factor for site specific use. Accounting for radon escape in the top 15.24 cm of soil will allow the system to better estimate Ra-226 concentrations in soil. This is a significant advancement over the gross counting statistics used in most studies to satisfy MARSSIM survey requirements, because it will increase stakeholder confidence. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗