Search NASA⌕ Search

SEARCH · Search NASA

Results for “THEORIES - FLOW LINEAR”

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

Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fermionic shift symmetries in (anti) de Sitter space

We study extended shift symmetries that arise for fermionic fields on anti-de Sitter (AdS) space and de Sitter (dS) space for particular values of the mass relative to the curvature scale. We classify these symmetries for general mixed-symmetry fermionic fields in arbitrary dimension and describe how fields with these symmetries arise as the decoupled longitudinal modes of massive fermions as they approach partially massless points. For the particular case of AdS 4 , we look for non-trivial Lie superalgebras that can underly interacting theories that involve these fields. We study from this perspective the minimal such theory, the Akulov-Volkov theory on AdS 4 , which is a non-linear theory of a spin-1/2 Goldstino field that describes the spontaneous breaking of $\mathcal{N}$ = 1 supersymmetry on AdS 4 down to the isometries of AdS 4 . We show how to write the nonlinear supersymmetry transformation for this theory using the fermionic ambient space formalism. We also study the Lie superalgebras of candidate multi-field examples and rule out the existence of a supersymmetric special galileon on AdS 4 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Parallel Diffusion Coefficient of Energetic Charged Particles in the Inner Heliosphere from the Turbulent Magnetic Fields Measured by Parker Solar Probe

Diffusion coefficients of energetic charged particles in turbulent magnetic fields are a fundamental aspect of diffusive transport theory but remain incompletely understood. In this work, we use quasi-linear theory to evaluate the spatial variation of the parallel diffusion coefficient κ ∥ from the measured magnetic turbulence power spectra in the inner heliosphere. We consider the magnetic field and plasma velocity measurements from Parker Solar Probe made during Orbits 5–13. The parallel diffusion coefficient is calculated as a function of radial distance from 0.062 to 0.8 au, and the particle energy from 100 keV to 1 GeV. We find that κ ∥ increases exponentially with both heliocentric distance and energy of particles. The fluctuations in κ ∥ are related to the episodes of large-scale magnetic structures in the solar wind. By fitting the results, we also provide an empirical formula of κ ∥ = (5.16 ± 1.22) × 10 18 r 1.17 ± 0.08 E 0.71 ± 0.02 (cm 2 s -1 ) in the inner heliosphere, which can be used as a reference in studying the transport and acceleration of solar energetic particles as well as the modulation of cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two-dimensional simulations of internal gravity waves in the radiation zones of intermediate-mass stars

Intermediate-mass main-sequence stars have large radiative envelopes overlying convective cores. This configuration allows internal gravity waves (IGWs) generated at the convective–radiative interface to propagate towards the stellar surface. The signatures of these waves can be observed in the photometric and spectroscopic data from stars. Here, we have studied the propagation of these IGWs using two-dimensional (2D) fully non-linear hydrodynamical simulations with realistic stellar reference states from the 1D stellar evolution code, Modules for Stellar Astrophysics (MESA). When a single wave is forced, we observe wave self-interaction. When two waves are forced, we observe non-linear interaction (i.e. triadic interaction) between these waves forming waves at different wavelengths and frequencies. When a spectrum of waves similar to that found in numerical simulations is forced, we find that the surface IGW frequency slope is consistent with recent observations. This power law is similar to that predicted by linear theory for the wave propagation, with small deviations that can be an effect of non-linearities. When the same generation spectrum is applied to 3 M ⊙ models at different stellar rotation and ages, the surface IGW spectrum slope is very similar to the generation spectrum slope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cosmological constraints from the cross-correlation of DESI Luminous Red Galaxies with CMB lensing from Planck PR4 and ACT DR6

