SEARCH · Search NASA

Results for “sample complexity analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces – part II

The modeling of the near-field interaction in the scattering-type scanning near-field optical microscope (s-SNOM) is rapidly advancing, although an accurate yet versatile modeling framework that can be easily adapted to various complex situations is still lacking. In this work, we propose a time-efficient numerical scheme in the quasi-electrostatic limit to capture the tip-sample interaction in the near field. This method considers an extended tip geometry, which is a significant advantage compared to the previously reported method based on the point-dipole approximation. Using this formalism, we investigate, among others, nontrivial questions such as uniaxial and biaxial anisotropy in the near-field interaction, the relationship between various experimental parameters (e.g. tip radius, tapping amplitude, etc.), and the tip-dependent spatial resolution. The demonstrated method further sheds light on the understanding of the contrast mechanism in s-SNOM imaging and spectroscopy, while also representing a valuable platform for future quantitative analysis of the experimental observations.

Chen, Xinzhong (ORCID:0000000261033848)↗

NASA NTRS2014

Radiation Quality Effects on Transcriptome Profiles in 3-d Cultures After Particle Irradiation

In this work, we evaluate the differential effects of low- and high-LET radiation on 3-D organotypic cultures in order to investigate radiation quality impacts on gene expression and cellular responses. Reducing uncertainties in current risk models requires new knowledge on the fundamental differences in biological responses (the so-called radiation quality effects) triggered by heavy ion particle radiation versus low-LET radiation associated with Earth-based exposures. We are utilizing novel 3-D organotypic human tissue models that provide a format for study of human cells within a realistic tissue framework, thereby bridging the gap between 2-D monolayer culture and animal models for risk extrapolation to humans. To identify biological pathway signatures unique to heavy ion particle exposure, functional gene set enrichment analysis (GSEA) was used with whole transcriptome profiling. GSEA has been used extensively as a method to garner biological information in a variety of model systems but has not been commonly used to analyze radiation effects. It is a powerful approach for assessing the functional significance of radiation quality-dependent changes from datasets where the changes are subtle but broad, and where single gene based analysis using rankings of fold-change may not reveal important biological information. We identified 45 statistically significant gene sets at 0.05 q-value cutoff, including 14 gene sets common to gamma and titanium irradiation, 19 gene sets specific to gamma irradiation, and 12 titanium-specific gene sets. Common gene sets largely align with DNA damage, cell cycle, early immune response, and inflammatory cytokine pathway activation. The top gene set enriched for the gamma- and titanium-irradiated samples involved KRAS pathway activation and genes activated in TNF-treated cells, respectively. Another difference noted for the high-LET samples was an apparent enrichment in gene sets involved in cycle cycle/mitotic control. It is plausible that the enrichment in these particular pathways results from the complex DNA damage resulting from high-LET exposure where repair processes are not completed during the same time scale as the less complex damage resulting from low-LET radiation.

Patel, Z. S.↗

NASA NTRS2020

Thermal system and environmental testing of the Mars Helicopter

The Mars Helicopter will be a technology demonstration conducted during the Mars 2020 mission. The primary mission objective is to achieve several 90-second flights demonstrating the feasibility of heavier than air flight on Mars and capture visible light images via forward and nadir mounted cameras. These flights could possibly provide reconnaissance data for sampling site selection for other Mars surface missions. A solar array and batteries for flight operations, imaging, communications, and survival heating power the Helicopter. The thermal design is driven by minimizing survival heater energy while maintaining compliance with allowable flight temperatures in a variable thermal environment. Due to the small size of the Helicopter and its complex geometries, along with the fact that it operates with very low power and small margins in the extreme Mars environment, additional care had to be paid while planning thermal tests and designing the thermal system. The first section of the paper describes the evolution of the thermal system of the Mars Helicopter. After the first thermal vacuum test of the engineering model, the thermal team has conducted a partial effect analysis on the thermal design components that had a major impact on the system performance. Several design choices derived by analysis and test have been made to meet the energy allocation and the temperature requirements. These changes included increased gas gaps to reduce gas conduction, low emissivity coatings for internal components, blanket implementation, optimization of wire routing and fine-tuning of surface operations to optimize waste heat recovery. The second part of the paper describes the flight model thermal vacuum test and the subsequent thermal model correlation necessary to confirm the fidelity of the analysis results.

Pauken, Michael↗

DOE OSTI2020

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Bayesian inference in band excitation scanning probe microscopy for optimal dynamic model selection in imaging

