Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Granular medium”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The (mis)use of subjective process measures in software engineering

A variety of measures are used in software engineering research to develop an understanding of the software process and product. These measures fall into three broad categories: quantitative, characteristics, and subjective. Quantitative measures are those to which a numerical value can be assigned, for example effort or lines of code (LOC). Characteristics describe the software process or product; they might include programming language or the type of application. While such factors do not provide a quantitative measurement of a process or product, they do help characterize them. Subjective measures (as defined in this study) are those that are based on the opinion or opinions of individuals; they are somewhat unique and difficult to quantify. Capturing of subjective measure data typically involves development of some type of scale. For example, 'team experience' is one of the subjective measures that were collected and studied by the Software Engineering Laboratory (SEL). Certainly, team experience could have an impact on the software process or product; actually measuring a team's experience, however, is not a strictly mathematical exercise. Simply adding up each team member's years of experience appears inadequate. In fact, most researchers would agree that 'years' do not directly translate into 'experience.' Team experience must be defined subjectively and then a scale must be developed e.g., high experience versus low experience; or high, medium, low experience; or a different or more granular scale. Using this type of scale, a particular team's overall experience can be compared with that of other teams in the development environment. Defining, collecting, and scaling subjective measures is difficult. First, precise definitions of the measures must be established. Next, choices must be made about whose opinions will be solicited to constitute the data. Finally, care must be given to defining the right scale and level of granularity for measurement.

Valett, Jon D.↗
NASA NTRS2000

Shear Band Formation in Plane Strain Experiments of Sand

A series of biaxial (plane strain) experiments were conducted on three sands under low (15 kPa) and high (100 kPa) confining pressure conditions to investigate the effects of specimen density, confining pressure, and sand grain size and shape on the constitutive and stability behavior of granular materials. The three sands used in the experiments were fine-, medium-, and coarse-grained uniform silica sands with rounded, subangular, and angular grains, respectively. Specimen deformation was readily monitored and analyzed with the help of a grid pattern imprinted on the latex membrane. The overall stress-strain behavior is strongly dependent on the specimen density, confining pressure, sand grain texture, and the resulting failure mode(s). That became evident in different degrees of softening responses at various axial strains. The relationship between the constitutive behavior and the specimens' modes of instability is presented. The failure in all specimens was characterized by two distinct and opposite shear bands. It was found that the measured dilatancy angles increase as the sand grains' angularities and sizes increase. The measured shear band inclination angles are also presented and compared with classical Coulomb and Roscoe solutions.

Alshibli, Khalid A.↗
NASA NTRS2000

Assessment of Constitutive and Stability Behavior of Sands Under Plane Strain Condition

A series of biaxial (plane strain) experiments were conducted on three sands under low (15 kPa) and high (100 kPa) confining pressure conditions to investigate the effects of specimen density, confining pressure, and sand grains size and shape on the constitutive and stability behavior of granular materials. The three sands used in the experiments were fine, medium, and coarse-grained uniform silica sands with rounded, sub-angular, and angular grains, respectively. Specimen deformation was readily monitored and analyzed with the help of a grid pattern imprinted on the latex membrane. The overall stress-strain behavior is strongly dependent on the specimen density, confining pressure, sand grain texture, and the resulting failure mode(s). That became evident in different degrees of softening responses at various axial strains. The relationship between the constitutive behavior and the specimens' modes of instability is presented. The failure in all specimens was characterized by two distinct and opposite shear bands. It was found that the measured dilatancy angles increase as the sand grains' angularity and size increase. The measured shear band inclination angles are also presented and compared with classical Coulomb and Roscoe solutions.

Alshibli, Khalid A.↗
NASA NTRS1996

Capillary movement of liquid in granular beds in microgravity

A more complete understanding of the dynamics of capillary flow through an unsaturated porous medium would be useful for the development of an effective water and nutrient delivery system for the growth of plants in space. An experiment was conducted on the Mir Space Station that used an experimental cuvette called "Capillary Test Bed" to compare fluid migration under terrestrial laboratory conditions by positioning the cuvette such that the hydrostatic force is negated and on Mir under microgravity conditions. Differences in fluid migration in the cuvette were observed with migration being slower in microgravity compared with some ground control experiments.

Mir Project↗
NASA NTRS1994

Capillary Movement of Liquid in Granular Beds in Microgravity

A more complete understanding of the dynamics of capillary flow through an g unsaturated porous medium would be useful for a number of space and terrestrial applications. One such application involves the development of an effective water and nutrient delivery system for the growth of plants in space. An experiment was conducted on the Mir Space Station that used an experimental cuvette called "Capillary Test Bed". An objective of this experiment was to compare fluid migration in the "Capillary Test Bed" under terrestrial laboratory simulated microgravity conditions by positioning the cuvette such that the hydrostatic force is negated and on Mir at microgravity conditions. Difference in fluid migration in the cuvette were observed. Based on these observations, improvements in the design of the cuvette are under-way that should allow for collection of more precise data in future flight experiments. The results of the MIR experiment and modifications being made to the cuvette for improved data collection will be discussed.

Yendeler, Boris S.↗
DOE OSTI2025

Building a Healthy Urban Design Index (HUDI): how to promote health and sustainability in European cities

As global urbanisation accelerates, alongside declining environmental quality and increasing climate challenges, it is increasingly vital for urban planners and policy makers to integrate health and wellbeing considerations into urban planning. This study introduces the Healthy Urban Design Index (HUDI), a high-resolution spatial index developed for European cities. HUDI combines policy-relevant indicators related to urban design, sustainable transportation, environmental quality, and greenspace accessibility—key factors influencing human health and well-being. Unlike existing indices, which often focus on few or large metropolitan cities and lack spatial granularity, HUDI offers high resolution and extends its scope to small-sized and medium-sized cities, home to over 50% of Europe's population.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2016

Modeling and Simulation of Trapping Mechanisms of Granular Media in Space

This paper describes the modeling and simulation of trapped granular media, within the context of the Granular Imager project. We describe the physics of trapped granular media in space, and the methodologies used to stably confine and shape such a medium using electromagnetic fields. The numerical models have also been validated with results in the literature, obtaining excellent agreement. The results of the numerical tests indicate that it is possible, with structural arrangements of rings and plates at different levels of electrostatic potential, to stably confine one or more charged particles, when driven by voltages that can be modulated in time and space.

Quadrelli, Marco B.↗
NASA NTRS1996

Capillary Movement in Substrates in Microgravity

A more complete understanding of the dynamics of capillary flow through an unsaturated porous medium would be useful for a number of space and terrestrial applications. Knowledge of capillary migration of liquids in granular beds in microgravity would significantly enhance the development and understanding of how a matrix based nutrient delivery system for the growth of plants would function in a microgravity environment. Thus, such information is of interest from the theoretical as well as practical point of view.

Bula, R. J.↗
NASA NTRSDate not supplied

Effect of Vacuum on Force Response of an Ultrasonic Penetrator

Introduction: The Apollo astronauts encountered higher than expected resistances when interacting with the lunar soil via the Apollo Lunar Surface Drill (ALSD) and the trenching tool. Reducing the force required to move tools or other mechanical components through regolith will impact many steps of the resource extraction process. Force reduction has been achieved in soil materials by imparting vibration to tooling interfaces such as a vibratory farming cultivator, a percussive scoop, and ultrasonically resonant penetrators. Vibration-assisted tools in granular media reduce interaction forces by fluidizing a volume around the tool, allowing the tool to progress through a dynamic (fluid) medium instead of a static (solid) medium. This work seeks to quantify ultrasonic vibration’s effect on the force response of a penetrator in lunar soil simulant in vacuum sufficient to be within the molecular flow regime of any disturbed gases. Methods: A custom vacuum chamber setup, CUBEvac, was designed and built to facilitate penetration testing in a high vacuum environment, for comparison to penetration behavior in ambient terrestrial environment. A two-stage pumping system (Agilent Triscroll 600 roughing pump, Agilent VHS-6 oil diffusion pump) reached chamber pressures of about 5x10-6 Torr with regolith simulant in place. Figure 1 is a schematic of the heart of the assembly (note the penetration drive mechanisms above the chamber feedthrough and the regolith simulant sample in the bottom are not shown). The penetration actuation stack was comprised of a stepper motor driving a lead screw to move the ultra-sonic probe vertically inside the chamber. Motion was coordinated with an Arduino Uno. GRC-3 lunar simulant was used for this set of experiments. Samples were prepared in a four-liter stain-less steel, cylindrical pot with an internal diameter of approximately 15.56 cm (6.125 in) and a depth of 19.37 cm (7 5/8 in) for testing. The soil was baked out prior to compaction preparation as a measure to reduce soil moisture which interfered with pump down capacity. The soil was not baked again if it was removed from the vacuum chamber, prepped, and immediately returned to the vacuum chamber for pump down. The soil was compacted using a 60 Hz vibration table with a surcharge of 34 kg place on top of the soil in the container. Prepared soil samples weighed approximately 6.5 kg (bulk density 1.895-1.934 g/cm3). Two probe end effectors were tested: A cone penetrometer (static only) with a nominal diameter of 12.7 mm (0.5 in) and a nominal height of 28.6 mm (1.125 in); and a vibrating cylindrical probe measuring 12.7 mm in diameter and 50.8 mm in effective length from the tip (Figure 2). The cylindrical probe vibrated resonantly at 20 kHz with an amplitude 23 μm. The cone penetration tests were conducted to assess potential soil behavior differences in vacuum. The cylinder probe tests were conducted as the primary subject of this investigation to assess force response in vacuum. For each test, a regolith simulant sample was load-ed and compacted in the chamber, which was then evacuated for roughly 18 hours to reach the lowest possible pressure (approximately 5x10-6 Torr for most tests). The probe was then moved to about 10 mm above the soil surface before being pushed to a depth of 50 mm for the cylinder probe tests and to a depth of 100 mm for the cone penetrometer tests, both at a speed of 2 mm/s. The simulant samples were prepared the same for all tests. Ideally, they would respond consistently to probe penetration under ambient and vacuum conditions. This was evaluated by measuring the resistance of representative prepared simulant beds with a standard cone penetrometer in both environments. Results and Discussion: The resistance of the simulant samples in the vacuum tests was consistently lower than in the ambient tests as determined by the cone penetration tests. Thus, the ambient and vacuum results cannot be compared directly; work is underway to de-confound and better correlate the data. Still, figures 3-6 show that probe penetration forces are lower overall in the vacuum environment. In both environments, resonant vibration of the probe provides two useful effects: It reduces the probe penetration force and smooths the force-depth curve, significantly reducing local maxima. These effects have implications for various potential applications, such as astronaut hand-tools, where benefits (reducing astronaut effort) outweigh the cost of the additional energy re-quired to generate vibration. These results demonstrate that resonantly vibrating tools can meaningfully reduce the penetration force required for excavation, probing, and drilling tools in simulated lunar regolith deposits under vacuum levels approaching those that will be experience on the Moon’s surface. Lunar-gravity, ambient environment tests are scheduled soon. The effects of realistic temperatures and temperature gradients and deeper vacuum remain to be tested.

E Rezich↗
DOE OSTI2024

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2020

Spatial U-Pb Age Distribution in Shock-Recrystallized Zircon – A Case Study from the Rochechouart Impact Structure, France

Age determination of impact structures via the zircon U–Pb system remains challenging and often ambiguous due to highly variable effects of shock metamorphism on U-Pb geochronology. It is, therefore, crucial to link the observed zircon microtextures, including their temperature and pressure conditions associated with their formation, directly to the U–Pb ages preserved. Here, we analyzed three recrystallized zircon grains and one plastically deformed zircon crystal from the medium-sized Rochechouart impact structure in the northwestern Massif Central of France. For the Rochechouart impact structure the impact age (206.92 ± 0.32 Ma [40Ar/39Ar]) as well as the 2 tectono-themal history is well established making this study site ideal to test concepts about U-Pb systematics in shocked zircon and to differentiate between shock-driven age resetting and preimpact crystallization and metamorphic overprinting. Zircon microstructures were studied using scanning electron imaging, cathodoluminescence imaging, and electron backscatter diffraction (EBSD) mapping. Further, we conducted U–Pb Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) depth-profiling analysis, allowing us to interpret the resultant age data in discrete steps with increasing ablation time/depth. The U–Pb depth profiling data demonstrate that plastically strained grains are incompletely reset and preserve rim and interior age domains reflecting typical pre-impact (pre-Variscan and Variscan) regional tectonic ages. Our results also reveal that the granular crystals encountered contain microstructural evidence for “former reidite in granular neoblastic” (FRIGN) zircon, reflecting both high-pressure (⩾30GPa) and high temperature (⩾1200℃) conditions. This signifies that FRIGN zircon is now known from an additional, medium sized, impact structure further supporting the hypothesis that this impact induced microstructure is commonly preserved. In addition, FRIGN zircon has a high potential to preserve completely impact-reset crystal domains (~204 to 207 Ma) that can be identified by our combined analytical approach of U–Pb depth profiling and EBSD mapping, and thus are suitable for determining reliable impact ages

U-Pb depth profiling↗
NASA NTRS1978

Balloon-borne imagery of the solar granulation. II - The lifetime of solar granulation

Phenomenological aspects of the temporal evolution of photospheric granulation are reported as derived from time series of granulation photographs obtained during a flight of a balloon-borne telescope. The distribution of granule lifetime probabilities is determined, and it is found that the data can be represented by an exponential decrease with a 'decay constant' of 5.9 min. The general properties of granular evolution are described along with the way individual granules evolve with time. The most common type of granule is shown to be a medium-sized or small fragment, and it is suggested that all granules are produced by fragmentation of preexisting granules. The relative frequencies of granule destruction by fragmentation, fading, and merging are determined to be 51%, 21%, and 28%, respectively. An average radial velocity of 0.8 km/s is computed for conglomerates with an average diameter of 2.25 arcsec.

Mehltretter, J. P.↗
NASA NTRS2007

Ejecta- and Size-Scaling Considerations from Impacts of Glass Projectiles into Sand

One of the most promising means of learning how initial impact conditions are related to the processes leading to the formation of a planetary-scale crater is through scaling relationships.1,2,3 The first phase of deriving such relationships has led to great insight into the cratering process and has yielded predictive capabilities that are mathematically rigorous and internally consistent. Such derivations typically have treated targets as continuous media; in many, cases, however, planetary materials represent irregular and discontinuous targets, the effects of which on the scaling relationships are still poorly understood.4,5 We continue to examine the effects of varying impact conditions on the excavation and final dimensions of craters formed in sand. Along with the more commonly treated variables such as impact speed, projectile size and material, and impact angle,6 such experiments also permit the study of changing granularity and friction angle of the target materials. This contribution presents some of the data collected during and after the impact of glass spheres into a medium-grained sand.

Anderson J. L. B.↗
DOE OSTI2023

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2023

On the applicability of various levels of detail for occupant behavior representation and modeling in building performance simulation

Occupant behavior (OB) is one of the significant sources of uncertainty in building performance simulation. While OB modeling has received increased attention in the past decade, research on the degree of granularity or level of detail (LoD) required for representing occupants is still in the nascent stages. This paper analyzes the modeling and applicability of three LoDs to represent occupants in building performance assessment. A medium-sized prototype office building located in Chicago, Illinois is used as the simulation case study. Ten occupant-centric attributes are adopted to develop the LoDs for OB representation. We first demonstrate the different modeling approaches required for simulating the three fidelity levels. Later, we illustrate the suitability of the developed LoDs in supporting six building performance use cases across different lifecycle stages. Furthermore, this study intends to provide guidance for the building simulation community on appropriate OB representation to support various use cases.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Performance of the CMS High Granularity Calorimeter prototype to charged pion beams of 20–300 GeV/c

The upgrade of the CMS experiment for the high luminosity operation of the LHC comprises the replacement of the current endcap calorimeter by a high granularity sampling calorimeter (HGCAL). The electromagnetic section of the HGCAL is based on silicon sensors interspersed between lead and copper (or copper tungsten) absorbers. The hadronic section uses layers of stainless steel as an absorbing medium and silicon sensors as an active medium in the regions of high radiation exposure, and scintillator tiles directly read out by silicon photomultipliers in the remaining regions. As part of the development of the detector and its readout electronic components, a section of a silicon-based HGCAL prototype detector along with a section of the CALICE AHCAL prototype was exposed to muons, electrons and charged pions in beam test experiments at the H2 beamline at the CERN SPS in October 2018. The AHCAL uses the same technology as foreseen for the HGCAL but with much finer longitudinal segmentation. The performance of the calorimeters in terms of energy response and resolution, longitudinal and transverse shower profiles is studied using negatively charged pions, and is compared to GEANT4 predictions. This is the first report summarizing results of hadronic showers measured by the HGCAL prototype using beam test data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