Search NASA⌕ Search

SEARCH · Search NASA

Results for “specific yield”

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 163 records · Page 9

Calibration, measurement, and characterization of soil moisture dynamics in a central Amazonian tropical forest

Soil moisture plays a key role in hydrological, biogeochemical, and energy budgets of terrestrial ecosystems. Accurate soil moisture measurements in remote ecosystems such as the Amazon are difficult and limited because of logistical constraints. Time domain reflectometry (TDR) sensors are widely used to monitor soil moisture and require calibration to convert the TDR's dielectric permittivity measurement (K a ) to volumetric water content (θ v ). In this study, our objectives were to develop a field-based calibration of TDR sensors in an old-growth upland forest in the central Amazon, to evaluate the performance of the calibration, and then to apply the calibration to determine the dynamics of soil moisture content within a 14.2 -m-deep vertical soil profile. Depth-specific TDR calibration using local soils in a controlled laboratory setting yielded a novel K a –θ v third-degree polynomial calibration. The sensors were later installed to their specific calibration depth in a 14.2-m pit. The widely used K a –θ v relationship (Topp model) underestimated the site-specific θ v by 22–42%, indicating significant error in the model when applied to these well-structured, clay-rich tropical forest soils. The calibrated wet- and dry-season θ v data showed a variety of depth and temporal variations highlighting the importance of soil textural differentiation, root uptake depths, as well as event to seasonal precipitation effects. Data such as these are greatly needed for improving our understanding of ecohydrological processes within tropical forests and for improving models of these systems in the face of changing environmental conditions.

54 ENVIRONMENTAL SCIENCES↗

Superconducting properties of the spin Hall candidate Ta 3 Sb with eightfold degeneracy

We report the synthesis and characterization of phase pure Ta 3 Sb, a material predicted to be topological with eightfold degenerate fermionic states [Bradlyn et al., Science 353, aaf5037 (2016)] and to exhibit a large spin Hall effect [Derunova et al., Sci. Adv. 5, eaav8575 (2019)]. We observe superconductivity in Ta 3 Sb with T c ~0.67 K in both electrical resistivity ρ(T) and specific heat c(T) measurements. Field-dependent measurements yield the superconducting phase diagram with an upper critical field of H c2 (0) ~0.95 T, corresponding to a superconducting coherence length of ξ ≈ 18.6 nm. The gap ratio deduced from specific heat anomaly, 2Δ 0 /k B T c is 3.46, a value close to the Bardeen-Cooper-Schrieffer value of 3.53. From a detailed analysis of both the transport and thermodynamic data within the Ginsburg-Landau (GL) framework, a GL parameter of κ ≈ 90 is obtained, identifying Ta 3 Sb as an extreme type-II superconductor. The observation of superconductivity in an eightfold degenerate fermionic compound with topological surface states and predicted large spin Hall conductance positions Ta 3 Sb as an appealing platform to further explore exotic quantum states in multifold degenerate systems.

36 MATERIALS SCIENCE↗

Pure quark and gluon observables in collinear drop

We construct a class of pure quark and gluon observables by using the collinear drop grooming technique. The construction is based on linear combinations of multiple cumulative distributions of the jet mass in collinear drop, whose specific weights are fully predicted perturbatively. This yields observables which obtain their values purely from quarks (or purely from gluons) in a wide region of phase space. We demonstrate this by showing that these observables are effective in two phase space regions, one dominated by perturbative resummation and one dominated by nonperturbative effects. The nonperturbative effects are included using shape functions which only appear as a common factor in the linear combinations constructed. We test this construction using a numerical analysis with next-to-leading logarithmic resummation and various shape function models, as well as analyzing these observables with PYTHIA and VINCIA. Choices for the collinear drop parameters are optimized for experimental use.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Energetic Basis and Design of Enzyme Function Demonstrated Using GFP, an Excited-State Enzyme

We report the past decades have witnessed an explosion of de novo protein designs with a remarkable range of scaffolds. It remains challenging, however, to design catalytic functions that are competitive with naturally occurring counterparts as well as biomimetic or nonbiological catalysts. Although directed evolution often offers efficient solutions, the fitness landscape remains opaque. Green fluorescent protein (GFP), which has revolutionized biological imaging and assays, is one of the most redesigned proteins. While not an enzyme in the conventional sense, GFPs feature competing excited-state decay pathways with the same steric and electrostatic origins as conventional ground-state catalysts, and they exert exquisite control over multiple reaction outcomes through the same principles. Thus, GFP is an “excited-state enzyme”. Herein we show that rationally designed mutants and hybrids that contain environmental mutations and substituted chromophores provide the basis for a quantitative model and prediction that describes the influence of sterics and electrostatics on excited-state catalysis of GFPs. As both perturbations can selectively bias photoisomerization pathways, GFPs with fluorescence quantum yields (FQYs) and photoswitching characteristics tailored for specific applications could be predicted and then demonstrated. The underlying energetic landscape, readily accessible via spectroscopy for GFPs, offers an important missing link in the design of protein function that is generalizable to catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scaling High-Resolution Soil Organic Matter Composition to Improve Predictions of Potential Soil Respiration Across the Continental United States

Despite the importance of microbial soil organic matter (SOM) respiration in regulating the flux of carbon between soils and the atmosphere, soil carbon cycling models remain primarily based on climate and soil properties, leading to large uncertainty in predictions. To address this knowledge gap, we analyzed high-resolution water-extractable SOM profiles from soil cores collected across the United States by the 1,000 Soils Pilot of the Molecular Observation Network. Our innovation lies in using machine learning to distill thousands of SOM formula into tractable units; and it enables integrating data from molecular measurements into soil respiration models. In surface soils, SOM chemistry provided better estimates of potential soil respiration than soil physicochemistry, and using them combined yielded the best prediction. Overall, we identify specific subsets of organic molecules that may improve predictions of global soil respiration and create a strong basis for developing new representations in process-based models.

54 ENVIRONMENTAL SCIENCES↗

Inverse design of equilibrium cluster fluids applied to a physically informed model

Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. Here, we then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Metabolic rewiring and biomass redistribution enable optimized mixotrophic growth in Chlamydomonas

Aquatic photosynthetic systems account for approximately one-half of all global carbon assimilation and could be a significant source of renewable fuels and feedstocks. However, rapid growth and biomass production in algae have not always translated into high product yields, partly because central metabolism is context specific, with metabolic fluxes being influenced by nutrient conditions and other environmental factors. In the green microalga Chlamydomonas reinhardtii (Chlamydomonas), mixotrophic cultures (acetate + light) grow far faster than phototrophic (light only) or heterotrophic (acetate + dark) cultures, even though acetate partially suppresses photosynthesis. Here, an isotopic dilution strategy with unlabeled acetate was combined with 13 CO 2 transient labeling to perform isotopically nonstationary metabolic flux analysis (INST-MFA) and to directly compare autotrophic and mixotrophic metabolism in Chlamydomonas supported by data from transcriptomics, proteomics, and metabolomics. INST-MFA indicated that acetate induces a synergistic rewiring of metabolism, conserving carbon by using the glyoxylate cycle and suppressing gluconeogenesis, the latter of which was discordant with omics results and prior models. Additionally, our data provide a plausible rationale for the well-known suppression of photosynthesis by acetate. We propose that reduced total protein content in mixotrophic versus phototrophic cells, much of which is attributed to reduced levels of photosynthetic proteins, decreases the costly metabolic burden of protein synthesis and represents a growth rate optimization strategy.

59 BASIC BIOLOGICAL SCIENCES↗

Cation Dynamics in Hybrid Halide Perovskites

Hybrid halide perovskite semiconductors exhibit complex, dynamical disorder while also harboring properties ideal for optoelectronic applications that include photovoltaics. However, these materials are structurally and compositionally distinct from traditional compound semiconductors composed of tetrahedrally coordinated elements with an average valence electron count of silicon. The additional dynamic degrees of freedom of hybrid halide perovskites underlie many of their potentially transformative physical properties. Neutron scattering and spectroscopy studies of the atomic dynamics of these materials have yielded significant insights into their functional properties. Specifically, inelastic neutron scattering has been used to elucidate the phonon band structure, and quasi-elastic neutron scattering has revealed the nature of the uncorrelated dynamics pertaining to molecular reorientations. Understanding the dynamics of these complex semiconductors has elucidated the temperature-dependent phase stability and origins of defect-tolerant electronic transport from the highly polarizable dielectric response. Furthermore, the dynamic degrees of freedom of the hybrid perovskites provide additional opportunities for application engineering and innovation.

Materials Science↗

Variation in morpho-physiological and metabolic responses to low nitrogen stress across the sorghum association panel

