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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Coronal Heating, Weak MHD Turbulence, and Scaling Laws

Long-time high-resolution simulations of the dynamics of a coronal loop in Cartesian geometry are carried out, within the framework of reduced magnetohydrodynamics (RMHD), to understand coronal heating driven by the motion of field lines anchored in the photosphere. We unambiguously identify MHD anisotropic turbulence as the physical mechanism responsible for the transport of energy from the large scales, where energy is injected by photospheric motions, to the small scales, where it is dissipated. As the loop parameters vary, different regimes of turbulence develop: strong turbulence is found for weak axial magnetic fields and long loops, leading to Kolmogorov-like spectra in the perpendicular direction, while weaker and weaker regimes (steeper spectral slopes of total energy) are found for strong axial magnetic fields and short loops. As a consequence we predict that the scaling of the heating rate with axial magnetic field intensity B, which depends on the spectral index of total energy for given loop parameters, must vary from B3/2 for weak fields to B2 for strong fields at a given aspect ratio. The predicted heating rate is within the lower range of observed active region and quiet-Sun coronal energy losses.

corona↗
NASA NTRS2004

Narrow Scale Flow and a Weak Field by the Top of Earth's Core: Evidence from Orsted, Magsat and Secular Variation

As Earth's main magnetic field weakens, our magnetic shield against the onslaught of the solar wind thins. And the field strength needed to fend off battering by solar coronal mass ejections is decreasing, just when the delicate complexity of modem, vulnerable, electro-technological systems is increasing at an unprecedented rate. Recently, a working group of distinguished scientist from across the nation has asked NASA's Solid Earth and Natural Hazards program a key question: What are the dynamics of Earth s magnetic field and its interactions with the Earth system? Paleomagnetic studies of crustal rocks magnetized in the geologic past reveal that polarity reversals have occurred many times during Earth s history. Networked super-computer simulations of core field and flow, including effects of gravitational, pressure, rotational Coriolis, magnetic and viscous forces, suggest how this might happen in detail. And space-based measurements of the real, time-varying magnetic field help constrain estimates of the speed and direction of fluid iron flowing near the top of the core and enable tests of some hypotheses about such flow. Now scientists at NASA s Goddard Space Flight Center have developed and applied methods to test the hypotheses of narrow scale flow and of a dynamically weak magnetic field near the top of Earth s core. Using two completely different methods, C. V. Voorhies has shown these hypotheses lead to specific theoretical forms for the "spectrum" of Earth s main magnetic field and the spectrum of its rate of change. Much as solar physicists use a prism to separate sunlight into its spectrum, from long wavelength red to short wavelength blue light, geophysicists use a digital prism, spherical harmonic analysis, to separate the measured geomagnetic field into its spectrum, from long to short wavelength fields. They do this for the rate of change of the field as well.

Voorhies, Coerte V.↗
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

ALEGRA Parallel Scaling for Shock in a Heterogeneous Structure

We investigate the strong and weak parallel scaling performance of the ALEGRA multiphysics finite element program when solving a problem involving shock propagation through a heterogeneous material. We determine that ALEGRA scales well over a wide range of problem sizes, cores, and element sizes, and that scaling generally improves as the minimum element size in the mesh increases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Parthenon—a performance portable block-structured adaptive mesh refinement framework

On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multidimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of 1.7 × 10 13 zone-cycles/s on 9216 nodes (73,728 logical GPUs) at [Formula: see text] weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. Here, two-level GDSW (Generalized Dryja–Smith–Widlund) type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the MPI-parallel implementation of multi-level Schwarz preconditioners provided by the package FROSch (Fast and Robust Schwarz)from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To our knowledge, this is the first time two-level Schwarz preconditioners are applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The pre-conditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the sub-domains. Several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored as well. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as non uniform meshes for the Greenland ice sheet are considered. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32 K processor cores (8 K MPI-ranks and 4 OpenMP threads) and 566 M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4 K processor cores (and MPI-ranks) and 68 M degrees of freedom for the coupled problem.

58 GEOSCIENCES↗
DOE OSTI2024

