Search NASA⌕ Search

SEARCH · Search NASA

Results for “Carbon Sequestration”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evaluation of accuracy and convergence of numerical coupling approaches for poroelasticity benchmark problems

Accurate modeling of subsurface flow and transport processes is vital as the prevalence of subsurface activities such as carbon sequestration, geothermal recovery, and nuclear waste disposal increases. Computational modeling of these problems leverages poroelasticity theory, which describes coupled fluid flow and mechanical deformation. Although fully coupled monolithic schemes are accurate for coupled problems, they can demand significant computational resources for large problems. In this work, a fixed stress scheme is implemented into the Sandia Sierra Multiphysics toolkit. Two implementation methods, along with the fully coupled method, are verified with one-dimensional (1D) Terzaghi, 2D Mandel, and 3D Cryer sphere benchmark problems. The impact of a range of material parameters and convergence tolerances on numerical accuracy and efficiency was evaluated. Overall the fixed stress schemes achieved acceptable numerical accuracy and efficiency compared to the fully coupled scheme. However, the accuracy of the fixed stress scheme tends to decrease with low permeable cases, requiring the finer tolerance to achieve a desired numerical accuracy. For the fully coupled scheme, high numerical accuracy was observed in most of cases except a low permeability case where an order of magnitude finer tolerance was required for accurate results. Finally, a two-layer Terzaghi problem and an injection–production well system were used to demonstrate the applicability of findings from the benchmark problems for more realistic conditions over a range of permeability. Simulation results suggest that the fixed stress scheme provides accurate solutions for all cases considered with the proper adjustment of the tolerance. Furthermore, this work clearly demonstrates the robustness of the fixed stress scheme for coupled poroelastic problems, while a cautious selection of numerical tolerance may be required under certain conditions with low permeable materials.

58 GEOSCIENCES↗
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
DOE OSTI2020

Mechanistic modeling of CO 2 well leakage in a generic abandoned well through a bridge plug cement-casing gap

Both known and unmapped plugged and abandoned wells are potential leakage pathways for CO 2 from geologic carbon sequestration (GCS) sites. Although many abandoned wells have cement bridge plugs installed to prevent leakage, the seal between the cement and the inner casing wall is subject to failure. In this study, we carried out detailed T2Well simulations of cases of sudden non-Darcy flow of CO 2 and brine leakage up the gap between a cement plug and the inner steel casing wall that becomes a fully connected flow path during the post-injection period. The goal of our study was two-fold: (1) to understand the dynamics, rates, and the characteristic temporal signals associated with the onset of leakage through various gap-aperture sizes, and (2) to suggest potential monitoring strategies based on the findings. Simulation results show that the leakage of CO 2 and brine upward is transient with interesting phase interference behavior. Time-dependent oscillatory flows with varying pressure, temperature, and flow rates of CO 2 and brine show strong dependence on gap aperture. Phase-change and decompression lead to very low temperatures at the top of the well for gap apertures larger than 4 mm suggesting that remote thermal monitoring at the ground surface may be an effective way of monitoring even if well locations are not known a priori. Pressure in the well is also indicative of CO 2 leakage. The temporal patterns of changing temperature and pressure may be useful diagnostic signals for leakage detection. Finally, these transient leakage signals may provide information on the cause of leakage and/or characteristics of the flow path that could inform effective remediation design and execution approaches.

58 GEOSCIENCES↗
DOE OSTI2022

Development of lean, efficient, and fast physics-framed deep-learning-based proxy models for subsurface carbon storage

