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 541 records · Page 30

Quasi-steady solar wind dynamics

Progress in understanding the large scale dynamics of quasisteady, corotating solar wind structure was reviewed. The nature of the solar wind at large heliocentric distances preliminary calculations from a 2-D MHD model are used to demonstrate theoretical expectations of corotating structure out to 30 AU. It is found that the forward and reverse shocks from adjacent CIR's begin to interact at about 10 AU, producing new shock pairs flanking secondary CIR's. These sawtooth secondary CIR's interact again at about 20 AU and survive as visible entities to 30 AU. The model predicts the velocity jumps at the leading edge of the secondary CIR's at 30 AU should be very small but there should still be sizable variations in the thermodynamic and magnetic parameters. The driving dynamic mechanism in the distant solar wind is the relaxation of pressure gradients. The second topic is the influence of weak, nonimpulsive time dependence in quasisteady dynamics. It is suggested that modest large scale variations in the coronal flow speed on periods of several hours to a day may be responsible for many of the remaining discrepancies between theory and observation. Effects offer a ready explanation for the apparent rounding of stream fronts between 0.3 and 1.0 AU discovered by Helios.

Pizzo, V. J.↗

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗

From Weibel seeds to dynamo beyond pair-plasmas

Bridging the spatiotemporal scales of magnetic seed field generation and subsequent dynamo amplification in the weakly collisional intracluster medium presents an extreme numerical challenge. We perform collisionless turbulence simulations with initially unmagnetised electrons that capture both magnetic seed generation via the electron Weibel instability and the ensuing dynamo amplification. Going beyond existing pair-plasma studies, we use an ion-to-electron mass ratio of 100 for which we find electron and ion dynamics are sufficiently decoupled. These simulations are enabled by the 10-moment collisionless fluid solver of Gkeyll, which evolves the full pressure tensor for all species. The electron heat-flux closure regulates pressure isotropisation and effectively sets the magnetic Reynolds number. We investigate how the strength of the closure influences the transition between a regime reminiscent of previous kinetic pair-plasma simulations and a regime exhibiting dynamo behaviour qualitatively similar to magnetohydrodynamics.

astrophysical plasmas↗

Imaging mechanism for hyperspectral scanning probe microscopy via Gaussian process modelling

We investigate the ability to reconstruct and derive spatial structure from sparsely sampled 3D piezoresponse force microcopy data, captured using the band-excitation (BE) technique, via Gaussian Process (GP) methods. Even for weakly informative priors, GP methods allow unambiguous determination of the characteristic length scales of the imaging process both in spatial and frequency domains. We further show that BE data set tends to be oversampled in the spatial domains, with ~30% of original data set sufficient for high-quality reconstruction, potentially enabling faster BE imaging. At the same time, reliable reconstruction along the frequency domain requires the resonance peak to be within the measured band. This behavior suggests the optimal strategy for the BE imaging on unknown samples. Finally, we discuss how GP can be used for automated experimentation in SPM, by combining GP regression with non-rectangular scans.

47 OTHER INSTRUMENTATION↗

Mass calibration of DES Year-3 clusters via SPT-3G CMB cluster lensing

We measure the stacked lensing signal in the direction of galaxy clusters in the DarkEnergy Survey Year 3 (DES Y3) redMaPPer sample, using cosmic microwave background (CMB)temperature data from SPT-3G, the third-generation CMB camera on the South PoleTelescope (SPT). Here, we estimate the lensing signal using temperature maps constructed fromthe initial 2 years of data from the SPT-3G 'Main' survey, covering 1500 deg$^{2}$ of the Southernsky. We then use this lensing signal as a proxy for the mean cluster mass of the DES sample. Thethermal Sunyaev-Zel'dovich (tSZ) signal, which can contaminate the lensing signal if notaddressed, is isolated and removed from the data before obtaining the mass measurement. In thiswork, we employ three versions of the redMaPPer catalogue: a Flux-Limited sample containing 8865clusters, a Volume-Limited sample with 5391 clusters, and a Volume&Redshift-Limited sample with4450 clusters. For the three samples, we detect the CMB lensing signal at a significance of12.4σ, 10.5σ and 10.2σ and find the mean cluster masses to be M$_{200m}$ = 1.66±0.13 [stat.]± 0.03 [sys.], 1.97±0.18 [stat.]± 0.05 [sys.],and 2.11±0.20 [stat.]± 0.05 [sys.]×10$^{14}$ M$_{⊙}$, respectively. Thisis a factor of ∼ 2 improvement relative to the precision of measurements with previousgenerations of SPT surveys and the most constraining cluster mass measurements using CMB clusterlensing to date. Overall, we find no significant tensions between our results and masses given byredMaPPer mass-richness scaling relations of previous works, which were calibrated using CMBcluster lensing, optical weak lensing, and velocity dispersion measurements from variouscombinations of DES, SDSS and Planck data. We then divide our sample into 3 redshift and 3richness bins, finding no significant discrepancies with optical weak-lensing calibrated masses inthese bins. We forecast a 5.7% constraint on the mean cluster mass of the DES Y3 sample withthe complete SPT-3G surveys when using both temperature and polarization data and including anadditional ∼ 1400 deg$^{2}$ of observations from the 'Extended' SPT-3G survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxy zoo: stronger bars facilitate quenching in star-forming galaxies

