Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SANDSTONES”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Manganese exposure from spring and well waters in the Shenandoah Valley: interplay of aquifer lithology, soil composition, and redox conditions

Manganese (Mn) is of particular concern in groundwater, as low-level chronic exposure to aqueous Mn concentrations in drinking water can result in a variety of health and neurodevelopmental effects. Much of the global population relies on drinking water sourced from karst aquifers. Thus, we seek to assess the relative risk of Mn contamination in karst by investigating the Shenandoah Valley, VA region, as it is underlain by both karst and non-karst aquifers and much of the population relies on water wells and spring water. Water and soil samples were collected throughout the Shenandoah Valley, to supplement pre-existing well water and spring data from the National Water Information System and the Virginia Household Water Quality Program, totaling 1815 wells and 119 springs. Soils were analyzed using X-ray fluorescence and Mn K-Edge X-ray absorption near-edge structure spectroscopy. Factors such as soil type, soil geochemistry, and aquifer lithology were linked with each location to determine if correlations exist with aqueous Mn concentrations. Analyzing the distribution of Mn in drinking water sources suggests that water wells and springs within karst aquifers are preferable with respect to chronic Mn exposure, with < 4.9% of wells and springs in dolostone and limestone aquifers exceeding 100 ppb Mn, while sandstone and shale aquifers have a heightened risk, with > 20% of wells exceeding 100 ppb Mn. The geochemistry of associated soils and spatial relationships to various hydrologic and geologic features indicates that water interactions with aquifer lithology and soils contribute to aqueous Mn concentrations. Relationships between aqueous Mn in spring waters and Mn in soils indicate that increasing aqueous Mn is correlated with decreasing soil Mn(IV). These results point to redox conditions exerting a dominant control on Mn in this region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

The Geomechanics of CO 2 Storage in Deep Sedimentary Formations

This study provides a review of the geomechanics and modeling of geomechanics associated with geologic carbon storage (GCS), focusing on storage in deep sedimentary formations, in particular saline aquifers. The paper first introduces the concept of storage in deep sedimentary formations, the geomechanical processes and issues related with such an operation, and the relevant geomechanical modeling tools. This is followed by a more detailed review of geomechanical aspects, including reservoir stress-strain and microseismicity, well integrity, caprock sealing performance, and the potential for fault reactivation and notable (felt) seismic events. Geomechanical observations at current GCS field deployments, mainly at the In Salah CO 2 storage project in Algeria, are also integrated into the review. The In Salah project, with its injection into a relatively thin, low-permeability sandstone is an excellent analogue to the saline aquifers that might be used for large scale GCS in parts of Northwest Europe, the U.S. Midwest, and China. Some of the lessons learned at In Salah related to geomechanics are discussed, including how monitoring of geomechanical responses is used for detecting subsurface geomechanical changes and tracking fluid movements, and how such monitoring and geomechanical analyses have led to preventative changes in the injection parameters. Recently, the importance of geomechanics has become more widely recognized among GCS stakeholders, especially with respect to the potential for triggering notable (felt) seismic events and how such events could impact the long-term integrity of a CO 2 repository (as well as how it could impact the public perception of GCS). As described in the paper, to date, no notable seismic event has been reported from any of the current CO 2 storage projects, although some unfelt microseismic activities have been detected by geophones. However, potential future commercial GCS operations from large power plants will require injection at a much larger scale. In conclusion, for such large-scale injections, a staged, learn-as-you-go approach is recommended, involving a gradual increase of injection rates combined with continuous monitoring of geomechanical changes, as well as siting beneath a multiple layered overburden for multiple flow barrier protection, should an unexpected deep fault reactivation occur.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

An Uncertainty Quantification Framework for Studying the Effect of Spatial Heterogeneity in Reservoir Permeability on CO2 Sequestration

In this paper, we introduce an uncertainty quantification (UQ) software framework for carbon sequestration, focused on the effect of spatial heterogeneity of reservoir properties on CO2 migration. We use a sequential Gaussian method (SGSIM) to generate realizations of permeability fields with various spatial statistical attributes. To deal with the computational difficulties, we integrate the following ideas/approaches. First, we use three different sampling approaches (probabilistic collocation, quasi-Monte Carlo, and adaptive sampling) to reduce the number of forward calculations while trying to explore the parameter space and quantify the input uncertainty. Second, we use eSTOMP as the forward modeling simulator. eSTOMP is implemented with the Global Arrays toolkit that is based on one-sided inter-processor communication and supports a shared memory programming style on distributed memory platforms, providing a highly-scalable performance. Third, we built an adaptive system infrastructure to select the best possible data transfer mechanisms, to optimally allocate system resources to improve performance and to integrate software packages and data for composing carbon sequestration simulation, computation, analysis, estimation and visualization. We demonstrate the framework with a given CO2 injection scenario in heterogeneous sandstone reservoirs.

