Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice Integrable 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling of standing wave RF cavities for tracking through multi-pass energy recovery linac

Short bunches, high current and multiple linac pass are all characteristics of Energy Recovery Linacs (ERLs), which may result in collective effects. They in turn, may affect the beam, degrading its quality, or even yield to instabilities causing a beam loss. To study and mitigate these effects one needs a numerical simulation code, that can take into account both the collective effects, as well as, particular ERL features, such as a multi-turn design that does not reach a steady state or the multiple passages of the beam through Radio-Frequency (RF) cavities at different energies. CODAL [1], a code developed by SOLEIL in collaboration with IJCLab, enables such studies. It is a 6 dimensional (6D) tracking code applying 'kicks' based on the integration of the local Hamiltonian for each element of the lattice. It is also capable of simulating space charge, wakefields and coherent synchrotron radiation. However, to correctly take into account the ERL dynamics, an upgrade had to be made to include the effect of a standing wave RF cavity in 6D. In this paper, we will concentrate on the implementation and benchmarking (with DESY’s tracking code ASTRA [2]) of both the longitudinal and the transverse models (by J.B. Rosenzweig and L. Serafini [3]), which we use to carry out tracking of fully analytical 6D RF cavity.

Bogacz, A.↗
DOE OSTI2021

Metallization of colloidal crystals

Colloidal crystals formed by size-asymmetric binary particles coassemble into a wide variety of colloidal compounds with lattices akin to ionic crystals. Recently, a transition from a compound phase with a sublattice of small particles to a metal-like phase in which the small particles are delocalized has been predicted computationally and observed experimentally. In this colloidal metallic phase, the small particles roam the crystal maintaining the integrity of the lattice of large particles, as electrons do in metals. A similar transition also occurs in superionic crystals, termed sublattice melting. Here, we use energetic principles and a generalized molecular dynamics model of a binary system of functionalized nanoparticles to analyze the transition to sublattice delocalization in different coassembled crystal phases as a function of temperature (T), number of grafted chains on the small particles, and number ratio between the small and large particles n s : n l . Here we find that n s : n l is the primary determinant of crystal type due to energetic interactions and interstitial site filling, while the number of grafted chains per small particle determines the stability of these crystals. We observe first-order sublattice delocalization transitions as T increases, in which the host lattice transforms from low- to high-symmetry crystal structures, including A20 → bct → bcc, A d → bct → bcc, and bcc → bcc/fcc → fcc transitions and lattices. Analogous sublattice transitions driven primarily by lattice vibrations have been seen in some atomic materials exhibiting an insulator-metal transition also referred to as metallization. We also find minima in the lattice vibrations and diffusion coefficient of small particles as a function of n s : n l , indicating enhanced stability of certain crystal structures for n s : n l values that form compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling Transverse Space Charge effects in IOTA with pyORBIT

The role and mitigation of space charge effects are important aspects of the beam physics research to be performed at the Integrable Optics Test Accelerator (IOTA) at Fermilab. The impact of nonlinear integrability (partial and complete) on space charge driven incoherent and coherent resonances will be a new feature of this accelerator and has the potential to influence the design of future low energy proton synchrotrons. In this report we model transverse space charge effects using the PIC code pyORBIT. First we benchmark the single particle tracking results against MADX with checks on symplecticity, tune footprints, and dynamic aperture in a partially integrable lattice realized with a special octupole string insert. Our space charge calculations begin with an examination of the 4D symplecticity. Short term tracking is done first with the initial bare lattice and then with a matching of the rms values with space charge. Next, we explore slow initialization of charge as a technique to establish steady state and reduce emittance growth and beam loss following injection into a mismatched lattice. We establish values of space charge simulation parameters so as to ensure numerical convergence. Finally, we compare the simulated space charge induced tune shifts and footprints against theory.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Spin crossovers and superdiffusion in the one-dimensional Hubbard model

We use tools from integrability and generalized hydrodynamics to study finite-temperature dynamics in the one-dimensional Hubbard model. First, we examine charge, spin, and energy transport away from half-filling and zero magnetization, focusing on the strong coupling regime where we identify a rich interplay of temperature and energy scales, with crossovers between distinct dynamical regimes. We identify an intermediate-temperature regime analogous to the spin-incoherent Luttinger liquid, where spin degrees of freedom are hot but charge degrees of freedom are at low temperature. We demonstrate that the spin Drude weight exhibits sharp features at the crossover between this regime and the low-temperature Luttinger liquid regime, which are absent in the charge and energy response, and rationalize this behavior in terms of the properties of Bethe ansatz quasiparticles.We then turn to the dynamics along special lines in the phase diagram corresponding to half-filling and/or zero magnetization where on general grounds we anticipate that the transport is subballistic but superdiffusive. We provide analytical and numerical evidence for Kardar-Parisi-Zhang (KPZ) dynamical scaling (with length and time scales related via x ~ t 2/3 ) along both lines and at the SO(4)-symmetric point where they intersect. Furthermore, our results suggest that both spin-coherence crossovers and KPZ scaling may be accessed in near-term experiments with optical lattice Hubbard emulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Planar Systems for Quantum Information

