Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical mechanics models”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A multisite decomposition of the tensor network path integrals

Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.

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

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗
DOE OSTI2024

Isotope effects in supercooled H 2 O and D 2 O and a corresponding-states-like rescaling of the temperature and pressure

Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H 2 O and D 2 O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations—with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a “corresponding-states-like” rescaling for water properties in three classical water models that all exhibit a liquid–liquid transition and critical point, the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H 2 O and D 2 O is investigated here. Utilizing previously published data and equations-of-state for H 2 O and D 2 O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid–liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water’s anomalies.

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

On-the-fly autonomous control of neutron diffraction via physics-informed Bayesian active learning

We demonstrate the first live, autonomous control over neutron diffraction experiments by developing and deploying ANDiE: the autonomous neutron diffraction explorer. Neutron scattering is a unique and versatile characterization technique for probing the magnetic structure and behavior of materials. However, instruments at neutron scattering facilities in the world is limited, and instruments at such facilities are perennially oversubscribed. We demonstrate a significant reduction in experimental time required for neutron diffraction experiments by implementation of autonomous navigation of measurement parameter space through machine learning. Prior scientific knowledge and Bayesian active learning are used to dynamically steer the sequence of measurements. We show that ANDiE can experimentally determine the magnetic ordering transition of both MnO and Fe 1.09 Te all while providing a fivefold enhancement in measurement efficiency. Furthermore, in a hypothesis testing post-processing step, ANDiE can determine transition behavior from a set of possible physical models. ANDiE's active learning approach is broadly applicable to a variety of neutron-based experiments and can open the door for neutron scattering as a tool of accelerated materials discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Controlling morphology in hybrid isotropic/patchy particle assemblies

Brownian dynamics is used to study self-assembly in a hybrid system of isotropic particles (IPs), combined with anisotropic building blocks that represent special “designer particles.” Those are modeled as spherical patchy particles (PPs) with binding only allowed between their patches and IPs. In this study, two types of PPs are considered: Octahedral PPs (Oh-PPs) and Square PPs (Sq-PPs), with octahedral and square arrangements of patches, respectively. The self-assembly is additionally facilitated by the simulated annealing procedure. Here, the resultant structures are characterized by a combination of local correlations in cubatic ordering and a symmetry-specific variation of bond orientation order parameters (SymBOPs). By varying the PP/IP size ratio, we detected a sharp crossover between two distinct morphologies in both types of systems. High symmetry phases, NaCl crystal for Oh-PP and square lattice for Sq-PP, are observed for larger size ratios. For the smaller ones, the dominant morphologies are significantly different, e.g., Oh-PPs form a compact amorphous structure with predominantly face-to-face orientation of neighboring PPs. Unusually, for a morphology without a long-range order, it is still possible to identify well organized coherent clusters of this structure, thanks to the adoption of our SymBOP-based characterization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Constrained non-negative matrix factorization enabling real-time insights of in situ and high-throughput experiments

Non-negative matrix factorization (NMF) is an appealing class of methods for performing unsupervised learning on streaming spectral data, particularly in time-sensitive applications such as in situ characterization of materials. These methods seek to decompose a dataset into a small number of components and weights that can compactly represent the underlying signal while effectively reconstructing the observations with minimal error. However, canonical NMF methods have no underlying requirement that the reconstruction uses components or weights that are representative of the true physical processes. In this work, we demonstrate how constraining a subset of the NMF weights or components as rigid priors, provided as known or assumed values, can provide significant improvement in revealing true underlying phenomena. We present a PyTorch-based method for efficiently applying constrained NMF and demonstrate its application to several synthetic examples. Our implementation allows an expert researcher-in-the-loop to provide and dynamically adjust the constraints during a live experiment involving streaming spectral data. Such interactive priors allow researchers to specify known or identified independent components, as well as functional expectations about the mixing or transitions between the components. We further demonstrate the application of this method to measured synchrotron x-ray total scattering data from in situ beamline experiments. In such a context, constrained NMF can result in a more interpretive and scientifically relevant decomposition than canonical NMF or other decomposition techniques. As a result, the details of the method are provided, along with general guidance for employing constrained NMF in the extraction of critical information and insights during time-sensitive experimental applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gauge-free duality in pure square spin ice: Topological currents and monopoles

