Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Local electric field measurement in GaN diodes by exciton Franz–Keldysh photocurrent spectroscopy

The eXciton Franz–Keldysh (XFK) effect is observed in GaN p–n junction diodes via the spectral variation of photocurrent responsivity data that redshift and broaden with increasing reverse bias. Photocurrent spectra are quantitatively fit over a broad photon energy range to an XFK model using only a single fit parameter that determines the line shape and the local bias (V l ), uniquely determining the local electric field maximum and depletion widths. As expected, the spectrally determined values of V l vary linearly with the applied bias (V) and reveal a large reduction in the local electric field due to electrostatic non-uniformity. The built-in bias (V bi ) is estimated by extrapolating V l at V=0, which, when compared with independent C-V measurements, indicates an overall ±0.31 V accuracy of V l . This demonstrates sub-bandgap photocurrent spectroscopy as a local probe of electric field in wide bandgap diodes that can be used to map out regions of device breakdown (hot spots) for improving electrostatic design of high-voltage devices.

42 ENGINEERING↗
DOE OSTI2020

Viscous flow properties and hydrodynamic diameter of phenothiazine-based redox-active molecules in different supporting salt environments

We report viscous flow properties of a redox-active organic molecule, N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT), a candidate for non-aqueous redox flow batteries, and two of its radical cation salts. A microfluidic viscometer enabled the use of small sample volumes in determining viscosity as a function of shear rate and concentration in the non-aqueous solvent, acetonitrile, both with and without supporting salts. All solutions tested show Newtonian behavior over shear rates of up to 30 000 s−1, which was rationalized by scaling arguments for the diffusion-based relaxation time of a single MEEPT molecule without aggregation. Neat MEEPT is flowable but with a large viscosity (412 mPa⋅s at room temperature), which is ∼1000 times larger than that of acetonitrile. MEEPT solutions in acetonitrile have low viscosities; at concentrations up to 0.5 M, the viscosity increases by less than a factor of two. From concentration-dependent viscosity measurements, molecular information was inferred from intrinsic viscosity (hydrodynamic diameter) and the Huggins coefficient (interactions). Model fit credibility was assessed using the Bayesian Information Criterion. It is found that the MEEPT and its charged cations are “flowable” and do not flocculate at concentrations up to 0.5 M. MEEPT has a hydrodynamic diameter of around 8.5 Å, which is almost insensitive to supporting salt and state of charge. This size is comparable to molecular dimensions of single molecules obtained from optimized structures using density functional theory calculations. The results suggest that MEEPT is a promising candidate for redox flow batteries in terms of its viscous flow properties.

Wang, Yilin (ORCID:0000000316102995)↗
DOE OSTI2020

Isotope effects on the high pressure viscosity of liquid water measured by differential dynamic microscopy

In this paper, differential dynamic microscopy is performed in diamond anvil cells to measure the viscosity of water along the 24 °C isotherm to high-pressure by determination of the tracer diffusion coefficient of monodisperse silica spheres of known diameter and application of the Stokes-Einstein-Sutherland equation. This technique allows liquid samples to be compressed to greater pressure prior to freezing than with other viscometry methods. The highest-pressure measurement was made at 1.67 GPa, considerably deeper into the supercompressed regime than previously reported. The effect of isotopic composition is investigated with samples of normal water, heavy water, and partially deuterated water. When data below 0.25 GPa are excluded a free volume model fits the observed viscosities well yielding a theoretical glass transition density close to that observed in very-high-density amorphous ice. The improved fit above 0.25 GPa coincides with the loss of other anomalous behaviors in liquid water caused by hydrogen bonding and represents a transition to properties closer to those of a simple liquid.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Neural network potential from bispectrum components: A case study on crystalline silicon

In this article, we present a systematic study on developing machine learning force fields (MLFFs) for crystalline silicon. While the main-stream approach of fitting a MLFF is to use a small and localized training set from molecular dynamics simulations, it is unlikely to cover the global features of the potential energy surface. Additionally, to remedy this issue, we used randomly generated symmetrical crystal structures to train a more general Si-MLFF. Furthermore, we performed substantial benchmarks among different choices of material descriptors and regression techniques on two different sets of silicon data. Our results show that neural network potential fitting with bispectrum coefficients as descriptors is a feasible method for obtaining accurate and transferable MLFFs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spin dynamics and relaxation in 7.6 nm thin film of La 0.7 Sr 0.3 MnO 3 /SrTiO 3 : ac magnetic susceptibility and magnetic viscosity investigations

