Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Integrable field theories”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydrodynamic generators in relativistic kinetic theory

Here, we resum the nonequilibrium gradient corrections to a single-particle distribution function evolved by the Boltzmann equation in the relaxation time approximation (RTA). We first study a system undergoing Bjorken expansion and show that, for a constant relaxation time, the exact solution of the RTA Boltzmann equation at late times (i.e., after the decay of nonhydrodynamic modes) generates the Borel resummed Chapman-Enskog series. Extending this correspondence to systems without Bjorken symmetry, we construct a ( 3 + 1 )-dimensional hydrodynamic generator for RTA kinetic theory, which is an integral representation of the Chapman-Enskog series in the limit of vanishing nonhydrodynamic modes. Relaxing this limit, we discover at earlier times a set of nonhydrodynamic modes coupled to the RTA Chapman-Enskog expansion. Including the dynamics of these nonhydrodynamic modes is shown to control the emergence of hydrodynamics as an effective field theory description of nonequilibrium fluids, which works well even for far-off-equilibrium situations where the Knudsen number is initially large.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Post-Minkowskian effective field theory for conservative binary dynamics

We develop an Effective Field Theory (EFT) formalism to solve for the conservative dynamics of binary systems in gravity via Post-Minkowskian (PM) scattering data. Our framework combines a systematic EFT approach to compute the deflection angle in the PM expansion, together with the ‘Boundary-to-Bound’ (B2B) dictionary introduced. Due to the nature of scattering processes, a remarkable reduction of complexity occurs both in the number of Feynman diagrams and type of integrals, compared to a direct EFT computation of the potential in a PM scheme. We provide two illustrative examples. Firstly, we compute all the conservative gravitational observables for bound orbits to 2PM, which follow from only one topology beyond leading order. The results agree with those, obtained through the ‘impetus formula’ applied to the classical limit of the one loop amplitude. For the sake of comparison we reconstruct the conservative Hamiltonian to 2PM order, which is equivalent to the one derived from a matching calculation. Secondly, we compute the scattering angle due to tidal effects from the electric- and magnetic-type Love numbers at leading PM order. Using the B2B dictionary we then obtain the tidal contribution to the periastron advance. We also construct a Hamiltonian including tidal effects at leading PM order. Although relying on (relativistic) Feynman diagrams, the EFT formalism developed here does not involve taking the classical limit of a quantum amplitude, neither integrals with internal massive fields, nor additional matching calculations, nor spurious (‘super-classical’) infrared singularities. By construction, the EFT approach can be automatized to all PM orders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Acoustic scattering by smooth elastic cylinders insonified by directional transceivers: Monostatic theory and experiments

A theoretical model for predicting the acoustic field scattered by an elastic cylinder that is partially insonified by a directional transceiver is proposed in the form of a simple approximate one-dimensional integral. This model accounts for spherical spreading and directivity of the incident waves and extends the formulation used in a preceding article by including effects due to oblique insonification of a long cylinder assuming negligible end-contributions. The scattered field of an infinitely long cylinder for obliquely incident plane waves and point receivers is used to approximate the apparent volume flow of cylinders partially insonified by directional transceivers. The scattered pressure that is derived using the apparent volume flow, in contrast to the previous formulation, is capable of predicting axially propagating guided wave resonances; these natural modes are excited, in addition to circumferential ones, at off-normal incident angles. The model is compared with exact numerical simulations and with previously published as well as new laboratory data. The analysis illustrates the different realistic effects associated with scattering from elastic cylinders insonified by a directional transceiver both theoretically and experimentally.

42 ENGINEERING↗
NASA NTRS1972

Theoretical prediction of interference loading on aircraft stores. Part 1: Subsonic speeds

A method is developed for theoretically predicting the loading on pylon-mounted stores in subsonic compressible flow. Linear theory is used, without two-dimensional or slender body assumptions, to predict the flow field produced by the aircraft wing, nose, inlet, and pylons. The interference loading is integrated over the store length by considering the local crossflow, its axial and radial derivatives, and buoyancy. Store moment calculations under an F-4 aircraft at Mach 8 are compared to wind tunnel data. The method is computerized, and program user information is included.

Danfernandes, F.↗
NASA NTRS2015

Nonlinear Dynamic Modeling of a Supersonic Commercial Transport Turbo-Machinery Propulsion System for Aero-Propulso-Servo-Elasticity Research

