Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dynamics of poly[n]catenane melts

Inspired by advances in the chemical synthesis of interlocking polymer architectures, extensive molecular dynamics simulations have been conducted to study the dynamical properties of poly[n]catenanes-polymers composed entirely of interlocking rings-in the melt state. Both the degree of polymerization (number of links) and the number of beads per ring are systematically varied, and the results are compared to linear and ring polymers. Herein, a simple Rouse-like model is presented, and its analytical solution suggests a decomposition of the dynamics into "ring-like" and "linear-like" regimes at short and long times, respectively. In agreement with this picture, multiple sub-diffusive regimes are observed in the monomer mean-squared-displacements even though interchain entanglement is not prevalent in the system. However, the Rouse-type model does not account for the topological effects of the mechanical bonds, which significantly alter the dynamics at intermediate length scales both within the rings and at the chain segment scales. The stress relaxation in the system is extremely rapid and may be conveniently separated into ring-like and linear-like contributions, again in agreement with the Rouse picture. However, the viscosity has a non-monotonic dependence on the ring size for long chains, which disagrees strongly with theoretical predictions. This unexpected observation cannot be explained in terms of chain disentanglement and is inconsistent with other measures of polymer relaxation. Possible mechanisms for this behavior are proposed and implications for materials design are discussed.

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

Microscopic model of the doping dependence of linewidths in monolayer transition metal dichalcogenides

A fully microscopic model of the doping-dependent exciton and trion linewidths in the absorption spectra of monolayer transition metal dichalcogenides in the low temperature and low-doping regime is explored. The approach is based on perturbation theory and avoids the use of phenomenological parameters. In the low-doping regime, we find that the trion linewidth is relatively insensitive to doping levels, while the exciton linewidth increases monotonically with doping. On the other hand, we argue that the trion linewidth shows a somewhat stronger temperature dependence. The magnitudes of the linewidths are likely to be masked by phonon scattering for T ≥ 20 K in encapsulated samples in the low-doping regime. We discuss the breakdown of perturbation theory, which should occur at relatively low-doping levels and low temperatures. Our work also paves the way toward understanding a variety of related scattering processes, including impact ionization and Auger scattering in clean 2D samples

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

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

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

Pedestal magnetic turbulence measurements in ELMy H-mode DIII-D plasmas by Faraday-effect polarimetry

Internal magnetic fluctuation measurements are utilized to identify turbulence associated with micro-tearing modes (MTM) in the DIII-D Edge-Localized-Mode (ELM)-y H-mode pedestal. Using a Faraday-effect polarimeter, magnetic turbulence (150–500 kHz) is directly observed with a typical line-averaged fluctuation amplitude of ~0.8 G at peak frequency (250 kHz) and ~15 G integrated over the spectrum from 150 to 500 kHz. Frequency, poloidal wavenumber, and propagation direction of the magnetic turbulence all serve to identify as MTM. Magnetic turbulence amplitude non-monotonically correlates with collision frequency, peaks off mid-plane, and correlates with electron temperature gradient evolution between ELMs, consistent with MTM features identified from theory and gyro-kinetic simulation. The magnetic turbulence growth correlates with confinement degradation in ELMy H-mode plasmas during a slow density ramp. In conclusion, these internal measurements provide unique constraints toward developing physics understanding and validating models of the H-mode pedestal for future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure—An experimental and theoretical study

The 57Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α→ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth’s inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe2O3. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

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

Gold-tantalum alloy films deposited by high-density-plasma magnetron sputtering

Gold-tantalum alloy films are of interest for biomedical and magnetically-assisted inertial confinement fusion applications. Here, we systematically study the effects of substrate tilt (0°–80°) and negative substrate bias (0–100 V) on properties of ≲3-μm-thick films deposited by high-power impulse magnetron sputtering (HiPIMS) from a Au–Ta alloy target (with 80 at. % of Ta). Results reveal that, for all the substrate bias values studied, an increase in substrate tilt leads to a monotonic decrease in film thickness, density, residual compressive stress, and electrical conductivity. Larger substrate bias favors the formation of a body-centered cubic phase, with films exhibiting lower column tilt and higher density, electrical conductivity, and residual compressive stress. Furthermore, these changes are attributed to metal atom ionization effects, based on the lack of correlation with distributions of landing energies and incident angles of depositing species as calculated by Monte Carlo simulations of ballistic collisions and gas phase atomic transport. By varying substrate tilt and bias in HiPIMS deposition, properties of Au–Ta alloy films can be controlled in a very wide range, including residual stress from –2 to +0.5 GPa, density from 12 to 17 g/cm 3 , and the electrical resistivity from 50 to 4500 μΩ cm, enabling optimum deposition conditions to be selected for specific applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mode spectrum characteristics and onset of the low-shear MHD stability regime