Analysis of the results from our investigations of the temperature dependence of ac susceptibilities (χ' and χ'') and magnetic viscosity in a 7.6 nm thin film of La 0.7 Sr 0.3 MnO 3 grown on SrTiO 3 (001) substrate is presented. The Curie temperature (T C ) of this film is magnetic field (H) dependent, varying as [T C (H) - T C (0)] ~ H n with T C (0) = 270 K and n = 0.15. Temperature dependence of χ' and χ'' measured at ac frequencies f m in the range of 0.1–10 kHz shows a broad peak near 230 K associated with the blocking temperature (T B ) of spin clusters present in the 1.4 nm surface layer and a frequency dependent peak near 270 K associated with T C . The mean relaxation time τ determined from the Cole–Cole analysis of the temperature dependent χ' and χ'' is shown to fit the Vogel–Fulcher law: τ = τ 0 exp[ΔE/k B (T - T 0 )] with T 0 = 245 K, ΔE/k B = 270 K = T C and τ 0 = 1.2 x 10 -9 s. Magnetic viscosity S (measured in H = 0 for a sample cooled in H = 50 Oe) determined from the fit to the magnetization: M (t)= M (0) – S ln t, with time t up to 2 h, shows a peak at 230 K above which M (2 h) switches to negative values for temperatures up to T C (0). It is argued that this negative magnetization results from an interaction between the superparamagnetic spin clusters and the ferromagnetic phase, both being present in the temperature range of T B < T < T C .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

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

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vibrational spectroscopy of shock compressed condensed phase nitrous oxide: Frequency shifts in the ν 1 mode

Over many years, we have published vibrational spectra of molecular materials at extreme pressure and temperature obtained using shock compression and coherent Raman spectroscopy methods. For diatomic molecules, we were able to extract vibrational temperatures from the intensities of spectrally resolved vibrational hot bands. Larger molecules in the condensed phase suffer band broadening effects that obscure the vibrational hot bands as the anharmonicities are typically smaller than the widths of the bands. This inability to resolve the hot bands inhibits the extraction of vibrational temperature and the measurement of the vibrational frequency of the fundamental. Here, we use a hot band model based on gas-phase anharmonic coupling coefficients to fit coherent anti-Stokes Raman spectra of the ν 1 vibrational mode of shock compressed condensed phase N 2 O with shock pressures and temperatures estimated from literature equations of state and compare to fits from a model using a single Gaussian peak. We report the resulting vibrational frequency shifts with shock pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Response of CR-39 nuclear track detectors to protons with non-normal incidence

This paper presents data from experiments with protons at non-normal incidence to CR-39 nuclear track detectors, analyzing the properties of detection efficiency, proton track diameter, track contrast, and track eccentricity. Understanding the CR-39 response to protons incident at an angle is important for designing charged particle detectors for inertial confinement fusion (ICF) applications. This study considers protons with incident energies less than 3 MeV. In this regime, an incident angle of 10° has no effect on CR-39 detection efficiency, and >85% detection efficiency is preserved up through 25° in the range of 1.0 MeV–2.1 MeV. For ICF applications, incident angles above 30° are deemed impractical for detector design due to significant drops in proton detection at all energies. We observe significant reductions in detection efficiency compared to theoretical predictions, particularly at low energies where proton tracks are etched away. The proton track diameter measured by the scan system is observed to decrease with higher incident angles. The track diameters are analyzed with two fitting models, and it is shown that the diameter–energy relation can be fit with the existing models at angles up to 30° . The optical contrast of the tracks tends to increase with the angle, meaning that the tracks are fainter, and a larger increase is observed for higher energies. Eccentricity, a measure of how elongated proton tracks are, increases with the incident angle and drops after the critical angle. The lowest energy tracks remain nearly circular even at higher angles.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Magnetic and electronic structures of antiferromagnetic topological material candidate EuMg 2 Bi 2

In this work, EuMg 2 Bi 2 has been investigated to understand the electronic and magnetic behaviors as an antiferromagnetic (AFM) topological semimetal candidate. High-quality single crystals of EuMg 2 Bi 2 were grown via a Bi flux and, subsequently, characterized to be consistent with the previously reported bulk magnetic and resistivity properties. A ferromagnetic interaction is indicated by the positive Curie–Weiss temperature obtained through fitting the bulk magnetic susceptibility data. The bulk resistivity measurements reveal an interesting electronic behavior that is potentially influenced by a competing antiferromagnetic and ferromagnetic interaction in and out of the ab plane. From the resulting refinement of the neutron diffraction data, EuMg 2 Bi 2 was found to exhibit an A-type magnetic structure with Eu 2+ moments ferromagnetically aligned in the plane and antiferromagnetically stacked between neighbor ferromagnetic Eu layers. The power law fitting magnetic ordering parameter below TN~8 K agrees with the 2D Heisenberg model, indicating a weak interlayer antiferromagnetic interaction. Considering the magnetic structure determined by neutron diffraction, the surface state calculation suggests that EuMg 2 Bi 2 is an AFM topological insulator candidate. Linearly dispersed Dirac surface states were also observed in our angle-resolved photoemission spectroscopy measurements, consistent with the calculation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On wavelength-dependent exciton lifetime distributions in reconstituted CP29 antenna of the photosystem II and its site-directed mutants

To provide more insight on the excitonic structure and exciton lifetimes of the wild type (WT) CP29 complex of Photosystem II (PSII) we measured high-resolution (low temperature) absorption, emission, and hole burned spectra for the A2 and B3 mutants, which lack chlorophylls a612 and b614 (Chls), respectively. Experimental and modeling results obtained for the WT CP29 and A2/B3 mutants provide new insight on the mutation-induced changes at the molecular level and shed more light on energy transfer dynamics. Simulations of the A2 and B3 optical spectra, using the second-order non-Markovian theory, and comparison with improved fits of WT CP29 optical spectra, provide more insight on their excitonic structure, mutation induced changes, and frequency-dependent distributions of exciton lifetimes (T 1 ). A new Hamiltonian obtained for WT CP29 reveals that deletion of Chls a612 or b614 induces changes in the site-energies of all remaining Chls. Hamiltonians obtained for A2 and B3 mutants are discussed in the context of the energy landscape of chlorophylls, excitonic structure and the transfer kinetics. Our data suggest that the lowest exciton states in A2 and B3 mutants are contributed by a611(57%), a610(17%), a615(15%) and a615(58%), a611(20%), a612(15%) Chls, respectively, although other compositions of lowest energy states are also discussed. Lastly, we argue that the calculated exciton decay times are consistent with both the hole-burning and recent transient absorption measurements. Wavelength-dependent T 1 distributions offer more insight into the interpretation of kinetic traces commonly described by discrete exponentials in global analysis/global fitting of transient absorption experiments.

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

Collective optical Thomson scattering in pulsed-power driven high energy density physics experiments (invited)

Optical collective Thomson scattering is used to diagnose magnetized high energy density physics experiments at the Magpie pulsed-power generator at Imperial College London. The system uses an amplified pulse from the 2nd harmonic of a Nd:YAG laser (3 J, 8 ns, 532 nm) to probe a wide diversity of hightemperature plasma objects; with densities in the range of 10 17 -10 19 cm -3 and temperatures between 10 eV and a few keV. The scattered light is collected from 100 μm-scale volumes within the plasmas, which are imaged onto optical fiber arrays. Multiple collection systems observe these volumes from different directions, providing simultaneous probing with different scattering K-vectors (and different associated α-parameters, typically in the range 0.5 – 3) allowing independent measurements of separate velocity components of the bulk plasma flow. The fiber arrays are coupled to an imaging spectrometer with a gated ICCD. The spectrometer is configured to view the ion-acoustic waves (IAWs) of the collective Thomson scattered spectrum. Fits to the spectra with the theoretical spectral density function S(K,ω) yield measurements of the local plasma temperatures and velocities. Fitting is constrained by independent measurements of the electron density from laser interferometry, and the corresponding spectra for different scattering vectors. Furthermore, this TS diagnostic has been successfully implemented on a wide range of experiments, revealing temperature and flow velocity transitions across magnetized shocks, inside rotating plasma jets and imploding wire arrays, as well as providing direct measurements of drift velocities inside a magnetic reconnection current sheet.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observations of a plectonemic configuration in a stable magnetized plasma jet

Astrophysical jets are collimated high-speed outflows emerging from spinning and accreting matter around celestial objects and may spontaneously result from self-organized processes. Magnetic self-organization is commonly observed in laboratory plasma physics experiments; however, they require close-fitting flux conservers to constrain and stabilize the toroidal or cylindrical structures. Here we report the first observations of a long, stable, free-boundary plasma jet far from chamber walls, embedding a double-helix magnetic structure resembling a force-free plectonemic Taylor state. The jets arise from an experimental setup that mimics an accretion disk and has no close-fitting solid flux conserver. The results support the hypothesis that self-organization could be a universal, intrinsic explanation for jet formation, collimation, and stability and may help explain double-helix features in celestial observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Free space Thomson scattering to study high energy density shocks

A free space collective Thomson scattering system has been developed to study pulsed power produced plasmas. While most Thomson scattering diagnostics on pulsed power machines use a bundle of fibers to couple scattered light from the plasma to the spectrometer, this system used free space coupling of the light, which enabled a spatially continuous image of the plasma. Initial experiments with this diagnostic were performed on an inverse wire array generated by a 200 kA, 1100 ns rise time pulse power generator. The capabilities of this diagnostic were demonstrated by using the low frequency ion acoustic wave feature of the Thomson scattering spectra to measure the plasma flow velocity. The diagnostic was demonstrated to measure velocities between 20 and 40 km/s with an error of 4.7 km/s when fitting with a 600 μm spatial resolution or 8.9 km/s when fitting with a 150 μm spatial resolution. In some experiments, the diagnostic observed a bow shock in the plasma flow as the scattering intensity increased and flow velocity decreased.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

X-ray multi-probe data acquisition: A novel technique for laser pump x-ray transient absorption spectroscopy

We report the development and implementation of a novel data acquisition (DAQ) technique for synchrotron-based laser pump X-ray Transient Absorption (XTA) spectroscopy, called X-ray Multi-Probe DAQ (XMP DAQ). This technique utilizes high performance analog to digital converters and home-built software to efficiently measure and process the XTA signal from all x-ray pulses between laser excitations. XMP DAQ generates a set of time resolved x-ray absorption spectra at thousands of different pump–probe time delays simultaneously. Two distinct XMP DAQ schemes are deployed to accommodate different synchrotron storage ring filling patterns. Current Integration (CI) DAQ is a quasi-analog technique that implements a fitting procedure to extract the time resolved absorption intensity from the averaged fluorescence detector response. The fitting procedure eliminates issues associated with small drifts in the voltage baseline and greatly enhances the accuracy of the technique. Photon Counting (PC) DAQ is a binary technique that uses a time resolved histogram to calculate the XTA spectrum. While PC DAQ is suited to measure XTA data with closely spaced x-ray pulses (~10 ns) and a low count rate (<1 detected photon/pulse), CI DAQ works best for widely spaced pulses (tens of ns or greater) with a high count rate (>1 detected photon/pulse). XMP DAQ produces a two-dimensional XTA dataset, enabling efficient quantitative analysis of photophysical and photochemical processes from the sub-nanosecond timescale to 100 μs and longer.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

A dimensionless number for high-throughput design of multi-principal element alloys in directed energy deposition

The present study demonstrates the development of a dimensionless number to predict the build height in the additive manufacturing technique of directed energy deposition (DED). The build height can also be used to estimate the dendrite arm spacing and, thus, the cooling rate in the fabrication of samples. A baseline sample, 316L stainless steel, was used to fit the build height to the dimensionless number. A range of process parameters, including laser power, laser feed rate, powder flow rate, layer thickness, and hatch spacing, were varied. Based upon dendrite arm spacing, the estimated cooling rate varied between 102 and 104 K/s. Using the fitted relationship for the stainless steel, high-throughput (HT) processing of multi-principal element alloys (MPEAs) was performed. For this study, HT is the ability to fabricate a batch of 25 bulk samples (∼1 cm3) with different compositions within a 5-h period with ±10 at. % accuracy. A range of compositions using in situ alloying of elemental powders in the Fe–Ni–Cr–Mo system were made. The MPEAs' build height followed the same relationship to the dimensionless number as the 316L alloy. The dimensionless number predicts both macro and meso-scale features in HT processing, thus offering a design tool for choosing process parameters in DED additive manufacturing. Also, the ability to control or increase cooling rates can enhance the ability to promote metastability as well as control meso-scale chemical distributions of alloy samples.

Islam, Zahabul (ORCID:0000000169613209)↗