Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite temperature field theory”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Modeling of three-dimensional mixing and reacting ducted flows

A computer code, based upon a finite element solution algorithm, was developed to solve the governing equations for three-dimensional, reacting boundary region, and constant area ducted flow fields. Effective diffusion coefficients are employed to allow analyses of turbulent, transitional or laminar flows. The code was used to investigate mixing and reacting hydrogen jets injected from multiple orifices, transverse and parallel to a supersonic air stream. Computational results provide a three-dimensional description of velocity, temperature, and species-concentration fields downstream of injection. Experimental data for eight cases covering different injection conditions and geometries were modeled using mixing length theory (MLT). These results were used as a baseline for examining the relative merits of other mixing models. Calculations were made using a two-equation turbulence model (k+d) and comparisons were made between experiment and mixing length theory predictions. The k+d model shows only a slight improvement in predictive capability over MLT. Results of an examination of the effect of tensorial transport coefficients on mass and momentum field distribution are also presented. Solutions demonstrating the ability of the code to model ducted flows and parallel strut injection are presented and discussed.

Zelazny, S. W.↗
DOE OSTI2023

Multipartite entanglement in the one-dimensional spin- 1 2 Heisenberg antiferromagnet

Multipartite entanglement refers to the simultaneous entanglement between multiple subsystems of a many-body quantum system. While multipartite entanglement can be difficult to quantify analytically, it is known that it can be witnessed through the quantum Fisher information (QFI), a quantity that can also be related to dynamical Kubo response functions. In this work, we first show that the finite temperature QFI can generally be expressed in terms of a static structure factor of the system, plus a correction that vanishes as T→0. This implies that the static structure factor witnesses multipartite entanglement near quantum critical points at temperatures below a characteristic energy scale of the system. Therefore, in systems with a known static structure factor, we can deduce finite temperature scaling of multipartite entanglement and low temperature entanglement depth without knowledge of the full dynamical response function of the system. This is particularly useful to study 1D quantum critical systems in which sub-power-law divergences can dominate entanglement growth, where the conventional scaling theory of the QFI breaks down. Furthermore, the 1D spin- 1 2 antiferromagnetic Heisenberg model is an important example of such a system, and we show that multipartite entanglement in the Heisenberg chain diverges nontrivially as ~ln(1/T) 3/2 . We verify these predictions with calculations of the QFI using conformal field theory and matrix product state simulations. Finally, we discuss the implications of our results for experiments to probe entanglement in quantum materials, comparing to neutron scattering data in KCuF 3 , a material well described by the Heisenberg chain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamics of non-Gaussian fluctuations in model A

Motivated by the experimental search for the QCD critical point, we perform simulations of a stochastic field theory with purely relaxational dynamics (model A). We verify the expected dynamic scaling of correlation functions. Using a finite size scaling analysis, we obtain the dynamic critical exponent z = 2.026(56). We investigate time dependent correlation functions of higher moments M n (t) of the order parameter M(t) for n = 1, 2, 3, 4. We obtain dynamic scaling with the same critical exponent z for all n, but the relaxation constant depends on n. We also study the relaxation of M n (t) after a quench, where the simulation is initialized in the high temperature phase, and the dynamics is studied at the critical temperature T c . Finally, we find that the evolution does not follow simple scaling with the dynamic exponent z, and that it involves an early time rise followed by late stage relaxation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Classical analog of quantum models in synthetic dimensions

We introduce a classical analog of quantum matter in ultracold molecule-synthetic or Rydberg atom-synthetic dimensions, by extending the Potts model to include interactions J 1 between atoms adjacent in both real and synthetic space and studying its finite-temperature properties. For intermediate values of J 1 , the resulting phases and phase diagrams are similar to those of the clock and Villain models, in which three phases emerge. There exists a sheet phase analogous to that found in quantum synthetic dimension models between the high-temperature disordered phase and the low-temperature ferromagnetic phase. Furthermore, we also employ machine learning to uncover nontrivial features of the phase diagram using the learning by confusion approach, which is able to discern several successive phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Design and Analysis of an Advanced Three-Point Bend Test Approach for Miniature Irradiated Disk Specimens

Candidate tungsten armor materials in a magnetic confinement fusion device must be able to withstand thermal variation that leads to internal stresses caused by the impinging heat load. In addition, the thermomechanical properties of these materials are degraded by irradiation-induced defect accumulation. Fission reactor–based irradiation data are used to predict the fusion neutron damage and property change. This study examines the motivation and design of a custom-designed three-point bend test for neutron-irradiated disk specimens that are 3 mm in diameter to be able to define the flexural strength of advanced tungsten materials, alloys, and composites—and to the extent that embrittlement occurs after neutron irradiation. The theory provided shows a calculation for the flexural deflection and shear deflection due to the small-geometry constraints. A finite element deformation analysis is performed to evaluate the mechanical stress field of disk bend specimens. The stress values above 80% of the maximum stress are concentrated in 2.4 mm of the 3.0-mm length of the centerline across the tungsten disk diameter. A bend test fixture has been designed and fabricated to enable testing of these specimens with precisely engineered tolerance and minimal machine compliance. This fixture will be able to be placed inside a universal testing frame at elevated temperatures for the mechanical property evaluation of future neutron-irradiated disk specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Modeling prebiotic chemistries with quantum accuracy at classical costs

Molecular Dynamics (MD) simulations using classical force-fields are commonly employed in numerous scientific investigations. However, many natural processes involve bond breaking and quantum forces. This complexity is compounded by the presence of multiple competing length and timescales. For example, accurately modeling the thermodynamics and dynamics of a chemical reaction requires accounting for the concerted movements of numerous solvent molecules and ions with their own fast or slow timescales. While widely used static Density Functional Theory (DFT) calculations at 0 temperature can be beneficial for such investigations, they do not account for dynamics, and lack precision in describing the molecular environments. They particularly fail at correct, rigorous treatments of finite-temperature fluctuations, and thus generalization to experimentally relevant conditions. In PNAS Benayad et al develop a scalable, generalizable approach for designing Neural Network Potentials (NNPs) that can handle chemical reactivity in solvated systems with quantum accuracy at classical costs. Specifically, they study phosphoester bond formation and rupture, which is fundamentally relevant to the Phosphorus-Oxygen bond formation central to life, and especially for the RNA world hypothesis. The framework developed here has the potential to generalize to different chemical reactions of energy and biological relevance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

Modeling of three-dimensional mixing and reacting ducted flows

A computer code based on a finite-element solution algorithm is developed to perform an analytical investigation on the turbulent mixing and reaction of hydrogen jets injected from multiple orifices transverse and parallel to a supersonic airstream. A laser optical cavity flow field was also analyzed to demonstrate the generality of the proposed model. Computational results provide a three-dimensional description of velocity, temperature, and species-concentration fields downstream of injection. Major conclusions are that the analysis has immediate utility in evaluating the mixing effectiveness of transverse H2 injection data since it has been tested in its ability to model this type of data and that turbulent mixing length theory, constant effective Prandtl number, and a Lewis number of unity provide reasonable agreement with transverse H2 injection data downstream of the near-injection region. Efforts are presently being directed toward using the code in modeling laser and scramjet data for a wide range of flow conditions.

Zelazny, S. W.↗
DOE OSTI2024

The Schwinger-Keldysh coset construction