This paper covers the development of an integrated nonlinear dynamic model for a variable cycle turbofan engine, supersonic inlet, and convergent-divergent nozzle that can be integrated with an aeroelastic vehicle model to create an overall Aero-Propulso-Servo-Elastic (APSE) modeling tool. The primary focus of this study is to provide a means to capture relevant thrust dynamics of a full supersonic propulsion system by using relatively simple quasi-one dimensional computational fluid dynamics (CFD) methods that will allow for accurate control algorithm development and capture the key aspects of the thrust to feed into an APSE model. Previously, propulsion system component models have been developed and are used for this study of the fully integrated propulsion system. An overview of the methodology is presented for the modeling of each propulsion component, with a focus on its associated coupling for the overall model. To conduct APSE studies the described dynamic propulsion system model is integrated into a high fidelity CFD model of the full vehicle capable of conducting aero-elastic studies. Dynamic thrust analysis for the quasi-one dimensional dynamic propulsion system model is presented along with an initial three dimensional flow field model of the engine integrated into a supersonic commercial transport.

Control Theory↗
NASA NTRS2021

Boundary-layer Instability Measurements on a Cone at Freestream Mach 3.5

An experimental study was conducted in the NASA (National Aeronautics and Space Administration) Langley Supersonic Low-Disturbance Tunnel to investigate naturally-occurring instabilities in a supersonic boundary layer on a 7◦ half-angle cone at nominal freestream conditions: Mach 3.5, total temperature of 299.8K, and unit Reynolds numbers (millions per m) of 9.89, 13.85, 21.77, and 25.73. Instability measurements were acquired under noisy-flow and quiet-flow conditions. Pitot-pressure and calibrated hot-wire measurements were obtained using a model-integrated traverse system to document the model flow field. In noisy-flow conditions, growth rates and mode shapes achieved good agreement between the measured results and linear stability theory (LST). The corresponding N factor at transition from LST is N ≈ 3.9. Under quiet-flow conditions, the most unstable first-mode instabilities as predicted by LST were measured, but this mode was not the dominant instability measured in the boundary layer. Instead, the dominant instabilities were less-amplified, low-frequency disturbances predicted by LST, and grew according to linear theory. These low-frequency unstable disturbances were initiated by freestream acoustic disturbances through a receptivity process believed to occur near the branch I location of the cone. Under quiet-flow conditions, the boundary layer remained laminar up to the last measurement station for the largest unit Reynolds number, implying a transition N factor of N > 8.5.

Boundary Layer↗
DOE OSTI2022

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Radiative corrections to superallowed beta decays at $\mathcal{O}\left({\alpha}^2Z\right)$

We compute $\mathcal{O}\left({\alpha}^2Z\right)$ radiative corrections to superallowed β decays with a heavy-particle effective field theory that systematically describes the interactions of low-energy ultrasoft photons with nuclei. We calculate two-loop virtual and one-loop real-virtual amplitudes by reducing the Feynman integrals to a set of master integrals, which we solve analytically using a variety of techniques. These techniques can be applied to other phenomenologically interesting observables. The ultrasoft corrections can then be combined with contributions arising from the exchange of potential photons to obtain the complete $\mathcal{O}\left({\alpha}^2Z\right)$ correction to the decay rate, with resummation of large logarithms of the electron energy times the nuclear radius. We find that $\mathcal{O}\left({\alpha}^2Z\right)$ ultrasoft loops induce a relative correction to the decay rate that ranges from 0.7 ∙ 10 −3 in the decay of 10 C to 3.6 ∙ 10 −3 in the decay of 54 Co, and will thus impact the extraction of V ud at the permille level. We show that the inclusion of these corrections reduces the residual renormalization scale dependence of the decay rate to a negligible level, making missing ultrasoft perturbative corrections a subdominant source of theoretical uncertainty.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1999

Crack Turning and Arrest Mechanisms for Integral Structure

In the course of several years of research efforts to predict crack turning and flapping in aircraft fuselage structures and other problems related to crack turning, the 2nd order maximum tangential stress theory has been identified as the theory most capable of predicting the observed test results. This theory requires knowledge of a material specific characteristic length, and also a computation of the stress intensity factors and the T-stress, or second order term in the asymptotic stress field in the vicinity of the crack tip. A characteristic length, r(sub c), is proposed for ductile materials pertaining to the onset of plastic instability, as opposed to the void spacing theories espoused by previous investigators. For the plane stress case, an approximate estimate of r(sub c), is obtained from the asymptotic field for strain hardening materials given by Hutchinson, Rice and Rosengren (HRR). A previous study using of high order finite element methods to calculate T-stresses by contour integrals resulted in extremely high accuracy values obtained for selected test specimen geometries, and a theoretical error estimation parameter was defined. In the present study, it is shown that a large portion of the error in finite element computations of both K and T are systematic, and can be corrected after the initial solution if the finite element implementation utilizes a similar crack tip discretization scheme for all problems. This scheme is applied for two-dimensional problems to a both a p-version finite element code, showing that sufficiently accurate values of both K(sub I) and T can be obtained with fairly low order elements if correction is used. T-stress correction coefficients are also developed for the singular crack tip rosette utilized in the adaptive mesh finite element code FRANC2D, and shown to reduce the error in the computed T-stress significantly. Stress intensity factor correction was not attempted for FRANC2D because it employs a highly accurate quarter-point scheme to obtain stress intensity factors.

Pettit, Richard↗
DOE OSTI2021

Interstellar Extinction, Polarization, and Grain Alignment in the Sh 2-185 (IC 59 and IC 63) Region

Optical and infrared continuum polarization from the interstellar medium is driven by radiative processes aligning the grains with the magnetic field. While a quantitative, predictive theory of radiative alignment torques (RATs) exists and has been extensively tested, several parameters of the theory remain to be fully constrained. In a recent paper, Medan & Andersson showed that the polarization efficiency (and therefore grain alignment efficiency) at different locations in the wall of the Local Bubble (LB) could be modeled as proportional to the integrated light intensity from the surrounding stars and OB associations. Here we probe that relationship at high radiation field intensities by studying the extinction and polarization in the two reflection nebulae IC 59 and IC 63 in the Sh 2-185 H ii region, illuminated by the B0 IV star γ Cassiopeia. We combine archival visual polarimetry with new seven-band photometry in the Vilnius system, to derive the polarization efficiency from the material. We find that the same linear relationship seen in the LB wall also applies to the Sh 2-185 region, strengthening the conclusion from the earlier study.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A search for heavy Higgs bosons decaying into vector bosons in same-sign two-lepton final states in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for heavy Higgs bosons produced in association with a vector boson and decaying into a pair of vector bosons is performed in final states with two leptons (electrons or muons) of the same electric charge, missing transverse momentum and jets. A data sample of proton–proton collisions at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider between 2015 and 2018 is used. The data correspond to a total integrated luminosity of 139 fb –1 . The observed data are in agreement with Standard Model background expectations. The results are interpreted using higher-dimensional operators in an effective field theory. Upper limits on the production cross-section are calculated at 95% confidence level as a function of the heavy Higgs boson’s mass and coupling strengths to vector bosons. Limits are set in the Higgs boson mass range from 300 to 1500 GeV, and depend on the assumed couplings. The highest excluded mass for a heavy Higgs boson with the coupling combinations explored is 900 GeV. Limits on coupling strengths are also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Axions in neutron star mergers

Supernovae and cooling neutron stars have long been used to constrain the properties of axions, such as their mass and interactions with nucleons and other Standard Model particles. We investigate the prospects of using neutron star mergers as a similar location where axions can be probed in the future. We examine the impact axions would have on mergers, considering both the possibility that they free-stream through the dense nuclear matter and the case where they are trapped. We calculate the mean free path of axions in merger conditions, and find that they would free-stream through the merger in all thermodynamic conditions. In contrast to previous calculations, we integrate over the entire phase space while using a relativistic treatment of the nucleons, assuming the matrix element is momentum-independent. In particular, we use a relativistic mean field theory to describe the nucleons, taking into account the precipitous decrease in the effective mass of the nucleons as density increases above nuclear saturation density. We find that within current constraints on the axion-neutron coupling, axions could cool nuclear matter on timescales relevant to neutron star mergers. Our results may be regarded as first steps aimed at understanding how axions affect merger simulations and potentially interface with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Noncollinear ground states of solids with a source-free exchange correlation functional

In this paper, we expand upon the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and coworkers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+U+J functional. In the process, we compare the noncollinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground-state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin-polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analyses on probability currents, net torque, and constant terms draw additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and pave the way for future extensions and corrections to the SF corrected XC functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

Physics of Rotation Powered Pulsars and Their Nebulae

This letter is my progress report for the Astrophysics Theory grant. The first goal of the research supported by this grant is to produce a time dependent theory of the unsteady relativistic collisionless shock wave terminating the relativistic wind from a pulsar, and compare the predicted surface brightness fluctuations to Hubble Space Telescope observations of the wisps in the Crab Nebula. The second goal is to model the production of electron-positron pairs over the polar caps of rotation powered pulsars, and use the results to predict the heating of the surface due to particle trapping and bombardment of the atmosphere at the base of the polar field lines. We have succeeded in creating a one dimensional hybrid code, in which the electron-positron pairs incident on the shock structure are modeled as a relativistic, adiabatic, ideal MHD fluid, while the heavy ions are treated as particles using a particle-in-cell algorithm. The electromagnetic fields are evaluated from the currents and charge densities in the pairs and the ions, while the particles and the fluid accelerate in response to the computed self-consistent electromagnetic fields. The results are promising, in that the underlying ion cyclotron instability generates finite amplitude, propagating magnetosonic waves in the pairs, whose wavelengths and frequencies, when translated into physical units, are comparable to the observed running waves of brightness observed by HST near the Crab pulsar. The code is undergoing a number of tests, to assure us that this preliminary correspondence is not an artifact. In the coming year, the observational appearance of the models will be computed and compared to the HST observations of the Crab now in hand, and used to predict the HST results which will be obtained the year after next. WE also developed a one dimensional cascade theory for pair creation over pulsars' polar caps. A linear integral equation describing the synchrotron cascade has been derived and solved by iterative techniques, in the case when a high energy electron moving parallel to a star centered dipole magnetic field initiates the cascade through curvature gamma ray emission.

