SEARCH · Engineering Papers

Results for “sequential”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Data from a throughfall exclusion experiment: Fine root dynamics, morphology, chemistry, and AMF colonization across four lowland Panamanian forests

Fine roots regulate forest nutrient, carbon, and water cycling, yet their variation within and among tropical forests remains under-characterized. We quantified root productivity, disappearance, and stocks to 1 m using minirhizotron imaging, and we measured morphology, elemental composition [root carbon (C), root nitrogen (N), root phosphorus (P)], and arbuscular mycorrhizal fungi (AMF) colonization to 20 cm using ingrowth cores and sequential coring. Sampling took place in four distinct lowland Panamanian forests (32 plots; 8 per forest) from 2018 through 2022 under control and throughfall-exclusion (drought) treatments in the Panama Rainforest Changes with Experimental Drying (PARCHED) experiment.The dataset is presented as an Excel workbook with six tabs. The first tab is the data dictionary. Tab S1 contains ingrowth-core production and mortality, morphology and soil moisture. Tab S2 contains sequential-coring standing stocks with associated morphology and soil moisture. Tab S3 contains minirhizotron row data records to 1 m depth, including per-frame root length and diameter, normalized length metrics, and session timing. Tab S4 contains AMF colonization. Tab S5 contains fine-root chemistry at 0–10 cm, reporting %P, %C, %N, and C:N for samples collected via ingrowth cores and sequential-coring standing stocks. CSV mirrors for each tab are provided, and a KML file supplies coordinates for all 32 plots.Key variables span live and dead fine-root biomass (and coarse fractions where applicable), specific root length (SRL) and area (SRA), diameter, root tissue density (RTD), soil moisture, AMF colonization, root %N, %C, %P, and C:N, along with minirhizotron root length and diameter. Depth, season, treatment, and plot/site identifiers are included to support cross-tab integration and analysis from 0–100 cm (minirhizotron) and 0–20 cm (cores).Units are reported in-column and missing values are coded as NA. No special software is required to open or use the files (Excel, CSV, and KML compatible).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Control Co-Design Framework for Joint Optimization of PID Governor and Hydraulic Dynamics for a Hydropower Plant

This paper presents a control co-design (CCD) framework tailored to the hydraulic and mechanical subsystems of hydropower systems, aimed at enhancing their dynamic performance and responsiveness to grid demands. Conventional sequential design approaches often fall short in capturing the coupled interactions between plant dynamics and control objectives. The proposed CCD methodology enables simultaneous optimization of key hydro plant parameters - such as governor settings, penstock characteristics, and valve control - with the control system architecture. Results indicate that both CCD and sequential design approaches achieve comparable transient and steady-state performance, with minor discrepancies arising from differences in their optimized parameters. However, CCD demonstrates a distinct advantage by tripling the water time constant (Tw) while maintaining similar dynamic performance to the sequential approach. This increase in Tw enables optimization of penstock dimensions, as it is directly influenced by penstock design, thereby contributing to overall cost optimization.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗

DOE OSTI2021

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Iterative hybrid manufacture of a titanium alloy component

Here, this paper describes an iterative hybrid (additive + subtractive) manufacturing approach for a titanium alloy (Ti6Al4V) part using a laser hotwire directed energy deposition system (LHWDED) and a traditional four-axis milling machine tool. The term iterative hybrid manufacturing is used to described hybrid manufacturing where the additive and subtractive operations occur in multiple stages rather than sequentially. It is currently common to produce an entire part by sequential hybrid manufacturing by additively manufacturing (AM) an entire preform geometry that then requires post processing by another machine tool to create final part features. By contrast, a part produced by iterative hybrid manufacturing (IHM) does not produce the entire preform geometry in a single AM process. Instead, a portion of the entire preform geometry is manufactured by an AM process, then that portion is transferred to another machine tool which creates features in that portion, and then that machined portion is transferred back to the AM machine to complete another AM process. There is no limit to the number of iterations that an IHM process can have. IHM offers several advantages over sequential hybrid manufacturing such as the use of shorter and stiffer subtractive tooling, better access to part geometries that require subtractive processes, and the separation of the AM heat source from the subtractive machine tool. A titanium alloy demonstration part was successfully manufactured by IHM with three iterations using a shared pallet system between the AM machine tool and the subtractive machine tool.

