Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer programming integrity”

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,333 records · Page 74

Fast neutron leakage spectra of the EUCLID experiment

Special nuclear material in sub-critical and critical configurations measured in integral experiments are important for validation and adjustment of nuclear data. Many different evaluations of nuclear data exist, and these different evaluations can provide different values for individual cross sections that vary due to the uncertainties in differential experiments or lack of such data. For integral experiments, differences in these individual cross sections can have compensating errors, which lead to the same answer. One example of this is the Jezebel critical assembly, where k eff of the system is correctly computed by both ENDF/B-VIII.0 and JEFF-3.3, despite having substantially different underlying evaluated values for specific reactions (such as elastic and inelastic cross sections). To reduce compensating errors in nuclear data, the Experiments Underpinned by Computational Learning for Improvements in Nuclear Data (EUCLID) project has utilized machine learning to design a set of sub-critical and critical experiments. These experiments include slab- and cube-like configurations of 239 Pu in the form of the ZPPR plates. Six different responses were measured on a total of thirteen different configurations. One of these responses, the neutron leakage spectrum, was measured using an EJ301D detector. Finally, the results of the neutron leakage spectra show good agreement (within 1–2 σ ) with the expected spectrum from simulations and will be used in the subsequent nuclear data adjustment done by the EUCLID team.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multistage economic MPC for systems with a cyclic steady state: A gas network case study

Multistage model predictive control (MPC) provides a robust control strategy for dynamic systems with uncertainties and a setpoint tracking objective. Moreover, extending MPC to minimize an economic cost instead of tracking a pre-calculated optimal setpoint improves controller performance. This paper presents a novel multistage economic nonlinear model predictive control (E-NMPC) framework for dynamic systems operating under uncertainty, with specific application to natural gas transmission networks. A key innovation lies in the integration of cyclic steady-state (CSS) constraints within the multistage MPC formulation, enabling the controller to manage periodic operating conditions commonly observed in energy systems. A Lyapunov-based descent condition is enforced to ensure robust stability of the controller. The multistage economic MPC framework is validated on two gas pipeline case studies, where it successfully minimizes net energy consumption, respects operational constraints under uncertain demand profiles, and guides the network to optimal cyclic operation. The Lyapunov function remains bounded in both case studies, validating the robust stability of multistage E-NMPC.

03 NATURAL GAS↗

Phenomena Identification and Ranking Table (PIRT) for Heat Pipes

This Phenomena Identification and Ranking Table (PIRT) report provides an evaluation of key phenomena affecting the performance and operational regimes of heat pipes, particularly in the context of heat pipe microreactors (HPMRs). Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in the heat pipe, including phase change, turbulent transition, and compressibility effects, among others, there is high uncertainty in identifying and ranking the important phenomena affecting the operation of heat pipes and the current knowledge for their modeling and simulation and experimental measurements and instrumentation. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. The report analyzes phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics, discussing challenges and future research directions for improving their modeling and simulation and experimental measurements. Additionally, the report addresses phenomena with low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geysering, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. This comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Validation and Calibration of Nuclear Thermal Hydraulics Multiscale Multiphysics Models - Subcooled Flow Boiling Study

In addition to validation data plan, development of advanced techniques for calibration and validation of complex multiscale, multiphysics nuclear reactor simulation codes are a main objective of the CASL VUQ plan. Advanced modeling of LWR systems normally involves a range of physico-chemical models describing multiple interacting phenomena, such as thermal hydraulics, reactor physics, coolant chemistry, etc., which occur over a wide range of spatial and temporal scales. To a large extent, the accuracy of (and uncertainty in) overall model predictions is determined by the correctness of various sub-models, which are not conservation-laws based, but empirically derived from measurement data. Such sub-models normally require extensive calibration before the models can be applied to analysis of real reactor problems. This work demonstrates a case study of calibration of a common model of subcooled flow boiling, which is an important multiscale, multiphysics phenomenon in LWR thermal hydraulics. The calibration process is based on a new strategy of model-data integration, in which, all sub-models are simultaneously analyzed and calibrated using multiple sets of data of different types. Specifically, both data on large-scale distributions of void fraction and fluid temperature and data on small-scale physics of wall evaporation were simultaneously used in this work’s calibration. In a departure from traditional (or common-sense) practice of tuning/calibrating complex models, a modern calibration technique based on statistical modeling and Bayesian inference was employed, which allowed simultaneous calibration of multiple sub-models (and related parameters) using different datasets. Quality of data (relevancy, scalability, and uncertainty) could be taken into consideration in the calibration process. This work presents a step forward in the development and realization of the “CIPS Validation Data Plan” at the Consortium for Advanced Simulation of LWRs to enable quantitative assessment of the CASL modeling of Crud-Induced Power Shift (CIPS) phenomenon, in particular, and the CASL advanced predictive capabilities, in general. This report is prepared for the Department of Energy’s Consortium for Advanced Simulation of LWRs program’s VUQ Focus Area.