This project aims to develop two‐dimensional (2D) moiré materials as a quantum simulator to implement model Hamiltonians and their phase diagrams. Progress in quantum information science (QIS) requires the development of advanced quantum materials systems. The rich family of layered van der Waals materials and their heterostructures present opportunities to create previously unrealized types of applications for QIS. Specifically, when two layers of van der Waals materials are overlaid with a small twist angle or/and lattice mismatch, a moiré superlattice with a period of about ten nanometers is formed. This provides a periodic trapping potential for electrons. Electrons can tunnel between the traps and repel each other by their mutual Coulomb interactions. The platform of 2D moiré materials provides many attractive features, including tunability of length and energy scales, charge density, and even lattice symmetry. It presents new possibilities for realizing quantum simulation of the many-body physics in a solid-state platform. This integrated team of six investigators seeks to develop relevant theoretical treatments to link ab-initio studies of 2D moiré materials to model Hamiltonians and to evaluate correlated phases predicted by these model Hamiltonians in the relevant regimes. On the experimental side, the team aims to develop methods to realize a homogeneous and highly controlled potential landscape for the electrons and to initiate, protect, and measure their quantum many-body states.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

A Data-driven approach to Core Power distribution reconstruction in a Nuclear Reactor

This report presents the initial development of a data-driven approach for reconstructing the core power distribution in a nuclear reactor (power shape synthesis) using ex-core sensors. Traditional techniques rely on deploying a large number of detectors throughout the reactor core. However, this approach is not feasible for innovative reactor concepts like Advanced Reactors and Microreactors. First, the tight lattice pitch, designed to maximize power density, limits the space available for sensors. Secondly, the harsh operating conditions are not compatible with commercially available detectors. The method proposed in this work integrates high-fidelity modeling with data-driven techniques to accurately reconstruct power distribution across various reactor types, thereby reducing the reliance on in-core sensors. Purdue University Reactor One (PUR-1) was selected as the test case. The CAD model representing the latest configuration of the PUR-1 core was imported into the OpenMC simulation framework, and the model was built. Additionally, the previously developed MCNP6 model was updated. The two models were assessed against the data collected during an experimental campaign conducted in July 2024. Thirty gold foils were placed in three Irradiation Assemblies in PUR-1 core. Using the measured activity of the irradiated foils, the neutron flux at different core locations was estimated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Vibrational response and temperature of shock-compressed Pt: In situ extended x-ray absorption fine structure measurements to 325 GPa

Single pulse (~100 ps duration) extended x-ray absorption fine structure (EXAFS) measurements were obtained in laser shocked Pt to examine atomic vibrational properties and to determine temperatures for shock pressures from 72-325 GPa. Pt x-ray absorption spectra were recorded while a planar shock wave was propagating through the Pt sample. The absorption spectrum corresponding to the shocked Pt was obtained by subtracting the ambient Pt contribution from the measured absorption spectrum. Fits to the shocked state EXAFS data provided Pt lattice parameters consistent with the known Pt Hugoniot and the mean squared relative displacements (MSRDs) for nearest neighbor Pt atoms. Pt temperatures in the shocked state, estimated from the MSRDs using the correlated Debye model and correlated Einstein models, were consistent with both the Hugoniot temperatures calculated by integration along the Hugoniot and with the Hugoniot temperatures from published first-principles calculations. However, the agreement between the Hugoniot temperatures from EXAFS measurements and the calculated Hugoniot temperatures is not as good above 200 GPa, likely due to anharmonic effects. Furthermore, the present results demonstrate that single-pulse synchrotron EXAFS measurements in laser-shocked solids are useful for quantitative temperature determination and for examination of vibrational properties, including anharmonicity, at extreme pressure-temperature conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Random forest prediction of crystal structure from electron diffraction patterns incorporating multiple scattering

