Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FIT”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Overcoming significant challenges in extracting off-stoichiometric thermodynamics using the compound energy formalism through complementary use of experimental and first principles data: A case study of Ba 1-x Sr x FeO 3-δ

The compound energy formalism (CEF) is a powerful framework to describe the thermodynamics of metal oxides as a function of off-stoichiometry, temperature, and composition. The thermodynamic properties are crucial materials design attributes in metal oxide-based oxygen-exchange chemical processes. Despite the richness of information an accurate CEF model provides, a method to determine a unique and accurate fit for oxygen exchange materials remains elusive. This contribution details a method for fitting the CEF model that overcomes the current fitting challenges through three innovations: 1) the combination of density functional theory calculations with experimental data decorrelates excess terms and delineates the enthalpic/entropic contributions to the Gibbs free energy; 2) a systematic determination of the important CEF model terms, removing thermodynamic predetermining human intervention; 3) a self-consistent solution of the starting oxygen offstoichiometry (δ 0 ) of thermogravimetric measurements. Thus, our method enables the reliable extraction of off-stoichiometric metal oxide thermodynamic properties and facilitates rapid materials compositional screening, and reliable process design of systems dependent on off-stoichiometric redox-active metal oxides. We apply this method to a Ba x Sr 1-x FeO 3-δ test case. We find by systematically examining the performance of the CEF model fit with and without each innovation that all three innovations are necessary for an accurate fit. We determined that reduction enthalpy is higher and more sensitive to off-stoichiometry when the Sr fraction is large (139.5 and 185.3 kJ/mol O 2 for SrFeO 3 at δ = 0 and δ = 0.5, respectively vs. nearly constant 83 kJ/mol O 2 for BaFeO 3 ). However, the reduction entropy is mostly insensitive to Sr fraction, but highly dependent on δ suggesting larger contributions in the non-configurational entropy parameters.

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

A quantitative assessment of the model form error of friction models across different interface representations for jointed structures

Hysteretic models are widely used to model frictional interactions in joints to recreate experimental behavior. However, it is unclear which models are best suited for fitting or predicting the responses of structures. Here, the present study evaluates 26 friction model/interface representation combinations to quantify the model form error. A Quasi-Static Modal Analysis approach (termed Rayleigh Quotient Nonlinear Modal Analysis) is adopted to calculate the nonlinear system response, and a Multi-Objective Optimization is solved to fit experimental data of the first mode of the Brake-Reuß Beam. Optimized parameters from the first mode are applied to the second and third bending modes to quantify the predictive ability of the models. Formulations for both tracing full hysteresis loops and recreating hysteresis loops from a single loading curve (Masing assumptions) are considered. Smoothly varying models applied to a five patch representation showed the highest flexibility (for fitting mode 1) and good predictive potential (for modes 2 and 3). For a second formulation, which uses 152 frictional elements to represent the interface, the physically motivated spring in series with a Coulomb slip model (elastic dry friction) has high error for fitting mode 1 and performs near the middle for predicting higher modes. For both interface representation, the best fit models are not the most physical, but rather the ones with the most parameters (as expected); however, the more physical models perform somewhat better for predicting the higher modes.

42 ENGINEERING↗
DOE OSTI2018

Parameterized yields of semivolatile products from isoprene oxidation under different NOx levels: impacts of chemical aging and wall-loss of reactive gases

We developed a new parameterization for the yields of semi-volatile/intermediate volatility products from isoprene oxidation that took into account the chemical aging under different NOx levels, including functionalization and fragmentation reactions within the gas phase, oligomerization and photolysis within the condensed phase, as well as the loss of particles and gases to chamber walls. We used a box model to fit the dynamic evolution of measured secondary organic aerosol (SOA) concentrations from chamber experiments. The explicit inclusion of multigenerational gas phase chemistry and SOA oligomerization/photolysis during fitting itself had large impacts on the fitted yields that relate to the initial volatility distribution of isoprene oxidation products. During further multigenerational aging of these initial organics, we varied the branching ratio between functionalization and fragmentation. The best fits suggested increasing amounts of fragmentation with increasing initial NOx levels for mixed experiments. Using a box model to simulate evolution of SOA in the atmosphere at constant dilution rate (in absence of wall losses), we found that the yields determined by including different chemical scenarios during fitting changed predicted SOA mass concentrations substantially over aging timescale of 1-2 days. Our results demonstrated the need for better constraining smog-chamber-derived volatility distributions of both gas- and particle-phase SOA constituents with measurements.

Xing, Li↗
DOE OSTI2015

Molecular dynamics force-field refinement against quasi-elastic neutron scattering data

Quasi-elastic neutron scattering (QENS) is one of the experimental techniques of choice for probing the dynamics at length and time scales that are also in the realm of full-atom molecular dynamics (MD) simulations. This overlap enables extension of current fitting methods that use time-independent equilibrium measurements to new methods fitting against dynamics data. We present an algorithm that fits simulation-derived incoherent dynamical structure factors against QENS data probing the diffusive dynamics of the system. We showcase the difficulties inherent to this type of fitting problem, namely, the disparity between simulation and experiment environment, as well as limitations in the simulation due to incomplete sampling of phase space. We discuss a methodology to overcome these difficulties and apply it to a set of full-atom MD simulations for the purpose of refining the force-field parameter governing the activation energy of methyl rotation in the octa-methyl polyhedral oligomeric silsesquioxane molecule. Our optimal simulated activation energy agrees with the experimentally derived value up to a 5% difference, well within experimental error. We believe the method will find applicability to other types of diffusive motions and other representation of the systems such as coarse-grain models where empirical fitting is essential. In addition, the refinement method can be extended to the coherent dynamic structure factor with no additional effort.

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

Parameterization Strategies for Intermolecular Potentials for Predicting Trajectory-Based Collision Parameters

The accuracy of separable strategies for constructing full-dimensional potential energy surfaces for collisional energy transfer and collision rate calculations is studied systematically for three alcohols (A = methanol, ethanol, and butanol) and three bath gases (M = Ar, N 2 , and H 2 O). The fitting efficiency (defined as the number of ab initio data required to achieve parametrizations of a desired accuracy) is quantified for both pairwise (Buckingham or "exp6") and nonpairwise (permutationally invariant polynomials, PIPs, of Morse variables) functional forms, for four sampling strategies, and as a function of the complexity and anisotropy of the interaction potential. We find that convergence with respect to the number of sampled ab initio data is largely independent of the choice of functional form but instead varies nearly linearly with the number of adjustable parameters and depends strongly on the sampling strategy. Specifically, the use of biased Sobol quasirandom sampling is similar to 7X more efficient than using unbiased pseudorandom sampling, on average, requiring just similar to 3 computed ab initio energies per adjustable fitting parameter. The pairwise exp6 functional form is shown to provide accurate and transferable parametrizations for M = Ar but is unable to accurately describe alcohol interactions with M = N 2 and H 2 O. The nonpairwise PIP functional form, which is systematically improvable, can produce separable parametrizations with arbitrarily small fitting errors. However, these can suffer from overfitting, which is demonstrated using dynamics calculations of collision parameters for a large number of exp6 and PIP parametrizations. The tests described here validate a robust strategy for automatically generating A + M potential energy surfaces with minimal human intervention, including a quantifiable out of sample metric for judging the accuracy of the fitted surface. We further analyze this set of automatically generated potential energy surfaces to identify areas where more sophisticated fitting strategies may be desired, including pruning of the PIP expansions for large systems and improved sampling strategies more closely coupled with the description of the functional form.

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

Analysis of Individual Terrestrial Gamma‐Ray Flashes With Lightning Leader Models and Fermi Gamma‐Ray Burst Monitor Data

Abstract The Gamma‐ray Burst Monitor (GBM) onboard the Fermi spacecraft has observed many tens of sufficiently bright events, which are suitable for individual analysis. In our previous study, we fit individual, bright terrestrial gamma‐ray flashes (TGFs) with Relativistic Runaway Electron Avalanche (RREA) models for the first time. For relativistic‐feedback‐based models, the TGF‐producing electrons, which are seeded internally by a positive feedback effect, are usually accelerated in a large‐scale field with fully developed RREAs. Alternatively, lightning leader models may apply to either a large‐scale thunderstorm fields with fully developed RREAs or to inhomogeneous fields in front of lightning leaders where RREAs only develop partially. The predictions of the latter, inhomogeneous models for the TGF‐beaming geometry show some differences from estimations of the relativistic feedback models in homogeneous fields. In this work, we analyze a large sample of 66 bright Fermi GBM TGFs in the framework of lightning leader models, making comparisons with previous results from the homogeneous‐field RREA models. In most cases, the spectral analysis does not strongly favor one mechanism over the other, with 59% of the TGF events being best fit with the fully developed RREA mechanism, which corresponds to high‐potential leader models. The majority of the GBM‐measured TGFs can be best fit if the source altitude is below 15 km and 70% of events best fit by leader models cannot be satisfactorily modeled unless a tilted photon beam is used. For several spectrally soft TGFs, the tilted beam low‐potential leader model can best fit the data.

Mailyan, B. G.↗
DOE OSTI2013

Nanoparticle-Based Immunochromatographic Test Strip with Fluorescent Detector for Quantification of Phosphorylated Acetycholinesterase: An Exposure Biomarker of Organophosphorous Agents

A nanoparticle-based fluorescence immunochromatographic test strip (FITS) coupled with a hand-held detector for highly selective and sensitive detection of phosphorylated acetylcholinesterase (AChE), an exposure biomarker of organophosphate (OP) pesticides and nerve agents, is reported. In this approach, OP-AChE adducts were selectively captured by quantum dot-tagged anti-AChE antibodies (Qdot-anti-AChE) and zirconia nanoparticles (ZrO2 NPs). The sandwich-like immunoreactions were performed among the Qdot-anti-AChE, OP-AChE and ZrO2 NPs to form Qdot-anti-AChE/OP-AChE/ZrO2 complex, which was detected by recording the fluorescence intensity of Qdot captured on the test line. Paraoxon was used as the model OP pesticides. Under optimal conditions, this portable FITS immunosensor demonstrates a highly linear absorption response over the range of 0.01 nM to 10 nM OP-AChE, with a detection limit of 4 pM, coupled with a good reproducibility. Moreover, the FITS immunosensor has been validated with OP-AChE spiked human plasma samples. This is the first report on the development of ZrO2 NPs-based FITS for detection of OP-AChE adduct. The FITS immunosensor provides a sensitive and low-cost sensing platform for on-site screening/evaluating OP pesticides and nerve agents poisoning.

Acetylcholinesterase↗
DOE OSTI2025

Challenges of standard halo models in constraining galaxy properties from cosmic infrared background anisotropies

The halo model, combined with halo occupation distribution (HOD) prescriptions, is widely used to interpret cosmic infrared background (CIB) anisotropies and extract physical information about star-forming galaxies and their connection to large-scale structures. Recent CIB-specific implementations of the halo model have adopted more physical parameterizations. However, the extent to which these models can reliably recover meaningful physical parameters remains uncertain. We assessed whether the current parameterization of CIB halo models is sufficient to recover astrophysical quantities, such as star formation efficiency, η(M h , z), and halo mass at which the peak of star formation efficiency occurs, M max , when fit to mock data. We also assessed whether discrepancies arise from assumptions about galaxy emission (the HOD ingredients) or from more fundamental components in the halo model, such as bias and matter clustering. We fit the M21 CIB HOD model, implemented within the halo model framework, to mock CIB power spectra and star formation rate density (SFRD) data generated from the SIDES-Uchuu simulation, and compared the best-fit parameters to the known simulation inputs. We then repeated the analysis using a simplified version of the simulation (SSU), explicitly designed to match the HOD assumptions. A detailed comparison of model and simulation outputs was carried out to trace the origin of observed discrepancies. While the M21 HOD model provides a good fit to the mock data, it failed to recover the intrinsic parameters accurately, particularly the halo mass at which star formation efficiency peaks. This mismatch persists even when fitting data generated with the same model assumptions. We find strong agreement (within 5%) in the emission-related components (SFRD, emissivity), but observe a scale- and redshift-dependent offset exceeding 20% in the two-halo term of the CIB power spectrum. This likely arises from limitations in the treatment of halo bias and matter clustering within the linear approximation. Additionally, incorporating scatter in the SFR–halo mass relation and the spectral energy distribution (SED) templates significantly affects the shot noise (∼50%), but has only a modest impact (less than 10%) on the clustered component. These results suggest that recovering physical parameters from CIB clustering requires improvements to the cosmological ingredients of the halo model framework, such as adopting scale-dependent halo bias and nonlinear matter power spectra in addition to careful modeling of emission physics.

cosmic background radiation↗
DOE OSTI2013

Minimal SUSY SO(10) and Yukawa unification

The minimal supersymmetric (SUSY) SO(10) model, where only two Higgs multiplets {l_brace}10 Circled-Plus 126-bar{r_brace} are utilized for Yukawa couplings with matter fields, can nicely fit the neutrino oscillation parameters as well as charged fermion masses and mixing angles. In the fitting of the fermion mass matrix data, the largest element in the Yukawa coupling with the 126-bar -plet Higgs (Y{sup 126}) is found to be of order one, so that the right see-saw scale should be provided by Higgs vacuum expectation values (VEVs) of {beta}(10{sup 14}GeV). This fact causes a serious problem, namely, the gauge coupling unification is spoiled because of the presence of many exotic Higgs multiples emerging at the see-saw scale. In order to solve this problem, we consider a unification between bottom-quark and tau Yukawa couplings (b - {tau} Yukawa coupling unification) at the grand unified theory (GUT) scale, due to threshold corrections of superpartners to the Yukawa couplings at the 1 TeV scale. When the b - {tau} Yukawa coupling unification is very accurate, the largest element in Y{sub 126} can become {beta}(0.01), so that the right see-saw scale is realized by the GUT scale VEV and the usual gauge coupling unification is maintained. Since the b - {tau} Yukawa unification alters the Yukawa coupling data at the GUT scale, we re-analyze the fitting of the fermion mass matrix data by taking all the relevant free parameters into account. Unfortunately, we find that no parameter region shows up to give a nice fit for the current neutrino oscillation data and therefore, the usual picture of the gauge coupling unification cannot accommodate the fermion mass matrix data fitting in our procedure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

The eigenmode perspective of NMR spin relaxation in proteins

We developed in recent years the two-body (protein and probe) coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics from NMR spin relaxation. So far we used as descriptors the set of physical parameters that enter the SRLS model. They include the global (protein-related) diffusion tensor, D{sub 1}, the local (probe-related) diffusion tensor, D{sub 2}, and the local coupling/ordering potential, u. As common in analyzes based on mesoscopic dynamic models, these parameters have been determined with data-fitting techniques. In this study, we describe structural dynamics in terms of the eigenmodes comprising the SRLS time correlation functions (TCFs) generated by using the best-fit parameters as input to the Smoluchowski equation. An eigenmode is a weighted exponential with decay constant given by an eigenvalue of the Smoluchowski operator, and weighting factor determined by the corresponding eigenvector. Obviously, both quantities depend on the SRLS parameters as determined by the SRLS model. Unlike the set of best-fit parameters, the eigenmodes represent patterns of motion of the probe-protein system. The following new information is obtained for the typical probe, the {sup 15}N−{sup 1}H bond. Two eigenmodes, associated with the protein and the probe, dominate when the time scale separation is large (i.e., D{sub 2} ≫ D{sub 1}), the tensorial properties are simple, and the local potential is either very strong or very weak. When the potential exceeds these limits while the remaining conditions are preserved, new eigenmodes arise. The multi-exponentiality of the TCFs is associated in this case with the restricted nature of the local motion. When the time scale separation is no longer large, the rotational degrees of freedom of the protein and the probe become statistically dependent (coupled dynamically). The multi-exponentiality of the TCFs is associated in this case with the restricted nature of both the local and the global motion. The effects of local diffusion axiality, potential strength, and extent of mode-coupling on the eigenmode setup are investigated. We detect largely global motional or largely local motional eigenmodes. In addition, we detect mixed eigenmodes associated with correlated/prograde or anti-correlated/retrograde rotations of the global (D{sub 1}) and local (D{sub 2}) motional modes. The eigenmode paradigm is applied to N−H bond dynamics in the β-sheet residue K19, and the α-helix residue A34, of the third immunoglobulin-binding domain of streptococcal protein G. The largest contribution to the SRLS TCFs is made by mixed anti-correlated D{sub 1} and D{sub 2} eigenmodes. The next largest contribution is made by D{sub 1}-dominated eigenmodes. Eigenmodes dominated by the local motion contribute appreciably to A34 and marginally to K19. Correlated D{sub 1} and D{sub 2} eigenmodes contribute exclusively to K19 and do not contribute above 1% to A34. The differences between K19 and A34 are delineated and rationalized in terms of the best-fit SRLS parameters and mode-mixing. It may be concluded that eigenmode analysis is complementary and supplementary to data-fitting-based analysis.

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

