Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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 595 records · Page 33

Simulation of Boundary-Layer Cumulus and Stratocumulus Clouds using a Cloud-Resolving Model With Low- and Third-Order Turbulence Closures

The effects of subgrid-scale condensation and transport become more important as the grid spacings increase from those typically used in large-eddy simulation (LES) to those typically used in cloud-resolving models (CRMs). Incorporation of these effects can be achieved by a joint probability density function approach that utilizes higher-order moments of thermodynamic and dynamic variables. This study examines how well shallow cumulus and stratocumulus clouds are simulated by two versions of a CRM that is implemented with low-order and third-order turbulence closures (LOC and TOC) when a typical CRM horizontal resolution is used and what roles the subgrid-scale and resolved-scale processes play as the horizontal grid spacing of the CRM becomes finer. Cumulus clouds were mostly produced through subgrid-scale transport processes while stratocumulus clouds were produced through both subgrid-scale and resolved-scale processes in the TOC version of the CRM when a typical CRM grid spacing is used. The LOC version of the CRM relied upon resolved-scale circulations to produce both cumulus and stratocumulus clouds, due to small subgrid-scale transports. The mean profiles of thermodynamic variables, cloud fraction and liquid water content exhibit significant differences between the two versions of the CRM, with the TOC results agreeing better with the LES than the LOC results. The characteristics, temporal evolution and mean profiles of shallow cumulus and stratocumulus clouds are weakly dependent upon the horizontal grid spacing used in the TOC CRM. However, the ratio of the subgrid-scale to resolved-scale fluxes becomes smaller as the horizontal grid spacing decreases. The subcloud-layer fluxes are mostly due to the resolved scales when a grid spacing less than or equal to 1 km is used. The overall results of the TOC simulations suggest that a 1-km grid spacing is a good choice for CRM simulation of shallow cumulus and stratocumulus.

Xu, Kuan-Man↗

Finding evidence for inflation and the origin of galactic magnetic fields with CMB surveys

The origin of the μG magnetic fields observed in galaxies is unknown. One promising scenario is that magnetic fields generated during inflation, larger than 0.1 nG on Mpc scales, were adiabatically compressed to μG strengths in galaxies during structure formation. Thus, detecting a scale-invariant primordial magnetic field (PMF) above 0.1 nG on Mpc scales just after recombination would indicate an inflationary origin of galactic magnetic fields. Furthermore, this would also provide compelling evidence that inflation occurred since only an inflationary mechanism could generate such a strong, scale-invariant magnetic field on Mpc scales. In contrast, constraining the scale-invariant PMF strength to be below 0.1 nG would imply an inflationary scenario is not the primary origin, since such weak PMFs cannot be amplified enough via adiabatic compression to produce the strength of the galactic fields we observe today. We find that measurements of anisotropic birefringence by future CMB surveys will be able to improve the sensitivity to Mpc-scale inflationary PMFs by an order of magnitude, and, in particular, that CMB-HD would lower the upper bound to 0.072 nG at the 95% C.L., which is below the critical 0.1 nG threshold for ruling out a purely inflationary origin. If inflationary PMFs exist, we find that a CMB-HD survey would be able to detect them with about 3σ significance or higher, providing evidence for inflation itself.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

2D-FFTLog: efficient computation of real-space covariance matrices for galaxy clustering and weak lensing

ABSTRACT Accurate covariance matrices for two-point functions are critical for inferring cosmological parameters in likelihood analyses of large-scale structure surveys. Among various approaches to obtaining the covariance, analytic computation is much faster and less noisy than estimation from data or simulations. However, the transform of covariances from Fourier space to real space involves integrals with two Bessel integrals, which are numerically slow and easily affected by numerical uncertainties. Inaccurate covariances may lead to significant errors in the inference of the cosmological parameters. In this paper, we introduce a 2D-FFTLog algorithm for efficient, accurate, and numerically stable computation of non-Gaussian real-space covariances for both 3D and projected statistics. The 2D-FFTLog algorithm is easily extended to perform real-space bin-averaging. We apply the algorithm to the covariances for galaxy clustering and weak lensing for a Dark Energy Survey Year 3-like and a Rubin Observatory’s Legacy Survey of Space and Time Year 1-like survey, and demonstrate that for both surveys, our algorithm can produce numerically stable angular bin-averaged covariances with the flat sky approximation, which are sufficiently accurate for inferring cosmological parameters. The code CosmoCov for computing the real-space covariances with or without the flat-sky approximation is released along with this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Intruder configurations in 29 Ne at the transition into the island of inversion: Detailed structure study of 28 Ne

Detailed γ-ray spectroscopy of the exotic neon isotope 28 Ne has been performed for the first time using the one-neutron removal reaction from 29 Ne on a liquid hydrogen target at 240 MeV/nucleon. Based on an analysis of parallel momentum distributions, a level scheme with spin-parity assignments has been constructed for 28 Ne and the negative-parity states are identified for the first time. The measured partial cross sections and momentum distributions reveal a significant intruder p-wave strength providing evidence of the breakdown of the N=20 and N=28 shell gaps. Only a weak, possible f-wave strength was observed to bound final states. Large-scale shell-model calculations with different effective interactions do not reproduce the large p-wave and small f-wave strength observed experimentally, indicating an ongoing challenge for a complete theoretical description of the transition into the island of inversion along the Ne isotopic chain.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Misaligned Wind‐Waves Behind Atmospheric Cold Fronts

Abstract Atmospheric fronts embedded in extratropical cyclones are high‐impact weather phenomena, contributing significantly to mid‐latitude winter precipitation. The three vital characteristics of the atmospheric fronts, high wind speeds, abrupt change in wind direction, and rapid translation, force the induced surface waves to be misaligned with winds exclusively behind the cold fronts. The effects of the misaligned waves under atmospheric cold fronts on air‐sea fluxes remain undocumented. Using the multi‐year in situ near‐surface observations and direct covariance flux measurements from the Pioneer Array off the coast of New England, we find that the majority of the passing cold fronts generate misaligned waves behind the cold front. Once generated, the waves remain misaligned, on average, for about 8 hr. The parameterized effect of misaligned waves in a fully coupled model significantly increases the roughness length (185%), drag coefficient (19%), and air‐sea momentum flux (11%). The increased surface drag reduces the wind speeds in the surface layer. The upward turbulent heat flux is weakly decreased by the misaligned waves because of the decrease in temperature and humidity scaling parameters being greater than the increase in friction velocity. The misaligned wave effect is not accurately represented in a commonly used wave‐based bulk flux algorithm. Yet, considering this effect in the current formulation improves the overall accuracy of parameterized momentum flux estimates. The results imply that better representing a directional wind‐wave coupling in the bulk formula of the numerical models may help improve the air‐sea interaction simulations under the passing atmospheric fronts in the mid‐latitudes.

Oceanography↗

Detailed characterization of runaway electron driven whistler waves in low-density DIII-D discharges

RE-driven whistler waves during quiescent DIII-D shots have been investigated further. The waves are confirmed to be mostly perpendicularly propagating and are observed for the first time with frequencies up to 700 MHz. Phase-spectral analysis has been used to infer their toroidal mode numbers, n, which are expected to scale with the wavenumber, k, of the mode. Though we derive a theoretical scaling of k ≈ 4n, the measured mode numbers are found to exhibit a very weak dependence on k. In addition, increases in synchrotron emission have been found to consistently lag whistler wave bursts by roughly 3–5 ms, suggesting the waves are causing pitch-angle scattering, since the emitted synchrotron radiation is a strong function of the REs' perpendicular energy. The stronger the wave bursts, the greater the subsequent increase in synchrotron emission. A predator-prey model is used to describe these nonlinear wave-particle interactions, from which the wave damping rates and the loss parameter can be inferred. The damping rates are found to be of the order of (1.6 ± 0.8) × 10 4 /s, and the unitless loss parameter is found to be approximately 2, suggesting that the loss mechanism is diffusive. These observations will serve to validate models of RE-driven waves in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulence-level dependence of cosmic ray parallel diffusion

ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.

79 ASTRONOMY AND ASTROPHYSICS↗

AbacusHOD : a highly efficient extended multitracer HOD framework and its application to BOSS and eBOSS data

ABSTRACT We introduce the AbacusHOD model and present two applications of AbacusHOD and the AbacusSummit simulations to observations. AbacusHOD is a Halo Occupation Distribution (HOD) framework written in Python that is particle-based, multitracer, highly generalized, and highly efficient. It is designed specifically with multitracer/cosmology analyses for next-generation large-scale structure surveys in mind, and takes advantage of the volume and precision offered by the new state-of-the-art AbacusSummit cosmological simulations. The model is also highly customizable and should be broadly applicable to any upcoming surveys and a diverse range of cosmological analyses. In this paper, we demonstrate the capabilities of the AbacusHOD framework through two example applications. The first example demonstrates the high efficiency and the large HOD extension feature set through an analysis of full-shape redshift-space clustering of BOSS galaxies at intermediate to small scales ($\lt 30\, h^{-1}$ Mpc), assessing the necessity of introducing secondary galaxy biases (assembly bias). We find strong evidence for using halo environment instead of concentration to trace secondary galaxy bias, a result which also leads to a moderate reduction in the ‘lensing is low’ tension. The second example demonstrates the multitracer capabilities of the AbacusHOD package through an analysis of the extended Baryon Oscillation Spectroscopic Survey cross-correlation measurements between three different galaxy tracers: luminous red galaxies, emission-line galaxies, and quasi-stellar objects. We expect the AbacusHOD framework, in combination with the AbacusSummit simulation suite, to play an important role in a simulation-based analysis of the upcoming Dark Energy Spectroscopic Instrument data sets.

