Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow simulation”

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 1,081 records · Page 60

Development of a reduced tri-propylene glycol monomethyl ether– n -hexadecane–poly-aromatic hydrocarbon mechanism and its application for soot prediction

A reduced chemical kinetic mechanism for Tri-Propylene Glycol Monomethyl Ether (TPGME) has been developed and applied to computational fluid dynamics (CFD) calculations for predicting combustion and soot formation processes. The reduced TPGME mechanism was combined with a reduced n-hexadecane mechanism and a Poly-Aromatic Hydrocarbon (PAH) mechanism to investigate the effect of fuel oxygenation on combustion and soot emissions. The final version of the TPGME-n-hexadecane-PAH mechanism consists of 144 species and 730 reactions and was validated with experiments in shock tubes as well as in a constant volume spray combustion vessel (CVCV) from the Engine Combustion Network (ECN). The effects of ambient temperature, varying oxygen content in the tested fuels on ignition delay, spray liftoff length and soot formation under diesel-like conditions were analyzed and addressed using multidimensional reacting flow simulations and the reduced mechanism. Here, the results show that the present reduced mechanism gives reliable predictions of the combustion characteristics and soot formation processes. In the CVCV simulations, two important trends were identified. First, increasing the initial temperature in the CVCV shortens the ignition delay and lift-off length, reduces the fuel-air mixing, thereby increasing the soot levels. Secondly, fuel oxygenation introduces more oxygen into the central region of a fuel jet and reduces residence times of fuel rich area in active soot forming regions, thereby reducing soot levels.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modeling the potential impact on the US blood supply of transfusing critically ill patients with fresher stored red blood cells

Although some studies have suggested that transfusion recipients may have better medical outcomes if transfused with red blood cell units stored for a short time, the overall body of evidence shows mixed results. It is important to understand how using fresher stored red blood cell units for certain patient groups may affect blood availability. Based on the Stock-and-Flow simulation model of the US blood supply developed by Simonetti et al. 2014, we evaluated a newly implemented allocation method of preferentially transfusing fresher stored red blood cell units to a subset of high-risk group of critically ill patients and its potential impact on supply. Simulation results showed that, depending on the scenario, the US blood total supply might be reduced between 2-42%, when compared to the standard of care in transfusion medicine practice. Among our simulated scenarios, we observed that the number of expired red blood cell units modulated the supply levels. The age threshold of the required red blood cell units was inversely correlated with both the supply levels and the number of transfused units that failed to meet that age threshold. Conclusion To our knowledge, this study represents the first attempt to develop a comprehensive framework to evaluate the impact of preferentially transfusing fresher stored red blood cells to the higher-risk critically ill patients on supply. Model results show the difficulties to identify an optimal scenario.

59 BASIC BIOLOGICAL SCIENCES↗

Streamline-Based Time-Lapse-Seismic-Data Integration Incorporating Pressure and Saturation Effects

In this paper, we present an efficient history-matching technique that simultaneously integrates 4D repeat seismic surveys with well-production data. This approach is particularly well-suited for the calibration of the reservoir properties of high-resolution geologic models because the seismic data are a really dense but sparse in time, whereas the production data are finely sampled in time but spatially averaged. The joint history matching is performed by use of streamline-based sensitivities derived from either finite-difference or streamline-based flow simulation. For the most part, earlier approaches have focused on the role of saturation changes, but the effects of pressure have largely been ignored. Here, we present a streamline-based semianalytic approach for computing model-parameter sensitivities, accounting for both pressure and saturation effects. The novelty of the method lies in the semianalytic sensitivity computations, making it computationally efficient for high-resolution geologic models. The approach is implemented by use of a finite-difference simulator incorporating the detailed physics. Its efficacy is demonstrated by use of both synthetic and field applications. Lastly, for both the synthetic and the field cases, the advantages of incorporating the time-lapse variations are clear, seen through the improved estimation of the permeability distribution, the pressure profile, the evolution of the fluid saturation, and the swept volumes.

58 GEOSCIENCES↗

Relationship between Anisotropy in Soil Hydraulic Conductivity and Saturation

