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 199 records · Page 11

Tensor train continuous time solver for quantum impurity models

The simulation of strongly correlated quantum impurity models is a significant challenge in modern condensed matter physics that has multiple important applications. Thus far, the most successful methods for approaching this challenge involve Monte Carlo techniques that accurately and reliably sample perturbative expansions to any order. However, the cost of obtaining high precision through these methods is high. Recently, tensor train decomposition techniques have been developed as an alternative to Monte Carlo integration. In this study, we apply these techniques to the single-impurity Anderson model at equilibrium by calculating the systematic expansion in power of the hybridization of the impurity with the bath. Furthermore, we demonstrate the performance of the method in a paradigmatic application, examining the first-order phase transition on the infinite-dimensional Bethe lattice, which can be mapped to an impurity model through dynamical mean field theory. Our results indicate that using tensor train decomposition schemes allows the calculation of finite-temperature Green's functions and thermodynamic observables with unprecedented accuracy. The methodology holds promise for future applications to frustrated multiorbital systems, using a combination of partially summed series with other techniques pioneered in diagrammatic and continuous time quantum Monte Carlo.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parametrization of Generalized Parton Distributions from 𝑡-Channel String Exchange in AdS Spaces

We introduce a string-based parametrization for nucleon quark and gluon generalized parton distributions (GPDs) that is valid for all skewness. Our approach leverages conformal moments, representing them as the sum of spin-𝑗 nucleon 𝐴-form factor and skewness-dependent spin-𝑗 nucleon 𝐷-form factor, derived from 𝑡-channel string exchange in AdS spaces consistent with Lorentz invariance and unitarity. This model-independent framework, satisfying the polynomiality condition due to Lorentz invariance, uses Mellin moments from empirical data to estimate these form factors. With just five Regge slope parameters, our method accurately produces various nucleon quark GPD types and symmetric nucleon gluon GPDs through pertinent Mellin-Barnes integrals. Our isovector nucleon quark GPD is in agreement with existing lattice data, promising to improve the empirical extraction and global analysis of nucleon GPDs in exclusive processes, by avoiding the deconvolution problem at any skewness, for the first time.

QCD phenomenology↗

Low-G Slosh Workshop Results From 2023: State of the Art, Gaps and Forward Work

The NASA Engineering and Safety Center (NESC) sponsored a workshop to dis-cuss low-g slosh. Over a dozen different NASA centers and industry partners participated in this workshop. The main topics were to identify the state of the art regarding test data, Computational Fluid Dynamics (CFD) analysis resources, and GN&C models that are available for this problem. The primary intent was to determine deficiencies in the collective knowledge base and to identify priorities that should be established for proactive work in the areas of model development, validation, and model integration. Recommendations from workshop participants include establishing an online slosh modeling resource, conducting ground/flight testing, developing surrogate models, and establishing co-simulation techniques. This paper summarizes the workshop findings and recommendations for the larger community.

Lattice-Boltzmann↗

Tight-Binding Model Study of Anti-ferromagnetic Order in AA-Stacked Bi-layer Graphene

We address here the anti-ferromagnetic order present in AA-stacked bi-layer graphene in a transversely applied electric field. The system is described by kinetic energy with nearest-neighbor electron hopping with same hopping integral t{sub 1} for both the layers. Besides this, Coulomb interaction exists at A and B sub-lattices with same Coulomb correlation energy. The electron Green’s functions are calculated by Zubarev’s Green’s technique. The temperature-dependent anti-ferromagnetic magnetization is calculated from the Green’s function and is computed numerically and self-consistently. The strong on-site Coulomb interaction stabilizes the anti-ferromagnetic order in graphene. We assume that the electron spin at A site in the first layer is directed in the opposite direction to that of A site electron in the second layer. Similar spin order is observed for electrons in B site atom in reversed order. It is observed that anti-ferromagnetic (AFM) magnetization in the first layer nearly remains constant up to certain temperature and then increases with temperature, while the AFM magnetization in the second layer remains nearly constant and then rapidly decreases with temperature. The net AFM magnetization in bi-layer graphene remains constant and then rapidly increases with temperature. The evolution of the AFM magnetization is studied by varying transverse electric field, Coulomb energy, and temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗

An experimental and theoretical analysis of the aerodynamic characteristics of a biplane-winglet configuration

Improving the aerodynamic characteristics of an airplane with respect to maximizing lift and minimizing induced and parasite drag are of primary importance in designing lighter, faster, and more efficient aircraft. Previous research has shown that a properly designed biplane wing system can perform superiorly to an equivalent monoplane system with regard to maximizing the lift-to-drag ratio and efficiency factor. Biplanes offer several potential advantages over equivalent monoplanes, such as a 60-percent reduction in weight, greater structural integrity, and increased roll response. The purpose of this research is to examine, both theoretically and experimentally, the possibility of further improving the aerodynamic characteristics of the biplanes configuration by adding winglets. Theoretical predictions were carried out utilizing vortex-lattice theory, which is a numerical method based on potential flow theory. Experimental data were obtained by testing a model in the Pennsylvania State University's subsonic wind tunnel at a Reynolds number of 510,000. The results showed that the addition of winglets improved the performance of the biplane with respect to increasing the lift-curve slope, increasing the maximum lift coefficient, increasing the efficiency factor, and decreasing the induced drag. A listing of the program is included in the Appendix.

Gall, P. D.↗

Prediction of Broadband Blade-Wake Back-Scatter Noise from a Hovering Ideally Twisted Rotor using OVERFLOW2-ANOPP2

This work illustrates the use of a hybrid RANS/LES framework for the prediction of broadband noise via integral formulations from an ideally twisted hovering rotor. Four simulations were conducted, which investigated the effect of different turbulence modeling and spatiotemporal resolutions. The aerodynamic performance and aeroacoustics predicted by these four computational cases were compared against data acquired in the Small Hover Anechoic Chamber at the NASA Langley Research Center and against previous simulations conducted using the lattice-Boltzmann solver, PowerFLOW. The simulation case with the finest spatial resolution had the most favorable aerodynamic performance and broadband noise comparison with measured data; however, the simulation case utilizing the SST-DDES turbulence model predicted the most accurate tonal noise signature. Aerodynamic flowfield visualizations identified near-wake trailing edge vortex shedding, perpendicular blade-vortex interactions, and tip vortex formation, which are known to cause broadband noise. The fluid dynamics associated with a new broadband noise source, blade-wake back-scatter, were also identified. Lastly, a broadband noise source visualization technique was utilized to confirm the presence of blade-wake back-scatter as a dominant broadband noise source.

blade-wake back-scatter↗

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

Perturbative Stability and Error-Correction Thresholds of Quantum Codes

Topologically ordered phases are stable to local perturbations, and topological quantum error-correcting codes enjoy thresholds to local errors. We connect the two notions of stability by constructing classical statistical mechanics models for decoding general Calderbank-Shor-Steane codes and classical linear codes. Our construction encodes correction success probabilities under uncorrelated bit-flip and phase-flip errors, and simultaneously describes a generalized ℤ 2 lattice-gauge theory with quenched disorder. We observe that the clean limit of the latter is precisely the discretized imaginary-time path integral of the corresponding quantum code Hamiltonian when the errors are turned into a perturbative 𝑋 or 𝑍 magnetic field. Motivated by error-correction considerations, we define general order parameters for all such generalized ℤ 2 lattice-gauge theories, and show that they are generally lower bounded by success probabilities of error correction. For CSS codes satisfying the low-density parity-check condition and with a sufficiently large code distance, we prove the existence of a low-temperature ordered phase of the corresponding lattice-gauge theories, particularly for those lacking Euclidean spatial locality and/or when there is a nonzero code rate. We further argue that these results provide evidence for stable phases in the corresponding perturbed quantum Hamiltonians, obtained in the limit of continuous imaginary time. To do so, we distinguish space- and timelike defects in the lattice-gauge theory. A high free-energy cost of spacelike defects corresponds to a successful “memory experiment” and suppresses the energy splitting among the ground states, while a high free-energy cost of timelike defects corresponds to a successful “stability experiment” and points to a nonzero gap to local excitations.

