Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Wealth inequality in the prehispanic northern US Southwest: from Malthus to Tyche

Persistent differences in wealth and power among prehispanic Pueblo societies are visible from the late AD 800s through the late 1200s, after which large portions of the northern US Southwest were depopulated. In this paper we measure these differences in wealth using Gini coefficients based on house size, and show that high Ginis (large wealth differences) are positively related to persistence in settlements and inversely related to an annual measure of the size of the unoccupied dry-farming niche. We argue that wealth inequality in this record is due first to processes inherent in village life which have internally different distributions of the most productive maize fields, exacerbated by the dynamics of systems of balanced reciprocity; and second to decreasing ability to escape village life owing to shrinking availability of unoccupied places within the maize dry-farming niche as villages get enmeshed in regional systems of tribute or taxation. We embed this analytical reconstruction in the model of an ‘Abrupt imposition of Malthusian equilibrium in a natural-fertility, agrarian society’ proposed by Puleston et al. (Puleston C, Tuljapurkar S, Winterhalder B. 2014 PLoS ONE 9, e87541 (doi:10.1371/journal.pone.0087541)), but show that the transition to Malthusian dynamics in this area is not abrupt but extends over centuries. This article is part of the theme issue ‘Evolutionary ecology of inequality’.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Model Scripts for "Old-Aged Groundwater Contributes to Mountain Hillslope Hydrologic Dynamics"

The partitioning of water inputs between deep and shallow groundwater flow paths is a fundamental processes, yet is challenging to observe. Numerical models provide a valuable tool to further develop insights on these groundwater mixing processes. This repository contains scripts to run the ParFlow-CLM and EcoSLIM integrated hydrologic models along the Pumphouse Hillslope in the East River Watershed and associated python scripts to process model outputs. Files includes input decks for the associated models, python scripts, and text and css files to support the model runs and processing. The 2-D hillslope model simulates from 2000 to 2021 using transient forcing conditions from the publicly available NLDAS-2 dataset. We develop an ensemble of models with variable hydrogeologic and soil parameters. Model outputs from the ensemble of runs are compared to to water level observations at the PLM1 and PLM6 wells published at (https://data.ess-dive.lbl.gov/view/doi:10.15485/1866836) and groundwater mean ages published at (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1960042. The model results are used to evaluate hillslope water mass-balance transience and mixing dynamics between groundwater flow paths with young (<10 year) and old (>10 year) ages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Facile Measurement of the Rotation of a Single Optically Trapped Nanoparticle Using the Diagonal Ratio of a Quadrant Photodiode

Optical tweezers are a powerful tool for exploring physical properties of particles in various environments through analysis of their dynamics in a trapping potential. Analysis of the trapped particle's position is often done using interferometric back-focal-plane detection microscopy in which interference between the trapping laser and the light interacting with the trapped particle is projected on a quadrant photo-diode (QPD). The QPD measurements are almost universally the voltage ratios between the transverse components of the detector (left/right; up/down) that are calibrated to ascertain particle positions and orientations (i.e., rotation). In this paper, we demonstrate how measuring the diagonal ratio of the QPD allows tracking the spinning and orbital motion of optically trapped objects with enhanced sensitivity compared to the conventional transverse ratio and reveals additional information about particle dynamics. The diagonal ratio "balances" the signal from opposite sides of the QPD, allowing common mode noise reduction and increased measurement sensitivity of the rotating nanospheres at frequencies from below 200 Hz to above 15 kHz. We show that the variability in particle rotation rates is due to slight differences in the aspect ratio of the nanoparticles. We also measure in situ a gradual acceleration in the spinning frequency of trapped nanoparticles due to photothermal shedding of surface ligands and the resulting decrease of rotational friction in solution. The accelerating particles spin at rates exceeding 30 kHz and are eventually "optically printed" onto the glass coverslip on top of the liquid cell enclosure. The diagonal ratio is a simple, readily available, and powerful enabler for measuring rotational dynamics with a 5-fold sensitivity increase over conventional QPD measurements. This affords opportunities for controlling the orientational dynamics of trapped nanoparticles, which is crucial for various applications, including nanoscale mechanical motors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗
DOE OSTI2025

Distinctive features of structural evolution and thermodynamic response in wide-bandgap semiconductors driven by intense electronic excitation

Radiation-tolerant material selection requires balancing lattice rigidity, defect dynamics, and electronic stability, as shown by covalent SiC outperforming ionic Ga 2 O 3 and GaN under extreme environments. Responding to intense electronic excitation, irradiation-driven phase segregation (β → δ/κ in Ga 2 O 3 ) and core–shell track (disordered structure in GaN), accompanied by elemental redistribution, contrastingly, exceptional radiation tolerance manifested by comparatively minimal lattice distortion (0.17 % strain variation) was demonstrated in SiC. These differential responses are primarily attributed to two fundamental mechanisms: (i) thermodynamic driving forces governing defect migration and phase separation, and (ii) the synergistic effects of robust covalent bonding composition coupled with efficient defect recombination processes. Here, the stronger electron–phonon (e-ph) coupling in Ga 2 O 3 (4.34 × 1018 W m −3 K −1 ) and GaN (3.55 × 10 18 W m −3 K −1 ) enhances lattice energy deposition, triggering thermal spikes (ΔT ≫ T m ) and structural transition behaviors, whereas weaker e-ph coupling in SiC (3.69 × 10 18 W m −3 K −1 ), relatively high thermodynamic parameters and efficient energy dissipation suppress thermal spikes to maintaining lattice integrity. The photoresponse degradation driven by enhanced radiative recombination is dominant in N-doped SiC, while V-doped systems achieve defect-mediated photoconduction optimization characterized by abrupt current transitions, matching fluorescence yield evolutions, and directly connecting defect engineering to optoelectronic performance.

Intense electronic excitation↗
DOE OSTI2020

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluating Methods for Measuring Grid Frequency in Low-Inertia Power Systems: Preprint

Accurate measurement of grid frequency is a critical component of reliable grid control. Traditionally, methods such as phase locked loops (PLLs) and discrete Fourier transforms (DFTs) have been used in inverters and phasor measurement units (PMUs) to measure frequency. However, as the percentage of inverter-based resources (IBRs) such as solar and wind has increased, these conventional frequency measurement methods are proving unable to guarantee reliable control in some cases. One challenge is measuring frequency during transient events, where there is a disruption in the steady state sinusoidal voltage. During these events, the underlying frequency of the grid may barely change, but measurement methods report a large spike in frequency due to the disrupted waveform. New methods must balance between suppressing spikes in frequency during faults, and providing fast, accurate, measurements in all other grid operation conditions, especially during events with high rate-of-change-of frequency (ROCOF), which are more prevalent in high-IBR power systems. This paper first surveys frequency measurement methods that have been proposed to reduce measurement errors during transient events. Then, both conventional and more novel frequency measurement methods are tested against an IEEE standard and industry recommendations, and their performance is evaluated for events simulated in PSCAD. Results quantify the trade-offs in performance during different grid conditions and lead to suggestions for the most appropriate frequency and ROCOF measurement methods for low inertia grids.

frequency↗
DOE OSTI2024

Tropical root responses to global changes: A synthesis

Tropical ecosystems face escalating global change. These shifts can disrupt tropical forests' carbon (C) balance and impact root dynamics. Since roots perform essential functions such as resource acquisition and tissue protection, root responses can inform about the strategies and vulnerabilities of ecosystems facing present and future global changes. However, root trait dynamics are poorly understood, especially in tropical ecosystems. We analyzed existing research on tropical root responses to key global change drivers: warming, drought, flooding, cyclones, nitrogen (N) deposition, elevated (e) CO 2 , and fires. Based on tree species- and community-level literature, we obtained 266 root trait observations from 93 studies across 24 tropical countries. We found differences in the proportion of root responsiveness to global change among different global change drivers but not among root categories. In particular, we observed that tropical root systems responded to warming and eCO 2 by increasing root biomass in species-scale studies. Drought increased the root: shoot ratio with no change in root biomass, indicating a decline in aboveground biomass. Despite N deposition being the most studied global change driver, it had some of the most variable effects on root characteristics, with few predictable responses. Episodic disturbances such as cyclones, fires, and flooding consistently resulted in a change in root trait expressions, with cyclones and fires increasing root production, potentially due to shifts in plant community and nutrient inputs, while flooding changed plant regulatory metabolisms due to low oxygen conditions. The data available to date clearly show that tropical forest root characteristics and dynamics are responding to global change, although in ways that are not always predictable. This synthesis indicates the need for replicated studies across root characteristics at species and community scales under different global change factors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling Growth Kinetics, Interspecies Cell Fusion, and Metabolism of a Clostridium acetobutylicum / Clostridium ljungdahlii Syntrophic Coculture

ABSTRACT Clostridium acetobutylicum and Clostridium ljungdahlii grown in a syntrophic culture were recently shown to fuse membranes and exchange cytosolic contents, yielding hybrid cells with significant shifts in gene expression and growth phenotypes. Here, we introduce a dynamic genome-scale metabolic modeling framework to explore how cell fusion alters the growth phenotype and panel of metabolites produced by this binary community. Computational results indicate C. ljungdahlii persists in the coculture through proteome exchange during fusing events, which endow C. ljungdahlii cells with expanded substrate utilization, and access to additional reducing equivalents from C. acetobutylicum -evolved H 2 and through acquisition of C. acetobutylicum -native cofactor-reducing enzymes. Simulations predict maximum theoretical ethanol and isopropanol yields that are increased by 0.64 mmol and 0.39 mmol per mmol hexose sugar consumed, respectively, during exponential growth when cell fusion is active. This modeling effort provides a mechanistic explanation for the metabolic outcome of cellular fusion and altered homeostasis achieved in this syntrophic clostridial community. IMPORTANCE Widespread cell fusion and protein exchange between microbial organisms as observed in synthetic C. acetobutylicum / C. ljungdahlii culture is a novel observation that has not been explored in silico . The mechanisms responsible for the observed cell fusion events in this culture are still unknown. In this work, we develop a modeling framework that captures the observed culture composition and metabolic phenotype, use it to offer a mechanistic explanation for how the culture achieves homeostasis, and identify C. ljungdahlii as primary beneficiary of fusion events. The implications for the events described in this study are far reaching, with potential to reshape our understanding of microbial community behavior synthetically and in nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Xyce(™) Parallel Electronic Simulator v.7.5

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗
DOE OSTI2021

yuwangcn/C4_dynamic_model

A dynamic systems model of C4 photosynthesis was developed based on the previous NADP-ME metabolic model for maize (Wang et al., 2014 ab). The NADP-ME metabolic model is an ordinary differential equation model including all individual steps in C4 photosynthetic carbon metabolism. Here, the model is extended to include posttranslational regulation and temperature response of enzyme activities, dynamic stomata conductance, and leaf energy balance. This model is written in Matlab (R2019a) Steady-state and dynamic gas exchange data for maize (B73), sugarcane (CP88-1762) and sorghum (Tx430) were measured in a greenhouse in Urbana, IL from July 25, 2019 through August 8, 2019. Data include: CO2 response curves Light response curves Photosynthetic induction curves measured in the transition from darkness to high light (1800 μmol m-2 s-1), data logged every 1 min. Photosynthetic induction in the transition from darkness to high light (1800 μmol m-2 s-1) to determine the kinetics of rubisco activation in these C4 crops (τ_Rubisco), data logged every 10 s. Gas exchange under fluctuating light. After dark adaptation, the leaves undergo three light change steps, light intensity was set as 1800 µmol m-2 s-1, 200 µmol m-2 s-1 and 1800 µmol m-2 s-1 for each 1800 s step.

Wang, Yu↗
DOE OSTI2026

Discerning Influences from Enthalpy and Entropy at Aqueous Interfaces Involved in Biomass Conversions in Porous Catalysts

Project Summary: The goal in this proposal is to learn how solvent influences the enthalpies and entropies of catalytic species in zeolite pores. Specifically, the physical, chemical, and structural features of solvent, catalytic species, and zeolite pores that determine enthalpies and entropies of solvation in solution phase biomass conversions will be interrogated using multiscale simulations and machine learning. Multiscale simulations are based off of existing strategies and employ quantum mechanics and classical molecular dynamics, providing an excellent balance between chemical accuracy and computational expense. They are capable of calculating enthalpies and entropies of solvation separately and have been validated in prior work to achieve high accuracy compared to their parent methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

diffReplication - An Energy-Aware Fault Tolerance Model for Silent Error Detection and Mitigation in Heterogeneous Extreme-scale Computing Environment

At extreme scale, the frequency of silent errors – a class of errors that remain undetected by low-level error detection mechanisms – increases significantly with the computational complexity of the application and the scale of the computing infrastructure. As hardware and software advances are made to usher in the next scientific era of computing, developing new approaches to mitigate the impact of silent errors remains a challenging problem. In this work, we propose an energy-aware fault-tolerance model, referred to diffReplication to overcome silent errors. In the proposed model, the main process is associated with one replica that executes at the same rate as the main process, and one diffReplica that is executed at a fraction of the main process' execution rate. If the main and its replica reach consensus at the end of a computation phase, the state of the diffReplica is updated and computation is resumed. If the synchronization attempt results in a disagreement, however, the diffReplica increases its execution speed to complete the computation and quickly reach the synchronization barrier. Assuming a single error over any given synchronization interval, a majority voting is used to reach consensus and tolerate silent errors. To further enhance its performance, diffReplication is augmented with speculative execution, whereby the main or its fast replica is selected to continue execution without waiting for the diffReplica. The selection process is based on the previous behaviour of the main and its replica. A performance analysis study is carried out to assess the performance of diffReplication, in terms of the energy saving and time-to-completion reduction achieved by the diffReplication scheme. The experiment shows that speculative execution reduces the time to completion with additional energy, and dynamic decision-making balances the energy consumption and time to completion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Insights into the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty

Abstract. The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Exploration of Defect Dynamics and Color Center Qubit Synthesis with Pulsed Ion Beams

Short-pulse ion beams have been developed in recent years and now enable applications in materials science. A tunable flux of selected ions delivered in pulses of a few nanoseconds can affect the balance of defect formation and dynamic annealing in materials. We report results from color center formation in silicon with pulses of 900 keV protons. G-centers in silicon are near-infrared photon emitters with emerging applications as single-photon sources and for spin-photon qubit integration. G-centers consist of a pair of substitutional carbon atoms and one silicon interstitial atom and are often formed by carbon ion implantation and thermal annealing. Here, we report on G-center formation with proton pulses in silicon samples that already contained carbon, without carbon ion implantation or thermal annealing. The number of G-centers formed per proton increased when we increased the pulse intensity from 6.9 × 10 9 to 7.9 × 10 10 protons/cm 2 /pulse, demonstrating a flux effect on G-center formation efficiency. We observe a G-center ensemble linewidth of 0.1 nm (full width half maximum), narrower than previously reported. Pulsed ion beams can extend the parameter range available for fundamental studies of radiation-induced defects and the formation of color centers for spin-photon qubit applications.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Reinforcement Learning for Load-balanced Parallel Particle Tracing

We explore an online reinforcement learning (RL) paradigm to dynamically optimize parallel particle tracing performance in distributed-memory systems. Our method combines three novel components: (1) a work donation algorithm, (2) a high-order workload estimation model, and (3) a communication cost model. First, we design an RL-based work donation algorithm. Our algorithm monitors workloads of processes and creates RL agents to donate data blocks and particles from high-workload processes to low-workload processes to minimize program execution time. The agents learn the donation strategy on the fly based on reward and cost functions designed to consider processes' workload changes and data transfer costs of donation actions. Second, we propose a workload estimation model, helping RL agents estimate the workload distribution of processes in future computations. Third, we design a communication cost model that considers both block and particle data exchange costs, helping RL agents make effective decisions with minimized communication costs. We demonstrate that our algorithm adapts to different flow behaviors in large-scale fluid dynamics, ocean, and weather simulation data. Our algorithm improves parallel particle tracing performance in terms of parallel efficiency, load balance, and costs of I/O and communication for evaluations with up to 16,384 processors.

Distributed and parallel particle tracing↗