Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow Visualization”

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 325 records · Page 18

Quantitatively Monitoring Bubble-Flow at a Seep Site Offshore Oregon: Field Trials and Methodological Advances for Parallel Optical and Hydroacoustical Measurements

Two lander-based devices, the Bubble-Box and GasQuant-II, were used to investigate the spatial and temporal variability and total gas flow rates of a seep area offshore Oregon, United States. The Bubble-Box is a stereo camera–equipped lander that records bubbles inside a rising corridor with 80 Hz, allowing for automated image analyses of bubble size distributions and rising speeds. GasQuant is a hydroacoustic lander using a horizontally oriented multibeam swath to record the backscatter intensity of bubble streams passing the swath plain. The experimental set up at the Astoria Canyon site at a water depth of about 500 m aimed at calibrating the hydroacoustic GasQuant data with the visual Bubble-Box data for a spatial and temporal flow rate quantification of the site. For about 90 h in total, both systems were deployed simultaneously and pressure and temperature data were recorded using a CTD as well. Detailed image analyses show a Gaussian-like bubble size distribution of bubbles with a radius of 0.6–6 mm (mean 2.5 mm, std. dev. 0.25 mm); this is very similar to other measurements reported in the literature. Rising speeds ranged from 15 to 37 cm/s between 1- and 5-mm bubble sizes and are thus, in parts, slightly faster than reported elsewhere. Bubble sizes and calculated flow rates are rather constant over time at the two monitored bubble streams. Flow rates of these individual bubble streams are in the range of 544–1,278 mm 3 /s. One Bubble-Box data set was used to calibrate the acoustic backscatter response of the GasQuant data, enabling us to calculate a flow rate of the ensonified seep area (~1,700 m 2 ) that ranged from 4.98 to 8.33 L/min (5.38 × 10 6 to 9.01 × 10 6 CH 4 mol/year). Such flow rates are common for seep areas of similar size, and as such, this location is classified as a normally active seep area. For deriving these acoustically based flow rates, the detailed data pre-processing considered echogram gridding methods of the swath data and bubble responses at the respective water depth. The described method uses the inverse gas flow quantification approach and gives an in-depth example of the benefits of using acoustic and optical methods in tandem.

54 ENVIRONMENTAL SCIENCES↗

CFD-grade Experiments Addressing Thermal Fatigue in Nuclear Reactor Branch Lines

Thermal fatigue has been identified as one of the main degradation mechanisms contributing to Loss of Coolant Accident frequencies in Light Water Reactors (LWRs). Thermal fatigue in stagnant branch lines departing from hot or cold legs of LWRs has become of high interest to the nuclear industry, as the current management program has had mixed effectiveness in preventing loss of power incidents. In fact, at least 8 such events involving branch lines have been registered at US plants over the past 2 years. Thermal cycling in a branch line occurs at the fluctuating boundary created by the penetration of hot fluid from the main flow into the cold stagnant fluid of the branch line. Penetration occurs as the fast moving flow of the main line crosses the opening of a branch line, causing cavity flow. Turbulent cavity flow dissipates to a disturbed flow further into the branch line and then to a laminar spiral flow. At the extent of the spiral flow, hot and cold fluids stratify to create the fluctuating boundary of interest. Qualitative experimental techniques, such as visualization of tracer particles, have been employed in past research efforts to capture the large scale evolution from cavity flow to the formation of the fluctuating stratified layer. The insight gained has been crucial in understanding the fundamentals of thermal cycling and stratification in branch lines; however, this qualitative information is insufficient for improving the current methods of predicting thermal fatigue in branch lines. High resolution quantitative measurements of the flow, capturing the small scale motions that affect the thermal cycling at the pipe wall, are needed to shed more light on the phenomenon and improve the industry fatigue management program, especially in terms of reliability. In the work presented here, time-resolved particle image velocimetry (PIV) is employed to investigate the flow in a stagnant branch line. Both time-averaged and time-dependent velocity fields of interest were captured, providing the groundwork for the experimental campaign that will be translated into model enhancement. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Using ASCEM modeling and visualization to inform stakeholders of contaminant plume evolution and remediation efficacy at F-Basin Savannah river, SC - 15156

Communication with stakeholders, regulatory agencies, and the public is an essential part of implementing different remediation and monitoring activities, and developing site closure strategies at contaminated sites. Modeling of contaminant plume evolution plays a critical role in estimating the benefit, cost, and risk of particular options. At the same time, effective visualization of monitoring data and modeling results are particularly important for conveying the significance of the results and observations. In this paper, we present the results of the Advanced Simulation Capability for Environmental Management (ASCEM) project, including the discussion of the capabilities of newly developed ASCEM software package, along the its application to the F-Area Seepage Basins located in the U.S. Department of Energy Savannah River Site (SRS). ASCEM software includes state-of-the-art numerical methods for simulating complex flow and reactive transport, as well as various tool-sets such as a graphical user interface (GUI), visualization, data management, uncertainty quantification, and parameter estimation. Using this software, we have developed an advanced visualization of tritium plume migration coupled with a data management system, and simulated a three-dimensional model of flow and plume evolution on a high-performance computing platform. We evaluated the effect of engineered flow barriers on a nonreactive tritium plume, through advanced plume visualization and modeling of tritium plume migration. In addition, we developed a geochemical reaction network to describe complex geochemical processes at the site, and evaluated the impact of coupled hydrological and geochemical heterogeneity. These results are expected to support SRS's monitoring activities and operational decisions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Quantifying Dissolution Dynamics in Porous Media Using a Spatial Flow Focusing Profile

Abstract The diverse range of patterns in porous media formed by dissolution processes depends on the relative magnitude of flow, transport, and chemical reactions at pore surfaces. However, distinguishing between regimes often relies solely on qualitative, visual comparisons of emergent structures. Here, we propose a quantitative measure capable of identifying different regimes using the concept of the spatial flow focusing profile, which segments the medium into cross sections along the flow direction to calculate the flow focusing index for each section. We employ this measure in numerical simulations of a dissolving porous medium using a pore network model. We obtain a morphological phase diagram of dissolution patterns, which we characterize using the flow focusing profile. In particular, we demonstrate that analyzing the temporal changes in the profile allows one to quantitatively distinguish between wormholing and channeling. The transition between them is shown to be affected by the heterogeneity of the system.

58 GEOSCIENCES↗

Remote visual analysis of large turbulence databases at multiple scales

The remote analysis and visualization of raw large turbulence datasets is challenging. Current accurate direct numerical simulations (DNS) of turbulent flows generate datasets with billions of points per time-step and several thousand time-steps per simulation. Until recently, the analysis and visualization of such datasets was restricted to scientists with access to large supercomputers. The public Johns Hopkins Turbulence database simplifies access to multi-terabyte turbulence datasets and facilitates the computation of statistics and extraction of features through the use of commodity hardware. In this paper, we present a framework designed around wavelet-based compression for high-speed visualization of large datasets and methods supporting multi-resolution analysis of turbulence. By integrating common technologies, this framework enables remote access to tools available on supercomputers and over 230 terabytes of DNS data over the Web. Finally, the database toolset is expanded by providing access to exploratory data analysis tools, such as wavelet decomposition capabilities and coherent feature extraction.

97 MATHEMATICS AND COMPUTING↗

In situ feature analysis for large-scale multiphase flow simulations

The study of multiphase flow is essential for designing chemical reactors such as fluidized bed reactors (FBR), as a detailed understanding of hydrodynamics is critical for optimizing reactor performance and stability. An FBR allows scientists to conduct different types of chemical reactions involving multiphase materials, especially interaction between gas and solids. During such complex chemical processes, the formation of void regions in the reactor, generally termed as bubbles, is an important phenomenon. The study of these bubbles has a deep implication in predicting the reactor’s overall efficiency. But physical experiments needed to understand bubble dynamics are costly and non-trivial due to the technical difficulties involved and harsh working conditions of the reactors. Therefore, to study such chemical processes and bubble dynamics, a state-of-the-art computational simulation MFIX-Exa is being developed. Despite the proven accuracy of MFIX-Exa in modeling bubbling phenomena, the large-scale output data prohibits the use of traditional post hoc analysis capabilities in both storage and I/O time. Herein, to address these issues and allow the application scientists to explore the bubble dynamics in an efficient and timely manner, we have developed an end-to-end analytics pipeline that enables in situ detection of bubbles, followed by a flexible post hoc visual exploration methodology of bubble dynamics. The proposed method enables interactive analysis of bubbles, along with quantification of several bubble characteristics, enabling experts to understand the bubble interactions in detail. Positive feedback from the experts has indicated the efficacy of the proposed approach for exploring bubble dynamics in very-large-scale multiphase flow simulations.

97 MATHEMATICS AND COMPUTING↗

Water Observations of Flow/No-Flow for the East-Taylor Watershed, Colorado (June-July 2025 and 2026)

