Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model predictive control scheme”

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 91 records · Page 5

Resonance Suppression based on Predictive Control of Grid-following Inverters with LCL Filter in Weak Grid Condition

This paper presents an effective resonance suppression control when grid interactive inverter with LCL filter experience weak grid conditions. The resonance suppression mechanism is based on finite-set model predictive control (FSMPC) with adaptive cost function to alter the controller objectives seamlessly. The large parasitic impedance and low short circuit ratio (SCR) of weak grid challenges the operation of gridconnected inverters. Specifically, the LCL filter resonance may get excited, resulting in collapse of the inverter operation. To address this issue, an effective active damping approach is implemented in the proposed control. In the proposed active damping scheme for resonance suppression, the feedback currents are switched based on the operating conditions of inverter. During the stiff grid conditions, grid current serves as feedback to the controller. While in the weak grid conditions, the inverter current is selected as feedback signal. The toggling action between these two feedback currents is determined by comparing moving RMS of grid current with threshold current as constraint in the cost function of the proposed MPC scheme. The threshold current for toggling action is based on the reference active and reactive power setpoints plus an acceptable band of variation. The theoretical analysis is verified by several case studies.

Umar, Muhammad Farooq↗

A Triple-Moment Representation of Ice in the Predicted Particle Properties (P3) Microphysics Scheme

In the original Predicted Particle Properties (P3) bulk microphysics scheme, all ice-phase hydrometeors are represented by one or more “free” ice categories, where the physical properties evolve smoothly through changes to the four prognostic variables (per category), and with a two-moment representation of the particle size distribution. As such, the spectral dispersion cannot evolve independently, which thus results in limitations in representation of ice—in particular, hail—due to necessary constraints in the scheme to prevent excessive gravitational size sorting. To overcome this, P3 has been modified to include a three-moment representation of the size distribution of each ice category through the addition of a fifth prognostic variable, the sixth moment of the size distribution. The details of the three-moment ice parameterization in P3 are provided. The behavior of the modified scheme, with the single-ice-category configuration, is illustrated through simulations in a simple 1D kinematic model framework as well as with near large-eddy-resolving (250-m grid spacing) 3D simulations of a hail-producing supercell. It is shown that the three-moment ice configuration controls size sorting in a physically based way and leads to an improved capacity to simulate large, heavily rimed ice (hail), including mean and maximum sizes and reflectivity, and thus an overall improvement in the representation of ice-phase particles in the P3 scheme.

54 ENVIRONMENTAL SCIENCES↗

Assessing turbulence and mixing parameterizations in the gray-zone of multiscale simulations over mountainous terrain during the METEX21 field experiment

Multiscale numerical weather prediction models transition from mesoscale, where turbulence is fully parameterized, to microscale, where the majority of highly energetic scales of turbulence are resolved. The turbulence gray-zone is situated between these two regimes and multiscale models must downscale through these resolutions. Here, we compare three multiscale simulations which vary by the parameterization used for turbulence and mixing within the gray-zone. The three parameterizations analyzed are the Mellor-Yamada Nakanishi and Niino (MYNN) Level 2.5 planetary boundary layer scheme, the TKE-1.5 large eddy simulation (LES) closure scheme, and a recently developed three-dimensional planetary boundary layer scheme based on the Mellor-Yamada model. The simulation domain includes complex (i.e., mountainous) terrain in Nevada that was instrumented with meteorological towers, profiling and scanning lidars, a tethered balloon, and a surface flux tower. Simulations are compared to each other and to observations, with assessment of model skill at predicting wind speed, wind direction and TKE, and qualitative evaluations of transport and dispersion of smoke from controlled releases. This analysis demonstrates that microscale predictions of transport and dispersion can be significantly influenced by the choice of turbulence and mixing parameterization in the terra incognita , particularly over regions of complex terrain and with strong local forcing. This influence may not be apparent in the analysis of model skill, and motivates future field campaigns involving controlled tracer releases and corresponding modeling studies of the turbulence gray-zone.

54 ENVIRONMENTAL SCIENCES↗

Tendencies of Soil Microbial NO Emissions During HI‐SCALE as Predicted by a Nitrification/Denitrification Scheme

