Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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,549 records · Page 86

Development and Implementation of 3-D, High Speed Capacitance Tomography for Imaging Large-Scale, Cold-Flow Circulating Fluidized Bed

A detailed understanding of multiphase flow behavior inside a Circulating Fluidized Bed (CFB) requires a 3-D technique capable of visualizing the flow field in real-time. Electrical Capacitance Volume Tomography (ECVT) is a newly developed technique that can provide such measurements. The attractiveness of the technique is in its low profile sensors, fast imaging speed and scalability to different section sizes, low operating cost, and safety. Moreover, the flexibility of ECVT sensors enable them to be designed around virtually any geometry, rendering them suitable to be used for measurement of solid flows in exit regions of the CFB. Tech4Imaging LLC has worked under contract with the U.S. Department of Energy's National Energy Technology Laboratory (DOE NETL) to develop an ECVT system for cold flow visualization and install it on a 12 inch ID circulating fluidized bed. The objective of this project was to help advance multi-phase flow science through implementation of an ECVT system on a cold flow model at DOE NETL. This project has responded to multi-phase community and industry needs of developing a tool that can be used to develop flow models, validate computational fluid dynamics simulations, provide detailed real-time feedback of process variables, and provide a comprehensive understating of multi-phase flow behavior. In this project, a complete ECVT system was successfully developed after considering different potential electronics and sensor designs. The system was tested at various flow conditions and with different materials, yielding real-time images of flow interaction in a gas-solid flow system. The system was installed on a 12 inch ID CFB of the US Department of Energy, Morgantown Labs. Technical and economic assessment of Scale-up and Commercialization of ECVT was also conducted. Experiments conducted with larger sensors in conditions similar to industrial settings are very promising. ECVT has also the potential to be developed for imaging multi-phase flow systems in high temperature and high pressure conditions, typical in many industrial applications.

42 ENGINEERING↗

he Impact of Primary Marine Aerosol on Atmospheric Chemistry, Radiation and Climate: A CCSM Model Development Study

This project examined the potential large-scale influence of marine aerosol cycling on atmospheric chemistry, physics and radiative transfer. Measurements indicate that the size-dependent generation of marine aerosols by wind waves at the ocean surface and the subsequent production and cycling of halogen-radicals are important but poorly constrained processes that influence climate regionally and globally. A reliable capacity to examine the role of marine aerosol in the global-scale atmospheric system requires that the important size-resolved chemical processes be treated explicitly. But the treatment of multiphase chemistry across the breadth of chemical scenarios encountered throughout the atmosphere is sensitive to the initial conditions and the precision of the solution method. This study examined this sensitivity, constrained it using high-resolution laboratory and field measurements, and deployed it in a coupled chemical-microphysical 3-D atmosphere model.

54 ENVIRONMENTAL SCIENCES↗

Physiochemical Evidence of Faulting Processes and Modeling of Fluid in Evolving Fault Systems in Southern California

Our study targets recent (Plio-Pleistocene) faults and young (Tertiary) petroleum fields in southern California. Faults include the Refugio Fault in the Transverse Ranges, the Ellwood Fault in the Santa Barbara Channel, and most recently the Newport- Inglewood in the Los Angeles Basin. Subsurface core and tubing scale samples, outcrop samples, well logs, reservoir properties, pore pressures, fluid compositions, and published structural-seismic sections have been used to characterize the tectonic/diagenetic history of the faults. As part of the effort to understand the diagenetic processes within these fault zones, we have studied analogous processes of rapid carbonate precipitation (scaling) in petroleum reservoir tubing and manmade tunnels. From this, we have identified geochemical signatures in carbonate that characterize rapid CO 2 degassing. These data provide constraints for finite element models that predict fluid pressures, multiphase flow patterns, rates and patterns of deformation, subsurface temperatures and heat flow, and geochemistry associated with large fault systems.

01 COAL, LIGNITE, AND PEAT↗

Computer Modeling and Analysis of Truck Generated Salt-Spray Transport Near Bridges

Bridge painting is very expensive. Over the lifetime of a bridge, painting to protect the steel from the elements can cost a large fraction of the original construction cost. In addition to monetary costs, maintenance workers face health risks if strict procedures are not in place to avoid exposure to hazardous dust when old paint is removed and the fumes of hazardous chemicals when new paint is applied. Use of weathering steel in bridge construction avoids the large cost of periodic painting and the hazards of paint removal and repainting because this type of steel is never painted. When exposed to weather, the initial corrosion of weathering steel forms a natural corrosion resistant layer that prevents further corrosion under a wide variety of conditions. Given the advantages of using weathering steel, the choice might appear to be obvious. However, the large potential savings in bridge maintenance costs in using weathering steel are not always realized because the protective patina layer fails to adhere when the steel is exposed to excessively salt-laden and moisture-laden environments. A multi-year study is underway by the FHWA to better quantify the conditions that lead to excessive corrosion in weathering steel bridges. This report documents work using computational fluid dynamics (CFD) analysis to study the conditions and mechanisms that lead from the salt spray thrown from truck tires to salt water droplets reaching bridge girders. Computer simulations use the advanced motion modeling capability called a “sliding mesh” in the CFD software to move one or more large trucks under bridges at 60 mph. The computations include multiphase spray droplet tracking of salt-laden water droplets coming off the tires.

Lottes, S. A.↗

Crystalline Ceramic Waste Forms: Comparison Of Reference Process For Ceramic Waste Form Fabrication

The research conducted in this work package is aimed at taking advantage of the long term thermodynamic stability of crystalline ceramics to create more durable waste forms (as compared to high level waste glass) in order to reduce the reliance on engineered and natural barrier systems. Durable ceramic waste forms that incorporate a wide range of radionuclides have the potential to broaden the available disposal options and to lower the storage and disposal costs associated with advanced fuel cycles. Assemblages of several titanate phases have been successfully demonstrated to incorporate radioactive waste elements, and the multiphase nature of these materials allows them to accommodate variation in the waste composition. Recent work has shown that they can be produced from a melting and crystallization process. The objective of this report is to explore the phase formation and microstructural differences between lab scale melt processing in varying gas environments with alternative densification processes such as Hot Pressing (HP) and Spark Plasma Sintering (SPS). The waste stream used as the basis for the development and testing is a simulant derived from a combination of the projected Cs/Sr separated stream, the Trivalent Actinide - Lanthanide Separation by Phosphorous reagent Extraction from Aqueous Komplexes (TALSPEAK) waste stream consisting of lanthanide fission products, the transition metal fission product waste stream resulting from the transuranic extraction (TRUEX) process, and a high molybdenum concentration with relatively low noble metal concentrations. Melt processing as well as solid state sintering routes SPS and HP demonstrated the formation of the targeted phases; however differences in microstructure and elemental partitioning were observed. In SPS and HP samples, hollandite, pervoskite/pyrochlore, zirconolite, metallic alloy and TiO 2 and Al 2 O 3 were observed distributed in a network of fine grains with small residual pores. The titanate phases that incorporate M +3 rare earth elements were observed to be distinct phases (ex. Nd 2 Ti 2 O 7 with less degree of substitution as compared to the more homogeneous melt processed samples where a high degree of substitution and variation of composition within grains was observed. Liquid phase sintering was enhanced in reducing gas environments and resulted in large (10-200 microns) irregular shaped grains along with large voids associated with the melt process; SPS and HP samples exhibited finer grain size with smaller voids. Metallic alloys were observed in the bulk of the sample for SPS and HP samples, but were found at the bottom of the crucible in melt processed trials. These results indicate that for a first melter trial, the targeted phases can be formed in air by utilizing Ti/TiO 2 additives which aid phase formation and improve the electrical conductivity. Ultimately, a melter run in reducing gas environments would be beneficial to study differences in phase formation and elemental partitioning.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Continued Need for Modeling and Scaled Testing to Advance the Hanford Tank Waste Mission

Hanford tank wastes are chemically complex slurries of liquids and solids that can exhibit changes in rheological behavior during retrieval and processing. The Hanford Waste Treatment and Immobilization Plant (WTP) recently abandoned its planned approach to use computational fluid dynamics (CFD) supported by testing at less than full scale to verify the design of vessels that process these wastes within the plant. The commercial CFD tool selected was deemed too difficult to validate to the degree necessary for use in the design of a nuclear facility. Alternative, but somewhat immature, CFD tools are available that can simulate multiphase flow of non-Newtonian fluids. Yet both CFD and scaled testing can play an important role in advancing the Hanford tank waste mission—in supporting the new verification approach, which is to conduct testing in actual plant vessels; in supporting waste feed delivery, where scaled testing is ongoing; as a fallback approach to design verification if the Full Scale Vessel Testing Program is deemed too costly and time-consuming; to troubleshoot problems during commissioning and operation of the plant; and to evaluate the effects of any proposed changes in operating conditions in the future to optimize plant performance.

Fluid Dynamics↗

FY10 Report on Multi-scale Simulation of Solvent Extraction Processes: Molecular-scale and Continuum-scale Studies

This task is aimed at predictive modeling of solvent extraction processes in typical extraction equipment through multiple simulation methods at various scales of resolution. We have conducted detailed continuum fluid dynamics simulation on the process unit level as well as simulations of the molecular-level physical interactions which govern extraction chemistry. Through combination of information gained through simulations at each of these two tiers along with advanced techniques such as the Lattice Boltzmann Method (LBM) which can bridge these two scales, we can develop the tools to work towards predictive simulation for solvent extraction on the equipment scale (Figure 1). The goal of such a tool-along with enabling optimized design and operation of extraction units-would be to allow prediction of stage extraction effrciency under specified conditions. Simulation efforts on each of the two scales will be described below. As the initial application of FELBM in the work performed during FYl0 has been on annular mixing it will be discussed in context of the continuum-scale. In the future, however, it is anticipated that the real value of FELBM will be in its use as a tool for sub-grid model development through highly refined DNS-like multiphase simulations facilitating exploration and development of droplet models including breakup and coalescence which will be needed for the large-scale simulations where droplet level physics cannot be resolved. In this area, it can have a significant advantage over traditional CFD methods as its high computational efficiency allows exploration of significantly greater physical detail especially as computational resources increase in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Developments in digital in-line holography enable validated measurement of 3D particle field dynamics

Digital in-line holography is an optical technique which can be applied to measure the size, three-dimensional position, and three-component velocity of disperse particle fields. This work summarizes recent developments at Sandia National Laboratories focused on improvement in measurement accuracy, experimental validation, and applications to multiphase flows. New routines are presented which reduce the uncertainty in measured position along the optical axis to a fraction of the particle diameter. Furthermore, application to liquid atomization highlights the ability to measure complex, three-dimensional structures. Finally, investigation of particles traveling at near sonic conditions prove accuracy despite significant experimental noise due to shock-waves.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Annual Report: Carbon Capture Simulation Initiative (CCSI) (30 September 2013)

The Carbon Capture Simulation Initiative (CCSI) is a partnership among national laboratories, industry and academic institutions that is developing and deploying state-of-the-art computational modeling and simulation tools to accelerate the commercialization of carbon capture technologies from discovery to development, demonstration, and ultimately the widespread deployment to hundreds of power plants. The CCSI Toolset will provide end users in industry with a comprehensive, integrated suite of scientifically validated models, with uncertainty quantification (UQ), optimization, risk analysis and decision making capabilities. The CCSI Toolset incorporates commercial and open-source software currently in use by industry and is also developing new software tools as necessary to fill technology gaps identified during execution of the project. Ultimately, the CCSI Toolset will (1) enable promising concepts to be more quickly identified through rapid computational screening of devices and processes; (2) reduce the time to design and troubleshoot new devices and processes; (3) quantify the technical risk in taking technology from laboratory-scale to commercial-scale; and (4) stabilize deployment costs more quickly by replacing some of the physical operational tests with virtual power plant simulations. CCSI is led by the National Energy Technology Laboratory (NETL) and leverages the Department of Energy (DOE) national laboratories’ core strengths in modeling and simulation, bringing together the best capabilities at NETL, Los Alamos National Laboratory (LANL), Lawrence Berkeley National Laboratory (LBNL), Lawrence Livermore National Laboratory (LLNL), and Pacific Northwest National Laboratory (PNNL). The CCSI’s industrial partners provide representation from the power generation industry, equipment manufacturers, technology providers and engineering and construction firms. The CCSI’s academic participants (Carnegie Mellon University, Princeton University, West Virginia University, Boston University and the University of Texas at Austin) bring unparalleled expertise in multiphase flow reactors, combustion, process synthesis and optimization, planning and scheduling, and process control techniques for energy processes. During Fiscal Year (FY) 13, CCSI announced the initial release of its first set of computational tools and models during the October 2012 meeting of its Industry Advisory Board. This initial release led to five companies licensing the CCSI Toolset under a Test and Evaluation Agreement this year. By the end of FY13, the CCSI Technical Team had completed development of an updated suite of computational tools and models. The list below summarizes the new and enhanced toolset components that were released following comprehensive testing during October 2013. 1. FOQUS. Framework for Optimization and Quantification of Uncertainty and Sensitivity. Package includes: FOQUS Graphic User Interface (GUI), simulation-based optimization engine, Turbine Client, and heat integration capabilities. There is also an updated simulation interface and new configuration GUI for connecting Aspen Plus or Aspen Custom Modeler (ACM) simulations to FOQUS and the Turbine Science Gateway. 2. A new MFIX-based Computational Fluid Dynamics (CFD) model to predict particle attrition. 3. A new dynamic reduced model (RM) builder, which generates computationally efficient RMs of the behavior of a dynamic system. 4. A completely re-written version of the algebraic surrogate model builder for optimization (ALAMO). The new version is several orders of magnitude faster than the initial release and eliminates the MATLAB dependency. 5. A new suite of high resolution filtered models for the hydrodynamics associated with horizontal cylindrical objects in a flow path. 6. The new Turbine Science Gateway (Cluster), which supports FOQUS for running multiple simulations for optimization or UQ using a local computer or cluster. 7. A new statistical tool (BSS-ANOVA-UQ) for calibration and validation of CFD models. 8. A new basic data submodel in Aspen Plus format for a representative high viscosity capture solvent, 2-MPZ system. 9. An updated RM tool for CFD (REVEAL) that can create a RM from MFIX. A new lightweight, stand-alone version will be available in late 2013. 10. An updated RM integration tool to convert the RM from REVEAL into a CAPE-OPEN or ACM model for use in a process simulator. 11. An updated suite of unified steady-state and dynamic process models for solid sorbent carbon capture included bubbling fluidized bed and moving bed reactors. 12. An updated and unified set of compressor models including steady-state design point model and dynamic model with surge detection. 13. A new framework for the synthesis and optimization of coal oxycombustion power plants using advanced optimization algorithms. This release focuses on modeling and optimization of a cryogenic air separation unit (ASU). 14. A new technical risk model in spreadsheet format. 15. An updated version of the sorbent kinetic/equilibrium model for parameter estimation for the 1st generation sorbent model. 16. An updated process synthesis superstructure model to determine optimal process configurations utilizing surrogate models from ALAMO for adsorption and regeneration in a solid sorbent process. 17. Validation models for NETL Carbon Capture Unit utilizing sorbent AX. Additional validation models will be available for sorbent 32D in 2014. 18. An updated hollow fiber membrane model and system example for carbon capture. 19. An updated reference power plant model in Thermoflex that includes additional steam extraction and reinjection points to enable heat integration module. 20. An updated financial risk model in spreadsheet format.

01 COAL, LIGNITE, AND PEAT↗

Gasification Studies: Task 4 Topical Report, Utah Clean Coal Program

A key objective of the Task 4 activities has been to develop simulation tools to support development, troubleshooting and optimization of pressurized entrained-flow coal gasifiers. The overall gasifier models (Subtask 4.1) combine submodels for fluid flow (Subtask 4.2) and heat transfer (Subtask 4.3) with fundamental understanding of the chemical processes (Subtask 4.4) processes that take place as coal particles are converted to synthesis gas and slag. However, it is important to be able to compare predictions from the models against data obtained from actual operating coal gasifiers, and Subtask 4.6 aims to provide an accessible, non-proprietary system, which can be operated over a wide range of conditions to provide well-characterized data for model validation. Highlights of this work include: (1) Verification and validation activities performed with the Arches coal gasification simulation tool on experimental data from the CANMET gasifier (Subtask 4.1). (2) The simulation of multiphase reacting flows with coal particles including detailed gas-phase chemistry calculations using an extension of the one-dimensional turbulence model’s capability (Subtask 4.2). (3) The demonstration and implementation of the Reverse Monte Carlo ray tracing (RMCRT) radiation algorithm in the ARCHES code (Subtask 4.3). (4) Determination of steam and CO 2 gasification kinetics of bituminous coal chars at high temperature and elevated pressure under entrained-flow conditions (Subtask 4.4). In addition, attempts were made to gain insight into the chemical structure differences between young and mature coal soot, but both NMR and TEM characterization efforts were hampered by the highly reacted nature of the soot. (5) The development, operation, and demonstration of in-situ gas phase measurements from the University of Utah’s pilot-scale entrained-flow coal gasifier (EFG) (Subtask 4.6). This subtask aimed at acquiring predictable, consistent performance and characterizing the environment within the gasifier.

01 COAL, LIGNITE, AND PEAT↗

SESAME 96171, a three-phase equation of state for CeO2

We describe an equation of state (EOS) for cerium (IV) oxide, CeO 2 , that includes two solid phases and the liquid. The models and parameters used to calculate the EOS are presented in detail, and we compare with data for both crystal phases and ambient melt. We discuss complications that arise when making multiphase EOS.

36 MATERIALS SCIENCE↗

Modeling and Field Test Planning Activities in Support of Disposal of Heat-Generating Waste in Salt

The modeling efforts in support of the field test planning conducted at LBNL leverage on recent developments of tools for modeling coupled thermal-hydrological-mechanical-chemical (THMC) processes in salt and their effect on brine migration at high temperatures. This work includes development related to, and implementation of, essential capabilities, as well as testing the model against relevant information and published experimental data related to the fate and transport of water. These are modeling capabilities that will be suitable for assisting in the design of field experiment, especially related to multiphase flow processes coupled with mechanical deformations, at high temperature. In this report, we first examine previous generic repository modeling results, focusing on the first 20 years to investigate the expected evolution of the different processes that could be monitored in a full-scale heater experiment, and then present new results from ongoing modeling of the Thermal Simulation for Drift Emplacement (TSDE) experiment, a heater experiment on the in-drift emplacement concept at the Asse Mine, Germany, and provide an update on the ongoing model developments for modeling brine migration. LBNL also supported field test planning activities via contributions to and technical review of framework documents and test plans, as well as participation in workshops associated with field test planning.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SRNL CRP progress report [Development of Melt Processed Ceramics for Nuclear Waste Immobilization]

A multi-phase ceramic waste form is being developed at the Savannah River National Laboratory (SRNL) for treatment of secondary waste streams generated by reprocessing commercial spent nuclear. The envisioned waste stream contains a mixture of transition, alkali, alkaline earth, and lanthanide metals. Ceramic waste forms are tailored (engineered) to incorporate waste components as part of their crystal structure based on knowledge from naturally found minerals containing radioactive and non-radioactive species similar to the radionuclides of concern in wastes from fuel reprocessing. The ability to tailor ceramics to mimic naturally occurring crystals substantiates the long term stability of such crystals (ceramics) over geologic timescales of interest for nuclear waste immobilization [1]. A durable multiphase ceramic waste form tailored to incorporate all the waste components has the potential to broaden the available disposal options and thus minimize the storage and disposal costs associated with aqueous reprocessing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation

A multiscale linear-solver framework for the pressure equation associated with flow in highly heterogeneous porous formations was developed. The multiscale based approach is cast in a general algebraic form, which facilitates integration of the new scalable linear solver in existing flow simulators. The Algebraic Multiscale Solver (AMS) is employed as a preconditioner within a multi-stage strategy. The formulations investigated include the standard MultiScale Finite-Element (MSFE) andMultiScale Finite-Volume (MSFV) methods. The local-stage solvers include incomplete factorization and the so-called Correction Functions (CF) associated with the MSFV approach. Extensive testing of AMS, as an iterative linear solver, indicate excellent convergence rates and computational scalability. AMS compares favorably with advanced Algebraic MultiGrid (AMG) solvers for highly detailed three-dimensional heterogeneous models. Moreover, AMS is expected to be especially beneficial in solving time-dependent problems of coupled multiphase flow and transport in large-scale subsurface formations.

02 PETROLEUM↗

User's Manual for the Hydrate v1.5 Option of TOUGH+ v1.5: A Code for the Simulation of System Behavior in Hydrate-Bearing Geologic Media

HYDRATE v1.5 is a numerical code that for the simulation of the behavior of hydrate-bearing geologic systems, and represents the third update of the code since its first release [Moridis et al., 2008]. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. HYDRATE v1.5 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. By solving the coupled equations of mass and heat balance, the fully operational TOUGH+HYDRATE code can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH4-hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy’s law is valid. TOUGH+HYDRATE v1.5 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH4, hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.

54 ENVIRONMENTAL SCIENCES↗

User's Manual for the RealGasBrine v1.0 Option of TOUGH+ v1.5: A Code for the Simulation of System Behavior in Gas-Bearing Geologic Media

REALGASBRINE v1.0 is a numerical code that for the simulation of the behavior of gas-bearing porous and/fractured geologic media. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. REALGASBRINE v1.0 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. REALGASBRINE v1.0 describes the non-isothermal two- (for pure water) or three-phase (for brine) flow of an aqueous phase and a real gas mixture in a gas-bearing medium, with a particular focus in ultra-tight (such as tight-sand and shale gas) systems. Up to 12 individual real gases can be tracked, and salt can precipitate as solid halite. The capabilities of the code include coupled flow and thermal effects, real gas behavior, Darcy and non-Darcy flow, several isotherm options of gas sorption onto the grains of the porous media, complex fracture descriptions, gas solubility into water, and geomechanical effects on flow properties. REALGASBRINE v1.0 allows the study of flow and transport of fluids and heat over a wide range of time frames and spatial scales not only in gas reservoirs, but also in any problem involving the flow of gases in geologic media, including the geologic storage of greenhouse gas mixtures, the behavior of geothermal reservoirs with multi-component condensable (H 2 O and CO 2 ) and non-condensable gas mixtures, the transport of water and released H 2 in nuclear waste storage applications, etc.

15 GEOTHERMAL ENERGY↗

User's Manual of the TOUGH+ Core Code v1.5: A General-Purpose Simulator of Non-Isothermal Flow and Transport through Porous and Fractured Media

TOUGH+ v1.5 is a numerical code for the simulation of multi-phase, multi-component flow and transport of mass and heat through porous and fractured media, and represents the third update of the code since its first release [Moridis et al., 2008]. TOUGH+ is a successor to the TOUGH2 [Pruess et al., 1991] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstations, PC, Macintosh). TOUGH+ v1.5 employs dynamic memory allocation, thus minimizing storage requirements. It has a completely modular structure, follows the tenets of Object-Oriented Programming (OOP), and involves the advanced features of FORTRAN 95/2003, i.e., modules, derived data types, the use of pointers, lists and trees, data encapsulation, defined operators and assignments, operator extension and overloading, use of generic procedures, and maximum use of the powerful intrinsic vector and matrix processing operations. This report presents the core TOUGH+ v1.5 code, i.e., the part of the code that is common to all its applications. It provides a description of the underlying physics and thermodynamics of non-isothermal flow, of the mathematical and numerical approaches, as well as a detailed explanation of the general (common to all applications) input requirements, options, capabilities and output specifications. The core code cannot run by itself: it needs to be coupled with the code for the specific TOUGH+ application option that describes a particular type of problem. The additional input requirements specific to a particular TOUGH+ application options and related illustrative examples can be found in the corresponding User’s Manual.