quantum error correction↗

Domain epitaxial growth of ferroelectric films of barium strontium titanate on sapphire

A model of the epitaxial growth of crystalline multicomponent films on single-crystal substrates with a domain correspondence is presented using a solid solution of barium strontium titanate on sapphire substrates (r cut). The domain epitaxial growth suggests the matching of the lattice planes of the film and the substrate having similar structures by comparison of domain multiple of an integral number of the interplanar spacings. Variation of the component composition of the solid solution enables changes in the domain size in the range sufficient for epitaxial growth. This method can be used to project the epitaxial growth of films of various solid solutions on single-crystal substrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Aeroelastic modeling for the FIT team F/A-18 simulation

Some details of the aeroelastic modeling of the F/A-18 aircraft done for the Functional Integration Technology (FIT) team's research in integrated dynamics modeling and how these are combined with the FIT team's integrated dynamics model are described. Also described are mean axis corrections to elastic modes, the addition of nonlinear inertial coupling terms into the equations of motion, and the calculation of internal loads time histories using the integrated dynamics model in a batch simulation program. A video tape made of a loads time history animation was included as a part of the oral presentation. Also discussed is work done in one of the areas of unsteady aerodynamic modeling identified as needing improvement, specifically, in correction factor methodologies for improving the accuracy of stability derivatives calculated with a doublet lattice code.

Zeiler, Thomas A.↗

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗

Effects of 3d electron configurations on helium bubble formation and void swelling in concentrated solid-solution alloys

Currently, elemental specific chemical complexity is known to play a vital role in microstructure development in single-phase concentrated solid-solution alloys (SP-CSAs), including both He bubble formation and irradiation-induced void swelling. While cavity formation and evolution under ion irradiation at elevated temperature are complex nonequilibrium processes, chemical effects are revealed at the level of electrons and atoms herein in a simplified picture, using Ni and a special set of Ni-based SP-CSAs composed of 3d transition metals as model alloys. Based on Ni and the model alloys with minimized variables (e.g., atomic mass, size, and lattice structure), we report on the effects of chemically-biased energy dissipation, defect energetics, sluggish diffusion, and atomic transport on cavity formation and evolution under both self-ion Ni irradiation and He implantation. The observed difference in microstructure evolution is attributed to the effects of d electron interactions in their integrated ability to dissipate radiation energy. The exhibited impact of alloying 3d transition metals with larger differences in the outermost electron counts suggests a simple design strategy for tuning defect properties to improve radiation tolerance in structural alloys.

36 MATERIALS SCIENCE↗

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗

Coherent, atomically thin transition-metal dichalcogenide superlattices with engineered strain

Epitaxy forms the basis of modern electronics and optoelectronics. We report coherent atomically thin superlattices in which different transition metal dichalcogenide monolayers—despite large lattice mismatches—are repeated and laterally integrated without dislocations within the monolayer plane. Grown by an omnidirectional epitaxy, these superlattices display fully matched lattice constants across heterointerfaces while maintaining an isotropic lattice structure and triangular symmetry. This strong epitaxial strain is precisely engineered via the nanoscale supercell dimensions, thereby enabling broad tuning of the optical properties and producing photoluminescence peak shifts as large as 250 millielectron volts. We present theoretical models to explain this coherent growth and the energetic interplay governing the ripple formation in these strained monolayers. Such coherent superlattices provide building blocks with targeted functionalities at the atomically thin limit.

Science & Technology - Other Topics↗

Theory of d-Wave High Temperature Superconductivity in the Cuprates Involving Non-linear Lattice Modes