Carbon sequestration, uncertainty quantification, ↗
DOE OSTI2018

Methane Diffusion and Adsorption in Shale Rocks: A Numerical Study Using the Dusty Gas Model in TOUGH2/EOS7C-ECBM

Gas production from shale gas reservoirs plays a significant role in satisfying increasing energy demands. Compared with conventional sandstone and carbonate reservoirs, shale gas reservoirs are characterized by extremely low porosity, ultra-low permeability and high clay content. Slip flow, diffusion, adsorption and desorption are the primary gas transport processes in shale matrix, while Darcy flow is restricted to fractures. Understanding methane diffusion and adsorption, and gas flow and equilibrium in the low-permeability matrix of shale is crucial for shale formation evaluation and for predicting gas production. Modeling of diffusion in low-permeability shale rocks requires use of the Dusty gas model (DGM) rather than Fick’s law. The DGM is incorporated in the TOUGH2 module EOS7C-ECBM, a modified version of EOS7C that simulates multicomponent gas mixture transport in porous media. Also included in EOS7C-ECBM is the extended Langmuir model for adsorption and desorption of gases. In this study, a column shale model was constructed to simulate methane diffusion and adsorption through shale rocks. The process of binary CH 4 - N 2 diffusion and adsorption was analyzed. A sensitivity study was performed to investigate the effects of pressure, temperature and permeability on diffusion and adsorption in shale rocks. The results show that methane gas diffusion and adsorption in shale is a slow process of dynamic equilibrium, which can be illustrated by the slope of a curve in CH4mass variation. The amount of adsorption increases with the pressure increase at the low pressure, and the mass change by gas diffusion will decrease due to the decrease in the compressibility factor of the gas. With the elevated temperature, the gas molecules move faster and then the greater gas diffusion rates make the process duration shorter. The gas diffusion rate decreases with the permeability decrease, and there is a limit of gas diffusion if the permeability is less than 1.0×10 -15 m 2 . In conclusion, the results can provide insights for a better understanding of methane diffusion and adsorption in the shale rocks so as to optimize gas production performance of shale gas reservoirs.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2014

Exploring the effects of data quality, data worth, and redundancy of CO2 gas pressure and saturation data on reservoir characterization through PEST Inversion

This study examined the impacts of reservoir properties on CO2 migration after subsurface injection and evaluated the possibility of characterizing reservoir properties using CO2 monitoring data such as saturation distribution. The injection reservoir was assumed to be located 1400-1500 m below the ground surface such that CO2 remained in the supercritical state. The reservoir was assumed to contain layers with alternating conductive and resistive properties, which is analogous to actual geological formations such as the Mount Simon Sandstone unit. The CO2 injection simulation used a cylindrical grid setting in which the injection well was situated at the center of the domain, which extended up to 8000 m from the injection well. The CO2 migration was simulated using the PNNL-developed simulator STOMP-CO2e (the water-salt-CO2 module). We adopted a nonlinear parameter estimation and optimization modeling software package, PEST, for automated reservoir parameter estimation. We explored the effects of data quality, data worth, and data redundancy on the detectability of reservoir parameters using CO2 saturation monitoring data, by comparing PEST inversion results using data with different levels of noises, various numbers of monitoring wells and locations, and different data collection spacing and temporal sampling intervals. This study yielded insight into the use of CO2 saturation monitoring data for reservoir characterization and how to design the monitoring system to optimize data worth and reduce data redundancy.

Carbon sequestration↗
DOE OSTI2015

Rapid imbibition of water in fractures within unsaturated sedimentary rock