97 MATHEMATICS AND COMPUTING↗

MCNP6 Results for the Phase III Sensitivity Benchmark of the OCED/NEA Expert Group on Uncertainty Analysis for Criticality Safety Assessment

Within the last decade, there has been increasing interest in the calculation of cross section sensitivity coefficients of k eff for integral experiment design and uncertainty analysis. The OECD/NEA has an Expert Group devoted to Sensitivity and Uncertainty Analysis within the Working Party for Nuclear Criticality Safety. This expert group has developed benchmarks to assess code capabilities and performance for doing sensitivity and uncertainty analysis. Phase III of a set of sensitivity benchmarks evaluates capabilities for computing sensitivity coefficients. MCNP6 has the capability to compute cross section sensitivities for k eff using continuous-energy physics. To help verify this capability, results for the Phase III benchmark cases are generated and submitted to the Expert Group for comparison. The Phase III benchmark has three cases: III.1, an array of MOX fuel pins, III.2, a series of infinite lattices of MOX fuel pins with varying pitches, and III.3 two spheres with homogeneous mixtures of UF 4 and polyethylene with different enrichments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Status of PyGriffin Development for Integration into the NEAMS Workbench

The integration of Griffin into Workbench, by way of the PyGriffin code package, was initiated by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY-2022. PyGriffin was developed as a Python wrapper to Griffin to streamline the complex workflow involving mesh, cross section, and Griffin input generation, code execution, and results postprocessing, by improving user experience with high fidelity neutronics analysis through the Workbench GUI interface. PyGriffin can be used as a standalone application or through the PyARC code suite. FY-2023 saw the continued development of PyGriffin to expand its capabilities including advanced post-processing of simulation results and initial development of the Monte Carlo (MC) cross section generation workflow with the Shift MC code. PyGriffin was also approved for open-source software status in FY-2023, initiating plans to move the code repository into an open location. In addition to PyGriffin development, there were many improvements and developments made in PyARC leading to several software releases (latest version is v2.3.0), but the focus of this report is PyGriffin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Residuals-based distributionally robust optimization with covariate information

We consider data-driven approaches that integrate a machine learning prediction model within distributionally robust optimization (DRO) given limited joint observations of uncertain parameters and covariates. Our framework is flexible in the sense that it can accommodate a variety of regression setups and DRO ambiguity sets. We investigate asymptotic and finite sample properties of solutions obtained using Wasserstein, sample robust optimization, and phi-divergence-based ambiguity sets within our DRO formulations, and explore cross-validation approaches for sizing these ambiguity sets. Through numerical experiments, we validate our theoretical results, study the effectiveness of our approaches for sizing ambiguity sets, and illustrate the benefits of our DRO formulations in the limited data regime even when the prediction model is misspecified.

97 MATHEMATICS AND COMPUTING↗

Optimal Control of SOEC-Based Hydrogen Production Systems for Demand Response Using Deep Reinforcement Learning in Smart Grids

Solid oxide electrolysis cell (SOEC) hydrogen production technology can range in size from small, appliance-size equipment to large-scale, central production facilities that can be tied directly to renewable or non-greenhouse-gas-emitting forms of electricity production, making it an ideal resource for demand response (DR). The SOEC hydrogen production system is a complex integrated system that encompasses fluid dynamics, electrical dynamics, and electrochemical and thermal dynamics, all of which involve non-linearity and non-convexity. Proper control of the SOEC hydrogen production system is crucial to enable its participation in the DR program. Here, to overcome the difficulty of designing an explicit control law for such nonlinear systems with nonconvex optimization features in DR applications, deep reinforcement learning (DRL) is explored to achieve the optimal control of the SOEC system for DR participation. Specifically, a twin delayed deterministic policy gradient (TD3) control framework is applied to achieve optimal response performance during DR events by considering power tracking error and hydrogen production efficiency with a suitable reward function. Two case studies with grid connections for tracking different DR commands were investigated. The first case study involved operating conditions reaching the boundaries, while the second involved operating conditions within the boundaries. The results showed that the proposed DRL-based control for SOEC can track the DR signal in a timely manner while maintaining high energy efficiency.

