Loading [MathJax]/jax/output/SVG/jax.js

Methods for the localization of a leak in open water channels

  • Received: 01 November 2008 Revised: 01 January 2009
  • Primary: 35L40, 35B37; Secondary: 93C95.

  • In this paper, we present two methods for determining the position of a leak in an open water channel. The available measurements are the water level and the gate position at the upstream and downstream end of a channel reach. We assume that the size of the leak and the time it started are already estimated by a leak-detection method. Both of the proposed methods make use of a nonlinear Saint-Venant equation model of the channel where the leak is modelled as a lateral outflow. The first method makes use of a bank of N models corresponding to N possible positions of the leak along the channel. The estimated position of the leak is determined by the model which minimizes a quadratic cost function. The second method is based on the same principle except that it uses observers instead of pure models. The methods are tested on both real and simulated data from the Coleambally Channel 6 in Australia. It is further shown that the determination of the position of a leak is an inherently difficult problem.

    Citation: Nadia Bedjaoui, Erik Weyer, Georges Bastin. Methods for the localization of a leak in open water channels[J]. Networks and Heterogeneous Media, 2009, 4(2): 189-210. doi: 10.3934/nhm.2009.4.189

    Related Papers:

    [1] Huiting He, Mohamed Ousbika, Zakaria El Allali, Jiabin Zuo . Non-trivial solutions for a partial discrete Dirichlet nonlinear problem with p-Laplacian. Communications in Analysis and Mechanics, 2023, 15(4): 598-610. doi: 10.3934/cam.2023030
    [2] Dan Li, Yuhua Long . On periodic solutions of second-order partial difference equations involving p-Laplacian. Communications in Analysis and Mechanics, 2025, 17(1): 128-144. doi: 10.3934/cam.2025006
    [3] Zhen Wang, Luhan Sun . The Allen-Cahn equation with a time Caputo-Hadamard derivative: Mathematical and Numerical Analysis. Communications in Analysis and Mechanics, 2023, 15(4): 611-637. doi: 10.3934/cam.2023031
    [4] Leandro Tavares . Solutions for a class of problems driven by an anisotropic (p,q)-Laplacian type operator. Communications in Analysis and Mechanics, 2023, 15(3): 533-550. doi: 10.3934/cam.2023026
    [5] Tianfang Wang, Wen Zhang . Existence and concentration of homoclinic orbits for first order Hamiltonian systems. Communications in Analysis and Mechanics, 2024, 16(1): 121-146. doi: 10.3934/cam.2024006
    [6] Xiulan Wu, Yaxin Zhao, Xiaoxin Yang . On a singular parabolic p-Laplacian equation with logarithmic nonlinearity. Communications in Analysis and Mechanics, 2024, 16(3): 528-553. doi: 10.3934/cam.2024025
    [7] Wang Xiao, Xuehua Yang, Ziyi Zhou . Pointwise-in-time α-robust error estimate of the ADI difference scheme for three-dimensional fractional subdiffusion equations with variable coefficients. Communications in Analysis and Mechanics, 2024, 16(1): 53-70. doi: 10.3934/cam.2024003
    [8] Katica R. (Stevanović) Hedrih, Gradimir V. Milovanović . Elements of mathematical phenomenology and analogies of electrical and mechanical oscillators of the fractional type with finite number of degrees of freedom of oscillations: linear and nonlinear modes. Communications in Analysis and Mechanics, 2024, 16(4): 738-785. doi: 10.3934/cam.2024033
    [9] Caojie Li, Haixiang Zhang, Xuehua Yang . A new α-robust nonlinear numerical algorithm for the time fractional nonlinear KdV equation. Communications in Analysis and Mechanics, 2024, 16(1): 147-168. doi: 10.3934/cam.2024007
    [10] Mohamed Karim Hamdani, Lamine Mbarki, Mostafa Allaoui . A new class of multiple nonlocal problems with two parameters and variable-order fractional p()-Laplacian. Communications in Analysis and Mechanics, 2023, 15(3): 551-574. doi: 10.3934/cam.2023027
  • In this paper, we present two methods for determining the position of a leak in an open water channel. The available measurements are the water level and the gate position at the upstream and downstream end of a channel reach. We assume that the size of the leak and the time it started are already estimated by a leak-detection method. Both of the proposed methods make use of a nonlinear Saint-Venant equation model of the channel where the leak is modelled as a lateral outflow. The first method makes use of a bank of N models corresponding to N possible positions of the leak along the channel. The estimated position of the leak is determined by the model which minimizes a quadratic cost function. The second method is based on the same principle except that it uses observers instead of pure models. The methods are tested on both real and simulated data from the Coleambally Channel 6 in Australia. It is further shown that the determination of the position of a leak is an inherently difficult problem.


    In this article, we study the existence of infinitely many homoclinic solutions of the following fractional discrete p-Laplacian equations:

    (Δ1)spu(a)+V(a)|u(a)|p2u(a)=λf(a,u(a)),in Z, (1.1)

    where s(0,1) and p(1,) are fixed constants, V(a)R+, λ is a positive parameter, f(a,) is a continuous function for all aZ and (Δ1)sp is the fractional discrete p-Laplacian given by

    (Δ1)spu(a)=2bZ,ba|u(a)u(b)|p2(u(a)u(b))Ks,p(ab),in Z,

    where the discrete kernel Ks,p has the following property: There exist two constants 0<cs,pCs,p, such that

    {cs,p|d|1+psKs,p(d)Cs,p|d|1+ps,for alldZ{0};Ks,p(0)=0. (1.2)

    The fractional operator has received more attention recent decades because of its many applications in the real world. Many scholars have paid attention to this kind of problem, and have produced a lot of classical works, see for example [1,2,3]. As a classical fractional operator, the fractional Laplacian has wide applications in various fields such as optimization, population dynamics and so on. The fractional Laplacian on R can be defined for 0<s<1 and νC0(RN) as

    (Δ)sν(x)=CN,sP.V.RNν(x)ν(y)|xy|N+2sdy,xRN,

    where CN,s is a positive constant and P.V. denotes the Cauchy principle value. In various cases involving differential equations, the Laplace operator is replaced by either the fractional Laplace operator or other more general operators, and hence the existence results have been obtained by employing variational approaches, see for instance [4,5,6,7]. These alternative approaches have been found to offer improved descriptions of numerous phenomena observed in the natural world. Correspondingly, it is necessary to give some qualitative results by employing numerical analysis. The nonlocal feature of the fractional Laplacian is one of the important aspects to be considered in numerical methods, which makes it necessary to study the existence of solutions.

    Let ZH denote a grid of fixed size H>0 on R, i.e., ZH={Ha|aZ}. In [8], the definition of the fractional discrete Laplacian on ZH is given by

    (ΔH)sν(a)=bZ,ba(ν(a)ν(b))KHs(ab),

    where s(0,1), νs={ν:ZHR|ωZ|ν(ω)|(1+|ω|)1+2s<} and

    KHs(a)={4sΓ(1/2+s)π|Γ(s)|Γ(|a|s)H2sΓ(|a|+1+s),aZ{0},0,a=0.

    The above discrete kernel KHs has the following property: There exist two constants 0<csCs, such that for all aZ{0} there holds

    csH2s|a|1+2sKHs(a)CsH2s|a|1+2s.

    In [8], Ciaurri et al. also proved that if ν is bounded then lims1(ΔH)sν(a)=ΔHν(a), where ΔH is there discrete Laplacian on ZH, i.e.

    ΔHν(a)=1H2(ν(a+1)2ν(a)+ν(a1)).

    Moreover, under some suitable conditions and H0, the fractional Laplacian can be approximated by the fractional discrete Laplacian.

    Subsequently, let us give some existence results on the fractional difference equations. Xiang et al. [9] first investigated the fractional discrete Laplacian equations based on variational methods:

    {(Δ1)sν(a)+V(a)ν(a)=λf(a,ν(a)),foraZ,ν(a)0,as|a|, (1.3)

    where f(a,) is a continuous function for all aZ, λ>0, V(a)R+ and

    (Δ1)sν(a)=2bZ,ba(ν(a)ν(b))Ks(ab),in Z.

    Using the mountain pass theorem and Ekeland's variational principle under some suitable conditions, they obtained two homoclinic solutions for problem (1.3). It is evident that when p=2, the discrete fractional p-Laplace operator corresponds to the discrete fractional Laplace operator. After that, Ju et al. [10] studied the following fractional discrete p-Laplacian equations

    {(Δ1)spν(d)+V(d)|ν(d)|p2ν(d)=λa(d)|ν(d)|q2ν(d)+b(d)|ν(d)|r2ν(d),fordZ,ν(d)0,as|d|, (1.4)

    where appq, b, 1<q<p<r<, λ>0, V(a)R+, (Δ1)sp is the fractional discrete p-Laplacian. Under certain conditions, they employed the Nehari manifold method to achieve the existence of at least two homoclinic solutions for problem (1.4). In [11], Ju et al. investigated the existence of multiple solutions for the fractional discrete p-Laplacian equations with various nonlinear terms via different Clark's theorems. In a recent study conducted by Ju et al. in [12], it was demonstrated that using the fountain theorem and the dual fountain theorem under the same hypotheses, two separate sequences of homoclinic solutions were derived for the fractional discrete Kirchhoff-Schrödinger equations. Based on the findings from [8], it could be deduced that Eq.(1.3) can be reformulated as the renowned discrete version of the Schrödinger equation

    Δμ(ξ)+V(ξ)μ(ξ)=λf(ξ,μ(ξ)),in Z. (1.5)

    It is worth mentioning that in [13], Agarwal et al. first employed the variational methods to analyze Eq.(1.5). Here, We give some literature on the study of difference equations using the critical point theory, see [14,15,16].

    In particular, we observe that both the nonlinear terms in [11,12] have the following symmetry condition:

    (S) f(a,ν) is odd in ν.

    Therefore, in this paper, we consider the nonlinear term without condition (S), and study the existence of multiple homoclinic solutions of problem (1.1). For this, let us first recall if the solution ν of Eq.(1.1) satisfies ν(d)0 as |d|, then ν is called a homoclinic solution. Suppose that V(a) and f(a,u(a)) in problem (1.1) satisfy the following assumptions:

    (V) V1 and there is a constant V0(0,infaZV(a)]; (1 is defined in next section)

    (F) |f(a,u)|C(|u|p1+|u|t1) for any aZ and uR, where p<t< and C>0 is a constant.

    Set

    A:=lim infτ0+max|ζ|τaZF(a,ζ(a))τp,B:=lim supτ0+aZF(a,τ(a))τp,θ:=CsV0CbV1,

    where F(a,u)=u0f(a,η)dη and Cs, Cb will appear in next section. Here we give the main conclusion of our paper as follows.

    Theorem 1.1. Suppose that (V) and (F) are satisfied. Furthermore, the following inequality holds: A<θB. Then, for every λ(CbV1pB,CsV0pA), problem (1.1) possesses infinitely many nontrivial homoclinic solutions. In addition, their critical values and their -norms tend to zero.

    The rest of this article is arranged as follows: In Section 2, we introduce some definitions and give some preliminary results. In Section 3, we give the proof of Theorem 1.1. In Section 4, we give an example to demonstrate the main result.

    Here we illustrate some notations used in this paper:

    C, Cs,ω, Cs, Cb and C are diverse positive constants.

    denotes the embedding.

    denotes the strong convergence.

    First we give some basic definitions.

    Let 1ω, we give the definition of the space (ω,ω) as follows:

    ω:={{μ:ZR|aZ|μ(a)|ω<},if1ω<;{μ:ZR|supaZ|μ(a)|<},ifω=;
    μω:={(aZ|μ(a)|ω)1/ω,if1ω<;supaZ|μ(a)|,ifω=.

    Through the corresponding conclusions in [17], we know that ω is a Banach space. Moreover, ω1ω2 and μω2μω1 if 1ω1ω2.

    Next, we give the variational framework and some lemmas of this paper.

    The space (Q,Q) is defined by

    Q={σ:ZR|aZbZ|σ(a)σ(b)|pKs,p(ab)+dZV(d)|σ(d)|p<};
    σpQ=[σ]ps,p+dZV(d)|σ(d)|p=aZbZ|σ(a)σ(b)|pKs,p(ab)+dZV(d)|σ(d)|p.

    Lemma 2.1. (see [10, Lemma 2.1]) If ξω, then [ξ]s,ωCs,ωξω<.

    Lemma 2.2. (see [11, Lemma 2.2]) Under the hypothesis (V), (Q,Q) is a reflexive Banach space, and

    σ:=(aZV(a)|σ(a)|p)1/p.

    is an equivalent norm of Q.

    Through Lemma 2.2, we obtain that there exist 0CsCb such that

    CsμpμpQCbμp. (2.1)

    Lemma 2.3. Under the hypothesis (V), Qr is continuous for all pr.

    Proof. Using the above conclusions and (V), we can deduce that

    σrσp=(aZ|σ(a)|p)1pV1p0(aZV(a)|σ(a)|p)1p,σQ.

    As desired.

    Lemma 2.4. (see [10, Lemma 2.4]) If WQ is a compact subset, then for ι>0, a0N such that

    [|a|>a0V(a)|ξ(a)|p]1/p<ι,for eachξW.

    For all uQ, we define

    K(u)=D(u)λE(u)

    where

    D(u)=1paZbZ|u(a)u(b)|pKs,p(ab)+1pdZV(d)|u(d)|p=1pupQ

    and

    E(u)=dZF(d,u(d)).

    Clearly

    infQD(μ)=infQ1pμpQ=D(0)=0. (2.2)

    Lemma 2.5. (see [10, Lemma 2.5]) Under the hypothesis (V), then D(σ)C1(Q,R) with

    D(σ),ξ=aZbZ|σ(a)σ(b)|p2(σ(a)σ(b))(ξ(a)ξ(b))Ks,p(ab)+dZV(d)|σ(d)|p2σ(d)ξ(d),

    for all σ,ξQ.

    Lemma 2.6. (see [12, Lemma 2.6]) Under the hypotheses (V) and (F), then E(σ)C1(Q,R) with

    E(σ),ξ=dZf(d,σ(d))ξ(d)

    for all σ,ξQ.

    Combining Lemma 2.5 and Lemma 2.6, we know that K(σ)C1(Q,R).

    Lemma 2.7. Under the hypotheses (V) and (F), then for λ>0, every critical point of K is a homoclinic solution of problem (1.1).

    Proof. Assume σ be a critical point of K, we get for ξQ

    aZbZ|σ(a)σ(b)|p2(σ(a)σ(b))(ξ(a)ξ(b))Ks,p(ab)+aZV(a)|σ(a)|p2σ(a)ξ(a)=λaZf(a,σ(a))ξ(a). (2.3)

    For each aZ, we define edQ as follows:

    ed(a):={0,ifad;1,ifa=d.

    Taking ξ=ed in (2.3), we have

    2bZ,ba|u(a)u(b)|p2(u(a)u(b))Ks,p(ab)+V(a)|u(a)|p2u(a)=λf(a,σ(a)).

    So σ is a solution of problem (1.1). Moreover, by Lemma 2.3 and σQ, we know σ(a)0 as |a|. Thus, σ is a homoclinic solution of problem (1.1).

    In this section, we shall use the following Thoerem 3.1 to prove our main result. In fact, this theorem is a special version of Ricceri's variational principle [18, Lemma 2.5].

    Theorem 3.1. (see [19, Lemma 2.1]) Let Q be a reflexive Banach space, K(μ):=D(μ)+λE(μ) for each μQ, where D,EC1(Q,R), D is coercive, and λ is a real positive parameter. For every γ>infQD(μ), let

    η(γ):=infμD1((,γ))(supνD1((,γ))E(ν))E(μ)γD(μ),

    and

    ρ:=lim infγ(infQD(μ))+η(γ)

    If ρ<+, then for every λ(0,1ρ), the following conclusions holds only one:

    (a) there exists a global minimum of D which is a local minimum of K.

    (b) there exists a sequence {μm} of pairwise distinct critical points (local minima) of K, with limmD(μm)=infQD(μ), which converges to a global minimum of D.

    Remark 3.1. Obviously, ρ0. In addition, when ρ=0, we think that 1ρ=+.

    Proof of Theorem 1.1. Let us recall

    A=lim infτ0+max|ζ|τaZF(a,ζ(a))τp,B=lim supτ0+aZF(a,τ(a))τp,θ=CsV0CbV1,

    where F(a,u)=u0f(a,ω)dω. Fix λ(CbV1pB,CsV0pA) and set K,D,E as in Section 2. By Lemma 2.2, Lemma 2.5 and Lemma 2.6, we know Q be a reflexive Banach space and D,EC1(Q,R). Because of

    D(μ)=1paZbZ|μ(a)μ(b)|pKs,p(ab)+1pdZV(d)|μ(d)|p=1pμpQ+

    as μQ+, i.e. D is coercive. Now, we show that ρ<+. For this purpose, let {δn} be a positive sequence such that limnδn=0 and

    limnmax|ζ|δnaZF(a,ζ(a))δpn=A.

    Put

    γn:=CsV0pδpn,

    for all nN. Clearly, limnγn=0. For n>0 is big enough, by Lemma 2.2 and (2.1), we can derive that

    D1((,γn)){νQ:|ν(d)|δn,dZ}. (3.1)

    Since D(0)=E(0)=0, for each n large enough, by (3.1), we get

    η(γn)=infμD1((,γn))(supνD1((,γn))aZF(a,ν(a)))E(μ)γnD(μ)(supνD1((,γn))aZF(a,ν(a)))E(0)γnD(0)=supνD1((,γn))aZF(a,ν(a))γnmax|w|δnaZF(a,w(a))γn=pmax|w|δnaZF(a,w(a))CsV0δpn.

    Therefore, by (2.2), we acquire that

    ρ=lim infγ(infQD(μ))+η(γ)=lim infγn0+η(γn)limnη(γn)limnpmax|w|δnaZF(a,w(a))CsV0δpn=pACsV0. (3.2)

    From (3.2), we get

    λ(CbV1pB,CsV0pA)(0,1ρ).

    Next, we verify that 0 is not a local minimum of K. First, suppose that B=+. Choosing M such that M>CbV1pλ and let {hn} be a sequence of positive numbers, with limnhn=0, there exists n1N such that for all nn1

    aZF(a,hn)>Mhpn. (3.3)

    Therefore, let {ln} be a sequence in Q defined by

    ln(a):=hn,for allaZ.

    It is easy to infer that ln(a)Q0 as n. By (V), (2.1) and (3.3), we obtain

    K(ln)=D(ln)λE(ln)=1plnpQλaZF(a,ln)=1plnpQλaZF(a,hn)CbplnpλaZF(a,hn)=CbpaZV(a)|ln(a)|pλaZF(a,hn)=Cbp(aZV(a))|hn|pλaZF(a,hn)<CbpV1hpnλMhpn=(CbpV1λM)hpn.

    So, K(ln)<0=K(0) for each nn1 big enough. Next, suppose that B<+. Since λ>CbV1pB, there exists ε>0 such that ε<BCbV1pλ. Hence, also choosing {hn} be a sequence of positive numbers, with limnhn=0, there is n2N such that for all nn2

    aZF(a,hn)>(Bε)hpn. (3.4)

    Arguing as before and by choosing {ln} in Q as above, we get

    K(ln)=D(ln)λE(ln)=1plnpQλaZF(a,ln)=1plnpQλaZF(a,hn)CbplnpλaZF(a,hn)=CbpaZV(a)|ln(a)|pλaZF(a,hn)=Cbp(aZV(a))|hn|pλaZF(a,hn)<CbpV1hpnλ(Bε)hpn=(CbpV1λ(Bε))hpn.

    So, K(ln)<0=K(0) for each nn2 big enough. In general, 0 is not a local minimum of K. By Theorem 3.1, (a) is not valid, then we have a sequence {μn}Q of critical points of K such that

    limnD(μn)=limn1pμnQ=infQD(μ)=0

    and

    limnK(μn)=infQD(μ)=0.

    By Lemma 2.3, we gain

    μnCμnQ0

    as n. By Lemma 2.7, the problem (1.1) admits infinitely many nontrivial homoclinic solutions. In addition, their critical values and their -norms tend to zero. This completes the proof.

    Here, we give an example of a nonlinear term which can apply Theorem 1.1.

    Example 4.1. We define

    Ψ(n):=1333n,fornN+;
    Φ(n):=1333n1,fornN+;
    χ(n):=13(p+1)33n3,fornN+.

    Obviously, we know that Φ(n+1)<Ψ(n)<Φ(n) for all nN+ and limnΨ(n)=limnΦ(n)=0. Set

    f(a,u)=0,aZN+.

    And for each aN+, let f(a,) is a nonnegative continuous function such that

    f(a,u)=0,uR(Ψ(a),Φ(a))andΦ(a)Ψ(a)f(a,η)dη=χ(a).

    There are many nonlinear terms that satisfy the above conditions. Here we give one of them as an example.

    f(a,u)=nN+χ(n)2((u+Φ(n)+Ψ(n)2)2Φ2(n)Ψ2(n))e{n}×[Ψ(n),Φ(n)](a,u)

    where eM×N is the indicator function on M×N. Then

    A=lim infτ0+max|ζ|τaZF(a,ζ(a))τplimnmax|ζ|Ψ(n)aZF(a,ζ(a))Ψp(n)=limna=n+1F(a,ζ(a))Ψp(n)=limna=n+1χ(a)Ψp(n)limn3χ(n+1)Ψp(n)=0

    and

    B=lim supτ0+aZF(a,ζ(a))τplimnaZF(a,ζ(a))Φp(n)limnF(n,ζ(n))Φp(n)limnχ(n)Φp(n)=+.

    Now it is easy to see that all the assumptions of Theorem 1.1 are satisfied, hence the corresponding conclusion can be delivered by Theorem 1.1.

    B. Zhang was supported by the National Natural Science Foundation of China (No. 12171152), the Shandong Provincial Natural Science Foundation, PR China (No. ZR2023MA090) and the Cultivation Project of Young and Innovative Talents in Universities of Shandong Province.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

  • This article has been cited by:

    1. Meiqiang Feng, Nontrivial p-convex solutions to singular p-Monge-Ampère problems: Existence, Multiplicity and Nonexistence, 2024, 16, 2836-3310, 71, 10.3934/cam.2024004
    2. Caihong Gu, Yanbin Tang, Well-posedness of Cauchy problem of fractional drift diffusion system in non-critical spaces with power-law nonlinearity, 2024, 13, 2191-950X, 10.1515/anona-2024-0023
    3. Eleonora Amoroso, Angela Sciammetta, Patrick Winkert, Anisotropic (p,q)-Laplacian problems with superlinear nonlinearities, 2024, 16, 2836-3310, 1, 10.3934/cam.2024001
  • Reader Comments
  • © 2009 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4157) PDF downloads(54) Cited by(21)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog