Applications

Storyboard

The porosity of the soil allows rain or irrigation water to penetrate the soil and reach the napa. Therefore we must study how it can be modeled based on our geometric model as the water moves.

>Model

ID:(367, 0)



Mechanisms

Concept


ID:(15205, 0)



Flow density and hydraulic conductivity

Concept

The flux density ($j_s$) can be expressed in terms of the hydraulic conductivity ($K_s$), in the infinitesimal limit with the column height differential ($dh$) and the distance differential ($dx$), as follows:

$ j_s = - K_s \displaystyle\frac{ dh }{ dx }$



This means that the steeper the gradient or the steeper the terrain, the larger the flux density ($j_s$) will be, as illustrated in the graph:

None



The graph shows how bars with equal ERROR:10142,0 values have progressively smaller ERROR:10141,0 values, resulting in a decreasing ERROR:7220,0. Since the volume of the liquid is conserved, this can only be possible if there is another flow that compensates for this reduction in ERROR:7220,0. This could be a flow perpendicular to the one shown, for example, if the shorter bars are wider in a direction perpendicular to the graph.

This issue leads to the following:

The height $h$ of the liquid can only be calculated as a result of solving a differential equation, as it must meet the requirement that volume is conserved throughout the entire area where flow occurs.



Additionally, it is important to keep in mind that:

The negative sign reflects the fact that flow always goes from the higher to the lower elevation zone. If the slope is negative, the negative sign results in positive flow (from left to right), and conversely, if the slope is positive, the flow is negative (from right to left).

ID:(930, 0)



Flow equation in one dimension

Concept

If we study the one-dimensional case, describing the process along the $x$-axis, we can observe how the height of the column $\Delta h$ varies over a time interval $\Delta t$. In this case, a column with width $\Delta x$ will change its volume per unit length over time as $\Delta x \Delta h/\Delta t$. On the other hand, the amount of liquid entering along the column at $x$ is $h(x) j_s(x)$, while at $x+\Delta x$ it exits as $h(x+\Delta x) j_s(x+\Delta x)$:

None



Therefore, the variation of the height of the water column on the ground ($h$) over time is equal to the variation of the product of the height of the water column on the ground ($h$) and the flux density ($j_s$) at position:

$\displaystyle\frac{\partial h}{\partial t} = - \displaystyle\frac{\partial}{\partial x}( h j_s )$



Partial derivatives are similar to ordinary derivatives, with the difference that they are applied to functions that depend on more than one variable. In these cases, the partial derivative, denoted by the symbol $\partial$, reminds us of the typical derivative denoted by the letter $d$, but with the peculiarity that the variables not mentioned in the denominator are held constant.

ID:(2290, 0)



Flow into a channel

Concept

In the case of flow towards a channel, the system can be modeled in a one-dimensional manner, where the height of the water column on the ground ($h$) is a function of the position of the water column on the ground ($x$) representing the flux density ($j_s$), and it satisfies the condition

$ h j_s = h_0 j_{s0} $



with the flow at a reference point ($j_{s0}$) and the reference height of the water column ($h_0$) defining the water profile in the soil:

None



The key to this equation is that the product of the height of the water column on the ground ($h$) and the flux density ($j_s$) must always remain constant. In this sense, if the height of the water column on the ground ($h$) increases, the flux density ($j_s$) decreases, and vice versa. Moreover, the sign remains the same; hence, flow towards the channel, i.e., negative flow, will occur only when the groundwater level is higher than that of the channel. As the liquid approaches the channel, the groundwater level decreases, leading to an increase in flow density.

ID:(15104, 0)



Flow height solution towards a channel

Concept

The solution to the one-dimensional flow equation towards a channel, where the height of the water column on the ground ($h$) is calculated as a function of the reference height of the water column ($h_0$) and the position of the water column on the ground ($x$) at the channel's edge, along with the characteristic length of the flow in the ground ($s_0$), takes the following form:

$ \displaystyle\frac{ h }{ h_0 } = \sqrt{1 + \displaystyle\frac{ 2 x }{ s_0 }} $



This solution is graphically represented in terms of the additional factors $h/h_0$ and $x/s_0$ as follows:

None



The profile reveals that, away from the channel, the height of the water column is significantly high. However, due to water extraction by the channel, this height begins to decrease until it reaches the channel's edge. Dynamically, the flux density ($j_s$) determines the amount of water flowing into the channel, while the reference height of the water column ($h_0$) gradually adjusts until it reaches an equilibrium state. In other words, if the value of the reference height of the water column ($h_0$) is too low relative to the total amount of water arriving at the channel, it increases; and if it is too high, it decreases. In this way, the reference height of the water column ($h_0$) acquires the value that balances the amount of incoming water with the amount of water flowing through the channel.

ID:(15109, 0)



Flux density solution towards a channel

Concept