Many atmospheric chemical processes, including the formation of secondary organic aerosol (SOA), are strongly modulated by the reactions of NO and NO 2 (NO x ). Though NO x is controlled by anthropogenic emissions near urban areas, in rural areas soil microbes can be a significant contribution to NO emission globally. The relative rates of emissions of different nitrogen-containing species (e.g., NO, N 2 O, HONO, and N 2 ) are strong functions of soil properties such as temperature, moisture content, pH, and soil carbon and nitrogen pools. However, typical large-scale biogeochemical models either express these emissions simplistically, or not at all. Here we investigate the potential impact of soil NO emissions on atmospheric chemistry and SOA formation over regional and monthly time scales, specifically the 2016 spring and summer Intensive Observational Periods (IOPs) of the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign. We implement the soil NO nitrification/denitrification parameterization of Rasool et al. (2019) into the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem), supplemented by a 1-km soil moisture analysis. We then simulate both IOPs over the U.S. Great Plains, evaluating against ground stations and flight data. We show that soil NO emissions can account for a large fraction of total NO x and locally increase O 3 concentrations by up to 25%, while alleviating negative biases of gases and aerosols toward observations. Soil moisture and temperature changes between IOP1 and IOP2 lead to overall differences in emissions, but with large regional variability due to heterogeneous surface characteristics.

Atmoshpheric Chemistry↗

The use of idealised experiments in testing a new convective parametrization: Performance of CoMorph-A

CoMorph is a new mass-flux convection parametrization under development at the Met Office designed for use within the Unified Model and its successor model, LFRic. Use of a three-dimensional idealised model enables controlled tests of the performance of the scheme across different regimes. This includes the interaction between the physical parametrizations and the resolved dynamics, allowing study of the emergent organisation of convection on the resolved scale. A selection of well-known cases is revisited here, with the purpose of documenting the extent to which CoMorph captures a range of important, but challenging, behaviour such as the diurnal cycle and sensitivity to tropospheric moisture. Simulations using CoMorph-A, a new physics package, that has been demonstrated to perform well at numerical weather prediction (NWP) and climate scales, are compared against the current global atmosphere configuration and high-resolution results. In addition to an entirely new convection scheme, the package of changes includes significant changes to the cloud, microphysics, and boundary-layer parametrizations. Recognising that CoMorph-A is the first version of a scheme that will continue to be substantially developed and to obtain good performance, compromises in tuning have had to be made. These idealised tests therefore show what works well in this configuration, and what areas will require further work. As such, it is quite a demanding testbed and could be viewed as some of the equipment required for a “convective playground”.

54 ENVIRONMENTAL SCIENCES↗

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗

Selective acoustic control of photon-mediated qubit-qubit interactions

Quantum technologies such as quantum sensing, quantum imaging, quantum communications, and quantum computing rely on the ability to actively manipulate the quantum state of light and matter. Quantum emitters, such as color centers trapped in solids, are a useful platform for the realization of elementary building blocks (qubits) of quantum information systems. In particular, the modular nature of such solid-state devices opens up the possibility to connect them into quantum networks and create nonclassical states of light shared among many qubits. The function of a quantum network relies on efficient and controllable interactions among individual qubits. In this context, we present a scheme where optically active qubits of differing excitation energies are mutually coupled via a dispersive interaction with a shared mode of an optical cavity. This generally off-resonant interaction is prohibitive of direct exchange of information among the qubits. However, we propose a scheme in which by acoustically modulating the qubit excitation energies it is, in fact, possible to tune to resonance a preselected pair of qubits and thus open a communication channel between them. This method potentially enables fast (approximately nanoseconds) and parallelizable on-demand control of a large number of physical qubits. In this work, we develop an analytical and a numerical theoretical model demonstrating this principle and suggest feasible experimental scenarios to test the theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling active control of resistive wall mode with power saturation and sensor noise on HL-2M