Equilibria with extended regions of weak magnetic shear, including some tokamak scenarios and stellarators, can be susceptible to pressure-driven internal MHD instabilities even though there is no mode rational surface in the plasma. Nonresonant modes, in particular, can have properties that are unattractive for confinement, including displacing substantial volumes of the plasma and leading to more efficient pressure gradient flattening in the nonlinear regime. The onset and linear properties of the low shear stability regime are examined using the initial-value, extended-MHD code M3D-C. For monotonic q-profiles, we demonstrate a clear correlation between the convergents associated with the continued fraction representation of q 0 and the spectrum of unstable modes. Nonresonant modes are observed to be destabilized preferentially to any other resonant instability with the same toroidal mode number when n > 1. Using the observed connection between the spectrum associated with q 0 and the overall equilibrium stability properties, we suggest a technique for reducing the uncertainty on both q 0 and magnetic shear in the core region, obtained either via measurement or through the analysis and reconstruction of experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A theoretical approach for transient shock strengthening in high-energy-density laser compression experiments

In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dynamic strain evolution in an optically excited Pt thin film

The structural evolution of a Pt thin film following photo-thermal excitation by 1 ps optical laser pulses was studied with a time resolution of 100 ps over a total time period of 1 ms. Laser pulse fluences below 50 mJ/cm 2 were insufficient to relax the residual stress state of the as-prepared film even after 10 000 pulses. In this fluence regime, a rapid initial lattice expansion and a decrease in the lattice coherence length due to ultrafast photo-thermal heating were observed. The lattice expansion reached a maximum, and the coherence length reached a minimum, 100–200 ps after excitation before monotonically decaying back to their initial values in about 1 µs. Laser pulse fluences greater than 50 mJ/cm 2 produced irreversible stress relaxation within the first 10 optical pulses. In this regime, the lattice expansion was qualitatively similar to that in the low fluence regime, except that the initial structural state was not recovered. The evolution in the coherence length, however, was more complex. Following an initial decrease similar to that observed at low fluence, the coherence length then increased to a broad maximum greater than the initial value, before recovery.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effects of strain, disorder, and Coulomb screening on free-carrier mobility in doped cadmium oxide

Here, the interplay of stress, disorder, and Coulomb screening dictating the mobility of doped cadmium oxide (CdO) is examined using Raman spectroscopy to identify the mechanisms driving dopant incorporation and scattering within this emerging infrared optical material. Specifically, multi-wavelength Raman and UV-vis spectroscopies are combined with electrical Hall measurements on a series of yttrium (X = Y) and indium (X = In) doped X:CdO thin-films. Hall measurements confirm n-type doping and establish carrier concentrations and mobilities. Spectral fitting along the low-frequency Raman combination bands, especially the TA+TO(X) mode, reveals that the evolution of strain and disorder within the lattice as a function of dopant concentration is strongly correlated with mobility. Coupling between the electronic and lattice environments was examined through analysis of first- and second-order longitudinal–optical phonon–plasmon coupled modes that monotonically decrease in energy and asymmetrically broaden with increasing dopant concentration. By fitting these trends to an impurity-induced Fröhlich model for the Raman scattering intensity, exciton–phonon and exciton–impurity coupling factors are quantified. These coupling factors indicate a continual decrease in the amount of ionized impurity scattering with increasing dopant concentration and are not as well correlated with mobility. This shows that lattice strain and disorder are the primary determining factors for mobility in donor-doped CdO. In aggregate, the study confirms previously postulated defect equilibrium arguments for dopant incorporation in CdO while at the same time identifying paths for its further refinement.

42 ENGINEERING↗
DOE OSTI2022

A robust three-parameter reference curve for condensed phase materials

Many isothermal equations of state lead to numerical problems under tension conditions due to negative bulk modulus. We summarize the properties of four existing models and devise a continuous non-piecewise reference curve to address this issue. The final functional form contains three fitting parameters that are directly related to the bulk modulus and its pressure derivative at the reference state and at infinite pressures. We also include bounds that guarantee the monotonicity of its pressure derivative in all volume space. To demonstrate the qualitative features of each isotherm, we conduct calibrations against experimental data for cyclotetramethylene tetranitramine, magnesium oxide, gold, and platinum. Furthermore, pressure residuals and root mean square deviations show that the proposed functional form is capable of adequately fitting the data compared to other three-parameter models.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
DOE OSTI2022

The effect of polymer stiffness on magnetization reversal of magnetorheological elastomers

Ultrasoft magnetorheological elastomers (MREs) offer convenient real-time magnetic field control of mechanical properties that provides a means to mimic mechanical cues and regulators of cells in vitro. Here, we systematically investigate the effect of polymer stiffness on magnetization reversal of MREs using a combination of magnetometry measurements and computational modeling. Poly-dimethylsiloxane-based MREs with Young’s moduli that range over two orders of magnitude were synthesized using commercial polymers Sylgard™ 527, Sylgard 184, and carbonyl iron powder. The magnetic hysteresis loops of the softer MREs exhibit a characteristic pinched loop shape with almost zero remanence and loop widening at intermediate fields that monotonically decreases with increasing polymer stiffness. A simple two-dipole model that incorporates magneto-mechanical coupling not only confirms that micrometer-scale particle motion along the applied magnetic field direction plays a defining role in the magnetic hysteresis of ultrasoft MREs but also reproduces the observed loop shapes and widening trends for MREs with varying polymer stiffnesses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Metriplectic foundations of gyrokinetic Vlasov–Maxwell–Landau theory

Here, this Letter reports on a metriplectic formulation of a collisional, nonlinear full-f electromagnetic gyrokinetic theory compliant with energy conservation and monotonic entropy production. In an axisymmetric background magnetic field, the toroidal angular momentum is also conserved. Notably, a new collisional current, contributing to the gyrokinetic Maxwell–Ampère equation and the gyrokinetic charge conservation law, is discovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermal activation of low-density Ga implanted in Ge

The nuclear spins of low-density implanted Ga atoms in Ge are interesting candidates for solid state-based qubits. To date, activation studies of implanted Ga in Ge have focused on high densities. Here, we extend activation studies into the low-density regime. We use spreading resistance profiling and secondary ion mass spectrometry to derive electrical activation of Ga ions implanted into Ge as a function of the rapid thermal anneal temperature and implant density. We show that for our implant conditions, the activation is best for anneal temperatures between 400 and 650 °C with a maximum activation of 69% at the highest fluence. Below 400 °C, remaining implant damage results in defects that act as superfluous carriers, and above 650 °C, surface roughening and loss of Ga ions are observed. The activation increased monotonically from 10% to 69% as the implant fluence increased from 6×1010 to 6×1012 cm−2. The results provide thermal anneal conditions to be used for initial studies of using low-density Ga atoms in Ge as nuclear spin qubits.

Foster, Natalie D. (ORCID:0000000235712054)↗
DOE OSTI2022

Diagnosing the pedestal magnetic field and magnetohydrodynamics radial structure with pedestal–scrape of layer electron cyclotron emission radiation inversion in H-mode plasma (invited)

Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Linear stability analysis via simulated annealing and accelerated relaxation

Simulated annealing (SA) is a kind of relaxation method for finding equilibria of Hamiltonian systems. A set of evolution equations is solved with SA, which is derived from the original Hamiltonian system so that the energy of the system changes monotonically while preserving Casimir invariants inherent to noncanonical Hamiltonian systems. The energy extremum reached by SA is an equilibrium. Since SA searches for an energy extremum, it can also be used for stability analysis when initiated from a state where a perturbation is added to an equilibrium. The procedure of the stability analysis is explained, and some examples are shown. Because the time evolution is computationally time consuming, efficient relaxation is necessary for SA to be practically useful. An acceleration method is developed by introducing time dependence in the symmetric kernel used in the double bracket, which is part of the SA formulation described here. An explicit formulation for low-beta reduced magnetohydrodynamics (MHD) in cylindrical geometry is presented. In conclusion, since SA for low-beta reduced MHD has two advection fields that relax, it is important to balance the orders of magnitude of these advection fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