The solution obtained for the height and the parameters the flow at a reference point ($j_{s0}$) and the reference height of the water column ($h_0$) reveals that the flux density ($j_s$) is given by:

$ \displaystyle\frac{ j_s }{ j_{s0} } = \displaystyle\frac{1}{\sqrt{1 + \displaystyle\frac{ 2 x }{ s_0 }}} $



We can graphically represent the flux density ($j_s$) as a function of the additional factors $j_s/j_{s0}$ and $x/s_0$ as follows:

None



It is noticeable that the flux density ($j_s$) continues to increase as we approach the channel, as the height of the water column on the ground ($h$) decreases. This increase is necessary to maintain the flow velocity in the flux density ($j_s$) or, alternatively, to increase it.

ID:(15110, 0)



Flow from a channel

Concept

In the case where the flow emerges from the channel, a situation arises where the level of the height of the water column on the ground ($h$) must decrease as we move away from the channel, ensuring the existence of the pressure gradient that drives the flow. The problem is that if the flow rapidly moves within the medium, the height will tend to zero, and as a result, the flow will approach infinity, which doesn't make sense.



This means that there is no stationary solution in such a scenario, and the only solution is for the medium to fill up until it reaches the height of the channel, effectively becoming constant.

The question is whether there exists a non-trivial stationary situation that represents a real and interesting scenario. One possible case is when the level of the medium decreases to the point where it becomes lower than the column before the solution diverges. This case corresponds to the situation where the flow emerges at the surface, and there is no divergence in the solution. This would imply that a flow is generated that exits to the exterior at a certain point, with the risk of weakening the foundation and thereby destabilizing the medium, which acts as a dam.

ID:(4746, 0)



Situation that meets boundary conditions

Concept

If we consider a situation where the flow from the channel can emerge at the surface, we have a scenario where the flow enters and then exits the medium, making the solution viable.

Emerging at the surface simply implies that the height of the liquid column becomes higher than that of the surrounding medium. In fact, similar to the case of flow towards a channel, this would generate water on the surface, which, if not allowed to flow, would actually form a new channel.

In the case of flow from a channel, it is possible to model the system in a one-dimensional manner, where the height of the water column on the ground ($h$) is a function of the position of the water column on the ground ($x$) representing the flux density ($j_s$) and satisfying the condition:

$ h j_s = h_0 j_{s0} $



With the flow at a reference point ($j_{s0}$) and the reference height of the water column ($h_0$) defining the water profile in the soil, as shown in the following image:



The key to the equation lies in the fact that the product of the height of the water column on the ground ($h$) and the flux density ($j_s$) must remain constant at all times. In this regard, if the height of the water column on the ground ($h$) increases, the flux density ($j_s$) will decrease, and vice versa. Furthermore, the sign remains the same. Therefore, flow from the channel, i.e., positive flow, will occur only if the height of the channel is greater than that of the point where the flow emerges. As the liquid moves away from the channel, the height will decrease, and the flow density will increase.

ID:(4370, 0)



Flow height solution from a channel

Concept

The solution to the one-dimensional flow equation from a channel, in which the value of the height of the water column on the ground ($h$) is calculated as a function of the reference height of the water column ($h_0$) and the position of the water column on the ground ($x$) at the channel's edge, along with the characteristic length of the flow in the ground ($s_0$), takes the following form:

$ \displaystyle\frac{ h }{ h_0 } = \sqrt{1 - \displaystyle\frac{ 2 x }{ s_0 }} $



This solution is graphically represented in terms of the additional factors $h/h_0$ and $x/x_0$ as follows:



The profile reveals that the height decreases as one moves away from the channel to maintain a pressure gradient. However, a problem arises when the distance reaches half of the characteristic length of the flow in the ground ($s_0$), as the height of the column reaches zero, and there is no solution for greater distances (the argument of the square root becomes negative). In other words, for the solution to make sense, there must be a mechanism that removes liquid before reaching this critical distance.

ID:(4374, 0)



Flux density solution from a channel

Concept

The solution obtained for the height and the parameters the flow at a reference point ($j_{s0}$) and the reference height of the water column ($h_0$) shows us that the flux density ($j_s$) is equal to:

$ \displaystyle\frac{ j_s }{ j_{s0} } = \displaystyle\frac{1}{\sqrt{1 - \displaystyle\frac{ 2 x }{ s_0 }}} $



We can graphically represent the flux density ($j_s$) in terms of the additional factors $j_s/j_{s0}$ and $x/x_0$ as follows:



the flux density ($j_s$) continues to increase as we approach the channel, as the height of the water column on the ground ($h$) decreases. This increase is necessary to maintain the flow velocity in the flux density ($j_s$) or, alternatively, to increase it.

ID:(7827, 0)



Dam I - Mina Córrego do Feijão

Concept

An example that illustrates the effect of flow through the base in the case of a dam occurred at Dam 1 of the 'Córrego do Feijão' mine in Brumadinho, Minas Gerais, Brazil.