The resistive wall mode (RWM) control on the HL-2M tokamak is simulated with the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681), aiming at quantifying control current and voltage requirements when more realistic issues are taken into account, i.e. the control power saturation and the sensor signal noise. The fluid model predicts a narrow stability region for the n = 1 RWM without magnetic feedback, in the 2D parameter space of the plasma pressure versus the toroidal flow speed. Magnetic feedback can fully stabilize the RWM on HL-2M. Without considering the voltage limitation and the sensor signal noise, it is found that plasma flow helps active control of the mode, by reducing the required critical feedback gain for both flux-to-current and flux-to-voltage control schemes. In the absence of the sensor signal noise, the lowest control voltage saturation level, below which the RWM control is lost, is found to roughly satisfy a linear relation to the plasma flow frequency, indicating that subsonic plasma flow is effective in relaxing the control power requirement for the RWM feedback stabilization. The presence of the sensor signal noise substantially modifies the feedback results. A statistical study finds that the sensor signal noise, with the standard deviation of 0.1 G on HL-2M, roughly doubles the required control voltage for successful mode control. The synergistic stabilization effect due to plasma flow is somewhat weakened by the presence of the sensor signal noise. At a given rotation, the tolerable voltage limit generally increases with increasing feedback gain due to the sensor signal noise.

Physics↗

Distributed Load Shedding Application Architecture and Bi-Level Predictive Estimator Algorithm

Increasing penetrations of distributed renewables are decreasing the effectiveness of traditional decentralized under-frequency load shedding (UFLS) schemes. As more distribution circuits begin to back-feed the transmission system, operation of traditional UFLS may exacerbate frequency instability. This paper presents the conceptual framework for a data-rich environment to coordinate UFLS across multiple distribution providers based on the laminar coordination framework in order to ensure optimal adaptive setting of UFLS relays. Communication and control are enabled through a distributed implementation of the IEC 61968-1 Common Information Model message bus structure. In addition to the proposed architecture, a novel adaptive UFLS scheme informed by a bi-level state estimator to create optimal relay setpoints is introduced. Initial simulation results are presented for the IEEE 14-bus test system on scenarios leading to mis-operation of traditional UFLS.

Anderson, Alexander A.↗

Software Validation Work With The ZPPR-15 Data

The analysis activities for fast reactors involve using many different pieces of software that are relied upon for their predictive capabilities. For this software to be considered reliable, documented proof that the predictions of the software are accurate is required. In this manuscript, the validation work that covers some of the Argonne software used in fast reactor design activities is discussed and displayed. This validation work includes neutron and gamma flux distributions, reaction rate distributions, and reactivity worth. In an ideal world, a reactor development program would have access to a comprehensive set of experimental facilities to help inform the design aspects of the reactor itself. While thermal-hydraulics experiments, and to a limited degree mechanical experiments, can be carried out today for validation needs, neutronics related experimental facilities are rather impractical because of the lack of experimental facilities. Given the desired time table for construction of new reactors, the reconstitution or creation of new neutronic experimental facilities is untenable and thus those reactor development programs must rely upon any available experimental measurements that are qualitatively similar to the design. While a methodology has been proposed to assess the similarity between the past experimental measurements and the reactor itself, that aspect is beyond the scope of this manuscript. In this manuscript, the focus is entirely placed on the analysis results for a series of experiments carried out at the ZPPR facility in Idaho in the mid-1980s. In this regard, this manuscript only shows the validation of the stated neutronics software for specific loadings of the ZPPR reactor. Because of the fuel form, its proposed enrichment, and the material content of the reactor core, the ZPPR-15 experiments were identified as potential validation data for the reactor. The ZPPR-15 experiments were intended as mockups of a 330 MWe Integral Fast Reactor program which was a follow on program to the Clinch River Breeder Reactor. In the ZPPR-15 series of experiments, measurements of the neutron spectrum, control rod worth, sodium void worth, foil reaction rate distributions, Doppler worth of heated samples, gamma dose, and axial expansion worth were all carried out and published. In many cases, these reactivity coefficients are good candidates to validate the reactivity coefficient calculation scheme used by the analysis software and included in the safety analysis activities of fast reactor development projects today. This manuscript discusses the modeling methodology and accuracy of the calculated experimental results using the LANL software MCNP and the ANL software package ARC (Argonne Reactor Codes). As will be shown, for many of the experimental measurements, the two software packages are found to be good predictive analysis tools for those experiments. In other cases, problems with the analysis methodology or underlying cross section data are exposed which indicates where predictive analysis is not as reliable. Finally, in some of the measurements the conclusion is reached that the experimental measurement cannot be reproduced with the analysis software as it is simply too difficult.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Software Validation Work With The ZPPR-15 Data

The analysis activities for fast reactors involve using many different pieces of software that are relied upon for their predictive capabilities. For this software to be considered reliable, documented proof that the predictions of the software are accurate is required. In this manuscript, the validation work that covers some of the Argonne software used in fast reactor design activities is discussed and displayed. This validation work includes neutron and gamma flux distributions, reaction rate distributions, and reactivity worth. In an ideal world, a reactor development program would have access to a comprehensive set of experimental facilities to help inform the design aspects of the reactor itself. While thermal-hydraulics experiments, and to a limited degree mechanical experiments, can be carried out today for validation needs, neutronics related experimental facilities are rather impractical because of the lack of experimental facilities. Given the desired time table for construction of new reactors, the reconstitution or creation of new neutronic experimental facilities is untenable and thus those reactor development programs must rely upon any available experimental measurements that are qualitatively similar to the design. While a methodology has been proposed to assess the similarity between the past experimental measurements and the reactor itself, that aspect is beyond the scope of this manuscript. In this manuscript, the focus is entirely placed on the analysis results for a series of experiments carried out at the ZPPR facility in Idaho in the mid-1980s. In this regard, this manuscript only shows the validation of the stated neutronics software for specific loadings of the ZPPR reactor. Because of the fuel form, its proposed enrichment, and the material content of the reactor core, the ZPPR-15 experiments were identified as potential validation data for the reactor. The ZPPR-15 experiments were intended as mockups of a 330 MWe Integral Fast Reactor program which was a follow on program to the Clinch River Breeder Reactor. In the ZPPR-15 series of experiments, measurements of the neutron spectrum, control rod worth, sodium void worth, foil reaction rate distributions, Doppler worth of heated samples, gamma dose, and axial expansion worth were all carried out and published. In many cases, these reactivity coefficients are good candidates to validate the reactivity coefficient calculation scheme used by the analysis software and included in the safety analysis activities of fast reactor development projects today. This manuscript discusses the modeling methodology and accuracy of the calculated experimental results using the LANL software MCNP and the ANL software package ARC (Argonne Reactor Codes). As will be shown, for many of the experimental measurements, the two software packages are found to be good predictive analysis tools for those experiments. In other cases, problems with the analysis methodology or underlying cross section data are exposed which indicates where predictive analysis is not as reliable. Finally, in some of the measurements the conclusion is reached that the experimental measurement cannot be reproduced with the analysis software as it is simply too difficult.

Aliberti, Gerardo↗

Hybrid model predictive control techniques for safety factor profile and stored energy regulation while incorporating NBI constraints

Abstract A novel hybrid Model Predictive Control (MPC) algorithm has been designed for simultaneous safety factor ( q ) profile and stored energy ( w ) control while incorporating the pulse-width-modulation constraints associated with the neutral beam injection (NBI) system. Regulation of the q -profile has been extensively shown to be a key factor for improved confinement as well as non-inductive sustainment of the plasma current. Simultaneous control of w is necessary to prevent the triggering of pressure-driven magnetohydrodynamic instabilities as the controller shapes the q profile. Conventional MPC schemes proposed for q -profile control have considered the NBI powers as continuous-time signals, ignoring the discrete-time nature of these actuators and leading in some cases to performance loss. The hybrid MPC scheme in this work has the capability of incorporating the discrete-time actuator dynamics as additional constraints. In nonlinear simulations, the proposed hybrid MPC scheme demonstrates improved q -profile+ w control performance for NSTX-U operating scenarios.

Physics↗

Resolution requirements for high-speed droplet impact along a rigid wall

The mechanical and thermal loads that are produced on high-speed projectiles from a droplet are key in the projectile development.To accurately predict the pressure produced along the wall surface during impact, it is important to correctly model the air/water interface of the droplet. Novel modeling approaches such as the Phase-Field enable control of the (numerical) interface thickness. However, the interplay between the grid resolution (including potential adaptive mesh refinement) and the interface thickness is unclear. In this work, the resolution requirements for numerical simulations of a two-dimensional cylindrical water droplet impacting a rigid wall at Mach four are computationally investigated. Numerical simulations are performed using a second-order accurate scheme with adaptive mesh refinement and a Phase-Field interface model. The choice of numerical interface thickness and resolution parameters affects the maximum pressure and location of the maximum pressure along the wall. The overall maximum pressure does not converge even at high resolutions. The number of levels and thickness (number of cells across the interface) play a role in obtaining a well behaved solution.

Burrell, E.↗

Event-Driven Predictive Approach for Real-Time Volt/VAR control with CVR in solar PV rich Active Distribution Network

The focus of this paper is on analyzing the impact of conservation voltage reduction in the presence of active devices such as solar photovoltaic (PV) and developing controls that leverage these distributed energy resources. An event-driven predictive approach for real-time volt/volt-ampere reactive (VAR) optimization, along with local two-level adaptive volt/VAR droop-based control algorithm for advanced distribution management systems, is introduced. The methodology covers aggregated and autonomous controls under different timescale operations, including the impact and effect of unpredicted events such as cloud transients on PV power production. In addition, the control schemes include the uncertainties in PV power generation and load power demand. The proposed methodology is validated in a real-time framework using the real-time digital simulator platform through co-simulation with models based on Python and OpenDSS (Open Distribution System Simulator). The developed methodology is tested on the modified IEEE 123-feeder test system. The results reveal that the proposed methodology works well in the presence of high penetrations of PV power, produces significant energy savings, and mitigates over-/undervoltage problems.

14 SOLAR ENERGY↗

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗

Hierarchical Model-Free Transactive Control of Building Loads to Support Grid Services

Residential buildings consume 4.4 quads of electricity annually, approximately 37% of the total electricity consumption in the United States. This represents a vast resource that can be used for demand management and other ancillary services. This project aims to develop a robust, scalable hierarchical transactional control mechanism incorporating elements of model-free control (MFC) and game theory to harness buildings to provide ancillary services to the grid. This approach is being taken to address the challenges of incorporating traditional transactional control schemes into existing buildings. The challenges include small individual building sizes requiring aggregation of many buildings, unpredictable energy usage that makes model identification difficult, and satisfying the sensitive occupant comfort constraints. In the proposed approach, by separating the control mechanism into two layers above and below the load aggregator, MFC can be used below the aggregator to modulate flexible building loads in response to pricing signals with guaranteed performance. This allows the burden of identifying an accurate model of the system to be shifted to the above-aggregator layer, where fluctuations in individual building usage have less impact on predicted building system behavior. Game theory concepts can then be used to determine pricing curves and control signals among regional aggregators. Managing this control in a game-theoretic approach will allow us to build in financial incentives that increase customer engagement. Additionally, the usage of MFC necessitates less burdensome computational and communication requirements, thus, it is easily deployable on small, embedded devices. In a broader sense, developing a strategy capable of effectively incorporating residential and small commercial buildings will allow greater throughput of existing and emerging grid services in addition to future transactive energy grid management methods. Using MFC within a hierarchical control architecture will allow the shifting of existing forecasting challenges to an aggregate level, where dynamics are slower and more predictable. This will enable a smooth interface between the grid services requests of utilities and the reliable control required by participating buildings. MFC, which supports distributed control architecture, permits a scalable solution that can be deployed to neighborhood-size systems as well as individual buildings. This project focuses on three objectives: (1) developing the mathematical framework, algorithm toolkit, and software toolset of the two-layer transactive control testbed; (2) developing a scalable solution for application over many residential and small-size commercial buildings with sparse distributed communication; and (3) field testing and implementation on hardware of the control strategies developed in the previous two objectives. The research and development activities are focused and designed to be impactful within the relevant 2025 targets timeframe. An open-source control framework for exploiting variability and dispatchability of building loads will be delivered as the outcome of the project. This capability enables greater participation of loads in electricity markets and ancillary services that are both useful for the utility and financially beneficial for building owners.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Identification and Online Updating of Dynamic Models for Demand Response of an Industrial Air Separation Unit

Demand-response operation of air separation units requires frequent changes in production rate(s), and scheduling calculations must explicitly consider process dynamics to ensure feasibility of the solutions. To this end, scale-bridging models (SBMs) approximate the scheduling-relevant dynamics of a process and its controller in a low-order representation. In contrast to previous works that have employed nonlinear SBMs, this paper proposes linear SBMs, developed using time-series analysis, to facilitate online scheduling computations. Using a year-long industrial dataset, we nd that compact linear SBMs are suitable approximations over typical scheduling horizons, but that their accuracies are unpredictable over time. We introduce a strategy for online updating of the SBMs, based on Kalman ltering schemes for online parameter estimation. The approach greatly improves the accuracy of SBM predictions and will enable the use of linear SBM-based demand-response scheduling in the future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