Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 415 records · Page 23

Hidden sector monopole dark matter with matter domination

The thermal freeze-out mechanism for relic dark matter heavier than O (10 – 100 TeV) requires cross-sections that violate perturbative unitarity. Yet the existence of dark matter heavier than these scales is certainly plausible from a particle physics perspective, pointing to the need for a non-thermal cosmological history for such theories. Topological dark matter is a well-motivated scenario of this kind. Here the hidden-sector dark matter can be produced in abundance through the Kibble-Zurek mechanism describing the non-equilibrium dynamics of defects produced in a second order phase transition. We revisit the original topological dark matter scenario, focusing on hidden-sector magnetic monopoles, and consider more general cosmological histories. We find that a monopole mass of order (1–10 5 ) PeV is generic for the thermal histories considered here, if monopoles are to entirely reproduce the current abundance of dark matter. In particular, in a scenario involving an early era of matter domination, the monopole number density is always less than or equal to that in a pure radiation dominated equivalent provided a certain condition on critical exponents is satisfied. This results in a larger monopole mass needed to account for a fixed relic abundance in such cosmologies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigation of desorption of hydrogen gas from polymer matrix using thermal desorption analysis and finite element modeling

Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.

36 MATERIALS SCIENCE↗

Topology optimization of compliant mechanisms under transient thermal conditions

This work considers multi-material topology optimization of compliant mechanisms under transient thermal and quasi-static mechanical conditions wherein thermally actuated devices are optimized for different operating conditions. The materials are modeled using finite strain thermo-hyperelasticity and a two way coupling between the energy balance and equilibrium equations is investigated. The design updates are generated from the gradient-based method of moving asymptotes optimizer and the sensitivities are computed using the time dependent adjoint sensitivity analysis. Results show the impact of designing for short versus long actuation times.

42 ENGINEERING↗

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗

Quantifying Energy Conversion in Higher-Order Phase Space Density Moments in Plasmas

Weakly collisional and collisionless plasmas are typically far from local thermodynamic equilibrium (LTE), and understanding energy conversion in such systems is a forefront research problem. The standard approach is to investigate changes in internal (thermal) energy and density, but this omits energy conversion that changes any higher-order moments of the phase space density. In this Letter, we calculate from first principles the energy conversion associated with all higher moments of the phase space density for systems not in LTE. Particle-in-cell simulations of collisionless magnetic reconnection reveal that energy conversion associated with higher-order moments can be locally significant. Furthermore, the results may be useful in numerous plasma settings, such as reconnection, turbulence, shocks, and wave-particle interactions in heliospheric, planetary, and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Competing Easy-Axis Anisotropies Impacting Magnetic Tunnel Junction-Based Molecular Spintronics Devices (MTJMSDs)

Molecular spintronics devices (MSDs) attempt to harness molecules’ quantum state, size, and configurable attributes for application in computer devices—a quest that began more than 70 years ago. In the vast number of theoretical studies and limited experimental attempts, MSDs have been found to be suitable for application in memory devices and futuristic quantum computers. MSDs have recently also exhibited intriguing spin photovoltaic-like phenomena, signaling their potential application in cost-effective and novel solar cell technologies. The molecular spintronics field’s major challenge is the lack of mass-fabrication methods producing robust magnetic molecule connections with magnetic electrodes of different anisotropies. Another main challenge is the limitations of conventional theoretical methods for understanding experimental results and designing new devices. Magnetic tunnel junction-based molecular spintronics devices (MTJMSDs) are designed by covalently connecting paramagnetic molecules across an insulating tunneling barrier. The insulating tunneling barrier serves as a mechanical spacer between two ferromagnetic (FM) electrodes of tailorable magnetic anisotropies to allow molecules to undergo many intriguing phenomena. Our experimental studies showed that the paramagnetic molecules could produce strong antiferromagnetic coupling between two FM electrodes, leading to a dramatic large-scale impact on the magnetic electrode itself. Recently, we showed that the Monte Carlo Simulation (MCS) was effective in providing plausible insights into the observation of unusual magnetic domains based on the role of single easy-axis magnetic anisotropy. Here, we experimentally show that the response of a paramagnetic molecule is dramatically different when connected to FM electrodes of different easy-axis anisotropies. Motivated by our experimental studies, here, we report on an MCS study investigating the impact of the simultaneous presence of two easy-axis anisotropies on MTJMSD equilibrium properties. In-plane easy-axis anisotropy produced multiple magnetic phases of opposite spins. The multiple magnetic phases vanished at higher thermal energy, but the MTJMSD still maintained a higher magnetic moment because of anisotropy. The out-of-plane easy-axis anisotropy caused a dominant magnetic phase in the FM electrode rather than multiple magnetic phases. The simultaneous application of equal-magnitude in-plane and out-of-plane easy-axis anisotropies on the same electrode negated the anisotropy effect. Our experimental and MCS study provides insights for designing and understanding new spintronics-based devices.

42 ENGINEERING↗

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

H- and Dissociation in Ultra-hot Jupiters: A Retrieval Case Study of WASP-18b

Atmospheres of a number of ultra-hot Jupiters (UHJs) with temperatures ≳2000 K have been observed recently. Many of these planets show largely featureless thermal spectra in the near-infrared observed with the HST WFC3 spectrograph (1.1–1.7 μm) even though this spectral range contains strong H{sub 2}O opacity. Recent works have proposed the possibility of H- opacity masking the H{sub 2}O feature and/or thermal dissociation of H{sub 2}O causing its apparent depletion at the high temperatures of UHJs. In this work, we test these hypotheses using observations of the exoplanet WASP-18b as a case study. We report detailed atmospheric retrievals of the planet using the HyDRA retrieval code, extended to include the effects of H- opacity and thermal dissociation. We report constraints on the H{sub 2}O, CO, and H- abundances as well as the pressure–temperature profile of the dayside atmosphere for retrievals with and without H-/dissociation for each data set. We find that the H{sub 2}O and H- abundances are relatively unconstrained given the featureless WFC3 spectra. We do not conclusively detect H- in the planet, contrary to previous studies that used equilibrium models to infer its presence. The constraint on the CO abundance depends on the combination of WFC3 and Spitzer data, ranging from solar to super-solar CO values. We additionally see signs of a thermal inversion from two of the data sets. Our study demonstrates the potential of atmospheric retrievals of UHJs, including the effects of H- and thermal dissociation of molecules.

79 ASTRONOMY AND ASTROPHYSICS↗

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗

How does negative triangularity mitigate ITG turbulence and transport?

Improved confinement in negative triangularity (NT) experiments is attributed to reduced fluxes driven by micro-turbulence. The physical mechanism of why thermal confinement improves in NT relative to PT is unknown. This study employs gyrokinetic flux tube simulations using the GENE code with local Miller equilibrium to elucidate the physical mechanisms behind the beneficial effects of NT flux surface shapes. The focus is on collisionless ion temperature gradient (ITG) driven turbulence with adiabatic electrons. The kinetic profiles are held fixed across a scan of triangularity values, thus enabling comparisons on a level playing field. The reduced linear growth rates for NT is shown to be due to a reduced eigenmode averaged magnetic drift frequency and a wider, stronger negative local magnetic shear region about the outboard mid-plane. The nonlinear heat flux is lower for NT than that for PT, due to reduced radial correlation length and increased correlation time (τ c ) of fluctuations. These, in turn, are due to a comparatively higher level of self-generated zero-frequency E × B zonal shearing rate ω E in NT as compared to PT. Though the linear zonal potential residual is lower for NT, the nonlinearly generated E × B zonal shearing rate is higher for NT than for PT. This outcome is linked to the distinctive features of the radial wavenumber spectra of the zonal potential and the zonal shearing rate. The dimensionless parameter ω E τ c is suggested as a figure of merit. This is higher for NT than for PT. Thus, the reduced heat diffusivity for NT is linked to increased ω E τ c . Self-generated temperature corrugations (i.e. zonal temperature gradients) are much weaker than the background mean temperature gradient. Nevertheless, temperature corrugations are more pronounced in NT than in PT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tuning Epoxy Thermomechanics via Thermal Isomerization: A Route to Negative Coefficient of Thermal Expansion Materials

Fine control over the thermal expansion and contraction behavior of polymer materials is challenging. Most polymers have large coefficients of thermal expansion (CTEs), which preclude long performance lifetimes of composite materials. Herein, we report the design and synthesis of epoxy thermosets with low CTE values below their T g and large contraction behavior above T g by incorporating thermally contractile dibenzocyclooctane (DBCO) motifs within the thermoset network. This atypical thermomechanical behavior was rationalized in terms of a twist-boat to chair conformational equilibrium of the DBCO linkages. We anticipate these findings to be generally useful in the preparation of materials with designed CTE values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Models of resistive wall tearing mode disruptions