The universal tendency in scanning probe microscopy (SPM) over the last two decades is to transition from simple 2D imaging to complex detection and spectroscopic imaging modes. The emergence of complex SPM engines brings forth the challenge of reliable data interpretation, i.e., conversion from detected signals to descriptors specific to tip–surface interactions and subsequently to material’s properties. In this work, we implemented a Bayesian inference approach for the analysis of the image formation mechanisms in band excitation SPM. Compared to the point estimates in classical functional fit approaches, Bayesian inference allows for the incorporation of extant knowledge of materials and probe behavior in the form of corresponding prior distribution and return the information on the material functionality in the form of readily interpretable posterior distributions. We explore the nonlinear mechanical behaviors spatially in a classical ferroelectric material, PbTiO 3 . We observe the non-trivial evolution of the Duffing stiffness term and the nonlinearity of the sample surface, determine spatial clustering of the nonlinear response, and perform a Landau analysis on predicting the nonlinear coefficient, which indicates that ferroelectric behavior can be a cause of the observed results. These observations suggest that the spectrum of anomalous behaviors at the ferroelectric domain walls may be broader than previously believed and can extend to non-conventional mechanical properties in addition to static and microwave conductance.

36 MATERIALS SCIENCE↗

DOE OSTI2022

First comprehensive identification of cardiac proteins with putative increased O-GlcNAc levels during pressure overload hypertrophy

Protein posttranslational modifications (PTMs) by O-GlcNAc globally rise during pressure-overload hypertrophy (POH). However, a major knowledge gap exists on the specific proteins undergoing changes in O-GlcNAc levels during POH primarily because this PTM is low abundance and easily lost during standard mass spectrometry (MS) conditions used for protein identification. Methodologies have emerged to enrich samples for O-GlcNAcylated proteins prior to MS analysis. Accordingly, our goal was to identify the specific proteins undergoing changes in O-GlcNAc levels during POH. We used C57/Bl6 mice subjected to Sham or transverse aortic constriction (TAC) to create POH. From the hearts, we labelled the O-GlcNAc moiety with tetramethylrhodamine azide (TAMRA) before sample enrichment by TAMRA immunoprecipitation (IP). We used LC-MS/MS to identify and quantify the captured putative O-GlcNAcylated proteins. We identified a total of 700 putative O-GlcNAcylated proteins in Sham and POH. Two hundred thirty-three of these proteins had significantly increased enrichment in POH over Sham suggesting higher O-GlcNAc levels whereas no proteins were significantly decreased by POH. We examined two MS identified metabolic enzymes, CPT1B and the PDH complex, to validate by immunoprecipitation. We corroborated increased O-GlcNAc levels during POH for CPT1B and the PDH complex. Enzyme activity assays suggests higher O-GlcNAcylation increases CPT1 activity and decreases PDH activity during POH. In summary, we generated the first comprehensive list of proteins with putative changes in O-GlcNAc levels during POH. Our results demonstrate the large number of potential proteins and cellular processes affected by O-GlcNAc and serve as a guide for testing specific O-GlcNAc-regulated mechanisms during POH.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Multi-Fidelity Stochastic Economic Dispatch for Operating Low Carbon Power Grids

Grid operators can address the inherently stochastic nature of renewables by solving a two-stage stochastic programming model that minimizes the cost of dispatch decisions while accounting for the complex grid dynamics. It is common to use a sample average approximation to estimate the expectation of the second stage costs in this model. However, the large sample count needed for numerical accuracy makes effective modeling large-scale electric grids computationally intractable. We introduce a control variate multi-fidelity estimator for the second-stage recourse that enables high quality dispatch decisions in real-time with a reduced computational burden. We obtain a hierarchy of model fidelities by linearizing the AC power flow system representation to DC power flow, and by relaxing transmission and voltage network constraints. We evaluate the performance of our proposed method on a synthetic grid with 73 buses against a deterministic baseline with persistence forecast and a high-fidelity reference. Our analysis shows a computational speed-up of 7.62x with a minimal loss in accuracy. The multi-fidelity method is well suited to fidelity combinations that use a simplified network topology in their lower fidelity model and is an attractive option for applications where accurate grid modeling needed on a limited computational budget.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS2016

Science Enabling Exploration: Using LRO to Prepare for Future Missions