Hybrid manufacturing↗

DOE OSTI2023

Commissioning of a two-target DC cylindrical magnetron sputter coater for depositing Nb 3 Sn film on Nb superconducting radiofrequency cavities

A DC cylindrical magnetron sputter coater was commissioned and used to coat Nb 2.6 GHz superconducting radiofrequency (SRF) cavity with Nb 3 Sn. The sputter coater has two identical cylindrical magnetrons that can move, with a controlled speed, along the axis of the SRF cavity to coat the inside surface of the cavity. The design of the sputter coater is discussed, along with its performance. Initially, a sample holder that allows coating on flat substrates at positions similar to the equator and beam tubes of a 2.6 GHz SRF cavity was used to test conditions for fabricating Nb 3 Sn layers. Multilayers of Nb and Sn were sequentially sputtered on flat Nb and sapphire substrates mounted on the equivalent positions of the cavity's beam tubes and the equator using the two identical cylindrical magnetrons. Then, the Nb/Sn multilayers were annealed at 950 °C for 3 h. Here, the ~1.2 μm thick Nb 3 Sn film did not show any other Nb–Sn compounds and had a superconducting transition temperature of 17.61–17.76 K. The 2.6 GHz SRF cavity was coated using similar conditions as flat samples. Cryogenic RF testing of the Nb 3 Sn-coated cavity demonstrated a quality factor of 3.2 × 10 8 at an accelerating gradient of 5 MV/m at 4.4 K.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Nonthermal Plasma-Stimulated C–N Coupling from CH 4 and N 2 Depends on the Presence of Surface CH x and Plasma-Phase CN Species

Formation of C–N containing compounds from plasma-catalytic coupling of CH 4 and N 2 over various transition metals (Ni, Pd, Cu, Ag, and Au) is investigated using a multimodal spectroscopic approach, combining polarization-modulation infrared reflection–absorption spectroscopy (PM-IRAS) and optical emission spectroscopy (OES). Through sequential experiments utilizing CH 4 and N 2 nonthermal plasmas, we minimize plasma-phase reactions and identify key intermediates for C–N coupling on metal surfaces. Results show that simultaneous CH 4 and N 2 exposure with plasma stimulation produces surface C–N species. However, N 2 –CH 4 sequential exposure does not lead to C–N species formation, while CH 4 –N 2 sequential exposure reveals the presence of CH x surface species and CN radical species as key precursors to C–N species formation. From further analysis using X-ray photoelectron spectroscopy and liquid chromatography–mass spectrometry, the influence of exposure conditions on the degree of nitrogen incorporation and the nature of C–N species formed were revealed. The work highlights the importance of surface chemistry and exposure conditions in surface C–N coupling with plasma stimulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Antimicrobial Properties of Corn Stover Lignin Fractions Derived from Catalytic Transfer Hydrogenolysis in Supercritical Ethanol with a Ru/C Catalyst

Converting lignin to value-added products at high yields provides an avenue for making ethanol biorefineries more profitable while reducing the carbon footprint of products generally derived from petroleum. In this study, corn stover lignin was depolymerized by catalytic transfer hydrogenolysis (CTH) in supercritical ethanol with a Ru/C catalyst. The lignin-derived bio-oil was then sequentially extracted utilizing hexane, petroleum ether, chloroform, and ethyl acetate as solvents in order of less polar to polar, and the subsequent bio-oils were characterized using GPC, GC/MS, and HSQC NMR. In this work, results show that the monomers in the bio-oil fractions contained primarily alkylated phenols, hydrogenated hydroxycinnamic acid derivatives, syringol and guaiacol-type lignins created from reductive cleavages of ether linkages, which were sequentially extracted into groups depending on the solvent polarity. The antimicrobial properties of the bio-oils were screened against Gram-positive (Bacillus subtilis, Lactobacillus amylovorus, and Staphylococcus epidermidis) and Gram-negative (Escherichia coli) bacteria and yeast (Saccharomyces cerevisiae) by examining microbial growth inhibition. Results show that CTH-derived bio-oils inhibited all tested organisms at concentrations less than 3 mg/mL. Total monomer concentration and the presence of specific monomers (i.e., syringyl propane) showed correlations to antimicrobial activity, likely due to cell death or membrane damage. This study provides insights into using sequential extraction to fractionate lignin-derived compounds and correlations between the properties of the extracted compounds and their antimicrobial activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Optimized active layer morphology toward efficient and polymer batch insensitive organic solar cells

