Application of Newton–Raphson Method for Computing the Final Air–Water Interface Location in a Pipe Water Filling

The estimation of thermodynamic behavior during filling processes with entrapped air in water pipelines is a complex task as it requires solving a system of algebraic-differential equations. A lot of different numerical methods have been used for this purpose in literature including the rigid water...

Full description

Autores:
Bonilla-Correa, Dalia M.
Coronado-Hernández, Óscar E.
Fuertes-Miquel, Vicente S.
Besharat, Mohsen
Ramos, Helena M.
Tipo de recurso:
Fecha de publicación:
2023
Institución:
Universidad Tecnológica de Bolívar
Repositorio:
Repositorio Institucional UTB
Idioma:
eng
OAI Identifier:
oai:repositorio.utb.edu.co:20.500.12585/12271
Acceso en línea:
https://hdl.handle.net/20.500.12585/12271
https://doi.org/10.3390/w15071304
Palabra clave:
Air; Geysers;
Emptying
LEMB
Rights
openAccess
License
http://creativecommons.org/licenses/by-nc-nd/4.0/
Description
Summary:The estimation of thermodynamic behavior during filling processes with entrapped air in water pipelines is a complex task as it requires solving a system of algebraic-differential equations. A lot of different numerical methods have been used for this purpose in literature including the rigid water column (RWC) model. The main advantage of the RWC model is its acceptable accuracy with very low computational load. In that context, this research presents the computation of critical points of the physical equations that describe the phenomenon. These points provide information about the final position of the air–water interface. The Newton–Raphson method was then applied to obtain a unique equation that can be used by engineers to directly compute variables such as air pocket pressure and water column length at the end of the hydraulic event. A case study was analyzed to compare the results of the mathematical model with the obtained equation for computing critical points. Both methods provided the same values for the water column length at the end of the hydraulic event. A sensitivity analysis was conducted to identify dependent and non-dependent parameters for evaluating the critical points. The proposed formulation was validated through an experimental set of data. © 2023 by the authors.