Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONSERVATION”

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
DOE OSTI2018

Comparison of Conservation Behavior in Several Variationally Derived Implementations for Simulating Electromagnetic Plasmas

Variational techniques based on Low's Lagrangian formulation have demonstrated significant improvement in conservation properties as compared to traditional particle-in-cell (PIC) algorithms. Previous work showed that equations derived from a discretized Lagrangian preserve a discretized version of the connection between Lagrangian symmetries and conservation laws, i.e., Noether's theorem. Specifically, a thorough analysis of energy conservation within the new system explained the absence of “grid heating” phenomena. In the present work, variationally derived implementations in a gridded periodic domain were compared to analogous implementations using Fourier bases. Since the Fourier basis allows for exact momentum conservation, this approach enables us to compare the relative merits of exact and approximate momentum conservation. Furthermore, to establish a tangible comparison for numerical methods, we investigate two simple physical examples; (1) the coupling of an electrostatic and electromagnetic plasma wave and (2) an equilibrium supporting Weibel instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Phylogenetic and amino acid conservation analyses of bacterial l-aspartate-α-decarboxylase and of its zymogen-maturation protein reveal a putative interaction domain

Background: All organisms must synthesize the enzymatic cofactor coenzyme A (CoA) from the precursor pantothenate. Most bacteria can synthesize pantothenate de novo by the condensation of pantoate and β-alanine. The synthesis of β-alanine is catalyzed by l-aspartate-α-decarboxylase (PanD), a pyruvoyl enzyme that is initially synthesized as a zymogen (pro-PanD). Active PanD is generated by self-cleavage of pro-PanD at Gly24-Ser25 creating the activesite pyruvoyl moiety. In Salmonella enterica, this cleavage requires PanM, an acetyl-CoA sensor related to the Gcn5-like N-acetyltransferases. PanM does not acetylate pro-PanD, but the recent publication of the three-dimensional crystal structure of the PanM homologue PanZ in complex with the PanD zymogen of Escherichia coli provides validation to our predictions and provides a framework in which to further examine the cleavage mechanism. In contrast, PanD from bacteria lacking PanM efficiently cleaved in the absence of PanM in vivo. Results: Using phylogenetic analyses combined with in vivo phenotypic investigations, we showed that two classes of bacterial l-aspartate-α-decarboxylases exist. This classification is based on their posttranslational activation by self-cleavage of its zymogen. Class I l-aspartate-α-decarboxylase zymogens require the acetyl-CoA sensor PanM to be cleaved into active PanD. This class is found exclusively in the Gammaproteobacteria. Class II l-aspartate-αdecarboxylase zymogens self cleave efficiently in the absence of PanM, and are found in a wide number of bacterial phyla. Several members of the Euryarchaeota and Crenarchaeota also contain Class II l-aspartate-α-decarboxylases. Phylogenetic and amino acid conservation analyses of PanM revealed a conserved region of PanM distinct from conserved regions found in related Gcn5-related acetyltransferase enzymes (Pfam00583). This conserved region represents a putative domain for interactions with l-aspartate-α-decarboxylase zymogens. This work may inform future biochemical and structural studies of pro-PanD-PanM interactions. Conclusions: Experimental results indicate that S. enterica and C. glutamicum l-aspartate-α-decarboxylases represent two different classes of homologues of these enzymes. Class I homologues require PanM for activation, while Class II self cleave in the absence of PanM. Computer modeling of conserved amino acids using structure coordinates of PanM and l-aspartate-α-decarboxylase available in the protein data bank (RCSB PDB) revealed a putative site of interactions, which may help generate models to help understand the molecular details of the self-cleavage mechanism of l-aspartate-α-decarboxylases.

Aspartate decarboxylase↗
DOE OSTI2024

Exploiting bacterial effector proteins to uncover evolutionarily conserved antiviral host machinery

Arboviruses are a diverse group of insect-transmitted pathogens that pose global public health challenges. Identifying evolutionarily conserved host factors that combat arbovirus replication in disparate eukaryotic hosts is important as they may tip the balance between productive and abortive viral replication, and thus determine virus host range. Here, we exploit naturally abortive arbovirus infections that we identified in lepidopteran cells and use bacterial effector proteins to uncover host factors restricting arbovirus replication. Bacterial effectors are proteins secreted by pathogenic bacteria into eukaryotic hosts cells that can inhibit antimicrobial defenses. Since bacteria and viruses can encounter common host defenses, we hypothesized that some bacterial effectors may inhibit host factors that restrict arbovirus replication in lepidopteran cells. Thus, we used bacterial effectors as molecular tools to identify host factors that restrict four distinct arboviruses in lepidopteran cells. By screening 210 effectors encoded by seven different bacterial pathogens, we identify several effectors that individually rescue the replication of all four arboviruses. We show that these effectors encode diverse enzymatic activities that are required to break arbovirus restriction. We further characterize Shigella flexneri-encoded IpaH4 as an E3 ubiquitin ligase that directly ubiquitinates two evolutionarily conserved proteins, SHOC2 and PSMC1, promoting their degradation in insect and human cells. We show that depletion of either SHOC2 or PSMC1 in insect or human cells promotes arbovirus replication, indicating that these are ancient virus restriction factors conserved across invertebrate and vertebrate hosts. Collectively, our study reveals a novel pathogen-guided approach to identify conserved antimicrobial machinery, new effector functions, and conserved roles for SHOC2 and PSMC1 in virus restriction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Genome-Wide Identification of Evolutionarily Conserved Alternative Splicing Events in Flowering Plants

Alternative splicing (AS) plays important roles in many plant functions, but its conservation across the plant kingdom is not known. We describe a methodology to identify AS events and identify conserved AS events across large phylogenetic distances using RNASeq datasets. We applied this methodology to transcriptome data from nine angiosperms including Amborella, the single sister species to all other extant flowering plants. AS events within 40–70% of the expressed multi-exonic genes per species were found, 27,120 of which are conserved among two or more of the taxa studied. While many events are species specific, many others are shared across long evolutionary distances suggesting they have functional significance. Conservation of AS event data provides an estimate of the number of ancestral AS events present at each node of the tree representing the nine species studied. Furthermore, the presence or absence of AS isoforms between species with different whole genome duplication (WGD) histories provides the opportunity to examine the impact of WDG on AS potential. Examining AS in gene families identifies those with high rates of AS, and conservation can distinguish ancient events vs. recent or species specific adaptations. The MADS-box and SR protein families are found to represent families with low and high occurrences of AS, respectively, yet their AS events were likely present in the MRCA of angiosperms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fluctuations and Conservation Laws

We report on recent progress concerning effects of global conservation laws on cumulants of conserved quantities. Specifically, we will relate for an arbitrary equation of state cumulants of a conserved charge measured in a subvolume of a thermal system with the corresponding grandcanonical susceptibilities, taking into account exact global conservation of that charge. Applications to actual measurement at the RHIC and LHC as well as extensions to multiple conserved charges will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Draft Environmental Impact Statement for Proposed Energy Conservation Standards for Manufactured Housing

This draft environmental impact statement (EIS) analyzes the impacts related to DOE’s proposed energy conservation standards for manufactured homes. DOE’s proposed energy conservation standards are based on the 2021 International Energy Conservation Code (IECC) and relate to the building thermal envelope; air sealing; installation of insulation; duct sealing; heating, ventilation and air conditioning (HVAC); service hot water systems; mechanical ventilation fan efficacy; and heating and cooling equipment sizing. The draft EIS evaluates DOE’s proposed action (energy conservation standards tiered based the manufacturers retail list price), alternatives to the proposed action (including energy conservation standards tiered based on size of the manufactured home and untiered standards), and a no action alternative. The draft EIS presents analysis of impacts related to (1) energy resources; (2) air resources (including greenhouse gases, climate change, ambient air quality and indoor air quality); (3) health and safety; (4) socioeconomics; and (5) environmental justice.

MacDonell, Margaret↗
DOE OSTI2026

Entropy Stable Conservative Flux Form Neural Networks

We propose an entropy-stable conservative flux form neural network (CFN) to predict the dynamics of unknown governing conservation laws. The design of the network is based on the entropy-stable, second-order, and non-oscillatory Kurganov-Tadmor (KT) scheme. The proposed entropy-stable CFN, hereafter referred to as ESCFN, uses slope limiting as a denoising mechanism, ensuring accurate predictions in both noisy and sparse observation environments, as well as in both smooth and discontinuous regions. Importantly, our method is designed to predict long term dynamics of the unknown conservation law exclusively from a short temporal window of observed data, that is, without oracle knowledge of the PDE or later-time solution profiles. Numerical experiments demonstrate that the ESCFN achieves both stability and conservation while maintaining accuracy over extended time domains, and successfully predicts shock propagation speeds in long-term simulations. Furthermore, it is also robust to both noisy and sparse data environments.

Hyperbolic conservation laws↗
DOE OSTI2019

Conservation agriculture based on diversified and high-performance production system leads to soil carbon sequestration in subtropical environments

