Next Article in Journal
Stability and Threshold Dynamics in a Seasonal Mathematical Model for Measles Outbreaks with Double-Dose Vaccination
Next Article in Special Issue
Generalized Halanay Inequalities and Relative Application to Time-Delay Dynamical Systems
Previous Article in Journal
Research on Relative Threshold of Abnormal Travel in Subway Based on Bilateral Curve Fitting
Previous Article in Special Issue
Hybrid Impulsive Pinning Control for Mean Square Synchronization of Uncertain Multi-Link Complex Networks with Stochastic Characteristics and Hybrid Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Positive Periodic Solution for Pipe/Tank Flow Configurations with Friction

School of Mathematics and Statistics, Huaiyin Normal University, Huai’an 223300, China
*
Author to whom correspondence should be addressed.
Submission received: 7 March 2023 / Revised: 3 April 2023 / Accepted: 7 April 2023 / Published: 9 April 2023

Abstract

:
In this article, we study a periodic boundary value problem related to valveless pumping. The valveless pumping is described by the unidirectional flow of liquid in a system. We establish some conditions for globally asymptotic stability and the existence of a positive periodic solution to the considered equation. Finally, a numerical example shows that the theoretical results in this paper are feasible.

1. Introduction

In 2006, Propst [1] presented a detailed explanation of the pumping effect for flow configurations of one to three rigid tanks connected by rigid pipes. A pump can be considered as a periodic differential system. Propst found the existence of periodic solutions to the corresponding differential equation for a system with two or three tanks. For the simplest structure of one pipe and one tank, Propst proved the existence result of a periodic solution in particular conditions. In general, the structure of one pipe and one tank can be expressed by the following boundary value problem
x + b x = 1 x ( γ ( t ) c ( u ) 2 ) d ,
x ( 0 ) = x ( T ) , x ( 0 ) = x ( T ) ,
where  b 0 , c > 1 , d > 0  and  γ ( t )  is a continuous periodic function. The specific meaning of the involved coefficients in (1) can be found in [1,2]. Equation (1) is a singular equation. Hakl, Torres and Zamora [3,4] studied the following singular equation:
y + u ( y ) y + v ( t , y ) = w ( t , y ) ,
where both u and v may have singularity at the origin. Since the study of singular equations is more complex and difficult than the study of regular equations, we attempt to find a more simple and effective method to study Equation (1). In fact, Equation (1) can be transformed into an equivalent regular equation through appropriate variable substitution. Let  x = y μ  with  μ = 1 c + 1 , and then we transform boundary value problem (1) and (2) to the following boundary value problem
y + a y + δ ( t ) y β ρ ( t ) y α = 0 ,
y ( 0 ) = y ( T ) , y ( 0 ) = y ( T ) ,
where  ρ ( t ) = γ ( t ) μ , δ ( t ) = d μ , α = 1 2 μ , β = 1 μ , γ ( t )  is a continuous periodic function. Since Equation (3) is not a singular equation, many classic research methods, such as the upper and lower solution method and topological degree theory, can be used for studying Equation (3). In [2], Cid, Propst and Tvrdý obtained the following theorems.
Theorem 1. 
Assume that  b > 0 , c > 1 , d > 0 , γ > 0  and that  γ ( t )  is continuous periodic on  R . Then, for BVP (1) and (2), there exists a positive solution if the following condition is satisfied:
( c + 1 ) d 2 4 γ < π 2 T 2 + b 2 4 ,
where  γ = min t [ 0 , T ] γ ( t ) .
Theorem 2. 
Assume that  b > 0 , c > 1 , d > 0 , γ > 0  and that  γ ( t )  is continuous periodic on  R . Then, for BVP (1) and (2), there exists at least one stable positive solution if the following conditions are satisfied:
d 2 γ < π 2 T 2 + b 2 4
and
( b 1 ) e < b e ,
where  γ = min t [ 0 , T ] γ ( t ) , γ = max t [ 0 , T ] γ ( t ) .
After that, Dorociaková et al. [5] considered the existence and stability of a periodic solution for valveless pumping using some fixed point theorems and the inequality technique. However, when using the lower and upper functions method, constructing proper lower and upper functions for BVP (1) and (2) is difficult in [2]. Furthermore, in order to obtain the existence and stability of the periodic solution for BVP (1) and (2), Lemmas 2.2, 2.4 and 2.5 in [2] were used, which greatly increases the difficulty of proof. In [5], sufficient conditions for the existence and exponential stability of a positive periodic solution to BVP (3) and (4) are complicated and are difficult to verify. In Remark 3 of this paper, we find that conditions (39)–(42) are sharp. Therefore, it is important to achieve simple and easily verifiable conditions.
Encouraged by the above work, we study periodic boundary value problem (3) and (4) by using coincidence degree theory. Topological degree theory is often used to study periodic solutions of differential equations; see, for example, quaternion-valued inertial memristor-based neural networks [6]; higher-order delay differential equations [7,8,9]; and fractional multi-point boundary value problems [10]. In this paper, we use the coincidence degree theory in the study of periodic solutions of a strong role combined with the appropriate variable substitution and mathematical analysis to study the periodic solution of pipe/tank flow configurations with friction.
With regard to some other important results about differential equations and systems, see, e.g., [11,12,13,14,15,16,17]. In particular, we transform Equation (3) into a equivalent two-dimensional system through appropriate variable substitution and can then easily study a periodic solution for the above two-dimensional system. Furthermore, we obtain that the periodic solution is stable through using the Lyapunov function method. The main contributions are listed as follows:
(1)
Using appropriate variable substitution, the second-order equation can be transformed into a low-order two-dimensional system so that the dynamic properties of the second-order equation can be conveniently studied.
(2)
Compared with the article [5], the existence and stability conditions of periodic solutions obtained in this paper are easier to verify.
(3)
In [2], strong mathematical analysis skills are required, while the method used in this paper does not require strong mathematical analysis skills.
We organize the remaining research content of this article as follows: Section 2 lists some existing research results and lemmas needed in this paper. In Section 3 and Section 4, we give the main results of this paper. In Section 5, we give a numerical example. Finally, we summarize the research content of the full text and provide a discussion of some issues.