08 HYDROGEN↗

Collaborative Proposal: Improving Decadal Prediction of Arctic Climate Variability and Change Using a Regional Arctic System Model (RASM)

RASM is a multi-disciplinary project, which brings together researchers from six state universities, one military postgraduate school, and one DoE laboratory to address the core modeling objectives of the arctic research community articulated in the Arctic System Modeling report by Roberts et al. (2010b). This report advocates the construction of a regional downscaling tool to generate probabilistic decadal projections of Greenland ice sheet retreat, evolution of arctic sea ice cover, changes in land surface vegetation, and regional processes leading to arctic amplification. Unified coupled models such as RASM are ideal for this purpose because they simulate fine-scale physics, essential for the realistic representation of intra-annual variability, in addition to processes fundamental to long term climatic shifts (Hurrell et al. 2009). By using RASM with boundary conditions from a global model, we can generate many-member ensembles essential for understanding uncertainty in regional climate projections (Hawkins and Sutton 2009). This probabilistic approach is computationally prohibitive for high-resolution global models in the foreseeable future, and also for regional models interactively nested within global simulations. Yet it is fundamental for quantifying uncertainty in decadal forecasts to make them useful for decision makers (Doherty et al. 2009). For this reason, we have targeted development of ensemble generation techniques as a core project task (Task 4.5). Environmental impact assessment specialists need high-fidelity regional ensemble projections to improve the accuracy of their work (Challinor et al. 2009; Moss et al. 2010). This is especially true of the Arctic, where economic, social and national interests are rapidly reshaping the high north in step with regional climate change. During the next decade, considerable oil and gas discoveries are expected across many parts of the marine and terrestrial Arctic (Gautier et al. 2009), the economics of the Northern Sea Route will steadily improve (Arctic Council 2009), and sovereign claims over the Arctic Ocean will increasingly be subject to international negotiations (Proelss 2009). Issues such as these have led to an expanding demand for Arctic climate projections to aid national and commercial decisions. However, detailed information from existing models is lacking. RASM will enhance the existing Arctic system modeling capabilities and align them with the scientific and societal needs outlined above. The science involved in the development of RASM will be integrated with teaching and training at a few different levels. RASM PIs will supervise at least six postdoctoral and doctoral students with additional involvement of undergraduate research assistants. RASM postdoctoral fellows and graduate trainees will benefit from being a part of a large, collaborative, and multi-disciplinary research projects employing state-of-the-art modeling and computational tools. They will also have the opportunity to contribute to the UTEP-led education and outreach effort developed specifically for RASM. The UTEP team will use RASM as a platform for education and outreach. First, we will develop products for public dissemination, such as curriculum units, lesson plans and other materials (simulations, movies, and images) for use by students and teachers in high school and university classrooms. Second, we will facilitate bringing RASM PIs, postdocs, and graduate students into the classroom, through electronic mentorship and by contributing content via online lectures and presentations. Third, by the use of problem-based learning (PBL) approaches, we will provide real-world scenarios and problems enabling students to do research and develop position papers or presentations on topics related to RASM. (Problem-based learning is a student-centered, inquiry-based approach in which students work in teams to solve challenging, open-ended problems.) Fourth, we will develop teacher-training materials that will be developed into workshops (face-to-face or online) to help teachers understand how to use the materials we supply. Overall, we will provide a broad, and long-term RASM program for education and outreach. RASM education and outreach activities will target underrepresented populations - in particular, Hispanic students and teachers. For example, the UTEP student body is more than 70 percent Hispanic and 55 percent female, and the majority of students are the first in their families to attend college. Students from underrepresented groups will be given the chance to expand their understanding of the different scientific knowledge and processes associated with RASM. We will develop an eye-catching RASM website to include information about the project, model results and educational materials, and lists of RASM publications and presentations.

54 ENVIRONMENTAL SCIENCES↗

CRADA Final Report: Open Microgrid Platform