We consider pure square spin ice, that is, square ice, where only nearest neighbors are coupled. The gauge-free duality between the perpendicular and collinear structure leads to a natural description in terms of topological currents and charges as the relevant degrees of freedom. That, in turn, can be expressed via a continuous field theory where the discrete spins are subsumed into entropic interactions among charges and currents. This approach produces structure factors, correlations, and susceptibilities for spins, monopoles, and currents. It also generalizes the height formalism of the disordered ground state to non-zero temperature. The framework can be applied to the zoology of recent experimental results, especially realizations on quantum annealers, and can be expanded to include longer range interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗
NASA NTRS1994

Damage Mechanisms and Controlled Crack Propagation in a Hot Pressed Silicon Nitride Ceramic

The subcritical growth of cracks from pre-existing flaws in ceramics can severely affect the structural reliability of a material. The ability to directly observe subcritical crack growth and rigorously analyze its influence on fracture behavior is important for an accurate assessment of material performance. A Mode I fracture specimen and loading method has been developed which permits the observation of stable, subcritical crack extension in monolithic and toughened ceramics. The test specimen and procedure has demonstrated its ability to generate and stably propagate sharp, through-thickness cracks in brittle high modulus materials. Crack growth for an aluminum oxide ceramic was observed to be continuously stable throughout testing. Conversely, the fracture behavior of a silicon nitride ceramic exhibited crack growth as a series of subcritical extensions which are interrupted by dynamic propagation. Dynamic initiation and arrest fracture resistance measurements for the silicon nitride averaged 67 and 48 J/sq m, respectively. The dynamic initiation event was observed to be sudden and explosive. Increments of subcritical crack growth contributed to a 40 percent increase in fracture resistance before dynamic initiation. Subcritical crack growth visibly marked the fracture surface with an increase in surface roughness. Increments of subcritical crack growth loosen ceramic material near the fracture surface and the fracture debris is easily removed by a replication technique. Fracture debris is viewed as evidence that both crack bridging and subsurface microcracking may be some of the mechanisms contributing to the increase in fracture resistance. A Statistical Fracture Mechanics model specifically developed to address subcritical crack growth and fracture reliability is used together with a damaged zone of material at the crack tip to model experimental results. A Monte Carlo simulation of the actual experiments was used to establish a set of modeling input parameters. It was demonstrated that a single critical parameter does not characterize the conditions required for dynamic initiation. Experimental measurements for critical crack lengths, and the energy release rates exhibit significant scatter. The resulting output of the model produces good agreement with both the average values and scatter of experimental measurements.

Calomino, Anthony Martin↗
NASA NTRS1996

There aren't Non-Standard Solutions for the Braid Group Representations of the QYBE Associated with 10-D Representations of SU(4)

It is well known that the quantum Yang-Baxter equations (QYBE) play an important role in various theoretical and mathematical physics, such as completely integrable system in (1 + 1)-dimensions, exactly solvable models in statistical mechanics, the quantum inverse scattering method and the conformal field theories in 2-dimensions. Recently, much remarkable progress has been made in constructing the solutions of the QYBE associated with the representations of lie algebras. It is shown that for some cases except the standard solutions, there also exist new solutions, but the others have not non-standard solutions. In this paper by employing the weight conservation and the diagrammatic techniques we show that the solution associated with the 10-D representations of SU (4) are standard alone.

Yijun, Huang↗
NASA NTRS2000

Scale-Invariant Forms of Conservation Equations in Reactive Fields and a Modified Hydro-Thermo-Diffusive Theory of Laminar Flames

A scale-invariant model of statistical mechanics is applied to present invariant forms of mass, energy, linear, and angular momentum conservation equations in reactive fields. The resulting conservation equations at molecular-dynamic scale are solved by the method of large activation energy asymptotics to describe the hydro-thermo-diffusive structure of laminar premixed flames. The predicted temperature and velocity profiles are in agreement with the observations. Also, with realistic physico-chemical properties and chemical-kinetic parameters for a single-step overall combustion of stoichiometric methane-air premixed flame, the laminar flame propagation velocity of 42.1 cm/s is calculated in agreement with the experimental value.

Sohrab, Siavash H.↗
NASA NTRSDate not supplied

Colloids in Microgravity

Colloidal particles are important in a broad range of applications in technology and materials processing. Due to high sensitivity to thermal fluctuations, a suspension of colloidal particles is considered as an assembly of super-atoms in condensed matter physics. This feature provides an excellent experimental model in statistical mechanics for fundamental studies of microstructural evolution and phase transformation at the individual particle level. However, using colloids to probe three-dimensional phase transitions on Earth remains a challenge as they can be influenced by the gravity force causing the particles to sediment and the solvent to flow. To suppress the undesirable gravity effects, experiments on concentrated colloids were carried out in the International Space Station (ISS) Light Microscopy Module. We will present results of these experiments.