Edible oil structures at low and intermediate concentrations. II. Ultra-small angle X-ray scattering of in situ tristearin solids in triolein

Ultra-small angle X-ray scattering has been used for the first time to elucidate, in situ, the aggregation structure of a model edible oil system. The three-dimensional nano- to micro-structure of tristearin solid particles in triolein solvent was investigated using 5, 10, 15, and 20% solids. Three different sample preparation procedures were investigated: two slow cooling rates of 0.5°/min, case 1 (22 days of storage at room temperature) and case 2 (no storage), and one fast cooling of 30°/min, case 3 (no storage). The length scale investigated, by using the Bonse-Hart camera at beamline ID-15D at the Advanced Photon Source, Argonne National Laboratory, covered the range from 300 Å to 10 μm. The unified fit and the Guinier-Porod models in the Irena software were used to fit the data. The former was used to fit 3 structural levels. Level 1 structures showed that the primary scatterers were essentially 2-dimensional objects for the three cases. The scatterers possessed lateral dimensions between 1000 and 4300 Å. This is consistent with the sizes of crystalline nanoplatelets present which were observed using cryo-TEM. Level 2 structures were aggregates possessing radii of gyration, R{sub g2} between 1800 Å and 12000 Å and fractal dimensions of either D{sub 2}=1 for case 3 or 1.8≤D{sub 2}≤2.1 for case 1 and case 2. D{sub 2} = 1 is consistent with unaggregated 1-dimensional objects. 1.8 ≤ D{sub 2} ≤ 2.1 is consistent with these 1-dimensional objects (below) forming structures characteristic of diffusion or reaction limited cluster-cluster aggregation. Level 3 structures showed that the spatial distribution of the level 2 structures was uniform, on the average, for case 1, with fractal dimension D{sub 3}≈3 while for case 2 and case 3 the fractal dimension was D{sub 3}≈2.2, which suggested that the large-scale distribution had not come to equilibrium. The Guinier-Porod model showed that the structures giving rise to the aggregates with a fractal dimension given by D{sub 2} in the unified fit level 2 model were cylinders described by the parameter s≈1 in the Guinier-Porod model. The size of the base of these cylinders was in agreement with the cryo-TEM observations as well as with the results of the level 1 unified fit model. By estimating the size of the nanoplatelets and understanding the structures formed via their aggregation, it will be possible to engineer novel lipids systems that embody desired functional characteristics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Multi-layer Potfit: An accurate potential representation for efficient high-dimensional quantum dynamics