Diffraction is the most common method to solve for unknown or partially known crystal structures. However, it remains a challenge to determine the crystal structure of a new material that may have nanoscale size or heterogeneities. Here, in this study, we train an architecture of hierarchical random forest models capable of predicting the crystal system, space group, and lattice parameters from one or more unknown two-dimensional electron diffraction patterns. Our initial model correctly identifies the crystal system of a simulated electron diffraction pattern from a 20-nm-thick specimen of arbitrary orientation 67% of the time. We achieve a topline accuracy of 79% when aggregating predictions from ten patterns of the same material but different zone axes. The space group and lattice predictions range from 70% to 90% accuracy and median errors of 0.01-0.5Å, respectively, for cubic, hexagonal, trigonal, and tetragonal crystal systems while being less reliable on orthorhombic and monoclinic systems. We apply this architecture to a four-dimensional scanning transmission electron microscopy scan of gold nanoparticles, where it accurately predicts the crystal structure and lattice constants. These random forest models can be used to significantly accelerate the analysis of electron diffraction patterns, particularly in the case of unknown crystal structures. Additionally, due to the speed of inference, these models could be integrated into live transmission electron microscopy experiments, allowing real-Time labeling of a specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calculation of neutron flux spectra of the VVER-1000 mock-up shielding benchmark with Monte Carlo code MCS utilizing mesh-based weight window

The measurements of neutron spectra compiled inside the NEA-1517/82 package from the Shielding Integral Benchmark Archive and Database (SINBAD) are chosen as benchmark cases to validate the variance reduction technique based on weight window in Monte Carlo Code MCS. A full 3D model for fixed source mode calculation with hexagonal lattice source definition is developed to simulate total of 6 points of measurements at the vicinity of the reactor and the reactor pressure vessel region. A code/code comparison against MCNP6 code is first conducted as verification element for the mesh-based weight window capability in MCS. Finally, the validation results are presented against measurements. The verification against MCNP6 code gives good agreement in addition of the insight to the importance of user understanding to determine proper reference point and reference lower weight bound for scaling which is not required in MCS code due to its capability of automatic scaling. The comparison of neutron spectra between MCS and measurements shows good agreement within 3 standard deviations for all of six detector positions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Renormalization of excitonic properties by polar phonons

Here we employ quasiparticle path integral molecular dynamics to study how the excitonic properties of model semiconductors are altered by electron–phonon coupling. We describe ways within a path integral representation of the system to evaluate the renormalized mass, binding energy, and radiative recombination rate of excitons in the presence of a fluctuating lattice. To illustrate this approach, we consider Fröhlich-type electron–phonon interactions and employ an imaginary time influence functional to incorporate phonon-induced effects nonperturbatively. The effective mass and binding energies are compared with perturbative and variational approaches, which provide qualitatively consistent trends. We evaluate electron-hole recombination rates as mediated through both trap-assisted and bimolecular processes, developing a consistent statistical mechanical approach valid in the reaction limited regime. These calculations demonstrate how phonons screen electron–hole interactions, generically reducing exciton binding energies and increasing their radiative lifetimes.

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

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2022

Digital Quantum Simulation of the Schwinger Model and Symmetry Protection with Trapped Ions

Tracking the dynamics of physical systems in real time is a prime application of quantum computers. Using a trapped-ion system with up to six qubits, we simulate the real-time dynamics of a lattice gauge theory in 1+1 dimensions, i.e., the lattice Schwinger model, and demonstrate nonperturbative effects such as pair creation for times much longer than previously accessible. We study the gate requirement of two formulations of the model using the Suzuki-Trotter product formula, as well as the trade-off between errors from the ordering of the Hamiltonian terms, the Trotter step size, and experimental imperfections. To mitigate experimental errors, a recent symmetry-protection protocol for suppressing coherent errors and a symmetry-inspired postselection scheme are applied. This work demonstrates the integrated theoretical, algorithmic, and experimental approach that is essential for efficient simulation of lattice gauge theories and other complex physical systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coherent electrical control of a single high-spin nucleus in silicon

Nuclear spins are highly coherent quantum objects. In large ensembles, their control and detection via magnetic resonance is widely exploited, for example, in chemistry, medicine, materials science and mining. Nuclear spins also featured in early proposals for solid-state quantum computers and demonstrations of quantum search and factoring algorithms. Scaling up such concepts requires controlling individual nuclei, which can be detected when coupled to an electron. However, the need to address the nuclei via oscillating magnetic fields complicates their integration in multi-spin nanoscale devices, because the field cannot be localized or screened. Control via electric fields would resolve this problem, but previous methods relied on transducing electric signals into magnetic fields via the electron–nuclear hyperfine interaction, which severely affects nuclear coherence. In this work, we demonstrate the coherent quantum control of a single 123 Sb (spin-7/2) nucleus using localized electric fields produced within a silicon nanoelectronic device. The method exploits an idea proposed in 1961 but not previously realized experimentally with a single nucleus. Our results are quantitatively supported by a microscopic theoretical model that reveals how the purely electrical modulation of the nuclear electric quadrupole interaction results in coherent nuclear spin transitions that are uniquely addressable owing to lattice strain. The spin dephasing time, 0.1 seconds, is orders of magnitude longer than those obtained by methods that require a coupled electron spin to achieve electrical driving. These results show that high-spin quadrupolar nuclei could be deployed as chaotic models, strain sensors and hybrid spin-mechanical quantum systems using all-electrical controls. Integrating electrically controllable nuclei with quantum dots could pave the way to scalable, nuclear- and electron-spin-based quantum computers in silicon that operate without the need for oscillating magnetic fields.