In this work, we present deep-learning-based surrogate models for CCUS developed with four different algorithms and a physics-framed two-phase flow problem involving displacement of water by CO 2 . The deep-learning models were trained using 3D datasets describing the pressure plume, CO 2 saturation plume, and water extraction rate generated by numerical simulation. The hyperparameters defining the architecture of the neural networks were optimized to determine the slimmest network size and training parameters that give the most efficient performance at the least training cost. To develop a robust model that closely mimics the governing physical laws, the discretized form of the two-phase fluid transport equation was used to formulate the supervised deep-learning task. The algorithms investigated in this study predicted the data to above 95% accuracy, with the multi-layer perceptron model demonstrating the best performance by balancing training speed, prediction time, and prediction accuracy with lean network capacity. Furthermore, the surrogate models simultaneously predict reservoir pressure and CO 2 saturation in every grid block, including the surface well extraction rate and bottomhole pressure, at all simulation times for a given static model realization in just a few seconds on a standard desktop computer. A key outcome of this study is that limits can be placed on network design parameters to avoid over designing neural networks, with associated efficiencies in training and prediction times. This is very useful because large volumes of data may be generated in CCUS projects and over-design of neural network architectures imposes penalties that are antithetical to the goal of near-real time forecasting.

58 GEOSCIENCES↗
DOE OSTI2024

An assessment of controlled source EM for monitoring subsurface CO 2 injection at the wyoming carbonSAFE geologic carbon storage site

Here we evaluate if electromagnetic (EM) geophysical methods for monitoring geologic carbon storage (GCS) efforts at the Wyoming CarbonSAFE project adjacent to the Dry Fork Station power plant near Gillette, Wyoming. This first involved acquiring both electric and magnetic fields at eleven different locations ranging in distance from immediately adjacent to 4 km from the plant. Passive EM measurements were made to provide spectral EM noise measurements generated by electricity production at the plant and to determine if useful magnetotelluric (MT) data can be successfully collected in the region. The processed data indicate that useful MT data can be collected as long as the site is located more than 2km away from the power plant as well as active roads and rail lines. Controlled source EM data were collected using three different source configurations, two of which connected to steel casings used to complete the injection wells. Comparing the EM noise measurements to the CSEM data show measurable electric and magnetic field signals at all sites. Next a series of three-dimensional (3D) numerical models were built that simulate resistivity changes caused by the proposed CO2 injection at depths ranging from 2.4 to 3.0km. These models were used to simulate various EM measurement configurations. The modeling shows that casing-source CSEM monitoring can provide sensitivity to the injected CO 2 if source electrodes are connected to the bottom of one or both of the injection wells.

58 GEOSCIENCES↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗
DOE OSTI2026

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2021

Deep learning inversion of gravity data for detection of CO 2 plumes in overlying aquifers

In this work, we developed an effective U-Net based deep learning (DL) model for inversion of surface gravity data on a rectangular grid to predict 2-D high-resolution subsurface CO 2 distribution along a vertical cross-section due to CO 2 leakage through a wellbore within a deep CO 2 storage reservoir. We used synthetic data to model two types of CO 2 leakage scenarios: one CO 2 plume in a shallow aquifer (single plume case), and two plumes present at different depths (double plume case). The 3-D synthetic plume samples were created by sampling among predetermined CO 2 plume depths, saturations, and volumes. The corresponding surface gravity data on a rectangular grid were generated by a 3-D forward model. The U-Net model detected 72% of single-plume samples, and one or both plumes in 75% of double-plume samples. Most of the undetected single plumes have small gravity field strengths below the typical noise level of 5 μGal. This model generated reproducible, reliable predictions with acceptable errors and demonstrated improved spatial resolution over the conventional least-squares inversion. In contrast to the conventional least-squares inversion, which often overestimates the size of its target and underestimates its density, this U-Net model accurately delineated the boundary of a target. Furthermore, this DL inversion detected deep, small, or low saturation CO 2 plumes that are often more difficult to resolve with conventional gravity inversion methods. We note the limitations of this feasibility study, including the use of synthetic data with regular CO 2 plume shapes, and the prediction of a 2-D plume cross-section rather than the full 3-D plume, as well, we recognize the lower detection fraction for double-plume scenarios. Nevertheless, this study demonstrates that DL gravity inversion is a promising and potentially superior method to conventional least-squares inversion. Our U-Net based deep learning inversion approach may be adapted for inversion of other types of geophysical data. DL inversion can facilitate near real-time monitoring of geologic carbon sequestration to provide site operators with prompt information about subsurface CO 2 distribution for risk management and mitigation.