Current microgrid control technology is mostly proprietary, expensive, and inflexible. An open-source, vendor-neutral, publish-subscribe controller will significantly lower the barrier to entry for new technologies and market entrants: inverter, battery, and load control manufacturers; cloud and services providers; and energy aggregators. An OpenFMB controller can significantly reduce the acquisition, integration, and ongoing operating and maintenance costs and increase revenue opportunities for energy producers and consumers. It allows distributed energy resource (DER) owners to replace individual DER components as richer-function, lower-cost devices become available or to incorporate new forecast and optimization algorithms as they are developed. A standard, secure, full-function DER field controller will also provide a low-cost means for grid operators to effectively manage the variability and uncertainty of solar photovoltaics (PV). Under the terms of this CRADA, ORNL has worked with Open Energy Solutions (OES) to advance the transition of ORNL’s existing open-source microgrid controller into a platform geared toward mainstream implementation. Using a more prevalent and memory-safe programming language, Rust, OES has converted ORNL’s current generation on-grid optimization into an application more suitable for industry use. The results of the ported optimization were analyzed by inputting the same arguments in both the current-generation and next-generation optimizers and verifying that the results were equivalent.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Risk-Informed Safety Margins Characterization (RISMC) Pathway Technical Program Plan

Safety is central to the design, licensing, operation, and economics of nuclear power plants (NPPs). As the current Light Water Reactor (LWR) NPPs age beyond 60 years, there are possibilities for increased frequency of systems, structures, and components (SSCs) degradations or failures that initiate safety-significant events, reduce existing accident mitigation capabilities, or create new failure modes. Plant designers commonly “over-design” portions of NPPs and provide robustness in the form of redundant and diverse engineered safety features to ensure that, even in the case of well-beyond design basis scenarios, public health and safety will be protected with a very high degree of assurance. This form of defense-in-depth is a reasoned response to uncertainties and is often referred to generically as “safety margin.” Historically, specific safety margin provisions have been formulated, primarily based on “engineering judgment.” The ability to better characterize (i.e., describe and quantify) safety margin is important to improved decision making about LWR design and long-term operation. A systematic approach to characterizing safety margins and the subsequent risk informed margins management options represents a vital input to the licensee and regulatory analysis and decision making that will be involved. In addition, as research and development (R&D) in the LWRS Program and other collaborative efforts yield new data and improved scientific understanding of physical processes that govern the aging and degradation of plant SSCs (and concurrently support technological advances in nuclear reactor fuels and plant instrumentation and control systems), needs and opportunities to better optimize plant safety, economics, and performance will become known. The purpose of the Risk-Informed Safety Margins Characterization (RISMC) Pathway R&D is to support plant decisions for risk-informed margins management with the aim to improve economics, reliability, and sustain safety of current NPPs over periods of extended plant operations. The goals of the RISMC Pathway are twofold: (1) develop and demonstrate a risk-assessment method that is coupled to safety margin quantification that can be used by NPP decision makers as part of risk-informed margin management strategies; (2) create an advanced RISMC Toolkit that enables more accurate representation of NPP safety margins. The methods and tools provided by RISMC are essential to a comprehensive and integrated risk-informed margin management approach that supports effective preservation of margin for both active and passive SSCs.

97 MATHEMATICS AND COMPUTING↗

Waveform processing using neural network algorithms on the front-end electronics

In a multi-channel radiation detector readout system, waveform sampling, digitization, and raw data transmission to the data acquisition system constitute a conventional processing chain. The deposited energy on the sensor is estimated by extracting peak amplitudes, area under pulse envelopes from the raw data, and starting times of signals or time of arrivals. However, such quantities can be estimated using machine learning algorithms on the front-end Application-Specific Integrated Circuits (ASICs), often termed as “edge computing”. Edge computation offers enormous benefits, especially when the analytical forms are not fully known or the registered waveform suffers from noise and imperfections of practical implementations. In this work, we aim to predict peak amplitude from a single waveform snippet whose rising and falling edges containing only 3 to 4 samples. We thoroughly studied two well-accepted neural network algorithms, Multi-Layer Perceptron (MLP) and Convolutional Neural Network (CNN) by varying their model sizes. Further, to better fit front-end electronics, neural network model reduction techniques, such as network pruning methods and variable-bit quantization approaches, were also studied. By combining pruning and quantization, our best performing model has the size of 1.5 KB, reduced from 16.6 KB of its full model counterpart. It can reach mean absolute error of 0.034 comparing to that of a naive baseline of 0.135. Such parameter-efficient and predictive neural network models established feasibility and practicality of their deployment on front-end ASICs.