This dataset provides multi-year, ground-truth visual observations of surface water flow/no-flow conditions within the East-Taylor Watershed, Colorado, collected during June and July of 2025 and 2026. In June and July 2025, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function Scientific Focus Area (SFA) and Rocky Mountain Biological Laboratory (RMBL) Colorado Headwaters Ecological Spectroscopy Study (CHESS) campaign (further details are provided within the CHESS Project Description). We obtained 377 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. These ground-truth observations were collected to validate classification maps from remote sensing data and model results within the East-Taylor Watershed. In 2025, flow/no-flow measurements were collected using a field-based app for the CHESS Campaign (Zerion iForm). Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included information about visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, and beaver activity. For 2025 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2025_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2025_Water_Observations_Locations.kmz); (3) photos (.jpg and .jpeg) of the water observation points, organized by location, contained within 2025_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2025_Water_Observation_Protocols.pdf). In June and July 2026, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function SFA project. We obtained 365 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. The 2026 observations focused on collecting repeat measurements at the 2025 flow/no-flow observation locations conducted as part of the CHESS campaign. These ground-truth observations were collected to understand differences in flow/no-flow in 2026, given the unprecedented 2026 drought in Colorado. In 2026, flow/no-flow measurements were collected using ArcGIS (Geographic Information System) Survey123. Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included repeat information from the 2025 water observation effort, including visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, beaver activity, and a new metadata component of estimated stream depth (for select locations). For 2026 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2026_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2026_Water_Observations_Locations.kmz); (3) photos (.jpg) of the water observation points, organized by location, contained within 2026_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2026_Water_Observation_Protocols.pdf). For 2025 and 2026 water observations, this dataset contains: (1) a location metadata file (locations.csv); (6) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and (7) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. 2026-09-02: This dataset was updated to include 2026 water observation measurements. The 2025 observation files were also updated to ensure a consistent file naming convention across water observation years.

2018 NEON and 2025 CHESS Campaigns↗

Precipitation of calcium carbonate and calcium phosphate under diffusion controlled mixing

Multi-component mineral precipitation in porous, subsurface environments is challenging to simulate or engineer when in situ reactant mixing is controlled by diffusion. In contrast to well-mixed systems, the conditions that favor mineral precipitation in porous media are distributed along chemical gradients, which evolve spatially due to concurrent mineral precipitation and modification of solute transport in the media. The resulting physical and chemical characteristics of a mixing/precipitation zone are a consequence of coupling between transport and chemical processes, and the distinctive properties of individual chemical systems. We examined the spatial distribution of precipitates formed in “double diffusion” columns for two chemical systems, calcium carbonate and calcium phosphate. Polyacrylamide hydrogel was used as a low permeability, high porosity medium to maximize diffusive mixing and minimize pressure- and density-driven flow between reactant solutions. In the calcium phosphate system, multiple, visually dense and narrow bands of precipitates were observed that were reminiscent of previously reported Liesegang patterns. In the calcium carbonate system, wider precipitation zones characterized by more sparse distributions of precipitates and a more open channel structure were observed. In both cases, formation of precipitates inhibited, but did not necessarily eliminate, continued transport and mixing of the reactants. A reactive transport model with fully implicit coupling between diffusion, chemical speciation and precipitation kinetics, but where explicit details of nucleation processes were neglected, was able to qualitatively simulate properties of the precipitation zones. The results help to illustrate how changes in the physical properties of a precipitation zone depend on coupling between diffusion-controlled reactant mixing and chemistry-specific details of precipitation kinetics.

54 ENVIRONMENTAL SCIENCES↗

Grain Size Variations in the Murray Formation: Stratigraphic Evidence for Changing Depositional Environments in Gale Crater, Mars

The lowermost exposure of the Murray formation in Gale crater, Mars, was interpreted as sediment deposited in an ancient lake based on data collected by the Curiosity rover. Constraining the stratigraphic extent and duration of this environment has important implications for the paleohydrology of Gale. Insights into early Martian environments and paleofluid flow velocity can be obtained from grain size in rocks. Visual inspection of grain size is not always available for rocks investigated at field sites on Mars due to limited image coverage. But grain sizes can also be estimated from the Gini Index Mean Score, a grain-size proxy that uses point-to-point chemical variations in ChemCam Laser Induced Breakdown Spectroscopy data. Here, new Gini Index Mean Score results indicate that the Murray formation is dominated by mudstones with grains smaller than the spatial resolution of all rover cameras. Intervals of fine to coarse sandstone also are present, some of which are verified using observations of grain size and sedimentary structures in associated images. Overall, results demonstrate that most of the Murray consists of mudstone, suggesting settling of grains from suspension in low energy depositional environments such as lakes. Some of the mudstones contain desiccation cracks indicating periods of drying with a lowering of lake water level. However, beds and lenses of cross-bedded sandstones are common at specific intervals, suggesting episodes of fluvial and possibly eolian deposition. The persistence of lacustrine deposits interspersed with fluvial deposits suggests that liquid water was sustained on the Martian surface for tens of thousands to millions of years.