79 ASTRONOMY AND ASTROPHYSICS↗

Criteria for fragmentation in a collapsing rotating cloud

The isothermal collapse of a rotating protostar with pressure and self-gravity is calculated by a hydrodynamic computer code that treats three space dimensions. A number of initial conditions are tested to determine the conditions under which a cloud is unstable to fragmentation. It is shown that fragmentation does not proceed to a significant extent during the first free-fall time, but that in most cases fragmentation has occurred by 1.5 to 2 times the initial free-fall time. The extent to which the density in the fragments is enhanced over that in their surroundings on this time scale depends on the initial ratio of thermal to gravitational energy, but is only weakly dependent on the type and amplitude of the perturbation and on the initial rotational energy.

Bodenheimer, P.↗

Photometric Studies of Heat Flow at the Photosphere

Continuum photometry is carried out and the results of comparing these observations with models of photospheric heat flow are described. The main results are: (1) a possible detection of weak bright rings around some spot penumbrae; (2) no evidence is found for large scale photospheric brightness inhomogeneities exceeding 2-3 K which places tighter constraints on models of global scale convection; (3) supergranular scale continuum structures observed across the photosphere appear mainly due to random clumping of granules; (4) the one case observed of a sunspot emergence shows no thermal shadow exceeding 1.5 K rms one day prior to umbra appearance; (5) network and faculae are found to show a small excess brightness even at mu = 1, so detection of faculae at mu = 1 by differential photometry indicates a gentler temperature gradient near tau = 1 in the facular (relative to cell) atmosphere; (6) the limb darkening study shows no significant global variations to within 0.1% rms.

Foukal, P.↗

Numerical simulations of turbulent trapping in the weak beam-plasma instability

Numerical simulations of the weak beam-plasma instability were done in the turbulent regime where small-scale trapping is a dominant feature of the instability, a regime with behavior not predicted by quasi-linear theory. The results of the simulations were compared with those of a specific model of the turbulence, the so-called 'turbulent trapping' model, which gives precise formulas for the particle correlation functions, and predicts a growth rate well enhanced over the quasi-linear value. It was found that the model gives accurate predictions for the correlation functions. On the other hand, while growth rates were enhanced over the quasi-linear values, the enhancements observed are smaller than expected from the quantitative predictions of the model.

Theilhaber, K.↗

Strong plasma turbulence in the earth's electron foreshock

A quantitative model is developed to account for the distribution in magnitude and location of the intense plasma waves observed in the earth's electron foreshock given the observed rms levels of waves. In this model, nonlinear strong-turbulence effects cause solitonlike coherent wave packets to form and decouple from incoherent background beam-excited weak turbulence, after which they convect downstream with the solar wind while collapsing to scales as short as 100 m and fields as high as 2 V/m. The existence of waves with energy densities above the strong-turbulence wave-collapse threshold is inferred from observations from IMP 6 and ISEE 1 and quantitative agreement is found between the predicted distribution of fields in an ensemble of such wave packets and the actual field distribution observed in situ by IMP 6. Predictions for the polarization of plasma waves and the bandwidth of ion-sound waves are also consistent with the observations. It is shown that strong-turbulence effects must be incorporated in any comprehensive theory of the propagation and evolution of electron beams in the foreshock. Previous arguments against the existence of strong turbulence in the foreshock are refuted.

Robinson, P. A.↗

Wave diffraction in weak cosmic-ray-modified shocks

Weakly multidirectional, long-wavelength cosmic-ray-modified shocks are studied via multiple scale perturbation techniques. The effects of diffraction are discussed in terms of Green's function solutions of the linearized 1 + 3D Burgers and 1 + 3D KdVB equations, and also in terms of solutions with singular Dirac delta initial distributions. The solutions show a monotonic decrease of the wave-front curvature with increasing time owing to the effects of wave diffraction. The shape of the wave surface is discussed in terms of solutions S to the wave eikonal equation corresponding to singular initial conditions. For the fast magnetosonic wave propagating in the positive x-direction, the wave phase surface S = 0 has elliptic cross sections with the planes x = constant and has a convex paraboloidal shape. Plane-wave solutions of the 1 + 3D KdVB equation are discussed.

Webb, G. M.↗

Weakly inhomogeneous MHD turbulence and transport of solar wind fluctuations

An evaluation is conducted of recent theories of small-scale MHD turbulence transport in an inhomogeneous background that are pertinent to the evolution of solar wind turbulence. Attention is given to the WKB formalism that has been used in many solar wind-related physics applications, with a view to its shortcomings. Also discussed are the structure of two-scale transport theories, and their relationship to WKB theory in light of multiple-scales analysis.

Matthaeus, W. H.↗

Mixing in a stratified shear flow: Energetics and sampling

Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed for Richardson numbers from 0 to 1 and for Prandtl numbers between 0.1 and 2. The results indicate that mixing efficiency R(sub f) varies with turbulent Froude number in a manner consistent with laboratory experiments performed with Prandtl numbers of 0.7 and 700. However, unlike the laboratory results, for a particular Froude number, the simulations do not show a clear dependence on the magnitude of R(sub f) on Pr. The observed maximum value of R(sub f) is 0.25. When averaged over vertical length scales of an order of magnitude greater than either the overturning or Ozmidov scales of the flow, the simulations indicate that the dissipation rate epsilon is only weakly lognormally distributed with an intermittency of about 0.01 whereas estimated values in the ocean are 3 to 7.

Ivey, G. N.↗

Coupling between middle atmosphere trend estimates and solar effects in ozone vertical structure

Studies are performed on both the Stratospheric Aerosol and Gas Experiment (SAGE) II (1985-1989) and Solar Backscatter Ultraviolet Spectrometer (SBUV)-Version 6 (1979-1986) global ozone vertical structure satellite data sets to determine the long-term trends in ozone as a function of altitude (pressure) and latitude. SAGE II data are only available during the period of increasing solar activity and show increases in ozone with time in the upper stratosphere which are attributed largely to rising solar activity. Looking at this data set independently, the solar effects and trends are highly coupled and cannot be clearly separated. However, a study of combined SBUV and SAGE II data over the 11-year solar cycle shows a clear response of ozone to 11-year solar variations and allows a decoupling of solar effects, quasibiennial oscillations (QBO), and trends. The detailed pattern of long-term ozone trends become clear using this approach. In the upper stratosphere, ozone depletion increases sharply with latitude. Global trends are fairly symmetric about the equator but are somewhat stronger in the Southern Hemisphere. Near the equator, some layers of ozone are decreasing with time while others appear to be increasing. Near 30 mb, there is evidence of intrusion to mid latitudes of high latitude negative trends. Near 15 mb, trends appear to be very weak. Near the tropopause there appears to be strong ozone depletion on a global scale. Two regions of unexpectedly strong ozone response to 11-year solar variations were detected: the first near 2 mb and the second near 30 mb at low latitudes and near 15 mb at mid latitudes.

Keating, G. M.↗

X-Ray Emission from the Unidentified Gamma-Ray Transient GRO J1838-0415

The gamma-ray transient GRO J1838-04 is one of the most enigmatic sources in the sky. Despite its closeness to the Galactic plane, it showed a strong gamma-ray flare in 1995 typical of blazars. However, no blazar is detected in its approx. 1 deg. error box, and a Galactic source is suspected. Other time variable gamma-ray sources are known in the Galactic disk, and GRO J1838-04 may belong to a new class of sources (isolated young pulsars, the only proven Galactic gamma-ray sources, do not show variability). The Advanced Satellite for Cosmology and Astrophysics (ASCA) observed the centroid of the Energetic Gamma Ray Experiment Telescope (EGRET) error box of GRO J1838-04 in April 1997. Five weak X-ray sources are detected, one of which appears to be diffuse (3 arcmins). No known radio source is coincident with the X-ray sources. A preliminary report of our results appeared and a more systematic report is being written . The fact that no prominent X-ray source appears in the field, excludes an X-ray transient source with strong persistent emission is a possible counterpart of GRO JI838- 04 (such as superluminal transients). To uncover the nature of GRO J1838-04 requires more study of the weak X-ray sources in its error box, and a long-time scale monitoring of their variability.

Tavani, Marco↗

Light scalars at the cosmological collider

We study the self-energies of weakly interacting scalar fields in de Sitter space with one field much lighter than the Hubble scale. We argue that self-energies drastically simplify in this light limit. We illustrate this in theories with two scalar fields, one heavy and one light, interacting with one another through either cubic or quartic interactions. To regulate infrared divergences, we compute these self-energies in Euclidean de Sitter space and then carefully analytically continue to Lorentzian signature. In particular, we do this for the most general renormalizable theory of two scalar fields with even interactions to leading order in the coupling and the mass of the light field. These self-energies are determined by de Sitter sunset diagrams, whose analytic structure and UV divergences we derive. Even at very weak couplings, the light field can substantially change how the heavy field propagates over long distances. The light field’s existence may then be inferred from how it modifies the heavy field’s oscillatory contribution to the primordial bispectrum in the squeezed limit, i.e. its cosmological collider signal.

79 ASTRONOMY AND ASTROPHYSICS↗