47 OTHER INSTRUMENTATION↗

Flexible Integration of Diverse HVAC Technologies in EnergyPlus via Python-Enabled Workflows

Analysis of advanced controls and novel system types is often not directly feasible in building energy simulation tools. Various techniques extend building energy simulation tool capabilities to allow the use of user-defined scripts and programs, but these approaches have limitations. The EnergyPlus Python plugin offers users new flexibility to use EnergyPlus to call an external Python module at specific points in the simulation, as well as to use Python to call EnergyPlus functionality through an application programming interface (API). This paper presents four case studies leveraging the EnergyPlus Python plugin to facilitate analysis of advanced controls and system types. The use of the Python plugin offers greater modularity and flexibility relative to previous approaches, is less error prone, and is simpler for users to adopt. The Python plugin allows EnergyPlus to be used in a more flexible manner and to accommodate the expanding realm of energy modeling applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Liquid Metals as Plasma-facing Materials for Fusion Energy Systems: From Atoms to Tokamaks

The objective of our studies was to advance our fundamental understanding of liquid metals as plasma-facing materials for fusion energy systems, with a broad scope: from atoms to tokamaks. The flow of liquid metals offers solutions to significant problems of the plasma-facing materials for fusion energy systems. Candidate metals include lithium, tin, gallium, and their eutectic combinations. However, such liquid metal solutions can only be designed efficiently if a range of scientific and engineering issues are resolved that require advances in fundamental fluid dynamics, materials science and surface science. In our research we investigated a range of significant and timely problems relevant to current and proposed engineering designs for fusion reactors, including high-heat flux configurations that are being considered by leading fusion energy groups world-wide. Using experimental and theoretical tools spanning atomistic to continuum descriptions of liquid metals, and bridging surface chemistry, wetting/dewetting and flow, our research has advanced the science and engineering of fusion energy materials and systems. Specifically, we developed a combined experimental and theoretical program to investigate flows of liquid metals in fusion-relevant geometries, including equilibrium and stability of thin-film flows, e.g. wetting and dewetting, effects of electromagnetic and thermocapillary fields on liquid metal thin-film flows, and how chemical interactions and the properties of the surface are influenced by impurities and in turn affect the surface wetting characteristics, the surface tension, and its gradients. Because high-heat flux configurations produce evaporation and sputtering, which forces rearrangement of the liquid, and any dewetting exposes the substrate to damage from the plasma, our studies addressed such evaporatively driven liquid flows and measured and simulated properties of the different bulk phases and material interfaces. The range of our studies included (i) quantum mechanical calculations that allow inclusion of many thousands of atoms for the characterization of the interface of liquid metals exposed to continuous bombardment by deuterium and tritium as expected in fusion, (ii) molecular dynamics studies of the phase behavior of liquid metals, which (a) utilize thermodynamic properties computed using our quantum mechanical calculations and (b) establish material and wetting properties of the liquid metals, including relevant eutectics, (iii) experimental investigations of the surface science of liquid metals, interacting both with the solid substrate as well as gaseous species, and (iv) fluid dynamical studies that incorporate the material and surface science results of (ii) and (iii) in order to characterize flow in capillary porous materials and the thin-film flow along curved boundaries, both of which are potentially major components of plasma-facing materials. The outcome of these integrated studies was new understanding that enables developing design rules useful for future developments of the plasma-facing components critical to the success of fusion energy systems.

36 MATERIALS SCIENCE↗

Preface to the Special Issue on TOUGH Symposium 2015