ABSTRACT We have used Galaxy Zoo DECaLS (GZD) to study strong and weak bars in disc galaxies. Out of the 314 000 galaxies in GZD, we created a volume-limited sample (0.01 < z < 0.05, Mr < − 18.96) which contains 1867 galaxies with reliable volunteer bar classifications in the ALFALFA footprint. In keeping with previous Galaxy Zoo surveys (such as GZ2), the morphological classifications from GZD agree well with previous morphological surveys. GZD considers galaxies to either have a strong bar (15.5 per cent), a weak bar (28.1 per cent) or no bar (56.4 per cent), based on volunteer classifications on images obtained from the DECaLS survey. This places GZD in a unique position to assess differences between strong and weak bars. We find that the strong bar fraction is typically higher in quiescent galaxies than in star-forming galaxies, while the weak bar fraction is similar. Moreover, we have found that strong bars facilitate the quenching process in star-forming galaxies, finding higher fibre star formation rates (SFRs), lower gas masses, and shorter depletion time-scales in these galaxies compared to unbarred galaxies. However, we also found that any differences between strong and weak bars disappear when controlling for bar length. Based on this, we conclude that weak and strong bars are not fundamentally different phenomena. Instead, we propose that there is a continuum of bar types, which varies from ‘weakest’ to ‘strongest’.

79 ASTRONOMY AND ASTROPHYSICS↗

Advances in constraining intrinsic alignment models with hydrodynamic simulations

We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans four redshift slices over the approximate range of contemporary lensing surveys z = 0–1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ~2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A 2 , in TNG, above statistical uncertainties, although MBII favours a moderately negative A 2 ~ –2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.

79 ASTRONOMY AND ASTROPHYSICS↗

Low-temperature T -linear resistivity due to umklapp scattering from a critical mode

In this study, we consider the transport properties of a model of fermions scattered by a critical bosonic mode. The mode is overdamped and scattering is mainly in the forward direction. Such a mode appears at the quantum critical point for an electronic nematic phase transition and in gauge theories for a U(1) spin liquid. It leads to a short fermion lifetime, violating Landau's criterion for a Fermi liquid. In spite of this, transport can be described by a Boltzmann equation. We include momentum relaxation by umklapp scattering, supplemented by weak impurity scattering. We find that above a very low temperature which scales with Δ q 3 , where Δ q is the minimum umklapp scattering vector, the resistivity is linear in T with a coefficient which is independent of the amount of disorder. This picture holds until the temperature dependent part of the resistivity exceeds that due to impurity scattering. We compare the relaxation time approximation with an exact numerical solution of the Boltzmann equation. Surprisingly we find that, unlike the resistivity, the Hall coefficient strongly deviates from the relaxation time approximation and shows a strong reduction with increasing temperature. We comment on possible comparisons with experiments on high T c cuprates.

36 MATERIALS SCIENCE↗

Quantum Computing for Neutrino-Nucleus Scattering

Neutrino-nucleus cross section uncertainties are expected to be a dominant systematic in future accelerator neutrino experiments. The cross sections are determined by the linear response of the nucleus to the weak interactions of the neutrino, and are dominated by energy and distance scales of the order of the separation between nucleons in the nucleus. These response functions are potentially an important early physics application of quantum computers. Here we present an analysis of the resources required and their expected scaling for scattering cross section calculations. The current estimates of Trotter steps needed to achieve an energy resolution of 10 MeV and the number of CNOT gates for analyzing Ar40 highlights the need for significant improvements in algorithms. We also examine simple small-scale neutrino-nucleus models on modern quantum hardware. In this paper, we use variational methods to obtain the ground state of a three nucleon system (the triton) and then implement the relevant time evolution. To tame the errors in present-day NISQ devices, we explore the use of different error-mitigation techniques to increase the fidelity of the calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Estimating the energy threshold of phonon-mediated superconducting qubit detectors operated in an energy-relaxation sensing scheme