Abstract Background Access to biologically available nitrogen is a key constraint on plant growth in both natural and agricultural settings. Variation in tolerance to nitrogen deficit stress and productivity in nitrogen limited conditions exists both within and between plant species. However, our understanding of changes in different phenotypes under long term low nitrogen stress and their impact on important agronomic traits, such as yield, is still limited. Results Here we quantified variation in the metabolic, physiological, and morphological responses of a sorghum association panel assembled to represent global genetic diversity to long term, nitrogen deficit stress and the relationship of these responses to grain yield under both conditions. Grain yield exhibits substantial genotype by environment interaction while many other morphological and physiological traits exhibited consistent responses to nitrogen stress across the population. Large scale nontargeted metabolic profiling for a subset of lines in both conditions identified a range of metabolic responses to long term nitrogen deficit stress. Several metabolites were associated with yield under high and low nitrogen conditions. Conclusion Our results highlight that grain yield in sorghum, unlike many morpho-physiological traits, exhibits substantial variability of genotype specific responses to long term low severity nitrogen deficit stress. Metabolic response to long term nitrogen stress shown higher proportion of variability explained by genotype specific responses than did morpho-pysiological traits and several metabolites were correlated with yield. This suggest, that it might be possible to build predictive models using metabolite abundance to estimate which sorghum genotypes will exhibit greater or lesser decreases in yield in response to nitrogen deficit, however further research needs to be done to evaluate such model.

59 BASIC BIOLOGICAL SCIENCES↗

Magnetic Microcalorimeter (MMC) Gamma-Ray Detectors with Ultra-High Energy Resolution. Comprehensive Technology Readiness Assessment

Magnetic Microcalorimeter (MMC) gamma ray detectors have demonstrated very high energy resolution of 38 eV FWHM for gamma ray energies below 100 keV. This is an order or magnitude better than conventional HPGe detector and can remove most line-overlap problems in NDA by gamma spectroscopy below ~130 keV. Unlike superconducting transition edge sensors (TESs), MMCs have a reproducible and almost perfectly linear response, which greatly simplifies adding spectra from different detector pixels and reduces systematic errors in quantitative assays. The instrument is currently at TRL-5 but could advance to TRL-6 by improving the signal processing software and using the system for the quantitative analysis of a safeguards sample provided by e.g. the IAEA. Making MMC detectors commercially viable and used more widely then depends on making the technology more reliable, specifically by increasing the device uniformity and yield, and by developing commercial software. This situation is often only half-jokingly referred to as the “valley of death” and could be addressed though the SBIR program. Our collaborators at STAR Cryoelectronics LLC are well qualified to further advance MMC development, and XIA LLC is well-positioned to develop the data acquisition system and write the software for MMC readout. In addition, the advance to TRL-6 would be desirable to increase the visibility and showcase the value of MMC spectrometers, e.g. by placing one of them into the Safeguards Analytical Laboratory (SAL) at the IAEA, whose NDA experts have recently started to show interest in cryogenic detector technology. As a result of this LCP, we also now have an instrument at LLNL that we will use for various applications of MMCs in nuclear safeguards and science. We are confident that as the number of high-profile results increases, so will the interest in and commercial viability of this technology.

42 ENGINEERING↗

Magnetic Microcalorimeter (MMC) Gamma-Ray Detectors with Ultra-High Energy Resolution. Comprehensive Technology Readiness Assessment

Magnetic Microcalorimeter (MMC) gamma ray detectors have demonstrated very high energy resolution of 38 eV FWHM for gamma ray energies below 100 keV. This is an order or magnitude better than conventional HPGe detector and can remove most line-overlap problems in NDA by gamma spectroscopy below ~130 keV. Unlike superconducting transition edge sensors (TESs), MMCs have a reproducible and almost perfectly linear response, which greatly simplifies adding spectra from different detector pixels and reduces systematic errors in quantitative assays. The instrument is currently at TRL-5 but could advance to TRL-6 by improving the signal processing software and using the system for the quantitative analysis of a safeguards sample provided by e.g. the IAEA. Making MMC detectors commercially viable and used more widely then depends on making the technology more reliable, specifically by increasing the device uniformity and yield, and by developing commercial software. This situation is often only half-jokingly referred to as the “valley of death” and could be addressed though the SBIR program. Our collaborators at STAR Cryoelectronics LLC are well qualified to further advance MMC development, and XIA LLC is well-positioned to develop the data acquisition system and write the software for MMC readout. In addition, the advance to TRL-6 would be desirable to increase the visibility and showcase the value of MMC spectrometers, e.g. by placing one of them into the Safeguards Analytical Laboratory (SAL) at the IAEA, whose NDA experts have recently started to show interest in cryogenic detector technology. As a result of this LCP, we also now have an instrument at LLNL that we will use for various applications of MMCs in nuclear safeguards and science. We are confident that as the number of high-profile results increases, so will the interest in and commercial viability of this technology.

42 ENGINEERING↗