Here, we infer the growth of large scale structure over the redshift range 0.4 ≲ z ≲ 1 from the cross-correlation of spectroscopically calibrated Luminous Red Galaxies (LRGs) selected from the Dark Energy Spectroscopic Instrument (DESI) legacy imaging survey with CMB lensing maps reconstructed from the latest Planck and ACT data. We adopt a hybrid effective field theory (HEFT) model that robustly regulates the cosmological information obtainable from smaller scales, such that our cosmological constraints are reliably derived from the (predominantly) linear regime. We perform an extensive set of bandpower- and parameter-level systematics checks to ensure the robustness of our results and to characterize the uniformity of the LRG sample. We demonstrate that our results are stable to a wide range of modeling assumptions, finding excellent agreement with a linear theory analysis performed on a restricted range of scales. From a tomographic analysis of the four LRG photometric redshift bins we find that the rate of structure growth is consistent with ΛCDM with an overall amplitude that is ≃ 5-7% lower than predicted by primary CMB measurements with modest (∼ 2σ) statistical significance. From the combined analysis of all four bins and their cross-correlations with Planck we obtain S 8 = 0.765 ± 0.023, which is less discrepant with primary CMB measurements than previous DESI LRG cross Planck CMB lensing results. From the cross-correlation with ACT we obtain S 8 = 0.790 +0.024 -0.027 , while when jointly analyzing Planck and ACT we find S 8 = 0.775 +0.019 -0.022 from our data alone and σ 8 = 0.772 +0.020 -0.023 with the addition of BAO data. These constraints are consistent with the latest Planck primary CMB analyses at the ≃ 1.6-2.2σ level, and are in excellent agreement with galaxy lensing surveys.

cosmological parameters from LSS↗
DOE OSTI2024

Impossible symmetries and conformal gravity

We explore the physics of relativistic gapless phases defined by a mixed anomaly between two generalized conserved currents. The gapless modes can be understood as Goldstone modes arising from the nonlinear realization of (generically higher-form) symmetries arising from these currents. In some cases, the anomaly cannot be reproduced by any local and unitary theory, indicating that the corresponding symmetries are impossible, in the sense that they cannot appear in a Lorentzian physical system. We give a general construction and illustrate it with several examples. Most notably, we study conformal gravity from this perspective, describing the higher-form symmetries of the linear theory and showing how it can be understood in terms of anomalies. Along the way we clarify some aspects of electric-magnetic duality in linear conformal gravity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bi-Kappa Proton Mirror and Cyclotron Instabilities in the Solar Wind

Abstract The charged particles in the solar wind are often observed to possess a nonthermal tail in the velocity distribution function, a feature that can be fitted with the Kappa model. The anisotropic, or bi-Kappa, model of protons, electrons, and other charged particles is thus adopted in the literature for interpreting the data as well as in the context of the analysis of wave–particle interactions. The present paper develops an approximate but efficient theory of the mirror and cyclotron instabilities excited by the bi-Kappa protons in the solar wind. A velocity moment-based quasi-linear theory of these instabilities is also formulated in order to investigate the saturation behavior. Applications of the formalism are made for instabilities close to the marginally unstable state, which is typical of the solar wind near 1 au.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Improving initialization and evolution accuracy of cosmological neutrino simulations

Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations — the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. Here, a recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart Script $\mathcal{O}$(1%) numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the 15% level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of m ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Consistent linearization of micromorphic continuum theories

We outline the procedure of consistent linearization and apply it to the micromorphic, microstretch, and micropolar theories of continua. This yields tractable linear theories for nonlinearly elastic microstructured materials undergoing finite deformations. The results may be readily utilized in computational mechanics, stability and bifurcation analyses, and small-deformation problems in the context of these types of continua. Our results generalize those existing in current literature and facilitate their recovery upon incorporating appropriate kinematic and constitutive assumptions.

McAvoy, Ryan C.↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Design of Cylindrical Implosion Experiments to Demonstrate Scale-Invariant Rayleigh-Taylor Instability Growth