In recent years, the lack of a conclusive detection of weakly interacting massive particle dark matter at the 10 GeV/c 2 mass scale and above has encouraged development of low-threshold detector technology aimed at probing lighter dark matter candidates. Detectors based on Cooper-pair-breaking sensors have emerged as a promising avenue for this detection due to the low (meV-scale) energy required for breaking a Cooper pair in most superconductors. Among them, devices based on superconducting qubits are interesting candidates for sensing due to their observed sensitivity to broken Cooper pairs. Here, we have developed an end-to-end G 4 CMP -based simulation framework and have used it to evaluate performance metrics of qubit-based devices operating in a gate-based “energy relaxation” readout scheme, akin to those used in recent studies of qubit sensitivity to ionizing radiation. We find that for this readout scheme, the qubit acts as a phonon sensor with an energy threshold ranging down to ≃ 0.4 eV for near-term performance parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the Jet–Ejecta Interaction in 3D GRMHD Simulations of a Binary Neutron Star Merger Aftermath

Short γ-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta's isotropic equivalent mass along the pole at the time of the BH formation is ≲10 -4 M ⊙ , setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳1015 G; and (iii) if the jets are weak (≲10 50 erg), the ejecta's isotropic equivalent mass along the pole must be small (≲10 -2 M ⊙ ). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ~1.5).

Jets↗

Qualitative properties of mathematical model for data flow

In this paper, properties of a recently proposed mathematical model for data flow in large-scale asynchronous computer systems are analyzed. In particular, the existence of special weak solutions based on propagating fronts is established. Qualitative properties of these solutions are investigated, both theoretically and numerically.

97 MATHEMATICS AND COMPUTING↗

Spectroscopic observations of winds on Venus

Measurements of the differential Doppler shift between various points on Venus using a high-resolution PEPSIOS interferometer. Using both a CO2 line and a Fraunhofer line a mean zonal-wind velocity of about -83 m/s is found near the equator; the velocity appears to vary from about -2 to -125 m/s with a time scale of greater than one week. Meridional velocities are measured to be weak. The equatorial zonal velocity appears to be smaller in the 'morning' than in the 'afternoon'. A comparison with reported velocities of the ultraviolet dark markings reveals general agreement in that both find the motion to be retrograde, variable, and accelerated during the day. A new potential source of systematic error in all spectroscopic determinations of the differential Doppler shift of nonuniformly illuminated objects is pointed out.

Traub, W. A.↗

The scattering of cosmic rays by magnetic bubbles

During large flares, the sun ejects magnetic bubbles which are then convected by the solar wind into the interstellar medium. The bubbles partake of the spherical expansion of the solar wind and grow to spatial scales comparable to the gyroradii of 100- to 1000-GeV cosmic rays. Assuming that most stars leave behind a wake of magnetic bubbles in the interstellar medium and that the rate of bubble production is comparable to the solar rate, 100- to 1000-GeV and possibly higher-energy cosmic rays will be well scattered and have mean free paths consistent with the compound diffusion theory. Particular choices for the scale-length spectrum of bubbles yield mean free paths which depend only weakly on cosmic-ray energy.

Flewelling, R. F.↗

Eolian sedimentation on earth and Mars - Some comparisons

Eolian sediments on earth are mostly formed from quartz. The quartz particles originally came from a granitic source. With respect to eolian sediments on Mars, it appears that an entirely different set of criteria must apply, but some critical parameters can be usefully compared. Impact experiments with basalt in eolian abrasion devices suggest that basalt sand-sized particles fragment rapidly to produce silt and clay-sized detritus. Cohesive forces must be more effective on Mars since the gravitational contribution to the bond/weight ratio is lower. Compared to the terrestrial situation, both larger and smaller particles can be expected to make significant contributions to eolian sediments on Mars. The low gravity and the high speed of moving particles and the relative weak rock material of which they are composed will allow large-scale fine particle production.

Smalley, I. J.↗

Galactic distribution models of gamma-ray burst sources

The log N-log S relationships for several galactic distribution models of gamma-ray burst sources are derived. The resulting curves are substantially different from the galactic curves usually assumed (slope = -1). There is a relatively weak dependence on the radial distribution assumed but a strong dependence on scale height. Available observational data from several experiments appear to rule out both thin disk distributions and extended halo distributions which have a strong radial density dependence. The prospects are good for obtaining self-consistent distribution data from future large-area detectors.

Fishman, G. J.↗

Strong turbulence and atmospheric waves in stellar occultations

General techniques for producing model lightcurves for a variety of realistic atmospheric irregularities, including turbulence and inertia-gravity waves, are presented and applied. The restrictions of weak scintillation theory are relaxed and model lightcurves are constructed using wave optics for atmospheres with strong, anisotropic turbulence. This is accomplished by numerical simulations which model the propagation of a wave through a phase-changing screen while maintaining complete amplitude and phase information from the wave. The results are compared with available weak scintillation theory and with recent occultation data. The effects of large scale atmospheric waves with realistic horizontal structure are examined, and the reliability of the numerical inversion method in retrieving the true atmospheric vertical structure under conditions of strong ray crossing and horizontal inhomogeneities is assessed. The nature of model lightcurve spikes generated using geometric optics and wave optics are compared.

French, R. G.↗