The multi-layer multi-configuration time-dependent Hartree method (ML-MCTDH) is a highly efficient scheme for studying the dynamics of high-dimensional quantum systems. Its use is greatly facilitated if the Hamiltonian of the system possesses a particular structure through which the multi-dimensional matrix elements can be computed efficiently. In the field of quantum molecular dynamics, the effective interaction between the atoms is often described by potential energy surfaces (PES), and it is necessary to fit such PES into the desired structure. For high-dimensional systems, the current approaches for this fitting process either lead to fits that are too large to be practical, or their accuracy is difficult to predict and control. This article introduces multi-layer Potfit (MLPF), a novel fitting scheme that results in a PES representation in the hierarchical tensor (HT) format. The scheme is based on the hierarchical singular value decomposition, which can yield a near-optimal fit and give strict bounds for the obtained accuracy. Here, a recursive scheme for using the HT-format PES within ML-MCTDH is derived, and theoretical estimates as well as a computational example show that the use of MLPF can reduce the numerical effort for ML-MCTDH by orders of magnitude, compared to the traditionally used POTFIT representation of the PES. Moreover, it is shown that MLPF is especially beneficial for high-accuracy PES representations, and it turns out that MLPF leads to computational savings already for comparatively small systems with just four modes.

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

Constructing diabatic representations using adiabatic and approximate diabatic data – Coping with diabolical singularities

We have recently introduced a diabatization scheme, which simultaneously fits and diabatizes adiabatic ab initio electronic wave functions, Zhu and Yarkony J. Chem. Phys. 140, 024112 (2014). The algorithm uses derivative couplings in the defining equations for the diabatic Hamiltonian, H{sup d}, and fits all its matrix elements simultaneously to adiabatic state data. This procedure ultimately provides an accurate, quantifiably diabatic, representation of the adiabatic electronic structure data. However, optimizing the large number of nonlinear parameters in the basis functions and adjusting the number and kind of basis functions from which the fit is built, which provide the essential flexibility, has proved challenging. In this work, we introduce a procedure that combines adiabatic state and diabatic state data to efficiently optimize the nonlinear parameters and basis function expansion. Further, we consider using direct properties based diabatizations to initialize the fitting procedure. To address this issue, we introduce a systematic method for eliminating the debilitating (diabolical) singularities in the defining equations of properties based diabatizations. We exploit the observation that if approximate diabatic data are available, the commonly used approach of fitting each matrix element of H{sup d} individually provides a starting point (seed) from which convergence of the full H{sup d} construction algorithm is rapid. The optimization of nonlinear parameters and basis functions and the elimination of debilitating singularities are, respectively, illustrated using the 1,2,3,4{sup 1}A states of phenol and the 1,2{sup 1}A states of NH{sub 3}, states which are coupled by conical intersections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Hyperfine interactions of constituent quarks and the mass spectrum of tetraquark states

The simplest form of the Hyperfine Interaction (HFI) of quarks is a color-spin interaction, sometimes called Fermi-Breit (FB) HFI. For conventional hadrons, this HFI for quark-antiquark (qq) pairs is twice as strong as the same HFI for quark-quark (qq or qq) pairs, but this difference is absorbed into the different HFI constants for mesons and baryons, which are free parameters of the fits. For tetraquarks, however, qqHFI is only half as strong as the qq/qqHFI and both kinds of pairs coexist within the same state. Hence it is worth investigating a model with different fit constants for qq and qq/q pairs. We have investigated 3 models of spin-color HFI for light scalar tetraquarks: 2 models with single constants for both kinds of pairs and one model with different fit constants. The light scalar nonet is chosen because of a lack of other experimentally observed multiplets interpretable as pure exotic states. Since even the nonet contains only 3 particles with distinct masses (i. e. 3 fit points), the constituent quark masses had to be taken from conventional hadron fits obtained in our previous work. The two-constant model did not prove to be a viable alternative to the one-constant models. © 2019 Author(s).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Transport coefficients for magnetic-field evolution in inviscid magnetohydrodynamics