58 GEOSCIENCES↗
DOE OSTI2025

Comparative life cycle assessment of a modular cross-laminated timber residential building designed for disassembly and reuse versus traditional wood frame construction

There is a need for affordable housing across the U.S., with high-performance modular and prefabricated buildings providing a logical avenue for meeting some of this demand. However, there is a need to balance high performance construction – including low emissions – with affordability. To provide a proof-of-concept in meeting these goals, the Circular Home is a cross-laminated timber (CLT)-based deconstructible and reconfigurable single-family residence that meets high performance targets in moisture, energy, design, economics, and life cycle assessment (LCA). This study focuses on the LCA, presenting a cradle-to-cradle whole-building life cycle assessment (WBLCA) for the Circular Home and a functionally equivalent Baseline Home constructed with traditional materials and methods. The functional unit is 1 m 2 of gross floor area across 60 years. Revit building information models (BIM) provided material quantities and Tally LCA was utilized for impact data (inclusive of biogenic carbon sequestration), supplemented with manufacturer environmental product declarations (EPDs). The Circular Home outperforms the baseline residence in most measured impact categories, including global warming potential (GWP), producing −2.73 kgCO 2 eq/m 2 in embodied emissions, whereas the modeled baseline has an embodied GWP of 428 kgCO 2 eq/m 2 . The careful material selection and advanced building design optimizes performance, with the Circular Home containing only −0.006 times the embodied emissions and −0.02 times the operational emissions of its traditional counterpart. Finally, the unique contribution of this work is in the environmental impact comparison of a high-performance modular CLT structure that can be affordably scaled and mass produced in a U.S. market, compared to typical single family home construction.

Circularity↗
DOE OSTI2021

Restoring wetlands on intensive agricultural lands modifies nitrogen cycling microbial communities and reduces N 2 O production potential

The concentration of nitrous oxide (N 2 O), an ozone-depleting greenhouse gas, is rapidly increasing in the atmosphere. Most atmospheric N 2 O originates in terrestrial ecosystems, of which the majority can be attributed to microbial cycling of nitrogen in agricultural soils. Here, in this study, we demonstrate how the abundance of nitrogen cycling genes vary across intensively managed agricultural fields and adjacent restored wetlands in the Sacramento-San Joaquin Delta in California, USA. We found that the abundances of nirS and nirK genes were highest at the intensively managed organic-rich cornfield and significantly outnumber any other gene abundances, suggesting very high N 2 O production potential. The quantity of nitrogen transforming genes, particularly those responsible for denitrification, nitrification and DNRA, were highest in the agricultural sites, whereas nitrogen fixation and ANAMMOX was strongly associated with the wetland sites. Although the abundance of nosZ genes was also high at the agricultural sites, the ratio of nosZ genes to nir genes was significantly higher in wetland sites indicating that these sites could act as a sink of N 2 O. These findings suggest that wetland restoration could be a promising natural climate solution not only for carbon sequestration but also for reduced N 2 O emissions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Damage detection in a laboratory-scale wellbore applying Time Reversal and Nonlinear Elastic Wave Spectroscopy (TR NEWS)

We report Time Reversal and Nonlinear Elastic Wave Spectroscopy (TR-NEWS) has been used to focus acoustic energy, and make measurements correlated with damage in a variety of industrial materials. Most studies using TR-NEWS have focused on relatively small objects and may have multiple acoustic sources. In Earth, within energy extraction and carbon sequestration settings, the structure and scale of wellbores makes acoustic focusing challenging, as has been demonstrated through numerical simulations. This paper explores the feasibility of applying TR NEWS for damage detection in wellbores by constructing a laboratory-scale wellbore, and using TR to focus and make dynamic linear and nonlinear elastic measurements on it. During successive cycles of localized mechanical damage induced in the sample, the hysteretic nonlinear elastic parameter α increased with progressive mechanical damage. In addition to these strain-dependent changes in α, the width and peak amplitude of the time reversal focused energy changed near the damage site.

58 GEOSCIENCES↗
DOE OSTI2020

Distinct Assembly Processes and Microbial Communities Constrain Soil Organic Carbon Formation

Soil stores more carbon (C) than all vegetation and the atmosphere combined. Soil C stocks are broadly shaped by temperature, moisture, soil physical characteristics, vegetation, and microbial-mediated metabolic processes. The efficiency with which microorganisms use soil C regulates the balance between C storage in soil and the atmosphere. In this review, we discuss how microbial physiology and community assembly processes determine microbial growth rate and efficiency and, in turn, soil organic C cycling through the lens of community ecology. We introduce a conceptual framework cataloging life history (i.e., growth rate, resource acquisition, and stress tolerance) and assembly traits (i.e., competition, facilitation, and dispersal) that correspond with different growth efficiencies. We also compare how dominant mycorrhizal fungal type affects growth efficiency. We propose that further development and inclusion of specific community parameters in microbial-explicit Earth system models are needed for accurately predicting soil organic C responses to global change.

carbon sequestration↗
DOE OSTI2021

Estimating induced land use change emissions for sustainable aviation biofuel pathways

Sustainable aviation fuels (SAFs) are expected to play an essential role in achieving the aviation industries’ goal of carbon-neutral growth. However, producing biomass-based SAFs may induce changes in global land use and the associated carbon stock. The induced land use change (ILUC) emissions, as a part of the full life-cycle emissions for SAF pathways, will affect whether and to what extent SAFs reduce emissions compared with petroleum-based jet fuels. Here, we estimate the ILUC emission intensity for seventeen SAF pathways considered by the International Civil Aviation Organization (ICAO), covering five ASTM-certified technologies, nine biomass-based feedstocks, and four geographical regions. We introduce the SAF pathways into a well-established computable general equilibrium (CGE) model, GTAP-BIO, and its coupled emission accounting model, AEZ-EF, to study economy-wide implications of SAF production and estimate ILUC emissions intensity for each pathway. The estimated SAF ILUC emission intensities, using a 25-year amortization period, range from -58.5 g CO 2 e MJ -1 for the USA miscanthus alcohol (isobutanol)-to-jet (ATJ) pathway to 34.6 g CO 2 e MJ -1 for the Malaysia & Indonesia palm oil Hydrotreated Esters of Fatty Acids (HEFA) pathway. Notably, the vegetable oil pathways tend to have higher ILUC emission intensities due to their linkage to palm expansion and peatland oxidation in Southeast Asia. The cellulosic pathways studied provide negative ILUC emissions, mainly driven by the high carbon sequestrations in crop biomass and soil. Using the core life-cycle emissions established by ICAO, we show that fifteen of the assessed pathways have a lower full life-cycle emission intensity than petroleum-based jet fuels (89 g CO 2 e MJ -1 ), offering promising options to reduce aviation emissions.

09 BIOMASS FUELS↗
DOE OSTI2022

Achieving superior carbon transfer efficiency and pH control using membrane carbonation with a wide range of CO 2 contents for the coccolithophore Emiliania huxleyi

