American Institute of Mathematical Sciences

November  2020, 16(6): 3011-3033. doi: 10.3934/jimo.2019091

Numerical solution to an inverse problem on a determination of places and capacities of sources in the hyperbolic systems

 1 Baku State University, Z. Khalilov, 23, Baku, Azerbaijan 2 Institute of Control Systems of Azerbaijan, National Academy of Sciences, B. Vahabzade, 9, AZ1141 Baku, Azerbaijan

* Corresponding author

Received  September 2018 Revised  April 2019 Published  July 2019

In the work we investigate the numerical solution to a class of inverse problems with respect to the system of differential equations of hyperbolic type. The specialties of considered problems are: 1) the impulse impacts are present in the system and it is necessary to determine the capacities and the place of their location; 2) the differential equations of the system are only related to boundary values, and arbitrarily; 3) because of the long duration of the object functioning, the exact values of the initial conditions are not known, but a set of possible values is given. The inverse problem under consideration is reduced to the problem of parametric optimal control without initial conditions with non-separated boundary conditions. For the solution it is proposed to use first-order optimization methods. The results of numerical experiments are given on the example of the inverse problem of fluid transportation in the pipeline networks of complex structure. The problem is to determine the locations and the volume of leakage of raw materials based on the results of additional observations of the state of the transportation process at internal points or at the ends of sections of the pipeline network.

Citation: Yegana Ashrafova, Kamil Aida-Zade. Numerical solution to an inverse problem on a determination of places and capacities of sources in the hyperbolic systems. Journal of Industrial & Management Optimization, 2020, 16 (6) : 3011-3033. doi: 10.3934/jimo.2019091
References:

show all references

References:
The scheme of pipeline network in the problem 1
Plots of the exact ($q_{loss} = q^{(1,2)^\ast }(t) = 50 - 10e^{- 0.0003 t} {m^3} /{hour}$ and obtained leak volume functions $q^{(1,2)^\ast }(t)$ for five various initial approximations $(\xi ^{(1,2)},\,\,q^{(1,2)}(t))_0^i ,i = 1,..,5.$ for the iteration process (10) in problem 1
Plots of exact ($q^{loss*}$) and computed leak volume functions ($q^*(t)$)for various suspected segments
Results of the solution to the problem 1 regarding leakage in the segment (1, 2) in the presence of noise for various initial states $(\xi ^{(1,2)},\,\,q^{(1,2)}(t))_0^i ,i = 1,..,5.$
 $\xi _0^{(1,2)}$ 60 20 90 10 45,685 $\eta$ $q_0^{(1,2)} (t)$ $90 - 10e^{3\alpha t}$ $20 -10e^{3\alpha t}$ $30-10e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $\tilde {\Phi _0}$ 76.104 16.704 11.664 42.48 57.636 $\tilde {\Phi }$ 5.73$\cdot 10^{ - 7}$ 1.26$\cdot 10^{ - 7}$ 3.19$\cdot 10^{ - 6}$ 1.85$\cdot 10^{ - 6}$ 7.43$\cdot 10^{ - 7}$ 0% $\tilde {\xi }^{(1,2)}$ 30.003 29.998 30.008 29.994 29.998 $N_{iter}$ 6 5 16 14 8 $\delta \xi ^{(1,2)}$ 0.00009 0.00006 0.0003 0.0002 0.00006 $\delta q^{(1,2)}$ 0.0003 0. 0006 0.0003 0.0002 0.0003 $\tilde {\Phi _0}$ 76.352 16.837 11.683 42.148 57.732 $\tilde {\Phi }$ 0.023 0.014 0.024 0.017 0.020 0.5% $\tilde {\xi }^{(1,2)}$ 29.841 30.332 30.068 29.654 29.796 $N_{iter}$ 6 5 14 12 7 $\delta \xi ^{(1,2)}$ –0.005 0.011 0.002 –0.011 –0.006 $\delta q^{(1,2)}$ 0.043 0.030 0.049 0.041 0.042 $\tilde{\Phi _0}$ 77.119 16.924 11.832 43.744 57.413 $\tilde {\Phi }$ 0.067 0.071 0.073 0.062 0.065 1% $\tilde {\xi }^{(1,2)}$ 28.527 29.392 29.923 30.597 29.839 $N_{iter}$ 6 5 14 12 7 $\delta \xi ^{(1,2)}$ –0.049 –0.020 –0.002 0.020 –0.005 $\delta q^{(1,2)}$ 0.109 0.092 0.107 0.094 0.093
 $\xi _0^{(1,2)}$ 60 20 90 10 45,685 $\eta$ $q_0^{(1,2)} (t)$ $90 - 10e^{3\alpha t}$ $20 -10e^{3\alpha t}$ $30-10e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $\tilde {\Phi _0}$ 76.104 16.704 11.664 42.48 57.636 $\tilde {\Phi }$ 5.73$\cdot 10^{ - 7}$ 1.26$\cdot 10^{ - 7}$ 3.19$\cdot 10^{ - 6}$ 1.85$\cdot 10^{ - 6}$ 7.43$\cdot 10^{ - 7}$ 0% $\tilde {\xi }^{(1,2)}$ 30.003 29.998 30.008 29.994 29.998 $N_{iter}$ 6 5 16 14 8 $\delta \xi ^{(1,2)}$ 0.00009 0.00006 0.0003 0.0002 0.00006 $\delta q^{(1,2)}$ 0.0003 0. 0006 0.0003 0.0002 0.0003 $\tilde {\Phi _0}$ 76.352 16.837 11.683 42.148 57.732 $\tilde {\Phi }$ 0.023 0.014 0.024 0.017 0.020 0.5% $\tilde {\xi }^{(1,2)}$ 29.841 30.332 30.068 29.654 29.796 $N_{iter}$ 6 5 14 12 7 $\delta \xi ^{(1,2)}$ –0.005 0.011 0.002 –0.011 –0.006 $\delta q^{(1,2)}$ 0.043 0.030 0.049 0.041 0.042 $\tilde{\Phi _0}$ 77.119 16.924 11.832 43.744 57.413 $\tilde {\Phi }$ 0.067 0.071 0.073 0.062 0.065 1% $\tilde {\xi }^{(1,2)}$ 28.527 29.392 29.923 30.597 29.839 $N_{iter}$ 6 5 14 12 7 $\delta \xi ^{(1,2)}$ –0.049 –0.020 –0.002 0.020 –0.005 $\delta q^{(1,2)}$ 0.109 0.092 0.107 0.094 0.093
The results of the solution to the problem 1 regarding leakage in the section (1, 2) for different $\nu$ and for different initial states $(\xi ^{(1,2)},\,\,q^{(1,2)}(t))_0^i ,i = 1,..,5.$
 $\xi _0^{(1,2)}$ 60 20 90 10 45.685 $\nu$ $q_0^{(1,2)} (t)$ $90-10e^{3\alpha t}$ $20-10e^{3\alpha t}$ $30-10e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $\tilde{\Phi _0}$ 155.248 22.323 22.323 84.303 118.809 $\tilde {\Phi}$ 8.2$\cdot 10^{ - 4}$ 2.7$\cdot 10^{ - 4}$ 1.0$\cdot 10^{ - 3}$ 2.2$\cdot 10^{ - 4}$ 7.4$\cdot 10^{ - 5}$ 3 $\tilde {\xi }^{(1,2)}$ 30.145 29.996 30.190 30.022 30.029 $\delta \xi ^{(1,2)}$ 0.004 0.0001 0.006 0.0007 0.0005 $\delta q^{(1,2)}$ 0.0005 0.0004 0.0001 0.0009 0.0005 $\tilde{\Phi _0}$ 155.390 33.608 22.443 84.086 118.854 $\tilde {\Phi }$ 1.6$\cdot 10^{ - 4}$ 3.9$\cdot 10^{ - 4}$ 4.3$\cdot 10^{ - 4}$ 2.8$\cdot 10^{ - 4}$ 6.4$\cdot 10^{ - 5}$ 5 $\tilde {\xi }^{(1,2)}$ 29.963 30.093 30.100 29.985 29.995 $\delta \xi ^{(1,2)}$ 0.001 0.003 0.003 0.0005 0.0001 $\delta q^{(1,2)}$ 0.0006 0.0003 0.0005 0.0004 0.0002
 $\xi _0^{(1,2)}$ 60 20 90 10 45.685 $\nu$ $q_0^{(1,2)} (t)$ $90-10e^{3\alpha t}$ $20-10e^{3\alpha t}$ $30-10e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $66 + 20e^{3\alpha t}$ $\tilde{\Phi _0}$ 155.248 22.323 22.323 84.303 118.809 $\tilde {\Phi}$ 8.2$\cdot 10^{ - 4}$ 2.7$\cdot 10^{ - 4}$ 1.0$\cdot 10^{ - 3}$ 2.2$\cdot 10^{ - 4}$ 7.4$\cdot 10^{ - 5}$ 3 $\tilde {\xi }^{(1,2)}$ 30.145 29.996 30.190 30.022 30.029 $\delta \xi ^{(1,2)}$ 0.004 0.0001 0.006 0.0007 0.0005 $\delta q^{(1,2)}$ 0.0005 0.0004 0.0001 0.0009 0.0005 $\tilde{\Phi _0}$ 155.390 33.608 22.443 84.086 118.854 $\tilde {\Phi }$ 1.6$\cdot 10^{ - 4}$ 3.9$\cdot 10^{ - 4}$ 4.3$\cdot 10^{ - 4}$ 2.8$\cdot 10^{ - 4}$ 6.4$\cdot 10^{ - 5}$ 5 $\tilde {\xi }^{(1,2)}$ 29.963 30.093 30.100 29.985 29.995 $\delta \xi ^{(1,2)}$ 0.001 0.003 0.003 0.0005 0.0001 $\delta q^{(1,2)}$ 0.0006 0.0003 0.0005 0.0004 0.0002
The results of the solution to the problem 1 regarding leakage in the section (5, 6) for different initial states
 $\xi _0^{(5,6)}$ $q_0^{(5,6)} (t)$ $\tilde {\xi }^{(5,6)}$ $\tilde{\Phi_0}$ $\tilde {\Phi}$ $N_{iter}$ 25 $5+10e^{3\alpha t}$ 15.219 2.556 3.55 $\cdot 10^{ - 5}$ 165 10 $30 +10e^{3\alpha t}$ 14.781 44.615 3.55 $\cdot 10^{ - 5}$ 204 22.56 $15+10e^{3\alpha t}$ 15.222 3.467 3.59 $\cdot 10^{ - 5}$ 197
 $\xi _0^{(5,6)}$ $q_0^{(5,6)} (t)$ $\tilde {\xi }^{(5,6)}$ $\tilde{\Phi_0}$ $\tilde {\Phi}$ $N_{iter}$ 25 $5+10e^{3\alpha t}$ 15.219 2.556 3.55 $\cdot 10^{ - 5}$ 165 10 $30 +10e^{3\alpha t}$ 14.781 44.615 3.55 $\cdot 10^{ - 5}$ 204 22.56 $15+10e^{3\alpha t}$ 15.222 3.467 3.59 $\cdot 10^{ - 5}$ 197
Results obtained by solving Problem 2
 $\xi ^{0}$ $q^{0} (t)$ Segment $\tilde{\xi }^{ks}$ $\Im _{0}$ $\tilde{\Im }$ Number of indices $(k,s)$ iterations (1, 2) 5.991 43.272 39.614 6 (5, 2) 2.279 139.536 60.386 58 25 $15+10e^{-0.0003t}$ (3, 2) 29.974 10.188 1.22$\cdot 10^{-5}$ 12 (5, 4) 2.244 54.185 37.958 79 (5, 6) 2.995 72.396 20.621 185 (1, 2) 5.986 46.836 39.614 5 (5, 2) 2.305 84.132 60.408 5 5 $5+10e^{-0.0003t}$ (3, 2) 29.935 32.148 2.44$\cdot 10^{-5}$ 12 (5, 4) 1.898 38.628 37.728 150 (5, 6) 2.947 21.786 20.621 72
 $\xi ^{0}$ $q^{0} (t)$ Segment $\tilde{\xi }^{ks}$ $\Im _{0}$ $\tilde{\Im }$ Number of indices $(k,s)$ iterations (1, 2) 5.991 43.272 39.614 6 (5, 2) 2.279 139.536 60.386 58 25 $15+10e^{-0.0003t}$ (3, 2) 29.974 10.188 1.22$\cdot 10^{-5}$ 12 (5, 4) 2.244 54.185 37.958 79 (5, 6) 2.995 72.396 20.621 185 (1, 2) 5.986 46.836 39.614 5 (5, 2) 2.305 84.132 60.408 5 5 $5+10e^{-0.0003t}$ (3, 2) 29.935 32.148 2.44$\cdot 10^{-5}$ 12 (5, 4) 1.898 38.628 37.728 150 (5, 6) 2.947 21.786 20.621 72
 [1] Pengyu Chen. Non-autonomous stochastic evolution equations with nonlinear noise and nonlocal conditions governed by noncompact evolution families. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020383 [2] Mengni Li. Global regularity for a class of Monge-Ampère type equations with nonzero boundary conditions. Communications on Pure & Applied Analysis, 2021, 20 (1) : 301-317. doi: 10.3934/cpaa.2020267 [3] Gunther Uhlmann, Jian Zhai. Inverse problems for nonlinear hyperbolic equations. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 455-469. doi: 10.3934/dcds.2020380 [4] Yangrong Li, Shuang Yang, Qiangheng Zhang. Odd random attractors for stochastic non-autonomous Kuramoto-Sivashinsky equations without dissipation. Electronic Research Archive, 2020, 28 (4) : 1529-1544. doi: 10.3934/era.2020080 [5] Lin Shi, Xuemin Wang, Dingshi Li. Limiting behavior of non-autonomous stochastic reaction-diffusion equations with colored noise on unbounded thin domains. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5367-5386. doi: 10.3934/cpaa.2020242 [6] Manil T. Mohan. First order necessary conditions of optimality for the two dimensional tidal dynamics system. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020045 [7] Nguyen Huy Tuan. On an initial and final value problem for fractional nonclassical diffusion equations of Kirchhoff type. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020354 [8] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020449 [9] Héctor Barge. Čech cohomology, homoclinic trajectories and robustness of non-saddle sets. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020381 [10] Ying Lin, Qi Ye. Support vector machine classifiers by non-Euclidean margins. Mathematical Foundations of Computing, 2020, 3 (4) : 279-300. doi: 10.3934/mfc.2020018 [11] Sergey Rashkovskiy. Hamilton-Jacobi theory for Hamiltonian and non-Hamiltonian systems. Journal of Geometric Mechanics, 2020, 12 (4) : 563-583. doi: 10.3934/jgm.2020024 [12] Noufel Frikha, Valentin Konakov, Stéphane Menozzi. Well-posedness of some non-linear stable driven SDEs. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 849-898. doi: 10.3934/dcds.2020302 [13] Wenqiang Zhao, Yijin Zhang. High-order Wong-Zakai approximations for non-autonomous stochastic $p$-Laplacian equations on $\mathbb{R}^N$. Communications on Pure & Applied Analysis, 2021, 20 (1) : 243-280. doi: 10.3934/cpaa.2020265 [14] Dan Zhu, Rosemary A. Renaut, Hongwei Li, Tianyou Liu. Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020076 [15] Bernard Bonnard, Jérémy Rouot. Geometric optimal techniques to control the muscular force response to functional electrical stimulation using a non-isometric force-fatigue model. Journal of Geometric Mechanics, 2020  doi: 10.3934/jgm.2020032 [16] João Marcos do Ó, Bruno Ribeiro, Bernhard Ruf. Hamiltonian elliptic systems in dimension two with arbitrary and double exponential growth conditions. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 277-296. doi: 10.3934/dcds.2020138 [17] Karoline Disser. Global existence and uniqueness for a volume-surface reaction-nonlinear-diffusion system. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 321-330. doi: 10.3934/dcdss.2020326 [18] H. M. Srivastava, H. I. Abdel-Gawad, Khaled Mohammed Saad. Oscillatory states and patterns formation in a two-cell cubic autocatalytic reaction-diffusion model subjected to the Dirichlet conditions. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020433 [19] Marco Ghimenti, Anna Maria Micheletti. Compactness results for linearly perturbed Yamabe problem on manifolds with boundary. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020453 [20] Antoine Benoit. Weak well-posedness of hyperbolic boundary value problems in a strip: when instabilities do not reflect the geometry. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5475-5486. doi: 10.3934/cpaa.2020248

2019 Impact Factor: 1.366