Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical mechanics”

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

Experimental validation of simplified radionuclide transport bubble scrubbing code in sodium coolant pool

Solid radionuclides released following fuel pin failure may become entrained in gaseous fission products and be rapidly transported through the coolant pool to the cover gas region. Sensitivity studies on radionuclide transport have identified this potential pathway as one of both high concern and high uncertainty. The simplified radionuclide transport scrubbing code utilizes classical aerosol scrubbing mechanisms to model this phenomenon and predict aerosol masses reaching the cover gas region. This paper serves to validate this code in a sodium environment while conducting a parametric study to analyze the effects of aerosol size, bubble size, pool temperature, pool depth, aerosol density, and aerosol concentration. Through a series of experimental tests, it was determined that aerosol sizes ranging from 0.001 to 1.0 µm are of primary concern as aerosols in this range have a ratio of aerosol mass entering the sodium pool from the fuel pin to aerosol mass exiting the sodium pool to the cover gas region of less than 10. Further, the experimental results were found to match the trends found in the scrubbing model closely, but significantly more scrubbing was seen experimentally. Decreasing bubble size and increasing pool depth and aerosol density were all found to increase scrubbing both experimentally and theoretically. Pool temperature was found to have a negligible effect on scrubbing amounts; however, this was largely due to a subsequent increase in bubble size corresponding to increasing temperatures which offset the increase in scrubbing due to increased temperature. Varying aerosol concentration was found to have no effect on scrubbing ratios. A final series of tests was conducted for a more prototypic fuel pin failure with a heterogenous bubble swarm. From these tests, it was found that the experimental scrubbing quantities were larger than for the single bubble case. Overall, it was found that the simplified bubble transport scrubbing code accurately models the trends of the bubble scrubbing but provides a conservative estimate of scrubbing quantities. Model limitations fail to model the complex phenomena present for fission product scrubbing via bubble transport, but match trends seen experimentally.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Development and assessment of models for turbulent Rayleigh-Taylor mixing using the macroscopic forcing method

Reynolds-Averaged Navier Stokes (RANS) simulations are a popular method for designing ICF experiments, and accurate mixing models are crucial for these simulations to give good predictions. To this end, the present work seeks to demonstrate the Macroscopic Forcing Method (MFM) as a tool for both improving existing RANS models as well as assessing RANS model forms. First, MFM analysis from Lavacot et al. (Phys. Rev. Fluids, 2025) is used to develop the k–L–F model, an extension of the k–L model of Dimonte and Tipton (Phys. Fluids, 2006) that incorporates nonlocality through addition of a turbulent species flux transport equation. MFM is then applied to the k–L–F model along with the k–L and BHR–4 models to assess their forms and compare the model-implied eddy diffusivity moments to those measured from high-fidelity simulations. Furthermore, the analysis reveals that models incorporating nonlocality (k–L–F and BHR–4) match the high-fidelity simulation data better than purely local models (k–L), both in terms of mean fields and eddy diffusivity moments. However, all of the considered RANS models struggle to match temporal moments at high Atwood numbers, highlighting the importance of temporal nonlocality in these regimes and the need for additional improvement even among models incorporating nonlocality.

general physics↗
DOE OSTI2021

Density Functional Tight-Binding Simulations Reveal the Presence of Surface Defects on the Quartz (101)–Water Interface

Understanding the structure and reactivity of quartz–water interfaces is critical for numerous applications in the geological, environmental, and biological sciences. However, disagreements on the atomic-level structure of the interfaces between experiments and simulations are hampering our ability to predict the surface reactivity. Here, we used density functional tight-binding (DFTB)-based molecular dynamics to simulate a series of quartz (101) surfaces having different types and densities of surface defects in water and compared them with the structures determined by X-ray reflectivity measurements. The DFTB simulations are able to reproduce previous classical and quantum mechanical predictions of the pristine quartz (101)–water interface that disagree with experimental observations. To remedy this situation, a set of defective quartz surfaces having various surface silicon (Si) vacancies were built as indicated by recent experimental studies. We found that the rotation of surface [SiO 4 ] tetrahedra near Si vacancies can lead to outward displacements of Si atoms similar to those observed in the experiments. The presence of additional surface Si vacancies caused inward relaxations of terminal oxygens through the formation of hydrogen bonds. Finally, the overall results indicate that the quartz (101)–water interface may include a mixture of geminal (≡Si–(OH) 2 )- and vicinal (≡Si–OH)-type silanol groups together with the presence of surface Si vacancies.

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

Design Principles for the Synthesis of Self-Pillared ZSM-5 Zeolite Nanosheets

