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 487 records · Page 27

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The Hubble PanCET Program: Transit and Eclipse Spectroscopy of the Strongly Irradiated Giant Exoplanet WASP-76b

Ultra-hot Jupiters with equilibrium temperatures greater than 2000 K are uniquely interesting targets as they provide us crucial insights into how atmospheres behave under extreme conditions. This class of giant planets receives intense radiation from their host star and usually has strongly irradiated and highly inflated atmospheres. At such a high temperature, cloud formation is expected to be suppressed and thermal dissociation of water vapor could occur. We observed the ultra-hot Jupiter WASP-76b with seven transits and five eclipses using the Hubble Space Telescope and the Spitzer Space Telescope (Spitzer) for a comprehensive study of its atmospheric chemical and physical processes. We detected TiO and H{sub 2}O absorption in the optical and near-infrared transit spectrum. Additional absorption by a number of neutral and ionized heavy metals like Fe, Ni, Ti, and SiO help explain the short-wavelength transit spectrum. The secondary eclipse spectrum shows muted water feature but a strong CO emission feature in Spitzer’s 4.5 μm band indicating an inverted temperature pressure profile. We analyzed both the transit and eclipse spectra with a combination of self-consistent PHOENIX models and atmospheric retrieval. Both spectra were well fitted by the self-consistent PHOENIX forward atmosphere model in chemical and radiative equilibrium at solar metallicity, adding to the growing evidence that both TiO/VO and NUV heavy metals opacity are prominent NUV-optical opacity sources in the stratospheres of ultra-hot Jupiters.

47 OTHER INSTRUMENTATION↗

Fluctuation–dissipation relation for a quantum Brownian oscillator in a parametrically squeezed thermal field

In this paper we study the nonequilibrium evolution of a quantum Brownian oscillator, modeling the internal degree of freedom of a harmonic atom or an Unruh–DeWitt detector, coupled to a nonequilibrium and nonstationary quantum field bath and inquire whether a fluctuation–dissipation relation (FDR) can exist after/if it approaches equilibration. This is a nontrivial issue because a squeezed field bath cannot reach equilibration and yet, as this work shows, the system oscillator indeed can, which is a necessary condition for FDRs. We discuss three different settings: (A) The bath field essentially remains in a squeezed thermal state throughout, whose squeeze parameter is a mode- and time-independent constant. This situation is often encountered in quantum optics and quantum thermodynamics. (B) The bath field is initially in a thermal state, but is subjected to a parametric process leading to mode- and time-dependent squeezing. This scenario is encountered in cosmology and dynamical Casimir effects. The squeezing in the bath in both types of processes will affect the oscillator’s nonequilibrium evolution. We show that at late times it approaches equilibration and this stationarity condition warrants the existence of a FDR. The trait of squeezing is marked by the oscillator’s effective equilibrium temperature, and the proportionality factor in the FDR is only related to the stationary component of the noise kernel of the bath field. Setting (C) is more subtle: A finite system–bath coupling strength can set the oscillator in a squeezed state even though the bath field is stationary and does not engage in any parametric process. The squeezing of the system in this case is in general time-dependent but becomes constant when the internal dynamics is fully relaxed. We begin with comments on the broad range of physical processes involving squeezed thermal baths and end with some remarks on the significance of FDRs in capturing the essence of quantum backreaction in nonequilibrium and stochastic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development of high-performance long-pulse discharge in KSTAR

