Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Extreme scale”

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 289 records · Page 16

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

Hakobyan, Hayk [Flatiron Institute, New York, NY (↗

Flow instabilities in transonic small disturbance theory

The dynamics of unsteady transonic small disturbance flows about two-dimensional airfoils is examined, with emphasis on the behavior in the region where the steady state flow is nonunique. It is shown that nonuniqueness results from an extremely long time scale instability which occurs in a finite Mach number and angle of attack range. The similarity scaling rules for the instability are presented and the possibility of similar behavior in the Euler equations is discussed.

Williams, M. H.↗

Flow instabilities in transonic small-disturbance theory

The dynamics of unsteady transonic small disturbance flows about two-dimensional airfoils is examined, with emphasis on the behavior in the region where the steady state flow is nonunique. It is shown that nonuniqueness results from an extremely long time scale instability which occurs in a finite Mach number and angle of attack range. The similarity scaling rules for the instability are presented and the possibility of similar behvior in the Euler equations is discussed.

Williams, M. H.↗

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗

The solar origins of two high-latitude interplanetary disturbances

Two extremely similar interplanetary forward/reverse shock events, with bidirectional electron streaming were detected by Ulysses in 1994. Ground-based and Yohkoh/SXT observations show two strikingly different solar events that could be associated with them: an LDE flare on 20 Feb. 1994, and a extremely large-scale eruptive event on 14 April 1994. Both events resulted in geomagnetic storms and presumably were associated with coronal mass ejections. The sharply contrasting nature of these solar events argues against an energetic causal relationship between them and the bidirectional streaming events observed by Ulysses during its S polar passage. We suggest instead that for each pair of events. a common solar trigger may have caused independent instabilities leading to the solar and interplanetary phenomena.

Hudson, H. S.↗

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗

A Normal Incidence X-ray Telescope (NIXT) Sounding Rocket Payload

The solar corona, and the coronae of solar-type stars, consist of a low-density magnetized plasma at temperatures exceeding 10(exp 6) K. The primary coronal emission is therefore in the UV and soft X-ray range. The observed close connection between solar magnetic fields and the physical parameters of the corona implies a fundamental role for the magnetic field in coronal structuring and dynamics. Variability of the corona occurs on all temporal and spatial scales - at one extreme, as the result of plasma instabilities, and at the other extreme driven by the global magnetic flux emergence patterns of the solar cycle.

Golub, Leon↗

Climate Scenarios for the NASA / USAID SERVIR Project: Challenges for Multiple Planning Horizons

SERVIR, an acronym meaning "to serve" in Spanish, is a joint venture between NASA and the U.S. Agency for International Development (USAID) which provides satellite-based Earth observation data, modeling, and science applications to help developing nations in Central America, East Africa and the Himalayas improve environmental decision making. Anticipating climate variability / climate change impacts has now become an important component of the SERVIR efforts to build capacity in these regions. Uncertainty in hydrometeorological components of climate variations and exposure to extreme events across scales from weather to climate are of particular concern. We report here on work to construct scenarios or outlooks that are being developed as input drivers for decision support systems (DSSs) in a variety of settings. These DSSs are being developed jointly by a broad array NASA Applied Science Team (AST) Investigations and user communities in the three SERVIR Hub Regions, Central America, East Africa and the Himalayas. Issues span hydrologic / water resources modeling, agricultural productivity, and forest carbon reserves. The scenarios needed for these efforts encompass seasonal forecasts, interannual outlooks, and likely decadal / multi-decadal trends. Providing these scenarios across the different AST efforts enables some level of integration in considering regional responses to climate events. We will discuss a number of challenges in developing this continuum of scenarios including the identification and "mining" of predictability, addressing multiple continental regions, issues of downscaling global model integrations to regional / local applications (i.e. hydrologic and crop modeling). We compare / contrast the role of the U.S. National Multi- Model Experiment initiative in seasonal forecasts and the CMIP-5 climate model experiments in supporting these efforts. Examples of these scenarios, their use, and an assessment of their utility as well as limitations will be presented.

Robertson, Franklin R.↗

Dislocation Content Measured Via 3D HR-EBSD Near a Grain Boundary in an AlCu Oligocrystal

Interactions between dislocations and grain boundaries are poorly understood and crucial to mesoscale plasticity modeling. Much of our understanding of dislocation-grain boundary interaction comes from atomistic simulations and TEM studies, both of which are extremely limited in scale. High angular resolution EBSD-based continuum dislocation microscopy provides a way of measuring dislocation activity at length scales and accuracies relevant to crystal plasticity, but it is limited as a two-dimensional technique, meaning the character of the grain boundary and the complete dislocation activity is difficult to recover. However, the commercialization of plasma FIB dual-beam microscopes have made 3D EBSD studies all the more feasible. The objective of this work is to apply high angular resolution cross correlation EBSD to a 3D EBSD data set collected by serial sectioning in a FIB to characterize dislocation interaction with a grain boundary. Three dimensional high angular resolution cross correlation EBSD analysis was applied to an AlCu oligocrystal to measure dislocation densities around a grain boundary. Distortion derivatives associated with the plasma FIB serial sectioning were higher than expected, possibly due to geometric uncertainty between layers. Future work will focus on mitigating the geometric uncertainty and examining more regions of interest along the grain boundary to glean information on dislocation-grain boundary interaction.

Ruggles, Timothy↗

Advanced Energy Scale Correction Techniques for the X-ray Transition Edge Sensors of the Athena mission

The X-ray Integral Field Unit (X-IFU) onboard the future European X-ray telescope Athena will be the first space instrument carrying an array of more than a thousand transition edge sensors. One of the key challenges of the X-IFU is the measurement of narrow X-ray atomic lines to determine velocity shifts at an unprecedented level of accuracy. For this reason, the energy scale of the instrument needs to be known with extreme accuracy, of 0.4 eV (1σ) up to 7 keV. The energy scale will be measured on the ground through a dedicated calibration campaign using fiducial X-ray sources. Though calibrated, the energy scale is extremely sensitive to the environmental conditions around the TES array, and drifts in the readout chain electronics. Uncorrected, the energy scale can naturally drift up to hundreds of eVs. Changes of the TES gain will be monitored via onboard X-ray calibration sources, and the energy scale will be corrected either per pixel, or within a small groups of pixels. Although simulations show that a 0.4 eV level can be achieved, the very high accuracy required by the X-IFU calls for experimental validation. A dedicated measurement campaign has been performed by NASA Goddard Space Flight Center to characterize the energy scale of a prototype kilo-pixel array of X-IFU-representative TESs. The analysis of the data demonstrated the ability to correct for various drifts using two fiducial lines to track the temporal gain variation. In this paper, we propose to extend this study on the same data set by investigating multi-parameter correction techniques based on both the pulse-height of the fiducial line and the prepulse baseline level, using the knowledge of the TES energy scale at reference temperature/magnetic field set points acquired on the ground. Investigations on the co-adding of pixels to perform a joint correction over pools of pixels is also explored.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Origins of Extreme Velocity Stars as Revealed by Large-scale Galactic Surveys

We assemble a large sample of 12,784 high-velocity stars with total velocity V GSR ≥ 300 km s -1 , selected from RAVE DR5, SDSS DR12, LAMOST DR8, APOGEE DR16, GALAH DR2, and Gaia EDR3. In this sample, 52 are marginally hypervelocity star (HVS) candidates that have VGSR exceeding their local escape velocities within 2σ confidence levels, 40 of which are discovered for the first time. All of the candidates are metal-poor, late-type halo stars, which are significantly different from the previously identified HVSs, which are largely massive early-type stars, discovered by extreme radial velocity. This finding suggests that our newly identified HVS candidates are ejected by different mechanisms from the previous population. To investigate their origins, for 547 extreme velocity stars with V GSR ≥ 0.8V esc , we reconstruct their backward-integrated trajectories in the Galactic potential. According to the orbital analysis, no candidates are found to be definitely ejected from the Galactic-center (GC), while eight metal-poor extreme velocity stars are found to have a closest distance to the GC within 1 kpc. Intriguingly, 15 extreme velocity stars (including 2 HVS candidates) are found to have experienced close encounters with the Sagittarius dSph, suggesting that they originated from this dSph. This hypothesis is supported by an analysis of the [α/Fe]–[Fe/H] diagram. From a preliminary analysis of all of the 547 extreme velocity stars, we propose a general picture–star ejection from Galactic subsystems such as dwarf galaxies and globular clusters can be an important channel to produce extreme velocity stars or even HVSs, particularly the metal-poor late-type halo population.

79 ASTRONOMY AND ASTROPHYSICS↗

Local and Remote Mean and Extreme Temperature Response to Regional Aerosol Emissions Reductions

The climatic implications of regional aerosol and precursor emissions reductions implemented to protect human health are poorly understood. We investigate the mean and extreme temperature response to regional changes in aerosol emissions using three coupled chemistry–climate models: NOAA GFDL CM3, NCAR CESM1, and NASA GISS-E2. Our approach contrasts a long present-day control simulation from each model (up to 400 years with perpetual year 2000 or 2005 emissions) with 14 individual aerosol emissions perturbation simulations (160–240 years each). We perturb emissions of sulfur dioxide (SO2) and/or carbonaceous aerosol within six world regions and assess the statistical significance of mean and extreme temperature responses relative to internal variability determined by the control simulation and across the models. In all models, the global mean surface temperature response (perturbation minus control) to SO2 and/or carbonaceous aerosol is mostly positive (warming) and statistically significant and ranges from +0.17 K (Europe SO2) to -0.06 K (US BC). The warming response to SO2 reductions is strongest in the US and Europe perturbation simulations, both globally and regionally, with Arctic warming up to 1 K due to a removal of European anthropogenic SO2 emissions alone; however, even emissions from regions remote to the Arctic, such as SO2 from India, significantly warm the Arctic by up to 0.5 K. Arctic warming is the most robust response across each model and several aerosol emissions perturbations. The temperature response in the Northern Hemisphere midlatitudes is most sensitive to emissions perturbations within that region. In the tropics, however, the temperature response to emissions perturbations is roughly the same in magnitude as emissions perturbations either within or outside of the tropics. We find that climate sensitivity to regional aerosol perturbations ranges from 0.5 to 1.0 K (W m(exp -2))(exp -1) depending on the region and aerosol composition and is larger than the climate sensitivity to a doubling of CO2 in two of three models. We update previous estimates of regional temperature potential (RTP), a metric for estimating the regional temperature responses to a regional emissions perturbation that can facilitate assessment of climate impacts with integrated assessment models without requiring computationally demanding coupled climate model simulations. These calculations indicate a robust regional response to aerosol forcing within the Northern Hemisphere midlatitudes, regardless of where the aerosol forcing is located longitudinally. We show that regional aerosol perturbations can significantly increase extreme temperatures on the regional scale. Except in the Arctic in the summer, extreme temperature responses largely mirror mean temperature responses to regional aerosol perturbations through a shift of the temperature distributions and are mostly dominated by local rather than remote aerosol forcing.

Westervelt, Daniel M.↗

The Minos Computing Library: Efficient Parallel Programming for Extremely Heterogeneous Systems

Hardware specialization has become the silver bullet to achieve efficient high performance, from Systems-on-Chip systems, where hardware specialization can be ``extreme'', to large-scale HPC systems. As the complexity of the systems increases, so does the complexity of programming such architectures in a portable way. This work introduces the Minos Computing Library (MCL), as system software, programming model, and programming model runtime that facilitate programming extremely heterogeneous systems. MCL supports the execution of several multi-threaded applications within the same compute node, performs asynchronous execution of application tasks, efficiently balances computation across hardware resources, and provides performance portability. We show that code developed on a personal desktop automatically scales up to fully utilize powerful workstations with 8 GPUs and down to power-efficient embedded systems. MCL provides up to 17.5x speedup over OpenCL on NVIDIA DGX-1 systems and up to 1.88x speedup on single-GPU systems. In multi-application workloads, MCL dynamically resource allocation provides up to 2.43x performance improvement over manual, static allocation of computing resources.

Gioiosa, Roberto↗

Characterization of Sound Radiation by Unresolved Scales of Motion in Computational Aeroacoustics

Evaluation of the sound sources in a high Reynolds number turbulent flow requires time-accurate resolution of an extremely large number of scales of motion. Direct numerical simulations will therefore remain infeasible for the forseeable future: although current large eddy simulation methods can resolve the largest scales of motion accurately the, they must leave some scales of motion unresolved. A priori studies show that acoustic power can be underestimated significantly if the contribution of these unresolved scales is simply neglected. In this paper, the problem of evaluating the sound radiation properties of the unresolved, subgrid-scale motions is approached in the spirit of the simplest subgrid stress models: the unresolved velocity field is treated as isotropic turbulence with statistical descriptors, evaluated from the resolved field. The theory of isotropic turbulence is applied to derive formulas for the total power and the power spectral density of the sound radiated by a filtered velocity field. These quantities are compared with the corresponding quantities for the unfiltered field for a range of filter widths and Reynolds numbers.

Rubinstein, Robert↗

The Relationships Between the Trends of Mean and Extreme Precipitation

This study provides a better understanding of the relationships between the trends of mean and extreme precipitation in two observed precipitation data sets: the Climate Prediction Center Unified daily precipitation data set and the Global Precipitation Climatology Program (GPCP) pentad data set. The study employs three kinds of definitions of extreme precipitation: (1) percentile, (2) standard deviation and (3) generalize extreme value (GEV) distribution analysis for extreme events based on local statistics. Relationship between trends in the mean and extreme precipitation is identified with a novel metric, i.e. area aggregated matching ratio (AAMR) computed on regional and global scales. Generally, more (less) extreme events are likely to occur in regions with a positive (negative) mean trend. The match between the mean and extreme trends deteriorates for increasingly heavy precipitation events. The AAMR is higher in regions with negative mean trends than in regions with positive mean trends, suggesting a higher likelihood of severe dry events, compared with heavy rain events in a warming climate. AAMR is found to be higher in tropics and oceans than in the extratropics and land regions, reflecting a higher degree of randomness and more important dynamical rather than thermodynamical contributions of extreme events in the latter regions.

global warming↗

Multi-scale Multi-physics Scientific Machine Learning for Water Cycle Extreme Events Identification, Labelling, Representation, and Characterization

Impacts of climate are usually felt through extreme events such as droughts, floods, thunderstorms, windstorms, wildfires, and so on, that are intimately tied to the water cycle. Predicting the frequency and severity of extreme events under climate change remains a significant challenge; meanwhile, the mechanisms and impacts of these extremes are far from well understood. There are several major science challenges: (1) Lack of labelled extreme events data and missing standards in defining extremes; (2) Computational demand of high-resolution ensemble climate modeling; (3) Modeling the multiscale multi-physics hierarchical structure of compound extremes; (4) Lack of understanding of mechanisms of extreme events; (5) Large uncertainty in extreme events impacts on infrastructure; (6) Subjective assessment of weather-related risk from seasonal to multi-decadal time scales and lack of metrics for risk assessment and mitigation control.

54 ENVIRONMENTAL SCIENCES↗

Drought and Deluge: the Recurrence of Hydroclimate Extremes During the Past 600 Years in Eastern Australia’s Natural Resource Management (NRM) Clusters

Recent extremes of flood and drought across Australia have raised questions about the recurrence of such rare events and highlighted the importance of understanding multi-decadal climate variability. However, instrumental records over the past century are too short to adequately characterise climate variability on multi-decadal and longer timescales or robustly determine extreme event frequencies and their duration. Palaeoclimate reconstructions can provide much-needed information to help address this problem. Here, we use the 600-year hydroclimate record captured in the eastern Australian and New Zealand Drought Atlas (ANZDA) to analyse drought and pluvial frequency trends for East Australian Natural Resource Management (NRM) clusters. This partitioning of the drought atlas grid points into recognised biophysical areas (i.e. NRM clusters) revealed their differences and similarities in drought intensity and pluvial events over time. We find sustained multi-decadal periods of a wet–dry geographic 'seesaw' between eastern to central and southern NRMs (e.g. 1550–1600 CE and 1700–1750 CE). In contrast, other periods reveal spatially consistent wetting (e.g. 1500–1550 CE) or drying (e.g. 1750–1800 CE). Emerging hot spot analysis further shows that some areas that appear naturally buffered from severe drought during the instrumental period have a greater exposure risk when the longer 600-year record is considered. These findings are particularly relevant to management plans when dealing with the impacts of climate extremes developed at regional scales. Our results demonstrate that integrating and extending instrumental records with palaeoclimate datasets will become increasingly important for developing robust and locally specific extreme event frequency information for managing water resources.

Australian hydroclimate↗