Experimental Validation and Continuous Testing of an On-Purpose High-Yield Pitch Synthesis Process for Producing Carbon Fiber from US Domestic Coal

This project aims to develop technology that converts domestic United States (US) raw coal to high quality, high value, and marketable carbon fiber. More specifically, the project aims to significantly improve the selectivity and yield of carbon fiber produced per ton of coal over conventional coal pitch-based production by using low-severity direct coal conversion technology to maximize the yield of pitch from coal, suitable for production of carbon fiber. To meet the stated objective of developing a technology platform, capable of producing high quality, high-value and marketable carbon fiber from domestic US coal, the proposed scope of work involves testing of a low-severity direct coal liquefaction (LS-DCL) process approach. The isotropic pitch produced from the LS-DCL process will be thermally treated using conventional processes to convert to a liquid crystal or “mesophase” pitch as a precursor for structural carbon fibers. This mesophase pitch will then be melt-spun into fiber, oxygen stabilized, and carbonized using conventional processes. The resulting carbon fiber will be evaluated for mechanical properties and suitability for structural applications such as automotive parts and spars for wind turbine blades. The overall process is expected to significantly lower the cost of fiber compared to state of the art fibers produced from polyacrylonitrile (PAN) precursors.

01 COAL, LIGNITE, AND PEAT↗

Leveraging large language models to address data scarcity in machine learning for graphene synthesis

Machine learning in experimental materials science faces significant challenges due to the scarcity of data, which are costly and time-consuming to generate, particularly when relying on in-house experiments. Literature data mining offers a potential solution but introduces issues like mixed data quality, inconsistent formats, and non-uniform reporting of synthesis parameters, resulting in partially missing and heterogeneous features across the dataset. Here, we propose data imputation and feature engineering methods that employ pre-trained large language models (LLMs) to enhance machine learning performance on scarce, heterogeneous datasets, demonstrated on graphene CVD synthesis data and the ML-HydPARK hydrogen storage dataset. GPT models perform data imputation via tailored prompting and semantic normalization of inconsistently reported features through embeddings, for example, to harmonize the complex nomenclature of CVD substrates. Beyond yielding more diverse and richer feature representations than traditional methods such as K-nearest neighbors (KNN) and Multivariate Imputation by Chained Equations (MICE), LLM-based data imputation is evaluated against dataset characteristics and prompting strategies. We vary the level of autonomy granted to the LLM, from generic prompting that leverages pre-trained knowledge for autonomous data generation to data-informed prompting that constrains outputs using target-specific information, and demonstrate which level of autonomy yields superior imputation performance across datasets and feature types. The proposed data engineering methods markedly improve downstream performance; for example, in graphene layer number classification using a support vector machine (SVM), binary accuracy increases from 39% to 65% and ternary accuracy from 52% to 72%. Fine-tuning experiments on both datasets show that combining our proposed LLM-based data imputation and feature encoding methods with numerical machine learning predictors outperforms standalone fine-tuned LLM predictors in data-scarce settings. The proposed strategies emphasize data enhancement techniques rather than refining learning architectures or regularizing loss functions, offering a broadly applicable framework for improving machine learning performance on scarce, inhomogeneous datasets.

Chemical vapor deposition↗

Coproduction of Light and Heavy r-process Elements via Fission Deposition

We apply for the first time fission yields determined across the chart of nuclides from the macroscopic-microscopic theory of the Finite Range Liquid Drop Model to simulations of rapid neutron capture (r-process) nucleosynthesis. With the fission rates and yields derived within the same theoretical framework utilized for other relevant nuclear data, our results represent an important step toward self-consistent applications of macroscopic-microscopic models in r-process calculations. The yields from this model are wide for nuclei with extreme neutron excess. We show that these wide distributions of neutron-rich nuclei, and particularly the asymmetric yields for key species that fission at late times in the r process, can contribute significantly to the abundances of the lighter heavy elements, specifically the light precious metals palladium and silver. Since these asymmetric yields correspondingly also deposit into the lanthanide region, we consider the possible evidence for coproduction by comparing our nucleosynthesis results directly with the trends in the elemental ratios of metal-poor stars rich in r-process material. We show that for r-process enhanced stars palladium over europium and silver over europium display mostly flat trends suggestive of coproduction and compare to the lanthanum over europium trend which is often used to justify robustness arguments in the lanthanide region. We find that such robustness arguments may be extendable down to palladium and heavier and demonstrate that fission deposition is a mechanism by which such a universality or robustness can be achieved.

79 ASTRONOMY AND ASTROPHYSICS↗

Varied farm-level carbon intensities of corn feedstock help reduce corn ethanol greenhouse gas emissions

Abstract A reduction in the overall carbon intensity (CI) of a crop-based biofuel can be achieved by cutting down the CI of the biofuel’s feedstock, which in turn correlates significantly to agricultural management practices. Proposals are being made to incentivize low-carbon biofuel feedstocks under U.S. fuel regulatory programs to promote sustainable farming practices by individual farms. For such an incentive scheme to function properly, robust data collection and verification are needed at the farm level. This study presents our collaboration with U.S. private sector companies to collect and verify the corn production data necessary for feedstock-specific CI calculation at the farm level, through a carefully designed questionnaire, to demonstrate the practicality and feasibility of data collection at scale. We surveyed 71 farms that produced 0.2 million metric tons of corn grain in 2018 in a Midwestern U.S. state to obtain information on key parameters affecting corn ethanol feedstock CI, such as grain yields, fertilizer/chemical application rates, and agronomic practices. Feedstock-specific CI was calculated in the unit of grams (g) CO 2 equivalent (CO 2 e) of greenhouse gases per kilogram (kg) of corn produced. Results showed large CI variations—from 119 to 407 g CO 2 e kg −1 of corn—due to the farm-level inventory, while the production-weighted average CI for all surveyed farms was 210 g CO 2 e kg −1 , comparable to the national average CI of 204 g CO 2 e kg −1 . The nitrogen fertilizer type applied and rate were identified as key factors contributing most to CI variations at the farm level. The estimated N 2 O emissions from fertilizer and biomass nitrogen inputs to soil accounted for 51% of the overall farm-level CI and therefore need to be better monitored at farm level with high resolution. We concluded that this feedstock-specific, farm-level CI evaluation has the potential to be used to incentivize low-carbon feedstock for biofuel production.

54 ENVIRONMENTAL SCIENCES↗

Enzymatic depolymerization of polyester: Foaming as a pretreatment to increase specific surface area

Abstract Poly(ethylene terephthalate) (PET) is widely used for its high strength‐to‐weight ratio, gas barrier properties, and chemical resistance. The growing PET use highlights the demand for a better recycling system. Enzymatic recycling, alongside mechanical and chemical methods, is eco‐friendly and yields properties similar to virgin PET. Substrate properties ( T g , crystallinity, and specific surface area [SSA]) and enzyme stability significantly impact conversion efficiency. Higher SSA and lower crystallinity tend to yield improved depolymerization when employing leaf compost‐cutinase (LCC‐ICCG) enzymes. This study explored melt extrusion and foaming as pretreatment techniques to modify PET structural properties, using a low‐cost chemical foaming agent (CFA). The monomer conversion rate and efficiency during depolymerization were measured and related to the processing, extrudate micro‐ and meso‐structure, and polyester type. Pretreated PET substrates showed reduced T g , crystallinity, density, and enhanced SSA, resulting in a 90% mass loss for foamed RPET and VPET substrates within 2 days. In contrast, PET with ~30% of cyclohexanedimethanol comonomer exhibited a nearly 50% lower depolymerization rate, with zero BHET production. It indicates that the combination of low crystallinity, low T g , and high SSA leads to improved monomer conversion. These findings emphasize the significance of amorphization and foaming in enhancing PET enzymatic depolymerization.

42 ENGINEERING↗

Geometric prior of multi-resolution yielding manifolds and the local closest point projection for nearly non-smooth plasticity

Elastoplasticity models often introduce a scalar-valued yield function to implicitly represent the boundary between elastic and plastic material states. Here this paper introduces a new alternative where the yield envelope is represented by a manifold of which the topology and the geometry are learned from a set of data points in a parametric space (e.g. principal stress space, π-plane). Here, deep geometric learning enables us to reconstruct a highly complex yield envelope by breaking it down into multiple coordinate charts. The global atlas that consists of these coordinate charts in return allows us to represent the yield surface via multiple overlapping patches, each with a specific local parametrization. This setup provides several advantages over the classical implicit function representation approach. For instance, the availability of coordinate charts enables us to introduce an alternative stress integration algorithm where the trial stress may project directly on a local patch and hence circumvent the issues related to non-smoothness and the lack of convexity of yield surfaces. Meanwhile, the local parametric approach also enables us to predict hardening/softening locally in the parametric space, even without complete knowledge of the yield surface. Comparisons between the classical yield function approach on the non-smooth plasticity and anisotropic cam-clay plasticity model are provided to demonstrate the capacity of the models for highly precise yield surface and the feasibility of the implementation of the learned model in the local stress integration algorithm.

42 ENGINEERING↗

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