High-performance long-pulse plasma operation is essential for producing economically viable fusion energy in tokamak devices. To achieve such discharges in KSTAR, firstly, the rapid increase in the temperature of plasma-facing components was mitigated. The temperature increase of the poloidal limiter, especially, was associated with beam-driven fast ion orbit loss and the discrepancy of the equilibrium reconstructed with heated magnetic probes of signal drift. The fast ions lost to the poloidal limiter were reduced by optimizing the plasma shape and the composition of neutral beam injection (NBI). This nonlinear signal drift was successfully reduced by a new thermal shielding protector on the magnetic probes. Secondly, a lower loop voltage approach was implemented to reduce a poloidal flux consumption rate. A plasma current of 400 kA and a line-averaged electron density of ~2.0 × 10 19 m –3 were chosen by considering the L–H power threshold, fast ion orbit loss, and beam shine-through power loss for low loop voltage in KSTAR. In addition, the application of electron cyclotron heating also helped maintain the plasma with low loop voltage (~25 mV) by enhancing the NBI-driven current and achieving a high poloidal beta (β P ) state. KSTAR has achieved a long pulse (~90 s) operation with the high performance of β P $\leqslant$ 2.7, thermal energy confinement enhancement factor (H 98y2 ) ~ 1.1, and fraction of non-inductive current (f NI ) $\leqslant$ 0.96. Still, gradual degradation of the plasma performance has been observed over time in the discharges. In one of the long-pulse discharges, β P reduced by ~18% over the time of ~8τ R (current relaxation time, τ R ~5 s) and ~1067τ E,th (thermal energy confinement time, τ E,th ~45 ms). The degradation may be closely associated with weak, yet growing, and persistent toroidal Alfvén eigenmodes and their effect on fast ion confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vibrational spectroscopy analysis of silica and silicate glass networks

Vibrational spectroscopy has been widely used to investigate various structural aspects of the glass network, and there are a plethora of papers reporting subtle but consistent changes in infrared and Raman spectral features of glass upon alterations of glass compositions, thermal histories, mechanical stresses, or surface treatments. However, interpretations of such spectral features are still obscured due to the lack of well-established physical principles accurately describing vibrational modes of the non-crystalline glass network. Due to the non-equilibrium nature of the glass network, three-dimensionally connected without any long-range orders, vibrational spectral features of glass cannot be interpreted using the analogy with those of isolated molecular moieties or crystalline counterparts. Here, this feature article explains why such comparisons are outdated and describes the recent advances made from theoretical calculations of vibrational spectral features of amorphous networks or comparisons of computational results with experimental data. For the interpretation of vibrational spectral features of silica and silicate glasses, the following empirical relationships are suggested: (i) the intensity-weighted peak position of the Si-O-Si stretch mode negatively correlates with the weighted average of the Si-O bond length distribution, and (ii) the broad band of the Si-O-Si bending mode negatively correlates with the Si-O-Si bond angle distribution. Selected examples of vibrational spectroscopic imaging of surface defects are discussed to deliberate the implication of these findings in the structure-property relationship of silica and silicate glass materials. Unanswered questions and continuing research challenges are identified.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Improved treatment of multi-material cells in thermal radiation transport codes

High-energy-density physics simulations with non-conformal meshes of materials require a multi-material (MM) closure that affects the thermal radiation transport (TRT). We propose a set of novel closures that work for an arbitrary number of materials, both grey and multigroup energy discretizations, and any angular discretization (such as Sn, IMC, or diffusion). For each spatial cell, our closures let each species (ion and electron) of each material have its own temperature, density, and internal energy, but use a single radiation distribution that interacts with all materials within the cell. Our closures maintain energy conservation, do not incur increased computational cost in the TRT solve itself, are compatible with single-material TRT solvers, and do not make any temperature-equilibrium assumptions. Here we test our closures on a wide range of increasingly realistic problems and find them to be robust.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Revealing nanoscale dynamics during an epoxy curing reaction with x-ray photon correlation spectroscopy

In this work the evolution of nanoscale properties is measured during the thermally triggered curing of an industrial epoxy adhesive. We use x-ray photon correlation spectroscopy (XPCS) to track the progression of the curing reaction through the local dynamics of filler particles that reflect the formation of a thermoset network. Out-of-equilibrium dynamics are resolved through identification and analysis of the intensity– intensity autocorrelation functions obtained from XPCS. The characteristic time scale and local velocity of the filler is calculated as functions of time and temperature. We find that the dynamics speed up when approaching the curing temperature (T cure ), and decay rapidly once T cure is reached. We compare the results from XPCS to conventional macroscale characterization by differential scanning calorimetry (DSC). The demonstration and implementation of nanoscale characterization of curing reactions by XPCS proves useful for future development and optimization of epoxy thermoset materials and other industrial adhesive systems.

36 MATERIALS SCIENCE↗

Direct imaging of electrical switching of antiferromagnetic Néel order in α-Fe 2 O 3 epitaxial films

In this work, we report the direct observation of switching of the Néel vector of antiferromagnetic (AFM) domains in response to electrical pulses in micron-scale Pt/α-Fe 2 O 3 Hall bars using photoemission electron microscopy. Current pulses lead to reversible and repeatable switching with the current direction determining the final state, consistent with Hall effect experiments that probe only the spatially averaged response. Current pulses also produce irreversible changes in domain structure, in and even outside the current path. In both cases only a fraction of the domains switch in response to pulses. Furthermore, the analysis of images taken with different x-ray polarizations shows that the AFM Néel order has an out-of-plane component in equilibrium that is important to consider in analyzing the switching data. These results show that—in addition to effects associated with spin-orbit torques from the Pt layer, which can produce reversible switching—changes in AFM order can be induced by purely thermal effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of Equiatomic CrMnFeCoNiCu System and Subsequent Derivation of a Non-Equiatomic MnFeCoNiCu Alloy

Investigation into non-equiatomic high-entropy alloys has grown in recent years due to questions about the role of entropy stabilization in forming single-phase solid solutions. Non-equiatomic alloys have been shown to retain the outstanding mechanical properties exhibited by their equiatomic counterparts and even improve electrical, thermal, and magnetic properties, albeit with relaxed composition bounds. However, much remains to understand the processing–structure–property relationships in all classes of so-called high-entropy alloys (HEAs). Here, we are motivated by the natural phenomena of crystal growth and equilibrium conditions to introduce a method of HEA development where controlled processing conditions determine the most probable and stable composition. This is demonstrated by cooling an equiatomic CrMnFeCoNiCu alloy from the melt steadily over 3 days (cooling rate ~4 °C/h). The result is an alloy containing large Cr-rich precipitates and an almost Cr-free matrix exhibiting compositions within the MnFeCoNiCu system (with trace amounts of Cr). From this juncture, it is argued that the most stable composition is within the CrMnFeCoNiCu system rather than the CrMnFeCoNi system. With further optimization and evaluation, a unique non-equiatomic alloy, Mn 17 Fe 21 Co 24 Ni 24 Cu 14 , is derived. The alloy solidifies and recrystallizes into a single-phase face-centered cubic (FCC) polycrystal. In addition to possible applications where Invar is currently utilized, this alloy can be used in fundamental studies that contrast its behavior with its equiatomic counterpart and shed light on the development of HEAs.

36 MATERIALS SCIENCE↗

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomistic modeling of metal–nonmetal interphase boundary diffusion

Atomistic computer simulations are applied to investigate the atomic structure, thermal stability, and diffusion processes in Al–Si interphase boundaries as a prototype of metal–ceramic interfaces in composite materials. Some of the most stable orientation relationships between the phases found in this work were previously observed in epitaxy experiments. Here, a non-equilibrium interface can transform to a more stable state by a mechanism that we call interface-induced recrystallization. Diffusion of both Al and Si atoms in stable Al–Si interfaces is surprisingly slow compared with diffusion of both elements in Al grain boundaries but can be accelerated in the presence of interface disconnections. A qualitative explanation of the sluggish interphase boundary diffusion is proposed. Atomic mechanisms of interphase boundary diffusion are similar to those in metallic grain boundaries and are dominated by correlated atomic rearrangements in the form of strings and rings of collectively moving atoms.

36 MATERIALS SCIENCE↗

Defect engineering in nickel via electrodeposition following low temperature annealing

Electrodeposited nickel coatings are characterized by high densities of crystalline defects including dislocations, growth twins, and hydrogen bubbles/voids, arising from the non-equilibrium nature of deposition. While these metastable features influence as-deposited properties, their evolution under low temperature annealing remains unexplored. Here, in this study, we demonstrate that low temperature annealing (200 °C) induces significant microstructural rearrangement without grain growth, enabling defect engineering through thermally assisted dislocation motion and twin boundary migration. Notably, we observe faceted twin boundaries and dislocation organization into extended networks, which have been rarely reported under such annealing conditions. These transformations are attributed to hydrogen-assisted defect mobility, facilitated by the release or redistribution of hydrogen trapped during deposition. These structural transformations correlate with a significant enhancement in mechanical properties, including a twofold increase in yield strength and improved ductility. Our findings highlight the role of trapped hydrogen in mediating low-temperature defect mobility and twin boundary evolution, offering a unique pathway for microstructural tuning of metallic coatings through controlled annealing.

