Search NASA⌕ Search

SEARCH · Search NASA

Results for “asymptotic homogenization”

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.

34 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Insights on the peak in the speed of sound of ultradense matter

In this work we investigate the minimal physical requirements needed for generating a speed of sound that surpasses its asymptotic conformal limit. It is shown that a peak in the speed of sound of homogeneous matter naturally emerges in the transition from a phase with broken chiral symmetry to one with a gapped Fermi surface. We argue that this could be relevant for understanding the peak in the speed of sound displayed by some of the current models for cold ultradense matter. A minimal model implementation of this mechanism is presented, based on the spontaneous breakdown of an approximate particle-antiparticle symmetry, and its thermodynamic properties are determined.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Inflation alternative via the gravitational field of a singularity

Abstract We explore the scenario that the observable Universe emerged from the vicinity of a negative mass ring singularity, and all content of the Universe travels at the same group velocity close to the speed of light on a geodesic trajectory along the axis of rotation of the singularity. In appropriate coordinate parametrization and evaluated on the trajectory, we find that the metric tensor in the vicinity of the trajectory exhibits a conformal scale factor a ( η ) with contraction and subsequent expansion properties that solve the horizon problem. We then introduce a static flow of gravitating radiation along the trajectory (perturbatively with respect to the mass scale of the singularity) to model a homogeneous radiation dominated Universe. Solving the Einstein field equations with a physically motivated ansatz of metric perturbation then reveals that the effective conformal scale factor indeed grows asymptotically with the same power law as expected in a conventional radiation dominated Universe.

Astronomy & Astrophysics↗
DOE OSTI2026

Is there an exact magnetic moment for charged particle motion in a time-dependent, homogeneous magnetic field?

The non-perturbative guiding-centre model provides an exact alternative to full-orbit simulations of charged particle dynamics in situations where traditional guiding-centre theory may fail. We demonstrate that the charged particle motion in a homogeneous, time-varying magnetic field is a solvable example of the non-perturbative guiding-centre model. This entails showing that the exact magnetic moment of Qin and Davidson can be constructed to be asymptotic to the adiabatic invariant series of Kruskal. In contrast to the perturbative invariant, the exact invariant contains information about parametric resonances. These resonances destroy the conservation of the usual magnetic moment over very long times. This refutes some previous claims about the all-time invariance of the magnetic moment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
DOE OSTI2023

Magnetic potential based formulation for linear and non-linear 3D RF sheath simulation

Abstract This paper reports a new numerical scheme to simulate the radio-frequency (RF) induced RF sheath, which is suitable for a large 3D simulation. In the RF sheath boundary model, the tangential component of the electric field ( E t ) is given by the gradient of a scalar electric field potential. We introduce two additional scalar potentials for the tangential components of the magnetic field, which effectively impose the normal electric displacement ( D n ) on the plasma sheath boundary condition via in-homogeneous Neumann boundary condition and constrain the tangential electric field on the surface as curl-free ( ∇ × E t = 0 ). In our approach, the non-linear sheath impedance is formulated as a natural extension of the large thickness (or asymptotic) sheath limit ( D n = 0 ), allowing for handling both asymptotic and non-linear regimes seamlessly. The new scheme is implemented using the Petra-M finite element method analysis framework and is verified with simulations in the literature. The significance of non-linearity is discussed in various plasma conditions. An application of this scheme to asymptotic RF sheath simulation on the WEST ICRF antenna side limiters is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effect of boundary conditions on fluctuations of the Bose condensate of interacting atoms

For a gas of trapped interacting atoms forming a Bose condensate, we demonstrate a possibility of a significant influence of the boundary conditions on the statistics of the number of particles in the fundamental mode of the system. The analysis is carried out in the Bogoliubov – Popov low-temperature approximation for model homogeneous cubic traps with periodic or combined periodic and zero boundary conditions. It is shown that the effect of the boundary conditions does not weaken even in the region of parameters corresponding to the relatively strong interaction in the Thomas – Fermi asymptotic case, and does not disappear when proceeding to the thermodynamic limit. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A Stochastic Transport Model for the Cumulative Number of Fissions and Deposited Fission Energy