The spontaneous imbibition of water and other liquids into gas-filled fractures in variably-saturated porous media is important in a variety of engineering and geological contexts. However, surprisingly few studies have investigated this phenomenon. In this paper, we present a theoretical framework for predicting the 1-dimensional movement of water into air-filled fractures within a porous medium based on early-time capillary dynamics and spreading over the rough surfaces of fracture faces. The theory permits estimation of sorptivity values for the matrix and fracture zone, as well as a dispersion parameter which quantifies the extent of spreading of the wetting front. Quantitative data on spontaneous imbibition of water in unsaturated Berea sandstone cores were acquired to evaluate the proposed model. The cores with different permeability classes ranging from 50 to 500 mD and were fractured using the Brazilian method. Spontaneous imbibition in the fractured cores was measured by dynamic neutron radiography at the Neutron Imaging Prototype Facility (beam line CG-1D, HFIR), Oak Ridge National Laboratory. Water uptake into both the matrix and the fracture zone exhibited square-root-of-time behavior. The matrix sorptivities ranged from 2.9 to 4.6 mm s -0.5 , and increased linearly as the permeability class increased. The sorptivities of the fracture zones ranged from 17.9 to 27.1 mm s -0.5 , and increased linearly with increasing fracture aperture width. The dispersion coefficients ranged from 23.7 to 66.7 mm 2 s -1 and increased linearly with increasing fracture aperture width and damage zone width. Both theory and observations indicate that fractures can significantly increase spontaneous imbibition in unsaturated sedimentary rock by capillary action and surface spreading on rough fracture faces. Fractures also increase the dispersion of the wetting front. In conclusion, further research is needed to investigate this phenomenon in other natural and engineered porous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Mixing and trapping of dissolved CO 2 in deep geologic formations with shale layers

For dissolution trapping, the spatial variability of the geologic properties of naturally complex storage formations can significantly impact flow patterns and storage mechanisms of dissolved CO 2 . The significance of diffusive mixing that occurs in low permeability layers embedded between relatively higher permeability materials was highlighted by Agartan et al. (2015) using a highly controlled laboratory experimental study on trapping of dissolved CO 2 in multilayered systems. In this paper, we present a numerical modeling study on the impacts of low permeability layers on flow and storage of dissolved CO 2 in realistic field-scale settings. The simulator of variable-density flow used in this study was first verified using the experimental data in Agartan et al. (2015) to capture the observed processes. The simulator was then applied to a synthetic, field-scale multilayered system, with 19 sensitivity cases having variable permeability and thickness of the shale layers as well as the source strength and geometry of the source zone of dissolved CO 2 . Simulation results showed that the presence of continuous shale layers in the storage system disrupts the convective mixing by enhancing lateral spreading of dissolved CO 2 in sandstone layers and retarding the vertical mixing of dissolved CO 2 . The effectiveness of trapping of dissolved CO 2 depends on the physical properties of the shale layers and configurations of the source zone. The comparison to homogeneous cases with effective vertical permeability shows that it is important to capture these continuous thin shale layers in a storage formation and include them in the models to enhance dissolution trapping.

58 GEOSCIENCES↗
DOE OSTI2018

Transient analysis of advancing contact angle measurements on polished rock surfaces

We report that contact angle measurements for gas-liquid-rock systems are important for modeling multi-phase flow and transport in the subsurface. These data are needed in applications such as the extraction of oil and gas resources, geologic sequestration of carbon dioxide, contaminant fate and transport, and aquifer recharge through the vadose zone. Contact angles are frequently measured with the sessile drop method. Previous research has largely ignored the dynamic behavior of sessile droplets on geologic materials. This study investigates the dynamic behavior of sessile water droplets on prepared rock surfaces in the presence of air. Droplet diameter and advancing contact angle were determined at 0.5 s intervals for ~90 s on flat polished disks of Burlington limestone, Crossville sandstone, Mancos shale, Sierra White granite, Vermilion Bay granite, and Westerly granite using a Krüss DSA 30 Drop Shape Analyzer. The droplet diameter and advancing contact angle data sets were nonlinearly regressed against time using two different two-parameter models. The median coefficients of determination for the fits were 0.85 and 0.96, respectively. The resulting parameter estimates were used to compute the apparent equilibrium contact angle, θ e , for each disk following droplet diameter stabilization. Estimates of θ e ranged from 37.2° for Mancos shale to 75.6° for Burlington limestone. Analysis of variance indicated statistically significant differences in θ e between the rock types at the 95% confidence level. The variability of θe on the polished rock surfaces, as quantified by the coefficient of variation (CV) for θ e , varied between ~3 and ~9%; there were no significant differences in CV between the rock types. Neutron radiography indicated changes droplet morphology over time were due to the spontaneous imbibition of water into the rock matrix. Finally, the transient analysis employed in this study permits a more meaningful estimate of the equilibrium contact angle for rocks than taking the initial value or averaging over time as is frequently done.

