Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional performance testing”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Studying excited-state-specific perturbation theory on the Thiel set

Here, we explore the performance of a recently introduced N 5 -scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set’s inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles norm offers a relatively low-cost way to test for doubly excited character without the need to define an active space.

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

Deep Reinforcement Learning Based Volt-VAR Optimization in Smart Distribution Systems

This paper develops a model-free volt-VAR optimization (VVO) algorithm via multi-agent deep reinforcement learning (DRL) in unbalanced distribution systems. This method is novel since we cast the VVO problem in distribution networks to an intelligent deep Q-network (DQN) framework, which avoids solving a specific optimization model directly when facing time-varying operating conditions in the systems. We consider statuses/ratios of switchable capacitors, voltage regulators, and smart inverters installed at distributed generators as the action variables of the agents. A delicately designed reward function guides these agents to interact with the distribution system, in the direction of reinforcing voltage regulation and power loss reduction simultaneously. The forward-backward sweep method for radial three-phase distribution systems provides accurate power flow results within a few iterations to the DRL environment. The proposed method realizes the dual goals for VVO. We test this algorithm on the unbalanced IEEE 13-bus and 123-bus systems. Numerical simulations validate the excellent performance of this method in voltage regulation and power loss reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2020

Pressurized Gas Beam Monitor for Extremely High Intensities – Phase II

A novel pressurized gas-filled multi-RF-cavity beam profile monitor has been studied that is simple and robust in high-radiation environments. Charged particles passing through each RF-cavity in the monitor produce intensity-dependent ionized plasma, which changes the gas permittivity. Standard RF techniques to measure the change in quality factor (Q) and frequency (f) as a function of time are then used to determine the change in permittivity and corresponding beam intensity in each cavity in the profile monitor. The sensitivity to beam intensity is adjustable using gas pressure and RF gradient. The performance of the gas-filled beam profile monitor has been numerically simulated to evaluate the sensitivity of permittivity measurements and the effectiveness of calibration strategies. An entire multi-RF-cavity system has been designed and a demonstration test has been prepared.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Digital Twin Development of FASTR

This work package is focused on creating a digital twin (DT) of the Facility to Alleviate Salt Technology Risks (FASTR) experiment. FASTR is a high-profile facility tasked with maturing molten salt applications such as liquid salt energy storage, concentrated solar power facilities, and molten salt cooling systems supporting other advanced reactor technologies, all of which may be incorporated as components within the Integrated Energy Systems (IES) program. As part of creating a digital twin (DT), verification and validation (V&V) has been performed for the physics-based model, and gaps have been identified. Additionally, a preliminary integration of the FASTR model has been completed with the existing physics-based model of the Thermal Energy Distribution System (TEDS)/Microreactor AGile Non-nuclear Experimental Testbed (MAGNET) facility at Idaho National Laboratory (INL). A Functional Mock-up Unit (FMU) of the physics-based system model and a reduced-order model (ROM) were both used to test for a preliminary hardware in the loop operation when the physical loop was offline. Future tasks could include improving the physics-based model to better match the available experimental data, improving the hardware in the loop integration, and improving the integration with the TEDS model.

14 SOLAR ENERGY↗
DOE OSTI2024

Celeritas R&D Report: Accelerating Geant4

Celeritas is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a Graphics Processing Unit (GPU)-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.4, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling straightforward integration of Celeritas into the high energy physics (HEP) experimental frameworks CMSSW and ATLAS FullSimLight. On the Perlmutter supercomputer, Celeritas performs EM physics between 3× and 18× faster using the machine’s Nvidia GPUs compared to using only CPUs, corresponding to an electrical power efficiency up to a factor of 5. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.7–2.2× on GPU and 1.2× on CPU. In a CMS test application using tt¯ events and the prototype Run 4 configuration, compared to Geant4 CPU, Celeritas with a Nvidia A100 improves overall throughput up to a factor of 2.7× but cannot be efficiently shared with more than 8 cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Model-Free Primal-Dual Methods for Network Optimization with Application to Real-Time Optimal Power Flow

