SEARCH · Search NASA
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.
A mass shift method using ICP-MS/MS with O 2 reaction gas is used for ultra-trace determinations of U and Th in samples containing W, Ir, Pt, and Au derived polyatomic interferants without the need for extensive chemical sample preprocessing.
With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.
The application scope of the Poincare-Lighthill-Kuo (PLK) method is suggested by using the Particle-in-cell (PIC) numerical method to study head-on collision of two solitary waves. Comparisons between the numerical results from PIC simulations and the analytical ones from the PLK method indicate that the both are in good agreement with each other. The dependence of the phase shifts after the head-on collision on both amplitudes of two solitary waves is given from our PIC method. It is found that the phase shifts depended on the amplitude of both waves. The maximum amplitude during the colliding process is approximately equal to the sum of both amplitudes for the small amplitude solitary waves.
An improved method for measuring the cross sections for carrier trapping at defects in semiconductors is described. This method, a variation of deep level transient spectroscopy (DLTS) used with bipolar transistors, is applied to hot carrier trapping at vacancy-oxygen, carbon-oxygen, and three charge states of divacancy centers (V{sub 2}) in n- and p-type silicon. Unlike standard DLTS, we fill traps by injecting carriers into the depletion region of a bipolar transistor diode using a pulse of forward bias current applied to the adjacent diode. We show that this technique is capable of accurately measuring a wide range of capture cross sections at varying electric fields due to the control of the carrier density it provides. Because this technique can be applied to a variety of carrier energy distributions, it should be valuable in modeling the effect of radiation-induced generation-recombination currents in bipolar devices.
We extend the solid-state nudged elastic band method to handle a non-conserved order parameter, in particular, magnetization, that couples to volume and leads to many observed effects in magnetic systems. We apply this formalism to the well-studied magneto-volume collapse during the pressure-induced transformation in iron—from ferromagnetic body-centered cubic (bcc) austenite to hexagonal close-packed (hcp) martensite. We find a bcc-hcp equilibrium coexistence pressure of 8.4 GPa, with the transition-state enthalpy of 156 meV/Fe at this pressure. A discontinuity in magnetization and coherent stress occurs at the transition state, which has a form of a cusp on the potential-energy surface (yet all the atomic and cell degrees of freedom are continuous); the calculated pressure jump of 25 GPa is related to the observed 25 GPa spread in measured coexistence pressures arising from martensitic and coherency stresses in samples. Our results agree with experiments, but necessarily differ from those arising from drag and restricted parametrization methods having improperly constrained or uncontrolled degrees of freedom.
We present a time-implicit numerical method to solve the relativistic Vlasov–Ampere system of equations on a two dimensional phase space grid. The time-splitting algorithm we use allows the generalization of the work presented here to higher dimensions keeping the linear aspect of the resulting discrete set of equations. The implicit method is benchmarked against linear theory results for the relativistic Landau damping for which analytical expressions using the Maxwell-Jüttner distribution function are derived. We note that, independently from the shape of the distribution function, the relativistic treatment features collective behaviours that do not exist in the nonrelativistic case. The numerical study of the relativistic two-stream instability completes the set of benchmarking tests.
Electromagnetic interaction between a sub-wavelength particle (the “probe”) and a material surface (the “sample”) is studied theoretically. The interaction is shown to be governed by a series of resonances corresponding to surface polariton modes localized near the probe. The resonance parameters depend on the dielectric function and geometry of the probe as well as on the surface reflectivity of the material. Calculation of such resonances is carried out for several types of axisymmetric probes: spherical, spheroidal, and pear-shaped. For spheroids, an efficient numerical method is developed, capable of handling cases of large or strongly momentum-dependent surface reflectivity. Application of the method to highly resonant materials, such as aluminum oxide (by itself or covered with graphene), reveals a rich structure of multi-peak spectra and nonmonotonic approach curves, i.e., the probe-sample distance dependence. These features also strongly depend on the probe shape and optical constants of the model. For less resonant materials such as silicon oxide, the dependence is weak, so that the spheroidal model is reliable. The calculations are done within the quasistatic approximation with radiative damping included perturbatively.
We develop a new halo-independent strategy for analyzing emerging DM hints, utilizing the method of extended maximum likelihood. This approach does not require the binning of events, making it uniquely suited to the analysis of emerging DM direct detection hints. It determines a preferred envelope, at a given confidence level, for the DM velocity integral which best fits the data using all available information and can be used even in the case of a single anomalous scattering event. All of the halo-independent information from a direct detection result may then be presented in a single plot, allowing simple comparisons between multiple experiments. This results in the halo-independent analogue of the usual mass and cross-section plots found in typical direct detection analyses, where limit curves may be compared with best-fit regions in halo-space. The method is straightforward to implement, using already-established techniques, and its utility is demonstrated through the first unbinned halo-independent comparison of the three anomalous events observed in the CDMS-Si detector with recent limits from the LUX experiment.
Here, we demonstrate a multipole expansion method to calculate collective neutrino oscillations in supernovae using the neutrino bulb model. We show that it is much more efficient to solve multi-angle neutrino oscillations in multipole basis than in angle basis. The multipole expansion method also provides interesting insights into multi-angle calculations that were accomplished previously in angle basis.
This paper proposes a new formulation and solution algorithm that uses transmission level solar inverters to address the security-constrained optimal power flow (SCOPF) problem. The goal is to stabilize voltage fluctuations in transmission networks in base case and contingency scenarios, by using bulk solar power plant with a minimal number of post-contingency corrections. To achieve this goal, a two-stage volt/var optimization method is proposed to first correct all voltage violations with the volt-var alternating current optimal power flow (ACOPF) algorithm for a base case. Then a linearized SCOPF volt-var control algorithm is proposed to identify the corrective actions for all potential voltage violations in all contingency scenarios. The proposed method was tested and validated on a modified IEEE 118-bus system with solar photovoltaic (PV) data.
While commonly used in Si solar cell characterization, the Suns–Voc method of determining series resistance (Rs) is rarely used for CdTe and its alloys. However, it is advantageous relative to the slope method—commonly used in the CdTe community—because it measures series resistance at the maximum power point, at which a solar cell operates. Using the self-consistent predecessor to Suns–Voc, Jsc –Voc, , we determine Rs of two different CdSeTe/CdTe solar cell structures to be 0.96 and 2.72 Ωcm2, with repeated measurements showing little deviation
To accurately characterize shielded special nuclear materials (SNM) using passive gamma-ray spectroscopy measurement techniques, the effective atomic number and the thickness of shielding materials must be measured. Intervening materials between the source and detector may affect the estimated source isotopics (uranium enrichment and plutonium grade) for techniques which rely on raw count rates or photopeak ratios of gamma-ray lines separated in energy. Furthermore, knowledge of the surrounding materials can provide insight regarding the configuration of a device containing SNM. The described method was developed using spectra recorded using high energy resolution CdZnTe detectors, but can be expanded to any gamma-ray spectrometers with energy resolution of better than 1% FWHM at 662 keV. The effective atomic number, Z, and mass thickness of the intervening shielding material are identified by comparing the relative attenuation of different gamma-ray lines and estimating the proportion of Compton scattering interactions to photoelectric absorptions within the shield. While characteristic Kα x-rays can be used to identify shielding materials made of high Z elements, this method can be applied to all shielding materials. This algorithm has adequately estimated the effective atomic number for shields made of iron, aluminum, and polyethylene surrounding uranium samples using experimental data. The mass thicknesses of shielding materials have been estimated with a standard error of less than 1.3 g/cm2 for iron shields up to 2.5 cm thick. The effective atomic number was accurately estimated to 26 ± 5 for all iron thicknesses.
Over the past decades, the increasing penetration of distributed energy resources (DERs) has dramatically changed the power load composition in the distribution networks. The traditional static and dynamic load models can hardly capture the dynamic behavior of modern loads especially for fault-induced delayed voltage recovery (FIDVR) events. Thus, a more comprehensive composite load model with combination of static load, different types of induction motors, single-phase A/C motor, electronic load and DERs has been proposed by Western Electricity Coordinating Council (WECC). However, due to the large number of parameters and model complexity, the WECC composite load model (WECC CMLD) raises new challenges to power system studies. To overcome these challenges, in this paper, a cutting-edge parameter reduction (PR) approach for WECC CMLD based on active subspace method (ASM) is proposed. Firstly, the WECC CMLD is parameterized in a discrete-time manner for the application of the proposed method. Then, parameter sensitivities are calculated by discovering the active subspace, which is a lower-dimensional linear subspace of the parameter space of WECC CMLD in which the dynamic response is most sensitive. The interdependency among parameters can be taken into consideration by our approach. Finally, the numerical experiments validate the effectiveness and advantages of the proposed approach for WECC CMLD model.
In the last decade, zeolitic imidazolate frameworks (ZIFs) have been studied extensively for their potential as selective separation membranes. In this review, we highlight unique structural properties of ZIFs that allow them to achieve certain important separations, like that of propylene from propane, and summarize the state of the art in ZIF thin-film deposition on porous substrates and their modification by postsynthesis treatments. We also review the reported membrane performance for representative membrane synthesis approaches and attempt to rank the synthesis methods with respect to potential for scalability. To compare the dependence of membrane performance on membrane synthesis methods and operating conditions, we map out fluxes and separation factors of selected ZIF-8 membranes for propylene/propane separation. Finally, we provide future directions considering the importance of further improvements in scalability, cost effectiveness, and stable performance under industrially relevant conditions.
A major bottleneck to fully understanding the functional aspects of lysine acetylation is the lack of stoichiometry information. Here we describe a mass spectrometry method using a combination of isotope labeling and detection of a diagnostic fragment ion to determine the stoichiometry of lysine acetylation on proteins globally. Using this technique, we determined the modification occupancy on hundreds of acetylated peptides from cell lysates and cross-validated the measurements via immunoblotting.
The sensitivity and uncertainty analysis tools of the ORNL SCALE nuclear modeling and simulation code system that have been developed over the last decade have proven indispensable for numerous application and design studies for nuclear criticality safety and reactor physics. SCALE contains tools for analyzing the uncertainty in the eigenvalue of critical systems, but cannot quantify uncertainty in important neutronic parameters such as multigroup cross sections, fuel fission rates, activation rates, and neutron fluence rates with realistic three-dimensional Monte Carlo simulations. A more complete understanding of the sources of uncertainty in these design-limiting parameters could lead to improvements in process optimization, reactor safety, and help inform regulators when setting operational safety margins. A novel approach for calculating eigenvalue sensitivity coefficients, known as the CLUTCH method, was recently explored as academic research and has been found to accurately and rapidly calculate sensitivity coefficients in criticality safety applications. The work presented here describes a new method, known as the GEAR-MC method, which extends the CLUTCH theory for calculating eigenvalue sensitivity coefficients to enable sensitivity coefficient calculations and uncertainty analysis for a generalized set of neutronic responses using high-fidelity continuous-energy Monte Carlo calculations. Here, several criticality safety systems were examined to demonstrate proof of principle for the GEAR-MC method, and GEAR-MC was seen to produce response sensitivity coefficients that agreed well with reference direct perturbation sensitivity coefficients.