The coset construction is a tool for systematically building low energy effective actions for Nambu-Goldstone modes. This technique is typically used to compute time-ordered correlators appropriate for S-matrix computations for systems in their ground state. In this paper, we extend this technique to the Schwinger-Keldysh formalism, which enables one to calculate a wider variety of correlators and applies also to systems in a mixed state. We focus our attention on internal symmetries and demonstrate that, after identifying the appropriate symmetry breaking pattern, Schwinger-Keldysh effective actions for Nambu-Goldstone modes can be constructed using the standard rules of the coset construction. Particular emphasis is placed on the thermal state and ensuring that correlators satisfy the KMS relation. We also discuss explicitly the power counting scheme underlying our effective actions. We comment on the similarities and differences between our approach and others that have previously appeared in the literature. In particular, our prescription does not require the introduction of additional “diffusive” symmetries and retains the full non-linear structure generated by the coset construction. We conclude with a series of explicit examples, including a computation of the finite-temperature two-point functions of conserved spin currents in non-relativistic paramagnets, antiferromagnets, and ferromagnets. Along the way, we also clarify the discrete symmetries that set antiferromagnets apart from ferromagnets, and point out that the dynamical KMS symmetry must be implemented in different ways in these two systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

Electromagnetic fluctuations for anisotropic media and the generalized Kirchhoff's law

In this paper the polarimetric emission parameters for anisotropic media are derived using the generalized Kirchhoff's law for media with a uniform temperature and the fluctuation-dissipation theory for media with a temperature profile. Both finite-size objects and half-space media are considered. When the object has a uniform temperature across its body, the Kirchhoff's law, based on the condition of energy conservation in thermal equilibrium is generalized to obtain the emission parameters of an anisotropic medium, which can be interpreted as the absorptivity or the absorption cross section of the complementary object with a permittivity that is the transpose of the original object. When the medium has a nonuniform temperature distribution, the fluctuation-dissipation theory is applied for deriving the covariances between vector components of the thermal currents and, consequently, the covariances of the polarizations of electric fields radiated by the thermal currents. To verify the formulas derived from the fluctuation-dissipation theory, we let the temperature of the object be a constant and show that the results reduce to those obtained from the generalized Kirchhoff's law.

Yueh, Simon H.↗
DOE OSTI2022

Characters, quasinormal modes, and Schwinger pairs in dS 2 with flux

An integral representation of the 1-loop partition function for charged scalars and spinors, minimally coupled to a uniform U(1) field on S 2 , is given in terms of SO(1, 2) Harish-Chandra group characters and evaluated exactly in terms of Hurwitz ζ-functions. Analytically continuing the U(1) field, we interpret the path integrals as quasicanonical partition functions in dS 2 with an electric field. The character itself is obtained as a trace over states living at the future boundary of de Sitter and has a quasinormal mode expansion. The imaginary part of the partition function captures Schwinger pair creation in the static patch at finite temperature. The thermal enhancement is most noticeable for scalar masses below Hubble and leads to non-monotonicity of the current as a function of the field. This parameter range, when dimensionally reducing from a charged or rotating Nariai spacetime, is excluded by Swampland-inspired bounds. Around the AdS 2 black hole, in contrast to dS 2 , there is a threshold to pair creation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Relativistic mean-field theories for neutron-star physics based on chiral effective field theory

We describe and implement a procedure for determining the couplings of a relativistic mean-field theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars which consist of neutron-rich matter at densities ranging up to several times saturation density, which allows for additional astrophysical constraints. In our approach, rather than using the RMFT to extrapolate from symmetric to neutron-rich matter and from finite-sized nuclei to uniform matter, we fit the RMFT to properties of uniform pure neutron matter obtained from chiral effective field theory. Chiral effective field theory incorporates the experimental data for nuclei in the framework of a controlled expansion for nuclear forces valid at nuclear densities and enables us to account for theoretical uncertainties when fitting the RMFT. We construct four simple RMFTs that span the uncertainties provided by chiral effective field theory for neutron matter, and are consistent with current astrophysical constraints on the equation of state. Lastly, our RMFTs can be used to model the properties of neutron-rich matter across the vast range of densities and temperatures encountered in neutron stars and their mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 1: Formulation

