SEARCH · Engineering Papers

Results for “Models based on symmetries”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Meta-shell Approach for Constructing Lightweight and High Resolution X-Ray Optics

Lightweight and high resolution optics are needed for future space-based x-ray telescopes to achieve advances in high-energy astrophysics. Past missions such as Chandra and XMM-Newton have achieved excellent angular resolution using a full shell mirror approach. Other missions such as Suzaku and NuSTAR have achieved lightweight mirrors using a segmented approach. This paper describes a new approach, called meta-shells, which combines the fabrication advantages of segmented optics with the alignment advantages of full shell optics. Meta-shells are built by layering overlapping mirror segments onto a central structural shell. The resulting optic has the stiffness and rotational symmetry of a full shell, but with an order of magnitude greater collecting area. Several meta-shells so constructed can be integrated into a large x-ray mirror assembly by proven methods used for Chandra and XMM-Newton. The mirror segments are mounted to the meta-shell using a novel four point semi-kinematic mount. The four point mount deterministically locates the segment in its most performance sensitive degrees of freedom. Extensive analysis has been performed to demonstrate the feasibility of the four point mount and meta-shell approach. A mathematical model of a meta-shell constructed with mirror segments bonded at four points and subject to launch loads has been developed to determine the optimal design parameters, namely bond size, mirror segment span, and number of layers per meta-shell. The parameters of an example 1.3 m diameter mirror assembly are given including the predicted effective area. To verify the mathematical model and support opto-mechanical analysis, a detailed finite element model of a meta-shell was created. Finite element analysis predicts low gravity distortion and low thermal distortion. Recent results are discussed including Structural Thermal Optical Performance (STOP) analysis as well as vibration and shock testing of prototype meta-shells.

Lightweight and high resolution optics↗

DOE OSTI2020

Magnetic ground state of La 2 LiMoO 6 : A comparison with other Mo 5+ ( S=1/2 ) double perovskites

La 2 LiMoO 6 is a double perovskite (DP) with P2 1 /n symmetry based on the Mo5+ ion, 4d1, t2g1, S=1/2. It is isostructural with Sr 2 YMoO 6 , the magnetic ground state of which is apparently a very unusual collective spin singlet or valence-bond glass state as is the case for cubic (Fm-3m) Ba 2 YMoO 6 . Initial studies of La 2 LiMoO 6 suggested a different ground state from the other DPs but no clear conclusions could be drawn. A more detailed study is presented here including magnetic susceptibility, heat capacity, and elastic neutron-scattering results. This DP is now well characterized as an antiferromagnet, T N =18K, via observation of magnetic Bragg peaks in neutron scattering and an anomaly in the magnetic heat capacity. The ordering wave vector is k=( 1/21/2 0), consistent with a type I face-centered-cubic magnetic structure, and the ordered moment on Mo 5+ is 0.32(11)μ B , much reduced from the spin-only value of 1μ B . The index, f=|v c |/T N ~3, indicates a low level of frustration. The heat-capacity data above T N can be interpreted in terms of a one-dimensional spin-correlation model, as can the low-temperature data which follow a T 1 power law. This is consistent with an earlier suggestion. The difference with isostructural Sr 2 YMoO 6 is attributed to differences in the local distortion of the Mo–O octahedron and the resulting orbital ordering.

36 MATERIALS SCIENCE↗

DOE OSTI2022

1 H chemical shift anisotropy: a high sensitivity solid-state NMR dynamics probe for surface studies?

We report dynamics play significant roles in chemistry and biochemistry-molecular motions impact both large- and small-scale chemical reactions in addition to biochemical processes. In many systems, including heterogeneous catalysts, the characterization of dynamics remains a challenge. The most common approaches involve the solid-state NMR measurement of anisotropic interactions, in particular 2 H quadrupolar coupling and 1 H–X dipolar coupling, which generally require isotope enrichment. Due to the high sensitivity of 1 H NMR, 1 H chemical shift anisotropy (CSA) is a particularly enticing, and underexplored, dynamics probe. We carried out 1 H CSA and 1 H– 13 C dipolar coupling measurements in a series of model supported complexes to understand how 1 H CSA can be leveraged to gain dynamic information for heterogeneous catalysts. Mathematical descriptions are given for the dynamic averaging of the CSA tensor, and its dependence on orientation and asymmetry. The variability of the orientation of the tensor in the molecular frame, in addition to its magnitude and asymmetry, negatively impacts attempts to extract quantitative dynamic information. Nevertheless, 1 H CSA measurements can reveal useful qualitative insights into the motions of a particularly dilute site, such as from a surface species.

1H chemical shift anisotropy↗

DOE OSTI2021

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Indirect tunneling enabled spontaneous time-reversal symmetry breaking and Josephson diode effect in TiN/Al 2 ⁢O 3 /Hf 0.8 ⁢Zr 0.2 ⁢O 2 /Nb tunnel junctions

Josephson diode (JD) effect in Josephson tunnel junctions (JTJs) has attracted a great deal of attention due to its importance for developing superconducting-circuitry-based quantum technologies. Even though the preparation of high-quality JTJs by techniques employed in the semiconductor industry has been demonstrated, which was an important milestone because JTJs are the building blocks of superconducting electronics even before the quantum era, the JD effect has not been accomplished in them, nor has the highly desirable electrical control of the effect. We report here the fabrication of JTJs featuring a composite tunnel barrier of Al 2 ⁢O 3 and Hf 0.8 ⁢Zr 0.2 ⁢O 2 using complementary-metal-oxide-semiconductor compatible atomic layer deposition. These JTJs were found to show the JD effect in nominally zero magnetic fields with nonreciprocity controllable via an electric training current, yielding a surprisingly large diode efficiency. The quasiparticle tunneling, through which the Josephson coupling in a JTJ is established, was found to show theoretically expected gap features but no nonreciprocity. We attribute these observations to the simultaneous presence of positive and negative local Josephson couplings in the JTJs, with the negative Josephson coupling originating from indirect tunneling, which results in spontaneous time-reversal symmetry breaking. Finally, the double-minima washboard potential for the ensemble-averaged phase difference in the resistively and capacitively shunted junction model is shown to fully account for the experimentally observed JD effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Control of spin dynamics in artificial honeycomb spin-ice-based nanodisks

In this paper, we report the experimental and theoretical characterization of the angular-dependent spin dynamics in arrays of ferromagnetic nanodisks arranged on a honeycomb lattice. The magnetic field and microwave frequency dependence, measured by broadband ferromagnetic resonance, reveal a rich spectrum of modes that is strongly affected by the microstate of the network. Based on symmetry arguments with respect to the external field, we show that certain parts of the ferromagnetic network contribute to the detected signal. A comparison of the experimental data with micromagnetic simulations reveals that different subsections of the lattice predominantly contribute to the high-frequency response of the array. This is confirmed by optical characterizations using microfocused Brillouin light scattering. Furthermore, we find indications that nucleation and annihilation of vortexlike magnetization configurations in the low-field range affect the dynamics, which is different from clusters of ferromagnetic nanoellipses. Our work opens up new perspectives for designing magnonic devices that combine geometric frustration in gyrotropic vortex crystals at low frequencies with magnonic crystals at high frequencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2021

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT↗

NASA NTRS2021

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

BOLT↗

NASA NTRS2017

Jet Surface Interaction Noise in a High Aspect Ratio Rectangular Exhaust

A physics-based prediction model is employed to simulate jet surface interaction (JSI) noise in a transversely sheared jet exhaust. The methodology finds application in jets with a high aspect ratio (AR) rectangular exhaust in the proximity of a flat surface. Two component spectra are simulated: (i) mixing/scrubbing noise; (ii) trailing edge noise--and are superimposed to obtain the far field exhaust noise on either side of a nearby surface. This document describes the necessary input parameters (including mean flow and turbulence information for the nozzle exhaust of interest) that should be prepared in order to initiate the simulation for each noise component. Sample input/output files in connection with an 8:1 aspect ratio rectangular exhaust at Mach 0.98 near a rigid surface are described. Jet noise spectra are examined below at operating conditions listed in Table IV. Individual noise components, designated as Scrubbing Noise and Trailing Edge Noise, are presented and their sum Total Noise (Analysis) is compared with Measurement (Refs. 8 and 9) at selective number of observer polar angles at azimuth φ = 90deg. Results are presented on an arc R = 17.80-ft (i.e., R = 100Deq) on both sides of a nearby surface. Although the predicted TE noise component is symmetric with respect to the edge due to symmetry in the propagator, measurements for the majority of cases are not quite symmetric and exhibit a slightly larger peak on the reflected side of the surface. Turbulent mixing/scrubbing noise component has a greater presence on the reflected side, as expected. Figure 13 to Figure 18 show that the peak in the predicted TE component could differ from measurements by as much as 4 dB due to lack of symmetry in measured data, however, the general trend is in agreement with data across the three Mach numbers. The overall sound pressure level (OASPL) associated with the TE noise component follows a U5 velocity scaling in the current modeling (Ref. 4). Directivity predictions for the TE noise component as well as the total noise are shown in Figure 19 (bottom)-and are compared with measurements (top figure) at conditions of Table IV. As anticipated, the TE noise component (dashed-line) overwhelms the directivity factor due to its dominant spectral peak level. Only at small angles to the jet axis the mixing noise component contributes significant enough to weight noticeably on the total noise.

propulsion noise↗

DOE OSTI2022

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. Furthermore, we compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Extracting the Quantum Hall Conductance from a Single Bulk Wave Function

In this work, we propose a new formula that extracts the quantum Hall conductance from a single (2+1)⁢ D gapped wave function. The formula applies to general many-body systems that conserve particle number, and is based on the concept of modular flow, i.e., unitary dynamics generated from the entanglement structure of the wave function. The formula is shown to satisfy all formal properties of the Hall conductance: it is odd under time reversal and reflection, even under charge conjugation, universal and topologically rigid in the thermodynamic limit. Further evidence for relating the formula to the Hall conductance is obtained from conformal field theory arguments. Finally, we numerically check the formula by applying it to a noninteracting Chern band where excellent agreement is obtained.

2-dimensional systems↗

DOE OSTI2021

Multi-spin soft bootstrap and scalar-vector Galileon

We use the amplitude soft bootstrap method to explore the space of effective field theories (EFT) of massless vectors and scalars. It is known that demanding vanishing soft limits fixes uniquely a special class of EFTs: non-linear sigma model, scalar Galileon and Born-Infeld theories. Based on the amplitudes analysis, we conjecture no-go theorems for higher-derivative vector theories and theories with coupled vectors and scalars. We then allow for more general soft theorems where the non-trivial part of the soft limit of the ( n +1)-pt amplitude is equal to a linear combination of n-pt amplitudes. We derive the form of these soft theorems for general power-counting and spins of particles and use it as an input into the soft bootstrap method in the case of Galileon power-counting and coupled scalar-vector theories. We show that this unifies the description of existing Galileon theories and leads us to the discovery of a new exceptional theory: Special scalar-vector Galileon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1998

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗

NASA NTRS2013

Tropical Ocean Surface Energy Balance Variability: Linking Weather to Climate Scales

Radiative and turbulent surface exchanges of heat and moisture across the atmosphere-ocean interface are fundamental components of the Earth s energy and water balance. Characterizing the spatiotemporal variability of these exchanges of heat and moisture is critical to understanding the global water and energy cycle variations, quantifying atmosphere-ocean feedbacks, and improving model predictability. These fluxes are integral components to tropical ocean-atmosphere variability; they can drive ocean mixed layer variations and modify the atmospheric boundary layer properties including moist static stability, thereby influencing larger-scale tropical dynamics. Non-parametric cluster-based classification of atmospheric and ocean surface properties has shown an ability to identify coherent weather regimes, each typically associated with similar properties and processes. Using satellite-based observational radiative and turbulent energy flux products, this study investigates the relationship between these weather states and surface energy processes within the context of tropical climate variability. Investigations of surface energy variations accompanying intraseasonal and interannual tropical variability often use composite-based analyses of the mean quantities of interest. Here, a similar compositing technique is employed, but the focus is on the distribution of the heat and moisture fluxes within their weather regimes. Are the observed changes in surface energy components dominated by changes in the frequency of the weather regimes or through changes in the associated fluxes within those regimes? It is this question that the presented work intends to address. The distribution of the surface heat and moisture fluxes is evaluated for both normal and non-normal states. By examining both phases of the climatic oscillations, the symmetry of energy and water cycle responses are considered.

Roberts, J. Brent↗

DOE OSTI2022

Gravitational waves from global cosmic strings and cosmic archaeology

Global cosmic strings are predicted in many motivated extensions to the Standard Model of particle physics, with close connections to axion dark matter physics. Recent studies suggest that, although subdominant relative to Goldstone emission, gravitational wave (GW) signals from global strings can be detectable with current and planned GW detectors such as LIGO, LISA, DECIGO/BBO, ET/CE and AEDGE/AION, as well as pulsar timing arrays such as PPTA, NANOGrav and SKA. This work is an extensive, updated study on GWs from a global cosmic string network, taking into account of the most recent developments related to the subject. The main analysis is based on the analytical Velocity-dependent One-Scale (VOS) model calibrated with recent simulation results, which provides a generic protocol for such calculations with details given. We also demonstrate how the GW signal can be influenced with variations to the baseline model: this includes considering the uncertainties of model parameters and the potential deviation from the conventional VOS model prediction (i.e. the scaling behavior) as suggested by some of the recent simulation results. Furthermore, we investigated in detail the effect of a non-standard cosmology (e.g. early matter domination or kination) or new particle species on the GW signals from global strings. We demonstrate that the frequency spectrum of GW background from global cosmic strings can be used to probe the cosmic history prior to the Big Bang nucleosynthesis (BBN) (i.e. the primordial dark age) up to a temperature of T ~ 10 8 GeV.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Scattering wave packets of hadrons in gauge theories: Preparation on a quantum computer

Quantum simulation holds promise of enabling a complete description of high-energy scattering processes rooted in gauge theories of the Standard Model. A first step in such simulations is preparation of interacting hadronic wave packets. To create the wave packets, one typically resorts to adiabatic evolution to bridge between wave packets in the free theory and those in the interacting theory, rendering the simulation resource intensive. In this work, we construct a wave-packet creation operator directly in the interacting theory to circumvent adiabatic evolution, taking advantage of resource-efficient schemes for ground-state preparation, such as variational quantum eigensolvers. By means of an ansatz for bound mesonic excitations in confining gauge theories, which is subsequently optimized using classical or quantum methods, we show that interacting mesonic wave packets can be created efficiently and accurately using digital quantum algorithms that we develop. Specifically, we obtain high-fidelity mesonic wave packets in the Z 2 and U(1) lattice gauge theories coupled to fermionic matter in 1+1 dimensions. Our method is applicable to both perturbative and non-perturbative regimes of couplings. The wave-packet creation circuit for the case of the Z 2 lattice gauge theory is built and implemented on the Quantinuum H1-1 trapped-ion quantum computer using 13 qubits and up to 308 entangling gates. The fidelities agree well with classical benchmark calculations after employing a simple symmetry-based noise-mitigation technique. This work serves as a step toward quantum computing scattering processes in quantum chromodynamics.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1997

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Barth, Timothy↗

DOE OSTI2025

New state of matter between the hadronic phase and the quark-gluon plasma?

Lattice-quantum chromodynamics (QCD) simulations and theoretical arguments hint at the existence of an intermediate phase of strongly interacting matter between a confined hadron gas and a deconfined quark-gluon plasma (QGP). We qualitatively and semiquantitatively explore and differentiate the phase structures in the temperature window from the QCD pseudocritical temperature 𝑇 c ≃ 160 MeV to the pure gluonic deconfinement temperature 𝑇 d ≃ 285 MeV. We propose a three-regime picture using a hadron resonance gas description augmented with the exponential spectrum of strings, corresponding to highly excited mesons and glueballs, based on the analysis of a large number 𝑁 c of colors. We estimate the entropy density from our model to confirm that the lattice-QCD data are bracketed with three regimes, i.e., a hadron gas, a QGP, and a new phase for 𝑇 c ≲ 𝑇 ≲ 𝑇 d . In this new phase, which we name a spaghetti of quarks with glueballs (SQGBs), thermal degrees of freedom of quarks are liberated, yet gluons remain confined in glueballs. Since the Hagedorn temperature 𝑇 H ∼ 285 MeV is universal in the meson and the glueball sectors, in the infinite-𝑁 c limit, the phase diagram in the plane of the baryon chemical potential and the temperature is reduced to one with the confined and deconfined phases and quarkyonic matter at high density. At large but finite 𝑁 c , an SQGB window may open between these phases. We point out that the SQGB has interesting similarities with quarkyonic matter and that this matter in the large-𝑁 c limit is confined as measured by the interaction between heavy quarks, but behaves in other respects like a quasifree gas of quarks. As a result of the extrapolation to 𝑁 c = 3, we present a revised phase diagram with the SQGB phase bounded by thermal crossovers. Finally, we give a quantitative analysis of chiral-symmetry restoration in the SQGB phase.

Color confinement↗