Search NASA⌕ Search

SEARCH · Search NASA

Results for “CAPTURE EFFECT”

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 55 records · Page 3

Investigation of grain size and geometrically necessary dislocation density dependence of flow stress in Mg-4Al by using nanoindentation

Models for grain size dependence of flow stress typically ignore contributions from geometrically necessary dislocations (GNDs) in the grain boundary regions. Ashby had proposed the following equation for flow stress (τ) as a function of grain size (d) and effective plastic strain ($\bar{ε}$): $τ = τ_{0} + C'G\sqrt{^{b\bar{ε}}_{d}}$ using an unknown empirical coefficient, C’, In this study, the coefficient, C’, and equation for the grain size dependence of flow stress were investigated in Mg-4Al using nanoindentation. C’ ranges from 1.4 to 7.11 that is related to the hardness difference (∆H0) between the adjacent grains as $C' + 0.62\sqrt{ΔH_{0}}$. Ashby equation was used to predict uniaxial stress-strain curves in grain sizes of 55 µm, 187 µm, and 333 µm. The results showed lower flow stress and lower hardening rate compared to experimentally measured tensile stress-strain response. The discrepancy seems to be a result of hardening effect of SSDs. Here, we propose a modified Ashby equation that takes into account the hardening effect of SSDs and GNDs separately, resulting in a better agreement between calculated and experimentally measured values. This analysis demonstrates estimation of uniaxial stress-strain response using nanoindentation measurements that captures effects of grain size, and strain hardening from SSDs and GNDs. The combination of the revised equation and nanoindentation method will enable rapid exploration of stress-strain response of a diverse range of single-phase polycrystalline alloys.

36 MATERIALS SCIENCE↗

Modeling the Effects of Drift Orbit Bifurcation on the Magnetopause Shadowing Loss of Radiation Belt Electrons

Abstract Last closed drift shell (LCDS) has been identified as a crucial parameter for investigating the magnetopause shadowing loss of radiation belt electrons. However, drift orbit bifurcation (DOB) effects have not been physically incorporated into the LCDS calculation. Here we calculate event‐specific LCDS using different approaches to dealing with the DOB effects, that is, tracing field lines ignoring DOB, tracing test particles rejecting field lines with DOB, and tracing particles including field lines with DOB, and then incorporate them into a radial diffusion model to simulate the fast electron dropout observed by Van Allen Probes in May 2017. The model effectively captures the fast dropout at high L* and exhibits the best agreement with data when LCDS is calculated by tracing test particles with DOB more physically included. This study represents the first quantitative modeling of the DOB effects on radiation belt magnetopause shadowing loss via a more physical specification of LCDS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dual-Bed Radioiodine Capture from Complex Gas Streams with Zeolites: Regeneration and Reuse of Primary Sorbent Beds for Sustainable Waste Management

Dual-sorbent systems are proposed for radioiodine management with a regenerated primary bed for multiple cycles of use in complex conditions and a secondary bed for disposal with higher waste loadings. Sorbent approaches for the effective capture of gaseous radioiodine (isotopes 129 I and 131 I) produced from a range of nuclear processes have been studied for over half a century. (1−5) Whether or not a sorbent (e.g., molecular sieve) is required to physically screen/trap or chemically bind a radionuclide of interest through chemisorption, the complexity of the gas stream has a large impact on the performance (e.g., loading capacity, selectivity) and active life of a sorbent bed. (3) Silver mordenite (AgZ), the U.S. Department of Energy baseline sorbent for radioiodine capture from nuclear processes, performs well within acidic conditions and at elevated temperatures (6) and can be consolidated into a chemically durable waste form for long-term disposal. (7,8) However, new sorbents are being sought because optimal capture performance of AgZ significantly decreases in dynamic oxidizing environments with competing species, and it is expensive and it contains Ag (a toxic metal). (9) Until a new sorbent is found to replace AgZ, the regeneration and reuse of AgZ is an attractive alternative to a single-use primary sorbent bed. In this regard, a primary sorbent could be designed for enhanced capture in complex gas streams and the ability to be regenerated for reuse. Here, a secondary sorbent could then be tailored for maximum iodine loading in the gas stream and chemical durability within a disposal facility.

chemisorption↗

Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions

First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. In this work, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A simulation-based model studying monoethanolamine and aprotic heterocyclic anion ionic liquid (AHA-IL) mixtures for carbon capture