54 ENVIRONMENTAL SCIENCES↗

Area of Interest 1, CO 2 at the Interface: Nature and Dynamics of the Reservoir/Caprock Contact and Implications for Carbon Storage Performance

We examined the influence of geologic features present at the reservoir/caprock interface on the transmission of supercritical CO 2 into and through caprock. We focused on the case of deformation-band faults in reservoir lithologies that intersect the interface and transition to opening-mode fractures in caprock lithologies. Deformation-band faults are exceeding common in potential CO 2 injection units and our fieldwork in Utah indicates that this sort of transition is common. To quantify the impact of these interface features on flow and transport we first described the sedimentology and permeability characteristics of selected sites along the Navajo Sandstone (reservoir lithology) and Carmel Formation (caprock lithology) interface, and along the Slickrock Member (reservoir lithology) and Earthy Member (caprock lithology) of the Entrada Sandstone interface, and used this information to construct conceptual permeability models for numerical analysis. We then examined the impact of these structures on flow using single-phase and multiphase numerical flow models for these study sites. Key findings include: (1) Deformation-band faults strongly compartmentalize the reservoir and largely block cross-fault flow of supercritical CO 2 . (2) Significant flow of CO 2 through the fractures is possible, however, the magnitude is dependent on the small-scale geometry of the contact between the opening-mode fracture and the deformation band fault. (3) Due to the presence of permeable units in the caprock, caprock units are capable of storing significant volumes of CO 2 , particularly when the fracture network does not extend all the way through the caprock. The large-scale distribution of these deformation-bandfault-to-opening-mode-fractures is related to the curvature of the beds, with greater densities of fractures in high curvature regions. We also examined core and outcrops from the Mount Simon Sandstone and Eau Claire Formation reservoir/caprock interface in order to extend our work to a reservoir/caprock pair this is currently being assessed for long-term carbon storage. These analyses indicate that interface features similar to those observed at the Utah sites 3 were not observed. Although not directly related to our main study topic, one byproduct of our investigation is documentation of exceptionally high degrees of heterogeneity in the pore-size distribution of the Mount Simon Sandstone. This suggests that the unit has a greater-than-normal potential for residual trapping of supercritical CO 2 .

54 ENVIRONMENTAL SCIENCES↗