58 GEOSCIENCES↗

Experimental Plan for Crystal Accumulation Studies in the WTP Melter Riser

This experimental plan defines crystal settling experiments to be in support of the U.S. Department of Energy – Office of River Protection crystal tolerant glass program. The road map for development of crystal-tolerant high level waste glasses recommends that fluid dynamic modeling be used to better understand the accumulation of crystals in the melter riser and mechanisms of removal. A full-scale version of the Hanford Waste Treatment and Immobilization Plant (WTP) melter riser constructed with transparent material will be used to provide data in support of model development. The system will also provide a platform to demonstrate mitigation or recovery strategies in off-normal events where crystal accumulation impedes melter operation. Test conditions and material properties will be chosen to provide results over a variety of parameters, which can be used to guide validation experiments with the Research Scale Melter at the Pacific Northwest National Laboratory, and that will ultimately lead to the development of a process control strategy for the full scale WTP melter. The experiments described in this plan are divided into two phases. Bench scale tests will be used in Phase 1 (using the appropriate solid and fluid simulants to represent molten glass and spinel crystals) to verify the detection methods and analytical measurements prior to their use in a larger scale system. In Phase 2, a full scale, room temperature mockup of the WTP melter riser will be fabricated. The mockup will provide dynamic measurements of flow conditions, including resistance to pouring, as well as allow visual observation of crystal accumulation behavior.

36 MATERIALS SCIENCE↗

NRC HEAF Tests: Imaging & Measurement Methodology Report

The 2018 NRC HEAF tests were conducted in Chalfont Pennsylvania at KEMA High Power Laboratory during the week of September 10th. These scoping tests were executed to determine the most effective measurement methodologies for future tests. The goal of Sandia’s Photometrics group was to provide high-speed quantitative and qualitative imaging of the arcing fault tests for the Nuclear Regulatory Commission. The measurement methods included visible high-speed imaging, high-speed high-dynamic range visible imaging, thermal imaging, and quantitative flow imaging. In addition, data fusion products were generated to visualize instrumentation data and imaging measurements. All imaging has been time synchronized to the start of the arcing event.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

An Overview of the Circular Economy Lifecycle Assessment and VIsualization (CELAVI) Framework

A circular economy emphasizes the efficient use of all resources (e.g., materials, land, water). Despite anticipated overall benefits to society, the transition to a circular economy is likely to create regional differences in impacts. Current tools are unable to fully evaluate these potential externalities, which will be important for informing research prioritization and regional decision making. This presentation will review and contrast nine current impact assessment methods and describe their granularity, scope, data requirements, and capabilities within the context of identifying regional and sectoral transfers of impacts that could result from transitioning to a circular economy for energy systems. The review of current approaches highlights the need for a new hybrid circularity assessment framework that leverages multiple methods. We describe one such approach: the Circular Economy Lifecycle Assessment and VIsualization (CELAVI) framework. CELAVI uses system dynamics to model material flows for multiple circular economy pathways, network theory to track the spatial and sectoral flow of functional units across a graph, and discrete event simulation to to step through time and evaluate lifecycle assessment data at each time step. The framework is designed to be flexible and scalable enough to accommodate multiple energy materials and multiple energy technologies. CELAVI's hybridization of multiple existing methods yields additional capabilities, which may help answer questions about how flows between the technosphere and ecosphere may evolve if the circularity of energy systems is modified.

circular economy↗

Vortex breakdown in closed containers with polygonal cross sections

The vortex breakdown bubble in the confined flow generated by a rotating lid in closed containers with polygonal cross sections was analysed both experimentally and numerically for the height/radius aspect ratio equal to 2. The stagnation point locations of the breakdown bubble emergence and the corresponding Reynolds number were determined experimentally and in addition computed numerically by STAR-CCM+ CFD software for square, pentagonal, hexagonal, and octagonal cross section configurations. The flow pattern and the velocity were observed and measured by combining the seeding particle visualization and the temporal accuracy of laser Doppler anemometry. The vortex breakdown size and position on the container axis were determined for Reynolds numbers, ranging from 1450 to 2400. The obtained results were compared with the flow structure in the closed container of cubical and cylindrical configurations. It is shown that the measured evolution of steady vortex breakdown is in close agreement with the numerical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simplified Pore Network Model for Analysis of Flow Capacity and Residual Oil Saturation Distribution Based on Computerized Tomography