The TOUGH Symposium 2015 was held in Berkeley, California, September 28-30, 2015. The TOUGH family of codes, developed at the Energy Geosciences Division of Lawrence Berkeley National Laboratory (LBNL), is a suite of computer programs for the simulation of multiphase and multicomponent fluid and heat flows in porous and fractured media with applications in many geosciences fields, such as geothermal reservoir engineering, nuclear waste disposal, geological carbon sequestration, oil and gas reservoirs, gas hydrate research, vadose zone hydrology and environmental remediation. Since the first release in the 1980s, many modifications and enhancements have been continuously made to TOUGH and its various descendants (iTOUGH2, TOUGH+, TOUGH-MP, TOUGHREACT, TOUGH+HYDRATE, TMVOC...), at LBNL and elsewhere. Today, these codes are used worldwide in academia, government organizations and private companies in problems involving coupled hydrological, thermal, biogeochemical and geomechanical processes. The Symposia, organized every 2-3 years, bring together developers and users for an open exchange on recent code enhancements and applications. In 2015, the Symposium was attended by one hundred participants, representing thirty-four nationalities. This Special Issue in Computers & Geosciences gathers extended versions of selected Symposium proceedings related to (i) recent enhancements to the TOUGH family of codes and (ii) coupled flow and geomechanics processes modeling. In this special issue, a paper by Jung et al. provides a new base version of TOUGH, TOUGH3. This latest version addresses the computational robustness and efficiency needed for the increasing complexity of the simulated processes and the growing size of model domains. Moreover, it consolidates the serial and parallel versions into a single code, hence, inconsistencies between different versions are less likely to occur, and code maintenance and development become more efficient. Finsterle et al. present a review of the capabilities of iTOUGH2 and illustrate the workflow of an integrated simulation-optimization framework. Liu et al. present an approach to reduce the computational cost of performing Bayesian inverse modeling and demonstrate their approach using a synthetic experiment simulated with TOUGH2. Regarding the mesh generation process, the paper by Doughty presents a generalization of the grid-generation algorithm embedded in TOUGH2 (MESHMAKER) to create one-column grids with fractal dimension, useful for representing flow through fractures or highly heterogeneous media with fractal geometry, and also for inverse methods. Sentís and Gable present a module to interface the open source code Lagrit with TOUGH2, that assists in the processes of mesh generation and model setup, and that can be particularly useful for complex geometries and large models. Bonduà et al. present a pre-processing tool (VORO2MESH) that features a full 3D Voronoi approach for structured and unstructured grids, and an enhanced post-processor (TOUGH2Viewer) for the visualization of 3D Voronoi grids. The literature review shows that coupled fluid flow and geomechanics has gained increasing attention in the last years in many geosciences fields. In this context, several codes of the TOUGH family have been successfully coupled with geomechanics software. The paper by Rutqvist presents a review of some of these TOUGH-based geomechanics models, of which the TOUGH-FLAC simulator – presented for the first time in the TOUGH Symposium 2003 – is well-known worldwide and has been used in a large variety of applications, encompassing the application domains of the TOUGH codes. Blanco-Martín et al. present an extension of TOUGH-FLAC to the finite strain framework, motivated by the existence of large strains associated with rock salt deformation in the context of heat-generating nuclear waste disposal. Following a paper published in the Special Issue of the TOUGH Symposium 2012 in this journal, Kim et al. present the current status of the TOUGH-RSBN simulator for coupled hydro-mechanical modeling of hydraulic fracturing using a discrete fracture network approach, and present some validation and demonstration cases. Rinaldi and Nespoli present TOUGH-SEED, an injection-induced seismicity modeling approach using TOUGH2 and a geomechanics stochastic model, and apply it to the Basel EGS project. Finally, Rinaldi et al. use the PEST protocol implemented in iTOUGH2, combined with TOUGH-FLAC as the forward model, to perform inverse modeling for the In Salah CO2 demonstration site in Algeria. Using pore pressures and InSAR displacements as observations, this approach allows for a multi-physics joint inversion and takes advantage of the inverse capabilities of iTOUGH2 (parameter estimation, sensitivity analysis, uncertainty estimation …), applied to coupled fluid flow and geomechanics. On behalf of the members of the TOUGH Symposium 2015 Organizing Committee, I hope that this collection of papers will provide the readers of Computers & Geosciences with a recent view of the capabilities of the TOUGH family of codes, as well as with current modeling research domains in the geosciences. Last but not least, I would like to express my gratitude to all the reviewers who kindly contributed to this Special Issue, as well as the Journal staff who have made this publication possible.

99 GENERAL AND MISCELLANEOUS↗

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗

Countering Weapons of Mass Destruction (CWMD) Device Cybersecurity Characterization Process and Profile