58 GEOSCIENCES↗
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2023

A deep learning upscaling framework: Reactive transport and mineral precipitation in fracture-matrix systems

Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes for which effective medium approximations may become inaccurate, e.g., mineral precipitation and clogging during reactive transport. We have developed a deep learning upscaling framework, in which pore-scale modeling is directly employed in macroscale systems, without relying on effective medium approximations. The upscaling framework is first developed for general multiscale systems and then applied to modeling reactive transport with mineral precipitation in the altered layer in fracture-matrix structures. Solute transport from the fractures to the matrix is modeled as a wall boundary condition for the fractures, which, in turn, is predicted by recurrent neural networks using the concentration histories at the fracture-matrix boundary. Specifically, we consider a meter-scale fracture network embedded in sandstones, where the smallest feature is at the micron scale. Here the proposed framework allows us to span five orders of magnitude in length scales by capturing mineral precipitation in the altered layer of the rock matrix at the pore scale across the entire meter-scale fracture network.

42 ENGINEERING↗
DOE OSTI2013

Strontium isotope quantification of siderite, brine and acid mine drainage contributions to abandoned gas well discharges in the Appalachian Plateau

Unplugged abandoned oil and gas wells in the Appalachian region can serve as conduits for the movement of waters impacted by fossil fuel extraction. Strontium isotope and geochemical analysis indicate that artesian discharges of water with high total dissolved solids (TDS) from a series of gas wells in western Pennsylvania result from the infiltration of acidic, low Fe (Fe < 10 mg/L) coal mine drainage (AMD) into shallow, siderite (iron carbonate)-cemented sandstone aquifers. The acidity from the AMD promotes dissolution of the carbonate, and metal- and sulfate-contaminated waters rise to the surface through compromised abandoned gas well casings. Strontium isotope mixing models suggest that neither upward migration of oil and gas brines from Devonian reservoirs associated with the wells nor dissolution of abundant nodular siderite present in the mine spoil through which recharge water percolates contribute significantly to the artesian gas well discharges. Natural Sr isotope composition can be a sensitive tool in the characterization of complex groundwater interactions and can be used to distinguish between inputs from deep and shallow contamination sources, as well as between groundwater and mineralogically similar but stratigraphically distinct rock units. This is of particular relevance to regions such as the Appalachian Basin, where a legacy of coal, oil and gas exploration is coupled with ongoing and future natural gas drilling into deep reservoirs.

03 NATURAL GAS↗
DOE OSTI2016

Chemical and Sr isotopic characterization of North America uranium ores: Nuclear forensic applications

This study reports major, minor, and trace element data and Sr isotope ratios for 11 uranium ore (uraninite, UO 2+x ) samples and one processed uranium ore concentrate (UOC) from various U.S. deposits. The uraninite investigated represent ores formed via different modes of mineralization (e.g., high- and low-temperature) and within various geological contexts, which include magmatic pegmatites, metamorphic rocks, sandstone-hosted, and roll front deposits. In situ trace element data obtained by laser ablation-ICP-MS and bulk sample Sr isotopic ratios for uraninite samples investigated here indicate distinct signatures that are highly dependent on the mode of mineralization and host rock geology. Relative to their high-temperature counterparts, low-temperature uranium ores record high U/Th ratios (>1000), low total rare earth element (REE) abundances (<1 wt%), high contents (>300 ppm) of first row transition metals (Sc, Ti, V, Cr, Mn, Co, Ni), and radiogenic 87 Sr/ 86 Sr ratios (>0.7200). Comparison of chondrite normalized REE patterns between uraninite and corresponding processed UOC from the same locality indicates identical patterns at different absolute concentrations. Lastly, this result ultimately confirms the importance of establishing geochemical signatures of raw, uranium ore materials for attribution purposes in the forensic analysis of intercepted nuclear materials.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Enhancing predictive understanding and accuracy in geological carbon dioxide storage monitoring: Simulation and history matching of tracer transport dynamics

Co-injection of conservative tracers with carbon dioxide (CO 2 ) is a viable tool for monitoring subsurface processes during geological CO 2 storage (GCS). This research investigates the simulation and history-matching of a gas tracer (sulfur hexafluoride, SF 6 ) during CO 2 flooding, employing a core flooding result in Berea sandstone. Four extensively used saturation functions are assessed for their efficacy in history matching of CO 2 /SF 6 injection at the core scale. The history-matching process incorporates particle swarm optimization (PSO) to fine-tune constitutive relationships parameters. Next, employing transport models at the aquifer scale, we interrogate the impact on tracer transport and mixing of saturation function uncertainties, arising from the non-uniqueness of constitutive relationships parameters and saturation function types. To assess the effects of geological heterogeneity on behavior of tracer breakthrough curves (BTCs), we employ two normalized parameters assessing the degree of mixing and SF 6 breakthrough time. The aquifer-scale investigation encompasses both homogeneous and heterogeneous systems with and without capillary heterogeneity effects. Our findings underscore the critical importance of addressing saturation function uncertainties, emphasizing the significance of auxiliary experiments and innovative methodologies to enhance predictive accuracy. The findings highlight significant disparities in arrival times, BTC peaks, tails, and mixing levels, even under optimal conditions. Heterogeneity, with or without capillary heterogeneity, plays a crucial role in shaping BTC variations, resulting in accelerated SF 6 breakthrough times and reduced BTC peaks. Evaluation of monitoring points distant from the injector reveals a dampening effect on the SF6 BTC peak, particularly in heterogeneous systems with capillary heterogeneity, where the peak is halved. These insights underscore the challenges associated with tracer monitoring and the necessity for enhanced methodologies to improve predictive accuracy in subsurface processes during GCS.

58 GEOSCIENCES↗
DOE OSTI2020

Mobilization of trace metals from caprock and formation rocks at the Illinois Basin – Decatur Project demonstration site under geological carbon dioxide sequestration conditions

One concern for geologic CO 2 sequestration is the potential leakage of CO 2 or CO 2 -saturated brine containing trace metals into overlying aquifers, which poses the risk of adversely affecting underground sources of drinking water. In this work, rock and brine samples collected in the Illinois Basin – Decatur Project (IBDP), a large-scale geological CO 2 sequestration demonstration project in Decatur, Illinois, USA, as well as synthetic brine samples were used to understand the mobilization of trace metals as a result of CO 2 –rock–brine interactions under IBDP-specific conditions (50 °C and 20.6 MPa). Rock sample characterization indicated that, at the IBDP site, trace metal concentrations were greatest in cap rock samples from the Eau Claire Formation, compared with those from formations above and below this Formation. The natural brine samples collected from the Mt. Simon Sandstone reservoir were highly reducing and saline, with trace metal concentrations up to 680 times greater than the U.S. Environmental Protection Agency-prescribed drinking water standards. Batch leaching experiments indicated that both trace metal mobilization from rock and immobilization from the brine occurred when high-pressure CO 2 was introduced into the rock-brine system. The amounts of metals mobilized from the rock generally accounted for <5% of the total metals in the rock, but for some metals, including Ni, Pb, and Tl, up to 63% of the metals in the rock were mobilized to the brine. Leaching of trace metals into synthetic brines was different from that into the natural brine. The results of this study provided current information on the trace element sources and will help in risk simulations and experimental design to further evaluate potential impact of leakage on groundwater quality, which is important for future GCS projects to consider for monitoring and containment assurance purposes.

58 GEOSCIENCES↗
DOE OSTI2026

Geochemical evolution in Cacapon member: Fluid-rock interaction experiments and model insights for Appalachian Basin geothermal development

Here, this study combines recirculated flow-through experimental results conducted for 17 days at 90C and 200 PSI with reactive transport modeling to estimate fluid-rock interactions occurring in a sandy mudstone using an interbedded sandstone-shale sample from the Cacapon Member of the upper Tuscarora Sandstone/lower Rose Hill Formation for the purpose of geothermal exploration. Results suggest that the fluid and rock are likely to be in or near partial-equilibrium after approximately one year. In addition, after >400 h of continuous injection at 0.05 mL/min (5*10 −8 m 3 /min) the reactive front is restricted to the first ∼13% of the 4 cm*1.6 cm experimental rock length, whereas after >9100 h (∼1 year), the reactive front extends to 30% total length. The rate of changes in dissolution or precipitation are however, very minimal, with all major rock-forming minerals having rates <10 −11 mol/L porous media/s. Reservoir rock in the presence of dilute brine as may be the case during operation of an enhanced geothermal system would experience little alteration during the shut-in phase, and possibly up to one year. These results have utility in geothermal exploration for reservoirs at similar temperatures as well as general fluid-sandy mudstone rock interaction geochemistry.

Appalachian Basin↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