Arons, Jonathan↗
DOE OSTI2022

Giant gravitons, Harish-Chandra integrals, and BPS states in symplectic and orthogonal $\mathcal{N}$ = 4 SYM

We find generating functions for half BPS correlators in $\mathcal{N}$ = 4 SYM theories with gauge groups Sp(2N), SO(2N + 1), and SO(2N) by computing the norms of a class of BPS coherent states. These coherent states are built from operators involving Harish-Chandra integrals. Such operators have an interpretation as localized giant gravitons in the bulk of anti-de-Sitter space. This extends the analysis of [1] to Sp(2N), SO(2N + 1), and SO(2N) gauge theories. We show that we may use ordinary Schur functions as a basis for the sector of states with no cross-caps in these theories. This is consistent with the construction of these theories as orientifold projections of an SU(2N) theory. We make note of some relations between the symmetric functions that appear in the expansion of these coherent states and symplectic Schur functions. We also comment on some connections to Schubert calculus and Gromov-Witten invariants, which suggest that the Harish-Chandra integral may be extended to such problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1974

Subsonic potential aerodynamics for complex configurations - A general theory

A general theory of subsonic potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the velocity potential is obtained for both supersonic and subsonic flow. Under the small perturbation assumption, the potential at any point in the field depends only upon the values of the potential and its normal derivative on the surface of the body. On the surface of the body, this representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The theory is applied to finite-thickness wings in subsonic steady and oscillatory flows.

Morino, L.↗
DOE OSTI2025

High entropy powering green energy: hydrogen, batteries, electronics, and catalysis

A reformation in energy is underway to replace fossil fuels with renewable sources, driven by the development of new, robust, and multi-functional materials. High-entropy materials (HEMs) have emerged as promising candidates for various green energy applications, having unusual chemistries that give rise to remarkable functionalities. This review examines recent innovations in HEMs, focusing on hydrogen generation/storage, fuel cells, batteries, semiconductors/electronics, and catalysis—where HEMs have demonstrated the ability to outperform state-of-the-art materials. We present new master plots that illustrate the superior performance of HEMs compared to conventional systems for hydrogen generation/storage and heat-to-electricity conversion. We highlight the role of computational methods, such as density functional theory and machine learning, in accelerating the discovery and optimization of HEMs. The review also presents current challenges and proposes future directions for the field. We emphasize the need for continued integration of modeling, data, and experiments to investigate and leverage the underlying mechanisms of the HEMs that are powering progress in sustainable energy.

batteries↗
DOE OSTI2020

Insights into lithium ion deposition on lithium metal surfaces

Lithium metal is among the most promising anodes for the next generation of batteries due to its high theoretical energy density and high capacity. Challenges such as extreme reactivity and lithium dendrite formation have kept lithium metal anodes away from practical applications. However, the underlying mechanisms of Li ion deposition from the electrolyte solution onto the anode surface are still poorly understood due to their inherent complexity. In this work, density functional theory calculations and thermodynamic integration via constrained molecular dynamics simulations are conducted to study the electron and ion transfer between lithium metal slab and the electrolyte in absence of an external field. Here, we explore the effect of the solvent chemistry and structure, distance of the solvated complex from the surface, anion–cation separation, and concentration of Li-salts on the deposition of lithium ions from the electrolyte phase onto the surface. Ethylene carbonate (EC), 1,2-dimethoxyethane (DME), 1,3-dioxolane (DOL), and mixtures of them are used as solvents. These species compete with the salt anion and the Li cation for electron transfer from the surface. It is found that the structure and properties of the solvation shell around the lithium cation has a great influence on the ability of the cation to diffuse as well as on its surrounding electron environment. DME molecules allow easier motion of the lithium ion compared with EC and DOL molecules. The slow growth approach allows the study of energy barriers for the ion diffusion and desolvation during the deposition pathway. This method helps elucidating the underlying mechanisms on lithium-ion deposition and provides a better understanding of the early stages of Li nucleation.

25 ENERGY STORAGE↗