Radiation-hydrodynamics simulations are used to design laser-driven cylindrical implosion experiments to directly measure hydrodynamic instability growth in convergent geometry. Designs for two different size targets, varying in radial dimension by a factor of three, are presented. A set of beam pointings and powers are identified for each scale design that result in a nearly axially uniform implosion of an embedded marker layer. The implosion trajectories are shown to be scale-invariant between designs, with nearly identical scaled acceleration profiles. Finally, linear theory and radiation-hydrodynamics simulations predict that Rayleigh-Taylor instability growth of an azimuthal perturbation, machined on the inner surface of the embedded marker, is also scale-invariant between designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Forward and backward stimulated Raman scattering in multi-speckled beams: Density dependence and effects on cross-beam energy transfer

Particle-in-cell simulations of cross-beam energy transfer (CBET) between multi-speckled beams reveal a potentially prominent energy loss from forward stimulated Raman scattering (FSRS) in the amplified seed beam [Yin et al., Phys. Plasmas 26, 082708 (2019)]. In this work, an intensity scan is used to examine this energy loss in the CBET-amplified seed beam at intensities of a few times 1014 W/cm2. Electron tracers demonstrate both the trapping in the electron plasma waves that is sustained over multiple speckles and the subsequent acceleration to velocities well beyond the phase velocity. Motivated by expanding our understanding of this energy loss mechanism, we isolate the FSRS density dependence by using a 2×1014 W/cm2 pulse to scan densities between 0.06 ncr and 0.12 ncr; a spectral analysis characterizes the energy losses from scattering with the changing densities. At saturation, both greater beam attenuation is measured at the laser frequency and larger scattering FSRS angles are observed as density increases, along with a rapid transition to primarily backscatter SRS between 0.10 ncr and 0.12 ncr, earlier than predicted by linear theory. Density gradients are shown to effectively mitigate oblique FSRS and beam energy losses in the two lower density cases, whereas 0° forward scattering is primarily suppressed in the higher density cases. Finally returning to the two-beam CBET system of interest, the interplay and mitigation of the secondary instabilities at larger domain and beam spatial scales are discussed, here also examining backward stimulated Brillouin scattering.

Stark, D. J. (ORCID:0000000203642733)↗
DOE OSTI2023

Electronic density response of warm dense matter

Matter at extreme temperatures and pressures—commonly known as warm dense matter (WDM)—is ubiquitous throughout our Universe and occurs in astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are probed in x-ray Thomson scattering (XRTS) experiments and are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear response theory and nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques, including path-integral Monte Carlo simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrary complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations and models for the Richtmyer–Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P ∝ km, where k is the wavenumber. The hydrodynamic condition is the same as in the θ-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ∼1.86 shock and gamma-law = 5/3 fluids with Atwood number A = −0.5. The bubble amplitude hB is found to grow in two stages. Initially, hB undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, hB grows as a power law ∼tθ when k|hB| > O(1), where k is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

PIAFS: A 2D nonlinear hydrodynamics code to model gaseous optics

The survivability of final optics is expected to be a major challenge for all future inertial fusion energy concepts. Due to their higher damage threshold, gaseous optics have been identified as a promising solution to this problem. Gaseous optics can be created through the photoabsorption of spatially modulated UV light, which induces various chemical processes that heat the gas. This heating leads to a pressure perturbation, which in turn launches a density perturbation that can imprint a refractive index modulation such as a grating. In this article, we introduce a parallel C/C++ code to simulate gaseous optics. PIAFS2D is a high-order conservative finite-difference code to solve the compressible Navier–Stokes equations along with the photochemical heating sources on Cartesian grids. The simulations are validated by the linear theory derived in a previous paper [Michel et al., Phys. Rev. Appl. 22, 024014 (2024)]. For larger perturbations, the behavior of the system—particularly the evolution of the generated acoustic wave—demonstrates strong nonlinearity. PIAFS2D allows the study of nonlinear behaviors and can be used for the design of high-efficiency gaseous optics elements in realistic experimental conditions.

Oudin, A. [Lawrence Livermore National Laboratory ↗