Countering Weapons of Mass Destruction (CWMD) recognizes that threats in the cyberspace domain continue to grow, which requires CWMD devices and supporting systems to be both cybersecure (ability to protect or defend from cyber-attacks) and resilient (ability to maintain required capability in the face of adversity) to cyber threats. The CWMD cybersecurity characterization approach in this document supports existing cyber resilience activities within the Acquisition Lifecycle Framework. Similarly, this process supports existing Department of Homeland Security Cyber Resilience Test and Evaluation activities, which consist of iterative processes, starting at the initiation of system acquisition and continuing throughout the entire device and system life cycle. Cyber resilience is the ability of an information system to continue to operate while under attack, even if in a degraded or debilitated state, and to rapidly recover operational capabilities for essential functions after a successful attack. The goal of the security characterization task for CWMD is to support the development of a CBRN device-dependent profile that aligns with device network capabilities and maps to recommended security controls to create a characterization security profile impact levels. The impact levels for CWMD devices should be characterized as Low (L), Moderate (M), High (H) to align with the low, moderate, high control baselines. To estimate the impact levels, the device’s security-related attributes are translated into the security objectives: Confidentiality (C), Integrity (I), and Availability (A), known as the CIA triad. The potential impact for each device can be L, M, H, for devices that connect and transmit different types of data and may have different impact levels. National Institute of Standards and Technology Federal Information Processing Standards Publication 199 states, “the potential impact values assigned to the respective security objectives shall be the highest value from among those security categories that have been determined for each type of information resident on the information system.” As CWMD is determining the cybersecurity impact levels of CBRN devices based on network connections and data transfers, the impact levels are aligned with the associated attributes of network connections and communications. For example, if the device system is connected to a wireless network and transmits different data types based on the confidentiality of the data, the highest impact value for each security objective should represent the device’s CIA impact level. This document is intended to be used by test managers, test team, and program managers.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Mitigation of Syngas Cooler Plugging and Fouling

This Final Report summarizes research performed to develop a technology to mitigate the plugging and fouling that occurs in the syngas cooler used in many Integrated Gasification Combined Cycle (IGCC) plants. The syngas cooler is a firetube heat exchanger located downstream of the gasifier. It offers high thermal efficiency, but its’ reliability has generally been lower than other process equipment in the gasification island. The buildup of ash deposits that form on the fireside surfaces in the syngas cooler (i.e., fouling) lead to reduced equipment life and increased maintenance costs. Our approach to address this problem is that fouling of the syngas cooler cannot be eliminated, but it can be better managed. The research program was funded by DOE using two budget periods: Budget Period 1 (BP1) and Budget Period 2 (BP2). The project used a combination of laboratory scale experiments, analysis of syngas cooler deposits, modeling and guidance from industry to develop a better understanding of fouling mechanisms and to develop and evaluate strategies to mitigate syngas cooler fouling and thereby improve syngas cooler performance. The work effort in BP 1 and BP 2 focused on developing a better understanding of the mechanisms that lead to syngas cooler plugging and fouling and investigating promising concepts to mitigate syngas cooler plugging and fouling. The work effort focused on the following: • analysis of syngas cooler deposits and fuels provided by an IGCC plant collaborating with this project; • performing Jet cleaning tests in the University of Utah Laminar Entrained Flow Reactor to determine the bond strength between an ash deposit to a metal plate, as well as implementing planned equipment modifications to the University of Utah Laminar Entrained Flow Reactor and the one ton per day, pressurized Pilot Scale Gasifier; • performing Computational Fluid Dynamic modeling of industrially relevant syngas cooler configurations to develop a better understanding of deposit formation mechanisms; • performing Techno-Economic-Analysis for a representative IGCC plant to investigate the impact on plant economics, in particular the impacts on the Cost of Electricity (COE), due to plant shutdowns caused by syngas cooler plugging and fouling and potential benefits to plant economics of developing strategies to mitigate syngas cooler fouling; and • performing modeling and pilot scale tests to investigate the potential benefits of using a sorbent (fuel additive) to capture the vaporized metals that result in syngas cooler fouling. All project milestones for BP 1 and BP 2 were achieved. DOE was provided a briefing on our accomplishments in BP1 and BP2 and our proposed plans for Budget Period 3 (BP 3). Based on our research the mitigation technology selected to investigate in BP 3 was the use of a sorbent that can be injected into the gasifier with the fuel slurry to capture vaporized metals that lead to the deposit formation in the syngas cooler. The work effort proposed for BP 3 would have focused on addressing concerns raised by gasification industry personnel for the impacts on gasifier performance of sorbent injection, so that at the end of BP 3 the use of sorbent injection would be at “pre-commercial” stage and ready for use in a Field Demonstration that could be funded by industry or DOE. A Budget Continuation Application (BCA) was submitted to obtain funding for BP3 DOE but DOE chose to not fund the proposed BP3 effort.

47 OTHER INSTRUMENTATION↗