Disruptions are a serious issue in tokamaks. In a disruption, the thermal energy is lost by means of an instability which could be a resistive wall tearing mode (RWTM). During precursors to a disruption, the plasma edge region cools, causing the current to contract. Model sequences of contracted current equilibria are given, and their stability is calculated. A linear stability study shows that there is a maximum value of edge q a ≈ 3 for RWTMs to occur. Further, this also implies a minimum rational surface radius normalized to plasma radius from RWTMs to be unstable. Nonlinear simulations are performed using a similar model sequence derived from an equilibrium reconstruction. There is a striking difference in the results, depending on whether the wall is ideal or resistive. With an ideal wall, the perturbations saturate at moderate amplitude, causing a minor disruption without a thermal quench. With a resistive wall, there is a major disruption with a thermal quench, if the edge q a ≤ 3. There is a sharp transition in nonlinear behavior at q a = 3. This is consistent with the linear model and with experiments. If disruptions are caused by RWTMs, then devices with highly conducting walls, such as the International Tokamak Experimental Reactor will experience much milder, tolerable, disruptions than presently predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal conductivity and theory of inelastic scattering of phonons by collective fluctuations

We study the intrinsic scattering of phonons by a general quantum degree of freedom, i.e., a fluctuating “field” Q, which may have completely general correlations, restricted only by unitarity and translational invariance. From the induced scattering rates, generalizing the model studied in a companion paper [Mangeolle et al., Phys. Rev. X 12, 041031 (2022)], we obtain the consequences on the thermal conductivity tensor of the phonons. We confirm that, even within our generalized model, the off-diagonal scattering rates involve a minimum of three- to four-point correlation functions of the Q fields, and we discuss the “semiclassical” versus “quantum” nature of all contributions. We obtain general and explicit forms for these correlations which isolate the contributions to the Hall conductivity, and we provide a general discussion of the implications of symmetry and equilibrium; this elaborates on, and extends, the results of Mangeolle et al. [Phys. Rev. X 12, 041031 (2022)]. We also extend the discussion and evaluation of these two- (diagonal scattering) and four-point correlation functions, and hence the thermal transport, for the illustrative example of an ordered two-dimensional antiferromagnet, where the Q field is a composite of magnon operators arising from spin-lattice coupling, and confirm numerically that the results, while satisfying all the necessary symmetry restrictions, lead to nonvanishing scattering and Hall effects. In particular, we investigate, both analytically and numerically, the dependence of such intrinsic scattering on a crucial parameter—the magnon to phonon velocity ratio v. We in particular confirm that within some range of v of order 1 the skew-scattering mechanism leads to comparable thermal Hall conductivity for thermal currents within and normal to the plane of the antiferromagnetism, and we discover that the temperature scaling of the longitudinal conductivity displays a threshold effect and a nonuniversal, continuous variation of the scaling exponent with v.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of the 239 Pu$(n, f$) prompt fission neutron spectrum from 10 keV to 10 MeV induced by neutrons of energy 1–20 MeV

Although the prompt fission neutron spectrum (PFNS) is an essential component of neutron-driven systems that has been measured for decades, there are still multiple glaring unknowns regarding the PFNS of major actinides in the fission neutron incident energy range, specifically with regard to multichance fission and pre-equilibrium neutron emission processes. The only impactful experimental 239 Pu PFNS measurements included in recent nuclear data evaluations were measured over a limited outgoing neutron energy range at thermal and 1.5-MeV average incident neutron energy, while other potentially impactful measurements have been shown to contain errors that resulted in either large uncertainty increases or in complete exclusion from nuclear data evaluation. Here, we report a measurement of the 239 Pu PFNS over a wide range of incident neutron energy (1–20 MeV) and three orders of magnitude in outgoing neutron energy (0.01–10 MeV) resulting from the Chi-Nu experiment at the Los Alamos Neutron Science Center. These results are the combination of separate PFNS measurements in the same experimental area, one using a Li-glass and the other a liquid scintillator detector array. Covariances between all PFNS data points from each detector and within each incident energy range were generated between all other data in both detector arrays and within all other incident neutron energy bins, yielding a single covariance matrix for all 1300 PFNS data points reported here. These covariances are based on a thorough assessment of systematic bias and uncertainties associated with the measurement, PFNS extraction technique, combination of data from each detector type, and other aspects of the analysis. The existence of covariances between PFNS data points in different incident neutron energy ranges yielded covariances between average PFNS energy values at each incident energy to be reported here as well, which allowed for firm statements to be made regarding a shape of a purely experimental mean PFNS energy trend for the first time. Although minor PFNS shape differences exist between the results reported here and recent nuclear data evaluations, the ENDF/B-VIII.0 and JEFF-3.3 PFNS evaluations agree reasonably well with the present results from 1-to 10-MeV incident neutron energy, which spans the well-measured 1.5-MeV incident neutron energy PFNS from Lestone and Shores as well as the onset of second-chance fission. However, while the pre-equilibrium component of the PFNS above 12-MeV incident neutron energy roughly agrees in position and magnitude with ENDF/B-VIII.0 and JEFF-3.3, clear differences relating to the relative magnitude of third-chance fission PFNS features are present in the PFNS shape and in the mean PFNS energy trends.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Characterizing Temporal Heterogeneity by Quantifying Nanoscale Fluctuations in Amorphous Fe‐Ge Magnetic Films