The economic viability of microalgal-derived products relies on rapid CO 2 transfer in a cost-effective manner. Many industrial gas streams contain concentrated CO 2 that, if converted to useful products, would lower greenhouse gas emissions and valorize the wasted CO 2 . Membrane carbonation (MC) uses non-porous hollow-fiber gas-transfer membranes to deliver CO 2 without bubble formation, which makes it possible to achieve a high carbon-transfer efficiency (CTE). However, inert gasses in the industrial streams (e.g., N 2 , O 2 , and H 2 O) can significantly lower the CO 2 -delivery rate. The means to overcome the buildup of inert gases in the membrane lumen is to manage the distal end of the membranes to sweep out inert gases while not wasting significant CO 2 . A MC-venting strategy was evaluated for CO 2 inputs from 5% to 100%. Abiotic tests using a restricted exit flow could achieve >95% CTE abiotic for industrial CO 2 streams. When integrated with semi-continuous cultivation of a marine coccolithophore, Emiliania huxleyi , CO 2 delivery and venting were on-demand based on a pH set points and pH-actuated feed and venting valves. MC using the venting strategy achieved 100% CTE biotic when delivering 100% and 50% CO 2 , which was better than 50% CTE biotic obtained from pH-controlled sparging of 100% CO 2 -sparging. E. huxleyi consistently fixed ~80% of the delivered CO 2 into biomass, and the remaining ~20% to calcite coccoliths. Finally, the compact size of MC modules, stable pH control, and no shear forces from bubble agitation during the CO 2 delivery made MC an ideal match for cultivation of coccolithophores, which are sensitive to shear forces and pH fluctuations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

From pixels to patterns: Coupling Optical Coherence Tomography and machine learning for monitoring coastal wetland root systems

Coastal wetlands are crucial in shoreline stabilization, carbon sequestration, and storm protection. Yet, due to limitations in traditional destructive sampling techniques, the belowground biomass (live root mass) and necromass (dead and decaying roots) remain difficult to assess in coastal wetlands, limiting our understanding on coastal resilience, nutrient cycling, and soil structure. This study employs Optical Coherence Tomography (OCT) as a high-resolution imaging technique to analyze root biomass and necromass in the Terrebonne Basin, Louisiana. A Random Forest (RF) model was developed to classify root health states based on OCT-derived features, achieving an accuracy of 70% in distinguishing live from dead root segments. The results demonstrate that OCT, combined with ML, offers a promising novel approach to root analysis, providing fine-scale insights into root morphology and decay patterns that are not easily captured by conventional methods. This research lays the foundation for future integration of OCT with complementary imaging modalities such as X-ray Computed Tomography (XCT) and advanced ML algorithms to enhance classification accuracy and scalability. Future work aims to expand the dataset diversity across different wetland types and apply the methodology for large-scale, repeatable assessments of root biomass turnover and accumulation, with important implications for wetland monitoring, conservation, and restoration under changing environmental conditions.

AI/ML↗
DOE OSTI2025

Life Cycle Emissions and Health Cost Impacts of Producing Ethanol and Electricity from Willow and Switchgrass in the Riparian Buffers of Mid-Atlantic United States

The United States (US Mid-Atlantic Region (MAR) has the potential to grow a variety of perennial feedstocks such as switchgrass and shrub willow to increase domestic energy production. These cellulosic feedstocks have also shown improved ecosystem services, such as soil carbon sequestration, nitrate leaching reduction, and flood mitigation along rivers and streams as partially harvested riparian buffers. To examine the effects on greenhouse gases (GHGs) and criteria air pollutants (CAPs) from using these feedstocks to produce ethanol or electricity, we conducted a comprehensive life cycle assessment (LCA) and estimated the impact on human health costs when land use is changed from corn production for ethanol. Results indicate up to 54% reduced GHG per hectare from using willow and switchgrass feedstock sources to produce ethanol instead of corn. However, there was a trade-off in terms of CAP emission, as grass-based energy emitted more NO x and SO x compared to the corn ethanol pathway, except for SO X emissions from willow-based electricity. Electricity from cellulosic biomass had higher particulate matter (PM) emission compared to that from corn ethanol. Estimates for health cost to society ranged from $2498 ha –1 for electricity from switchgrass to a net benefit of $448 ha –1 for ethanol production from willow, depending on varying biomass yield under different market scenarios. Although using cellulosic feedstocks to produce bioenergy has great potential to reduce GHG emissions, CAP control measures are needed to manage CAP-induced health costs.

Bioenergy↗
DOE OSTI2020

Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

02 PETROLEUM↗