This paper examines the problem of real-time optimization of networked systems and develops online algorithms that steer the system towards the optimal trajectory without explicit knowledge of the system model. The problem is modeled as a dynamic optimization problem with time-varying performance objectives and engineering constraints. The design of the algorithms leverages the online zero-order primal-dual projected-gradient method. In particular, the primal step that involves the gradient of the objective function (and hence requires a networked systems model) is replaced by its zero-order approximation with two function evaluations using a deterministic perturbation signal. The evaluations are performed using the measurements of the system output, hence giving rise to a feedback interconnection, with the optimization algorithm serving as a feedback controller. The paper provides some insights on the stability and tracking properties of this interconnection. Finally, the paper applies this methodology to a real-time optimal power flow problem in power systems, and shows its efficacy on the IEEE 37-node distribution test feeder for reference power tracking and voltage regulation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Self-Cleaning Highly Porous TiO 2 Coating Designed by Swelling-Assisted Sequential Infiltration Synthesis (SIS) of a Block Copolymer Template

Photocatalytic self-cleaning coatings with a high surface area are important for a wide range of applications, including optical coatings, solar panels, mirrors, etc. Here, we designed a highly porous TiO 2 coating with photoinduced self-cleaning characteristics and very high hydrophilicity. This was achieved using the swelling-assisted sequential infiltration synthesis (SIS) of a block copolymer (BCP) template, which was followed by polymer removal via oxidative thermal annealing. The quartz crystal microbalance (QCM) was employed to optimize the infiltration process by estimating the mass of material infiltrated into the polymer template as a function of the number of SIS cycles. This adopted swelling-assisted SIS approach resulted in a smooth uniform TiO 2 film with an interconnected network of pores. The synthesized film exhibited good crystallinity in the anatase phase. The resulting nanoporous TiO 2 coatings were tested for their functional characteristics. Exposure to UV irradiation for 1 h induced an improvement in the hydrophilicity of coatings with wetting angle reducing to unmeasurable values upon contact with water droplets. Furthermore, their self-cleaning characteristics were tested by measuring the photocatalytic degradation of methylene blue (MB). The synthesized porous TiO 2 nanostructures displayed promising photocatalytic activity, demonstrating the degradation of approximately 92% of MB after 180 min under ultraviolet (UV) light irradiation. Thus, the level of performance was comparable to the photoactivity of commercial anatase TiO 2 nanoparticles of the same quantity. Our results highlight a new robust approach for designing hydrophilic self-cleaning coatings with controlled porosity and composition.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Model-Free Primal-Dual Methods for Network Optimization with Application to Real-Time Optimal Power Flow: Preprint

This paper examines the problem of real-time optimization of networked systems and develops online algorithms that steer the system towards the optimal trajectory without explicit knowledge of the system model. The problem is modeled as a dynamic optimization problem with time-varying performance objectives and engineering constraints. The design of the algorithms leverages the online zero-order primal-dual projected-gradient method. In particular, the primal step that involves the gradient of the objective function (and hence requires networked systems model) is replaced by its zero-order approximation with two function evaluations using a deterministic perturbation signal. The evaluations are performed using the measurements of the system output, hence giving rise to a feedback interconnection, with the optimization algorithm serving as a feedback controller. The paper provides some insights on the stability and tracking properties of this interconnection. Finally, the paper applies this methodology to a real-time optimal power flow problem in power systems, and shows its efficacy on the IEEE 37-node distribution test feeder for reference power tracking and voltage regulation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Loss Landscape Analysis for Reliable Quantized ML Models for Scientific Sensing

In this paper, we propose a method to perform empirical analysis of the loss landscape of machine learning (ML) models. The method is applied to two ML models for scientific sensing, which necessitates quantization to be deployed and are subject to noise and perturbations due to experimental conditions. Our method allows assessing the robustness of ML models to such effects as a function of quantization precision and under different regularization techniques -- two crucial concerns that remained underexplored so far. By investigating the interplay between performance, efficiency, and robustness by means of loss landscape analysis, we both established a strong correlation between gently-shaped landscapes and robustness to input and weight perturbations and observed other intriguing and non-obvious phenomena. Our method allows a systematic exploration of such trade-offs a priori, i.e., without training and testing multiple models, leading to more efficient development workflows. This work also highlights the importance of incorporating robustness into the Pareto optimization of ML models, enabling more reliable and adaptive scientific sensing systems.

Baldi, Tommaso [Pisa, Scuola Normale Superiore]↗
DOE OSTI2020

Mechanoadaptive strain and functional osseointegration of dental implants in rats

Spatiotemporal implant-bone biomechanics and mechanoadaptive strains in peri-implant tissue are poorly understood. Physical and chemical characteristics of an implant-bone complex (IBC) were correlated in three-dimensional space (along the length and around a dental implant) to gather insights into time related integration of the implant with the cortical portion of a jaw bone in a rat. Rats (N = 9) were divided into three experimental groups with three rats per time point; 3-, 11-, and 24-day. All rats were fed crumbled hard pellets mixed with water (soft-food diet) for the first 3 days followed by a hard-food diet with intact hard-food pellets (groups of 11- and 24-day only). Biomechanics of the IBCs harvested from rats at each time point was evaluated by performing mechanical testing in situ in tandem with X-ray imaging. The effect of physical association (contact area) of a loaded implant with adapting peri-implant tissue, and resulting strain within was mapped by using digital volume correlation (DVC) technique. The IBC stiffness at respective time points was correlated with mechanical strain in peri-implant tissue. Results illustrated that IBC stiffness at 11-day was lower than that observed at 3-day. However, at 24-day, IBC stiffness recovered to that which was observed at 3-day. Furthermore, correlative microscopy and spectroscopy illustrated that the lower IBC stiffness was constituted by softer and less mineralized peri-implant tissue that contained varying expressions of osteoconductive elements. Lower IBC stiffness observed at 11-day was constituted by less mineralized peri-implant tissue with osteoconductive elements that included phosphorus (P) which was co-localized with higher expression of zinc (Zn), and lower expression of calcium (Ca). Higher IBC stiffness at 24-day was constituted by mineralized peri-implant tissue with higher expressions of osteoconductive elements including Ca and P, and lower expressions of Zn. These spatiotemporal correlative maps of peri-implant tissue architecture, heterogeneous distribution of mineral density, and elemental colocalization underscore mechanoadaptive physicochemical properties of peri-implant tissue that facilitate functional osseointegration of an implant. These results provided insights into 1) plausible "prescription" of mechanical loads as an osteoinductive "therapeutic dose" to encourage osteoconductive elements in the peri-implant tissue that would facilitate functional osseointegration of the implant; 2) a "critical temporal window" between 3 and 11 days, and perhaps it is this acute phase during which key candidate regenerative molecules can be harnessed to accelerate osseointegration of an implant under load.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

A Low-Voltage Layered Na 2 TiGeO 5 Anode for Lithium-Ion Battery

Titanium-based anode materials have achieved much progress with the wide studies in lithium-ion batteries. However, these known materials usually possess high discharge voltage platforms and limited energy densities. In this work, a titanium-based oxide of Na 2 TiGeO 5 with layered structure, two-dimensional lamellar frame and exposed highly active (001) facet, exhibiting good electrochemical performance in terms of high capacity (410 mAh g –1 with a current density of 50 mA g –1 ), excellent rate capability and cycling stability with no obvious capacity attenuation after 4000 cycles, is reported. The appropriate discharge voltage plateau at around 0.2 V endows the Na 2 TiGeO 5 anode material high security compared with graphite and high energy density compared with spinel Li 4 Ti 5 O 12 . Combining the electrochemical tests and the density functional theory calculations, the Li + storage mechanism of Na 2 TiGeO 5 is elucidated and the conversion reaction process is revealed. More importantly, this study provides a way to develop low-voltage and high-capacity titanium-based anode materials for efficient energy storage.

25 ENERGY STORAGE↗
DOE OSTI2021

Light Propagation Considerations For Internally Clad Sapphire Optical Fiber Using the 6Li(n,)3H Reaction

Optical frequency domain reflectometry measurements in internally clad single crystal sapphire fiber have received attention in recent years due to their high temperature distributed sensing potential. As work with these fibers has proceeded, there have been some inconsistencies in the results. Deeper investigation and testing has identified two critical considerations for the proper functioning of these fibers. First, users must address the propagation of multimode light along the outer surface of the fiber. Here, by observing the far-field image of an internally clad sapphire fiber when adding index matching fluid to the outer fiber surface, we demonstrate the effects of removing the higher order modes from these fibers. The addition of index matching fluid resulted in nearly single mode performance where multimode performance was previously observed. Second, users must address the effect that coupling the fiber to the interrogator via silica based fiber has on the internally clad sapphire fibers performance. Direct fusion splicing of silica to sapphire, as has been used in the recent work with these fibers, has a mode filtering effect which can be beneficial towards the modal behavior of the fibers. However, in this paper we demonstrate that the splicing can cause a sensing failure due to little or no low order mode light, that is useful for sensing, returning to the detector. The positive results from recent years have demonstrated that optical frequency domain reflectometry sensing performance will be successful in clad sapphire fiber; but only when the considerations described herein are addressed properly.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Evaluating Space Object Conjunction Probabilities Using Characteristic Function Inversion

This report discusses an approach to computing the probability of a conjunction between two space objects in the short-term encounter scenario. A conjunction is defined here as an event where the miss distance between the objects is less than some specified value. The scenario assumptions are that the motion of the objects is linear, their positions are Gaussian distributed, and their velocities are known and constant. Under these assumptions, the squared-miss distance is shown to have the generalized chi-square distribution. An established statistical technique called characteristic function inversion is employed to evaluate the distribution and obtain conjunction probabilities. The method is closely related to a recent approach based on moment generating function inversion, and a qualitative comparison of the approaches is provided. Last, the method is tested on several benchmark test cases where it agrees with numerical integration on the cases with conjunction probabilities above 10 –12 . However, the exact probability in these cases is usually not needed and this probability can be bounded above using an independent Gaussian approximation. Overall, the report shows how to compute conjunction probabilities using a standard statistical method, though numerical integration seems to perform equally well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Novel niobium-doped titanium oxide towards electrochemical destruction of forever chemicals

Electrochemical advanced oxidative processes (EAOP) are a promising route to destroy recalcitrant organic contaminants such as per- and polyfluoroalkyl substances (PFAS) in drinking water. Central to EAOP are catalysis-induced reactive free radicals for breaking the carbon fluorine bonds in PFAS. Generating these reactive species electrochemically at electrodes provides an advantage over other oxidation processes that rely on chemicals or other harsh conditions. Herein, we report on the performance of niobium (Nb) doped rutile titanium oxide (TiO 2 ) as a novel EAOP catalytic material, combining theoretical modeling with experimental synthesis and characterization. Calculations based on density functional theory are used to predict the overpotential for oxygen evolution at these candidate electrodes, which must be high in order to oxidize PFAS. The results indicate a non-monotonic trend in which Nb doping below 6.25 at.% is expected to reduce performance relative to TiO 2 , while higher concentrations up to 12.5 at.% lead to increased performance, approaching that of state-of-the-art Magnéli Ti 4 O 7 . TiO 2 samples were synthesized with Nb doping concentration at 10 at.%, heat treated at temperatures from 800 to 1100 °C, and found to exhibit high oxidative stability and high generation of reactive oxygen free radical species. The capability of Nb-doped TiO 2 to destroy two common species of PFAS in challenge water was tested, and moderate reduction by ~ 30% was observed, comparable to that of Ti 4 O 7 using a simple three-electrode configuration. We conclude that Nb-doped TiO 2 is a promising alternative EAOP catalytic material with increased activity towards generating reactive oxygen species and warrants further development for electrochemically destroying PFAS contaminants.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Experimental study of compaction localization in carbonate rock and constitutive modeling of mechanical anisotropy