Abstract Morphology control in laboratory and industry setting remains as a major challenge for organic solar cells (OSCs) due to the difference in film-drying kinetics between spin coating and the printing process. A two-step sequential deposition method is developed to control the active layer morphology. A conjugated polymer that self-assembles into a well-defined fibril structure is used as the first layer, and then a non-fullerene acceptor is introduced into the fibril mesh as the second layer to form an optimal morphology. A benefit of the combined fibril network morphology and non-fullerene acceptor properties was that a high efficiency of 16.5% (certified as 16.1%) was achieved. The preformed fibril network layer and the sequentially deposited non-fullerene acceptor form a robust morphology that is insensitive to the polymer batches, solving a notorious issue in OSCs. Such progress demonstrates that the utilization of polymer fibril networks in a sequential deposition process is a promising approach towards the fabrication of high-efficiency OSCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Control Co-Design Studies for a 22 MW Semisubmersible Floating Wind Turbine Platform

We present a control co-design software framework that can be used to optimize floating wind turbines and their controllers. Because this framework has many options for design variables, constraints, and merit figures, along with modeling fidelity levels, we seek to demonstrate best practices for using the tool while designing a floating platform for the new 22 MW offshore reference wind turbine developed within the International Energy Agency Wind Technology Commercialization Programme 55 on Reference Wind Turbines and Farms. During these studies, we evaluate the use of different simulation fidelity levels, the effect of using different load cases for controller tuning, and the difference between sequential and simultaneous control co-design solutions. Based on these efforts, we suggest using an algorithm that performs an initial search of the design space before optimization. We find that solving smaller optimization problems, in a sequential manner, leads to more reliable outcomes in fewer iterations than larger, simultaneous control co-design solutions. However a simultaneous CCD solution produces a platform with a 2% lower mass than the sequential CCD outcome.

17 WIND ENERGY↗

DOE OSTI2023

The Physical and Engineering Limits of Coupled In situ TEM Experiments

Almost all future energy systems (advanced nuclear reactors, fusion energy system, concentrating solar-thermal power (CSP), and wind technologies) are limited by degradation of key material systems exposed to multiple environmental stressors. The degradation during exposure to high temperatures, radiation, mechanical loading, and chemical attack is often dictated by mechanisms active at the microstructural level. The nature of these mechanisms and the associated variations between sequential and concurrent interplay can be explored if the transmission electron microscope (TEM) is utilized as a toolbox for exploration [1]. One such tool developed at Sandia National Laboratories to couple several of these environments is the In situ Ion Irradiation TEM (I3TEM) [2]. Several studies over the last decade utilizing this tool have shown that the scientific intuition developed over decades of sequential experiments is not always a good indicator of concurrent mechanisms or failure routes. This presentation will highlight the recent and planned additions into the I3TEM facility of a Waviks gas injection system and Raman system, respectively, as can be seen in Fig. 1, permitting environmental degradation from gas species leaked into the pole piece region during quantitative mechanical loading (indentation, monotonic loading, high temperature creep, irradiation induced creep, and high-cycle fatigue); multi-beam ion irradiation (energies ranging from 1 keV to 48 MeV and species from H to Au); laser exposure (20 W and 1064 nm), or various combinations thereof that are already possible [3]. This addition permits the facility to explore both sequential or concurrently the four axes of stressors: thermal, mechanical, radiation, and chemical. This information can be directly coupled to modeling, expediting the refinement and validation of both atomistic and mesoscale models.

degradation↗

DOE OSTI2021

Galerkin Neural Networks: A Framework for Approximating Variational Equations with Error Control