Abstract Equilibrium phase transitions are influenced by fluctuations and often discussed within the framework of the Gibbs free energy, wherein the exchange of energy between system and thermal bath is stationary and all regions of the sample exhibit the same phase. Presence of spatial heterogeneity in the magnetic structures such as pinning centers, domain walls, topological defects, etc. may cause temporal heterogeneity that modifies the nature of the magnetic phase transition. This study reports that interplay of nanoscale thermodynamics with spatio‐temporal heterogeneity gives rise to complex phase transition pathways in amorphous Fe x Ge 1‐x thin films with temperature and Fe‐concentration ( x ). Coherent resonant soft X‐ray scattering experiments that have simultaneous spatial, temporal, and spectral sensitivity show that the origin of helical to paramagnetic phase transition in amorphous Fe‐Ge thin films lies in the appearance of enhanced‐fluctuation spots deep inside the ordered state. The fluctuations are heterogeneous, starting over a small fraction of the domains that increases and becomes isotropic over the entire film as the temperature increases or the Fe‐concentration decreases. The fluctuating‐fraction, when normalized to magnetization for different Fe‐concentrations, follows a single power law behavior, suggesting that the nature of the transition can be described in terms of the underlying spatio‐temporal fluctuations.

Singh, Arnab↗

Dark freeze-out cogenesis

We propose a new mechanism where a multi-component dark sector generates the observed dark matter abundance and baryon asymmetry and thus addresses the coincidence between the two. The thermal freeze-out of dark matter annihilating into meta-stable dark partners sets the dark matter relic abundance while providing the out-of-equilibrium condition for baryogenesis. The meta-stable state triggers baryon asymmetry production by its decay well after the freeze-out and potentially induces a period of early matter domination before its decay. The dark matter and baryon abundances are related through number conservation within the dark sector (cogenesis). The “coincidence” is a natural outcome with GeV- to TeV-scale symmetric dark matter and the dark sector’s interactions with the Standard Model quarks. We present a UV-complete model and explore its phenomenological predictions, including dark matter direct detection signals, LHC signatures of new massive particles with color charges and long-lived particles with displaced vertices, dark matter-induced nucleon conversions, (exotic) dark matter indirect detection signals, and effects on the cosmological matter power spectrum. As a side result, we provide a novel analytical treatment for dark sector freeze-out, which may prove useful in the study of related scenarios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detection of diffusion limited oxidation in thermally aged ethylene propylene rubber cable insulation material

The lifetime prediction of electrical cable insulation in nuclear power plants (NPPs) is primarily based on accelerated aging. Polymeric insulation exposed to accelerated aging at high temperatures may demonstrate diffusion limited oxidation (DLO) where oxygen within the material is consumed faster than can be supplied by diffusion from the surrounding air. This situation leads to equilibrium oxidation at the exposed surfaces and limited oxidation away from the surfaces. Therefore, DLO results in a heterogeneously aged material which can lead to underestimation of lifetime for in-service cable insulation. In this study, the presence of inhomogeneous aging caused by DLO was evaluated on cross-sectioned cable insulation specimens made of ethylene propylene rubber (EPR). The studies were performed on thermally aged specimens (165 °C, 1 day and 165 °C, 16 days) using local measurement techniques including nanoindentation, micro-Fourier transform infrared spectroscopy (FTIR), and scanning electron microscopy (SEM) - energy dispersive spectroscopy (EDS). All the results indicate some degree of heterogenous aging, with 165 °C, 16 days, aging demonstrating highest degree of oxidative degradation on the inside edges presumably caused due to copper catalyzed oxidation.

Pallaka, Madhusudhan Reddy↗