Search NASA⌕ Search

SEARCH · Search NASA

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermalization in the one-dimensional Salerno model lattice

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Continuous-Time Multidimensional Walks as an Integrable Model

We consider continuous-time random walks on multidimensional symplectic lattices. It is shown that the generating functions of random walks and the transition amplitudes of continuous-time quantum walks can be expressed through dynamical correlation functions of an exactly solvable model describing strongly correlated bosons on a chain, the so-called phase model. The number of random lattice paths with a fixed number of steps connecting the starting and ending points on the multidimensional lattice is expressed through solutions of the Bethe equations of the phase model. Its asymptotics is obtained in the limit of a large number of steps.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Low length scale modeling of Cr precipitation in FeCrAl under thermal aging conditions

This report summarizes the progress on lower length scale modeling of Cr precipitation in FeCrAl alloys, which has been identified as a possible cladding material for accident tolerant fuels due to its excellent corrosion resistance via the formation of a protective, slow growing alumina surface layer. The precipitation of Cr-rich α' clusters in FeCrAl can lead to degraded mechanical properties such as the well-known “475oC embrittlement”, and should therefore be avoided to ensure cladding integrity. In FY16, a lattice kinetic Monte Carlo (LKMC) model was developed for FeCr binary alloys. Due to the assumption of constant bond energies, the previous model cannot correctly reproduce the low-temperature part of the Fe-Cr phase diagram. It is also not possible to model the anomalous ordering tendency and repulsive interactions between Cr atoms in α Fe at the dilute concentration regime using constant bond energies. To address these issues, a new local-environment-dependent bond energy formalism has been developed for Fe-Cr in FY17 that can simultaneously reproduce the negative mixing energy at low Cr concentrations and the positive mixing energy at high Cr concentrations, and the consequent large solubility of Cr in bcc Fe at low temperatures. Furthermore, the LKMC model has been extended to include Al to enable the assessment of Cr precipitation tendency in ternary FeCrAl alloys under thermal aging. Results from the present KMC simulations show good agreement with available experiments in the literature. This study is supported by the DOE NEAMS program Accident Tolerant Fuel project.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fermionic approach to variational quantum simulation of Kitaev spin models

We use the variational quantum eigensolver (VQE) to simulate Kitaev spin models with and without integrability breaking perturbations, focusing in particular on the honeycomb and square-octagon lattices. These models are well known for being exactly solvable in a certain parameter regime via a mapping to free fermions. We use classical simulations to explore a novel variational ansatz that takes advantage of this fermionic representation and is capable of expressing the exact ground state in the solvable limit. We also demonstrate that this ansatz can be extended beyond this limit to provide excellent accuracy when compared to other VQE approaches. In certain cases, this fermionic representation is advantageous because it reduces by a factor of two the number of qubits required to perform the simulation. We also comment on the implications of our results for simulating non-Abelian anyons on quantum computers.

Jahin, Ammar↗
DOE OSTI2025

Quantum annealing-assisted lattice optimization

High Entropy Alloys (HEAs) have drawn great interest due to their exceptional properties compared to conventional materials. The configuration of HEA system is considered a key to their superior properties, but exhausting all possible configurations of atom coordinates and species to find the ground energy state is extremely challenging. In this work, we proposed a quantum annealing-assisted lattice optimization (QALO) algorithm, which is an active learning framework that integrates the Field-aware Factorization Machine (FFM) as the surrogate model for lattice energy prediction, Quantum Annealing (QA) as an optimizer and Machine Learning Potential (MLP) for ground truth energy calculation. By applying our algorithm to the NbMoTaW alloy, we reproduced the Nb depletion and W enrichment observed in bulk HEA. We found our optimized HEAs to have superior mechanical properties compared to the randomly generated alloy configurations. Our algorithm highlights the potential of quantum computing in materials design and discovery, laying a foundation for further exploring and optimizing structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hierarchical excitations from correlated spin tetrahedra on the breathing pyrochlore lattice

The hierarchy of the coupling strengths in a physical system often engenders an effective model at low energies where the decoupled high-energy modes are integrated out. Here, using neutron scattering, we show that the spin excitations in the breathing pyrochlore lattice compound CuInCr 4 S 8 are hierarchical and can be approximated by an effective model of correlated tetrahedra at low energies. At higher energies, intratetrahedron excitations together with strong magnon-phonon couplings are observed, which suggests the possible role of the lattice degree of freedom in stabilizing the spin tetrahedra. Finally, our work illustrates the spin dynamics in CuInCr 4 S 8 and demonstrates a general effective-cluster approach to understand the dynamics on the breathing-type lattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
DOE OSTI2019