The transition mechanism in high temperature cuprate superconductors is an outstanding puzzle. A previous suggestion on the role of non-linear local lattice instability modes on the microscopic pairing mechanism in high temperature cuprate superconductors (Lee, J. Supercond. Nov. Magn. 23(3), 333; 2009) is re-examined to provide a viable mechanism for superconductivity in these cuprates via an unusual lattice vibration in which an electron is predominantly interacting with a non-linear Q{sub 2} mode of the oxygen clusters in the CuO{sub 2} planes. It is shown that the interaction has explicit d-wave symmetry and leads to an indirect coupling of d-wave symmetry between electrons. As a follow-up of Lee (J. Supercond. Nov. Magn. 23(3), 333; 2009), in this paper, we report detailed derivation of the superconducting gap equation and numerical solutions for the transition temperature as inherently integrated into the so-called extended Hubbard model (EHM). A unique feature in the EHM is that the transition temperature has an inherent k-dependence. In addition, superconducting gap solutions are restrained to specific regions in the first Brillouin zone (1BZ). It is very feasible to expect that the EHM naturally inherits a huge parameter space in which experimentally measured results, such as the well-known superconducting dome and the phase diagram from electronic Raman scattering (Sacuto et al., Rep. Prog. Phys. 76(2), 022502; 2013) can be accommodated. The EHM model hence offers a viable venue to search for or confirm any signature in k-point-sensitive experimental measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modelling intergranular and transgranular micro-cracking in polycrystalline materials

Presented is a grain boundary formulation for intergranular and transgranular micro-cracking in three-dimensional polycrystalline aggregates.. The formulation is based on the displacement and stress boundary integral equations of solid mechanics and it has the advantage of expressing the polycrystalline problem in terms of grain boundary variables only. The individual grains within the polycrystalline morphology are modelled as generally anisotropic linear elastic domains with random spatial orientation. Transgranular micro-cracking is assumed to occur along specific cleavage planes, whose orientation in space within the grains depend upon the crystallographic lattice. Both intergranular and transgranular micro-cracking are modelled using suitably defined cohesive laws, whose parameters characterise the behaviour of the two mechanisms. The algorithm developed to track the inter/transgranular micro-cracking history is presented and discussed. Several numerical tests involving pseudo-3D and fully 3D morphologies are performed and analysed. The presented numerical results show that the developed formulation is capable of tracking the initiation and evolution of both intergranular and transgranular cracking as well as their competition, thus providing a useful tool for the study of damage micro-mechanics.

36 MATERIALS SCIENCE↗

Low-temperature grapho-epitaxial La-substituted BiFeO 3 on metallic perovskite

Bismuth ferrite has garnered considerable attention as a promising candidate for magnetoelectric spin-orbit coupled logic-in-memory. As model systems, epitaxial BiFeO 3 thin films have typically been deposited at relatively high temperatures (650–800 °C), higher than allowed for direct integration with silicon-CMOS platforms. Here, we circumvent this problem by growing lanthanum-substituted BiFeO 3 at 450 °C (which is reasonably compatible with silicon-CMOS integration) on epitaxial BaPb 0.75 Bi 0.25 O 3 electrodes. Notwithstanding the large lattice mismatch between the La-BiFeO 3 , BaPb 0.75 Bi 0.25 O 3 , and SrTiO 3 (001) substrates, all the layers in the heterostructures are well ordered with a [001] texture. Polarization mapping using atomic resolution STEM imaging and vector mapping established the short-range polarization ordering in the low temperature grown La-BiFeO 3 . Current-voltage, pulsed-switching, fatigue, and retention measurements follow the characteristic behavior of high-temperature grown La-BiFeO 3 , where SrRuO 3 typically serves as the metallic electrode. These results provide a possible route for realizing epitaxial multiferroics on complex-oxide buffer layers at low temperatures and opens the door for potential silicon-CMOS integration.

36 MATERIALS SCIENCE↗