2. Preliminaries

Set U and V be two Banach spaces. Set  M : D ( M ) U V  as a Fredholm operator with index zero, which means that  d i m K e r M = c o d i m I m M < +  and  I m M  is closed in V. If M is a Fredholm operator with index zero, then there are continuous projectors  G : U U , H : V V  such that  M D ( M ) K e r G : ( I G ) X I m M  is invertible, and  I m M = K e r H = I m ( I H ) , I m G = K e r M . Denote, by  F G , the inverse of  M D ( M ) K e r G .
Let  Θ  be a bounded subset of U. A map  R : Θ ¯ V  is said to be L-compact in  Θ ¯  if the operator  F G ( I H ) R ( Θ ¯ )  is relatively compact and  H R ( Θ ¯ )  is bounded.
Lemma 1 
([18]). Assume that U and V are two Banach spaces, and  M : D ( M ) U V  is a Fredholm operator with index zero.  R : Θ ¯ V  is L-compact on  Θ ¯ , and  Θ U  is an open bounded set. If the following conditions are satisfied:
(1) 
M x μ R x , x Θ D ( M ) , μ ( 0 , 1 ) ,
(2) 
R x I m M , x Θ K e r M ,
(3) 
d e g { Φ H R , Θ K e r M , 0 } 0 ,
where  Φ : I m H K e r M  is a isomorphic mapping, then the operator equation  M x = R x  has a solution on  Θ ¯ D ( M ) .
Let  u ( t ) = y ( t ) + ξ y ( t ) , where  ξ > 0  is a constant. We rewrite Equation (3) as follows:
y ( t ) = ξ y ( t ) + u ( t ) , u ( t ) = a u ( t ) + ( ξ + a ) y ( t ) δ ( t ) y β ( t ) + ρ ( t ) y α ( t ) .
Furthermore, let  y ( t ) = e v ( t ) , and we rewrite (5) as follows:
v ( t ) = ξ + u ( t ) e v ( t ) , u ( t ) = a u ( t ) + ( ξ + a ) e v ( t ) δ ( t ) e β v ( t ) + ρ ( t ) e α v ( t ) .
It is easy to see that, if  z ( t ) = ( v ( t ) , u ( t ) ) T  is a T-periodic solution of system (6), then  y ( t ) = e v ( t )  is a positive T-periodic solution of Equation (3). For studying the periodic solution of Equation (3), we only need to study the periodic solution of system (6).
Set  C T = { ϑ C ( R , R ) : ϑ ( t + T ) = ϑ ( t ) }  with norm  | ϑ | 0 = max t [ 0 , T ] | ϑ ( t ) | . Set
U = V = { w = ( u ( · ) , v ( · ) ) T C ( R , R 2 ) : ω ( t + T ) = w ( t ) }
with norm  | | w | | = max t [ 0 , T ] { | u | 0 , | v | 0 } . Clearly, U and V are Banach spaces.
We define a linear operator by
M : D ( M ) U V , M z = z ( t ) = v ( t ) u ( t ) .
Furthermore, We define a nonlinear operator by
R : U V , R z = ξ + u ( t ) e v ( t ) a u ( t ) + ( ξ + a ) e v ( t ) δ ( t ) e β v ( t ) + ρ ( t ) e α v ( t ) .
We find
K e r M = R 2 , I m M = { ν | ν V , 0 T ν ( s ) d s = 0 } .
Thus,  d i m K e r M = c o d i m I m M = 2 , and  I m M  is closed in V. Therefore, M is a Fredholm operator with index zero. We define a continuous projector G by
G : U K e r M , G z = 1 T 0 T z ( s ) d s
and define the other continuous projector H by
H : Y V / I m M , H ψ = 1 T 0 T ψ ( s ) d s .
Set
M G = M | D ( M ) K e r G : D ( M ) K e r G I m M ,
and then
M G 1 = F G : I m M D ( M ) K e r G .
We have
( F G z ) ( u ) = 0 T Γ ( u , v ) z ( v ) d v ,
where
Γ ( u , v ) = v T , 0 v u T v T T , 0 u v T
Therefore,  F G  is a completely continuous operator in  I m M . Clearly,  H R ( Θ ¯ )  is bounded on  Θ ¯ . Thus, the operator R is L-compact on  Θ ¯ .
Throughout this paper, we give the following assumptions:
Assumption 1 (H 1 ). 
a 0 , δ ( t )  and  ρ ( t )  are T-periodic continuous functions with  δ ( t ) 0  and  ρ ( t ) 0 .
Assumption 2 (H 2 ). 
β > 1  and  α > 0 .
Assumption 3 (H 3 ). 
There is a positive constant d such that
σ ( ξ + a ) e σ ( t ) + δ ( t ) e β σ ( t ) ρ ( t ) e α σ ( t ) > 0 f o r σ > d
and
σ ( ξ + a ) e σ ( t ) δ ( t ) e β σ ( t ) + ρ ( t ) e α σ ( t ) > 0 f o r σ < d .

3. Existence of a Positive Periodic Solution

Theorem 3. 
Suppose assumptions (H 1 ), (H 2 ) and (H 3 ) are satisfied. Then, there exists at least one positive T-periodic solution for Equation (3), provided that
ξ a a ξ + ξ 1 β 1 N 2 ,
and
ln N 2 ξ N 3 ,
where  ξ > 0  is a constant, and  N 2  and  N 3  are defined by (20) and (22), respectively.
Proof. 
Consider the following equation:
M z = μ R z , μ ( 0 , 1 ) .
Take  Ω 1 = { z D ( M ) : M z = μ R z , μ ( 0 , 1 ) } , where M and R are expressed by (7) and (8), respectively. For  z = ( v ( t ) , u ( t ) ) T Ω 1 , it follows, by (11), that
v ( t ) = μ ξ + μ u ( t ) e v ( t )
and
u ( t ) = μ a u ( t ) + μ ( ξ + a ) e v ( t ) μ δ ( t ) e β v ( t ) + μ ρ ( t ) e α v ( t ) .
Integrate (12) over  [ 0 , T ] , and then
0 T u ( t ) e v ( t ) d t = ξ T .
In view of (14), there is  t 0 [ 0 , T ]  such that
u ( t 0 ) = ξ T e v ( t 0 ) .
On the other hand, by (12), we find
e v ( t ) v ( t ) = μ ξ e v ( t ) + μ u ( t ) ,
i.e.,
e v ( t ) = μ ξ e v ( t ) + μ u ( t ) .
For the above  t 0 , by (15) and (16), we find
e v ( t ) = e v ( t 0 ) e t 0 t μ ξ d τ + t 0 t μ u ( s ) e t 0 t μ ξ d τ d s = e v ( t 0 ) e μ ξ ( t t 0 ) + μ e μ ξ ( t t 0 ) t 0 t μ u ( s ) d s e v ( t 0 ) + e ξ ( T t 0 ) ( T t 0 ) | u | 0 1 ξ T + e ξ T T | u | 0 = N 1 | u | 0 ,
where  N 1 = 1 ξ T + e ξ T T .  From (13), (17), (H 1 ) and (H 2 ), we have
u ( t ) = u ( t 1 ) e t 1 t μ a d τ + t 1 t μ [ ( ξ + a ) e v ( s ) δ ( s ) e β v ( s ) + ρ ( s ) e α v ( s ) ] e t 1 t μ a d τ d s = u ( t 1 ) e μ a ( t t 1 ) + μ e μ a ( t t 1 ) t 1 t μ [ ( ξ + a ) e v ( s ) δ ( s ) e β v ( s ) + ρ ( s ) e α v ( s ) ] d s | u | 0 + e a ( T t 1 ) ( T t 1 ) ( ξ + a ) N 1 | u | 0 δ ˇ N 1 β | u | 0 β ,
where  δ ˇ = min t [ 0 , T ] | δ ( t ) | . By (18), we have
δ ˇ N 1 β | u | 0 β ( ξ + a ) N 1 | u | 0 .
By (19), we have
| u | 0 ( ξ + a ) N 1 δ ˇ N 1 β 1 β 1 = N 2 .
From (17) and (20), we have
e v ( t ) N 1 | u | 0 N 1 N 2 .
By (21), we have
| v | 0 ln ( N 1 N 2 ) = N 3 .
Let  Ω 2 = { z K e r M : R z I m M } . If  z Ω 2 ,  and then  z K e r M  and  H R z = 0 , which results in
u e v = ξ
and
a u T + ( ξ + a ) e v T = 0 T [ δ ( t ) e β v ρ ( t ) e α v ] d t .
From (23) and (24), we have
a u + ( ξ + a ) e v δ ˇ e β v = δ ˇ u β ξ β .
Thus, by (9) and (25), we find
u ξ β β 1 a ξ a + 1 1 β 1 N 2 .
It follows by (10), (11) and (26) that
v = ln u ξ ln N 2 ξ N 3 .
Let  Ω = { z = ( v , u ) T : | v | 0 M 1 , | u | 0 M 2 } , where  M 1 > N 3 , M 2 > N 2 . Clearly,  Ω 1 , Ω 2 Ω .  In view of (20), (22), (26) and (27), we verified the conditions (1) and (2) in the Lemma. Now, we claim that the condition (3) of Lemma 1 holds. Let the isomorphic mapping  Φ : I m H K e r M  be
Φ ( v , u ) = ( u , v ) for v > d , ( u , v ) for v < d .
and
H ( ω , λ ) = λ ω + 1 λ T Φ H R ω for ( ω , λ ) Ω × [ 0 , 1 ] .
When  v > d , for  ( z , λ ) ( Ω K e r M ) , by assumption (H 3 ), we have
z T H ( z , λ ) = λ ( v 2 + u 2 ) + ( 1 λ ) T [ a u v + u ξ u 2 e v ] + ( 1 λ ) T 2 0 T v ( ξ + a ) e v ( t ) + δ ( t ) e β v ( t ) ρ ( t ) e α v ( t ) d t > 0 .
Furthermore, when  v < d , for  ( z , λ ( Ω K e r M ) , by assumption (H 3 ), we find
z T H ( z , λ ) = λ ( v 2 + u 2 ) + ( 1 λ ) T [ a u v + u ξ u 2 e v ] + ( 1 λ ) T 2 0 T v ( ξ + a ) e v ( t ) δ ( t ) e β v ( t ) + ρ ( t ) e α v ( t ) d t > 0 .
For  λ [ 0 , 1 ]  and  z Ω K e r M , we find  z T H ( z , μ ) 0 .  Thus,
d e g Φ H R , Ω K e r M , 0 = d e g H ( z , 0 ) , Ω K e r M , 0 = d e g H ( z , 1 ) , Ω K e r M , 0 0 .
Hence, the condition (3) of Lemma 1 is satisfied. Applying Lemma 1, there exists a T-periodic solution  ( v , u ) T  for system (6). Therefore, there exists a positive T-periodic solution  x ( t ) = e v ( t )  for Equation (3). □
Remark 1. 
The coincidence degree method is one of the main methods to investigate differential equations. When using this method, we should first estimate the prior bound of the solution. In many cases, estimating the prior bound of the solution is the primary difficulty in using this method. For more results about the use of the coincidence degree method, see, e.g., [19,20,21,22,23] and related references.

4. Dynamic Properties of a Positive Periodic Solution

In this section, we consider the dynamic properties of system (5) due to the equivalence of (3) and (5).
Definition 1. 
If  u ( t ) = ( u 1 ( t ) , v 1 ( t ) ) T  is a periodic solution of system (5) and  u ( t ) = ( u 1 ( t ) , v 1 ( t ) ) T  is any solution of system (5) satisfying
lim t + | u 1 ( t ) u 1 ( t ) | + | v 1 ( t ) v 1 ( t ) | = 0 ,
then we call  u ( t )  globally asymptotic stable.
We give the following assumption:
Assumption 4 (H 4 ). 
Set  f ( u ) = u γ ,  where  u R , γ > 0 . There exists  L > 0  such that
| f ( u 1 ) f ( v 1 ) | L | u 1 v 1 | for each u 1 , v 1 R .
Theorem 4. 
Assume that all conditions of Theorem 3 and (H 4 ) are satisfied. Then, Equation (3) has a globally asymptotic stable positive periodic solution, provided that
1 1 ξ 2 ( ξ + a + δ ^ L + ρ ^ L ) > 0 ,
where  δ ^ = max t [ 0 , T ] | δ ( t ) |  and  ρ ^ = max t [ 0 , T ] | ρ ( t ) | .
Proof. 
Using the conditions of Theorem 3, we obtain that that system (5) has a bounded T-periodic solution  u ( t ) = ( u 1 ( t ) , v 1 ( t ) ) T . Assume  u ( t ) = ( u 1 ( t ) , v 1 ( t ) ) T  to be any solution of system (5). By (5), we find
( u 1 ( t ) u 1 ( t ) = ξ ( u 1 ( t ) u 1 ( t ) + ( v 1 ( t ) v 1 ( t ) , ( v 1 ( t ) v 1 ( t ) = a ( v 1 ( t ) v 1 ( t ) + ( ξ + a ) ( u 1 ( t ) u 1 ( t ) δ ( t ) ( u 1 β ( t ) ( u 1 ) β ( t ) + ρ ( t ) [ ( u 1 α ( t ) ( u 1 ) α ( t ) ] .
From the first equation of system (29), we have
u 1 ( t ) u 1 ( t ) = u 1 ( t 0 ) u 1 ( t 0 ) e t 0 t ξ d s + t 0 t e t s ξ d s ( v 1 ( s ) v 1 ( s ) d s = u 1 ( t 0 ) u 1 ( t 0 ) e ξ ( t t 0 ) + t 0 t e ξ ( s t ) ( v 1 ( s ) v 1 ( s ) d s ,
where  t t 0  and  t 0  is a given constant. By (30), we have
| u 1 ( t ) u 1 ( t ) | | u 1 ( t 0 ) u 1 ( t 0 ) | e ξ ( t t 0 ) + 1 ξ 1 e ξ ( t t 0 ) | v 1 v 1 | 0 .
By the use of the second equation in system (29), we obtain
v 1 ( t ) v 1 ( t ) = v 1 ( t 0 ) v 1 ( t 0 ) e t 0 t a d s + t 0 t e t s a d s [ ( ξ + a ) ( u 1 ( s ) u 1 ( s ) δ ( s ) ( u 1 β ( s ) ( u 1 ) β ( s ) + ρ ( t ) [ ( u 1 α ( s ) ( u 1 ) α ( s ) ] ] d s = v 1 ( t 0 ) v 1 ( t 0 ) e a ( t t 0 ) + t 0 t e a ( s t ) [ ( ξ + a ) ( u 1 ( s ) u 1 ( s ) s ( s ) ( u 1 β ( s ) ( u 1 ) β ( s ) + ρ ( s ) [ ( u 1 α ( s ) ( u 1 ) α ( s ) ] ] d s
which, together with assumption (H 4 ), results in
| v 1 ( t ) v 1 ( t ) | | v 1 ( t 0 ) v 1 ( t 0 ) | e a ( t t 0 ) + 1 ξ 1 e a ( t t 0 ) ( ξ + a + δ ^ L + ρ ^ L ) | u 1 u 1 | 0 .
From (31) and (32), we have
1 1 ξ 2 ( 1 e ξ ( t t 0 ) ) ( 1 e a ( t t 0 ) ) ( ξ + a + δ ^ L + ρ ^ L ) | u 1 u 1 | 0 | u 1 ( t 0 ) u 1 ( t 0 ) | e ξ ( t t 0 ) + 1 ξ 1 e ξ ( t t 0 ) | v 1 ( t 0 ) v 1 ( t 0 ) | e a ( t t 0 ) .
Let  t +  for (33), and then
lim t + max t R | u 1 ( t ) u 1 ( t ) | = 0 .
By (34), we find
lim t + | u 1 ( t ) u 1 ( t ) | = 0 .
In view of (35), let  t +  for (32), and then
lim t + | v 1 ( t ) v 1 ( t ) | = 0 .
By (35) and (36), we have
lim t + [ | u 1 ( t ) u 1 ( t ) | + | v 1 ( t ) v 1 ( t ) | ] = 0 .
In view of (37), there exists a globally asymptotic stable solution  u ( t ) = ( u 1 ( t ) , v 1 ( t ) ) T  for system (5). Thus, there exists a globally asymptotic stable solution  u 1 ( t )  for Equation (3). □
Remark 2. 
In [2], the authors indicated that a pump belongs to a periodically forced differential system. Hence, studying the asymptotic properties of periodic solutions has important significance to pipe/tank flow configurations. In [2], to obtain an asymptotically stable positive solution of Equation (3), constructing upper and lower and functions is very difficult. In [5], exponential stability of a periodic solution for (3) was obtained. However, the obtained conditions in [5] are strong and are not easy to verify. In this article, we obtain conditions for the existence of periodic solutions that are easy to verify.

5. Example

Equation (3) has wide applications in pipe/tank configurations. Consider the following nonlinear differential system:
y ( t ) = ξ y ( t ) + u ( t ) , u ( t ) = a u ( t ) + ( ξ + a ) y ( t ) δ ( t ) y β ( t ) + ρ ( t ) y α ( t ) ,
where  T = 2 π , ξ = 2 , a = 1.9 , δ ( t ) = 0.75 × 10 3 > 0 , ρ ( t ) = sin t 2 < 0 , α = 0.1  and  β = 1.5 .  Clearly, assumptions (H 1 ) and (H 2 ) hold. From  δ ( t ) > 0  and  ρ ( t ) < 0 , assumption (H 3 ) holds. After a simple calculation, we have
N 1 = 1 ξ T + e ξ T T 1.79 × 10 6 , N 2 = ( ξ + a ) N 1 δ ˇ N 1 β 1 β 1 16 , N 3 = ln ( N 1 N 2 ) 17.17
Thus,
ξ β β 1 a ξ a + 1 1 β 1 0.02 N 2 = 16 ,
and
ln N 2 ξ 5.05 17.17 = N 3 .
Thus, all conditions of Theorem 3 are satisfied. Hence, there is a positive periodic solution for (3). Furthermore, choosing  L = 1 × 10 3 , we have
1 1 ξ 2 ( ξ + a + s ^ L + r ^ L ) = 0.03 > 0
and condition (28) holds. Hence, all conditions of Theorem 4 hold. From Theorem 4, there is a globally asymptotic stable positive periodic solution  ( y , u ) T  for (38). Hence, there is globally asymptotic stable positive periodic solution x.
System (38) has a globally asymptotic stable positive periodic solution  ( y , u ) T , i.e., Equation (3) has a globally asymptotic stable positive periodic solution. Figure 1 supports this conclusion.
Remark 3. 
In [5], the authors obtained the following results for (3).
Theorem 5. 
Assume that  ρ C ( [ t 0 , ) , ( 0 , ) ) , a > 0 , 0 < α < β < 1 , c > 0 , and there is function  ϑ C ( [ t 0 , ) , R )  and constants  L , l  such that the following conditions are satisfied:
0 < l exp t 0 t [ a + ϑ ( s ) ] d s L for t t 0 ,
t t + T [ a + ϑ ( s ) ] d s = 0 for t t 0 ,
ϑ ( t ) exp t 0 t [ a + ϑ ( s ) ] d s = t 0 t [ δ ( s ) exp α t 0 s [ a + ϑ ( v ) ] d v ρ ( s ) t 0 t exp β t 0 s [ a + ϑ ( v ) ] d v ] d s for t t 0 .
In addition, there are constants  L , l > 0  such that  l l , L L  and
α l α 1 δ ( t ) β L β 1 ρ ( t ) 0 for t t 0 ,
where  δ ( t ) = c μ , ρ ( t ) = e ( t ) μ . Then, there is an exponentially stable positive periodic solution for (3).
Clearly, conditions (39)–(41) are strong and difficult to verify. In [5], the authors obtained the existence of positive periodic solution of Equation (3) by using Schauder’s fixed point theorem. Furthermore, they obtained exponential stability of a positive periodic solution of Equation (3) by using the Lyapunov function method. The numerical simulation in [5] supports these conclusions. In the present paper, under more relaxed conditions, the numerical simulation supports our conclusions. Therefore, the results of this paper improve and enhance the results of [5].

6. Conclusions and Discussion

This paper considered a positive periodic solution of a nonlinear equation related to valveless pumping. The second-order equation reduced to a two-dimensional system via a linear variable transformation. We obtained some results for a periodic solution to (3). A simulation example demonstrated the validity of our results.
Since Equation (3) is a singular second-order equation, in future research work, we will focus on studying how the singularity of the equation affects the existence of the solution. In recent years, stochastic differential systems have received great attention. Hu, Zhu and Karimi [24] studied the Razumikhin stability theorem for a class of impulsive stochastic delay differential systems by using the stochastic analysis technique and Razumikhin approach. In [14], the authors dealt with the exponential stability of an impulsive stochastic food chain system with time-varying delays. The stability problems of stochastic nonlinear delay systems with exogenous disturbances and event-triggered feedback control were investigated in [17]. In the future, we will study stochastic pipe/tank flow configurations with friction.

Author Contributions

Methodology, H.D.; Formal analysis, X.W.; Writing—original draft, B.D. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was supported by the Natural Science Foundation of Huai’an (HAB202231).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Propst, G. Pumping effects in models of periodically forced flow configurations. Phys. D Nonlinear Phenom. 2006, 217, 193–201. [Google Scholar] [CrossRef]
  2. Cid, J.Á.; Propst, G.; Tvrdý, M. On the pumping effect in a pipe/tank flow configuration with friction. Phys. D Nonlinear Phenom. 2014, 273, 28–33. [Google Scholar] [CrossRef]
  3. Hakl, R.; Torres, P.J.; Zamora, M. Periodic solutions of singular second order differential equations: The repulsive case. Topol. Methods Nonlinear Anal. 2012, 39, 199–220. [Google Scholar]
  4. Hakl, R.; Torres, P.J.; Zamora, M. Periodic solutions of singular second order differential equations: Upper and lower functions. Nonlinear Anal. 2011, 74, 7078–7093. [Google Scholar] [CrossRef]
  5. Dorociaková, B.; Michalková, M.; Olach, R.; Sága, M. Existence and Stability of Periodic Solution Related to Valveless Pumping. Math. Probl. Eng. 2018, 2018, 3982432. [Google Scholar] [CrossRef] [Green Version]
  6. Liu, W.; Huang, J.; Yao, Q. Stability analysis for quaternion-valued inertial memristor-based neural networks with time delays. Neurocomputing 2021, 448, 67–81. [Google Scholar] [CrossRef]
  7. Novaes, D.; Silva, F. Higher Order Analysis on the Existence of Periodic Solutions in Continuous Differential Equations via Degree Theory. Siam J. Math. Anal. 2021, 53, 2476–2490. [Google Scholar] [CrossRef]
  8. Liu, Y.; Yang, P.; Ge, W. Periodic solutions of higher-order delay differential equations. Nonlinear Anal. TMA 2005, 63, 136–152. [Google Scholar] [CrossRef]
  9. Liu, J.; Liu, W.; Liu, B. Periodic solutions for fourth-order p-Laplacian functional differential equations with sign-variable coefficient. Electron. J. Differ. Equ. 2013, 205, 1–9. [Google Scholar]
  10. Djebali, I.; Guedda, L. Fractional multipoint boundary value problems at resonance with kernel dimension greater than one. Math. Methods Appl. Sci. 2021, 44, 2621–2636. [Google Scholar] [CrossRef]
  11. Rao, R.; Lin, Z.; Ai, X.; Wu, J. Synchronization of epidemic systems with Neumann boundary value under delayed impulse. Mathematics 2022, 10, 2064. [Google Scholar] [CrossRef]
  12. Li, G.; Zhang, Y.; Guan, Y.; Li, J. Stability analysis of multi-point boundary conditions for fractional differential equation with non-instantaneous integral impulse. Math. Biosci. Eng. 2023, 20, 7020–7041. [Google Scholar] [CrossRef]
  13. Zhu, Q.; Kong, F.; Cai, Z. Advanced Symmetry Methods for Dynamics, Control, Optimization and Applications. Symmetry 2023, 15, 26. [Google Scholar] [CrossRef]
  14. Zhao, Y.; Wang, L. Practical exponential stability of impulsive stochastic food chain system with time-varying delays. Mathematics 2023, 11, 147. [Google Scholar] [CrossRef]
  15. Li, K.; Li, R.; Cao, L.; Feng, Y.; Onasanya, B.O. Periodically intermittent control of Memristor-based hyper-chaotic bao-like system. Mathematics 2023, 11, 1264. [Google Scholar] [CrossRef]
  16. Xia, M.; Liu, L.; Fang, J.; Zhang, Y. Stability analysis for a class of stochastic differential equations with impulses. Mathematics 2023, 11, 1541. [Google Scholar] [CrossRef]
  17. Zhu, Q. Stabilization of stochastic nonlinear delay systems with exogenous disturbances and the event-triggered feedback control. IEEE Trans. Autom. Control 2019, 64, 3764–3771. [Google Scholar] [CrossRef]
  18. Gaines, R.; Mawhin, J. Coincidence Degree and Nonlinear Differential Equations; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
  19. Lu, S. On the existence of positive peroiodic solutions for neutral functional differential equation with multiple deviating arguments. J. Math. Anal. Appl. 2003, 280, 321–333. [Google Scholar] [CrossRef] [Green Version]
  20. Yang, H.; Han, X. Existence of periodic solutions for the forced pendulum equations of variable length. Qual. Theory Dyn. Syst. 2023, 22, 20. [Google Scholar] [CrossRef]
  21. Shi, L.; Qi, L.; Zhai, S. Periodic and almost periodic solutions for a non-autonomous respiratory disease model with a lag effect. Acta Math. Sci. 2022, 42, 187–211. [Google Scholar] [CrossRef]
  22. Gao, F.; Zhang, W. Periodic solutions for a p-Laplacian-like NFDE system. J. Frankl. Inst. 2011, 348, 1020–1034. [Google Scholar] [CrossRef]
  23. Ge, W.; Ren, J. An extension of Mawhins continuation theorem and its application to boundary value problems with a p-Laplacain. Nonlinear Anal. 2004, 58, 477–488. [Google Scholar] [CrossRef]
  24. Hu, W.; Zhu, Q.; Karimi, H. Some improved Razumikhin stability criteria for impulsive stochastic delay differential systems. IEEE Trans. Autom. Control 2019, 64, 5207–5213. [Google Scholar] [CrossRef]
Figure 1. Positive periodic solution  ( y ( t ) , u ( t ) )  of system (38).
Figure 1. Positive periodic solution  ( y ( t ) , u ( t ) )  of system (38).
Mathematics 11 01789 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Du, H.; Wang, X.; Du, B. Positive Periodic Solution for Pipe/Tank Flow Configurations with Friction. Mathematics 2023, 11, 1789. https://0-doi-org.brum.beds.ac.uk/10.3390/math11081789

AMA Style

Du H, Wang X, Du B. Positive Periodic Solution for Pipe/Tank Flow Configurations with Friction. Mathematics. 2023; 11(8):1789. https://0-doi-org.brum.beds.ac.uk/10.3390/math11081789

Chicago/Turabian Style

Du, Haiqing, Xiaojing Wang, and Bo Du. 2023. "Positive Periodic Solution for Pipe/Tank Flow Configurations with Friction" Mathematics 11, no. 8: 1789. https://0-doi-org.brum.beds.ac.uk/10.3390/math11081789

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop