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 955 records · Page 53

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 24C 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↗

An extended JBN method of determining unsteady-state two-phase relative permeability

Relative permeability is the reduction of permeability of porous media when subjected to multiphase flow and a key parameter in subsurface hydrology. The JBN method is a well-known method of obtaining relative permeability, which measures the overall pressure drop and the effluent phase ratio versus time during two-phase displacements. By assuming no capillary pressure or gravity, the JBN method obtains the relative permeabilities to both phases at the core outlet. Since data across a range of saturations are acquired in a relatively short time, this method is widely used. This work extends the JBN method by having (1) section-wise pressure drop measurements between the core inlet, four pressure taps on the core and the outlet, (2) local saturation measurements, and (3) local phase fluxes. With these data, the extended JBN method can determine relative permeabilities to both phases at each pressure tap of the core (not just at the core outlet). The JBN extension is shown using a data set where CO 2 invades a brine-filled core. From this it is found that the advantages of the extended JBN method over the regular JBN method are: (1) four times more data are obtained, and (2) data are more accurate because the capillary end effect is experimentally avoided. Avoiding the end effect results in tripling the saturation range, and obtaining relative permeabilities that are consistent with steady-state measurements and roughly 40% higher than those from the regular JBN method.

54 ENVIRONMENTAL SCIENCES↗

Pore-scale water dynamics during drying and the impacts of structure and surface wettability

Plants and microbes secrete mucilage into soil during dry conditions, which can alter soil structure and increase contact angle. Structured soils exhibit a broad pore size distribution with many small and many large pores, and strong capillary forces in narrow pores can retain moisture in soil aggregates. Meanwhile, contact angle determines the water repellency of soils, which can result in suppressed evaporation rates. Although they are often studied independently, both structure and contact angle influence water movement, distribution, and retention in soils. Here drying experiments were conducted using soil micromodels patterned to emulate different aggregation states of a sandy loam soil. Micromodels were treated to exhibit contact angles representative of those in bulk soil (8.4° ± 1.9°) and the rhizosphere (65° ± 9.2°). Drying was simulated using a lattice Boltzmann single-component, multiphase model. In our experiments, micromodels with higher contact angle surfaces took 4 times longer to completely dry versus micromodels with lower contact angle surfaces. Microstructure influenced drying rate as a function of saturation and controlled the spatial distribution of moisture within micromodels. In conclusion, lattice Boltzmann simulations accurately predicted pore-scale moisture retention patterns within micromodels with different structures and contact angles.

Flow, Transport, and Reaction in Porous Media: Percolation Scaling, Critical‐Path Analysis, and Effective Medium Approximation

Abstract We describe the most important developments in the application of three theoretical tools to modeling of the morphology of porous media and flow and transport processes in them. One tool is percolation theory. Although it was over 40 years ago that the possibility of using percolation theory to describe flow and transport processes in porous media was first raised, new models and concepts, as well as new variants of the original percolation model are still being developed for various applications to flow phenomena in porous media. The other two approaches, closely related to percolation theory, are the critical‐path analysis, which is applicable when porous media are highly heterogeneous, and the effective medium approximation— poor man's percolation —that provide a simple and, under certain conditions, quantitatively correct description of transport in porous media in which percolation‐type disorder is relevant. Applications to topics in geosciences include predictions of the hydraulic conductivity and air permeability, solute and gas diffusion that are particularly important in ecohydrological applications and land‐surface interactions, and multiphase flow in porous media, as well as non‐Gaussian solute transport, and flow morphologies associated with imbibition into unsaturated fractures. We describe new applications of percolation theory of solute transport to chemical weathering and soil formation, geomorphology, and elemental cycling through the terrestrial Earth surface. Wherever quantitatively accurate predictions of such quantities are relevant, so are the techniques presented here. Whenever possible, the theoretical predictions are compared with the relevant experimental data. In practically all the cases, the agreement between the theoretical predictions and the data is excellent. Also discussed are possible future directions in the application of such concepts to many other phenomena in geosciences.

Hunt, Allen G.↗

The Effect of Wettability Heterogeneity on Relative Permeability of Two-Phase Flow in Porous Media: A Lattice Boltzmann Study

Relative permeability is a critical parameter characterizing multiphase flow in porous media and it is strongly dependent on the wettability. In many situations, the porous media are nonuniformly wet. In this study, to investigate the effect of wettability heterogeneity on relative permeability of two-phase flow in porous media, a multi-relaxation-time color-gradient lattice Boltzmann model is adopted to simulate oil/water two-phase flow in porous media with different oil-wet solid fractions. For the water phase, when the water saturation is high, the relative permeability of water increases with the increase of oil-wet solid fraction under a constant water saturation. However, as the water saturation decreases to an intermediate value (about 0.4–0.7), the relative permeability of water in fractionally wet porous media could be lower than that in purely water-wet porous media, meaning additional flow resistance exists in the fractionally wet porous media. For the oil phase, similar phenomenon is observed. This phenomenon is mainly caused by the wettability-related microscale fluid distribution. According to both our simulation results and theoretical analysis, it is found that the relative permeability of two-phase flow in porous media is strongly related to three parameters: the fluid saturation, the specific interfacial length of fluid, and the fluid tortuosity in the flow direction. Lastly, the relationship between the relative permeability and these parameters under different capillary numbers is explored in this paper.

58 GEOSCIENCES↗

Assessment of CO 2 Storage Potential in Naturally Fractured Reservoirs With Dual–Porosity Models

Naturally Fractured Reservoirs (NFR’s) have received little attention as potential CO 2 storage sites. Two main facts deter from storage projects in fractured reservoirs: (1) CO 2 tends to be nonwetting in target formations and capillary forces will keep CO 2 in the fractures, which typically have low pore volume; and (2) the high conductivity of the fractures may lead to increased spatial spreading of the CO 2 plume. Numerical simulations are a powerful tool to understand the physics behind brine-CO 2 flow in NFR’s. Dual-porosity models are typically used to simulate multiphase flow in fractured formations. However, existing dual-porosity models are based on crude approximations of the matrix-fracture fluid transfer processes and often fail to capture the dynamics of fluid exchange accurately. Therefore, more accurate transfer functions are needed in order to evaluate the CO 2 transfer to the matrix. This work presents an assessment of CO 2 storage potential in NFR’s using dual-porosity models. We investigate the impact of a system of fractures on storage in a saline aquifer, by analyzing the time scales of brine drainage by CO 2 in the matrix blocks and the maximum CO 2 that can be stored in the rock matrix. A new model to estimate drainage time scales is developed and used in a transfer function for dual-porosity simulations. Furthermore, we then analyze how injection rates should be limited in order to avoid early spill of CO 2 (lost control of the plume) on a conceptual anticline model. Numerical simulations on the anticline show that naturally fractured reservoirs may be used to store CO 2 .

54 ENVIRONMENTAL SCIENCES↗

Fundamental Effects of Al and Ta on Microstructure and Phase Transformations in the Al–Cr–Mo–Ta–Ti Refractory Complex Concentrated Alloy System

The effect of aluminum and tantalum concentrations on a refractory metal complex concentrated alloy is reported, particularly with respect to their effect on microstructure and phase composition of the alloy in cast and annealed form. Alloys with an equiatomic composition, (AlCrMoTaTi), an aluminum-lean composition (Al 0.75 CrMoTaTi), and a tantalum-lean composition (AlCrMoTa 0.75 Ti) are produced via arc melting. The alloys exhibit multiphase structures, confirmed by X-ray diffraction, microstructural characterization, and thermal analysis. The minor off-equiatomic adjustments of aluminum and tantalum in this alloy system did not drastically alter the prevalence of the Cr–Ta-based Laves phase. Correlations between thermodynamic predictions and observed phase transformations via thermal analysis are improved upon refinement of calculations removing impractical intermediate phases. As a result, experimental findings provide information for the refinement of thermodynamic modeling and deliver additional insight into the optimization of alloy compositions within this five-component system.

36 MATERIALS SCIENCE↗

A new class of renewable thermoplastics with extraordinary performance from nanostructured lignin-elastomers

In this study, a new class of thermoplastic elastomers has been created by introducing nanoscale-dispersed lignin (a biomass-derived phenolic oligomer) into nitrile rubber. Temperature-induced controlled miscibility between the lignin and the rubber during high shear melt-phase synthesis allows tuning the material s morphology and performance. The sustainable product has unprecedented yield stress (15–45 MPa), strain hardens at large deformation, and has outstanding recyclability. The multiphase polymers developed from an equal-mass mixture of a melt-stable lignin fraction and nitrile rubber with optimal acrylonitrile content, using the method described here, show 5–100 nm lignin lamellae with a high-modulus rubbery interphase. Molded or printed elastomeric products prepared from the lignin-nitrile material offer an additional revenue stream to pulping mills and biorefineries. Research was sponsored by the Technology Innovation Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the U.S. Department of Energy. TEM (J.C. and C.D.T.) and SAXS (J.K.K.) experiments were conducted at the Center for Nanophase Materials Sciences, which is a DOE Office of Science User Facility.

36 MATERIALS SCIENCE↗