Rock mass response for lined rock caverns subjected to high internal gas pressure

The storage of hydrogen gas in underground lined rock caverns(LRCs)enables the implementation of the first fossil-free steelmaking process to meet the large demand for crude steel.Predicting the response of rock mass is important to ensure that gas leakage due to rupture of the steel lining does not occur.Analytical and numerical models can be used to estimate the rock mass response to high internal pressure;however,the fitness of these models under different in situ stress conditions and cavern shapes has not been studied.In this paper,the suitability of analytical and numerical models to estimate the maximum cavern wall tangential strain under high internal pressure is studied.The analytical model is derived in detail and finite element(FE)models considering both two-dimensional(2D)and three-dimensional(3D)geometries are presented.These models are verified with field measurements from the LRC in Skallen,southwestern Sweden.The analytical model is inexpensive to implement and gives good results for isotropic in situ stress conditions and large cavern heights.For the case of anisotropic horizontal in situ stresses,as the conditions in Skallen,the 3D FE model is the best approach.

Effect of rock joints on lined rock caverns subjected to high internal gaspressure

The storage of hydrogen gas in lined rock caverns(LRCs)may enable the implementation of the firstlarge-scale fossil-free steelmaking process in Sweden,but filling such storage causes joints in the rockmass to open,concentrating strains in the lining.The structural interaction between the LRC componentsmust be able to reduce the strain concentration in the sealing steel lining;however,this interaction iscomplex and difficult to predict with analytical methods.In this paper,the strain concentration in LRCsfrom the opening of rock joints is studied using finite element(FE)analyses,where the large-and small-scale deformation behaviors of the LRC are coupled.The model also includes concrete crack initiation anddevelopment with increasing gas pressure and rock joint width.The interaction between the jointed rockmass and the reinforced concrete,the sliding layer,and the steel lining is demonstrated.The results showthat the rock mass quality and the spacing of the rock joints have the greatest influence on the straindistributions in the steel lining.The largest effect of rock joints on the maximum strains in the steellining was observed for geological conditions of“good”quality rock masses.

Modelling erosion of a single rock block using a coupled CFD-DEM approach

Rock block removal is the prevalent physical mechanism for rock erosion and could affect the stability of dam foundations and spillways.Despite this,understanding of block removal is still inadequate because of the complex interactions among block characteristics,hydraulic forces,and erosive processes acting on the block.Herein,based on a previously conducted physical experiment of erosion of a single rock block,the removal processes of two different protruding blocks are represented by a coupled computational fluid dynamics-discrete element model(CFD-DEM)approach under varied flow conditions.Additionally,the blocks could be rotated with respect to the flow direction to consider the effect of the discontinuity orientation on the block removal process.Simulation results visualize the entire block removal process.The simulations reproduce the effects of the discontinuity orientation on the critical flow velocity inducing block incipient motion and the trajectory of the block motion observed in the physical experiments.The numerical results present a similar tendency of the critical velocities at different discontinuity orientations but have slightly lower values.The trajectory of the block in the simulations fits well with the experimental measurements.The relationship between the dimensionless critical shear stress and discontinuity orientation observed from the simulations shows that the effect of block protrusion becomes more dominant on the block incipient motion with the increase of relative protrusion height.To our knowledge,this present study is the first attempt to use the coupled finite volume method(FVM)-DEM approach for modelling the interaction behavior between the block and the flowing water so that the block removal process can be reproduced and analyzed.

Effects of spatial variation in cohesion over the concrete-rock interface on dam sliding stability

The limit equilibrium method （LEM） is widely used for sliding stability evaluation of concrete gravitydams. Failure is then commonly assumed to occur along the entire sliding surface simultaneously.However, the brittle behaviour of bonded concrete-rock contacts, in combination with the varying stressover the interface, implies that the failure of bonded dam-foundation interfaces occurs progressively. Inaddition, the spatial variation in cohesion may introduce weak spots where failure can be initiated.Nonetheless, the combined effect of brittle failure and spatial variation in cohesion on the overall shearstrength of the interface has not been studied previously. In this paper, numerical analyses are used toinvestigate the effect of brittle failure in combination with spatial variation in cohesion that is taken intoaccount by random fields with different correlation lengths. The study concludes that a possible existenceof weak spots along the interface has to be considered since it significantly reduces the overallshear strength of the interface, and implications for doing so are discussed.

Influence of location of large-scale asperity on shear strength of concrete-rock interface under eccentric load

The location and geometry of large-scale asperity present at the foundation of concrete gravity dams and buttress dams affect the shear resistance of the concrete-rock interface.However,the parameters describing the frictional resistance of the interface usually do not account for these asperities.This could result in an underestimate of the peak shear stre ngth,which leads to significantly conservative design for new dams or unnecessary stability enhancing measures for existing ones.The aim of this work was to investigate the effect of the location of first-order asperity on the peak shear strength of a concrete-rock interface under eccentric load and the model discrepancy associated with the commonly used rigid body methods for calculating the factor of safety(FS)against sliding.For this,a series of direct and eccentric shear tests under constant normal load(CNL)was carried out on concrete-rock samples.The peak shear strengths measured in the tests were compared in terms of asperity location and with the predicted values from analytical rigid body methods.The results showed that the large-scale asperity under eccentric load significantly affected the peak shear strength.Furthermore,unlike the conventional assumption of sliding or shear failure of an asperity in direct shear,under the effect of eccentric shear load,a tensile failure in the rock or in the concrete could occur,resulting in a lower shear strength compared with that of direct shear tests.These results could have important implications for assessment of the FS against sliding failure in the concrete-rock interface.

Tailor-Made Electrospun Culture Scaffolds Control Human Neural Progenitor Cell Behavior—Studies on Cellular Migration and Phenotypic Differentiation

In neuroscience research, cell culture systems are essential experimental platforms. It is of great interest to explore in vivo-like culture substrates. We explored how basic properties of neural cells, nuclei polarization, phenotypic differentiation and distribution/migration, were affected by the culture at poly-L-lactic acid (PLLA) fibrous scaffolds, using a multipotent mitogen-expanded human neural progenitor cell (HNPC) line. HNPCs were seeded, at four different surfaces: two different electrospun PLLA (d = 1.2 - 1.3 μm) substrates (parallel or random aligned fibers), and planar PLL- and PLLA surfaces. Nuclei analysis demonstrated a non-directed cellular migration at planar surfaces and random fibers, different from cultures at aligned fibers where HNPCs were oriented parallel with the fibers. At aligned fibers, HNPCs displayed the same capacity for phenotypic differentiation as after culture on the planar surfaces. However, at random fibers, HNPCs showed a significant lower level of phenotypic differentiation compared with cultures at the planar surfaces. A clear trend towards greater neuronal formation at aligned fibers, compared to cultures at random fibers, was noted. We demonstrated that the topography of in vivo-resembling PLLA scaffolds significantly influences HNPC behavior, proven by different migration behavior, phenotypic differentiation potential and nuclei polarization. This knowledge is useful in future exploration of in vivo-resembling neural cell system using electrospun scaffolds.