Anisotropy in unsaturated hydraulic conductivity is saturation-dependent. Accurate characterization of soil anisotropy is very important in simulating flow and contaminant (e.g., radioactive nuclides in Hanford) transport. A recently developed tensorial connectivity-tortuosity (TCT) concept describes the hydraulic conductivity tensor of the unsaturated anisotropic soils as the product of a scalar variable, the symmetric connectivity tortuosity tensor, and the hydraulic conductivity tensor at saturation. In this study, the TCT model is used to quantify soil anisotropy in unsaturated hydraulic conductivity. The TCT model can describe different types of soil anisotropy; e.g., the anisotropy coefficient, C, can be monotonically increase or decrease with saturation and can vary from greater than unity to less than unity and vice versa. Soil anisotropy is independent of soil water retention properties and can be characterized by the ratio of the saturated hydraulic conductivities and the difference of the tortuosity-connectivity coefficients in two directions. ln(C) is linearly proportional to ln(Se) with Se being the effective saturation. The log-linear relationship between C and Se allows the saturation-dependent anisotropy to be determined using linear regression with the measurements of the directional hydraulic conductivities at a minimum of two water content levels, of which one may be at full saturation. The model was tested using measurements of directional hydraulic conductivities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A High-Performance Embedded Hybrid Methodology for Uncertainty Quantification With Applications

Multiphysics processes modeled by a system of unsteady differential equations are naturally suited for partitioned (modular) solution strategies. We consider such a model where probabilistic uncertainties are present in each module of the system and represented as a set of random input parameters. A straightforward approach in quantifying uncertainties in the predicted solution would be to sample all the input parameters into a single set, and treat the full system as a black-box. Although this method is easily parallelizable and requires minimal modifications to deterministic solver, it is blind to the modular structure of the underlying multiphysical model. On the other hand, using spectral representations polynomial chaos expansions (PCE) can provide richer structural information regarding the dynamics of these uncertainties as they propagate from the inputs to the predicted output, but can be prohibitively expensive to implement in the high-dimensional global space of uncertain parameters. Therefore, we investigated hybrid methodologies wherein each module has the flexibility of using sampling or PCE based methods of capturing local uncertainties while maintaining accuracy in the global uncertainty analysis. For the latter case, we use a conditional PCE model which mitigates the curse of dimension associated with intrusive Galerkin or semi-intrusive Pseudospectral methods. After formalizing the theoretical framework, we demonstrate our proposed method using a numerical viscous flow simulation and benchmark the performance against a solely Monte-Carlo method and solely spectral method.

97 MATHEMATICS AND COMPUTING↗

An evaluation of the carbon sequestration potential of the Cambro-Ordovician Strata of the Illinois and Michigan basins