The design of next-generation materials for emerging energy and environmental applications heavily relies on empirical approaches to direct nonclassical nucleation and crystal growth pathways, polymorphism, intercrystalline transformations, and seed-assisted growth processes. A long-standing obstacle to nanoporous materials design is the complexity of their crystallization, which hinders the development of predictive models and/or physical descriptors that can guide their synthesis. In this study, we use a combination of state-of-the-art synthesis, characterization, and computational design to prepare hierarchical MFI-type zeolites, which we couple with benchmark catalytic testing to assess structure-property-performance relationships. These hierarchical materials are intergrowths of two commercially relevant zeolite frameworks, MFI and MEL, prepared as self-pillared pentasil (SPP) zeolites through seed-assisted, organic-free syntheses for which little theoretical guidance existed. Here, by comparing a large library of zeolite seeds with different pore sizes, dimensions, and structural composite building units, we determined the relative impact of seed and silica source selection, among other synthesis variables. Combined experimental and computational studies are used to test several hypotheses in literature to rationalize the choice of seed structure and establish a more robust selection criteria for seed-assisted synthesis of zeolites. Specifically, we show that a data-driven approach to develop structural descriptors correlates to new, facile routes to rationally design SPP zeolites, addressing knowledge gaps in the fundamental understanding of (non)classical crystal growth mechanisms that are characteristic of nanoporous aluminosilicates.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Substrate-controlled dynamics of spin qubits in low dimensional van der Waals materials

We report a theoretical study of the coherence dynamics of spin qubits in two-dimensional materials (2DMs) and van der Waals heterostructures, as a function of the host thickness and the composition of the surrounding environment. We focus on MoS2 and WS2, two promising systems for quantum technology applications, and we consider the decoherence arising from the interaction of the spin qubit with nuclear spins. We show that the Hahn-echo coherence time is determined by a complex interplay between the source of decoherence in the qubit host and in the environment, which in turn determines whether the noise evolution is in a classical or quantum mechanical regime. We suggest that the composition and thickness of van der Waals heterostructures encapsulating a qubit host can be engineered to maximize coherence times. Finally, we discuss how quantum sensors may be able to probe the dynamics of the nuclear bath in 2DMs.

Onizhuk, Mykyta (ORCID:0000000304344575)↗
DOE OSTI2025

Real classical shadows

Efficiently learning expectation values of a quantum state using classical shadow tomography has become a fundamental task in quantum information theory. In a classical shadows protocol, one measures a state in a chosen basis $\mathcal{W}$ after it has evolved under a unitary transformation randomly sampled from a chosen distribution $\mathcal{U}$. In this work we study the case where $\mathcal{U}$ corresponds to either local or global orthogonal Clifford gates, and $\mathcal{W}$ consists of real-valued vectors. Our results show that for various situations of interest, this ‘real’ classical shadow protocol improves the sample complexity over the standard scheme based on general Clifford unitaries. For example, when one is interested in estimating the expectation values of arbitrary real-valued observables, global orthogonal Cliffords typically decrease the required number of samples by a factor of two. More dramatically, for k-local observables composed only of real-valued Pauli operators, sampling local orthogonal Cliffords leads to a reduction by an exponential-in-k factor in the sample complexity over local unitary Cliffords. Finally, we show that by measuring in a basis containing complex-valued vectors, orthogonal shadows can, in the limit of large system size, exactly reproduce the original unitary shadows protocol.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Gravitational Contact Interactions and the Physical Equivalence of Weyl Transformations in Effective Field Theory

Theories of scalars and gravity, with non-minimal interactions, $\sim (M_P^2 +F(\phi) )R +L(\phi)$, have graviton exchange induced contact terms. These terms arise in single particle reducible diagrams with vertices $\propto q^2$ that cancel the Feynman propagator denominator $1/q^2$ and are familiar in various other physical contexts. In gravity these lead to additional terms in the action such as $\sim F(\phi) T_\mu^\mu(\phi)/M_P^2$ and $F(\phi)\partial^2 F(\phi)/M_P^2$. The contact terms are equivalent to induced operators obtained by a Weyl transformation that removes the non-minimal interactions, leaving a minimal Einstein-Hilbert gravitational action. This demonstrates explicitly the equivalence of different representations of the action under Weyl transformations, both classically and quantum mechanically. To avoid such "hidden contact terms" one is compelled to go to the minimal Einstein-Hilbert representation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Creating kinks with quantum mediation

We consider the creation of kink-antikink pairs of a scalar field ϕ by the scattering of classical wave packets of a second scalar field ψ when there are no direct interactions between ϕ and ψ . The creation becomes possible only due to a quantum field that interacts with both ϕ and ψ . We scan parameter space and find it favorable for kink production when the initial wave packets have large total energy and wide spatial extent but scatter at low velocities. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effective Field Theory for Extreme Mass Ratio Binaries

We derive an effective field theory describing a pair of gravitationally interacting point particles in an expansion in their mass ratio, also known as the self-force (SF) expansion. The 0SF dynamics are trivially obtained to all orders in Newton’s constant by the geodesic motion of the light body in a Schwarzschild background encoding the gravitational field of the heavy body. The corrections at 1SF and higher are generated by perturbations about this configuration—that is, the geodesic deviation of the light body and the fluctuation graviton—but crucially supplemented by an operator describing the recoil of the heavy body as it interacts with the smaller companion. Using this formalism we compute new results at third post-Minkowskian order for the conservative dynamics of a system of gravitationally interacting massive particles coupled to a set of additional scalar and vector fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Quasiclassical Gluon Fields and Low’s Soft Theorem at Small Momentum-Fraction x

