scholarly journals Assimilation of Dynamic Combined Finite Discrete Element Methods Using the Ensemble Kalman Filter

2021 ◽  
Vol 11 (7) ◽  
pp. 2898
Author(s):  
Humberto C. Godinez ◽  
Esteban Rougier

Simulation of fracture initiation, propagation, and arrest is a problem of interest for many applications in the scientific community. There are a number of numerical methods used for this purpose, and among the most widely accepted is the combined finite-discrete element method (FDEM). To model fracture with FDEM, material behavior is described by specifying a combination of elastic properties, strengths (in the normal and tangential directions), and energy dissipated in failure modes I and II, which are modeled by incorporating a parameterized softening curve defining a post-peak stress-displacement relationship unique to each material. In this work, we implement a data assimilation method to estimate key model parameter values with the objective of improving the calibration processes for FDEM fracture simulations. Specifically, we implement the ensemble Kalman filter assimilation method to the Hybrid Optimization Software Suite (HOSS), a FDEM-based code which was developed for the simulation of fracture and fragmentation behavior. We present a set of assimilation experiments to match the numerical results obtained for a Split Hopkinson Pressure Bar (SHPB) model with experimental observations for granite. We achieved this by calibrating a subset of model parameters. The results show a steady convergence of the assimilated parameter values towards observed time/stress curves from the SHPB observations. In particular, both tensile and shear strengths seem to be converging faster than the other parameters considered.

2016 ◽  
Vol 2016 ◽  
pp. 1-10 ◽  
Author(s):  
Tu-bing Yin ◽  
Kang Peng ◽  
Liang Wang ◽  
Pin Wang ◽  
Xu-yan Yin ◽  
...  

The dynamic failure characteristics of coal rock exposed to high temperatures were studied by using a split Hopkinson pressure bar (SHPB) system. The relationship between energy and time history under different temperature conditions was obtained. The energy evolution and the failure modes of specimens were analyzed. Results are as follows: during the test, more than 60% of the incident energy was not involved in the breaking of the sample, while it was reflected back. With the increase of temperature, the reflected energy increased continuously; transmitted and absorbed energy showed an opposite variation. At the temperature of 25 to 100°C, the absorbed energy was less than that transmitted, while this phenomenon was opposite after 100°C. The values of specific energy absorption (SEA) were distributed at 0.04 to 0.1 J·cm−3, and its evolution with temperature could be divided into four different stages. Under different temperature conditions, the failure modes and the broken blocks of the samples were obviously different, combining with the variation of microstructure characteristics of coal at high temperatures; the physical mechanism of damage and failure patterns of coal rock are explained from the viewpoint of energy.


2009 ◽  
Vol 6 (4) ◽  
pp. 8279-8309 ◽  
Author(s):  
W. Ju ◽  
S. Wang ◽  
G. Yu ◽  
Y. Zhou ◽  
H. Wang

Abstract. Soil and atmospheric water deficits have significant influences on CO2 and energy exchanges between the atmosphere and terrestrial ecosystems. Model parameterization significantly affects the ability of a model to simulate carbon, water, and energy fluxes. In this study, an ensemble Kalman filter (EnKF) and observations of gross primary productivity (GPP) and latent heat (LE) fluxes were used to optimize model parameters significantly affecting the calculation of these fluxes for a subtropical coniferous plantation in southeastern China. The optimized parameters include the maximum carboxylation rate (Vcmax), the Ball-Berry coefficient (m) and the coefficient determining the sensitivity of stomatal conductance to atmospheric water vapor deficit D0). Optimized Vcmax and m showed larger seasonal and interannual variations than D0. Seasonal variations of Vcmax and m are more pronounced than the interannual variations. Vcmax and m are associated with soil water content (SWC). During dry periods, SWC at the 20 cm depth can explain 61% and 64% of variations of Vcmax and m, respectively. EnKF parameter optimization improves the simulations of GPP, LE and sensible heat (SH), mainly during dry periods. After parameter optimization using EnKF, the variations of GPP, LE and SH explained by the model increased by 1% to 4% at half-hourly steps and by 3% to 5% at daily time steps. Efforts are needed to develop algorithms that can properly describe the variations of these parameters under different environmental conditions.


2014 ◽  
Vol 919-921 ◽  
pp. 1924-1929 ◽  
Author(s):  
Husain Abbas ◽  
Tarek Almusallam ◽  
Yousef Al-Salloum

The strategic concrete structures are often required to resist impact loads arising from the projectile strike, falling weight, blast generated missile etc. The existing structures found deficient in resisting these loads are required to be retrofitted whereas the upcoming structures are required to be designed for expected impact loads. This paper explores the ways of strengthening existing reinforced concrete (RC) structures using externally bonded carbon fiber reinforced polymer (CFRP) sheets and improving the impact resistance of concrete by mixing hybrid fibers in its production. The impact response of concrete structures is assessed using experiments involving the impact of projectiles of different nose shapes on slab specimens. The material behavior at high strain rate is established using split Hopkinson pressure bar (SHPB) testing at varying strain rates. Analytical models are developed for predicting penetration depth, scabbing thickness, ballistic limit velocity and ejected mass. The experimental results were also validated through numerical modeling using LS-DYNA.