42 ENGINEERING↗
DOE OSTI2025

Solvent Transport in Disordered and Dynamic Membrane Pores: Implications for Reverse Osmosis and Nanofiltration Membranes

Pressure-driven separations with nanoporous membranes, such as reverse osmosis and nanofiltration, play a vital role in addressing water scarcity and enabling resource recovery. Understanding water or solvent transport in membrane pores is essential for advancing membrane separation technologies. A key question in transport modeling is to establish a relationship between solvent permeability and membrane porous structure properties, such as porosity or pore size. The nano- and subnanometer pores in polymeric membranes such as reverse osmosis and nanofiltration membranes are highly tortuous and dynamically connected, which challenges the conventional methods of transport modeling. This study addresses this challenge by developing a theoretical framework to describe solvent transport through membranes with dynamic and disordered porous structure. Specifically, we propose a lattice model to describe the pore network, while preserving the viscous nature of solvent permeation. We further establish a relationship between solvent permeability and membrane porosity or pore size, which is validated by molecular dynamics simulations and experimental data. By integrating this relationship into the solution-friction model, we define pore connectivity and local friction coefficient to quantify the impact of pore structure on solvent permeability. Our analysis highlights the dominant influence of pore connectivity on the permeability of reverse osmosis and nanofiltration membranes, particularly when the pore size approaches the dimensions of solvent molecules. Overall, this study provides critical insights into water and solvent transport mechanisms in nanoporous membranes, opening the door for strategies to substantially enhance membrane performance.

membranes↗
DOE OSTI2022

Characterization of micro-sandwich structures via direct ink writing epoxy based cores

Sandwich structured (SS) composites demonstrate considerable flexural stiffness and high strength-to-weight ratios and can be tailored as functional materials. Historically they have been constrained to specific material types and geometry due to limitations in manufacturing methods. However, employing additive manufacturing (AM), specifically direct ink writing (DIW), can provide an alternative method for making SS composites with complex and controllable micro and mesostructures with multifunctionality targeted at desired mechanical, thermal, and electrical properties. DIW, an extrusion-based AM technique, uses a viscous and thixotropic ink with desired components that, once printed, is cured to obtain the final complex net shape parts. In this paper, a novel hybrid AM technique is employed to manufacture SS composite materials containing bisphenol A-based epoxy core and carbon fiber reinforced polymer (CFRP) face sheets that are fabricated via DIW and vacuum infusion process (VIP), respectfully. We demonstrate that the fabrication of these SS composites can be tailored from a thermosetting material, from which additives and/or various lattice structures can be manufactured to achieve enhanced and desirable mechanical integrity with functional properties. Surface topology and mechanical testing techniques are used to characterize the fabricated hybrid SS composites to study and assess mechanical stability. A rheo-kinetic cure model was developed for the core material to allow for additive manufacturing process requirements while ensuring complete cross-linking for the thermoset-based core material. Because of the ability to obtain relatively small core-thickness and controlled architecture, this method now allows for fabricating layered micro-sandwich structures for realizing further light-weighting in relevant applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Chromium segregation-induced oxide evolution in Ni-10Cr alloys during high-temperature oxidation

The oxidation behavior of a Ni-10(wt%)Cr alloy under high-temperature O 2 conditions is investigated using transmission electron microscopy and first-principles calculations. Results reveal that chromium segregation plays a central role in driving the evolution of complex oxide phase structures during oxidation. At low Cr concentrations, Cr preferentially segregates to NiO grain boundaries or internal pores, substituting for Ni atoms and forming Ni(Cr)O solid solutions. As Cr content increases, enhanced diffusion promotes Cr penetration into the NiO lattice, leading to the formation of multiphase oxide structures. First-principles modeling corroborates these findings: at low Cr concentrations, Cr atoms favor surface and grain-boundary segregation, while higher concentrations lead to Cr aggregation within the NiO bulk. Furthermore, the integrated experimental-theoretical approach provides atomistic insights into Cr-mediated mass transport mechanisms during alloy oxidation and offers valuable guidance for controlling oxide growth kinetics and phase stability in Ni-Cr alloys, with implications for improving oxidation resistance in high-temperature structural applications.

36 MATERIALS SCIENCE↗