Porous medium hydraulic resistance

Storyboard

Darcy's Law considers hydraulic resistance, which in its basic form corresponds to that of a tube with a given length and radius. However, in many situations, the fluid flows through a medium containing pores rather than a single cavity. These pores act as capillaries, and their hydraulic resistance can be modeled as that of small tubes. The sum of these multiple hydraulic resistances in parallel forms the total hydraulic resistance of a porous material.

>Model

ID:(2071, 0)



Porous medium hydraulic resistance

Storyboard

Darcy's Law considers hydraulic resistance, which in its basic form corresponds to that of a tube with a given length and radius. However, in many situations, the fluid flows through a medium containing pores rather than a single cavity. These pores act as capillaries, and their hydraulic resistance can be modeled as that of small tubes. The sum of these multiple hydraulic resistances in parallel forms the total hydraulic resistance of a porous material.

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$l$
l
Capillary length
m
$r$
r
Capillary radius
m
$\Delta s$
Ds
Distance traveled in a time
m
$j_s$
j_s
Flux density
m/s
$\Delta h$
Dh
Height of liquid column
m
$R_h$
R_h
Hydraulic resistance
kg/m^4s
$\rho_w$
rho_w
Liquid density
kg/m^3
$N_p$
N_p
Number of equal hydraulic resistances in parallel
-
$N_s$
N_s
Number of identical hydraulic resistors in series
-
$V_p$
V_p
Pore volume
m^3
$f$
f
Porosity
-
$S$
S
Section Tube
m^2
$\rho_s$
rho_s
Solid Density
kg/m^3
$V_s$
V_s
Solid volume of a component
m^3
$\Delta t$
Dt
Time elapsed
s
$M_s$
M_s
Total Dry Mass of Sample
kg
$R_t$
R_t
Total hydraulic resistance
kg/m^4s
$R_{pt}$
R_pt
Total hydraulic resistance in parallel
kg/m^4s
$V_t$
V_t
Total volume
m^3
$\Delta L$
DL
Tube length
m
$R$
R
Tube radius
m
$\Delta p$
Dp
Variación de la Presión
Pa
$\eta$
eta
Viscosity
Pa s
$J_V$
J_V
Volume flow
m^3/s

Calculations


First, select the equation:   to ,  then, select the variable:   to 

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

 Variable   Given   Calculate   Target :   Equation   To be used



Equations

The volume flow ($J_V$) can be calculated from the hydraulic conductance ($G_h$) and the pressure difference ($\Delta p$) using the following equation:

equation=14471

Furthermore, using the relationship for the hydraulic resistance ($R_h$):

equation=15092

results in:

equation

Since the hydraulic resistance ($R_h$) is equal to the hydraulic conductance ($G_h$) as per the following equation:

equation=15092

and since the hydraulic conductance ($G_h$) is expressed in terms of the viscosity ($\eta$), the tube radius ($R$), and the tube length ($\Delta L$) as follows:

equation=15102

we can conclude that:

equation

If there is the pressure difference ($\Delta p$) between two points, as determined by the equation:

equation=4252

we can utilize the water column pressure ($p$), which is defined as:

equation=4250

This results in:

$\Delta p=p_2-p_1=p_0+\rho_wh_2g-p_0-\rho_wh_1g=\rho_w(h_2-h_1)g$



As the height difference ($\Delta h$) is:

equation=4251

the pressure difference ($\Delta p$) can be expressed as:

equation

Flow is defined as the volume the volume element ($\Delta V$) divided by time the time elapsed ($\Delta t$), which is expressed in the following equation:

equation=4347

and the volume equals the cross-sectional area the section Tube ($S$) multiplied by the distance traveled the tube element ($\Delta s$):

equation=4346

Since the distance traveled the tube element ($\Delta s$) per unit time the time elapsed ($\Delta t$) corresponds to the velocity, it is represented by:

equation=4348

Thus, the flow is a flux density ($j_s$), which is calculated using:

equation


Examples


mechanisms