QuSpin: a Python package for dynamics and exact diagonalisation of quantum many body systems. Part II: bosons, fermions and higher spins

We present a major update to QuSpin, SciPostPhys.2.1.003 – an open-source Python package for exact diagonalization and quantum dynamics of arbitrary boson, fermion and spin many-body systems, supporting the use of various (user-defined) symmetries in one and higher dimension and (imaginary) time evolution following a user-specified driving protocol. We explain how to use the new features of QuSpin using seven detailed examples of various complexity: (i) the transverse-field Ising chain and the Jordan-Wigner transformation, (ii) free particle systems: the Su-Schrieffer-Heeger (SSH) model, (iii) the many-body localized 1D Fermi-Hubbard model, (iv) the Bose-Hubbard model in a ladder geometry, (v) nonlinear (imaginary) time evolution and the Gross-Pitaevskii equation on a 1D lattice, (vi) integrability breaking and thermalizing dynamics in the translationally-invariant 2D transverse-field Ising model, and (vii) out-of-equilibrium Bose-Fermi mixtures. This easily accessible and user-friendly package can serve various purposes, including educational and cutting-edge experimental and theoretical research. The complete package documentation is available under http://weinbe58.github.io/QuSpin/index.html.

Weinberg, Phillip↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2018

Separating the effects of nucleotide and EB binding on microtubule structure

Microtubules (MTs) are polymers assembled from αβ-tubulin heterodimers that display the hallmark behavior of dynamic instability. MT dynamics are driven by GTP hydrolysis within the MT lattice, and are highly regulated by a number of MT-associated proteins (MAPs). How MAPs affect MTs is still not fully understood, partly due to a lack of high-resolution structural data on undecorated MTs, which need to serve as a baseline for further comparisons. Here we report three structures of MTs in different nucleotide states (GMPCPP, GDP, and GTPγS) at near-atomic resolution and in the absence of any binding proteins. These structures allowed us to differentiate the effects of nucleotide state versus MAP binding on MT structure. Kinesin binding has a small effect on the extended, GMPCPP-bound lattice, but hardly affects the compacted GDP-MT lattice, while binding of end-binding (EB) proteins can induce lattice compaction (together with lattice twist) in MTs that were initially in an extended and more stable state. We propose a MT lattice-centric model in which the MT lattice serves as a platform that integrates internal tubulin signals, such as nucleotide state, with outside signals, such as binding of MAPs or mechanical forces, resulting in global lattice rearrangements that in turn affect the affinity of other MT partners and result in the exquisite regulation of MT dynamics.

Science & Technology - Other Topics↗
NASA NTRSDate not supplied

Aerodynamic Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft with Variable-Camber Continuous Trailing-Edge Flap

This paper presents an aerodynamic optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft with Variable-Camber Continuous Trailing-Edge Flap (VC-CTEF). The VCCTEF is a novel wing shaping control concept to improve aircraft aerodynamic efficiency. Drag reduction studies are conducted for three different VCCTEF configurations with 6-, 10-, and 21-spanwise sections, respectively. A simple VCCTEF actuator weight model is proposed to account the weight penalty of the actuator in the design. A vortex-lattice model of the TTBW model is developed with transonic small disturbance, integral boundary-layer, and wing-strut interference corrections for rapid aerodynamic performance evaluations. The VSPAERO model has been validated against wind tunnel test data. The optimization results show that the 21-spanwise sections VCCTEF has largest absolute drag coefficient reduction, however, the 6-spanwise sections VCCTEF provides a relatively better solution for drag reduction when the actuator weight penalty is considered. A high-fidelity CFD solver FUN3D is used to validate the VCCTEF optimization.

Aerodynamic Optimization↗
DOE OSTI2026

Parametric reduced order models for graded lattice structures

Graded lattice structures, characterized by smoothly varying mechanical properties, hold significant promise for optimizing material distribution in advanced engineering applications. However, accurately modeling these structures poses substantial computational challenges due to the continuous geometric variations within their unit cells. Here, to address these challenges, this paper introduces a novel Efficient Reduced Order Model (EROM) that integrates the Matrix Discrete Empirical Interpolation Method (MDEIM) and Discrete Empirical Interpolation Method (DEIM) with polynomial regression to efficiently manage geometric parametrization in lattice structures. Unlike traditional reduced order models (ROMs) that require extensive precomputed libraries for each geometric configuration, our approach enables continuous geometric variations through a flexible algebraic formulation, significantly reducing computational costs while preserving high accuracy. The method constructs projection matrices for individual unit cells that can be efficiently assembled into global systems, leveraging the repetitive nature of lattice structures. Numerical studies demonstrate that our EROM achieves displacement errors below 1% and von Mises stress prediction errors below 4%, coupled with computational speedups exceeding two orders of magnitude compared to full-order simulations. The proposed method's modularity and scalability make it particularly suitable for design optimization and real-time simulation of functionally graded lattice structures, with applications spanning aerospace to biomedical engineering.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Nonlinear sigma models with compact hyperbolic target spaces

We explore the phase structure of nonlinear sigma models with target spaces corresponding to compact quotients of hyperbolic space, focusing on the case of a hyperbolic genus-2 Riemann surface. The continuum theory of these models can be approximated by a lattice spin system which we simulate using Monte Carlo methods. The target space possesses interesting geometric and topological properties which are reflected in novel features of the sigma model. In particular, we observe a topological phase transition at a critical temperature, above which vortices proliferate, reminiscent of the Kosterlitz-Thouless phase transition in the O(2) model [1, 2]. Unlike in the O(2) case, there are many different types of vortices, suggesting a possible analogy to the Hagedorn treatment of statistical mechanics of a proliferating number of hadron species. Below the critical temperature the spins cluster around six special points in the target space known as Weierstrass points. In conclusion, the diversity of compact hyperbolic manifolds suggests that our model is only the simplest example of a broad class of statistical mechanical models whose main features can be understood essentially in geometric terms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Cumulant methods for electron-phonon problems. II. The self-consistent cumulant expansion

Here in this paper we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattice sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demonstration that the SC-CE corrects many deficits of the standard perturbative CE in the (nondispersive) Holstein model. The natural phonon damping that arises in the more realistic dispersed Holstein model renders the SC-CE stable, allowing for a complete assessment of the method. Here, we find that self-consistency dramatically corrects many of the failures found in the perturbative CE, but also introduces some unphysical features. Finally, we comment on the potential use of SC-CE as a tool for calculating Green's functions in generic many-body problems.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Long wavelength perfect fluidity from short distance jet transport in quark-gluon plasmas

Here, we build a new phenomenological framework that bridges the long wavelength bulk viscous transport properties of the strongly-coupled quark-gluon plasma (sQGP) and short distance hard jet transport properties in the QGP. The full nonperturbative chromo-electric (E) and chromo-magnetic (M) structure of the near “perfect fluid” like sQGP in the critical transition region are integrated into a semi-Quark-Gluon-Monopole Plasma (sQGMP) model lattice-compatibly and implemented into the new CUJET3.0 jet quenching framework. All observables computed from CUJET3.0 are found to be consistent with available data at RHIC and LHC simultaneously. Moreover, a quantitative connection between the shear viscosity and jet transport parameter is rigorously established within this framework. Finally, we deduce the T = 160-600 MeV dependence of the QGP’s η/s: its near vanishing value in the near T c regime is determined by the composition of E and M charges, it increases as T rises, and its high T limit is fixed by color screening scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRSDate not supplied

Study of Trailing Edge Mini Flaps Interactions on a Wind Tunnel Scale Aspect Ratio 13.5 Common Research Model

This paper presents a computation study of an aspect ratio 13.5 wind tunnel scale Common Research Model (CRM) with trailing edge mini-flaps. A vortex-lattice model of the CRM is developed with transonic small disturbance and integral boundary-layer corrections coupled with NASTRAN equivalent beam model for the rapid CRM aeroelastic simulation. The simulations are performed with individual trailing edge mini-flap deflected and a series of groups of two adjacent trailing edge mini-flaps deflected at the same time. A surrogate model is developed to estimate the aerodynamic interactions between the flaps. To validate the surrogate model, static aeroelastic simulations are performed using FUN3D CFD solver.

Interaction↗