Computational Modeling of Fluid and Particle Flows in Rock Fractures
Presented by NETL at 2016 NETL Workshop on Multiphase Flow Science; 8/16
Computational fluid dynamics↗
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.
Impact of Variable Gas Mixtures on Bubble Size Distribution and Mass Transfer in Gas Fermentation Reactors
Gas fermentation has emerged as a promising new technology for the generation of fuels and chemicals from mixtures of greenhouse and energy rich gas streams (CO2/CH4/H2/CO) via microbial bioreaction. Example pathways include biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Effective Gas-liquid mass-transfer is an important physical phenomenon that determines the design and scale-up of these systems. There is currently a knowledge-gap regarding bubble-size distributions when using a mixture of gases with vastly different properties, which can have a significant impact on overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a detailed multiphase computational-fluid-dynamics (CFD) model to study the impact of gas-mixtures on overall mass-transfer in bubble column and air-lift reactors. The CFD tool previously developed by the authors (1) for simulating aerobic fermentation reactors at scale is used in this study. The Reynolds-averaged mass, momentum, energy, and species transport equations are solved for interpenetrating gas and liquid phase in this model. We use a population balance-based bubble-size-distribution model that is validated against small-scale experiments in our solver. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient is analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors.
BIOMASS FUELS,MATHEMATICS AND COMPUTING↗
Performance Assessment of a Generic Nuclear Waste Repository in Shale - 18392
A performance assessment (PA) is an important component of a comprehensive safety analysis for a nuclear waste repository. In a PA, probabilistic simulations of the total repository system are performed, and results are evaluated against performance metrics. Uncertainty and sensitivity analyses help prioritize additional research and model development. The United States Department of Energy (DOE) has been developing a state of the art PA simulation software tool-kit, Geologic Disposal Safety Assessment (GDSA) Framework, that couples increasingly higher fidelity models of subsystem processes into total system PA simulations. Over the past several years, PAs of generic repositories in three geologic media (salt, shale, and crystalline rock) have demonstrated ongoing developments in capability. The current PA of a nuclear waste repository in a generic shale formation showcases GDSA Framework, including capabilities in domain discretization (Cubit), multi-physics simulations (PFLOTRAN), uncertainty and sensitivity analysis (Dakota), and visualization (Paraview). The generic shale reference case considers the disposal of 22,000 metric tons heavy metal of commercial spent nuclear fuel (SNF) in a generic shale formation. PA simulations account for the thermal load and radionuclide inventory of the waste form, components of the engineered barrier system including waste package and buffer, and components of the natural barrier system including the host rock shale and underlying and overlying stratigraphic units. Two repository layouts are considered, one for emplacement of waste packages containing 12 pressurized water reactor (PWR) assemblies, and one for 4-PWR waste package emplacement. Model domains are half-symmetry and contain between 7 and 22 million grid cells. Grid refinement captures the detail of individual waste packages, emplacement drifts, access drifts, and shafts. Simulations are run in a high-performance computing (HPC) environment on 512 or 2048 processes. The governing equations describing coupled heat and fluid flow and reactive transport are solved with PFLOTRAN, an open-source, massively parallel multiphase flow and reactive transport code. Additional simulated processes include waste package degradation; waste form dissolution; radioactive decay and ingrowth in aqueous, solid, and adsorbed phases; and a simple biosphere based on a pumping well. Simulations are run to 106 y; dose and radionuclide concentrations are observed within aquifers at a point approximately 5 km downgradient of the repository. Dakota is used to sample likely ranges of input parameters including waste form and waste package degradation rates and properties of engineered and natural materials to quantify uncertainty in predicted concentrations and sensitivity to input parameters. Given the assumptions of the reference case, results of 12-PWR and 4-PWR simulations are very similar. Dose at the pumping well 5-km downgradient of the repository does not exceed 10-9 Sv/y. I-129 is the largest contributor to dose. The dose due to Cl-36, the next largest contributor, is orders of magnitude smaller. At the pumping well, the concentration of I-129 5-km downgradient does not exceed 10-14 mol/L. Rank correlation coefficients indicate that the maximum concentration of I-129 within the aquifer is sensitive to shale (repository host rock) porosity and aquifer permeability. These results affirm that HPC-capable codes can be used to simulate important multi-physics couplings directly in a total system performance assessment of a geologic repository, and can be used in prioritization of future research and model development. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
LDRD 22A1059-068FP Tailoring the Properties of Multi-Phase Materials Through the Use of Correlative Microscopy and Machine Learning - Poster
High strength alloys with good ductility, hardness, and toughness are needed to meet stringent design requirements for extreme environments. One complication in this pursuit is the evidence that metals rarely exhibit both high strength and good fracture toughness as the underlying mechanisms work in opposition. An exception to this behavior is found in multiphase alloys that form complex microstructures of mixed phases with variable grain sizes and shapes that provide increased fracture toughness by the arrangement of their constituent elements. We propose to explore this phenomenon using state-of-the-art machine learning (ML) techniques in a new and novel manner to identify and correlate the critical microstructural features in a Titanium-10Vanadium-2Iron-3Aluminum (Ti-10V-2Fe-3Al) alloy that is reported to exhibit high strength and fracture toughness. Additionally, we will employ multiple, complementary characterization techniques such as optical microscopy, electron backscatter diffraction (EBSD), energy dispersive spectroscopy (EDS) and scanning electron microscopy to provide multi-layer, quantitative ground truth measures of the microstructures. This data will be used to train a Convolutional Neural Network (CNN) in a semi-supervised environment to identify key microstructural features such as ? platelet dimensions and locations and ?/? phase boundaries and correlate those features with the strength and toughness. Here the ? and ? nomenclature refers to hexagonal close pack (hcp) and body center cubic (bcc) crystal structures, respectively. Previous work has focused on popular alloys and typically used one characterization technique. This research is focused on a promising titanium alloy, uses multiple complimentary characterization tools to provide precise microstructural information and correlates to improved fracture toughness. The resulting ML tool can be trained for additional microstructural features, different alloy(s), and or target mechanical properties.
36 MATERIALS SCIENCE↗
Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation
Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Simulation of Boiling Flow in a Helical Coil Using a Homogeneous Equilibrium Model with the Nek-2P CFD Code
Recent developments in the field of nuclear energy have lead to an increased interest in the design of helically coiled steam generators. In particular, their compact size, and high heat transfer rates make them ideal for small modular reactors. This has lead to a renewed effort to investigate the characteristics of fluid flow and heat transfer in and around helical coils both numerically and experimentally. The focus of this work is on the development of a consistent modeling, formulation, and simulation framework for boiling flows in helical coils using the highly scalable Nek-2P CFD code. As a first step, two phase flow is accounted for using a homogeneous equilibrium model. Results for both single and multiphase flows are compared to results obtained from Star-CCM+. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Simulating Bubble Deformations Using the Intersection Marker (ISM) Front-Tracking method
Surface and volume tracking methods have been under continuous development over the last 50+ years since the first MAC (Marker and Cell) and SLIC (Simple Line Interface Calculation) methods appeared in the 1960's. The sustained interest in simulating multiphase flow stem from their many practical applications in the physical world; from biological systems, to chemical systems, to boiling dynamics and atmospheric heat and mass transfer. The ability to capture the details of interface evolution as well as the physics of interfacial heat and mass transfer has thus been a tempting goal for many numerical researchers. Numerical analysts, apart from the well-established and popular VOF PLIC (Volume of Fluid, Piecewise Linear Interface Construction) method, can now call upon even more computationally intensive options such as the LS (Level Set) method, the FT (Front Tracking) method, the MoF (Moment of Fluid) method and hybrid methods.
42 ENGINEERING↗
Numerical Analysis of Water Harvesting from Saturated Vapor with the Electrospray Water Droplet Injection in The Solar-Driven Desalination Applications
Solar-driven seawater desalination systems provide freshwater through environmentally friendly and carbon-neutral processes. Air humidification and dehumidification desalination (HDD) systems extract water directly from moist air. This paper focuses on a new solar-driven seawater desalination system that shares some common traits with HDD systems but in which the air is eliminated from the process. Seawater is vaporized by solar thermal radiation in high-performance solar panels. Water vapor flows to a new electro-condensation chamber and passes through a series of electrosprays, which inject small nuclei of freshwater droplets into the vapor clouds. Here, due to dielectrophoresis and electrohydrodynamic flows, the vapor molecules are captured by the charged droplets, and vapor is condensed at the droplet surface. The water vapor condensation was modeled by using the multiphase flow numerical model. The CFD model, which was also experimentally validated, used SprayFoam solver modules available in the Open FOAM® CFD open-source software. In preliminary results, 802 grams/(hr-m2) (grams per hour per unit area) of water was harvested directly from moist air at 24C and 90% relative humidity. The simulation results indicated that freshwater productivity could increase up to 4,000 grams/(hr-m2) for fully saturated air conditions. This represents about 5 times the water productivity rate of current conventional HDD systems today. In addition, energy consumption could be reduced by 80 percent compared to conventional HDD systems.
Electrospray, electro-condensation, solar, droplet↗
Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification
The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions
12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN
To better inform the subsurface scientist on the expected performance of parallel simulators, this work investigates performance of the reactive multiphase flow and multicomponent biogeochemical transport code PFLOTRAN as it is applied to several realistic modeling scenarios run on the Jaguar supercomputer. After a brief introduction to the code's parallel layout and code design, PFLOTRAN's parallel performance (measured through strong and weak scalability analyses) is evaluated in the context of conceptual model layout, software and algorithmic design, and known hardware limitations. PFLOTRAN scales well (with regard to strong scaling) for three realistic problem scenarios: (1) in situ leaching of copper from a mineral ore deposit within a 5-spot flow regime, (2) transient flow and solute transport within a regional doublet, and (3) a real-world problem involving uranium surface complexation within a heterogeneous and extremely dynamic variably saturated flow field. Weak scalability is discussed in detail for the regional doublet problem, and several difficulties with its interpretation are noted.
54 ENVIRONMENTAL SCIENCES↗
A harmonic pulse testing method for leakage detection in deep subsurface storage formations
Abstract Detection of leakage in deep geologic storage formations (e.g., carbon sequestration sites) is a challenging problem. This study investigates an easy‐to‐implement frequency domain leakage detection technology based on harmonic pulse testing (HPT). Unlike conventional constant‐rate pressure interference tests, HPT stimulates a reservoir using periodic injection rates. The fundamental principle underlying HPT‐based leakage detection is that leakage modifies a storage system's frequency response function, thus providing clues of system malfunction. During operations, routine HPTs can be conducted at multiple pulsing frequencies to obtain experimental frequency response functions, using which the possible time‐lapse changes are examined. In this work, a set of analytical frequency response solutions is derived for predicting system responses with and without leaks for single‐phase flow systems. Sensitivity studies show that HPT can effectively reveal the presence of leaks. A search procedure is then prescribed for locating the actual leaks using amplitude and phase information obtained from HPT, and the resulting optimization problem is solved using the genetic algorithm. For multiphase flows, the applicability of HPT‐based leakage detection procedure is exemplified numerically using a carbon sequestration problem. Results show that the detection procedure is applicable if the average reservoir conditions in the testing zone stay relatively constant during the tests, which is a working assumption under many other interpretation methods for pressure interference tests. HPT is a cost‐effective tool that only requires periodic modification of the nominal injection rate. Thus it can be incorporated into existing monitoring plans with little additional investment.
Sun, Alexander Y.↗
Transient seafloor venting on continental slopes from warming‐induced methane hydrate dissociation
Methane held in frozen hydrate cages within marine sediment comprises one of the largest carbon reservoirs on the planet. Recent submarine observations of widespread methane seepage may record hydrate dissociation due to oceanic warming, which consequently may further amplify climate change. Here, in this work, we simulate the effect of seafloor warming on marine hydrate deposits using a multiphase flow model. We show that hydrate dissociation, gas migration, and subsequent hydrate formation can generate temporary methane venting into the ocean through the hydrate stability zone. Methane seeps venting through the hydrate stability zone on the eastern Atlantic margin may record this process due to warming begun thousands of years ago. Our results contrast with the traditional view that venting occurs only updip of the hydrate stability zone.
hydrate↗
Quantifying hydrate solidification front advancing using method of characteristics
Here, we develop a one-dimensional analytical solution based on the method of characteristics to explore hydrate formation from gas injection into brine-saturated sediments within the hydrate stability zone. Our solution includes fully coupled multiphase and multicomponent flow and the associated advective transport in a homogeneous system. Our solution shows that hydrate saturation is controlled by the initial thermodynamic state of the system and changed by the gas fractional flow. Hydrate saturation in gas-rich systems can be estimated by 1 – cl 0 /cl e when Darcy flow dominates, where cl 0 is the initial mass fraction of salt in brine, and cl e is the mass fraction of salt in brine at three-phase (gas, liquid, and hydrate) equilibrium. Hydrate saturation is constant, gas saturation and gas flux decrease, and liquid saturation and liquid flux increase with the distance from the gas inlet to the hydrate solidification front. The total gas and liquid flux is constant from the gas inlet to the hydrate solidification front and decreases abruptly at the hydrate solidification front due to gas inclusion into the hydrate phase. The advancing velocity of the hydrate solidification front decreases with hydrate saturation at a fixed gas inflow rate. This analytical solution illuminates how hydrate is formed by gas injection (methane, CO 2 , ethane, propane) at both the laboratory and field scales.
58 GEOSCIENCES↗
Experimental analysis of spatial correlation effects on capillary trapping of supercritical CO 2 at the intermediate laboratory scale in heterogeneous porous media
Several numerical studies have demonstrated that the heterogeneous nature of typical sedimentary formations can favorably dampen the accumulation of mobile CO 2 phase underneath the caprock. Core flooding experiments have also shown that contrasts in capillary entry pressure can lead to buildup of nonwetting fluid phase (NWP) at interfaces between facies. Explicit representation of geological heterogeneity at the intermediate (cm-to-m) scale is a powerful approach to identify the key mechanisms that control multiphase flow dynamics in porous media. The ability to carefully control flow regime and permeability contrast at a scale that is relevant to CO 2 plume dynamics in saline formations offers valuable information to understand immiscible displacement processes and provides a benchmark for mathematical models. To provide insight into the impact of capillary heterogeneity on flow dynamics and trapping efficiency of supercritical CO 2 under successive drainage and imbibition conditions, we present an experimental investigation conducted in a synthetic sand reservoir. By mimicking the interplay of governing forces at reservoir conditions via application of surrogate fluids, we performed three immiscible displacement experiments to observe the entrapment of NWP in heterogeneous porous media. Capillary trapping performance is evaluated for each scenario through spatial and temporal variations of NWP saturation; for this reason we adopted X-ray attenuation to precisely measure phase saturation throughout the flow domain and apply spatial moment analysis. Furthermore, the sweeping performance of two different permeability fields with comparable variance but distinct spatial correlation was compared against a homogeneous base case with equivalent mean permeability by means of spatial moment analysis.
58 GEOSCIENCES↗
Cloud droplet activation through oxidation of organic aerosol influenced by temperature and particle phase state: CLOUD ACTIVATION BY AGED ORGANIC AEROSOL
Chemical aging of organic aerosol (OA) through multiphase oxidation reactions can alter their cloud condensation nuclei (CCN) activity and hygroscopicity. However, the oxidation kinetics and OA reactivity depend strongly on the particle phase state, potentially influencing the hydrophobic-to-hydrophilic conversion rate of carbonaceous aerosol. Here, amorphous Suwannee River fulvic acid (SRFA) aerosol particles, a surrogate humic-like substance (HULIS) that contributes substantially to global OA mass, are oxidized by OH radicals at different temperatures and phase states. When oxidized at low temperature in a glassy solid state, the hygroscopicity of SRFA particles increased by almost a factor of two, whereas oxidation of liquid-like SRFA particles at higher temperatures did not affect CCN activity. Low-temperature oxidation appears to promote the formation of highly-oxygenated particle-bound fragmentation products with lower molar mass and greater CCN activity, underscoring the importance of chemical aging in the free troposphere and its influence on the CCN activity of OA.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Dynamics of the Bingham Canyon rock avalanches (Utah, USA) resolved from topographic, seismic, and infrasound data: Bingham Canyon Rock Avalanches
The 2013 Bingham Canyon Mine rock avalanches represent one of the largest cumulative landslide events in recorded U.S. history and provide a unique opportunity to test remote analysis techniques for landslide characterization. We combine aerial photogrammetry surveying, topographic reconstruction, numerical runout modeling, and analysis of broadband seismic and infrasound data to extract salient details of the dynamics and evolution of the multiphase landslide event. Our results reveal a cumulative intact rock source volume of 52 Mm 3 , which mobilized in two main rock avalanche phases separated by 1.5 h. We estimate that the first rock avalanche had 1.5–2 times greater volume than the second. Each failure initiated by sliding along a gently dipping (21°), highly persistent basal fault before transitioning to a rock avalanche and spilling into the inner pit. The trajectory and duration of the two rock avalanches were reconstructed using runout modeling and independent force history inversion of intermediate-period (10–50 s) seismic data. Intermediate- and shorter-period (1–50 s) seismic data were sensitive to intervals of mass redirection and constrained finer details of the individual slide dynamics. Back projecting short-period (0.2–1 s) seismic energy, we located the two rock avalanches within 2 and 4 km of the mine. Further analysis of infrasound and seismic data revealed that the cumulative event included an additional 11 smaller landslides (volumes ~10 4 –10 5 m 3 ) and that a trailing signal following the second rock avalanche may result from an air-coupled Rayleigh wave. These results demonstrate new and refined techniques for detailed remote characterization of the dynamics and evolution of large landslides.
58 GEOSCIENCES↗
Recent advances in understanding secondary organic aerosol: Implications for global climate forcing: Advances in Secondary Organic Aerosol
Anthropogenic emissions and land use changes have modified atmospheric aerosol concentrations and size distributions over time. Understanding preindustrial conditions and changes in organic aerosol due to anthropogenic activities is important because these features (1) influence estimates of aerosol radiative forcing and (2) can confound estimates of the historical response of climate to increases in greenhouse gases. Secondary organic aerosol (SOA), formed in the atmosphere by oxidation of organic gases, represents a major fraction of global submicron-sized atmospheric organic aerosol. Over the past decade, significant advances in understanding SOA properties and formation mechanisms have occurred through measurements, yet current climate models typically do not comprehensively include all important processes. Our review summarizes some of the important developments during the past decade in understanding SOA formation. We also highlight the importance of some processes that influence the growth of SOA particles to sizes relevant for clouds and radiative forcing, including formation of extremely low volatility organics in the gas phase, acid-catalyzed multiphase chemistry of isoprene epoxydiols, particle-phase oligomerization, and physical properties such as volatility and viscosity. Several SOA processes highlighted in this review are complex and interdependent and have nonlinear effects on the properties, formation, and evolution of SOA. Current global models neglect this complexity and nonlinearity and thus are less likely to accurately predict the climate forcing of SOA and project future climate sensitivity to greenhouse gases. Efforts are also needed to rank the most influential processes and nonlinear process-related interactions, so that these processes can be accurately represented in atmospheric chemistry-climate models.
54 ENVIRONMENTAL SCIENCES↗
Capillary pressure-saturation relations in quartz and carbonate sands: Limitations for correlating capillary and wettability influences on air, oil, and supercritical CO 2 trapping
Capillary pressure (P c )-saturation (S w ) relations are essential for predicting equilibrium and flow of immiscible fluid pairs in soils and deeper geologic formations. In systems that are difficult to measure, behavior is often estimated based on capillary scaling of easily measured P c –S w relations (e.g., air-water, and oil-water), yet the reliability of such approximations needs to be examined. In this study, 17 sets of brine drainage and imbibition curves were measured with air-brine, decane-brine, and supercritical (sc) CO 2 -brine in homogeneous quartz and carbonate sands, using porous plate systems under ambient (0.1 MPa, 23°C) and reservoir (12.0 MPa, 45°C) conditions. Comparisons between these measurements showed significant differences in residual nonwetting phase saturation, S nw,r . Through applying capillary scaling, changes in interfacial properties were indicated, particularly wettability. With respect to the residual trapping of the nonwetting phases, S nwr, CO2 > S nwr, decane > S nwr, air . Decane-brine and scCO 2 -brine Pc–Sw curves deviated significantly from predictions assuming hydrophilic interactions. Moreover, neither the scaled capillary behavior nor S nw,r for scCO 2 -brine were well represented by decane-brine, apparently because of differences in wettability and viscosities, indicating limitations for using decane (and other organic liquids) as a surrogate fluid in studies intended to apply to geological carbon sequestration. Thus, challenges remain in applying scaling for predicting capillary trapping and multiphase displacement processes across such diverse fields as vadose zone hydrology, enhanced oil recovery, and geologic carbon sequestration.
58 GEOSCIENCES↗