Herein, we present a new approach to using neural networks to approximate the solutions of variational equations, based on the adaptive construction of a sequence of finite-dimensional sub-spaces whose basis functions are realizations of a sequence of neural networks. Here, the finite-dimensional subspaces are then used to define a standard Galerkin approximation of the variational equation. This approach enjoys a number of advantages, including: the sequential nature of the algorithm offers a systematic approach to enhancing the accuracy of a given approximation; the sequential enhancements provide a useful indicator for the error that can be used as a criterion for terminating the sequential updates; the basic approach is largely oblivious to the nature of the partial differential equation under consideration; and, some basic theoretical results are presented regarding the convergence (or otherwise) of the method which are used to formulate basic guidelines for applying the method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Particle Filter Based Inference Testing

The primary intent of PAR-FIT (Particle Filter based Inference Testing) is to provide hard inductive evidence that a machine learning model is capable and proven for an individual test input. By examining training data used to form the underlying model functional correlation, an estimate of the reliability that a model will make the correct prediction can be made. The Sequential Probability Ratio Test is used to derive a qualitative evaluation for reliability based on hypothesis testing. The PAR-FIT framework achieves this by implementing a particle filter and the sequential probability ratio test algorithms on the machine learning model training data to determine relevancy of new individual test samples to the training dataset. The kernel function evaluates the local proximity and density of training data used to derive a prediction outcome. Particles are used to probabilistically determine which training data to evaluate for proximity. For test samples that are within a close proximity to and surrounded by multiple training data points, the evaluated reliability of the prediction is high. For test samples that are anomalies not represented by the training dataset, in low density data clusters, or are far from existing data points, the evaluated reliability is low as insufficient training evidence exists to suggest the model is capable of making the correct prediction. Sequential Probability Ratio Test is further used to determine when a hypothesis on whether a signal can be rejected or accepted for use. The ratio test collects sequence information from the particle filter to test whether the signal is anomalous or normal via hypothesis testing of the underlying distributions.

Chen, Edward [Idaho National Laboratory (INL), Ida↗

DOE OSTI2020

Parallel-in-Time Multigrid Methods for Hyperbolic Problems, with a Focus on the Shallow Water Equations

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations, which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic problems of programmatic interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic problems on future architectures. The difficulty lies in the fact that hyperbolic problems are well-known to be difficult for parallel-in-time methods, with the most common parallel-in-time method, Parareal, diverging in many cases. Here, the chosen methodology for scalably solving these hyperbolic space-time equations is multigrid-reductionin-time (MGRIT), because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. To further research in this area, this project will explore a model hyperbolic problem (the shallow water equations) in the context of recent advances (e.g., by Wingate and Haut) in constructing improved coarse time-grid time-propagators by using the slow asymptotic structure of the equations. In particular, we will investigate if such coarse time-propagators offer significant advantages in an MGRIT setting, and explore any broader insights gained into parallel-in-time for hyperbolic problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Parallel Multigrid in Time and Space for Extreme-Scale Computational Science

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Parallel Multigrid in Time and Space for Extreme-Scale Computational Science: Chaotic and Hyperbolic Problems

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallelin-time methods, with the most common method, parareal, diverging in many cases. The current stateof-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, we take inspiration from the Least Squares Shadowing (LSS by Wang) approach for solving chaotic problems. Here, an optimization approach is able to find “well-conditioned” shadow trajectories/solutions to the original “ill-conditioned” chaotic problem. Thus, the new multigrid methods researched here also arise in an optimization context.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Multigrid Reduction in Time for Chaotic and Hyperbolic Problems (Final Report)

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallelin-time methods, with the most common method, parareal, diverging in many cases. The current state-of-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, this work shall continue researching the directions from the current collaboration with Dr. Falgout, which are laid out in the work Toward Parallel in Time for Chaotic Dynamical Systems and showed the first known results of a parallel-in-time speedup for a chaotic problem. This work outlines two key improvements to XBraid for chaotic problems, the so-called “theta” and “delta-correction” methods. Here, these two improvements will be implemented in a high-performance but general way in XBraid and explored for more complicated problems. We will additionally research, as time allows, improvements to these techniques, as well as multigrid relaxation techniques based on Least Squares Shadowing (LSS by Wang) and a nonintrusive block tridiagonal solver based on MGRIT, called TriMGRIT.

97 MATHEMATICS AND COMPUTING↗