Testing the scalar weak gravity conjecture in no-scale supergravity

We explore possible extensions of the Weak Gravity Conjecture (WGC) to scalar field theories. To avoid charged black hole remnants, the WGC requires the existence of a particle with a mass m ≤ gqM P , with charge q and U(1) gauge coupling g, allowing the decay to shed the black hole charge. Although there is no obvious problem that arises in the absence of a U(1) charge, it has been postulated that gravity must remain the weakest force even when extended to scalar interactions. Quantifying this conjecture may be done by comparing scalar and gravitational amplitudes, or as we advocate here by comparing scattering cross sections. In theories with non-trivial field space geometries, by working out examples with perturbation theory around arbitrary field values and performing tadpole resummations, we argue that the conjecture must be applied only at extrema of the scalar potential (when expressed in locally canonical coordinates). We consider several toy models in the context of no-scale supergravity and also consider examples of inflationary models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1984

Comments on 'Comparative diagnostic case study of east coast secondary cyclogenesis under weak versus strong synoptic-scale forcing'

A comment is presented on a recent paper by Pagnotti and Bosart (1984) which contrasted different forcing mechanisms for weak and strong cyclones along the east coast of the U.S., and Bosart replies. The comment argues that the paper did not give a clear definition of 'weak' upper-level forcing and did not provide sufficient evidence to support its claim that the return branch of transverse circulations may be coupled to the middle as opposed to the lower troposphere.

Uccellini, L. W.↗
NASA NTRS1989

The need for enhanced initial moisture information in simulations of a complex summertime precipitation event

Initial simulations of the June 28, 1986 Cooperative Huntsville Meteorological Experiment case illustrate the need for mesoscale moisture information in a summertime situation in which deep convection is organized by weak large scale forcing. A methodology is presented for enhancing the initial moisture field from a combination of IR satellite imagery, surface-based cloud observations, and manually digitized radar data. The Mesoscale Atmospheric Simulation Model is utilized to simulate the events of June 28-29. This procedure insures that areas known to have precipitation at the time of initialization will be nearly saturated on the grid scale, which should decrease the time needed by the model to produce the observed Bonnie (a relatively weak hurricane that moved on shore two days before) convection. This method will also result in an initial distribution of model cloudiness (transmissivity) that is very similar to that of the IR satellite image.

Waight, Kenneth T., III↗
NASA NTRS1994

Large-scale structure effects on the gravitational lens image positions and time delay

We compute the fluctuations in gravitational lens image positions and time delay caused by large-scale structure correlations. We show that these fluctuations can be expressed as a simple integral over the density power spectrum. Using the Cosmic Background Explorer (COBE) normalization we find that positions of objects at cosmological distances are expected to deviate from their true positions by few arcminutes. These deflections are not directly observable. The positions of the images relative to one another fluctuate by a few percent of the relative separation, implying that one does not expect multiple images to be produced by large-scale structure. Nevertheless, the fluctuations are larger than the observational errors on the positions and affect reconstructions of the lens potential. The time delay fluctuations have a geometrical and a gravitational contribution. Both are much larger than the expected time delay from the primary lens, but partially cancel each other. We find that large-scale structure weakly affects the time delay and time delay measurements can be used as a probe of the distance scale in the universe.

Seljak, Uros↗
NASA NTRS1991

Model studies on the role of moist convection as a mechanism for interaction between the mesoscales

A three year research effort is described which had as its goal the development of techniques to improve the numerical prediction of cumulus convection on the meso-beta and meso-gamma scales. Two MESO models are used, the MASS (mesoscale) and TASS (cloud scale) models. The primary meteorological situation studied is the 28-29 Jun. 1986 Cooperative Huntsville Meteorological Experiment (COHMEX) study area on a day with relatively weak large scale forcing. The problem of determining where and when convection should be initiated is considered to be a major problem of current approaches. Assimilation of moisture data from satellite, radar, and surface data is shown to significantly improve mesoscale simulations. The TASS model is shown to reproduce some observed mesoscale features when initialized with 3-D observational data. Convection evolution studies center on comparison of the Kuo and Fritsch-Chappell cumulus parameterization schemes to each other, and to cloud model results. The Fritsch-Chappell scheme is found to be superior at about 30 km resolution, while the Kuo scheme does surprisingly well in simulating convection down to 10 km in cases where convergence features are well-resolved by the model grid. Results from MASS-TASS interaction experiments are presented and discussed. A discussion of the future of convective simulation is given, with the conclusion that significant progress is possible on several fronts in the next few years.

Waight, Kenneth T., III↗
NASA NTRS1998

Internal Dynamics and Boundary Forcing Characteristics Associated with Interannual Variability of the Asian Summer Monsoon

In this paper, we present a description of the internal dynamics and boundary forcing characteristics of two major components of the Asian summer monsoon (ASM), i.e., the South Asian (SAM) and the Southeast-East Asian monsoon (SEAM). The description is based on a new monsoon-climate paradigm in which the variability of ASM is considered as the outcome of the interplay of a "fast" and an "intermediate" monsoon subsystem, under the influenced of the "slow" varying external forcings. Two sets of regional monsoon indices derived from dynamically consistent rainfall and wind data are used in this study. For SAM, the internal dynamics is represented by that of a "classical" monsoon system where the anomalous circulation is governed by Rossby-wave dynamics, i.e., generation of anomalous vorticity induced by an off-equatorial heat source is balanced by planetary vorticity advection. On the other hand, the internal dynamics of SEAM is characterized by a "hybrid" monsoon system featuring multi-cellular meridional circulation over the East Asian section, extending from the deep tropics to midlatitudes. These meridional-cells link tropical heating to extratropical circulation system via the East Asian jetstream, and are responsible for the characteristic occurrences of zonally oriented anomalous rainfall patterns over East Asian and the subtropical western Pacific. In the extratropical regions, the major upper level vorticity balance is by anomalous vorticity advection and generation by the anomalous divergent circulation. A consequence of this is that compared to SAM, the SEAM is associated with stronger teleconnection patterns to regions outside the ASM. A strong SAM is linked to basin-scale sea surface temperature (SST) fluctuation with significant signal in the equatorial eastern Pacific. During the boreal spring SST warming in the Arabian Sea and the subtropical western Pacific may lead to a strong SAM. For SEAM, interannual variability is tied to SSTA over the Sea of Japan and the South China Sea regions, while the linkage to equatorial basin-scale SSTA is weak at best. A large scale SSTA dipole with warming (cooling) in the subtropical central (eastern) Pacific foreshadows a strong SEAM.

Lau, K.- M.↗
DOE OSTI2023

Multi-scale Dynamics of Kinetic Turbulence in Weakly Collisional, High-Beta Plasmas (Final Report for DOE Grant DE-SC0019046)

This grant was a collaborative grant between Prof Matt Kunz (Princeton) and Prof. Eliot Quataert. Quataert was initially a faculty member at UC Berkeley when the grant was funded but moved to Princeton during the timeframe of this grant. Over the course of this grant we made major progress on understanding turbulence and multi-scale dynamics in weakly collisional high-beta plasmas. The most important contributions included: We developed one of the most compelling theoretical explanations for a decades-old puzzle at the heart of our understanding of the origin of the solar wind. The puzzle is that many of the observations of ion temperatures and distribution functions in the solar wind are consistent with cyclotron resonant heating. However, theoretical models of MHD turbulence in the solar wind show that most of the turbulent energy remains at low frequencies below the cyclotron frequency. In Squire et al. (2022), we showed, however, that in imbalanced turbulence (in which there is an asymmetry in the Alfven-wave flux in opposite directions along the magnetic field), which is the norm in the fast solar wind, the turbulent energy reaches a ‘bottleneck’ near the ion Larmor radius, and the amplitude and characteristic frequency of the fluctuations grow until the cyclotron frequency is reached.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The Aemulus Project. V. Cosmological Constraint from Small-scale Clustering of BOSS Galaxies