Soils can be a source or sink of atmospheric CO 2 , depending on the historic and existing land use and management. We used long term soil management database of a production farm that is based on principles of conservation agriculture such as: a) eliminate soil disturbance; b) maintain permanent soil surface cover; C) adopt crop diversity with high biomass-C input; for 30 years and agroecosystem models to study the potential of different management options to sequester C in soils. Using Century and Roth-C models we simulated the carbon stocks evolution in the farm and four subtropical soil management scenarios and studied C sequestration potential. The scenarios were: a) existing farm biomass input (14.5 Mg ha -1 year -1 ) or C input (6.5 Mg ha -1 year -1 ); b) 15% increase of farm biomass input (16.7 Mg ha -1 year -1 ) or C input (7.5 Mg ha -1 year -1 ); c) 15% decrease of farm biomass (12.3 Mg ha -1 year -1 ) or C input (5.5 Mg ha -1 year -1 ) and, d) 30% decrease of farm biomass input (10.1 Mg ha -1 year -1 ) or C input (4.5 Mg ha -1 year -1 ). Our results demonstrate that soil organic carbon continuously increased after conservation management practices adoption in 1985 until 2015, and currently soil organic carbon is in equilibrium. We found that an increase of 2.2 Mg ha -1 year -1 biomass-C input for 60 years resulted into increase of 12 Mg ha -1 soil organic carbon stocks. The same way, crop yields increased with time, and were more pronounced for maize compared to soybean and wheat. The scaling up of model results to similar climate and soil types indicated that conservation management practices has the potential to sequester 2.7 ± 0.02 Pg C at 0–20 cm and 4.8 ± 3 Pg C at 0–100 cm soil depth in 43 million ha area globally. In the 30% and 15% decrease scenarios the sequestration were 2.2 ± 0.02 and 2.4 ± 0.02 Pg C at 0–20 cm an in 15% increase scenario it goes to 3.2 ± 0.02 Pg C. Here, this equilibrium soil organic carbon stocks considering the currently adopted system are equivalent to 3.5–4.5% of the world SOC stocks in 3% of the world croplands and correspond to 6 years of global land use and land use change emissions, indicating that conservation management practices can lead the soil be a sink and a promising tool to promote C sequestration in subtropical soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Energy-conserving perfect-conductor boundary conditions for an implicit, curvilinear Darwin particle-in-cell algorithm

We propose here a conservative prescription for perfect-conductor boundary conditions for a Darwin particle-in-cell (PIC) model on curvilinear meshes and arbitrarily shaped boundaries. The Darwin model is a subset of Maxwell's equation in which the electromagnetic mode (light wave) has been analytically removed. This renders the Darwin formulation elliptic, instead of hyperbolic. Historically, Darwin-PIC practitioners have had difficulty implementing a well-posed set of boundary conditions for realistic applications. In this study, we demonstrate the well-posedness and effectiveness of a simple boundary-condition prescription for perfect conductors of specified potential and of arbitrary shape in multiple dimensions, using a recently developed conservative, implicit Darwin-PIC algorithm. The boundary conditions conserve the electromagnetic energy, and preserve the Coulomb gauge of the vector potential, ∇ · A = 0 exactly. We demonstrate that global energy is exactly conserved in curvilinear, simply-connected domains in the specific case of particles reflecting off perfect conductors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Using Change Models to Envision Better Applications of Animal Behavior Research in Conservation Management and Beyond

While most animal behavior researchers have mastered the process of knowledge creation, generating knowledge that can readily be applied requires a different set of skills. The process and timeframe of fundamental scientific knowledge production is often not relevant to those who might apply it, such as conservation or wildlife managers. Additionally, the complex challenges that policy makers, managers and practitioners face are often not adequately communicated to and among scientists. This mutual disconnect in discourse, relationships, common terms, and practices is especially apparent when animal behavior researchers seek to have applied impact. We argue that bridging the complex implementation gap in animal behavior requires a formalized vision for change. We turn to change model theory, a tool commonly used in other fields for identifying the links between actions and outcomes necessary for enacting large-scale change. We focus on the subfield of conservation behavior with a change model that outlines specific ways to improve collaboration and coordination between animal behavior science and conservation practice. We present this targeted change model, review each strategy the model outlines, and highlight pressing actions that people from various career stages and backgrounds can take. We encourage researchers to further the alignment of science with management needs by developing the proper communication mechanisms for improved cultural exchange and plan future change model efforts directly targeting managers. Beyond the conservation behavior change model we present, we also discuss the broad applicability of change models to enhance the application of academic research to other fields. Fundamental science researchers are increasingly required to show impact of their work on society; the change model process we describe here can enable further impact.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Path-Integrated Stereo X-Ray Digital Image Correlation: Resolving the Violation of Conservation of Intensity