The stochastic theory of neutron transport is extended to describe the cumulative distribution of fission numbers and deposited fission energy in a multiplying assembly. Solutions for the probability distributions are obtained using analytical approximations and Monte Carlo simulation in lumped geometry and in symmetric homogeneous and heterogeneous spheres. The results show the development of a power-law tail in the steady state fission number and deposited energy distributions when the medium is critical, independent of the fission neutron multiplicity distribution and domain heterogeneity. In contrast, the asymptotic decay is faster than exponential in subcritical media due to rapid chain extinction and in supercritical media due to the increasing probability of chain divergence. Here, a formal asymptotic analysis of the problem in lumped geometry with an arbitrary fission neutron multiplicity confirms the existence of power-law tails at critical.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗
DOE OSTI2022

Transient phenomena in the axion assisted Schwinger effect

Particle production induced by a time-dependent background is well understood as the projection of the time-evolved initial state onto a set of final states. While the asymptotic initial and final states are well defined in the usual way, the definition of particles and antiparticles at intermediate times in the presence of external fields is ambiguous. These external fields moreover induce divergences which require regularization. In this paper we clarify some subtleties in the computation of transient effects in physical quantities for fermions in a homogeneous axion background, including Schwinger production in background electromagnetic fields. The presence of the axion requires particular care as well as knowledge of the UV theory when regulating the theory and computing the vacuum contribution to the fermion energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Gravity loop integrands from the ultraviolet

We demonstrate that loop integrands of (super-)gravity scattering amplitudes possess surprising properties in the ultraviolet (UV) region. In particular, we study the scaling of multi-particle unitarity cuts for asymptotically large momenta and expose an improved UV behavior of four-dimensional cuts through seven loops as compared to standard expectations. For N=8 supergravity, we show that the improved large momentum scaling combined with the behavior of the integrand under BCFW deformations of external kinematics uniquely fixes the loop integrands in a number of non-trivial cases. In the integrand construction, all scaling conditions are homogeneous. Therefore, the only required information about the amplitude is its vanishing at particular points in momentum space. This homogeneous construction gives indirect evidence for a new geometric picture for graviton amplitudes similar to the one found for planar N=4 super Yang-Mills theory. We also show how the behavior at infinity is related to the scaling of tree-level amplitudes under certain multi-line chiral shifts which can be used to construct new recursion relations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effects of Atwood and Reynolds numbers on the evolution of buoyancy-driven homogeneous variable-density turbulence

The evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at Atwood numbers up to 0.75 and large Reynolds numbers is studied by using high-resolution direct numerical simulations. To help understand the highly non-equilibrium nature of buoyancy-driven HVDT, the flow evolution is divided into four different regimes based on the behaviour of turbulent kinetic energy derivatives. The results show that each regime has a unique type of dependence on both Atwood and Reynolds numbers. It is found that the local statistics of the flow based on the flow composition are more sensitive to Atwood and Reynolds numbers compared to those based on the entire flow. It is also observed that, at higher Atwood numbers, different flow features reach their asymptotic Reynolds-number behaviour at different times. The energy spectrum defined based on the Favre fluctuations momentum has less large-scale contamination from viscous effects for variable-density flows with constant properties, compared to other forms used previously. The evolution of the energy spectrum highlights distinct dynamical features of the four flow regimes. Thus, the slope of the energy spectrum at intermediate to large scales evolves from -7/3 to -1 , as a function of the production-to-dissipation ratio. The classical Kolmogorov spectrum emerges at intermediate to high scales at the highest Reynolds numbers examined, after the turbulence starts to decay. Finally, the similarities and differences between buoyancy-driven HVDT and the more conventional stationary turbulence are discussed and new strategies and tools for analysis are proposed.

42 ENGINEERING↗
DOE OSTI2020

Unitary quantum lattice simulations for Maxwell equations in vacuum and in dielectric media

Utilizing the similarity between the spinor representation of the Dirac and the Maxwell equations that has been recognized since the early days of relativistic quantum mechanics, a quantum lattice algorithm (QLA) representation of unitary collision-stream operators of Maxwell's equations is derived for both homogeneous and inhomogeneous media. A second-order accurate 4-spinor scheme is developed and tested successfully for two-dimensional (2-D) propagation of a Gaussian pulse in a uniform medium whereas for normal (1-D) incidence of an electromagnetic Gaussian wave packet onto a dielectric interface requires 8-component spinors because of the coupling between the two electromagnetic polarizations. In particular, the well-known phase change, field amplitudes and profile widths are recovered by the QLA asymptotic profiles without the imposition of electromagnetic boundary conditions at the interface. The QLA simulations yield the time-dependent electromagnetic fields as the wave packet enters and straddles the dielectric boundary. QLA involves unitary interleaved non-commuting collision and streaming operators that can be coded onto a quantum computer: the non-commutation being the very reason why one perturbatively recovers the Maxwell equations.

Physics↗
DOE OSTI2022

Spatially resolved properties of early-type group-dominant galaxies with MUSE: gas content, ionization mechanisms, and metallicity gradients

ABSTRACT With the goal of a thorough investigation of the ionized gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). These data allowed us to study the spatially resolved warm gas properties, including the morphology of the ionized gas, EW(H α), and kinematics as well as the gas-phase metallicity (12 + log(O/H)) of these systems. In order to distinguish between different ionization mechanisms, we used the emission-line ratios [O iii]/H β and [N ii]/H α in the BPT diagrams and EW(H α). We find that the ionization sources in our sample have variable impacts at different radii; central regions are more influenced by low-luminosity active galactic nuclei, while extended regions of low-ionization nuclear emission-line region-like emission are ionized by other mechanisms with post-asymptotic giant branch stars photoionization likely contributing significantly. We classified our sample into three H α +[N ii] emission morphology types. We calculate the gas-phase metallicity assuming several methods and ionization sources. In general, 12 + log(O/H) decreases with radius from the centre for all galaxies, independently of nebular morphology type, indicating a metallicity gradient in the abundance profiles. Interestingly, the more extended filamentary structures and all extranuclear star-forming regions present shallow metallicity gradients. Within the uncertainties these extended structures can be considered chemically homogeneous. We suggest that group-dominant galaxies in our sample likely acquired their cold gas in the past as a consequence of one or more mechanisms, e.g. gas-clouds or satellite mergers/accretion and/or cooling flows that contribute to the growth of the ionized gas structures.

Lagos, P. (ORCID:0000000223218657)↗
DOE OSTI2024

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗
DOE OSTI2025

Neutron-capture Element Abundances of 491 Stars in Milky Way Dwarf Satellite Galaxies from Medium-resolution Spectra

The chemical compositions of evolved stars in Local Group dwarf spheroidal galaxies (dSphs) provide insight into the galaxy’s past star formation and nucleosynthesis. Neutron-capture element abundances are especially interesting. In particular, s-process elements can provide a third chemical clock for resolving star formation histories in addition to core-collapse and Type Ia supernovae. Likewise, the primary sites of the r-process are still areas of extensive research. Until now, the number of stars with neutron-capture element abundances in dSphs has been limited by the need for stars bright enough for high-resolution spectroscopy. We present abundance measurements of the neutron-capture elements Sr, Y, Ba, and Eu with errors <0.4 dex—as well as new measurements of Mg—in 491 stars in Sculptor, Fornax, Draco, Sextans, and Ursa Minor. The large number of stars in our sample is possible because we used medium-resolution spectra from the DEIMOS spectrograph, assembling the largest homogeneous set of neutron-capture abundances in dSphs to date. By utilizing the abundances of both s- and r-process elements, we find evidence of an s-process contribution at early times in Sculptor from our measurements of [Ba/Fe]. This is a potential signature of s-process nucleosynthesis in fast-rotating massive stars. By comparing our measurements of [Eu/Fe] with [Mg/Fe], we show the need for an r-process source that has a short delay time to enrich stars in the dSphs. Thus, neutron star mergers are likely not the sole source of r-process material in dSphs.

79 ASTRONOMY AND ASTROPHYSICS↗