Register to enable "Calculate" button.
Units:
barrel=42 U.S. gallons, cm=centimeter, cP=centipoise, cSt=centistoke, darcy=permeability unit (1.01325x10^{12}
darcy=1 m^{2}), inch H_{2}O=inch water at 60F, inch Hg=inch mercury at 60F, ft=foot,
ft fluid=ft of whatever fluid is flowing through the bed, g=gram, gal=U.S. gallon, hr=hour, kg=kilogram,
kg/m^{2}s=kg/s of flow through 1 m^{2} of upstream area, lb=pound, lb/ft^{2}s=lb/s
of flow through 1 ft^{2} of upstream area, lb(f)=pound (force), m=meter, m fluid=meters of
whatever fluid is flowing through the bed, mbar=millibar, mm=millimeter, mm H_{2}O=mm
water at 4C, min=minute, N=Newton, psi=lb(f)/in^{2}, s=second
Introduction
Chemical engineers know this as a packed or porous bed. Groundwater hydrologists
call this a permeameter or flow through porous media. We provide two methods 
Idelchik and Darcy's law  since hydrogeologists and civil engineers are usually more
familiar with Darcy's law while chemical and mechanical engineers may be more familiar
with the Idelchik approach. The Idelchik method is valid for laminar or turbulent
flow through the bed while Darcy's law is valid only for laminar flow (Darcy Reynolds
Number<10). Another difference is that only Darcy's law requires entering the
permeability. Permeability is used rather than hydraulic conductivity since
permeability is a soil property independent of fluid while hydraulic conductivity is
defined for water only. If the Idelchik method is selected, permeability will be
backcalculated (based on Darcy's law) in case you want to compare a bed material with a
soil type. We also compute a minor loss coefficient (K_{m}) for flow through
the bed in case the bed is part of a longer pipeline that you are modeling. Any
liquid or gas can be used with either method so long as the fluid's density and viscosity
are known. Several fluids have properties built into the program.
Equations
Idelchik Method (Fried and Idelchik, 1989):
Valid for any porosity and Reynolds number. However, if porosity is
computed from packing angle, the packing angle equation is only valid for packing angles
between 60^{o} and 90^{o} (resulting in computed porosity between 0.2595
and 0.4764).
Darcy's Law Method (e.g. Freeze and Cherry, 1979; Bear, 1979):
Valid for Darcy Reynolds Number < 10 and any porosity. Porosity can
be computed from packing angle from Idelchik's equation above for packing angles between
60^{o} and 90^{o}.
For Both Methods:
For our calculation, all equations were solved analytically for the desired variables,
except when computing velocity and flowrate using the Idelchik method. This required
a numerical solution based on a cubic solver such as presented in Rao (1985).
Variables (Units: F=Force, L=Length,
M=Mass, T=Time)
A_{up}=Upstream area (column crosssectional area) [L^{2}].
A_{bed}=Bed area (column crosssectional area minus area of grains) [L^{2}].
d=Grain diameter [L].
dH/dx=Hydraulic gradient [unitless]. Change in head with respect to
distance.
k=Permeability of particles [L^{2}].
K=Hydraulic conductivity [L/T]. Not used in calculation but shown in
equations above. Instead of K, k is used in calculations since it is independent of
fluid type.
K_{m}=Minor loss coefficient of bed [unitless].
L=Bed length [L]. Also known as permeameter length.
M=Mass flow rate [M/T].
n=Porosity of bed. Enter as decimal. For example, enter 35% as 0.35.
Q=Volumetric flow rate [L^{3}/T].
Re(Idelchik)=Reynolds number defined by Idelchik method [unitless]. Based
on upstream velocity.
Re(Darcy)=Reynolds number for Darcy method [unitless]. Based on bed
velocity.
V_{up}=Upstream velocity [L/T]. Also known as Darcy velocity.
Upstream of porous bed.
V_{bed}=Bed velocity [L/T]. Velocity of fluid flowing through the
bed. Also known as pore velocity.
ρ=Fluid density [M/L^{3}].
μ=Fluid dynamic viscosity [FT/L^{2}] or [M/LT].
ΔP=Pressure loss [F/L^{2}]. Upstream Pressure  Downstream Pressure.
θ=Packing angle [radians or degrees]. Angle between successive rows of particles.
Soil Properties
The following table of soil properties was compiled from Freeze and Cherry
(1979), Bear (1979), Hillel (1982), and Sanders (1998).
Table of Soil Properties
Soil Type 
Permeability, k
(m^{2}) 
Porosity, n 
Grain Size, d (mm) 
Clayey 
10^{18}  10^{15} 
0.4  0.7 
< 0.002 
Silty 
10^{16}  10^{12} 
0.35  0.5 
0.002  0.05 
Sandy 
10^{14}  10^{10} 
0.25  0.5 
0.05  2 
Gravelly 
10^{10}  10^{7} 
0.25  0.4 
> 2 
Error Messages given by calculation
Input checks. If one of these messages appear, the
calculation will not proceed.
Lower limits on variables:
"A must be > 0", "Density must be > 0", "Viscosity must
be > 0", "Grain dia must be > 0", "L must be > 0",
"k must be > 0", "Pressure loss must be > 0", "M must be
> 0", "Q must be > 0", "V_{bed} must be > 0".
Upper limits on variables:
"A must be < 1e50 m^{2}", "Density must be < 1e6 kg/m^{3}",
"Viscosity must be < 1e9 m^{2}/s", "d must be < 1000 m",
"L must be < 1e10 m", "k must be < 100 m^{2}", "P
loss must be < 1e50 Pa", "M must be < 1e50 kg/s", "Q must be
< 1e50 m^{3}/s", "V_{bed} must be < 1e50 m/s".
Other limits on variables:
"Porosity out of range". Porosity must be between 0.00001 and
0.99999.
"Angle out of range". If packing angle is entered, it must be
between 60^{o} and 90^{o}.
"L must be ≥ d". Bed length cannot be less than grain size.
Runtime errors. The following messages may be
generated during execution:
"Infeasible input." May occur when using Idelchik method and
solving for flowrate using input values that approach machine precision (e.g. Pressure
loss of 10^{200} Pa) and are physically unrealistic.
"L will be < d." May appear if bed length is being computed,
and the length will be less than the grain diameter.
"Darcy Re >10." Will appear only if the Darcy method was
selected, and the Darcy Reynolds number is > 10. Darcy's law is only valid for
laminar flow which is defined as Re(Darcy)<10. This is only a warning message;
results will be computed and shown.
References
Bear, J. 1979. Hydraulics of Groundwater. McGrawHill Pub Co.
Freeze, R. A. and J. A. Cherry. 1979. Groundwater. Prentice Hall,
Inc.
Fried, E. and I. E. Idelchik. 1989. Flow Resistance: A Design Guide for
Engineers. Hemisphere Pub. Corp.
Hillel, D. 1982. Introduction to Soil Physics. Academic Press, Inc.
Rao, S. S. 1985. Optimization Theory and Applications. Wiley Eastern
Limited. 2ed.
Sanders, L. L. 1998. A Manual of Field Hydrogeology. Prentice Hall,
Inc.
© 20012015 LMNO Engineering, Research, and
Software, Ltd. All rights reserved.
Please contact us for consulting or other questions.
LMNO Engineering, Research, and Software, Ltd.
7860 Angel Ridge Rd. Athens, Ohio 45701 USA Phone: (740) 5921890
LMNO@LMNOeng.com http://www.LMNOeng.com

To:
LMNO Engineering home page (more calculations)
Related Calculations:
Compute gradient from well head measurements
1D Step injection with advection, dispersion, retardation
3D Pulse (slug) injection with advection, dispersion
Drawdown for unsteady groundwater flow to pumping well in confined aquifer (Theis calculator)
Transmissivity T=Kb (and K table)
Unit Conversions
Register