Discoveries from LRO have transformed our understanding of the Moon, but LRO's instruments were originally designed to collect the measurements required to enable future lunar surface exploration. A high lunar exploration priority is the collection of new samples and their return to Earth for comprehensive analysis. The importance of sample return from South Pole-Aitken is well-established [Jolliff et al., this conference], but there are numerous other locations where sample return will yield important advances in planetary science. Using new LRO data, we have defined an achievability envelope based on the physical characteristics of successful lunar landing sites. Those results were then used to define 1km x 1km regions of interest where sample return could be executed, including: the basalt flows in Oceanus Procellarum (22.1N, 53.9W), the Gruithuisen Domes (36.1N, 39.7W), the Dewar cryptomare (2.2S, 166.8E), the Aristarchus pyroclastic deposit (24.8N, 48.5W), the Sulpicius Gallus formation (19.9N, 10.3E), the Sinus Aestuum pyroclastic deposit (5.2N, 9.2W), the Compton-Belkovich volcanic complex (61.5N, 99.9E), the Ina Irregular Mare Patch (18.7N, 5.3E), and the Marius Hills volcanic complex (13.4N, 55.9W). All of these locations represent safe landing sites where sample returns are needed to advance our understanding of the evolution of the lunar interior and the timescales of lunar volcanism. If LRO is still active when any future mission reaches the surface, LRO's capability to rapidly place surface activities into broader geologic context will provide operational advantages. LRO remains a unique strategic asset that continues to address the needs of future missions.

Lawrence, S. J.↗

NASA NTRS2002

Rapid, specific determination of iodine and iodide by combined solid-phase extraction/diffuse reflectance spectroscopy

A new, rapid methodology for trace analysis using solid-phase extraction is described. The two-step methodology is based on the concentration of an analyte onto a membrane disk and on the determination by diffuse reflectance spectroscopy of the amount of analyte extracted on the disk surface. This method, which is adaptable to a wide range of analytes, has been used for monitoring ppm levels of iodine and iodide in spacecraft water. Iodine is used as a biocide in spacecraft water. For these determinations, a water sample is passed through a membrane disk by means of a 10-mL syringe that is attached to a disk holder assembly. The disk, which is a polystyrene-divinylbenzene composite, is impregnated with poly(vinylpyrrolidone) (PVP), which exhaustively concentrates iodine as a yellow iodine-PVP complex. The amount of concentrated iodine is then determined in only 2 s by using a hand-held diffuse reflectance spectrometer by comparing the result with a calibration curve based on the Kubelka-Munk function. The same general procedure can be used to determine iodide levels after its facile and exhaustive oxidation to iodine by peroxymonosulfate (i.e., Oxone reagent). For samples containing both analytes, a two-step procedure can be used in which the iodide concentration is calculated from the difference in iodine levels before and after treatment of the sample with peroxymonosulfate. With this methodology, iodine and iodide levels in the 0.1-5.0 ppm range can be determined with a total workup time of approximately 60 s with a RSD of approximately 6%.

Non-NASA Center↗

NASA NTRS2003

Syringe filtration methods for examining dissolved and colloidal trace element distributions in remote field locations

It is well-established that sampling and sample processing can easily introduce contamination into dissolved trace element samples if precautions are not taken. However, work in remote locations sometimes precludes bringing bulky clean lab equipment into the field and likewise may make timely transport of samples to the lab for processing impossible. Straightforward syringe filtration methods are described here for collecting small quantities (15 mL) of 0.45- and 0.02-microm filtered river water in an uncontaminated manner. These filtration methods take advantage of recent advances in analytical capabilities that require only small amounts of waterfor analysis of a suite of dissolved trace elements. Filter clogging and solute rejection artifacts appear to be minimal, although some adsorption of metals and organics does affect the first approximately 10 mL of water passing through the filters. Overall the methods are clean, easy to use, and provide reproducible representations of the dissolved and colloidal fractions of trace elements in river waters. Furthermore, sample processing materials can be prepared well in advance in a clean lab and transported cleanly and compactly to the field. Application of these methods is illustrated with data from remote locations in the Rocky Mountains and along the Yukon River. Evidence from field flow fractionation suggests that the 0.02-microm filters may provide a practical cutoff to distinguish metals associated with small inorganic and organic complexes from those associated with silicate and oxide colloids.

Colloids/analysis↗

NASA NTRS2023

Soluble Organic Matter Molecular Atlas of Ryugu Reveals Cold Hydrothermalism on C-Type Asteroid Parent Body

