Discrete Fracture in Rocks

April 5, 2013

Have you ever noticed how water can flow through rocks and leave a trace of its passage by covering the surface with a patina, like white stripes? The whole process is quite complex but can roughly be explained by breaking it down into two coupled effects: gravity and chemical reactions. Gravity causes the water to infiltrate through discrete fractures (hard rocks) or pores (sedimentary rocks), until it charges the groundwater supply below to a point where excess water will surface on the face of the rock, like an overflowing glass of water. Chemical reactions, on the other hand, allow water to dissolve carbonates and similar substances present in the rock. Such chemical compounds will start sticking and accumulating on the surface.

White stripes due to discrete fractures. Image attribution: Cinzia Iacovelli
White stripes due to discrete fractures in the rocks
at Halibut Point State Park, Cape Ann, MA, United States.

Water Flowing Through Rocks

The phenomena of water flowing through rocks can be witnessed in nature, such as in the Halibut Point State Park in Cape Ann, MA (shown above). These rocks are made of 400+ year-old layers of granite, one of the most common types of rock, named for its granular texture. As you can see in the picture, the flux of water in those rocks is visible on their surface, giving way to streaks of a lighter color on their surface. The free flow of fluids through a fracture, in a block of solid or porous media, is what’s called a discrete fracture.

Simulating Discrete Fracture

If you are interested in simulating discrete fractures, you can use COMSOL Multiphysics together with the Subsurface Flow Module. In fact, our Model Gallery features a great Discrete Fracture documented example. The image below to the left shows the geometry used in this model, which consists of two blocks sharing a common interior boundary representing the discrete fracture. To the right you can see the analysis of time-dependent flow in the fractured porous blocks after 1000 seconds.

Geometry of the discrete fracture Analysis of time-dependent flow
Geometry of the discrete fracture. Analysis of time-dependent flow after 1000 s.

The two blocks are made of porous media, while the shared interior boundary represents a discrete fracture. The Darcy’s law physics interface is used to solve both the porous media flow and the fracture flow. The fracture is represented as an interior boundary, meaning Darcy’s law is solved on the surface instead of in a volume by adding the Fracture Flow feature to the Darcy’s law interface. Doing so avoids creating a very thin domain representing the fracture that would lead to a dense mesh consisting of very small elements. In other words, this is a very efficient way of solving this type of subsurface flow problem.

The isosurfaces in the plot to the right show how pressure is continuous across the fracture, a characteristic of subsurface flows, while arrows indicate the velocity field along the fracture. The fluid is moving from the inlet (right) to the outlet (left) because of an imposed pressure difference.

Further Reading

Now we know that water will surface thanks to discrete fractures in the rock and white stripes will appear because of chemical reactions happening on its external surface.

You will find more details on how to run a discrete fracture simulation yourself in the documentation available in our Model Gallery.


Comments (2)

Leave a Comment
Log In | Registration
Loading...
Dennys Armando Lopez-Falcon
Dennys Armando Lopez-Falcon
May 27, 2019

Thank you for your interesting contribution about modeling fractures in porous media.
I followed the documentation in the Model Gallery but I suspect that I am doing something wrong because there is no response to fracture width (d_f). Have you got any suggestion?
Thanks in advance for your assistance.
Best regards.

Victoria McGuire
Victoria McGuire
June 29, 2023

Hi Dennys, I am currently trying to follow the documentation for this model from the Model Gallery but the .mph file containing this Discrete Fracture seems to be missing from the download options. If there is any possibly you have still this file on hand or know where to retrieve one I would greatly appreciate a response back or a copy that I can download. Apologizes for any inconvenience this brings.

EXPLORE COMSOL BLOG