Capturing of waste CO 2 , particularly from large industrial point sources, is necessary for either permanent storage or CO 2 utilization. Amine solvents used for carbon capture have certain drawbacks, particularly the high energy required by the reboiler for solvent regeneration and degradation of the amine with absorption/regeneration cycling. These challenges may be overcome using ionic liquids (ILs) that can be tailored for operating under typical flue gas stream conditions (8–10% CO 2 , 18–20% H 2 O, 2–3% O 2 , and 67–72% N 2 ). ILs, however, oftentimes suffer from high viscosities, leading to diffusion limitations when compared to amine solvents. This work models the combinations of aprotic heterocyclic anion ILs (AHA-IL) with monoethanolamine (MEA) and studies the carbon capture effectiveness of these mixtures. Furthermore, this study found that all AHA-IL/MEA mixtures exhibited 32 % lower enthalpy of reaction (compared to MEA) and 64 % lower viscosity (compared to AHA-IL), which are desirable characteristics for CO 2 capture solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Testing the parametric model for self-interacting dark matter using matched halos in cosmological simulations

Here, we systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in Ref.Yang et al. (2023) with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and subhalos. We focus on the accuracy of the model in predicting halo density profiles at z = 0 and the evolution of maximum circular velocity. We find the model predictions agree with the simulations within 10%–50% for most of the simulated (sub)halos, 50%–100% for extreme cases. This indicates that the model effectively captures the gravothermal evolution of the halos with very strong, velocity-dependent self-interactions. For an example application, we apply the model to study the impact of various SIDM scenarios on strong lensing perturber systems, demonstrating its utility in predicting SIDM effects for small-scale structure analyses. Our findings confirm that the model is an effective tool for mapping CDM halos into their SIDM counterparts.