The sample from the near-Earth carbonaceous asteroid (162173) Ryugu is analyzed in the context of carbonaceous meteorites soluble organic matter. The analysis of soluble molecules of samples collected by the Hayabusa2 spacecraft shines light on an extremely high molecular diversity on the C-type asteroid. Sequential solvent extracts of increasing polarity of Ryugu samples are analyzed using mass spectrometry with complementary ionization methods and structural information confirmed by nuclear magnetic resonance spectroscopy. Here we show a continuum in the molecular size and polarity, and no organomagnesium molecules are detected, reflecting a low temperature and water-rich environment on the parent body approving earlier mineralogical and chemical data. High abundance of sulfidic and nitrogen rich compounds as well as high abundance of ammonium ions confirm the water processing. Polycyclic aromatic hydrocarbons are also detected in a structural continuum of carbon saturations and oxidations, implying multiple origins of the observed organic complexity, thus involving generic processes such as earlier carbonization and serpentinization with successive low temperature aqueous alteration.

Asteroids, comets and Kuiper belt↗

DOE OSTI2019

Structural Characterization of a Series of N5-Ligated Mn IV -Oxo Species

Not provided.

Chemistry↗

DOE OSTI2022

Report on PTT Imaging of Defects in AM Metallic Materials-Part 2

36 MATERIALS SCIENCE↗

DOE OSTI2022

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials and Components (Final Technical Report)

Metal Additive Manufacturing (AM) is a promising method for cost-efficient fabrication of complex shape structures for applications in harsh environment, such as in a nuclear reactor. However, internal defects (pores) occur in high-strength AM alloys, which are manufactured with Laser Powder Bed Fusion (LPBF) AM method. Pulsed Infrared Thermography (PIT) is an efficient nondestructive evaluation (NDE) method to examine actual structures, because this method offers one-sided non-contact measurements, and fast processing of large sample areas. However, imaging of material defects, particularly defects with sizes at microscopic level, is challenging. In this report, we benchmark the performance of several Unsupervised Learning (UL) algorithms designed to enhance imaging of microscopic defects in metals with PIT. UL aims to learn the latent principal patterns (dictionaries) in PIT data to detect defects with minimal human supervision. Performance of Independent Component Analysis (ICA), Sparse Coding (SC), Principal Component Analysis (PCA) and Exploratory Factor Analysis (EFA) was compared using F-score, UL model training time and defects reconstruction time. We obtained the average F-score of 0.75, and a highest F-score of 0.89 for the EFA algorithm. Overall, EFA outperforms other UL algorithms considered in this study. In another approach, we investigate Thermal Tomography (TT), which is a computational method for reconstruction of depth profile of internal material defects from PIT nondestructive evaluation (NDE). TT algorithm obtains depth reconstructions of thermal effusivity, which has been shown to provide visualization of subsurface internals defects in metals. In many applications, one needs to determine the defect shape and orientation from reconstructed effusivity images. Interpretation of TT images is non-trivial because of blurring, which increases with depth due to heat diffusion-based nature of image formation. We have developed a deep learning convolutional neural network (CNN) to classify size and orientation of subsurface material defects in TT images. CNN was trained with TT images produced with computer simulations of 2D metallic structures (thin plates) containing elliptical subsurface voids. Performance of CNN was investigated using test TT images developed with computer simulations of plates containing elliptical defects, and defects with shape imported from scanning electron microscopy (SEM) images. CNN demonstrated the ability to classify radii and angular orientation of elliptical defects in previously unseen test TT images. We have also demonstrated that CNN trained on TT images of elliptical defects is capable of classifying shape and orientation of irregular defects. Training the CNN on irregular defect shapes instead of on elliptical shapes would make the resulting classifications more descriptive of actual defect shapes. However, this requires a much higher volume of SEM images of material defects, which are difficult to obtain because of random occurrence of defects in LPBF. To address this challenge, we developed a generative adversarial network (GAN) to augment the existing dataset of SEM defect images. The GAN model is demonstrated to create novel yet realistic defect shapes that can be used as input for simulated PTT images to train CNN. We also investigate several approaches based on Gaussian Random Circle and Bezier Curves for constructing parametric models of irregular-shape defects.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Analysis of Organic Electron Materials using Neutron Scattering (Final Report)

