InternalFluidFlow.jl

Installing and Loading InternalFluidFlow

InternalFluidFlow can be installed and loaded either from the JuliaHub repository (last released version) or from the maintainer's repository.

Last Released Version

The last version of InternalFluidFlow can be installed from JuliaHub repository:

using Pkg
Pkg.add("InternalFluidFlow")
using InternalFluidFlow

If InternalFluidFlow is already installed, it can be updated:

using Pkg
Pkg.update("InternalFluidFlow")
using InternalFluidFlow

Pre-Release (Under Construction) Version

The pre-release (under construction) version of InternalFluidFlow can be installed from the maintainer's repository.

using Pkg
Pkg.add(path="https://github.com/aumpierre-unb/InternalFluidFlow.jl")
using InternalFluidFlow

Citation of InternalFluidFlow

You can cite all versions (both released and pre-released), by using DOI 105281/zenodo.7019888. This DOI represents all versions, and will always resolve to the latest one.

The InternalFluidFlow Module for Julia

InternalFluidFlow provides the following functions:

Re2f
f2Re
h2fRe
doPlot

Re2f

Re2f computes the Darcy friction f factor given the Reynolds number Re and the relative roughness ε.

By default, pipe is assumed to be smooth (ε = 0). If ε > 0.05, relative roughness is reset to upper limit ε = 0.05.

If fig = true is given a schematic Moody diagram is plotted as a graphical representation of the solution.

If lam = false is given then Re2f disregards the laminar flow bounds (Re < 4e3).

If turb = false is given then Re2f disregards the turbulent flow bounds (Re > 2.3e3).

Syntax:

Re2f(; # Darcy friction factor
    Re::Number, # Reynolds number
    ε::Number=0, # relative roughness
    fig::Bool=false # show plot
    )

Examples:

Compute the Darcy friction factor f given the Reynolds number Re = 120,000 and the relative roughness ε = 3e-3.

Re2f( # Darcy friction factor
    Re=120e3, # Reynolds number
    ε=3e-3 # relative roughness
    )

Compute the Darcy friction factor f given the Reynolds number Re = 120,000 and the relative roughness ε = 6e-2. In this case, relative roughness is reassigned to ε = 5e-2 for turbulent flow.

Re2f( # Darcy friction factor
    Re=120e3, # Reynolds number
    ε=6e-2 # relative roughness
    )

Compute the Darcy friction factor f given the Reynolds number Re = 3,500 and the relative roughness ε = 6e-3 and show results on a schematic Moody diagram.

Re2f( # Darcy friction factor
    Re=3500, # Reynolds number
    ε=6e-3, # relative roughness
    fig=true # show plot
    )

f2Re

f2Re computes the Reynolds number Re given the Darcy friction factor f and the relative roughness ε for both laminar and turbulent regime, if possible.

By default, pipe is assumed to be smooth (ε = 0). If ε > 0.05, relative roughness is reset to upper limit ε = 0.05.

If fig = true is given a schematic Moody diagram is plotted as a graphical representation of the solution.

If lam = false is given then f2Re disregards the laminar flow bounds (Re < 4e3).

If turb = false is given then f2Re disregards the turbulent flow bounds (Re > 2.3e3).

Syntax:

f2Re(; # Reynolds number
    f::Number, # Darcy friction factor
    ε::Number=0, # relative roughness, default is smooth pipe
    fig::Bool=false, # default hide plot
    turbulent::Bool=false, # default disregard turbulent flow
    msgs::Bool=true # default not msgs messages
    )

Examples:

Compute the Reynolds number Re given the Darcy friction factor f = 2.8e-2 and the pipe relative roughness ε = 5e-3. In this case, only laminar solution is possible:

f2Re( # Reynolds number
    f=2.8e-2, # Darcy friction factor
    ε=5e-3 # relative roughness
    )

Compute the Reynolds number Re given the Darcy friction factor f = 1.8e-2 and the pipe relative roughness ε = 5e-3. In this case, only turbulent solution is possible:

f2Re( # Reynolds number
    f=1.8e-2, # Darcy friction factor
    ε=5e-3 # relative roughness
    )

Compute the Reynolds number Re given the Darcy friction factor f = 1.2e-2 and the pipe relative roughness ε = 9e-3. In this case, both laminar and turbulent solutions are impossible:

f2Re( # Reynolds number
    f=1.2e-2, # Darcy friction factor
    ε=9e-3 # relative roughness
    )

Compute the Reynolds number Re given the Darcy friction factor f = 0.028 for a smooth pipe and plot and show results on a schematic Moody diagram. In this case, both laminar and turbulent solutions are possible:

f2Re( # Reynolds number
    f=0.028, # Darcy friction factor
    fig=true # show plot
    )

h2fRe

h2fRe computes the Reynolds number Re and Darcy friction factor f given the head loss h, the pipe hydraulic diameter D or the flow speed v or the volumetric flow rate Q, the pipe length L (default L = 100 cm), the pipe roughness k (default k = 0 cm) or the pipe relative roughness ε (default ε = 0), the fluid density ρ (default ρ = 0.997 g/cc), the fluid dynamic viscosity μ (default μ = 0.0091 g/cm/s), and the gravitational accelaration g (default g = 981 cm/s/s).

By default, pipe is assumed to be 1 m long, L = 100 (in cm).

By default, pipe is assumed to be smooth (ε = 0). If ε > 0.05, relative roughness is reset to upper limit ε = 0.05.

Notice that default values are given in the cgs unit system and, if taken, all other parameters must as well be given in cgs units.

If fig = true is given a schematic Moody diagram is plotted as a graphical representation of the solution.

If flow speed is given, both laminar and turbulent flow bounds are considered for possible solutions.

Syntax:

h2fRe(; # Reynolds number Re and Darcy friction factor f
    h::Number, # head loss in cm
    L::Number=100, # pipe length in cm, default is 100 cm
    ε::Number=NaN, # pipe relative roughness
    k::Number=NaN, # pipe roughness in cm
    D::Number=NaN, # pipe hydraulic diameter in cm
    v::Number=NaN, # flow speed in cm/s
    Q::Number=NaN, # volumetric flow rate in cc/s
    ρ::Number=0.997, # fluid dynamic density in g/cc
    μ::Number=0.0091, # fluid dynamic viscosity in g/cm/s
    g::Number=981, # gravitational accelaration in cm/s/s
    fig::Bool=false # default is hide plot
    )

Examples:

Compute the Reynolds number Re and the Darcy friction factor f given the head loss h = 40 cm, the pipe hydraulic diameter D = 4 in, the pipe length L = 25 m and the pipe relative roughness ε = 0.0021 for water flow.

h2fRe( # Reynolds number Re and Darcy friction factor f
     h=40, # head loss in cm
     D=4*2.54, # pipe hyraulic diameter in cm
     L=25e2, # pipe length in cm
     ε=0.0027 # pipe relative roughness
     )

Compute the Reynolds number Re and the Darcy friction factor f given the head loss per meter h/L = 1.6 cm/m, the volumetric flow rate Q = 8.6 L/s, the pipe length L = 25 m, the pipe roughness k = 0.08 cm, the fluid density ρ = 0.989 g/cc and the fluid dynamic viscosity μ = 0.89 cP.

h2fRe( # Reynolds number Re and Darcy friction factor f
     h=1.6*25, # head loss in cm
     Q=8.6e3, # volumetric flow rate in cc/s
     L=25e2, # pipe length in cm
     k=0.08, # pipe relative roughness
     ρ=0.989, # fluid dynamic density in g/cc
     μ=8.9e-3 # fluid dynamic viscosity in g/cm/s
     )

Compute the Reynolds number Re and the Darcy friction factor f, given the head loss h = 0.30 m, the flow speed v = 25 cm/s, the pipe length L = 25 m, the pipe roughness 0.02 cm for water flow and show results on a schematic Moody diagram. In this case, both laminar and turbulent solutions are possible:

h2fRe( # Reynolds number Re and Darcy friction factor f
     h=0.30e2, # head loss in cm
     v=25, # flow speed in cm/s
     L=25e2, # pipe length in cm
     k=0.02, # pipe roughness in cm
     fig=true # show plot
     )

Compute the Reynolds number Re and the Darcy friction factor f, given the head loss h = 0.12 m, the flow speed v = 23 cm/s, the pipe length L = 25 m, the pipe roughness k = 0.3 cm for water flow and show results on a schematic Moody diagram. In this case, both laminar and turbulent solutions are possible, however laminar flow is extended to Re = 4e3 and relative roughness is reassigned to maximum ε = 5e-2 for tubulent flow:

h2fRe( # Reynolds number Re and Darcy friction factor f
     h=0.12e2, # head loss in cm
     v=23, # flow speed in cm/s
     L=25e2, # pipe length in cm
     k=0.3, # pipe roughness in cm
     fig=true # show plot
     )

doPlot

doPlot produces a schematic Moody diagram..

Syntax:

doPlot()

Examples:

Build a schematic Moody diagram.

doPlot()