SEARCH · Engineering Papers
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.
The flowfield of a rectangular jet with 2:1 aspect ratio was studied at an axial Reynolds number of 127,000, using a three-dimensional laser anemometer. The flowfield surveys resulted in mean velocity vector field plots and contour plots of the Reynolds stress tensor components for the major and minor axes. These data contribute substantially to currently available data of jet flowfields.
This work describes the development of a model that accounts for the additional total surface recession in Thermal Protection Systems materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics capability module was integrated within the Porous material Analysis Toolbox based on OpenFOAM, PATO. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver module computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model of the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.
This work describes the development of a model that accounts for the additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics module was integrated within the Porous material Analysis Toolbox (PATO) based on OpenFOAM. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model capable to predict the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.
This work describes the development of a model that accounts for the additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics module was integrated within the Porous material Analysis Toolbox (PATO) based on OpenFOAM. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model capable to predict the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.
Sparse tensor contraction sequence has been widely employed in many fields, such as chemistry and physics. However, how to efficiently implement the sequence faces multiple challenges, such as redundant computations and memory operations, massive memory consumption, and inefficient utilization of hardware. To address the above challenges, we introduce Athena, a high-performance framework for SpTC sequences. Athena introduces new data structures, leverages emerging Optane-based heterogeneous memory (HM) architecture, and stage parallelism. In particular, Athena introduces shared hash table-represented sparse accumulator to eliminate unnecessary input processing and data migration; Athena uses a novel data-semantic guided dynamic migration solution to make the best use of the Optane-based HM for high performance; Athena also co-runs execution phases with different characteristics to enable high hardware utilization. Evaluating with 12 datasets, we show that Athena brings 327-7362× speedup over the state-of-the-art SpTC algorithm. With the dynamic data placement guided by data semantics, Athena brings performance improvement on Optane-based HM over a state-of-the-art software-based data management solution, a hardware-based data management solution, and PMM-only by 1.58×, 1.82×, and 2.34× respectively.
Chiral Weyl fermions with linear energy-momentum dispersion in the bulk accompanied by Fermi-arc states on the surfaces prompt a host of enticing optical effects. While new Weyl semimetal materials keep emerging, the available optical probes are limited. In particular, isolating bulk and surface electrodynamics in Weyl conductors remains a challenge. We devised an approach to the problem based on near-field photocurrent imaging at the nanoscale and applied this technique to a prototypical Weyl semimetal TaIrTe 4 . As a first step, we visualized nano-photocurrent patterns in real space and demonstrated their connection to bulk nonlinear conductivity tensors through extensive modeling augmented with density functional theory calculations. Notably, our nanoscale probe gives access to not only the in-plane but also the out-of-plane electric fields so that it is feasible to interrogate all allowed nonlinear tensors including those that remained dormant in conventional far-field optics. Surface- and bulk-related nonlinear contributions are distinguished through their “symmetry fingerprints” in the photocurrent maps. Robust photocurrents also appear at mirror-symmetry breaking edges of TaIrTe 4 single crystals that we assign to nonlinear conductivity tensors forbidden in the bulk. Here, nano-photocurrent spectroscopy at the boundary reveals a strong resonance structure absent in the interior of the sample, providing evidence for elusive surface states.
The energy-momentum tensor (EMT) of the proton encodes fundamental information about its mass, pressure, and shear distributions. Using recent lattice QCD data for the gravitational form factors, we show that the Breit-frame Wigner EMT may be of Hawking-Ellis type IV in the proton’s core. Such EMT violates all pointwise energy conditions and lacks a causal rest frame so that the usual mechanical picture fails at short distances. We define the —a new hadronic observable—marking the scale where the EMT becomes ordinary (type I) and the classical interpretation is restored. We also derive from the averaged null energy condition nonperturbative, model-independent quantum field theory constraints on gravitational form factors.
By applying the stress-tensor-scalar operator product expansion (OPE) twice, we search for algebraic structures in d = 4 conformal field theories (CFTs) with a pure Einstein gravity dual. We find that a rescaled mode operator defined by an integral of the stress tensor T ++ on a d = 2 plane satisfies a Virasoro-like algebra when the dimension of the scalar is large. The structure is enhanced to include a Kac-Moody-type algebra if we incorporate the T -- component. In our scheme, the central terms are finite. It remains challenging to directly compute the stress-tensor sector of d = 4 scalar four-point functions at large central charge, which, based on holography and bootstrap methods, were recently shown to have a Virasoro/ $\mathcal{W}$-algebra vacuum block-like structure.
We construct “soft-collinear gravity”, the effective field theory which describes the interaction of collinear and soft gravitons with matter (and themselves), to all orders in the soft-collinear power expansion. Despite the absence of collinear divergences in gravity at leading power, the construction exhibits remarkable similarities with soft-collinear effective theory of QCD (gauge fields). It reveals an emergent soft background gauge symmetry, which allows for a manifestly gauge-invariant representation of the interactions in terms of a soft covariant derivative, the soft Riemann tensor, and a covariant generalisation of the collinear light-cone gauge metric field. The gauge symmetries control both the unsuppressed collinear field components and the inherent inhomogeneity in λ of the invariant objects to all orders, resulting in a consistent expansion.
Abstract We review the basic ideas of the tensor renormalization group method and show how they can be applied for lattice field theory models involving relativistic fermions and Grassmann variables in arbitrary dimensions. We discuss recent progress for entanglement filtering, loop optimization, bond-weighting techniques and matrix product decompositions for Grassmann tensor networks. The new methods are tested with two-dimensional Wilson–Majorana fermions and multi-flavor Gross–Neveu models. We show that the methods can also be applied to the fermionic Hubbard model in 1+1 and 2+1 dimensions.
Holographic quantum-error correcting codes are models of bulk/boundary dualities such as the anti-de Sitter/conformal field theory (AdS/CFT) correspondence, where a higher-dimensional bulk geometry is associated with the code’s logical degrees of freedom. Previous discrete holographic codes based on tensor networks have reproduced the general code properties expected from continuum AdS/CFT, such as complementary recovery. However, the boundary states of such tensor networks typically do not exhibit the expected correlation functions of CFT boundary states. In this work, we show that a new class of exact holographic codes, extending the previously proposed hyperinvariant tensor networks into quantum codes, produce the correct boundary correlation functions. This approach yields a dictionary between logical states in the bulk and the critical renormalization group flow of boundary states. Furthermore, these codes exhibit a state-dependent breakdown of complementary recovery as expected from AdS/CFT under small quantum gravity corrections.
Here, we present a detailed investigation of a subdominant oscillating scalar field [“early dark energy” (EDE)] in the context of resolving the Hubble tension. Consistent with earlier work, but without relying on fluid approximations, we find that a scalar field frozen due to Hubble friction until log 10 ($z_c$) ~ 3.5, reaching $ρ_{\text{EDE}}(z_c)/ρ_{\text{tot}}$ ~ 10 % and diluting faster than matter afterwards, can bring cosmic microwave background (CMB), baryonic acoustic oscillations, supernovae luminosity distances, and the late-time estimate of the Hubble constant from the SH0ES Collaboration into agreement. A scalar field potential that scales as $V (\phi) \propto \phi ^{2n}$ with 2 ≲ $\textit{n}$ ≲ 3.4 around the minimum is preferred at the 68% confidence level, and the Planck polarization places additional constraints on the dynamics of perturbations in the scalar field. In particular, the data prefer a potential that flattens at large field displacements. A Markov-chain Monte Carlo analysis of mock data shows that the next-generation CMB observations (i.e., CMB-S4) can unambiguously detect the presence of the EDE at a very high significance. This projected sensitivity to the EDE dynamics is mainly driven by improved measurements of the E-mode polarization. We also explore new observational signatures of EDE scalar field dynamics: (i) We find that depending on the strength of the tensor-to-scalar ratio, the presence of the EDE might imply the existence of isocurvature perturbations in the CMB. (ii) We show that a strikingly rapid, scale-dependent growth of EDE field perturbations can result from parametric resonance driven by the anharmonic oscillating field for $\textit{n}$ ≈ 2. This instability and ensuing potentially nonlinear, spatially inhomogeneous, dynamics may provide unique signatures of this scenario.
We derive a factorization theorem that allows for resummation of small-x logarithms by exploiting Glauber operators in the soft collinear effective field theory. Our analysis is carried out for the hadronic tensor W μν in deep inelastic scattering, and leads to the definition of a new gauge invariant soft function S μν that describes quark and gluon emission in the central region. This soft function provides a new framework for extending resummed calculations for coefficient functions to higher logarithmic orders. Our factorization also defines impact factors by universal collinear functions that are process independent, for instance being identical in small-x DIS and Drell-Yan.
The boundary conditions corresponding to the design problem when the blades being simulated by the bound vorticity distribution are presented. The 3D flow is analyzed by the two steps S2 - S1 approach. In the first step, the number of blades is supposed to be infinite, the vortex distribution is transformed into an axisymmetric one, so that the flow field can be analyzed in a meridional plane. The thickness distribution of the blade producing the flow channel striction is taken into account by the modification of metric tensor in the continuity equation. Using the meridional stream function to define the flow field, the mass conservation is satisfied automatically. The governing equation is deduced from the relation between the azimuthal component of the vorticity and the meridional velocity. The value of the azimuthal component of the vorticity is provided by the hub to shroud equilibrium condition. This step leads to the determination of the axisymmetric stream sheets as well as the approximate camber surface of the blade. In the second step, the finite number of blades is taken into account, the inverse problem corresponding to the blade to blade flow confined in each stream sheet is analyzed. The momentum equation implies that the free vortex of the absolute velocity must be tangential to the stream sheet. The governing equation for the blade to blade flow stream function is deduced from this condition. At the beginning, the upper and the lower surfaces of the blades are created from the camber surface obtained from the first step with the assigned thickness distribution. The bound vorticity distribution and the penetrating flux conservation applied on the presumed blade surface constitute the boundary conditions of the inverse problem. The detection of this flux leads to the rectification of the geometry of the blades.
Motivated by the ground state structure of quantum models with all-to-all interactions such as mean-field quantum spin glass models and the Sachdev-Ye-Kitaev (SYK) model, we propose a tensor network architecture which can accomodate volume law entanglement and a large ground state degeneracy. We call this architecture the non-local renormalization ansatz (NoRA) because it can be viewed as a generalization of MERA, DMERA, and branching MERA networks with the constraints of spatial locality removed. We argue that the architecture is potentially expressive enough to capture the entanglement and complexity of the ground space of the SYK model, thus making it a suitable variational ansatz, but we leave a detailed study of SYK to future work. We further explore the architecture in the special case in which the tensors are random Clifford gates. Here the architecture can be viewed as the encoding map of a random stabilizer code. We introduce a family of codes inspired by the SYK model which can be chosen to have constant rate and linear distance at the cost of some high weight stabilizers. We also comment on potential similarities between this code family and the approximate code formed from the SYK ground space.
The complex interaction of forebody and wing vortices significantly impacts missile aerodynamics. The formation of these vortices involves smooth regions of the geometry or geometric discontinuities like leading edges, trailing edges, tips, and corners. Regions of supersonic flow and complex shock topologies interact with boundary layers and vortices. Smooth-body separation and 3D viscous effects strain current Reynolds-averaged Navier-Stokes (RANS) techniques. The quantification and control of discretization error is critical to obtaining reliable simulation results and often turbulence model assessments are made in the presence of unquantified (and potentially large) discretization errors. Two mesh adaptation schemes are applied to steady RANS simulations. Multiscale unstructured mesh adaptation is applied to control interpolation error estimates of the Mach field, which resolves boundary layers, vortices, and shocks. A dual-mesh approach with overset communication is applied between an expert-crafted near-body unstructured mesh and an adaptive off-body Cartesian mesh refined with Q-criterion scaled by the strain tensor magnitude. A generic missile configuration is examined in a supersonic flow field to show the interaction of mesh adaptation and turbulence model. Turbulence model modifications for rotational correction and a quadratic constitutive relationship show a strong influence on adaptive mesh refinement and predicted rolling moment.
We discuss the physics of sound propagation and charge diffusion in a plasma with non-vanishing charge density. Our analysis culminates the program initiated in to construct an open effective field theory of low-lying modes of the stress tensor and charge current in such plasmas. We model the plasma holographically as a Reissner-Nordström-AdS d+1 black hole, and study linearized fluctuations of longitudinally polarized scalar gravitons and photons in this background. We demonstrate that the perturbations can be decoupled and repackaged into the dynamics of two designer scalars, whose gravitational coupling is modulated by a non-trivial dilatonic factor. The holographic analysis allows us to isolate the phonon mode from the charge diffusion mode, and identify the combination of currents that corresponds to each of them. We use these results to obtain the real-time Gaussian effective action, which includes both the retarded response and the associated stochastic (Hawking) fluctuations, accurate to quartic order in gradients.
Abstract We study the mixed scalar-tensor two-point function in icosahedral inflation. Within the regime of validity of the effective field theory, this has to be perturbatively small; in particular, much smaller than the scalar spectrum. However, it can be much bigger that the tensor spectrum itself. We discuss observational implications for the CMB temperature-polarization spectrum.