In the high-energy limit, soft gluons can be approximately described by quasiclassical gluon fields. It is well known that the gluon field is a pure gauge field on the transverse plane at eikonal order. We derived the complete next-to-eikonal order solutions of the classical Yang-Mills equations for soft gluons in the dense nuclear regime. Utilizing these solutions, it is shown that Low’s soft theorem at small x can be obtained by considering off-diagonal matrix elements of quasiclassical chromoelectric field between single-gluon states in the dilute regime. We further propose on extending Low’s soft theorem at small x to incorporate the effects of gluon saturation in the dense regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Sign Problem in Tensor-Network Contraction

We investigate how the computational difficulty of contracting tensor networks depends on the sign structure of the tensor entries. Using results from computational complexity, we observe that the approximate contraction of tensor networks with only positive entries has lower computational complexity as compared to tensor networks with general real or complex entries. This raises the question of how this transition in computational complexity manifests itself in the hardness of different tensor-network-contraction schemes. We pursue this question by studying random tensor networks with varying bias toward positive entries. First, we consider contraction via Monte Carlo sampling and find that the transition from hard to easy occurs when the tensor entries become predominantly positive; this can be understood as a tensor-network manifestation of the well-known negative-sign problem in quantum Monte Carlo. Second, we analyze the commonly used contraction based on boundary tensor networks. The performance of this scheme is governed by the number of correlations in contiguous parts of the tensor network (which by analogy can be thought of as entanglement). Remarkably, we find that the transition from hard to easy—i.e., from a volume-law to a boundary-law scaling of entanglement—already occurs for a slight bias of the tensor entries toward a positive mean, scaling inversely with the bond dimension D , and thus the problem becomes easy the earlier the larger D occurs. This is in contrast both to expectations and to the behavior found in Monte Carlo contraction, where the hardness at fixed bias increases with the bond dimension. To provide insight into this early breakdown of computational hardness and the accompanying entanglement transition, we construct an effective classical statistical-mechanical model that predicts a transition at a bias of the tensor entries of 1 / D , confirming our observations. We conclude by investigating the computational difficulty of computing expectation values of tensor-network wave functions (projected entangled-pair states, PEPSs) and find that in this setting, the complexity of entanglement-based contraction always remains low. We explain this by providing a local transformation that maps PEPS expectation values to a positive-valued tensor network. This not only provides insight into the origin of the observed boundary-law entanglement scaling but also suggests new approaches toward PEPS contraction based on positive decompositions. Published by the American Physical Society 2025

Chen, Jielun (ORCID:0000000178411545)↗
DOE OSTI2026

Energy-momentum tensor in a classical model of the electron

We show that the leading nonanalytic terms in the small-𝑡 expansion of the energy-momentum tensor form factors of an electrically charged particle in QED can be correctly derived in a classical model of the electron by Białynicki-Birula. Based on the lucidity of the employed exactly solvable model, we comment also on the recently proposed concept of a regularized proton 𝐷-term.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Demonstration of a third-order hierarchy of topological states in a three-dimensional acoustic metamaterial

Classical wave systems have constituted an excellent platform for emulating complex quantum phenomena, such as demonstrating topological phenomena in photonics and acoustics. Recently, a new class of topological states localized in more than one dimension of a D -dimensional system, referred to as higher-order topological (HOT) states, has been reported, offering an even more versatile platform to confine and control classical radiation and mechanical motion. Here, we design and experimentally study a 3D topological acoustic metamaterial supporting third-order (0D) topological corner states along with second-order (1D) edge states and first-order (2D) surface states within the same topological bandgap, thus establishing a full hierarchy of nontrivial bulk polarization–induced states in three dimensions. The assembled 3D topological metamaterial represents the acoustic analog of a pyrochlore lattice made of interconnected molecules, and is shown to exhibit topological bulk polarization, leading to the emergence of boundary states.

Science & Technology - Other Topics↗
DOE OSTI2024

Determining Exact RANS Operators with the Macroscopic Forcing Method (Final Report)

This report contains a compilation of key results and findings from the ACT project “Determining Exact RANS Operators with the Macroscopic Forcing Method.” The Macroscopic Forcing Method (MFM), a numerical tool for determining closure operators, is used to measure eddy diffusivity moments in Rayleigh-Taylor (RT) instability. It is first applied to low-Atwood 2D RT instability; that work is then extended to 3D RT at different finite Atwood numbers. It is found that nonlocality is important for modeling the mean scalar transport closure operator in RT mixing. Additionally, work is done to improve the statistical convergence of MFM for chaotic problems like RT mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

58 GEOSCIENCES↗