Yang, Daneng [Chinese Academy of Sciences, Nanjing↗

FEED Study of Carbon Capture Inc. DAC and CarbonCure Utilization Using United States Steel's Gary Works Waste Plant

The University of Illinois Urbana-Champaign (UIUC) is spearheading a front-end engineering design (FEED) study of a Direct Air Capture and Utilization System (DACUS) that can remove 5,000 tonnes/yr net CO2 from ambient air, utilizing it in the production of concrete. The designed system, if built, would be larger than any currently existing DAC collector. Such carbon negative technologies are critical to meeting the goals of the DOE’s program to accelerate climate-critical technology. DAC is a promising new technology for reducing the greenhouse gas in the atmosphere but is expensive, partly due to the energy required to capture and then release CO2 during the process. By integrating CarbonCapture Technologies Inc. DAC collector at United States Steel's Gary Works and utilizing steel plant's waste heat, energy, and location so energy & transportation costs can be minimized. CarbonCapture has developed an innovative DAC system using novel adsorbents to cost- effectively capture CO2. The captured, liquified CO2 will be trucked to ready-mix concrete plants across Illinois, Indiana, & Wisconsin, where CarbonCure will inject it into concrete during the mixing process. The carbon dioxide reacts with concrete, mineralizing into calcium carbonate (CaCO3), sequestering the CO2 into the concrete, and reducing the amount of cement, the manufacturing of which is among the major sources of CO2 emissions. This FEED study demonstrates a full CO2 value chain for DACUS from industrial facilities. It also provided an opportunity to assess the impact of this holistic approach on job creation, regional economic impact, and environmental justice.

Carbon management, ccus, carbon capture, utilizati↗

Coupled THM modeling of bentonite heating and hydration in tank tests with a new temperature-dependent water retention model

This study presents a coupled thermo-hydro-mechanical (THM) model for simulating the heating and hydration behavior of bentonite, a buffer material in deep geological repositories (DGRs). The model incorporates a new temperature-dependent soil water retention curve which captures the thermal-induced shift in water retention behavior. It also distinguishes between liquid and gas permeability, modeling intrinsic gas permeability as a function of accessible porosity to improve vapor transport and desaturation predictions. The model was validated against two large-scale tank tests, demonstrating good agreement with measured temperature, relative humidity, and water inflow data. It revealed a complex porosity evolution driven by thermal expansion, vapor movement, vapor condensation, and hydration-induced swelling during heating and hydration processes. The simulation results also suggest that the permeability of the hydration layer plays a critical role in controlling water intake. Clogging of this layer can significantly reduce the volume of water inflow during the hydration phase. Furthermore, while the model effectively captures key THM behavior, further development of the mechanical constitutive law is required to account for possible thermo-elasto-plastic volume changes and microstructural effects. Overall, the model provides a robust tool for evaluating the evolution of bentonite-based barrier material in DGRs.

Guo, Guanlong [Lawrence Berkeley National Laborato↗

Treyson Ricks - Intern Showcase Poster

Quinone-based sorbents offer a tunable, energy-efficient route to electrochemical CO2 capture, but systematic guidance for molecular design is lacking. Here, we report a high-throughput computational workflow that combines density functional theory (DFT) screening with machine-learning (ML) modeling to evaluate CO2 binding thermodynamics across several quinone derivatives, spanning benzoquinones, naphthoquinones, and anthraquinones. In addition to using solvents to stabilize the quinone anion and dianion, we studied the effect of ion-pairing on the reduction potentials and the CO2 binding energy. Automated Python scripts handled geometry optimizations and adduct-formation energies on an HPC cluster, reducing manual effort significantly. This integrated platform can uncover structure–property relationships and enables rapid in silico evaluation of untested candidates. We present one example from our workflow to showcase the capability of using quinones with ion-pairing to effectively capture CO2. Our approach paves the way for the rational selection of optimal quinone sorbents and can be extended with experimental thermochemical and kinetic data, alternative redox cycles, and stability assessments to accelerate development of next-generation electrochemical CO2 capture materials.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

High Selectivity Reactive Carbon Dioxide Capture over Zeolite Dual-Functional Materials

Reactive carbon dioxide capture (RCC) is a process where carbon dioxide (CO 2 ) is captured from a mixed gas stream (such as air) and converted to products without first performing a separation step to concentrate the CO 2 . Here, in this work, zeolite dual-functional materials (ZFMs) are introduced and evaluated for simulated RCC. The studied ZFMs feature high surface area, crystalline, microporous zeolite faujasite (FAU) as the support. Sodium oxide (“Na 2 O”) is impregnated as an effective capture agent capable of scavenging low concentration CO 2 (1,000 ppm). Exchanged and impregnated sodium on FAU chemisorbs CO 2 as carbonates and bicarbonates but does not promote the conversion of sorbed CO 2 to products when heated in hydrogen. The addition of Ru promotes the formation of formates, while the addition of Pt generates carbonyl surface species when heated in hydrogen. The active metal then promotes extremely high selectivity for CO 2 hydrogenation to either methane on Ru catalyst (~150 °C) or carbon monoxide on Pt catalyst (~200 °C) when heated in reducing atmospheres.

36 MATERIALS SCIENCE↗

Relaxation behavior in low-frequency complex conductivity of sands caused by bacterial growth and biofilm formation by Shewanella oneidensis under a high-salinity condition

Complex electrical conductivity is increasingly used to monitor subsurface processes associated with microbial activities because microbial cells mostly have surface charges and thus electrical double layers. Although highly saline environments are frequently encountered in coastal and marine sediments, there are limited data available on the complex conductivity associated with microbial activities under a high-salinity condition. Therefore, we have developed the spectral responses of complex conductivity of sand associated with bacterial growth and biofilm formation under a highly saline condition of approximately 1% salinity and approximately 2 S/m pore water conductivity with an emphasis on relaxation behavior. A column test is performed, in which the model bacteria Shewanella oneidensis MR-1 are stimulated for cell growth and biofilm formation in a sand pack, whereas the complex conductivity is monitored from 0.01 Hz to 10 kHz. The test results indicate that the real conductivity increases in the early stage due to the microbial metabolites and the increased surface conduction with cell growth but soon begin to decrease because of the reduction of charge passages due to bioclogging. However, the imaginary conductivity significantly increases with time, and clear bell-shaped relaxation behaviors are observed with the peak frequency of 0.1–1 Hz, associated with the double-layer polarization of cells and electrically conductive pili and biofilms. The Cole-Cole relaxation model appears to capture such relaxation behaviors well, and the modeling results indicate gradual increases in normalized chargeability and decreases in relaxation time during bacterial growth and biofilm formation in the highly saline condition. Comparison with previous literature confirms that the high-salinity condition further increases the normalized chargeability, whereas it suppresses the phase shift and thus the imaginary conductivity. Our results suggest that the complex conductivity can effectively capture microbial biomass formation in sands under a highly saline condition.

58 GEOSCIENCES↗

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗

Protection against Chemical Warfare Agents and Biological Threats Using Metal–Organic Frameworks as Active Layers

The SARS-CoV-2 pandemic outbreak and the unfortunate misuse of toxic chemical warfare agents (CWAs) highlight the importance of developing functional materials to protect against these chemical and pathogen threats. Metal–organic frameworks (MOFs), which comprise a tunable class of crystalline porous materials built from inorganic nodes and organic linkers, have emerged as a class of heterogeneous catalysts capable of rapid detoxification of multiple classes of these harmful chemical or biological hazards. In particular, zirconium-based MOFs (Zr-MOFs) feature Lewis acidic nodes that serve as active sites for a wide range of catalytic reactions, including the hydrolysis of organophosphorus nerve agents within seconds in basic aqueous solutions. In addition, postsynthetic modification of Zr-MOFs enables the release of active species capable of reacting with and deactivating harmful pathogens. Despite this impressive performance, utilizing Zr- MOFs in powder form is not practical for application in masks or protective uniforms. To address this challenge, our team sought to develop MOF/fiber composite systems that could be adapted for use under realistic operating conditions to protect civilians, military personnel, and first responders from harmful pathogens and chemical warfare agents. Over the last several years, our group has designed and fabricated reactive and biocidal MOF/fiber composites that effectively capture and deactivate these toxic species. In this Account, we describe the evolution of these porous and reactive MOF/fiber composites and focus on key design challenges and considerations. First, we devised a scalable method for the integration of Zr-MOFs onto textile substrates using aqueous precursor solutions and without using pretreated textiles, highlighting the potential scalability of this method. Moving beyond standard textiles, we also developed a microbial synthesis strategy to prepare hierarchically porous MOF/bacterial cellulose nanofiber composite sponges that can both capture and detoxify nerve agents when exposed to contaminated gas flows. The mass loading of the MOF in the nanofibrous composite sponge is up to 90%, affording higher work capacities compared to those of textile-fiber-based composites with relatively lower MOF loadings. Next, we demonstrated that heterogeneous polymeric bases are suitable replacements for volatile liquid bases typically used in solution-phase reactions, and we showed that these composite systems are capable of effectively hydrolyzing nerve agents in the solid state by using only water that is present as humidity. Moreover, incorporating a reactive dye precursor into the composite affords a dual function sensing and detoxifying material that changes color from white to orange upon reaction with the byproduct following nerve agent hydrolysis, demonstrating the versatility of this platform for use in decontamination applications. We then created chlorine-loaded MOF/fiber composites that act as biocidal and reactive textiles that are capable of not only detoxifying sulfur-mustard-based chemical warfare agents and simulants but also deactivating both bacteria and the SARS-CoV-2 virus within minutes of exposure. Lastly, we synthesized a mixed-metal Ti/Zr-MOF coating on cotton fibers to afford a photoactive biocidal cloth that shows fast and broad-spectrum biocidal performance against viruses and Gram-positive and Gram-negative bacteria under visible light irradiation. Given the tunable, multifunctional nature of these MOF/fiber composites, we believe that this Account will offer new insights for the rational design and preparation of functional MOF/fiber composites and pave the way toward the development of next-generation reactive and protective textiles.

36 MATERIALS SCIENCE↗

A hybrid CNN-LSTM surrogate model for hyper-resolution spatiotemporal flood forecasting in Norfolk, Virginia

Study region: Norfolk, Virginia, United States Study focus: Accurate and timely flood forecasting is essential for enhancing resilience in coastal urban areas in the context of increasing frequency and intensity of rainfall, sea level rise and urbanization. This study presents a hybrid deep learning-based surrogate model that integrates Convolutional Neural Networks (CNN) and Long Short-Term Memory (LSTM) networks to enable real-time spatiotemporal flood forecasting. The model leverages CNN to capture spatial features from inputs such as elevation and Topographic Wetness Index (TWI), while LSTM processes time-series inputs of rainfall and tide data to capture temporal features. New hydrologic insights for the region: The hybrid CNN-LSTM model was trained using the physics-based hydrodynamic model simulations obtained from the Two-dimensional Unsteady FLOW (TUFLOW) model for Norfolk, Virginia, and achieved high predictive accuracy across diverse flood-prone areas. The reduced computational time from four to six hours using TUFLOW to 3.2 min per event using CNN-LSTM enables rapid flood inundation mapping and early warning applications. The model effectively captured both spatial flood extents and their temporal evolution across different flooding scenarios, providing forecasts at a 2.5-m spatial resolution and 15-min temporal resolution and a one-hour-ahead prediction horizon. While challenges remain in terms of transferability to new regions and real-time data assimilation, this approach demonstrates strong potential for supporting operational flood risk management in coastal urban environments.

Coastal urban flooding↗

Efficient graph representation framework for chemical molecule similarity tasks

Graph data has emerged in numerous scientific domains and machine learning techniques have been widely used for analysis and learning of diverse data for prediction and decision. Machine learning techniques can readily address complex problems by leveraging their structural information. But graphs cannot be directly used for existing machine learning algorithms unless encoded as vectors. The problem of efficient representation of graphs is a substantial challenge in graph machine learning. In this paper, we propose a novel two-stage framework for the representation of chemical molecule graphs based on the strengths of Graph Isomorphism Networks (GINs) and Siamese autoencoders. In the first stage, the GIN model is constructed and trained using the structural information of chemical molecule graphs. Node attributes, edge attributes, and edge indices are used as input data, while graph attributes are used as labels. The GIN model effectively captures the structural characteristics of graphs and can accurately predict graph attributes, i.e., molecular properties. It also generates Graph Embeddings, represented as vectors that encode the structural information of graphs. In the second stage, Graph Embedding vectors are further optimized for downstream similarity tasks while preserving the graph structural information. The Siamese autoencoder is constructed and trained, which reduces the dimensionality of the Graph Embedding vectors, while maximizing the preservation of structural information in the original high-dimensional vectors. The resulting low-dimensional Graph Embeddings can be effectively utilized for tasks such as approximate nearest neighbor search. The experimental results demonstrate the effectiveness of our proposed framework in accurately predicting graph similarity.

Ma, Jiaji↗

Inferring Plant Acclimation and Improving Model Generalizability With Differentiable Physics‐Informed Machine Learning of Photosynthesis

Net photosynthesis (A N ) is a key component of the global carbon cycle influencing climate feedback over decadal scales. Although plant acclimation to environmental changes can modify A N , traditional vegetation models in Earth system models (ESMs) often rely on plant functional type (PFT)-specific parameterizations or simplified acclimation assumptions limiting generalizability across time, space, and PFTs. In this study, we developed a differentiable photosynthesis model to learn the environmental dependencies of V c,max25 (maximum carboxylation rate at 25°C, representing photosynthetic capacity), as this genre of hybrid physics-informed machine learning can seamlessly train neural networks and process-based equations together. Compared to PFT-specific parameterization of V c,max25 , learning the environment dependencies of key photosynthetic parameters improved model spatiotemporal generalizability. Applying environmental acclimation to V c,max25 led to substantial variations in global mean A N indicating the need to address acclimation in ESMs. The model effectively captured multivariate observations (V c,max25 , A N , and stomatal conductance (g s )) simultaneously with multivariate constraints, improving generalization across space and PFTs. It also learned sensible acclimation relationships of V c,max25 to different environmental conditions. The model explained more than 54%, 57%, and 62% of the variance of A N , g s , and V c,max25 , respectively, presenting a first global-scale spatial test benchmark of A N and g s . These results highlight the potential for differentiable modeling to enhance process-based modules in ESMs and effectively leverage information from large, multivariate data sets.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of a Process to Separate Ensiled Corn Grain and Stover for Improved Utilization as a Biomass Feedstock

Corn grain and chopped stover were stored anaerobically at 37% to 50% (w.b.) aggregate moisture content to create a unique biomass feedstock. However, after storage, the two fractions must be separated to accommodate different conversion pathways. In this work, a modified cleaning system from a grain combine harvester was evaluated to investigate its effectiveness in separating the grain from the stover. While this system has been purposefully built for this task, there is no previous literature on separating grain at high material other than grain (MOG) to grain ratios. Using this system, material was separated into four fractions: Grain, Tailings, Heavy MOG, and Light MOG. Subsamples were collected, oven dried, and then hydrodynamically separated to quantify the grain content in the four fractions. Several different configurations of air velocities, mass flow rate, material MOG: grain ratios, and sieve types and openings were investigated. Grain capture was defined as the fraction of the total grain dry mass collected in the grain and tailings fractions. In five experiments, the grain capture effectiveness varied from 89.2% to 99.4% on a dry basis. The dry basis foreign matter (FM) in the clean grain was 0.6% to 10.8%. As a result, the use of a modified combine cleaning system has the potential to be one step in a system to fractionate corn grain from stover (i.e., MOG) in a biorefinery that uses both starch and cellulose as biomass feedstocks.

09 BIOMASS FUELS↗

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-art “black-box” solvers such as algebraic multigrid (AMG) are ineffective because they cannot effectively capture the strong coupling between the mechanics and the flow subproblems, as well as the coupling inherent in the multiphase flow and transport process. In this work, we develop an algebraic framework based on multigrid reduction (MGR) that is suited for tightly coupled systems of PDEs. Using this framework, the decoupling between the equations is done algebraically through defining appropriate interpolation and restriction operators. One can then employ existing solvers for each of the decoupled blocks or design a new solver based on knowledge of the physics. We demonstrate the applicability of our framework when used as a “black-box” solver for multiphase poromechanics. Here, we show that the framework is flexible to accommodate a wide range of scenarios, as well as efficient and scalable for large problems.

97 MATHEMATICS AND COMPUTING↗