Abstract
Fast transients in pumps and valves may induce significant variations in pressures and flow rates throughout pipelines that can even cause structural damages. This computational work deals with parameter estimation of a water hammer model, with focus on a transient friction coefficient and an empirical parameter related to the pipeline elasticity. The hyperbolic water hammer model was solved with a total variation diminishing version of the weighted average flux finite volume scheme. Simulated pressure and flow rate measurements taken near the inlet and the outlet of the pipeline were used for the solution of the parameter estimation problem with Markov Chain Monte Carlo methods. This work aimed at the reduction of the computational time of the inverse problem solution with two different strategies: (1) the application of a recent parallel computation version of the Metropolis–Hastings algorithm; and (2) the use of a machine learning metamodel obtained with the evolutionary neural network algorithm and the approximation error model approach. These two approaches are compared in terms of the parameter estimation accuracies and associated computational times.
Similar content being viewed by others
Data availability
Authors will provided data upon request.
Abbreviations
- a :
-
Celerity, m/s
- A :
-
Area of the pipe cross section, m2
- AEM:
-
Approximation error model
- B :
-
Vector of source terms
- c k :
-
Courant number
- D :
-
Pipe diameter, m
- e :
-
Pipe wall thickness, m
- e(P):
-
Vector of modeling errors
- E :
-
Young’s modulus of the pipe, Pa
- EvoNN:
-
Evolutionary neural network
- f :
-
Darcy–Weisbach friction coefficient
- f u :
-
Transient friction coefficient
- F :
-
Vector of fluxes
- G :
-
Pipe shear modulus, Pa
- g :
-
Gravity acceleration, m/s2
- H :
-
Head, m
- K :
-
Elastic modulus of the fluid, Pa
- K p,b :
-
Localized head loss coefficient
- k :
-
Iteration index of each parallel Markov chain
- L :
-
Pipe length, m
- MCMC:
-
Markov Chain Monte Carlo
- I :
-
Number of measurements
- N :
-
Number of waves in the Riemann problem
- N tr :
-
Number of training data for the neural networks
- n :
-
Time step
- P :
-
Vector of parameters
- Pr(Y):
-
Prior distribution
- Pr(P|Y):
-
Posterior distribution
- Pr(Y|P):
-
Likelihood function
- p :
-
Pressure, Pa
- PPGA:
-
Predator–Prey genetic algorithm
- q :
-
Mass flow rate, kg/s
- R :
-
Total number of cores for parallel computation
- R in :
-
Pipe internal radius
- R out :
-
Pipe external radius
- Re:
-
Reynolds number
- r :
-
Processor index
- S LS :
-
Least squares norm
- T(P) :
-
Vector of dependent variables obtained with the high-fidelity model
- T app (P) :
-
Vector of dependent variables obtained with the low-fidelity model
- t :
-
Time, s
- TVD:
-
Total variation diminishing
- u :
-
Velocity in the parallel direction of the flow, m/s
- U :
-
Vector of state variables
- W :
-
Covariance matrix
- WAF:
-
Weighted average flux
- x :
-
Axial position, m
- x b :
-
Position of the localized pressure loss, m
- Y :
-
Vector of measurements
- α :
-
Acceptance rate
- δ(.):
-
Dirac delta function
- ∆p b :
-
Localized pressure loss
- ε :
-
Measurements errors
- θ :
-
Angle between the flow direction and a horizontal reference plane
- ρ :
-
Density, kg/m3
- ϕ :
-
Proposal distribution
- Φ:
-
Limiter function
- Φ( t ) :
-
Diagonal matrix with the standard deviations of the proposal for each parameter
- υ P :
-
Poisson’s ratio
- Ψ:
-
Dimensionless parameter related to the pipe elasticity
- ω :
-
Pipe roughness
- η(P):
-
Total error
- Γ η P :
-
Covariance matrix of η and P
- 0:
-
Steady-state local value
- i :
-
Finite volume position
- k :
-
Index of the core
- max:
-
Maximum
- n :
-
Time step
- –:
-
Mean value
- *:
-
Candidate parameter
References
Abreu LA, Orlande HRB, Colaço MJ, Kaipio J, Kolehmainen V, Pacheco CC, Cotta RM (2018) Detection of contact failures with the Markov chain Monte Carlo method by using integral transformed measurement. Int J Therm Sci 132:486–497. https://doi.org/10.1016/j.ijthermalsci.2018.06.006
Axworthy DH, Ghidaoui MS, Mcinnis DA (2000) Extended thermodynamics derivation of energy dissipation in unsteady pipe flow. J Hydraul Eng 126(4):276–287. https://doi.org/10.1061/(ASCE)0733-9429(2000)126:4(276)
Bansal A, Barman A, Ghosh S, Chakraborti N (2012) Designing Cu-Zr glass using multiobjective genetic algorithm and evolutionary neural network metamodels-based classical molecular dynamics simulations. Mater Manuf Process 28(7):733–740. https://doi.org/10.1080/10426914.2013.763961
Beck JV, Arnold KJ (1997) Parameter estimation in engineering and science. Wiley Interscience, New York
Bhattacharya B, Kumar GRD, Agarwal A, Erkoç S, Singh A, Chakraborti N (2009) Analyzing Fe-Zn system using molecular dynamics, evolutionary neural nets and multi-objective genetic algorithms. Comput Mater Sci 46(4):821–827. https://doi.org/10.1016/j.commatsci.2009.04.023
Boulos PF, Karney BW, Wood DJ, Lingireddy S (2005) Hydraulic transient guideless for protecting water distribution systems. JAWRA 97(5):111–124. https://doi.org/10.1002/j.1551-8833.2005.tb10892.x
Brunone B, Golia UM, Greco M (1995) Effects of two-dimensionality on pipe transients modeling. J Hydraul Eng 121(12):906–912. https://doi.org/10.1061/(ASCE)0733-9429(1995)121:12(906)
Carvalho RC, Orlande HRB, Colaço MJ, Madeira IM (2018) Application of WAF–TVD scheme for the water hammer equations, Brazilian congress of thermal sciences and engineering, Águas de Lindóia, Brazil
Carvalho RC, Orlande HRB, Colaço MJ, Madeira IM (2022) Characterization of pipelines during flow transients. In: 10th international conference on inverse problems in engineering, Francavilla al Mare (Chieti), Italy
Chakraborti N (2022) Data-driven evolutionary modeling in materials technology. CRC Press, Boca Raton
Chaudhry MH (2014) Applied hydraulic transients, 3rd edn. Springer, New York
Coello C, Lamont GB, Van Veldhuizen D (2007) Evolutionary algorithms for solving multi-objective problems, vol 5. Springer, New York
Cui T, Fox C, Nicholls GK, O’Sullivan MJ (2019) Using parallel Markov Chain Monte Carlo to quantify uncertainties in geothermal reservoir calibration. Int J Uncertain Quantif 9(3):295–310. https://doi.org/10.1615/Int.J.UncertaintyQuantification.2019029282
Datta S, Pettersson F, Ganguly S, Saxen H, Chakraborti N (2008) Identification of factors governing mechanical properties of TRIP-aided steel using genetic algorithms and neural networks. Mater Manuf Process 23(2):130–137. https://doi.org/10.1080/10426910701774528
Deb K (2001) Multi-objective optimization using evolutionary algorithms. Wiley, Chichester
Gammerman D, Lopes HF (2006) Markov Chain Monte Carlo: stochastic simulation for Bayesian inference. Chapman & Hall/CRC, Boca Raton
Geneva N, Zabaras N (2019) Quantifying model form uncertainty in Reynolds-averaged turbulence models with Bayesian deep neural networks. J Comput Phys 385:125–147. https://doi.org/10.1016/j.jcp.2019.01.021
Ghidaoui MS, Zhao M, Mcinnis DA, Axworthy DH (2005) A review of water hammer theory and practice. Appl Mech Rev 58:49–76. https://doi.org/10.1115/1.1828050
Giri BK, Hakanen J, Miettinen J, Chakraborti N (2013) Genetic programming through bi-objective genetic algorithms with a study of a simulated moving bed process involving multiple objectives. Appl Soft Comput 13(5):2613–2623. https://doi.org/10.1016/j.asoc.2012.11.025
Govindan D, Chakraborty S, Chakraborti N (2010) Analyzing the fluid flow in continuous casting through evolutionary neural nets and multi-objective genetic algorithms. Steel Res Int 81(3):197–203. https://doi.org/10.1002/srin.200900128
Hastings WK (1970) Monte Carlo sampling methods using Markov Chains and their applications. Biometrika 57:97–109. https://doi.org/10.2307/2334940
Helle M, Petterson F, Chakraborti N, Saxen H (2006) Modelling noisy blast furnace data using genetic algorithms and neural networks. Steel Res Int 77(2):75–81. https://doi.org/10.1002/srin.200606357
Kaipio J, Somersallo P (2004) Computation and statistical methods for inverse problems. Springer, Berlin
Kumar H, Kumar S, Gnanasekaran N, Balaji C (2017) A Markov chain Monte Carlo-metropolis hastings approach for the simultaneous estimation of heat generation and heat transfer coefficient from a Teflon cylinder. Heat Transf Eng. https://doi.org/10.1080/01457632.2017.1305823
Lamien B, Orlande HRB, Eliçabe G (2017) Particle filter and approximation error model for state estimation in hyperthermia. J Heat Transf 139(1):12001. https://doi.org/10.1115/1.4034064
Lamien B, Orlande HRB (2013) Approximation error model to account for convective effects in liquids characterized by the line heat source probe, 4th Inverse Problems, Design and Optimization Symposium, Albi, France
Lee PM (2004) Bayesian statistics. Oxford University Press, London
Lüdecke HJ, Kothe B (2006) KSB Know-how Water Hammer. Volume 1. KSB
Madeira IM, Lucumi MAR, Orlande HRB (2021) State estimation problem for the detection of valve closure in gas pipelines. Inverse Probl Sci Eng 29(12):2186–2206. https://doi.org/10.1080/17415977.2021.1910682
Metropolis N, Rosenbluth A, Rosenbluth M, Teller A, Teller E (1953) Equation of state calculation by fast computing machines. J Chem Phys 21:1087–1092. https://doi.org/10.1063/1.1699114
Miettinen K (2012) Nonlinear multiobjective optimization, vol 12. Springer Science & Business Media, Berlin
Mijhaylov K, Rigopoulos S, Papadakis G (2021) Reconstruction of large-scale flow structures in a stirred tank from limited sensor data. AIChE J. https://doi.org/10.1002/aic.17348
Nissinen A (2011) Modelling errors in electrical impedance tomography, Dissertation in Forestry and Natural Sciences, University of Eastern Finland
Nissinen A, Heikkinen L, Kaipio J (2008) The Bayesian approximation error approach for electrical impedance tomography—experimental results. Meas Sci Technol 19(1):015501. https://doi.org/10.1088/0957-0233/19/1/015501
Nissinen A, Heikkinen L, Kolehmainen V, Kaipio J (2009) Compensation of errors due to discretization, domain truncation and unknown contact impedances in electrical impedance tomography. Meas Sci Technol 20:105504
Nissinen A, Kolehmainen V, Kaipio J (2011a) Compensation of modeling errors due to unknown boundary domain in electrical impedance tomography. IEEE Trans Med Imaging 30:231–242. https://doi.org/10.1088/0957-0233/20/10/105504
Nissinen A, Kolehmainen V, Kaipio J (2011b) Reconstruction of domain boundary and conductivity in electrical impedance tomography using the approximation error approach. Int J Uncertain Quantif 1:203–222. https://doi.org/10.1615/Int.J.UncertaintyQuantification.v1.i3.20
Orlande HRB, Colaço MJ, Dulikravich GS (2008) Approximation of the likelihood function in the Bayesian technique for the solution of inverse problems. Inverse Probl Sci Eng 16(6):677–692. https://doi.org/10.1080/17415970802231677
Orlande HRB, Fudym F, Maillet D, Cotta RM (2011) Thermal measurements and inverse techniques. CRC Press, Boca Raton
Orlande HRB, Dulikravich GS, Neumayer M, Watzenig D, Colaço MJ (2014) Accelerated Bayesian inference for the estimation of spatially varying heat flux in a heat conduction problem. Numer Heat Transf Part Appl 65(1):1–25. https://doi.org/10.1080/10407782.2013.812008
Ozisik MN, Orlande HRB (2021) Inverse heat transfer: fundamentals and applications, 2nd edn. CRC Press, Boca Raton
Ozisik MN, Orlande HRB, Colaço MJ, Cotta RM (2017) Finite difference methods in heat transfer, 2nd edn. CRC Press, Boca Raton
Peng GCY, Alber M, Tepole AB, Cannon WR, De S, Dura-Bernal S, Garikipati K, Karniadakis G, Lytton WW, Perdikaris P, Petzold L, Kuhl E (2021) Multiscale modeling meets machine learning: what can we learn? Achiev Comput Methods Eng 28:1017–1037. https://doi.org/10.48550/arXiv.1911.11958
Pettersson F, Chakraborti N, Saxen H (2007) A genetic algorithms based multi-objective neural net applied to noisy blast furnace data. Appl Soft Comput 7(1):387–397. https://doi.org/10.1016/j.asoc.2005.09.001
PyRVEA (2019) Research Group in Industrial Optimization. https://github.com/industrial-optimzation-group/pyRVEA/tree/master/pyrvea. Accessed 3 Sept 2019.
Raissi M, Perdikaris P, Karniadakis GE (2018) Physics-informed neural networks: a deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. J Comput Phys 378:686–707. https://doi.org/10.1016/j.jcp.2018.10.045
Rajak P, Tewary U, Das S, Bhattacharya B, Chakraborti N (2011) Phases in Zn-coated Fe analyzed trough an evolutionary meta-model and multi-objective genetic algorithms. Comput Mater Sci 50:2502–2516. https://doi.org/10.1016/j.commatsci.2011.03.034
Roberts GO, Gelman A, Gilks WR (1997) Weak convergence and optimal scaling of random walk Metropolis algorithms. Ann Appl Probab 7(1):110–120
Rogers CA, Kassab AJ, Divo EA, Ostrowski Z, Bialecki RA (2012) An inverse POD-RBF network approach to parameter estimation in mechanics. Inverse Probl Sci Eng 20(5):746–767. https://doi.org/10.1080/17415977.2012.693080
Rosenthal J (2011) Chapter 4: optimal proposal distribution and adaptive MCMC. In: Brooks S, Gelman A, Jones G, Meng X (eds) Handbook of Markov Chain Monte Carlo. CRC Press, Boca Raton
Szydlowski M (2002) Finite volume method for water hammer simulation, I international scientific and technical conference on technology, automation and control of wastewater and drinking water systems—TiASWiK'02, pp 159–165
Tan S, Fox C, Nicholls G (2006) Inverse problems—course notes for Physics 707, University of Auckland. https://math.unm.edu/~vageli/courses/Ma375/literature/
Toro EF (1997) Riemann solvers and numerical methods for fluid dynamics, 2nd edn. Springer, Berlin
Winkler R (2003) An introduction to Bayesian inference and decision. Probabilistic Publishing, Gainsville
Wylie EB (1997) Frictional effects in unsteady turbulent pipe flows. Appl Mech Rev 50(11):241–244. https://doi.org/10.1115/1.3101843
Yadav R, Balaji C, Venkasteshan SP (2019) Inverse estimation of number and location of discrete heaters in radiant furnaces using artificial neural networks and genetic algorithm. J Quant Spect Rad Transf 226:127–137. https://doi.org/10.1016/j.jqsrt.2018.12.031
Zhu Y, Zabaras N, Koutsourelakis PS, Perdikaris P (2019) Physics-constrained deep learning for high-dimensional surrogate modeling and uncertainty quantification without labeled data. J Comput Phys 394:56–81. https://doi.org/10.1016/j.jcp.2019.05.024
Acknowledgements
This study is mainly supported by Petrobras S.A., through project PT-112.01.13161. The support provided by CNPq, FAPERJ and CAPES (Finance Code 001, PQ1A, Cientista do Estado) is also greatly appreciated.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
There are no relationships of all authors with any people or organizations that could inappropriately influence (bias) this work.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Carvalho, R.C., Herzog, I.L., Orlande, H.R.B. et al. Parameter estimation with the Markov Chain Monte Carlo method aided by evolutionary neural networks in a water hammer model. Comp. Appl. Math. 42, 35 (2023). https://doi.org/10.1007/s40314-022-02162-0
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s40314-022-02162-0
Keywords
- Bayesian statistics
- Metropolis–Hastings algorithm
- Parallel computation
- Machine learning
- EvoNN
- Approximation error model