Abstract Sedimentary rocks are inherently anisotropic and prone to strain localization. While the influence of rock anisotropy on the brittle/dilative regime has been studied extensively, its influence on the ductile/compactive regime is much less explored. This paper discusses the anisotropic behavior of a high‐porosity carbonate rock from central Europe (the Maastricht Tuffeau). A set of triaxial tests with concurrent x‐ray tomography has been performed at different confining pressures. The anisotropic characteristics of this rock have been investigated by testing samples cored at different inclinations of the bedding, thus revealing non‐negligible effects of the coring direction on yielding and compaction behavior. Specifically, samples cored perpendicular to bedding display higher strength and longer stages of post‐yielding deformation before manifesting re‐hardening. Despite such alterations of the inelastic response, Digital Image Correlation has revealed that the strain localization mode is independent of the coring direction, thus being primarily affected by the confinement level. To capture the observed interaction between material anisotropy and compaction behavior at the continuum‐scale, an elastoplastic constitutive law has been proposed. For this purpose, a set of tensorial bases has been introduced to replicate how the oriented rock fabric modulates the yielding and plastic flow characteristics of the material. The analyses show that the impact of the coring direction on yield function and plastic flow rule is fundamentally different, thus requiring the use of distinct projection strategies (a strategy here defined heterotopic mapping). The performance of the model, studied through parametric analyses and by calibrating the experimental results, illustrates the improved capability of the proposed constitutive approach when applied to strongly anisotropic porous rocks.

Shahin, Ghassan↗
DOE OSTI2024

Strategically Modified Ligand Incorporating Mixed Phosphonate and Carboxylate Groups to Enhance Performance in All-Iron Redox Flow Batteries

Iron redox flow batteries (Fe-RFBs) hold significant promise for achieving cost-effectiveness and utilizing abundant materials for stationary energy storage applications. Here, a design of a novel Fe complex utilizing a nitrogenous phosphonate/carboxylate mixed ligand, N,N-Bis(phosphonomethyl)glycine (BMPG), is presented to achieve high performance Fe anolyte. Compared to its all-phosphonate form, nitrilotri(methylphosphonic acid) (NTMPA), the new complex Fe(BPMG) 2 demonstrates a negatively shifted redox potential, resulting in ≈0.07 V (≈10%) increase in battery output voltage. Full battery testing paired with ferrocyanide catholyte demonstrates stable cycling (capacity degradation <0.0001%/cycle) over 730 consecutive charge/discharge cycles with Coulombic Efficiency of 100% at a current density of 20 mA cm -2 under near neutral pH (≈8). Of particular interest, density functional theory (DFT) studies and operando Raman measurements provide strong evidence supporting a molecular structure in BPMG, which reveals the mixed phosphonate/carboxylate groups in BPMG maintain the octahedral coordination of the Fe ion center with phosphonates exclusively, while leaving the carboxylate unbound for both Fe(II) and Fe(III) complexes. This structural similarity between BPMG-based Fe(II) and Fe(III) complexes effectively mitigates the slow redox reaction kinetics observed in Fe(NTMPA) 2 anolyte, where significant ligand reorientation occurs between Fe(II) and Fe(III) complexes.

25 ENERGY STORAGE↗
DOE OSTI2019

A visco-plastic framework for interface processes in sedimentary reservoir rocks at HPHT conditions

As energy operations face the challenge of reservoirs at ever-increasing depths, modelling the response of reservoir rocks at high pressure and high temperature (HPHT) conditions is a crucial step for successfully unlocking new resources. At these conditions, rocks can experience thermal and pressure sensitivity, as well as admit internal phase changes at the pore-solid interfaces that determine their macroscopic response to external loading. It is therefore expected that constitutive laws of Geomechanics should be enriched to accommodate such effects. In this work, a multi-physics constitutive theory for sedimentary rocks is proposed within the general framework of viscoplasticity. The viscosity of the material is assumed to be a function of the temperature, pore-pressure and energy required to alter the inter-granular interfaces. This energy is expressed through the chemical potentials of the phases involved during processes like chemical dissolution/precipitation of the interfaces or mechanical debonding. The resulting flow law and corresponding stress equilibrium are coupled to the energy and mass conservation laws, constituting a closed system of equations, which is solved using the Finite Element simulator REDBACK. A series of numerical tests for performance and calibration is then successfully performed against different types of reservoir rocks, confining pressures and temperatures (from room temperature to over 800K). Here, we show that the mechanical response of sedimentary porous rocks at strains usually achieved in laboratory testing can be explained by accounting for the interface processes taking place at the cementitious material bonding the grains, a process that can also determine the brittle-to-ductile transition of sedimentary rocks.

58 GEOSCIENCES↗