skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Modeling of Rock Scour using Coupled 3-D Discrete Element and Lattice Boltzmann Methods
Scour of rock in unlined rock spillway channels is a critical issue facing many of the world’s dams. From a modeling point of view this poses a challenging and interesting problem that combines rock mechanics and hydraulics of turbulent flow. We analyze this interaction between the blocky rock mass and water by directly modeling the solid and fluid phases—the individual polyhedral blocks are modeled using the Discrete Element Method (DEM) while the water is modeled using the Lattice Boltzmann Method (LBM). The LBM mesh is entirely independent of the DEM discretization, making it possible to refine the LBM mesh such that transient and varied fluid pressures acting of the rock surface are directly modeled. This provides the capability to investigate the effect of water pressure inside the fractured rock mass, along potential sliding planes, and can be extended to rock falls and slides into standing bodies of water such as lakes and reservoirs. Herein we present preliminary results to demonstrate the capabilities of the methodology.  more » « less
Award ID(s):
1363354
PAR ID:
10066919
Author(s) / Creator(s):
;
Date Published:
Journal Name:
52nd US Rock Mechanics/Geomechancs Symposium
ISSN:
0085-574X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Geophysical Journal International)
    null (Ed.)
    SUMMARY The lattice Boltzmann method (LBM) is a method to simulate fluid dynamics based on modelling distributions of particles moving and colliding on a lattice. The Python scripting language provides a clean programming paradigm to develop codes based on the LBM, however in order to reach performance comparable to compiled languages, it needs to be carefully implemented, maximizing its vectorized tools, mostly integrated in the NumPy module. We present here the details of a Python implementation of a concise LBM code, with the purpose of offering a pedagogical tool for students and professionals in the geosciences who are approaching this technique for the first time. The first half of the paper focuses on how to vectorize a 2-D LBM code and show how if carefully done, this allows performance close to a compiled code. In the second part of the paper, we use the vectorization described earlier to naturally write a parallel implementation using MPI and test both weak and hard scaling up to 1280 cores. One benchmark, Poiseuille flow and two applications, one on sound wave propagation and another on fluid-flow through a simplified model of a rock matrix are finally shown. 
    more » « less
  2. ; (, Journal of Fluid Mechanics)
    Modelling of fluid–particle interactions is a major area of research in many fields of science and engineering. There are several techniques that allow modelling of such interactions, among which the coupling of computational fluid dynamics (CFD) and the discrete element method (DEM) is one of the most convenient solutions due to the balance between accuracy and computational costs. However, the accuracy of this method is largely dependent upon mesh size, where obtaining realistic results always comes with the necessity of using a small mesh and thereby increasing computational intensity. To compensate for the inaccuracies of using a large mesh in such modelling, and still take advantage of rapid computations, we extended the classical modelling by combining it with a machine learning model. We have conducted seven simulations where the first one is a numerical model with a fine mesh (i.e. ground truth) with a very high computational time and accuracy, the next three models are constructed on coarse meshes with considerably less accuracy and computational burden and the last three models are assisted by machine learning, where we can obtain large improvements in terms of observing fine-scale features yet based on a coarse mesh. The results of this study show that there is a great opportunity in machine learning towards improving classical fluid–particle modelling approaches by producing highly accurate models for large-scale systems in a reasonable time. 
    more » « less
  3. (, Journal of the Society of Materials Science, Japan)
    抄録 It is important to understand the long-term migration of radionuclides when considering long-lasting rock engineering projects such as the geological disposal of radioactive waste. The network of fractures and pores in a rock mass plays a major role in fluid migration as it provides pathways for fluid flow. The geometry of such a network can change due to fracture sealing by fine-grained material over extended periods of time. Groundwater commonly contains fine-grained material such as clay minerals, and it is probable that such minerals accumulate within rock fractures during groundwater flow, thereby decreasing fracture apertures and bulk permeability. It is therefore essential to conduct permeability measurements using water that includes fine-grained minerals in order to understand the evolving permeability characteristics of rock. However, this has not been studied to date in in-situ rock mass. Therefore, in the present study, we perform permeability measurements in a granite rock mass to investigate the change of permeability that occurs under the flow of water that includes clays. Our findings show that clay particles accumulate in fractures and that the permeability (hydraulic conductivity) of the granite rock mass decreases over time. The decrease was more significant in the earlier time. We conclude that the accumulation of clay minerals in the fracture decreases the permeability of a rock mass. Furthermore, we consider that the filling and closure of fractures in rock is possible under the flow of groundwater that contains clay minerals. 
    more » « less
  4. (, Transactions of the American Mathematical Society)
    Abstract. We consider a nonlinear, moving boundary, fluid-structure interaction problem between a time dependent incompressible, viscous fluid flow, and an elastic structure composed of a cylindrical shell supported by a mesh of elastic rods. The fluid flow is modeled by the time-dependent Navier- Stokes equations in a three-dimensional cylindrical domain, while the lateral wall of the cylinder is modeled by the two-dimensional linearly elastic Koiter shell equations coupled to a one-dimensional system of conservation laws defined on a graph domain, describing a mesh of curved rods. The mesh supported shell allows displacements in all three spatial directions. Two-way coupling based on kinematic and dynamic coupling conditions is assumed between the fluid and composite structure, and between the mesh of curved rods and Koiter shell. Problems of this type arise in many ap- plications, including blood flow through arteries treated with vascular prostheses called stents. We prove the existence of a weak solution to this nonlinear, moving boundary problem by using the time discretization via Lie operator splitting method combined with an Arbitrary Lagrangian-Eulerian approach, and a non-trivial extension of the Aubin-Lions-Simon compactness result to problems on moving domains. 
    more » « less
  5. ; ; (, SEG Meeting)
    Hydrothermal ore deposits involve complex fluid-rock interactions, where alteration zones yield information on physicochemical conditions of ore-forming processes and fluid pathways. Lithogeochemical vectors can be used to delineate alteration zones and potentially yield information for targeting ore zones for the exploration of mineral deposits. To be able to recognize these geochemical and mineralogical signatures, it is imperative to understand the underlying ore-forming processes and link them to geochemical vectors. Numerical modeling permits tracking changes in rock chemistry and metal precipitation mechanisms by simulating various alteration and fluid evolution scenarios such as boiling, fluid mixing, fluid-rock interaction, and changes in redox, pH, and temperature. The open-access MINES thermodynamic database (http://tdb.mines.edu) has been recently launched for modeling fluid-rock interaction and the chemical changes related to ore-forming processes, using the program GEM-Selektor (http://gems.web.psi.ch). The philosophy behind the MINES database is to focus on testing a series of numerical modeling projects used to simulate various ore deposits. We present a number of case studies where GEM-Selektor has been successfully applied to simulate ore-forming processes. These include Cu transport and mineralization at the Kansanshi Cu-Au deposit hosted in metasedimentary rocks, and the metasomatism of pegmatites and mobilization of rare earth elements (REEs) in the Strange Lake REE-Zr-Nb mineral deposit. We also use this method to demonstrate the coupling of complex chemical equilibria with reactive mass transport along a fluid flow path for simulating the deposition of Pb and Zn in a Mississippi Valley-type (MVT) deposit and the alteration associated with volcanogenic massive sulfide (VMS) deposits. The long-term goal of this project is the generation of a geochemical vectors database for interpreting field data, using knowledge gained from numerical simulations combined with field observations. This link between fundamental and applied research is expected to significantly advance ore vectoring for the exploration industry. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email