36 - MATERIALS SCIENCE↗

Entangled core/shell magnetic structure driven by surface magnetic symmetry-breaking in Cr 2 O 3 nanoparticles

Bulk Cr 2 O 3 is an antiferromagnetic (AFM) oxide that exhibits the magnetoelectric effect at room temperature, with neither spontaneous magnetization nor net electric polarization. These physical properties stem from a subtle competition between exchange and crystal field interactions. In this article, we exploit the symmetry breaking at the surface of Cr 2 O 3 nanoparticles for unbalancing this delicate physical equilibrium. The emerging weak ferromagnetic signal we observe persists up to near room temperature (≈ 270 K) at which the antiferromagnetic order disappears. In addition, an exchange-bias effect, that rapidly decreases on heating from low temperature up to 30 K, is resistant to thermal disorder above 200 K. In conclusion, our findings point to the possible formation of an entangled core/shell magnetic structure, where pinned uncompensated spins at the shell are randomly distributed in a low-temperature spin-glass ordering, with low net magnetic moment and an ordering temperature governed by the AFM Néel temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of nonequilibrium magnons in gapped Heisenberg antiferromagnets

Nonequilibrium dynamics in spin systems is a topic currently under intense investigation as it provides fundamental insights into thermalization, universality, and exotic transport phenomena. While most of the studies have been focused on ideal closed quantum many-body systems such as ultracold atomic quantum gases and one-dimensional spin chains, driven-dissipative Bose gases in steady states away from equilibrium in classical systems also lead to intriguing nonequilibrium physics. Here, in this work, we theoretically investigate out-of-equilibrium dynamics of magnons in a gapped Heisenberg quantum antiferromagnet based on Boltzmann transport theory. We show that, by treating scattering terms beyond the relaxation-time approximation in the Boltzmann transport equation, energy and particle number conservation mandate that nonequilibrium magnons cannot relax to equilibrium, but decay to other nonequilibrium stationary states. The only decay channel for these stationary states back to equilibrium is through the nonconserving interactions (i.e., changing particle number and/or energy within the magnon system) such as boundary or magnon-phonon scattering. At low temperatures, these nonconserving interactions are much slower processes than intrinsic magnon-magnon interaction in a gapped spin system. Using magnon-phonon interaction as a quintessential type of nonconserving interaction, we then propose that nonequilibrium steady states of magnons can be maintained and tailored using periodic driving at frequencies faster than relaxation due to phonon interactions. These findings reveal a class of classical material systems that are suitable platforms to study nonequilibrium statistical physics and macroscopic phenomena such as classical Bose-Einstein condensation of quasiparticles and magnon supercurrents that are relevant for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Correlations Between Fission Fragment and Neutron Anisotropies in Neutron-Induced Fission

Several sources of angular anisotropy for fission fragments and prompt neutrons have been studied in neutron-induced fission reactions. These include kinematic recoils of the target from the incident neutron beam and the fragments from the emission of the prompt neutrons, preferential directions of the emission of the fission fragments with respect to the beam axis due to the population of particular transition states at the fission barrier, and forward-peaked angular distributions of pre-equilibrium neutrons which are emitted before the formation of a compound nucleus. In addition, there are several potential sources of angular anisotropies that are more difficult to disentangle: the angular distributions of prompt neutrons from fully accelerated fragments or from scission neutrons, and the emission of neutrons from fission fragments that are not fully accelerated. Here, we study the effects of the first group of anisotropy sources, particularly exploring the correlations between the fission-fragment anisotropy and the resulting neutron anisotropy. While kinematic effects were already accounted for in our Hauser-Feshbach Monte Carlo code cgmf, anisotropic angular distributions for the fission fragments and pre-equilibrium neutrons resulting from neutron-induced fission on 233,234,235,238 U, 239,241 Pu, and 237 Np have been introduced for the first time. The effects of these sources of anisotropy are examined over a range of incident neutron energies, from thermal to 20 MeV, and compared with experimental data from the Chi-Nu liquid scintillator array. The anisotropy of the fission fragments is reflected in the anisotropy of the prompt neutrons, especially as the outgoing energy of the prompt neutrons increases, allowing for an extraction of the fission-fragment anisotropy to be made from a measurement of the neutrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