The studies summarized herein were conducted during 2009–2014 to investigate the utility of the Knox Group and St. Peter Sandstone deeply buried geologic strata for underground storage of carbon dioxide (CO 2 ), a practice called CO 2 sequestration (CCS). In the subsurface of the midwestern United States, the Knox and associated strata extend continuously over an area approaching 500,000 sq. km, about three times as large as the State of Illinois. Although parts of this region are underlain by the deeper Mt. Simon Sandstone, which has been proven by other Department of Energy-funded research as a resource for CCS, the Knox strata may be an additional CCS resource for some parts of the Midwest and may be the sole geologic storage (GS) resource for other parts. One group of studies assembles, analyzes, and presents regional-scale and point-scale geologic information that bears on the suitability of the geologic formations of the Knox for a CCS project. New geologic and geo-engineering information was developed through a small-scale test of CO 2 injection into a part of the Knox, conducted in western Kentucky. These studies and tests establish the expectation that, at least in some locations, geologic formations within the Knox will (a) accept a commercial-scale flow rate of CO 2 injected through a drilled well; (b) hold a commercial-scale mass of CO 2 (at least 30 million tons) that is injected over decades; and (c) seal the injected CO 2 within the injection formations for hundreds to thousands of years. In CCS literature, these three key CCS-related attributes are called injectivity, capacity, and containment. The regional-scale studies show that reservoir and seal properties adequate for commercial-scale CCS in a Knox reservoir are likely to extend generally throughout the Illinois and Michigan Basins. Information distinguishing less prospective subregions from more prospective fairways is included in this report. Another group of studies report the results of reservoir flow simulations that estimate the progress and outcomes of hypothetical CCS projects carried out within the Knox (particularly within the Potosi Dolomite subunit, which, in places, is highly permeable) and within the overlying St. Peter Sandstone. In these studies, the regional-scale information and a limited amount of detailed data from specific boreholes is used as the basis for modeling the CO 2 injection process (dynamic modeling). The simulation studies were conducted progressively, with each successive study designed to refine the conclusions of the preceding one or to answer additional questions. The simulation studies conclude that at Decatur, Illinois or a geologically similar site, the Potosi Dolomite reservoir may provide adequate injectivity and capacity for commercial-scale injection through a single injection well. This conclusion depends on inferences from seismic-data attributes that certain highly permeable horizons observed in the wells represent laterally persistent, porous vuggy zones that are vertically more common than initially evident from wellbore data. Lateral persistence of vuggy zones is supported by isotopic evidence that the conditions that caused vug development (near-surface processes) were of regional rather than local scale. Other studies address aspects of executing and managing a CCS project that targets a Knox reservoir. These studies cover well drilling, public interactions, representation of datasets and conclusions using geographic information system (GIS) platforms, and risk management.

54 ENVIRONMENTAL SCIENCES↗

Gulf of Mexico miocene CO₂ site characterization mega transect

This project characterized the Miocene-age sub-seafloor stratigraphy in the near-offshore portion of the Gulf of Mexico adjacent to the Texas coast. The large number of industrial sources of carbon dioxide (CO₂) in coastal counties and the high density of onshore urbanization and environmentally sensitive areas make this offshore region extremely attractive for long-term storage of carbon dioxide emissions from industrial sources (CCS). The study leverages dense existing geologic data from decades of hydrocarbon exploration in and around the study area to characterize the regional geology for suitability and storage capacity. Primary products of the study include: regional static storage capacity estimates, sequestration “leads” and prospects with associated dynamic capacity estimates, experimental studies of CO₂-brine-rock interaction, best practices for site characterization, a large-format ‘Atlas’ of sequestration for the study area, and characterization of potential fluid migration pathways for reducing storage risks utilizing novel high-resolution 3D (HR3D) seismic surveys. In addition, three subcontracted studies address source-to-sink matching optimization, offshore well bore management and environmental aspects. The various geologic data and interpretations are integrated and summarized in a series of cross-sections and maps, which represent a primary resource for any near-term commercial deployment of CCS in the area. The regional study characterized and mapped important geologic features (e.g., Clemente-Tomas fault zone, the regionally extensive Marginulina A and Amphistegina B confining systems, etc.) that provided an important context for regional static capacity estimates and specific sequestration prospects of the study. A static capacity estimate of the majority of the Study area (14,467 mi 2 ) was estimated at 86 metric Gigatonnes. While local capacity estimates are likely to be lower due to reservoir-scale characteristics, the offshore Miocene interval is a storage resource of National interest for providing CO₂ storage as an atmospheric emissions abatement strategy. The natural petroleum system was used as an analog to infer seal quality and predict possible migration pathways of fluids in an engineered system of anthropogenic CO₂ injection and storage. The regional structural features (e.g., Clemente-Tomas fault zone) that exert primary control on the trapping and distribution of Miocene hydrocarbons are expected to perform similarly for CCS. Industrial-scale CCS will require storage capacity utilizing well-documented Miocene hydrocarbon (dominantly depleted gas) fields and their larger structural closures, as well as barren (unproductive, brine-filled) closures. No assessment was made of potential for CO₂ utilization for enhanced oil and gas recovery. The use of 3D numerical fluid flow simulations have been used in the study to greatly assist in characterizing the potential storage capacity of a specific reservoir. Due to the complexity of geologic systems (stratigraphic heterogeneity) and inherent limitations on producing a 3D geologic model, these simulations are typically simplified scenarios that explore the influence of model property variability (sensitivity study). A specific site offshore San Luis Pass (southern Galveston Island) was undertaken successfully, indicating stacked storage potential. Downscaling regional capacity estimates to the local scale (and the inverse) has proven challenging, and remains an outstanding gap in capacity assessments. In order to characterize regional seal performance and identify potential brine and CO₂ leakage pathways, results from three high-resolution 3D (HR3D) seismic datasets acquired by the study using novel HR3D (P-Cable) acquisition system showed steady and significant improvements in data quality because of improved acquisition and processing technique. Finely detailed faults and stratigraphy in the shallowest 1000 milliseconds (~800 m) of data allowed for the identification and mapping of unconformable surfaces including what is probably a surface associated with the last Pleistocene glacial lowstand. The identification of a previously unrecognized (in commercial seismic data) gas chimney that was clearly defined in the 2013 HR3D survey, indicates that HR3D surveys may be useful as both a characterization tool for the overburden of a potential carbon sequestration site and as an additional monitoring tool for future engineered injection sites. Geochemical modeling indicated that injection of CO₂ would result in minor dissolution of calcite, K-feldspar and albite. In addition, modeling of typical brines in Miocene age rocks indicate that approximately 5% of injection capacity would result from CO₂ dissolution into the brine. After extensive searches, no rock samples of the Marginulina A and Amphistegina B seals (“caprocks”) were obtained, but analyses of available core samples of other Miocene age mudrocks (seals or caprocks) indicate that they have sealing ability sufficient for potential CO 2 storage in underlying sandstone units.

01 COAL, LIGNITE, AND PEAT↗

sRecovery Act: Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO 2 Storage

This study focuses on evaluating the feasibility and suitability of using the Jurassic/Triassic (J/TR) sediments of the South Georgia Rift basin (SGR) for CO 2 storage in southern South Carolina and southern Georgia The SGR basin in South Carolina (SC), prior to this project, was one of the least understood rift basin along the east coast of the U.S. In the SC part of the basin there was only one well (Norris Lightsey #1) the penetrated into J/TR. Because of the scarcity of data, a scaled approach used to evaluate the feasibility of storing CO 2 in the SGR basin. In the SGR basin, 240 km (~149 mi) of 2-D seismic and 2.6 km2 3-D (1 mi2) seismic data was collected, process, and interpreted in SC. In southern Georgia 81.3 km (~50.5 mi) consisting of two 2-D seismic lines were acquired, process, and interpreted. Seismic analysis revealed that the SGR basin in SC has had a very complex structural history resulting the J/TR section being highly faulted. The seismic data is southern Georgia suggest SGR basin has not gone through a complex structural history as the study area in SC. The project drilled one characterization borehole (Rizer # 1) in SC. The Rizer #1 was drilled but due to geologic problems, the project team was only able to drill to 1890 meters (6200 feet) instead of the proposed final depth 2744 meters (9002 feet). The drilling goals outlined in the original scope of work were not met. The project was only able to obtain 18 meters (59 feet) of conventional core and 106 rotary sidewall cores. All the conventional core and sidewall cores were in sandstone. We were unable to core any potential igneous caprock. Petrographic analysis of the conventional core and sidewall cores determined that the average porosity of the sedimentary material was 3.4% and the average permeability was 0.065 millidarcy. Compaction and diagenetic studies of the samples determined there would not be any porosity or permeability at depth in SC. In Georgia there appears to be porosity in the J/TR section based on neutron log porosity values. The only zones in Rizer #1 that appear to be porous were fractured diabase units where saline formation water was flowing into the borehole. Two geocellular models were created for the SC and GA study area. Flow simulation modeling was performed on the SC data set. The injection simulation used the newly acquired basin data as well as the Petrel 3-D geologic model that included geologic structure. Due to the new basin findings as a result of the newly acquired data, during phase two of the modeling the diabase unit was used as reservoir and the sandstone units were used as caprock. Conclusion are: 1) the SGR basin is composed of numerous sub-basins, 2) this study only looked at portions of two sub-basins, 3) in SC, 30 million tonnes of CO 2 can be injected into the diabase units if the fracture network is continuous through the units, 4) due to the severity of the faulting there is no way of assuring the injected CO 2 will not migrate upward into the overlying Coastal Plain aquifers, 5) in Georgia there appears to porous zones in the J/TR sandstones, 6) as in SC there is faulting in the sub-basin and the seismic suggest the faulting extends upward into the Coastal Plain making that area not suitable for CO 2 sequestration, 7) the complex faulting observed at both study areas appear to be associated with transfer fault zones (Heffner 2013), if sub-basins in the Georgia portion of the SGR basin can be located that are far away from the transfer fault zones there is a strong possibility of sequestering CO 2 in these areas, and 9) the SGR basin covers area in three states and this project only studied two small areas so there is enormous potential for CO 2 sequestration in other portions the basin and further research needs to be done to find these areas.

58 GEOSCIENCES↗

Generalization of mixed multiscale finite element methods with applications

Many science and engineering problems exhibit scale disparity and high contrast. The small scale features cannot be omitted in the physical models because they can affect the macroscopic behavior of the problems. However, resolving all the scales in these problems can be prohibitively expensive. As a consequence, some types of model reduction techniques are required to design efficient solution algorithms. For practical purpose, we are interested in mixed finite element problems as they produce solutions with certain conservative properties. Existing multiscale methods for such problems include the mixed multiscale finite element methods. We show that for complicated problems, the mixed multiscale finite element methods may not be able to produce reliable approximations. This motivates the need of enrichment for coarse spaces. Two enrichment approaches are proposed, one is based on generalized multiscale finte element metthods (GMsFEM), while the other is based on spectral element-based algebraic multigrid (rAMGe). The former one, which is called mixed GMsFEM, is developed for both Darcy’s flow and linear elasticity. Application of the algorithm in two-phase flow simulations are demonstrated. For linear elasticity, the algorithm is subtly modified due to the symmetry requirement of the stress tensor. The latter enrichment approach is based on rAMGe. The algorithm differs from GMsFEM in that both of the velocity and pressure spaces are coarsened. Due the multigrid nature of the algorithm, recursive application is available, which results in an efficient multilevel construction of the coarse spaces. Stability, convergence analysis, and exhaustive numerical experiments are carried out to validate the proposed enrichment approaches. iii

97 MATHEMATICS AND COMPUTING↗

Advanced gasifier and water gas shift technologies for low cost coal conversion to high hydrogen syngas

The Gas Technology Institute (GTI) and team members RTI International (RTI), Coanda Research and Development, and Nexant, are developing and maturing a portfolio of technologies to meet the United States Department of Energy (DOE) goals for lowering the cost of producing high hydrogen syngas from coal for use in carbon capture power and coal-to-liquids/chemicals. This project matured an advanced pilot-scale gasifier, with scalable and commercially traceable components, to readiness for use in a first-of-a-kind commercially-relevant demonstration plant on the scale of 500-1,000 tons per day (TPD). This was accomplished through cold flow simulation of the gasifier quench zone transition region at Coanda and through an extensive hotfire gasifier test program on highly reactive coal and high ash/high ash fusion temperature coals at GTI. RTI matured an advanced water gas shift process and catalyst to readiness for testing at pilot plant scale through catalyst development and testing, and development of a preliminary design basis for a pilot scale reactor demonstrating the catalyst. A techno-economic analysis was performed by Nexant to assess the potential benefits of the gasifier and catalyst technologies in the context of power production and methanol production. This analysis showed an 18%reduction in cost of power and a 19%reduction in cost of methanol relative to DOE reference baseline cases.

01 COAL, LIGNITE, AND PEAT↗

Analysis of Factors that Influence Infiltration Rates using the HELP Model

The Hydrologic Evaluation of Landfill Performance (HELP) model is used by Savannah River National Laboratory (SRNL) in conjunction with PORFLOW groundwater flow simulation software to make longterm predictions of the fate and transport of radionuclides in the environment at radiological waste sites. The work summarized in this report supports preparation of the planned 2018 Performance Assessment for the E-Area Low-Level Waste Facility (LLWF) at the Savannah River Site (SRS). More specifically, this project focused on conducting a sensitivity analysis of infiltration (i.e., the rate at which water travels vertically in soil) through the proposed E-Area LLWF closure cap. A sensitivity analysis was completed using HELP v3.95D to identify the cap design and material property parameters that most impact infiltration rates through the proposed closure cap for a 10,000-year simulation period. The results of the sensitivity analysis indicate that saturated hydraulic conductivity (Ksat) for select cap layers, precipitation rate, surface vegetation type, and geomembrane layer defect density are dominant factors limiting infiltration rate. Interestingly, calculated infiltration rates were substantially influenced by changes in the saturated hydraulic conductivity of the Upper Foundation and Lateral Drainage layers. For example, an order-of-magnitude decrease in Ksat for the Upper Foundation layer lowered the maximum infiltration rate from a base-case 11 inches per year to only two inches per year. Conversely, an order-of-magnitude increase in Ksat led to an increase in infiltration rate from 11 to 15 inches per year. This work and its results provide a framework for quantifying uncertainty in the radionuclide transport and dose models for the planned 2018 E-Area Performance Assessment. Future work will focus on the development of a nonlinear regression model for infiltration rate using Minitab 17® to facilitate execution of probabilistic simulations in the GoldSim® overall system model for the E-Area LLWF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗

NRAP-Open-IAM: FutureGen2 Component Models

This report describes the development and testing of three component models for NRAP-Open-IAM, the National Risk Assessment Partnership’s open-source integrated assessment model. The FutureGen2 Lookup Table Reservoir component model is based on interpolation of data from a set of lookup tables. The lookup tables contain pressures and saturations predicted by multiphase flow simulations performed with a reservoir simulator. The FutureGen2 Above Zone Monitoring Interval (AZMI) component is a surrogate model that can be used to estimate the impact that carbon dioxide (CO 2 ) and brine leaks from the CO 2 storage reservoir at the FutureGen 2.0 site might have had on overlying aquifers or monitoring units were a leak to occur. The model estimates the size of “impact plumes” according to five metrics: pH, Total Dissolved Solids (TDS), pressure, dissolved CO 2 and temperature. The FutureGen2 Aquifer component is similar, but is limited to four metrics: pH, Total Dissolved Solids (TDS), pressure, and dissolved CO 2 . The input parameters for each model are the same, but the Aquifer component is applicable to depths between 100 m and 700 m and the AZMI component is applicable from depths between 700 m and 1050 m.

42 ENGINEERING↗

Performance tuning of CEED software and 1st and 2nd wave apps

The goal of this milestone was the performance tuning of the CEED software, as well as the use and tuning of CEED to accelerate the first and second wave of targeted ECP applications. In this milestone, the CEED team developed optimization techniques and tuned for performance the CEED software to accelerate the first and second wave target ECP applications. Specifically, the focus was on the following: (1) Efficient use of the memory sub-system for managing data motion and interactions among different physics packages and libraries; this included integration of the new developments in a libCEED 0.5 software release; (2) Enhancing the CEED libraries with GPU and AMD GPU support in close interaction with vendors; (3) Optimal data locality and motion, and enhanced scalability and parallelism; this included vendors interactions for improvements in data motion and making strong scaling easier and more efficient; and (4) Continue boosting performance for the first-wave and second-wave of ECP target applications, including ExaSMR, ExaWind, Urban, MARBL, E3SM, and ExaAM. A main part of this milestone was the performance tuning of the CEED first-wave and second-wave ECP applications. This included the ExaSMR application – Coupled Monte Carlo Neutronics and Fluid Flow Simulation of Small Modular Reactors (ORNL), the MARBL application – Next-Gen Multi-physics Simulation Code (LLNL), and the ExaAM ExaConstit – a miniapp for the Transforming Additive Manufacturing through Exascale Simulation applications project (ExaAM), as well as ExaWind, Urban, and E3SM. The artifacts delivered include performance improvements in CEED’s 1st and 2nd wave of applications, and tuned CEED software for various architectures through a number of backends, freely available in the CEED’s repository on GitHub. See the CEED website, http://ceed.exascaleproject.org and the CEED GitHub organization, http://github.com/ceed for more details. In addition to details and results from the above R&D efforts, in this document we are also reporting on other project-wide activities performed in Q4 of FY19 including: CEED’s third annual meeting, MFEM GPU kernels scaling improvements, libCEED 0.5 and NekRS releases, an overview paper covering latest results on a set of bake-off high-order finite element problems (BPs), and other outreach efforts.

97 MATHEMATICS AND COMPUTING↗

Performance tuning of CEED software and first wave apps

The goal of this milestone was the performance tuning of the CEED software and first wave apps. In this milestone, the CEED team developed optimization techniques and tuned for performance the CEED software. Specifically, the focus was on the following: (1) Fast finite element operator storage and evaluation using partial assembly/matrix-free algorithms that take advantage of tensor-product element structure; (2) Architecture optimizations, including the development of performant discretization libraries targeting heterogeneous systems with multi-core CPUs, GPUs, and/or many-core CPUs in single node, cluster, or large-scale ECP supercomputer configurations; and (3) Global kernels for finite element operators, including the development of algorithms and libraries for fast gather-scatter exchanges for action assembly on unstructured graphs and for components needed for global linear solvers. A main part of this milestone was the performance tuning of the CEED first-wave ECP applications. This included the ExaSMR application – Coupled Monte Carlo Neutronics and Fluid Flow Simulation of Small Modular Reactors (ORNL), and the MARBL application – Next-Gen Multi-physics Simulation Code (LLNL). The artifacts delivered include performance improvements in CEED’s 1st wave of applications, and tuned CEED software for various architectures through a number of backends, freely available in the CEED’s repository on GitHub. See the CEED website, http://ceed.exascaleproject.org and the CEED GitHub organization, http://github.com/ceed for more details. In addition to details and results from the above R&D efforts, in this document we are also reporting on other project-wide activities performed in Q4 of FY18 including: CEED’s second annual meeting, a successful minisymposium at the premier international conference on high-order methods (ICOSAHOM), making the CEED milestone reports publicly available, new Laghos release, and other outreach efforts.

97 MATHEMATICS AND COMPUTING↗

NRAP-Open-IAM: Generic Aquifer Component Development and Testing

The Generic Aquifer Model calculates the concentrations of dissolved salt and dissolved CO 2 surrounding a leaking legacy well. The Generic Aquifer model can also estimate the size of an “impact plume” where concentration changes exceed user-specified thresholds. The model is a component of NRAP-Open-IAM, an open-source Integrated Assessment Model (IAM) developed by the National Risk Assessment Partnership (NRAP) to perform risk assessment for geologic CO 2 storage. The input parameters were selected to cover a wide range of groundwater aquifers and leakage rates. The generic aquifer model was developed using a generative adversarial deep learning network, trained using a large synthetic dataset of STOMP multiphase flow simulations. The deep learning model predictions of dissolved salt and dissolved CO 2 in the aquifer compare well to the original STOMP simulation results. The extent of aquifer impacted by leaking CO 2 or brine is calculated using a user-defined mass fraction threshold. The aquifer impact volumes calculated based on STOMP simulation results compare well to those calculated based on the deep learning model. In a provided python script, gridded observation results from the generic aquifer component of NRAP-Open-IAM are converted to HDF5 format files for monitoring design with the DREAM code.

54 ENVIRONMENTAL SCIENCES↗

Composite Analysis Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Sources.

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The parallel exascale version of the Subsurface Transport Over Multiple Phases (eSTOMP1) simulator is used to simulate flow and transport for the CA VZ models. The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Composite Analysis Vadose Zone Transport of Selected Hanford Site Composite Analysis: Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Waste Sites

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The parallel exascale version of the Subsurface Transport Over Multiple Phases (eSTOMP) simulator is used to simulate flow and transport for the CA VZ models. The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