The magnetized resistivity and electrothermal tensors when substituted into the induction equation lead to electrothermal magnetic field generation, resistive magnetic diffusion, and magnetic field advection due to resistivity gradients, temperature gradients and currents. Here, the advection terms driven by temperature gradient and current have cross-field components (perpendicular to both the magnetic field and the driving term) that depend on significantly modified versions of Braginskii’s transport coefficients. The improved fits to Braginskii’s coefficients given by Epperlein and Haines and Ji and Held give physically incorrect results for cross-field advection at small Hall parameters (product of cyclotron frequency and collision time). The errors in Epperlein and Haines’ fits are particularly severe, giving increasing advection velocities below a Hall parameter of one when they should decrease linearly to zero. Epperlein and Haines’ fits can also give erroneous advection terms due to variations in effective atomic number. The only serious error in Braginskii’s fits is an overestimate in advection due to perpendicular resistivity. New fits for the cross-field advection terms are obtained from a direct numerical solution of the Fokker-Planck equation and Ji and Held’s higher order expansion approach that are continuous functions of the effective atomic number.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantifying electron temperature distributions from time-integrated x-ray emission spectra

K-shell x-ray emission spectroscopy is a standard tool used to diagnose the plasma conditions created in high-energy-density physics experiments. In the simplest approach, the emissivity-weighted average temperature of the plasma can be extracted by fitting an emission spectrum to a single temperature condition. It is known, however, that a range of plasma conditions can contribute to the measured spectra due to a combination of the evolution of the sample and spatial gradients. In this work, we define a parameterized model of the temperature distribution and use Markov Chain Monte Carlo sampling of the input parameters, yielding uncertainties in the fit parameters to assess the uniqueness of the inferred temperature distribution. Here we present the analysis of time-integrated S and Fe x-ray spectroscopic data from the Orion laser facility and demonstrate that while fitting each spectral region to a single temperature yields two different temperatures, both spectra can be fit simultaneously with a single temperature distribution. We find that fitting both spectral regions together requires a maximum temperature of $1310^{+90}_{-70}$ eV with significant contributions from temperatures down to 200 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Derivation of the Two-Exponential Probability Density Function for Rossi-Alpha Measurements of Reflected Assemblies and Validation for the Special Case of Shielded Measurements

Rossi-alpha measurements of fissionable assemblies are used to estimate the prompt neutron decay constant, α, of a system, which can be used to infer reactivity. Previous works have demonstrated that two-exponential fits for Rossi-alpha measurements of reflected assemblies are better than one exponential fits; however, the two-exponential probability density function that is needed to obtain α from the fit parameters has not been derived. This paper derives the two-region point kinetics model (two-exponential fit) for Rossi-alpha measurements, a generalization of the current, one region model (one-exponential fit). The new model is validated with measurements of reflected, fissionable assemblies (highly enriched uranium with k eff ≈ 0.95 and weapons-grade plutonium with k eff > 0.77) with α held constant. Here, the results show that the twoexponential model can (i) predict the constant α, and (ii) deconvolve α and the time a neutron spends in the reflector, neither of which is possible with the one-exponential model.

42 ENGINEERING↗
DOE OSTI2015

OPEN CLUSTERS AS PROBES OF THE GALACTIC MAGNETIC FIELD. I. CLUSTER PROPERTIES

Stars in open clusters are powerful probes of the intervening Galactic magnetic field via background starlight polarimetry because they provide constraints on the magnetic field distances. We use 2MASS photometric data for a sample of 31 clusters in the outer Galaxy for which near-IR polarimetric data were obtained to determine the cluster distances, ages, and reddenings via fitting theoretical isochrones to cluster color–magnitude diagrams. The fitting approach uses an objective χ{sup 2} minimization technique to derive the cluster properties and their uncertainties. We found the ages, distances, and reddenings for 24 of the clusters, and the distances and reddenings for 6 additional clusters that were either sparse or faint in the near-IR. The derived ranges of log(age), distance, and E(B−V) were 7.25–9.63, ∼670–6160 pc, and 0.02–1.46 mag, respectively. The distance uncertainties ranged from ∼8% to 20%. The derived parameters were compared to previous studies, and most cluster parameters agree within our uncertainties. To test the accuracy of the fitting technique, synthetic clusters with 50, 100, or 200 cluster members and a wide range of ages were fit. These tests recovered the input parameters within their uncertainties for more than 90% of the individual synthetic cluster parameters. These results indicate that the fitting technique likely provides reliable estimates of cluster properties. The distances derived will be used in an upcoming study of the Galactic magnetic field in the outer Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