In this study, x-ray imaging addresses many challenges with visible light imaging in extreme environments, where optical diagnostics such as digital image correlation (DIC) and particle image velocimetry (PIV) suffer biases from index of refraction changes and/or cannot penetrate visibly occluded objects. However, conservation of intensity—the fundamental principle of optical image correlation algorithms—may be violated if ancillary components in the X-ray path besides the surface or fluid of interest move during the test. The test series treated in this work sought to characterize the safe use of fiber-epoxy composites in aerospace and aviation industries during fire accident scenarios. Stereo X-ray DIC was employed to measure test unit deformation in an extreme thermal environment involving a visibly occluded test unit, incident surface heating to temperatures above 600°C, and flames and soot from combusting epoxy decomposition gasses. The objective of the current work is to evaluate two solutions to resolve the violation of conservation of intensity that resulted from both the test unit and the thermal chamber deforming during the test. The first solution recovered conservation of intensity by pre-processing the path-integrated X-ray images to isolate the DIC pattern of the test unit from the thermal chamber components. These images were then correlated with standard, optical DIC software. The second solution, called Path-Integrated Digital Image Correlation (PI-DIC), modified the fundamental matching criterion of DIC to account for multiple, independently-moving components contributing to the final image intensity. The PI-DIC algorithm was extended from a 2D framework to a stereo framework and implemented in a custom DIC software. Both solutions accurately measured the cylindrical shape of the undeformed test unit, recovering radii values of R = 76.20±0.12 mm compared to the theoretical radius of R theor = 76.20 mm. During the test, the test unit bulged asymmetrically as decomposition gasses pressurized the interior and eventually burned in a localized jet. Both solutions measured the heterogeneous radius of this bulge, which reached a maximum of approximately R = 91 mm, with a discrepancy of 2–3% between the two solutions. Two solutions that resolve the violation of conservation of intensity for path-integrated X-ray images were developed. Both were successfully employed in a stereo X-ray DIC configuration to measure deformation of an aluminum-skinned, fiber-epoxy composite test unit in a fire accident scenario.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

High-order, stable, and conservative boundary schemes for central and compact finite differences

Stable and conservative numerical boundary schemes are constructed such that they do not diminish the overall accuracy of the method for interior schemes of orders 4, 6, and 8 using both explicit (central) and compact finite differences. Previous attempts to develop stable numerical boundary schemes for non-linear problems have resulted in schemes which significantly reduced the global accuracy and/or required some form of artificial dissipation. Thus, the schemes developed in this paper are the first to not require this tradeoff, while also ensuring discrete conservation and allowing for direct boundary condition enforcement. After outlining a general procedure for the construction of conservative boundary schemes of any order, a simple, yet novel, optimization strategy which focuses directly on the compressible Euler equations is presented. Furthermore, the result of this non-linear optimization process is a set of high-order, stable, and conservative numerical boundary schemes which demonstrate excellent stability and convergence properties on an array of linear and non-linear hyperbolic problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A Legendre–Fourier spectral method with exact conservation laws for the Vlasov–Poisson system

In this study, we present the design and implementation of an L 2 -stable spectral method for the discretization of the Vlasov–Poisson model of a collisionless plasma in one space and velocity dimension. The velocity and space dependence of the Vlasov equation are resolved through a truncated spectral expansion based on Legendre and Fourier basis functions, respectively. The Poisson equation, which is coupled to the Vlasov equation, is also resolved through a Fourier expansion. The resulting system of ordinary differential equation is discretized by the implicit second-order accurate Crank–Nicolson time discretization. The non-linear dependence between the Vlasov and Poisson equations is iteratively solved at any time cycle by a Jacobian-Free Newton–Krylov method. In this work we analyze the structure of the main conservation laws of the resulting Legendre–Fourier model, e.g., mass, momentum, and energy, and prove that they are exactly satisfied in the semi-discrete and discrete setting. The L 2 -stability of the method is ensured by discretizing the boundary conditions of the distribution function at the boundaries of the velocity domain by a suitable penalty term. The impact of the penalty term on the conservation properties is investigated theoretically and numerically. An implementation of the penalty term that does not affect the conservation of mass, momentum and energy, is also proposed and studied. A collisional term is introduced in the discrete model to control the filamentation effect, but does not affect the conservation properties of the system. Numerical results on a set of standard test problems illustrate the performance of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm

The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. A novel, implicit, particle-in-cell based scheme for the hybrid model is derived for fully 3D electromagnetic problems with multiple ion species, which features global mass, momentum and energy conservation. The scheme includes sub-cycling and orbit-averaging for the ions, with cell-centered finite differences and implicit midpoint time advance. To reduce discrete particle noise, the scheme allows arbitrary-order shape functions for the particle-mesh interpolations and the application of conservative binomial smoothing. The algorithm is verified for a number of test problems to demonstrate the correctness of the implementation, the unique conservation properties, and the favorable stability properties of the new scheme. In particular, there is no indication of unstable growth of the finite-grid instability for a population of cold ions drifting through a uniform spatial mesh, in a set-up where several commonly used non-conservative schemes are highly unstable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High order explicit local time stepping methods for hyperbolic conservation laws

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RK-DG) methods and design LTS algorithms based on the strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integration in different regions of the computational domain. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of the interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

58 GEOSCIENCES↗
DOE OSTI2022

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