2016 ◽  
Vol 20 (12) ◽  
pp. 4949-4961 ◽  
Author(s):  
Chao Deng ◽  
Pan Liu ◽  
Shenglian Guo ◽  
Zejun Li ◽  
Dingbao Wang

Abstract. Hydrological model parameters play an important role in the ability of model prediction. In a stationary context, parameters of hydrological models are treated as constants; however, model parameters may vary with time under climate change and anthropogenic activities. The technique of ensemble Kalman filter (EnKF) is proposed to identify the temporal variation of parameters for a two-parameter monthly water balance model (TWBM) by assimilating the runoff observations. Through a synthetic experiment, the proposed method is evaluated with time-invariant (i.e., constant) parameters and different types of parameter variations, including trend, abrupt change and periodicity. Various levels of observation uncertainty are designed to examine the performance of the EnKF. The results show that the EnKF can successfully capture the temporal variations of the model parameters. The application to the Wudinghe basin shows that the water storage capacity (SC) of the TWBM model has an apparent increasing trend during the period from 1958 to 2000. The identified temporal variation of SC is explained by land use and land cover changes due to soil and water conservation measures. In contrast, the application to the Tongtianhe basin shows that the estimated SC has no significant variation during the simulation period of 1982–2013, corresponding to the relatively stationary catchment properties. The evapotranspiration parameter (C) has temporal variations while no obvious change patterns exist. The proposed method provides an effective tool for quantifying the temporal variations of the model parameters, thereby improving the accuracy and reliability of model simulations and forecasts.


2020 ◽  
Vol 2020 ◽  
pp. 1-13
Author(s):  
Rongrong Zhang ◽  
Dongdong Ma ◽  
Qingqing Su ◽  
Kun Huang

RMT-150B rock mechanics and split Hopkinson pressure bar (SHPB) devices were adopted to investigate the physical and mechanical properties, energy dissipation, and failure modes of argillaceous sandstone after different high temperatures under air-dried and saturation states. In addition, SEM and EDS tests were conducted to investigate its microstructure characteristics. Results showed that both the P-wave velocity and density of argillaceous sandstone specimen decreased with the increase of high temperature, while its porosity increased. Compared with static stress-strain curves, there was no obvious compaction stage for dynamic stress-strain curves, and the decrease rate of dynamic curves after peak strain was obviously slow compared with static curves. Both the static and dynamic strengths of argillaceous sandstone specimens decreased with increasing temperature, and the critical temperature point for the strength of argillaceous sandstone was 400°C. At the same temperature, the specific energy absorption under air-dried state was generally smaller compared with that under saturated state. Both the strain rate and temperature showed significant effect on the failure mode. After 100∼1000°C heat treatment, the granular crystals of the clastic structure gradually became larger, and both the number and average size of the original pores decreased, resulting in the deterioration of mechanical properties of argillaceous sandstone specimen.


SPE Journal ◽  
2010 ◽  
Vol 15 (02) ◽  
pp. 382-394 ◽  
Author(s):  
Haibin Chang ◽  
Yan Chen ◽  
Dongxiao Zhang

Summary In reservoir history matching or data assimilation, dynamic data, such as production rates and pressures, are used to constrain reservoir models and to update model parameters. As such, even if under certain conceptualization the model parameters do not vary with time, the estimate of such parameters may change with the available observations and, thus, with time. In reality, the production process may lead to changes in both the flow and geomechanics fields, which are dynamically coupled. For example, the variations in the stress/strain field lead to changes in porosity and permeability of the reservoir and, hence, in the flow field. In weak formations, such as the Lost Hills oil field, fluid extraction may cause a large compaction to the reservoir rock and a significant subsidence at the land surface, resulting in huge economic losses and detrimental environmental consequences. The strong nonlinear coupling between reservoir flow and geomechanics poses a challenge to constructing a reliable model for predicting oil recovery in such reservoirs. On the other hand, the subsidence and other geomechanics observations can provide additional insight into the nature of the reservoir rock and help constrain the reservoir model if used wisely. In this study, the ensemble-Kalman-filter (EnKF) approach is used to estimate reservoir flow and material properties by jointly assimilating dynamic flow and geomechanics observations. The resulting model can be used for managing and optimizing production operations and for mitigating the land subsidence. The use of surface displacement observations improves the match to both production and displacement data. Localization is used to facilitate the assimilation of a large amount of data and to mitigate the effect of spurious correlations resulting from small ensembles. Because the stress, strain, and displacement fields are updated together with the material properties in the EnKF, the issue of consistency at the analysis step of the EnKF is investigated. A 3D problem with reservoir fluid-flow and mechanical parameters close to those of the Lost Hills oil field is used to test the applicability.


Sign in / Sign up

Export Citation Format

Share Document