If the porous medium is modeled as a network of the hydraulic resistance ($R_h$) identical elements connected in parallel in groups of the number of equal hydraulic resistances in parallel ($N_p$), which are then summed in series as the number of identical hydraulic resistors in series ($N_s$):

In this way, the general parallel sum, which results in the total hydraulic resistance in parallel ($R_{pt}$) according to

equation=3181,

is transformed into

equation=15902.

Similarly, the general series sum, which yields the total hydraulic resistance in series ($R_{st}$) according to

equation=3180,

is transformed into

equation=15903.

By using the definition of the hydraulic resistance ($R_h$) with the values the viscosity ($\eta$), the tube length ($\Delta L$), and the tube radius ($R$) according to

equation=3629,

and calculating the number of equal hydraulic resistances in parallel ($N_p$) from the porosity ($f$) and the section Tube ($S$) using

equation=15905,

as well as the number of identical hydraulic resistors in series ($N_s$) with the tube length ($\Delta L$) and the capillary length ($l$) through

equation=15904,

the final result is obtained as

equation=15907.


model

The number of identical hydraulic resistors in series ($N_s$) can be obtained by dividing the tube length ($\Delta L$) by the capillary length ($l$):

kyon

The porosity ($f$) is the ratio of the empty section through which the liquid flows relative to the section Tube ($S$). The first is calculated by multiplying the section of each capillary, the tube radius ($R$), by the number of equal hydraulic resistances in parallel ($N_p$), so that:

kyon

The number of equal hydraulic resistances in parallel ($N_p$) is calculated as the fraction given by the porosity ($f$) of the section Tube ($S$), divided by the section of a capillary, the tube radius ($R$):

kyon

The total hydraulic resistance in parallel ($R_{pt}$) is the result for the hydraulic resistance ($R_h$) identical values, obtained by dividing this by the number of equal hydraulic resistances in parallel ($N_p$):

kyon

Since the hydraulic resistance ($R_h$) is equal to the inverse of the hydraulic conductance ($G_h$), it can be calculated from the expression of the latter. In this way, we can identify parameters related to geometry (the tube length ($\Delta L$) and the tube radius ($R$)) and the type of liquid (the viscosity ($\eta$)), which can be collectively referred to as a hydraulic resistance ($R_h$):

kyon

The total hydraulic resistance ($R_t$) is calculated from a type of hydraulic resistance density, which depends on the viscosity ($\eta$), the porosity ($f$), and the tube radius ($R$), as well as the geometric factors the tube length ($\Delta L$) and the section Tube ($S$):

kyon

The height difference, denoted by the height difference ($\Delta h$), implies that the pressure in both columns is distinct. In particular, the pressure difference ($\Delta p$) is a function of the liquid density ($\rho_w$), the gravitational Acceleration ($g$), and the height difference ($\Delta h$), as follows:

kyon

If we have a tube with a the section Tube ($S$) moving a distance the tube element ($\Delta s$) along its axis, having displaced the volume element ($\Delta V$), then it is equal to:

kyon

A flux density ($j_s$) can be expressed in terms of the volume flow ($J_V$) using the section or Area ($S$) through the following formula:

kyon

Darcy rewrites the Hagen Poiseuille equation so that the pressure difference ($\Delta p$) is equal to the hydraulic resistance ($R_h$) times the volume flow ($J_V$):

kyon

The porosity ($f$) expresses the relationship between the pore volume ($V_p$) and the total volume ($V_t$), allowing us to define the equation as follows:

kyon

The surface of a disk ($S$) of ERROR:5275.1 is calculated as follows:

kyon

The mean Speed ($\bar{v}$) can be calculated from the distance traveled in a time ($\Delta s$) and the time elapsed ($\Delta t$) using:

kyon

The total volume ($V_t$) is the sum of the pore volume ($V_p$), which includes both the micropores and the macropores in the soil, and the total Dry Mass of Sample ($M_s$), so that:

kyon

Since we already know the total Dry Mass of Sample ($M_s$) and the solid volume ($V_s$) from the sample, we can introduce the solid Density ($\rho_s$) and calculate it using the following equation:

kyon


>Model

ID:(2071, 0)