Neutron scattering is a powerful suite of measurement techniques that enable unique analysis of structure and dynamics of a wide variety of materials. Specifically, neutrons interact strongly with protons, which enables excellent contrast in light element organic materials that cannot be obtained with any other measurements. This series of three proposals focused on using inelastic neutron scattering (INS) to study the vibrations (dynamics) of the organic molecules involved in organic electronic devices. Our main goals were to demonstrate that INS is capable of elucidating new, useful, and unique knowledge about organic electronic samples that could be used to improve the function of devices or to solve outstanding problems that could not be solved using other measurement techniques. Since all of the goal involved using INS or other neutron techniques to do something “new” a large part of the effort was focused on complex data interpretation and development of modeling workflows that could be used by other members of the organic electronics community to study their materials. To facilitate the goal of making INS measurements more accessible to a broader community and for a wider range of materials, we optimized and engineered the computational workflows to be well documented, publicly available, and user friendly. The individual titles of the three sequential DOE-BES awards are: (1) Engineering doping profiles in organic semiconducting materials, (2) Measurement and modeling of molecular dopant structure and dynamics, and (3) Quantification of Dynamic Disorder in Electronic and Excitonic Organic Materials. For readability, this report is organized to highlight the progress of the research thematically. Many of the ideas started in one grant period and then spilled over into later awards and collaborative research.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Analysis of Solids Obtained from 201O, 201P, 202A, and 202B Contactors in the Salt Waste Processing Facility

In December 2022, Savannah River Mission Completion personnel recovered solids samples from four contractors (EXT-001, EXT-007, EXT-031, and EXT-040) employed at the Salt Waste Processing Facility (SWPF). These solids were submitted to Savannah River National Laboratory for analysis and characterization. This Technical Report outlines the analytical findings and observations associated with these samples and explores the potential source of solid generation within the SWPF. The following conclusions are offered as a result of this work: the light color solids recovered from EXT-007 are consistent with crystallized, soluble salt components (such as sodium nitrate, sodium hydroxide, and sodium carbonate), all of which can be easily re-dissolved in process water and do not represent a threat to SWPF CSSX processing. While the light solids from the top of EXT-040 were not analyzed, it is believed that they share similar characteristics to those recovered from EXT-007; The dark color solids recovered from EXT-001, EXT-031, and the inside of EXT-040 appear to be complex, consisting of several different types of solids: High concentrations of mercury (Hg, 40-75% by mass) partially present as elemental Hg are consistent with the precipitation of Hg 0 /Hg 2+ from disproportionation of Hg + species in pH swing conditions; Moderate concentrations of titanium (Ti) (2-7%, by mass) are also present, often co-located with Hg. This is consistent with co-precipitation, entrainment, and/or amalgamation of Ti and Hg; Moderate concentrations of iron (Fe) (1-3%, by mass) are observed in SWPF solids and are likely attributable to erosion of stainless-steel components within the CSSX process; Small concentrations of tungsten (W) and cobalt (Co) (<1%, by mass) are observed in SWPF solids. The presence of these solids is consistent with the erosion of components coated with W and Co (e.g., Stellite); Small concentrations of aluminum (Al) (<1%, by mass) are observed in the form of gibbsite and potentially sodium aluminosilicate. The presence of these compounds is likely attributable to the pH swing observed between the scrub and extraction cycles within the CSSX process. The following recommendations are made as a result of this work: Testing should be performed to determine the chemical drivers, process, and mechanism of Hg precipitation in CSSX processing. Special care should be taken to include titanium in tests to evaluate the potential for co-precipitation, entrainment, and amalgamation; Hg-mitigation options that would minimize or eliminate risks of Hg precipitation and solids accumulation (e.g., Hg-absorption techniques, ion exchange/adsorbents, flowsheet changes, etc.) should be assessed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1974

Stratigraphy of the Descartes region /Apollo 16/ - Implications for the origin of samples

Analysis of terrain in the Apollo 16 Descartes landing region shows a series of features that form a stratigraphic sequence which dominates the history and petrogenesis at the site. An ancient 150-km diam crater centered on the Apollo 16 site is one of the earliest recognizable major structures. Nectaris ejecta was concentrated in a regional low at the base of the back slope of the Nectaris basin to form the Descartes Mountains. Subsequently, a 60-km diam crater formed in the Descartes Mountains centered about 25 km to the west of the site. This crater dominates the geology and petrogenetic history of the site. Stone and Smoky Mountains represent the degraded terraced crater walls, and the dark matrix breccias and metaclastic rocks derived from North and South Ray craters represent floor fallback breccias from this cratering event. The interpretation is developed that the stratigraphy of the Cayley and Descartes, and thus the historical record of the Apollo 16 region, documents the complex interaction of deposits and morphology of local and regional impact cratering events. Large local 60- to 150-km diam craters have had a dramatic and previously unrecognized effect on the history and petrology of the Apollo 16 site.

Head, J. W.↗

DOE OSTI2020

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