The intention of this page is to provide Free / Open Source files for piping engineers
It consists of tools, standard reports and user developed functions (UDF), which are included in each application.
Versions for Excel (*.xls) are for the 2003 version. Few files are *.xlsm

2a. xls
Air receivers volume calculation.xls
(This file presents the calculation of an air receiver, and shows several references related to this theme.
Application example and derivation of equation to determine the receiver volume.
Rev. 07.08.2018)

2b. ods
Air receivers volume calculation.ods
(This file presents the calculation of an air receiver, and shows several references related to this theme.
Application example and derivation of equation to determine the receiver volume.
Rev. 07.08.2018)

3a. xls
Atmospheric temperature, pressure and density as function of the height above sea level.xls
(This file presents the calculation of atmospheric temperature, pressure and density as a function
of the heigth above sea level, according 1976 U.S. Standard Atmosphere. It includes also an
approximate method that can be applied for a range of heights 0 km.a.s.l. < H < 6 km.a.s.l.
with an error less than 0.1% . Also, it is presented an equation to calculate the water vapor.
pressure as a function of temperature. References included.
Rev. 09.08.2018)

3b. ods
Atmospheric temperature, pressure and density as function of the height above sea level.ods
(This file presents the calculation of atmospheric temperature, pressure and density as a function
of the heigth above sea level, according 1976 U.S. Standard Atmosphere. It includes also an
approximate method that can be applied for a range of heights 0 km.a.s.l. < H < 6 km.a.s.l.
with an error less than 0.1% . Also, it is presented an equation to calculate the water vapor.
pressure as a function of temperature. References included.
Rev. 09.08.2018)

4a. xls
Average particle size determination d50 from gravimetric analysis. Mesh-opening.xls.
(For a given granulometric analysis (mesh size vs. Retained percentage) a table of "Mesh vs. Particle
size (mesh opening)" is made. The plotted curve allows to find the particle size orresponding to
a 50% retained percentage: This is the average particle size or d50 value.
Evaluation of a filtrated sample properties.
Gravimetric composition of a flow resulting from the joint of two flows.
Rev. 09.08.2016)

4b. ods
Average particle size determination d50 from gravimetric analysis. Mesh-opening.ods.
(For a given granulometric analysis (mesh size vs. Retained percentage) a table of "Mesh vs. Particle
size (mesh opening)" is made. The plotted curve allows to find the particle size orresponding to
a 50% retained percentage: This is the average particle size or d50 value.
Evaluation of a filtrated sample properties.
Gravimetric composition of a flow resulting from the joint of two flows.
Rev. 09.08.2016)

6a. xls
Blower_Air_line.xls
(This file presents a calculation sheet for an air line. The pressure drop in each fitting and pipe
is calculated and the values in a new line are calculated based in the resulting pressure of the precedent line.
Rev. 09.08.2018)

6b. ods
Blower_Air_line.ods
(This file presents a calculation sheet for an air line. The pressure drop in each fitting and pipe
is calculated and the values in a new line are calculated based in the resulting pressure of the precedent line.
Rev. 09.08.2018)

15.
Combustion. Chemical reactions.xls
(Four examples of combustion analysis. Equations for stoichiometric combustion and combustion with excess air. Dew point temperature.
Rev. 04.12.2016)

21.
Compressor power and air discharge temperature.xls
(This file calculates the compressor power of an isentropic and of a real process.
Also, the exit temperature of an isentropic and of a real process is calculated.
Rev. 08.08.2016)

24.
Cooling tower. Application_ Treybal.xls
(This file is an application of the Merkel theory for cooling towers.
Some correction have been made, nomenclature reviewed and some literature added.
Still some explanations are required. Pending.
Re. 08.11.2014)

25.
Cooling tower. Merkel theory_Treybal.xls
(This file presents a resume of Merkel's theory for cooling tower, taken from Robert Treybal,
Operaciones de transferencia de masa.
Rev. 31.01.2014)

26.
Coupled water tanks. Stabilization time, oscilation amplitud.xls
(This file presents the solution of a system of ordinary differential equations, resulting from the water
movement between two tanks. The problem is solved using a finite differences method.
Rev. 08.08.2016)

27.
Dimensioning_compressed_air_installations_Atlas_Copco.xls
(This file uses an Atlas Copco reference for the dimensioning of a compressd air installation.
It includes compressor, aftercooler, receiver, dryer and pressure drop.
Rev. 15.04.2016 )

28.
Dryer with air impinging jets.xls
(Design of a strip dryer with air impinging jets. Air pressurized in a fan, heated in an heat exchanger and
impinging into a moving steel strip with a water film. Based on Prof. Martin Holger paper.
Rev. 13.01.2018)

32a.xls
Expansion_loop.xls
(This file derives the expansion loop equation an presents a calculation example. Spiraxsarco and Vitaulic references.
Rev. 01.08.2018)

32b.ods
Expansion_loop.ods
(This file derives the expansion loop equation an presents a calculation example. Spiraxsarco and Vitaulic references
Rev. 01.08.2018)

34.
Friction and singular pressure drop. Water and slurry.xls
(This file presents a routine for the calculation of pressure drops due to friction and singularities.
Several functions are included for the calculation of fittings and valves.
Rev. 08.08.2016)

47.
Gas. Viscosity of gases according Sutherland.xls
(Sutherland's formula. Dynamic viscosity of an ideal gas as a function of the temperature.
Valid for temperatures 0 < T < 555 K, with an error due to pressure less than 10%, below 3.45MPa
Rev. 27.02.2018)

48.
Gas. Wet air composition.xls
(Wet air for a given absolute air humidity. Wet air composition and molecular mass.
Rev. 11.08.2018)

49.
Gravitational adduction Water Hammer.xls
(Gravitational discharge of slurry with the option of three pipe sizes, under consideration of three flow rates..
Rev. 08.08.2016)

52.
Heat transfer. Annealing of a steel plate. Mills example 3.9.xls
(When steel plates are thinned by rolling, periodic reheating is required. A plain carbon steel plate of thickness "2*L", initially
at a temperature "tini", is to be reheated to a minimum temperature "tend" in a furnace mantained at "tfurn"
A one-term solution is implemented with VBA functions to replace the use of the Heisler charts.
Rev. 17.05.2018)

54.
Heat transfer. Cooling a pipe filled with water.xls
(This application calculates the cooling time of water in a carbon steel pipe.
The theory needed is deducted. Calculation examples of pipes with and without
insulation are presented. A comparison with the results presented in the
Mechanical Insulation Design Guide (NMIC) is included.
Rev. 26/04/2017)

58.
Heat transfer. Finned and bare pipes.xls.
(U-factor and heat flow of finned pipes. Fin efficiency according Schneider. Incropera and Schneider examples.
Rev. 27.02.2018)

71.
Heat transfer. Underwater pipe for effluent discharge.xls
(Underwater pipe for effluent discharge. Discharge temperature of effluent in the sea and
heat flow rate from the pipe into the sea. Exterior and interior convection coefficients.
Rev. 08.08.2016)

81.
Math. Runge-Kutta application to a tank problem.xls
(The Runge-Kutta Method is a commonly used numerical method for solving 1st-order ordinary
differential equations (ODEs) with a known initial condition. Jeff Munic application example to a tank .
Rev. 03.04.2018)

86.
Math. Tank Problem 1. Runge Kutta application.xls
(A brine solution of water is added at a given flow rate to a tank with pure water.
The volume in the tank is maintained at a constant volume with an overflow drain.
Required is the concentration change with time. Ref. Jeff Munic.
Rev. 17.04.2018)

87.
Math. Tank Problem 2. Runge Kutta application.xls
(A tank with a solution of caustic liquid is is inially at a given concentration and volume.
An upset occurs, and the supply flow rate and the inlet concentration drops down. Required
is the tank concentration. Ref. Jeff Munic.
Rev. 17.04.2018)

88.
Mc_Elvain_Cave_Durand_Bingham_fluids_HR_value.xls
(This file presents functions for the calculation of two slurry correction factors applied to the deposition velocity:
- Mc Elvain and Cave correction factor and Durand correction factor.
(This two function are a digitalization of the curves and therefore no equation is used).
- Also is presented a function for Weir - HR factor for estimating the head and the efficiency of slurries, based on the values for water
(Note. Weir, in later publications, is proposing a "HR-value" determination method that also requires the impeller diameter as input data.
Rev. 31.01.2014)

89.
Minimum_distance_between_pipes_with_ flanches.xls
(Flanches dimensions according ASME B16.5-2003.
Minimum distance between flanches and pipes: 30 mm.
Valid for pipes without insulation.
Distances to be verified if lateral movements or expansions could occur and
also if orifice plates or other elements are present.
Verify that there is not an occurrence of two flanges face to face.
Pipes according ASME B36.10M-1996.
Rev. 22.05.2016)

90.
Mollier_diagram.xls
(A Pressure-Enthalpy, mollier type diagram, built by means of Steamdat functions.
Rev. 04.04.2016)

91.
Moody diagram. Hagen Poiseuille, Colebrook and Churchill equations_Nikuradse data.xls
(For 0 < Re < 2300
Laminar region. Hagen - Poiseuille equation.
For 2300 =< Re =< 4000
Critical region. Churchill equation.
For 4000 < Re
Transition and turbulent regions. Colebrook equation.
There is not a theory describing the critical region.
Churchill equation describes relatively well this region, for smooth pipes with Rrel <= 0.01, giving conservative values,
when compared with Nikuradse experimental data.
Although Churchill equation describes also the transition and turbulente regions in accordance with Colebrook, this last equation is
used in these regions because its use is often required in certain design criteria.
Rev. 28.07.2016)

92.
Normal_to_real_flow_rate_and_FAD_flow rate.xls
(This file presents routines to transform Normal to Real flow rates and inversely,
Standard to real flow rates and inversely and FAD flow rates to real or Normal flow rates.
Rev. 21.04.2016)

93.
Orifice Plates.xls
(This file presents routines to calculate orifices plates with applications for air and water.
Also, Cameron eqautions for water are presented.
Rev. 26.09.2016)

94.a (xls)
>Pipe dimensions and friction factor.xls
(Visual Basic functions for pipe dimensions for carbon steel, stainless steel, HDPE PE100, HDPE PE80,
Fibre reinforced polyethylene, pipe friction factor for Darcy-Weisbach
equation and Manning's coefficient.
Rev. 06.07.2017)

94.b (ods)
Pipe dimensions and friction factor.ods
(LibreOffice Calc Basic functions for pipe dimensions for carbon steel, stainless steel, HDPE PE100, HDPE PE80,
pipe friction factor for Darcy-Weisbach equation and Manning's coefficient.
Rev. 06.07.2017)

95.a (xls)
>Pipe dimensions_CS_SS_HDPE100_HDPE80.xls
(Visual Basic functions for pipe dimensions for carbon steel, stainless steel, HDPE PE100, HDPE PE80,
Pipe friction factor for Darcy-Weisbach
equation and Manning's coefficient.
Rev. 04.07.2017)

95.b (ods)
>Pipe dimensions_CS_SS_HDPE100_HDPE80.ods
(LibreOffice Calc Basic functions for pipe dimensions for carbon steel, stainless steel, HDPE PE100, HDPE PE80,
Pipe friction factor for Darcy-Weisbach
equation and Manning's coefficient.
Rev. 04.07.2017)

101.
Pressure. Pressure loss in a steam pipe. Tabulated example.xls
(Pressure drop of a steam flow rate "m ton/h" in a carbon steel pipe with nominal diameter "dn", schedule "Sch"
and absolute rugosity "Rabs". The pipe is located at a hight above sea level "H m.a.s.l." The steam inlet pressure
is "pin_g bar (g)". Pipe lengths and fittings are shown in the calculation table.
Rev. 05.12.2017)

102.
Pressure. Pressure loss in a steam pipe. Tyler example.xls
(Pressure drop of a steam flow rate "m ton/h" in a carbon steel pipe with nominal diameter "dn", schedule "Sch"
and absolute rugosity "Rabs". Tyler Example with a pressure reducing valve.
Rev. 05.12.2017)

103.a (xls)
Pipes. Slope required for a pipe to avoid fluid accumulation.xls
(Slope of a pipe to avoid accumulation of fluid in case the pipe should be emptied.
To avoid the accumulation of fluid, one support shall be installed at a height lower than the other,
at a difference Dh [mm].
The tangent at the point of inflection (P) of the beam must become horizontal to get that
no fluid can remain stored.
Rev. 10.07.2017)

103.b (ods)
Pipe. Slope required for a pipe to avoid fluid accumulation.ods
(Slope of a pipe to avoid accumulation of fluid in case the pipe should be emptied.
To avoid the accumulation of fluid, one support shall be installed at a height lower than the other,
at a difference Dh [mm].
The tangent at the point of inflection (P) of the beam must become horizontal to get that
no fluid can remain stored.
Rev. 10.07.2017)

105.
Pneumatic transport in dilute phase. Rhodes example.xlsm
(Martin Rhodes, Introduction to particle technology. Example 8.1. Design calculation for dilute pneumatic transport.
Spreadsheet make use of some VBA functions. (Beta version)
Rev. 01.07.2018)

106.
Pressure. Maximum allowable pressure, ASME B31.3. Pipes A53, A106, API 5L (dn- Sch) at a given temperature.xls
(Maximum allowable pressure and temperature ratings for petroleum refinery piping and chemical plant piping systems
according ANSI/ASME B31.3 (2008) Process piping, materials grade B: A53, A106, API 5L, pipes with plane ends.
Allowable stresses from ASME B31.3, 2008, page 146) Maximum pressure calculated according Ec. 3a
Maximum temperature and pressure ratings of flanges conforming dimensions ASME B16.5 and materials specification ASTM A-105
Rev. 12.06.2015)

111.a (xls)
Psychrometric chart.xls
(Psychrometric charts: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
for heights above sea level of 0 m.a.s.l. and 5300 m.a.s.l.
Rev. 01.03.2018)

111.b (ods)
Psychrometric chart.ods
(Psychrometric charts: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
for heights above sea level of 0 m.a.s.l. and 5300 m.a.s.l.
Rev. 15.02.2017)

112.a (xls)
Psychrometric chart with process shown in diagram.xls
(Psychrometric charts: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
for heights above sea level of 0 m.a.s.l. and 5300 m.a.s.l.
Psychtometric functions for following input variable input groups:
1. tdb, f, H
2. tdw, twb, H
3. tdb, x, H
4. enthalpy, x, H
5. tdb, enthalpy, H
Rev. 01.03.2018)

112.b (ods)
Psychrometric chart with process shwon in diagram.ods
(Psychrometric charts: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
for heights above sea level of 0 m.a.s.l. and 5300 m.a.s.l.
Psychtometric functions for following input variable input groups:
1. tdb, f, H
2. tdw, twb, H
3. tdb, x, H
4. enthalpy, x, H
5. tdb, enthalpy, H
Rev. 22.05.2016)

113.a (xls)
Psychrometric functions_Deductions.xls
(Psychrometric functions: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
dew point temperature, specific volume and density, for heights above sea level til 5300 m.a.s.l.
Rev. 01.03.2018)

113.b (ods)
Psychrometric functions_Deductions.ods
(Psychrometric functions: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
dew point temperature, specific volume and density, for heights above sea level til 5300 m.a.s.l.
Rev. 15.02.2017)

114.a (xls)
Psychrometric functions_Resume.xls
(Psychrometric functions, only a resume: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
dew point temperature, specific volume and density, for heights above sea level til 5300 m.a.s.l.
Rev. 01.03.2018)

114.b (ods)
Psychrometric functions.ods
(Psychrometric functions, only a resume: Dry and wet bulb temperature, absolute humidity, relative humidity, enthalpy,
dew point temperature, specific volume and density, for heights above sea level til 5300 m.a.s.l.
Rev. 15.02.2017)

115.
Psychrometry. Heat recovery air handling unit (Ahu). By Ömer Faruk D.xls
(This spreadsheet calculates air flow and battery capacity for Air Handling Units.
The data used corresponds to a location in Turkey. In the example, data for the city of Bursa has been used.
You can change the data according to your city, in the Data page.
By Omer Faruk D., Makine Mühendisi , Mechanical Engineer
Rev. 01.10.2014)

117.
Pump. Detention time of a pump impulsion system.xls
(The routine calculates the time interval "t", from the de-energization of the pump, until the system comes to rest.
It is considered the inertia of the pump, motor and fluid and the friction between fluid and pipe.
An ascending pipe with constant slope is assumed. The friction factor is considered constant and with the value of
the steady state condition.
Rev. 29.07.2016)

119.
Pump. Heterogeneous_Slurries_Type_A_Warman.xls
(A froth pump calculation for heterogeneous slurry, according a Warman.
To calculate the pressure drop of a "Weir type A slurry", the system
is to be calculated as if the fluid were water.
The file presents a usual input data sheet a water pressure drop
calculation and finaly the calculation of the
pressure difference that in some cases has to be added to the
calculated pressure..
Ref. 14.06.2016)

126.
Pump. TDH , NPSH, Pump power.xls
(Pumping system between two water tanks.
Results are calculated in a spreadsheet and by means of "user defined Excel functions"
Ref. 06.07.2016)

130.
Sand_trap.xls
(Determination of the basic sand trap dimensions using a VBA function .
Rev. 31.01.2014)

131.
Settling_velocity_of_spherical_particles.xls
(Settling velocity of spherical particles as function of particle diameter,
solids density, liquid density and liquid absolute viscosity (VBA)
Rev. 05.06.2016)

138.a (xls)
Slurry. Settling velocity according JRI.xls
(JRI recommend three types of equations to calculated limit deposition velocities,
according the particle average size and pipe diameter.
Rev. 04.07.2017)

138.b (ods)
Slurry. Settling velocity according JRI.ods
(JRI recommend three types of equations to calculated limit deposition velocities,
according the particle average size and pipe diameter.
Rev. 04.07.2017)

140.
Steam. Applications using Steamdat_97.xls
(Steamdat function used to calculate a steam turbine stage and pressure
reducing valve "PRV" with desuperheating.
Rev. 06.04.2017)

141.a (xls)
Steam. Applications using Magnus Holmgren functions.xls
(Steam and water functions used to calculate a steam turbine stage and pressure
reducing valve "PRV" with desuperheating. The data used by the functions is included in the code.
Rev. 25.07.2017)

141.b (ods)
Steam. Applications using Magnus Holmgren functions.ods
(Steam and water functions used to calculate a steam turbine stage and pressure
reducing valve "PRV" with desuperheating. The data used by the functions is included in the code.
Rev. 25.07.2017)

148.
Tailings deposition_by Gordon McPhail_2008.pdf
(Prediction of the beach profile of high density thickened tailings from rheological
and small scale trial deposition data. By Gordon McPhail, from Metago Environmental
Engineers.
Rev. 13.10.2014)

153.
Viscosity of gases according Sutherland.ods
(Sutherland's formula. Dynamic viscosity of an ideal gas as a function of the temperature.
Valid for temperatures 0 < T < 555 K, with an error due to pressure less than 10%, below 3.45MPa
Rev. 15.02.2017)

154.
Viscosity of oils as a function of temperature.xls
(The viscosity of oils can be shown as straight lines in a Log-Nat Diagram
This concept is applied to the case of Rimula 15W-40 oil, where a paire
of points "viscosity - temperature" are known.
Rev. 26.09.2016)

158.
Water hammer. Method of characteristics. Equations, Boundary conditions and Initial conditions.xls
( The Method of Characteristics transform the system of Partial Differential Equations (PDE) in a system of Ordinary Differential Equations (ODE).
The system of ODE can be solved numerically by different mathematical methods . The Methods of Characteristics is used here.
The application is based on a Streeter example.
Rev. 05.06.2016)

159.
Water hammer. Method of characteristics. Example solved using Visual Basic and Finite Differences.xls
( An application for a simple case consisting in a reservoir, a horizontal pipe and a valve.
The water hammer problem is solved by Finite Differences and also using Visual Basic.
An application example is solved with input data from a Streeter example. Results of the solution by finite differences and V.B.
are in agreement with the results from Streeter example.
Rev. 05.06.2016)

163.
Water hammer. Water and slurry hammer.xls
( Water hammer examples:
Tyler. Water hammer in a carbon steel pipe
Pehmco: Water hammer in a HDPE pipe
Tsingua University: Slurry hammer in a HDPE PE80 pipe.
Rev. 28.05.2016)

Some Examples of Piping Tools Applications

Air and saturated water properties as a function of temperature [ºC]

The files are free of any virus.
As general good practice, it is recommended to check any type of file to be free of virus before using them.
The presented files contain Excel functions and it is known that Excel has some problems in recognizing
clean files. For this reason it is recommended that, once a file has been checked for viruses and has been
recognized as a clean file, to put it in a directory declared as a trusted directory.
To declare a directory or a file as a safe document, use the Excel trust Center.
This providence will avoid that under given circumstances Excel will erase the Visual Basic code.

DISCLAIMER

The author shall not be responsible for losses of any kind resulting from use of this program or of any
documentation and can in no way provide compensation for any losses sustained including but not limited
to any obligation, liability, right or remedy for tort nor any business expense, machine downtime or
damage caused to the user by any deficiency, defect or error in the program or in any such documentation
or any malfunction of the program or for any incidental or consequential losses, damages, or costs,
however caused.