On January 25, 2019, Dam 1, which is shown in the center of the image, collapsed as depicted in images 1 to 6. Initially, the base began to move while the top started to sink. Eventually, a torrent of water emerged from the base as the entire structure collapsed. In the lower central image, you can see the situation after the dam had completely emptied from the side that contained it ([1], [2]):



The upper-left image shows the dam before the collapse, and the diagram explains how the water pushes against the base surface (blue arrows) and causes the center to collapse (beige arrow). The images show the structure again before the collapse (top-right photo), when the base is being forced, causing the upper part to collapse (bottom-left photo), and the resulting water flow at the base (bottom-right photo) [3]:



The dynamics are driven by the high pressure and flow that exist at the base, explaining the emergence of water through this path.

In this case, there were multiple signs of danger, leading to detailed satellite monitoring of the movement of multiple points for over a year. The points are indicated in the upper photo, and in the lower left image, you can see a detail of the base. Specifically, the points that experienced the most total displacement (Bs and Bp) are highlighted, which are also shown in the graph on the right. The graph also shows the amount of rainfall, which contributes to some extent but is not necessarily a key factor [4]:

This example aims to demonstrate how high pressure at the base, combined with a high water flow, contributes to the observed dynamics, without necessarily explaining when or how it became unstable. This will be explored further.[1] Google Earth Pro for Brumadinho, Minas Gerais, Brazil, January 2019 and February 2019[2] Vale S.A. cameras[3] Criminal Investigative Procedure No. MPMG-0090.19.000013-4, Police Investigation No. PCMG-7977979, MINISTÉRIO PÚBLICO DO ESTADO DE MINAS GERAIS[4] Deformations Prior to the Brumadinho Dam Collapse Revealed by Sentinel-1 InSAR Data Using SBAS and PSI Techniques, Fábio F. Gama, José C. Mura, Waldir R. Paradella, and Cleber G. de Oliveira, MDPI, Remote Sens. 2020, 12, 3664.

ID:(4378, 0)



Flow into a well

Concept

In the case of groundwater flow towards a well, the height of the water column on the ground ($h$) as a function of the radius from center of well ($r$) with the water well radius ($r_0$), the characteristic length of the flow in the ground ($s_0$), and the reference height of the water column ($h_0$) is represented by

$ r \displaystyle\frac{ dh^2 }{ dr } = 2 h_0 ^2\displaystyle\frac{ r_0 }{ s_0 } $



which defines the water profile in the ground:

ID:(4371, 0)



Flow height solution towards a well

Concept

The solution to the one-dimensional flow equation towards a well, in which the value of the height of the water column on the ground ($h$) is calculated as a function of the radius from center of well ($r$), the reference height of the water column ($h_0$), and the water well radius ($r_0$) at the well's edge, along with the characteristic length of the flow in the ground ($s_0$), takes the following form:

$ \displaystyle\frac{ h }{ h_0 } = \sqrt{1 + \displaystyle\frac{ 2 r_0 }{ s_0 }\ln\left(\displaystyle\frac{ r }{ r_0 }\right)} $



This solution is graphically represented in terms of the additional factors $h/h_0$ and $r/r_0$ for various $r_0/s_0$ as follows:



The profile reveals that, far from the well, the height of the water column is significantly high. However, due to water extraction by the well, this height begins to decrease until it reaches the well's edge. Dynamically, the flux density ($j_s$) determines the amount of water flowing towards the well, while the reference height of the water column ($h_0$) gradually adjusts itself to reach an equilibrium state. In other words, if the value of the reference height of the water column ($h_0$) is too low relative to the total amount of water reaching the well, it increases, and if it is too high, it decreases. In this way, the reference height of the water column ($h_0$) acquires the value that balances the amount of water arriving with the amount of water being extracted through the well.

ID:(10591, 0)



Flux density solution towards a well

Concept

The solution obtained for the height and parameters the flow at a reference point ($j_{s0}$) and the radius from center of well ($r$), the water well radius ($r_0$), the characteristic length of the flow in the ground ($s_0$) shows us that the flux density ($j_s$) is equal to:

$ \displaystyle\frac{ j_s }{ j_{s0} } = \displaystyle\frac{1}{\displaystyle\frac{ r }{ r_0 }\sqrt{1 + \displaystyle\frac{ 2 r_0 }{ s_0 }\ln\left(\displaystyle\frac{r}{r_0}\right)}}$



This solution is graphically represented in terms of the additional factors $j_s/j_{s0}$ and $r/r_0$ for various values of $r_0/s_0$ as follows:

None



the flux density ($j_s$) continues to increase as we approach the channel, while the height of the water column on the ground ($h$) decreases. This increase is necessary to maintain the flow velocity in the flux density ($j_s$) or, alternatively, to increase it.

ID:(2209, 0)



Model

Concept


ID:(15224, 0)