Colloids↗
DOE OSTI2022

Understanding the Solid-State Electrode–Electrolyte Interface of a Model System Using First-Principles Statistical Mechanics and Thin-Film X-ray Characterization

Intermixing of atomic species at the electrode-electrolyte boundaries can impact the properties of the interfaces in solid-state batteries. Herein, this work uses first-principles statistical mechanics along with experimental characterization to understand intermixing at the electrode-electrolyte interface. For the model presented in this work, lithium manganese oxide (LiMn 2 O 4 , LMO) and lithium lanthanum titanate (Li 3x La 2/3-x TiO 3 , LLTO) are employed as the cathode and electrolyte, respectively. The results of the computational work show that Ti-Mn intermixing at the interface is significant at synthesis temperatures. The experimental results in this work find that, at some critical temperatures between 600 and 700 degrees C for material preparation, the interface of LLTO-LMO becomes blurred. Calculations predict that the interface is unstable with regard to Ti-Mn intermixing starting at 0 K, suggesting that the critical temperature found in the experiment is related to kinetics. The work overall suggests that, in designing a solid-state battery, the fundamental reactions such as intermixing need to be considered.

25 ENERGY STORAGE↗
NASA NTRS2010

Coherent Eigenmodes in Homogeneous MHD Turbulence

The statistical mechanics of Fourier models of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence is discussed, along with their relevance for dissipative magnetofluids. Although statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation, i.e., we have coherent structure. We use eigenanalysis of the modal covariance matrices in the probability density function to explain this phenomena in terms of `broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We provide examples from 2-D and 3-D magnetohydrodynamic simulations of homogeneous turbulence, and show new results from long-time simulations of MHD turbulence with and without a mean magnetic field

Shebalin, John V.↗
NASA NTRS2010

Broken Ergodicity in Ideal, Homogeneous, Incompressible Turbulence

We discuss the statistical mechanics of numerical models of ideal homogeneous, incompressible turbulence and their relevance for dissipative fluids and magnetofluids. These numerical models are based on Fourier series and the relevant statistical theory predicts that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. However, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We explain this phenomena in terms of broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We review the theoretical basis of broken ergodicity, apply it to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence, and show new results from simulations using GPU (graphical processing unit) computers.

Morin, Lee↗
DOE OSTI2021

Sparse expansions of multicomponent oxide configuration energy using coherency and redundancy

We report that compressed sensing has become a widely accepted paradigm to construct high dimensional cluster expansion models used for statistical mechanical studies of atomic configuration in complex multicomponent crystalline materials. However, strict sampling requirements necessary to obtain minimal coherence measurements for compressed sensing to guarantee accurate estimation of model parameters are difficult and in some cases impossible to satisfy due to the inability of physical systems to access certain configurations. Nevertheless, the dependence of energy on atomic configuration can still be adequately learned without these strict requirements by using compressed sensing by way of coherent measurements using redundant function sets known as frames. We develop a particular frame constructed from the union of all occupancy-based cluster expansion basis sets. We illustrate how using this highly redundant frame yields sparse expansions of the configuration energy of complex oxide materials that are competitive and often surpass the prediction accuracy and sparsity of models obtained from standard cluster expansions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Detecting Quantum Critical Points of Correlated Systems by Quantum Convolutional Neural Network Using Data from Variational Quantum Eigensolver

Machine learning has been applied to a wide variety of models, from classical statistical mechanics to quantum strongly correlated systems, for classifying phase transitions. The recently proposed quantum convolutional neural network (QCNN) provides a new framework for using quantum circuits instead of classical neural networks as the backbone of classification methods. We present the results from training the QCNN by the wavefunctions of the variational quantum eigensolver for the one-dimensional transverse field Ising model (TFIM). We demonstrate that the QCNN identifies wavefunctions corresponding to the paramagnetic and ferromagnetic phases of the TFIM with reasonable accuracy. The QCNN can be trained to predict the corresponding ‘phase’ of wavefunctions around the putative quantum critical point even though it is trained by wavefunctions far away. The paper provides a basis for exploiting the QCNN to identify the quantum critical point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