The present work is part of a research program for the numerical simulation of electromagnetic (EM) fields within conventional Ginzburg-Landau (GL) superconductors. The final goal of this research is to formulate, develop and validate finite element (FE) models that can accurately capture electromagnetic thermal and material phase changes in a superconductor. The formulations presented here are for a time-independent Ginzburg-Landau superconductor and are derived from a potential-based variational principle. We develop an appropriate variational formulation of time-independent supercontivity for the general three-dimensional case and specialize it to the one-dimensional case. Also developed are expressions for the material-dependent parameters alpha and beta of GL theory and their dependence upon the temperature T. The one-dimensional formulation is then discretized for finite element purposes and the first variation of these equations is obtained. The resultant Euler equations contain nonlinear terms in the primary variables. To solve these equations, an incremental-iterative solution method is used. Expressions for the internal force vector, external force vector, loading vector and tangent stiffness matrix are therefore developed for use with the solution procedure.

Schuler, James J.↗
NASA NTRS1978

Prediction of ground effects on aircraft noise

A unified method is recommended for predicting ground effects on noise. This method may be used in flyover noise predictions and in correcting static test-stand data to free-field conditions. The recommendation is based on a review of recent progress in the theory of ground effects and of the experimental evidence which supports this theory. It is shown that a surface wave must be included sometimes in the prediction method. Prediction equations are collected conveniently in a single section of the paper. Methods of measuring ground impedance and the resulting ground-impedance data are also reviewed because the recommended method is based on a locally reactive impedance boundary model. Current practice of estimating ground effects are reviewed and consideration is given to practical problems in applying the recommended method. These problems include finite frequency-band filters, finite source dimension, wind and temperature gradients, and signal incoherence.

Pao, S. P.↗
NASA NTRS1972

Thermal elastoplastic structural analysis of non-metallic thermal protection systems

An incremental theory and numerical procedure to analyze a three-dimensional thermoelastoplastic structure subjected to high temperature, surface heat flux, and volume heat supply as well as mechanical loadings are presented. Heat conduction equations and equilibrium equations are derived by assuming a specific form of incremental free energy, entropy, stresses and heat flux together with the first and second laws of thermodynamics, von Mises yield criteria and Prandtl-Reuss flow rule. The finite element discretization using the linear isotropic three-dimensional element for the space domain and a difference operator corresponding to a linear variation of temperature within a small time increment for the time domain lead to systematic solutions of temperature distribution and displacement and stress fields. Various boundary conditions such as insulated surfaces and convection through uninsulated surface can be easily treated. To demonstrate effectiveness of the present formulation a number of example problems are presented.

Chung, T. J.↗
DOE OSTI2024

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Signature of magnetoelectric coupling driven finite momentum pairing in 3D ising superconductor

The finite momentum superconducting paring states (FMPs) represent a forefront of condensed matter physics. Here we report experimental evidence of FMP in a locally noncentrosymmetric bulk superconductor 4H b -TaS 2 . Using hard X-ray diffraction and angle-resolved photoemission spectroscopy, we reveal unusual 2D ferro-rotational charge density wave (CDW) and weak interlayer hopping in 4H b -TaS 2 . The superconducting upper critical field, H c2 , linearly increases via decreasing temperature, and well exceeds the Pauli limit, suggesting the dominant orbital pair-breaking mechanism. Remarkably, we observed evidence of field-induced superconductivity-to-superconductivity transition that breaks continuous rotational symmetry of the s-wave uniform pairing in the Bardeen-Cooper-Schrieffer theory down to the six-fold rotation symmetry. Ginzburg-Landau free energy analysis shows that magnetoelectric coupling, induced by 2D ferro-rotational CDW, stabilizes FMP that provides an explanation of the lowering rotation symmetry. Our results provide a new understanding of unconventional superconducting behaviors of the bulk quantum heterostructure 4H b -TaS 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