The pore network (connected tunnels) can be used not only to describe pore geometry inside a rock but also to calculate various fluid flow properties on a micro scale. Based on a “digital core” extracted from micro computerized tomography (μ-CT), we have developed a simplified pore network model. This study simplifies the pore structure as a lattice of cubes. Pore size, pore distribution, tunnel connectivity, fluid viscosity, interfacial tension, and the external driving pressure gradient are considered for each adjacent cube to evaluate the flow capacity of the pore network. We calculated the permeability, effective permeability, and effective sweep efficiency. The results show that fluid flow initially starts in the largest tunnel, which corresponds to the threshold pressure of the rock sample (fluid flow onset pressure). Permeability changes nonlinearly before all the tunnels are involved in the flow. Changes in the microscopic residual oil saturation and the sweep efficiency were visualized based on the three-dimensional pore distribution. Visualization shows that more oil remains in the formation, the heterogeneity is greater, and the oil viscosity is higher.

58 GEOSCIENCES↗

Nozzle tip wetting in gasoline direct injection injector and its link with nozzle internal flow

Fuel film in the gasoline direct injection injector tip, or so-called nozzle tip wetting, has been found to be an important contributor of particle emissions. Attempts have been made to reduce the nozzle tip wetting by optimizing nozzle geometry designs. However, the inherent mechanism of the nozzle tip wetting formation and its link with nozzle internal flow is still unclear yet due to the lack of direct observations. To overcome this insufficiency, here the nozzle internal flow and the formation process of the nozzle tip wetting were visualized in the real-scale aluminum nozzles using the X-ray phase-contrast technique. Results showed that the needle bouncing, injection pressure, and hole configuration affect the formation of the nozzle tip wetting, while the influence of needle bouncing is the most critical. A further study was conducted to examine the effect of nozzle counterbore diameter on the nozzle tip wetting. It was found that with an increase in counterbore diameter, the nozzle tip wetting slightly increased first and then decreased sharply after the counterbore diameter exceeded 0.40 mm. The mechanisms of the aforementioned phenomena were discussed in detail, which can contribute to the better understandings and control strategies of nozzle tip wetting.

42 ENGINEERING↗

Visualization of Simulation and TAP Interaction

The Temporal Analysis of Products (TAP) Reactor applies transient kinetics to analyze sophisticated chemical reaction mechanisms, rapidly changing temperature, flow rate, pressure or concentration and then monitoring the system’s time dependent response. A simulation tool exists for the TAP reactor, but it requires a visualization platform to interface effectively with researchers. This graphical user interface remedies this problem, fixing many of the difficulties caused by users interacting directly with the simulation code by providing an intuitive web application that satisfies the simulation input parameters for fast generation of data. This tool will ensure robust and repeatable simulation results, expedite the data generating process and enhance the overall ability of TAP to produce significant findings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Deep learning for in situ data compression of large turbulent flow simulations

As the size of turbulent flow simulations continues to grow, in situ data compression is becoming increasingly important for visualization, analysis, and restart checkpointing. For these applications, single-pass compression techniques with low computational and communication overhead are crucial. In this paper we present a deep-learning approach to in situ compression using an autoencoder architecture that is customized for three-dimensional turbulent flows and is well suited for contemporary heterogeneous computing resources. The autoencoder is compared against a recently introduced randomized single-pass singular value decomposition (SVD) for three different canonical turbulent flows: decaying homogeneous isotropic turbulence, a Taylor-Green vortex, and turbulent channel flow. Our proposed fully convolutional autoencoder architecture compresses turbulent flow snapshots by a factor of 64 with a single pass, allows for arbitrarily sized input fields, is cheaper to compute than the randomized single-pass SVD for typical simulation sizes, performs well on unseen flow configurations, and has been made publicly available. The results reported here show that the autoencoder dramatically outperforms a randomized single-pass SVD with similar compression ratio and yields comparable performance to a higher-rank decomposition with an order of magnitude less compression in regard to preserving a number of important statistical quantities such as turbulent kinetic energy, enstrophy, and Reynolds stresses.

97 MATHEMATICS AND COMPUTING↗