We analyze clustering measurements of BOSS galaxies using a simulation-based emulator of two-point statistics. We focus on the monopole and quadrupole of the redshift-space correlation function, and the projected correlation function, at scales of 0.1 ~ 60 h -1 Mpc. Although our simulations are based on wCDM with general relativity (GR), we include a scaling parameter of the halo velocity field, γ f , defined as the amplitude of the halo velocity field relative to the GR prediction. We divide the BOSS data into three redshift bins. After marginalizing over other cosmological parameters, galaxy bias parameters, and the velocity scaling parameter, we find fσ 8 (z = 0.25) = 0.413 ± 0.031, fσ 8 (z = 0.4) = 0.470 ± 0.026, and fσ 8 (z = 0.55) = 0.396 ± 0.022. Compared with Planck observations using a flat Lambda cold dark matter model, our results are lower by 1.9σ, 0.3σ, and 3.4σ, respectively. These results are consistent with other recent simulation-based results at nonlinear scales, including weak lensing measurements of BOSS LOWZ galaxies, two-point clustering of eBOSS LRGs, and an independent clustering analysis of BOSS LOWZ. All these results are generally consistent with a combination of ${\gamma }_{f}^{1/2}{\sigma }_{8}\approx 0.75$. We note, however, that the BOSS data is well fit assuming GR, i.e., γf = 1. We cannot rule out an unknown systematic error in the galaxy bias model at nonlinear scales, but near-future data and modeling will enhance our understanding of the galaxy–halo connection, and provide a strong test of new physics beyond the standard model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Forecasting the potential of weak lensing magnification to enhance LSST large-scale structure analyses

ABSTRACT Recent works have shown that weak lensing magnification must be included in upcoming large-scale structure analyses, such as for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), to avoid biasing the cosmological results. In this work, we investigate whether including magnification has a positive impact on the precision of the cosmological constraints, as well as being necessary to avoid bias. We forecast this using an LSST mock catalogue and a halo model to calculate the galaxy power spectra. We find that including magnification has little effect on the precision of the cosmological parameter constraints for an LSST galaxy clustering analysis, where the halo model parameters are additionally constrained by the galaxy luminosity function. In particular, we find that for the LSST gold sample (i < 25.3) including weak lensing magnification only improves the galaxy clustering constraint on Ωm by a factor of 1.03, and when using a very deep LSST mock sample (i < 26.5) by a factor of 1.3. Since magnification predominantly contributes to the clustering measurement and provides similar information to that of cosmic shear, this improvement would be reduced for a combined galaxy clustering and shear analysis. We also confirm that not modelling weak lensing magnification will catastrophically bias the cosmological results from LSST. Magnification must therefore be included in LSST large-scale structure analyses even though it does not significantly enhance the precision of the cosmological constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Classical combinatorial optimization scaling for random Ising models on 2D heavy-hex graphs

Motivated by near term quantum computing hardware limitations, combinatorial optimization problems that can be addressed by current quantum algorithms and noisy hardware with little or no overhead are used to probe capabilities of quantum algorithms such as the quantum approximate optimization algorithm. In this study, a specific class of near term quantum computing hardware defined combinatorial optimization problems, Ising models on heavy-hex graphs both with and without geometrically local cubic terms, are examined for their classical computational hardness via empirical computation time scaling quantification. Specifically the time-to-solution (TTS) metric using the classical heuristic simulated annealing is measured for finding optimal variable assignments (ground states), as well as the time required for the optimization software Gurobi to find an optimal variable assignment. Because of the sparsity of these Ising models, the classical algorithms are able to find optimal solutions efficiently even for large instances (i.e. 100 000 spin variables). The Ising models both with and without geometrically local cubic terms exhibit average-case linear-time or weakly quadratic scaling when solved exactly using Gurobi, and the Ising models with no cubic terms show evidence of exponential-time TTS scaling when sampled using simulated annealing. These findings point to the necessity of developing and testing more complex, namely more densely connected, optimization problems in order for quantum computing to ever have a practical advantage over classical computing. Our results are another illustration that different classical algorithms can indeed have exponentially different running times, thus making the identification of the best practical classical technique important in any quantum computing vs. classical computing comparison.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