Next Article in Journal
On a Backward Problem for the Rayleigh–Stokes Equation with a Fractional Derivative
Next Article in Special Issue
Iteration with Bisection to Approximate the Solution of a Boundary Value Problem
Previous Article in Journal
Path Planning and Trajectory Tracking for Autonomous Obstacle Avoidance in Automated Guided Vehicles at Automated Terminals
Previous Article in Special Issue
On Boundary Controllability for the Higher-Order Nonlinear Schrödinger Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Behavior of Solutions for Nonlinear Klein–Gordon Wave Type Models with a Logarithmic Nonlinearity and Multiple Time-Varying Delays

Institut de Recherche Mathématique de Rennes, 20 Avenue des Buttes de Coësmes, CS 14315, 35043 Rennes Cédex, France
Submission received: 1 November 2023 / Revised: 20 December 2023 / Accepted: 26 December 2023 / Published: 30 December 2023
(This article belongs to the Special Issue Advances in Nonlinear Analysis and Boundary Value Problems)

Abstract

:
In this paper, we study the existence and exponential stability of solutions to a class of nonlinear delay Klein–Gordon wave type models on a bounded domain. Such models include multiple time-varying delays, frictional damping, and nonlinear logarithmic source terms. After showing the local existence result of the solutions using Faedo–Galerkin’s method and logarithmic Sobolev inequality, the global existence is analyzed. Then, under some appropriate conditions, energy decay estimates and exponential stability results of the global solutions are investigated.

1. Introduction and Mathematical Setting of the Problem

In this paper, we consider the following nonlocal initial-boundary value problem for a class of nonlinear Klein–Gordon wave type equations with frictional damping, logarithmic nonlinearity, and multiple time-varying delays in velocity:
2 u t 2 div ( K ( x ) u ) + d 0 ( t ) u t + Ψ ( x , t ; u t ) + M ( u L 2 ( Ω ) 2 ) u = u ln u θ , x Ω , t > 0 , u t ( x , t ) = h 0 ( x , t ) , ( x , t ) Q 0 = Ω × [ δ ( 0 ) , 0 ) , u ( x , 0 ) = ϕ 0 ( x ) , u t ( x , 0 ) = ϕ 1 ( x ) , x Ω , u ( x , t ) = 0 , x Γ , t > 0 ,
where  u ( x , t ) I R  is the variable state, the operator  M  is a continuous function on  [ 0 , + ) , the domain  Ω  is a boundary subset of  I R N N 1 , with smooth boundary  Γ = Ω , and the term  K  is the conductivity tensor function on  Ω ¯ , with  Ω ¯ = Ω Ω  the closure of  Ω . The functions  ϕ 0 ϕ 1 , and  h 0  are the initial/initial history data to be specified later. The operator  Ψ , which describes multiple time-varying delays related to velocity  u t , is defined as in [1]
Ψ ( x , t ; u t ) = i = 1 n d i ( t ) u t ( x , t e i ( t ) ) ,
where the non constant weights  ( d i , i = 1 , n ) : I R + I R  and  d 0 : I R + ] 0 , + [  are bounded and sufficiently regular functions,  ( e i , i = 1 , n ) : I R + I R +  are sufficiently regular functions representing multiple time-varying delays,  δ ( 0 ) = max i = 1 , n ( e i ( 0 ) ) , and the parameter  θ  (which measures the force of nonlinear interactions) is a real number that will be specified later.

1.1. Motivation and Outline of the Paper

A strong motivation for analyzing the behavior of the nonlinear problem with multiple time-varying delays (1) comes from its application in many branches of physics and other applied sciences.
The nonlinear Klein–Gordon equation, also known as the relativistic version of the Schrodinger equation, is an important class of partial differential equations and performs a significant role in mathematical physics and many other scientific applications such as relativistic quantum mechanics. They occur in various areas of physical sciences and engineering such as solid-state physics, nonlinear optics, quantum field theory, fluid dynamics, mathematical biology, chemical kinematics, propagation of dislocations in crystals, and the behavior of elementary particles. The Klein–Gordon equation, with and without damping terms, has been extensively studied, either from a theoretical (as global existence and nonexistence, exponential decay of energy, time blow-up, asymptotic behavior of solutions) or from a numerical point of view, for different nonlinear source terms  F ( u )  as  F ( u ) = a u u γ  or logarithmic nonlinearity  F ( u ) = a u l n u γ , where the parameter a measures the force of the nonlinear interactions. For the first type of nonlinearities, we can cite, e.g., [2,3,4,5] and the references therein. Problems with logarithmic nonlinearity arise naturally in many areas such as quantum optics and transport phenomena, via a logarithmic Schrodinger equation (see, e.g., [6,7]); fluid dynamics via a logarithmic Korteweg–de Vries equation or a logarithmic Kadomtsev–Petviashvili equation (see, e.g., [8]); or material sciences, with a Cahn–Hilliard equation (see, e.g., [9]). It also arises in nuclear physics, inflation cosmology, vibration, supersymmetric fields in quantum field theory, spinless particles, and viscoelastic mechanics where a logarithmic Klein–Gordon equation is considered (see, e.g., [7,10,11,12,13,14,15,16]). Such logarithmic Klein–Gordon problems have been the object of numerous studies either from a theoretical or from a numerical point of view (see, e.g., [17,18,19,20,21,22,23,24,25,26,27,28,29,30] and the references therein).
The introduction of retarded arguments is to reflect the different after-effects. Different time-varying delay configurations occur naturally in various areas of physics, biologics and engineering, such as in biochemical systems, in population dynamics, in quantum chaotic systems, in relativistic quantum waves, and in the area of plasma control (e.g., in the context of thermonuclear fusion with Tokamaks). Moreover, time-varying delays in signal transmission are inevitable in many applications and practical processes. A small delay can affect considerably dynamical behaviors of the system (e.g., destabilize the system which is asymptotically stable in the absence of time delays unless additional conditions, control functions, or stabilization mechanism functions have been used). Delay terms can lead to change in the stability of dynamics and give rise to highly complex behavior including instability, oscillations, and chaos (see, e.g., [1,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46] and the references therein). Therefore, these behaviors and aspects, by taking into account different sources of delays, motivate the study of multiple time-varying delays effects on properties of dynamical systems.
In this work, in order to take into account the influence of different sources of time delays in the velocity of signal transmission, nonlinear Klein–Gordon wave type models with logarithmic nonlinearity are modified by incorporating multiple time delays and a nonlocal operator. The proposed strategy consists in controlling these instabilities by imposing some suitable conditions involving different functions and parameters representing the multiple time-varying delays. The presence of the nonlocal term  M ( u L 2 ( Ω ) 2 )  and operator  Ψ , which describes multiple time-varying delays, makes the mathematical study of such a class of problems particularly interesting.
Remark 1.1.
1.
The Equation (1) is called nonlocal because of the presence of nonlocal term  M ( u L 2 ( Ω ) 2 ) .
2.
The functions  d k  are diffusion coefficients that represent the strength of each associated time delay. A zero coefficient means the associated previous state does not impact the system. □
The paper is organized as follows. In the next subsections, we give some necessary notations and preliminary results. In Section 2, we prove the local existence of solutions to problem (1) using Faedo–Galerkin’s approximation and logarithmic Sobolev inequality. Section 3 deals with the global existence and energy decay rate of the solutions to problem (1). In Section 4, under suitable conditions on data and involved functionals, the exponential decay of solutions, for initial data in a set of stability, is investigated. Finally, we present conclusions in Section 5.

1.2. Notations

Let  D  be a Banach space and  D  its dual Banach space. We denote the norm on  D  by  . D , the norm on  D  by  . D , and the duality pairing on  D  and  D  by  < , > D , D . We denote the norm and the scalar product in  L 2 ( Ω )  by  .  and  ( . , . ) , respectively. We say that a sequence  ( w n ) n  of  D  (respectively, of  D ) converges weakly (respectively, weakly*) to w if and only if  f D < f , w n w >  converges to 0 (respectively,  v D < w n w , v >  converges to 0) (see e.g., [33], Part I). The dual of  L 2 ( Ω )  is identified with itself, the dual of  H 0 1 ( Ω )  is  H 1 ( Ω )  and we have the following injections  H 0 1 ( Ω ) L 2 ( Ω ) H 1 ( Ω )  with continuous and dense embedding. Moreover we denote by  C Ω  (which depends on the geometry of domain  Ω ) the optimal constant of embedding inequality (this smallest possible  C Ω  is called the Poincaré constant)
v L 2 ( Ω ) C Ω v L 2 ( Ω ) , v H 0 1 ( Ω ) .
If  Ω  is a convex domain with diameter  d Ω , then  C Ω = d Ω π .
We can now introduce the following sets (for arbitrary final time  T > 0 ):
V 0 ( Ω ; δ ( 0 ) ) = H 0 1 ( Ω ) × L 2 ( Ω ) × L 2 ( Q 0 ) , H ( 0 , T ) = L ( 0 , T ; L 2 ( Ω ) ) , V ( 0 , T ) = L 2 ( 0 , T ; H 0 1 ( Ω ) ) , S ( 0 , T ) = C 0 ( [ 0 , T ] ; H 0 1 ( Ω ) ) C 1 ( [ 0 , T ] ; L 2 ( Ω ) ) C 2 ( [ 0 , T ) ; H 1 ( Ω ) ) .
Finally, for a bounded function f, we denote
f m i n = inf t 0 f ( t ) , f m a x = sup t 0 f ( t ) .
We now state some assumptions for the various functions and operators appearing in Equation (1).

1.3. Assumptions and Preliminaries

We start by assuming that  t I R + ( r i ( t ) = t e i ( t ) , i = 1 , n )  are strictly increasing functions (and consequently are bijective) and  ( e i , i = 1 , n )  are  C 1  non-negative and bounded functions on  I R + . So, we have the existence of inverse functions  ( f i , i = 1 , n )  of  ( r i , i = 1 , n ) . For simplicity, without loss of generality, we can suppose that  e i ( 0 ) = δ ( 0 )  for all  i = 1 , , n .
Remark 1.2.
We can define the following subdivision:  s 1 = δ ( 0 ) s 0 = 0  and  j I N s j = min i = 1 , n ( f i ( s j 1 ) ) , and we denote  τ j = s j s j 1 , and  Q j = Ω × ( s j 1 , s j )  for  j 0 . According to the hypotheses on functions  ( e i , i = 1 , n ) , we prove easily that:
(i)
the sequence  ( s j ) j I N  is strictly increasing  j 0 ,
(ii)
for  j 2 , if  t ( s j 1 , s j )  then  i = 1 , , n r i ( t ) s j 1 ,
(iii)
if  t ( s 0 , s 1 ) , then  i = 1 , , n r i ( t ) ( s 1 , s 0 ) .
In order to derive the solution of (1), we can use the following process: we solve the problem on  Q 1  by using the initial/initial history data and obtain the solution. Then, the solution on  Q 2  is obtained by using the solution on  Q 1  to generate the initial data at  s 1 . This advancing process is repeated for  Q 3 ,  Q 4 , …, until the final set is reached. 
For the tensor function  K , we suppose that the following assumptions hold.
(H1) We assume that the conductivity tensor function  K W 1 , ( Ω ¯ )  is a symmetric, positive definite matrix function and is uniformly elliptic, i.e., there exist constants  0 < ν μ  such that ( v I R N )
ν v 2 2 v T K v μ v 2 2 in Ω ¯ ,
where  . 2  is the Euclidean norm.
For the nonlinear operator  M , we set the following hypothesis.
(H2)  M C ( [ 0 , ) , I R )  such that  M ( ζ ) ϑ 0 0 ζ [ 0 , ) , with  ϑ 0  a constant. We denote  M p ( ζ ) = 0 ζ M ( η ) d η .
Finally, we impose the following assumption for the parameter  θ : there exists  α > 0   (depending on  θ ) such that
(AAA) (i)  ν θ α 2 2 π > 0  and (ii)  θ ( 1 + N ( 1 + ln α ) ) + ϑ 0 2 > 0 .
Remark 1.3.
-
If  θ < 0 , the relation (i) is always true for any  α > 0  and (ii) is true for all α such that  0 < α < exp ( N 1 ϑ 0 θ N ) . Consequently, it is always possible to find  α > 0  such that the relations (i) and (ii) hold.
-
If  θ > 0 , the relation (i) is true for all α such that  0 < α < 2 π ν θ  and (ii) is true for all α such that  α > exp ( N 1 ϑ 0 θ N ) . Consequently it is always possible to find α such that the relations (i) and (ii) hold, provided that  2 π ν θ > exp ( N 1 ϑ 0 θ N ) , i.e., that θ satisfies the following inequality
θ Θ N ( ν , ϑ 0 ) ,
with  Θ N ( ν , ϑ 0 )  satisfying  2 π ν Θ N ( ν , ϑ 0 ) = exp ( N 1 ϑ 0 Θ N ( ν , ϑ 0 ) N ) . We prove easily that the sequence  Θ N ( ν , ϑ 0 )  is decreasing with  Θ ( ν , ϑ 0 ) = lim N Θ N ( ν , ϑ 0 ) = 2 π e 2 ν  (the limit  Θ ( ν , ϑ 0 ) 46.4268 ν  is independent of  ϑ 0 ). In this case, α satisfies
α 0 = exp ( N 1 ϑ 0 θ N ) < α < α 1 = 2 π ν θ .
Lemma 1.1
(Logarithmic Sobolev inequality, see, e.g., [47]). Let  v H 0 1 ( Ω ) , then for any positive constant α, the following estimates hold:
Ω v 2 ln v d x 1 2 v L 2 ( Ω ) 2 ln v L 2 ( Ω ) 2 + α 2 2 π v L 2 ( Ω ) 2 N ( 1 + ln α ) 2 v L 2 ( Ω ) 2 .
Lemma 1.2
(Logarithmic Gronwall inequality, see [19] or, e.g., [25]). Given a positive time T and positive constants K and  A 0  with  A 0 1 , let l be a nonnegative function such that  l L 1 ( ( 0 , T ) ) . If a function  A : ( 0 , T ) [ 1 , + [  satisfies, for any  t ( 0 , T ) ,
A ( t ) A 0 + K 0 t l ( s ) A ( s ) l n A ( s ) d s ,
then  A ( t ) A 0 exp ( K 0 t l ( s ) d s ) , for any  t ( 0 , T ) .
Lemma 1.3.
Let  ϵ ( 0 , 1 ) , then there exists  d ϵ > 0  such that  r ln r r 2 + d ϵ r 1 ϵ , for any  r > 0 .
 
Proof. 
We prove easily that  d ϵ  is the maximum of the following function  g : r ] 0 , + [ g ( r ) = r ϵ ( ln r r ) .   □
In the sequel, we will always denote C (or  C i ) as some positive constant, which may be different at each occurrence.

2. Existence of Local Solution

In this section, we shall study the existence of local solutions for problem (1). For this, let  T > 0  be a fixed and but arbitrary real number, and we denote  Q = Ω × ( 0 , T )  and  Σ = Γ × ( 0 , T ) .
Introducing now the following new functions (for  i = 1 , , n )
w i ( x , t ; η ) = u t ( x , t η e i ( t ) ) , f o r ( x , t ; η ) Q × ( 0 , 1 ) .
Then, we have
w i ( x , t ; 0 ) = u t ( x , t ) , w i ( x , t ; 1 ) = u t ( x , t e i ( t ) ) , and e i ( t ) w i t ( x , t ; η ) + ( 1 η e i ( t ) ) w i η ( x , t ; η ) = 0 , for ( x , t ; η ) Q × ( 0 , 1 ) .
Consequently, problem (1) becomes (for  i = 1 , , n )
2 u t 2 div ( K ( x ) u ) + d 0 ( t ) u t + i = 1 n d i ( t ) w i ( x , t ; 1 ) + M ( u L 2 ( Ω ) 2 ) u = u ln u θ , f o r ( x , t ) Q , e i ( t ) w i t ( x , t ; η ) + ( 1 η e i ( t ) ) w i η ( x , t ; η ) = 0 , f o r ( x , t ; η ) Q × ( 0 , 1 ) , w i ( x , t ; 0 ) = u t ( x , t ) f o r ( x , t ) Q , w i ( x , 0 ; η ) = h 0 ( x , η δ ( 0 ) ) : = w 0 i ( x ; η ) , f o r ( x , η ) Ω × ( 0 , 1 ) , u ( x , 0 ) = ϕ 0 ( x ) , u t ( x , 0 ) = ϕ 1 ( x ) , f o r x Ω , u ( x , t ) = 0 , f o r ( x , t ) Σ .
For this solution u, the corresponding energy function  E u : I R + I R  is defined by (for all  t 0 )
E u ( t ) = u t L 2 ( Ω ) 2 + ( K u , u ) + M p ( u L 2 ( Ω ) 2 ) + θ 2 u L 2 ( Ω ) 2 Ω u 2 ln u θ d x + λ i = 1 n e i ( t ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 ,
where  λ > 0  is a fixed parameter that satisfies the following condition (H3)
D 1 = 1 n inf t 0 ( d 0 ( t ) 1 2 i = 1 n d i ( t ) ) > λ > D 2 = max i = 1 , n ( sup t 0 d i ( t ) 2 ( 1 e i ( t ) ) ) , under the following additional assumption 1 n inf t 0 ( d 0 ( t ) 1 2 i = 1 n d i ( t ) ) > max i = 1 , n ( sup t 0 d i ( t ) 2 ( 1 e i ( t ) ) ) .
For typographical convenience, we will denote the energy at time t by  E ( t )  in place of  E u ( t )  if no confusion arises.
Remark 2.1.
If the weight function  d 0  is a non-increasing function, we can replace the energy function E by
E ( t ) = u t L 2 ( Ω ) 2 + ( K u , u ) + M p ( u L 2 ( Ω ) 2 ) + θ 2 u L 2 ( Ω ) 2 Ω u 2 ln u θ d x + λ i = 1 n d 0 ( t ) e i ( t ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 .
By adapting easily the condition (12) and relation (15), and with very minimal modifications, the results of the paper remain valid.
 
Theorem 2.1.
Assume that the hypotheses (AAA) and (H1)(H3) hold. Then, for the initial/initial history conditions  ( ϕ 0 , ϕ 1 , h 0 ) V 0 ( Ω , δ ( 0 ) ) , there exists a local solution u of problem (1), with  u S ( 0 , T )  and  w i C 0 ( [ 0 , T ] × [ 0 , 1 ] ; L 2 ( Ω ) ) C 1 ( [ 0 , T ) × [ 0 , 1 ] ; H 1 ( Ω ) ) , for  i = 1 , , n , such that E is a nonincreasing function and satisfies (for  t ( 0 , T ) )
d E d t ( t ) 2 R 1 ( λ ) u t L 2 ( Ω ) 2 2 R 2 ( λ ) i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 0
where  R 1 ( λ )  and  R 2 ( λ )  are given by
R 1 ( λ ) = n ( D 1 λ ) and R 2 ( λ ) = n ( λ D 2 ) min i = 1 , n ( inf t 0 ( 1 e i ( t ) ) ) .
Proof. 
To establish the existence result of a weak solution to problem (1), we apply the Faedo–Galerkin method, derive a priori estimates, and then pass to the limit in the approximate solutions using compactness arguments. We approximate Equation (1) by projecting them onto finite m dimensional subspaces, and then we take the limit in m. For this, let  ( v k ) k I N  be the orthogonal basis of  H 0 1 ( Ω ) , which is orthogonal in  L 2 ( Ω ) , and define  V m = s p a n { v 1 , . . . , v m } , for  m 1 . From  ( v k ) k = 1 , m , we can find a sequence  ( ψ k ) k = 1 , m  with  ψ k ( x , 0 ) = v k ( x ) , for  k = 1 , m , such that  ( ψ k ) k = 1 , m  is orthogonal in  L 2 ( Ω × ( 0 , 1 ) )  and then define  Ψ m = s p a n { ψ 1 , . . . , ψ m } , for  m 1  (as, e.g., in [38]). Let  u m 0 u m 1  be sequences of  V m  and  w m 0  be a sequence of  Ψ m  such that  u m 0 ϕ 0  in  H 0 1 ( Ω ) u m 1 ϕ 1  in  L 2 ( Ω )  and  w m 0 i w 0 i  in  L 2 ( Ω × ( 0 , 1 ) )  as  m . For each  m I N , we would like to define the Faedo–Galerkin approximation solution  ( u m , ( w m i ) i = 1 , n )  of the problem (1). Setting
u m ( x , t ) = k = 1 m h k m ( t ) v k ( x ) and w m i ( x , t ; η ) = k = 1 m g k m i ( t ) ψ k ( x , η ) ,
where  ( h k m ) k = 1 , m  and  ( g k m i ) k = 1 , m  are unknown functions and replacing  ( u , ( w i ) i = 1 , n )  by  ( u m , ( w m i ) i = 1 , n )  in (1), we obtain a.e.  t ( 0 , T ) , the system of Galerkin equations (for all  ( v , ψ ) V m × Ψ m  and  i = 1 , , n )
( 2 u m t 2 , v ) + ( K u m , v ) + ( d 0 ( t ) u m t , v ) + i = 1 n ( d i ( t ) w m i ( x , t ; 1 ) , v ) + ( M ( u m L 2 ( Ω ) 2 ) u m , v ) = ( u m ln u m θ , v ) , 0 1 ( e i ( t ) w m i t ( . , t ; η ) , ψ ( η ) ) d η + 0 1 ( ( 1 η e i ( t ) ) w m i η ( . , t ; η ) , ψ ( η ) ) d η = 0 , u m ( x , 0 ) = u m 0 ( x ) , u m t ( x , 0 ) = u m 1 ( x ) , for x Ω w m i ( x , 0 ; η ) = w m 0 i ( x ; η ) for ( x , η ) Ω × ( 0 , 1 ) .
By virtue of the standard theory of ordinary differential equations, the problem (16) has a local solution that is extended to a maximal interval  [ 0 , T m )  (with  0 < T m < T , for any given  T > 0 ). The following estimate will give the local solution being extended to the whole interval  [ 0 , T ) .
Replacing  ( v , ψ )  by  ( u m t , λ w m i )  in (16) and using the relation
( u m ln u m θ , u m t ) = Ω θ 2 u m ln u m 2 u m t d x = d d t ( θ 4 Ω ( u m ) 2 ln u m 2 d x θ 4 u m L 2 ( Ω ) 2 ) , = d d t ( 1 2 Ω ( u m ) 2 ln u m θ d x θ 4 u m L 2 ( Ω ) 2 ) ,
we obtain
1 2 d d t ( u m t L 2 ( Ω ) 2 + ( K u m , u m ) + M p ( u m L 2 ( Ω ) 2 ) + θ 2 u m L 2 ( Ω ) 2 Ω ( u m ) 2 ln u m θ d x ) = d 0 ( t ) u m t L 2 ( Ω ) 2 i = 1 n d i ( t ) ( w m i ( . , t ; 1 ) , u m t ) , 1 2 d d t ( λ e i ( t ) w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 ) = 1 2 λ e i ( t ) w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 1 2 λ e i ( t ) w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 λ ( ( 1 e i ( t ) ) w m i ( . , t ; 1 ) L 2 ( Ω ) 2 w m i ( . , t ; 0 ) L 2 ( Ω ) 2 ) , = λ ( 1 e i ( t ) ) w m i ( . , t ; 1 ) L 2 ( Ω ) 2 w m i ( . , t ; 0 ) L 2 ( Ω ) 2 , u m ( x , 0 ) = u m 0 ( x ) , u m t ( x , 0 ) = u m 1 ( x ) , for x Ω w m i ( x , 0 ; η ) = w m 0 i ( x ; η ) for ( x , η ) Ω × ( 0 , 1 ) .
Let us introduce the energy  E m  of the solution  u m  of problem (16) (as (11))
E m ( t ) = u m t L 2 ( Ω ) 2 + ( K u m , u m ) + M p ( u m L 2 ( Ω ) 2 ) + θ 2 u m L 2 ( Ω ) 2 Ω ( u m ) 2 ln u m θ d x + λ i = 1 n e i ( t ) w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 .
Then, (since  w m i ( . , . ; 0 ) = u m t )
1 2 d E m d t ( t ) = d 0 ( t ) u m t L 2 ( Ω ) 2 λ i = 1 n u m t L 2 ( Ω ) 2 i = 1 n d i ( t ) ( w m i ( . , t ; 1 ) , u m t ) λ i = 1 n ( 1 e i ( t ) ) w m i ( . , t ; 1 ) L 2 ( Ω ) 2 .
By Young’s inequality, we obtain
1 2 d E m d t ( t ) d 0 ( t ) u m t L 2 ( Ω ) 2 λ i = 1 n u m t L 2 ( Ω ) 2 + 1 2 i = 1 n d i ( t ) u m t L 2 ( Ω ) 2 + 1 2 i = 1 n d i ( t ) w m i ( . , t ; 1 ) L 2 ( Ω ) 2 λ i = 1 n ( 1 e i ( t ) ) w m i ( . , t ; 1 ) L 2 ( Ω ) 2
and then, according to the result (12)
1 2 d E m d t ( t ) ( d 0 ( t ) i = 1 n ( λ + d i ( t ) 2 ) ) u m t L 2 ( Ω ) 2 i = 1 n ( λ ( 1 e i ( t ) ) d i ( t ) 2 ) w m i ( . , t ; 1 ) L 2 ( Ω ) 2 R 1 ( λ ) u m t L 2 ( Ω ) 2 R 2 ( λ ) i = 1 n w m i ( . , t ; 1 ) L 2 ( Ω ) 2 0 ,
where  R 1 ( λ )  and  R 2 ( λ )  are given by (15).
From (H2), (4), (7) and (21), we can deduce that
u m t L 2 ( Ω ) 2 + ( ν θ α 2 2 π ) u m L 2 ( Ω ) 2 + θ ( 1 + N ( 1 + ln α ) ) + ϑ 0 2 u m L 2 ( Ω ) 2 + 2 R 1 ( λ ) 0 t u m t L 2 ( Ω ) 2 d s + 2 R 2 ( λ ) i = 1 n 0 t w m i ( . , s ; 1 ) L 2 ( Ω ) 2 d s + i = 1 n w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 E m ( 0 ) + θ 2 u m L 2 ( Ω ) 2 ln u m L 2 ( Ω ) 2 .
According to assumption (AAA), we have that
1 K 1 = ν θ α 2 2 π > 0 and 1 K 2 = θ ( 1 + N ( 1 + ln α ) ) + ϑ 0 2 > 0
and then
u m t L 2 ( Ω ) 2 + 1 K 1 u m L 2 ( Ω ) 2 + 1 K 2 u m L 2 ( Ω ) 2 + 2 R 1 ( λ ) 0 t u m t L 2 ( Ω ) 2 d s + 2 R 2 ( λ ) i = 1 n 0 t w m i ( . , s ; 1 ) L 2 ( Ω ) 2 d s + i = 1 n w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 E m ( 0 ) + θ 2 u m L 2 ( Ω ) 2 ln u m L 2 ( Ω ) 2 .
Moreover, we have, according to the expression of  E m ( 0 ) , for large m E m ( 0 ) C I ( ϕ 0 , ϕ 1 , h 0 ) , where  C I  is a positive constant depending on the initial/initial history data  ( ϕ 0 , ϕ 1 , h 0 ) , and then
u m t L 2 ( Ω ) 2 + 1 K 1 u m L 2 ( Ω ) 2 + 1 K 2 u m L 2 ( Ω ) 2 + 2 R 1 ( λ ) 0 t u m t L 2 ( Ω ) 2 d s + 2 R 2 ( λ ) i = 1 n 0 t w m i ( . , s ; 1 ) L 2 ( Ω ) 2 d s + i = 1 n w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 C I + θ 2 u m L 2 ( Ω ) 2 ln u m L 2 ( Ω ) 2 .
Since  u m ( x , t ) = u m ( x , 0 ) + 0 t u m t ( x , s ) d s , then, for large m, (according to (24))
u m ( . , t ) L 2 ( Ω ) 2 2 u m ( 0 ) L 2 ( Ω ) 2 + C 1 0 t u m t L 2 ( Ω ) 2 d s C 0 + C 2 0 t u m L 2 ( Ω ) 2 ln u m L 2 ( Ω ) 2 d s
According to logarithmic Gronwall inequality, we can deduce the estimate (for all t)
u m ( . , t ) L 2 ( Ω ) 2 ( 1 + C 0 ) e x p ( C 2 t ) C T
and then from inequality (24) follows (for a.e.  t ( 0 , T ) )
u m t L 2 ( Ω ) 2 + u m L 2 ( Ω ) 2 + u m L 2 ( Ω ) 2 + 0 t u m t L 2 ( Ω ) 2 d s + i = 1 n 0 t w m i ( . , s ; 1 ) L 2 ( Ω ) 2 d s + i = 1 n w m i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 C .
Consequently,
( u m ) is bounded in L ( 0 , T ; H 0 1 ( Ω ) ) × L ( 0 , T ; L 2 ( Ω ) ) , ( u m t ) is bounded in L ( 0 , T ; L 2 ( Ω ) ) × L 2 ( 0 , T ; L 2 ( Ω ) ) , ( w m i ) is bounded in L ( 0 , T ; L 2 ( Ω × ( 0 , 1 ) ) ) , ( w m i ( . , . ; 1 ) ) is bounded in L 2 ( 0 , T ; L 2 ( Ω ) ) .
This result makes it possible to extract from  ( u m , w m i )  a subsequence also denoted by  ( u m , w m i )  and such that (by using Aubin–Lions compactness lemma, see [48] Theorem 5.1, p. 58)
u m u weakly star in L ( 0 , T ; H 0 1 ( Ω ) ) , u m t u t weakly star in L ( 0 , T ; L 2 ( Ω ) ) , u m t u t weakly in L 2 ( 0 , T ; L 2 ( Ω ) ) , w m i w i weakly star in L ( 0 , T ; L 2 ( Ω × ( 0 , 1 ) ) ) , w m i ( . , . ; 1 ) w i ( . , . ; 1 ) weakly in L 2 ( Q ) , u m u strongly in L 2 ( Q ) ( and so u m u a . e . in Q ) , ( ) , u m ln u m θ u ln u θ a . e . in Q ( from the continuity of y y l n y θ and ( ) ) , M ( u m L 2 ( Ω ) 2 ) M ( u L 2 ( Ω ) 2 ) strongly in L 2 ( 0 , T ) ( from the continuity of M and ( ) ) .
Moreover, since  sup 0 < r 1 r ln r = 1 e , we can deduce, for all  q > 2  if  N 2  and  2 < q 2 N N 2  if  N 3 , that
u m ln u m θ L 2 ( Ω ) θ 2 e ( 1 2 m e s ( Ω ) + 1 q 2 u m L q ( Ω ) q ) C ( 1 + u m L 2 ( Ω ) q ) .
Consequently, from (30) we obtain (according to the last result of (29) and Lebesgue dominated convergence arguments)
u m ln u m θ u ln u θ weakly in L 2 ( Q )
and (according to (29) and (16)) the boundedness of  2 u m t 2 in L ( 0 , T ; H 1 ( Ω ) )  and
2 u m t 2 2 u t 2 weakly star in L ( 0 , T ; H 1 ( Ω ) ) .
By using (29)–(32) and density properties of spaces spanned, respectively, by  ( v k ) k  and  ( ψ k ) k , we can pass to the limit ( m + ) in a standard way in (16). So we omit the details. The limit  ( u ; w i ) , for  i = 1 , n , then satisfies the following system (for all  ( v , ψ ) H 0 1 ( Ω ) × L 2 ( Ω × ( 0 , 1 ) ) )
< 2 u t 2 , v > H 1 , H 0 1 + ( K u , v ) + ( d 0 ( t ) u t , v ) + i = 1 n ( d i ( t ) w i ( x , t ; 1 ) , v ) + ( M ( u L 2 ( Ω ) 2 ) u , v ) = ( u ln u θ , v ) , 0 1 ( e i ( t ) w i t ( . , t ; η ) , ψ ( . , η ) ) d η + 0 1 ( ( 1 η e i ( t ) ) w i η ( . , t ; η ) , ψ ( . , η ) ) d η = 0 , u ( x , 0 ) = ϕ 0 ( x ) , u m t ( x , 0 ) = ϕ 1 ( x ) , for x Ω w i ( x , 0 ; η ) = w 0 i ( x ; η ) , for ( x , η ) Ω × ( 0 , 1 ) .
Finally, by a similar argument as to show (21), we can deduce
d E d t ( t ) 2 R 1 ( λ ) u t L 2 ( Ω ) 2 2 R 2 ( λ ) i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 0
where  R 1 ( λ )  and  R 2 ( λ )  are given in (15). Then, E is a nonincreasing function. This completes the proof. □
 
Remark 2.2.
If  T m a x = sup { T > 0 : the solution u existing on [ 0 , T ] }  is the maximal existence time of the weak solution to problem (1), and then if  T m a x = + , we say that the solution u is global and if  T m a x < + , the solution u blows up and  T m a x  is the blow-up time.

3. Global Existence and Energy Decay Estimate

In this section, we prove the global existence and energy decay rate of the solutions to problem (10).
We started by introducing the functions  I  and  J  such that (for  v H 0 1 ( Ω ) )
I ( v ) = ν v L 2 ( Ω ) 2 + ϑ 0 2 v L 2 ( Ω ) 2 θ Ω v 2 ln v d x and J ( v ) = I ( v ) + θ 2 v L 2 ( Ω ) 2 .
For  r > 0  and  v H 0 1 ( Ω ) , we obtain
J ( r v ) = I ( r v ) + r 2 θ 2 v L 2 ( Ω ) 2 = r 2 ( I ( v ) + θ ( 1 2 ln ( r ) ) 2 v L 2 ( Ω ) 2 ) .
Then
r ( J ( r v ) ) = 2 r ( I ( v ) + θ ( 1 2 ln ( r ) ) 2 v L 2 ( Ω ) 2 ) r θ v L 2 ( Ω ) 2 = 2 r ( I ( v ) θ ln ( r ) v L 2 ( Ω ) 2 ) .
Therefore, for an arbitrary element  v H 0 1 ( Ω ) / { 0 } , we have  r r ( J ( r v ) ) = 2 I ( r v ) lim r 0 + J ( r v ) = 0 lim r + J ( r v ) = , and  r J ( r v ) > 0  if  0 < r < r ¯ ( v )  and  r J ( r v ) < 0  if  r > r ¯ ( v ) , where the value  r ¯ ( v ) = exp ( I ( v ) θ v L 2 ( Ω ) 2 )  is the unique solution of  r J ( r v ) = 0 . Furthermore, we have that  2 r 2 J ( r v ) | r = r ¯ = 2 θ v L 2 ( Ω ) 2 < 0 arg max r > 0 J ( r v ) = r ¯ ( v )  and  I ( r ¯ ( v ) v ) = 0 .
Associated with  J , we have the well-known Nehari Manifold (the set of all nontrivial stationary solutions to the problem (10))
N = { v H 0 1 ( Ω ) / { 0 } : ( r J ( r v ) ) | r = 1 = 0 } .
Equivalently,
N = { v H 0 1 ( Ω ) / { 0 } : I ( v ) = 0 } .
We define the potential well depth d (also known as mountain pass level), as in the Mountain Pass theorem due to Ambrosetti and Rabinowitz [49], by
d = inf v H 0 1 ( Ω ) / { 0 } sup r > 0 J ( r v ) .
As Payne and Sattinger noted in [50], the potential well depth d can be also characterized as
d = inf v N J ( v ) .
Moreover, the weak solution of (10) blows up when  I ( u ( t ) ) < 0 . If  θ < 0 , then  I ( u ( t ) ) 0 , for all t and then the weak solution of (10) is global. Consequently, in the sequel, we assume that  θ > 0 .
Now, we can introduce the following spaces
N + = { v H 0 1 ( Ω ) / { 0 } : I ( v ) > 0 } { 0 } , N = { v H 0 1 ( Ω ) / { 0 } : I ( v ) < 0 } , W 1 = { v N + : J ( v ) < d } ( stable set ) , W 2 = { v N : J ( v ) < d } ( unstable set ) .
The space  W 1  is corresponding to the set of stability for the problem (10).
Theorem 3.1.
Assume that the hypotheses (AAA)(H2) hold, and that the initial condition  ( ϕ 0 , ϕ 1 ) H 0 1 ( Ω ) × L 2 ( Ω )  and initial history condition  h 0 L 2 ( Q 0 ) . Let u be a local weak solution to (10) and  E a [ 0 , d ] . If  ϕ 0 W 1  and  E ( 0 ) < E a , then u is a global solution.
 
Proof. 
First, from Lemma A3, we have that  u ( t ) W 1 , for every  t [ 0 , T m a x ) . Prove now that  T m a x = . For this purpose, it is sufficient to prove the boundedness (in time) of  u t L 2 ( Ω ) 2 + u L 2 ( Ω ) 2 . Since  0 < E ( 0 ) < E a d , we have for all  t [ 0 , T m a x )  (from the definition of  E ( t )  and the positivity of  I ( u ( t ) ) )
u t L 2 ( Ω ) 2 E ( t ) E ( 0 ) < E a < d , θ 2 u L 2 ( Ω ) 2 E ( t ) E ( 0 ) < E a < d .
Moreover, from Lemma 1.1, we can deduce
ν u L 2 ( Ω ) 2 = I ( u ( t ) ) ϑ 0 2 u L 2 ( Ω ) 2 + θ Ω u 2 ln u d x I ( u ( t ) ) + θ 2 u L 2 ( Ω ) 2 ( ln ( u L 2 ( Ω ) 2 ) N ( 1 + ln α ) ϑ 0 θ ) + θ α 2 2 π u L 2 ( Ω ) 2 E ( t ) + θ 2 u L 2 ( Ω ) 2 ( ln ( u L 2 ( Ω ) 2 ) N ( 1 + ln α ) ϑ 0 θ ) + θ α 2 2 π u L 2 ( Ω ) 2
and then
( ν θ α 2 2 π ) u L 2 ( Ω ) 2 E ( t ) + θ 2 u L 2 ( Ω ) 2 h α , θ , ϑ 0 ( u L 2 ( Ω ) ) ,
where
h α , θ , ϑ 0 ( u L 2 ( Ω ) ) = ln ( u L 2 ( Ω ) 2 ) N ( 1 + ln α ) ϑ 0 θ .
From (38), (A1) in Appendix A and (AAA), we can deduce that
h α , θ , ϑ 0 ( u L 2 ( Ω ) ) ln ( 2 d θ ) N ( 1 + ln α ) ϑ 0 θ , ln ( 2 d θ ) N ( 1 + ln α ) ϑ 0 θ ln ( 2 d I θ ) N ( 1 + ln α ) ϑ 0 θ 2 ( ν θ α 2 2 π ) C ˜ Ω 2 θ + N ( 1 + ln α ) N ( 1 + ln α ) ϑ 0 θ + ϑ 0 θ = 2 ( ν θ α 2 2 π ) C ˜ Ω 2 θ > 0
and then (since  E ( t ) < d )
u L 2 ( Ω ) 2 < d ( 1 + ln ( 2 d / θ ) N ( 1 + ln α ) ϑ 0 θ ) ( ν θ α 2 2 π ) 1 .
Thus, according to (38) and (43), we can conclude that the solution u is global (by the continue principle). This completes the proof. □
 
Remark 3.1.
1.
From (34), we can deduce the following boundedness results (since  0 < E ( t ) < d  for  t 0 )
0 t u t L 2 ( Ω ) 2 d s d 2 R 1 ( λ ) , 0 t w i ( s ; 1 ) L 2 ( Ω ) 2 d s d 2 R 2 ( λ ) , for i = 1 , , n .
2.
If we replace the conditions  ϕ 0 W 1  and  E ( 0 ) < E a , by the existence of a real number  t 1 [ 0 ; T m a x )  such that  u ( t 1 ) W 1  and  E ( t 1 ) < E a , the result of Theorem remains valid.
3.
From (A1) in Appendix A, (38) and (AAA), we can deduce the following relations.
(i)
If  E a d I
h α , θ , ϑ 0 ( u L 2 ( Ω ) ) ln ( 2 E a θ ) N ( 1 + ln α ) ϑ 0 θ , ln ( 2 E a θ ) N ( 1 + ln α ) ϑ 0 θ ln ( 2 d I θ ) N ( 1 + ln α ) ϑ 0 θ = 2 ( ν θ α 2 2 π ) C ˜ Ω 2 θ > 0 ,
(ii)
if  E a d I
h α , θ , ϑ 0 ( u L 2 ( Ω ) ) ln ( 2 E a θ ) N ( 1 + ln α ) ϑ 0 θ ln ( 2 d I θ ) N ( 1 + ln α ) ϑ 0 θ = 2 ( ν θ α 2 2 π ) C ˜ Ω 2 θ ,
(iii)
if  E a d A
h α , θ , ϑ 0 ( u L 2 ( Ω ) ) ln ( 2 E a θ ) N ( 1 + ln α ) ϑ 0 θ ln ( 2 d A θ ) N ( 1 + ln α ) ϑ 0 θ = 0 .

4. Asymptotic Behavior

In this section, we prove the exponential decay of solution of problem (34). We assume that the operator  M  satisfies the following condition (for all  v 0 )
M ( v ) v M p ( v ) 0
and we use the Nakao’s Lemma [51].
Theorem 4.1.
Let the assumptions of Theorem 3.1 hold. Assume that  E a d I  and  ν > θ α 2 π . Then, if  E a d A  or  2 ( ν θ α 2 2 π ) C ˜ Ω 2 ϑ 0 , there exist positive constants  C s  and δ such that the energy E associated to problem (10) satisfies
0 < E ( t ) C s e δ t , t 0 .
Proof. 
Let  v = u  and  ψ = exp ( η e i ) w i  in (33), we can deduce that
u t L 2 ( Ω ) 2 + d d t ( u t , u ) + ( K u , u ) + M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 + ( d 0 ( t ) u t , u ) + i = 1 n ( d i ( t ) w i ( . , t ; 1 ) , u ) = ( u ln u θ , u ) , 0 1 e i ( t ) exp ( η e i ( t ) ) t w i ( . , t ; η ) L 2 ( Ω ) 2 d η + 0 1 exp ( η e i ( t ) ) ( 1 η e i ( t ) ) η w i ( . , t ; η ) L 2 ( Ω ) 2 d η = 0 ,
then (since  w i ( . ; η = 0 ) = u t ),
d d t ( u t , u ) = u t L 2 ( Ω ) 2 ( K u , u ) M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 d 0 ( t ) ( u t , u ) i = 1 n ( d i ( t ) w i ( . , t ; 1 ) , u ) + Ω u 2 ln u θ d x , t ( e i ( t ) exp ( η e i ( t ) ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 ) 0 1 exp ( η e i ) e i ( 1 η e i ) ) w i ( . t , ; η ) L 2 ( Ω ) 2 ( t ; η ) d η + 0 1 exp ( η e i ( t ) ) ( e i + e i ( t ) ( 1 η e i ( t ) ) ) w i ( . , t ; η ) L 2 ( Ω ) 2 d η + e x p ( e i ( t ) ) ( 1 e i ( t ) ) w i ( . , t ; 1 ) L 2 ( Ω ) 2 u t L 2 ( Ω ) 2 = 0 .
Introduce now the following functions:
Φ ( t ) = ( u t , u ) and Ψ ( t ) = λ θ i = 1 n 0 1 e i ( t ) exp ( η e i ) w i L 2 ( Ω ) 2 ( t ; η ) d η .
Then,
Φ ( t ) = u t L 2 ( Ω ) 2 ( K u , u ) M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 d 0 ( t ) Φ ( t ) i = 1 n ( d i ( t ) w i ( x , t ; 1 ) , u ) + Ω u 2 ln u θ d x ( 1 + ( C Ω d 0 m a x ) 2 ν ) u t L 2 ( Ω ) 2 1 2 ( K u , u ) M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 + ( C Ω D ) 2 ν i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 + Ω u 2 ln u θ d x
and
Ψ ( t ) = λ θ i = 1 n 0 1 exp ( η e i ( t ) ) e i ( t ) w i ( . , t ; η ) L 2 ( Ω ) 2 ( t ; η ) d η λ θ i = 1 n exp ( e i ( t ) ) ( 1 e i ( t ) ) w i ( . , t ; 1 ) L 2 ( Ω ) 2 + λ θ u t L 2 ( Ω ) 2 λ θ i = 1 n e i ( t ) exp ( e i ( t ) ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 λ P θ i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 + λ θ u t L 2 ( Ω ) 2
where  P = min i ( ( 1 p i m a x ) exp ( e i m a x ) ) D = max i ( d i m a x ) . Consequently (according to (3) and (H1))
( Φ + Ψ ) ( t ) ( 1 + ( C Ω d 0 m a x ) 2 ν + λ θ ) u t L 2 ( Ω ) 2 1 2 ( K u , u ) M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 ( ( C Ω D ) 2 ν + λ P θ ) i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 + Ω u 2 ln u θ d x λ θ i = 1 n e i ( t ) exp ( e i ( t ) ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 .
Then, for all  ϵ , ξ > 0 , we have (according to (3) and (H1))
( E ( t ) + ϵ ( Φ + Ψ ) ) ( t ) + λ ξ θ E ( t ) ( 2 R 1 ( λ ) ϵ ( 1 + ( C Ω d 0 m a x ) 2 ν + λ θ ) λ ξ θ ) u t L 2 ( Ω ) 2 + ( C Ω 2 λ ξ 2 θ θ ν ( ϵ 2 λ ξ θ ) ) u L 2 ( Ω ) 2 ϵ M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 + λ ξ θ M p ( u L 2 ( Ω ) 2 ) ( 2 R 2 ( λ ) + ϵ ( ( C Ω D ) 2 ν + λ P θ ) ) i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 + ( ϵ λ ξ θ ) Ω u 2 ln u θ d x λ θ i = 1 n e i ( t ) ( ϵ exp ( e i ( t ) ) ξ λ ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 .
So, from the logarithmic Sobolev inequality (for  ϵ 2 > ρ = λ ξ θ > 0 )
( E ( t ) + ϵ ( Φ + Ψ ) ) ( t ) + ρ E ( t ) ( 2 R 1 ( λ ) ϵ ( 1 + ( C Ω d 0 m a x ) 2 ν + λ θ ) ρ ) u t L 2 ( Ω ) 2 ( C Ω 2 θ ρ 2 ( ϵ ρ ) θ α 2 2 π + ν ( ϵ 2 ρ ) ) u L 2 ( Ω ) 2 ( ϵ ρ ) M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 ρ ( M p ( u L 2 ( Ω ) 2 ) + M ( u L 2 ( Ω ) 2 ) u L 2 ( Ω ) 2 ) ( 2 R 2 ( λ ) + ϵ ( ( C Ω D ) 2 ν + λ P θ ) ) i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 + ( ϵ ρ ) θ ( ln u L 2 ( Ω ) 2 N ( 1 + ln α ) ) 1 2 u L 2 ( Ω ) 2 λ θ i = 1 n e i ( t ) ( ϵ exp ( e i ( t ) ) ρ θ ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 .
According to (H2) and (48), we can deduce that (according to the expression (41))
( E ( t ) + ϵ ( Φ + Ψ ) ) ( t ) + ρ E ( t ) ( 2 R 1 ( λ ) ϵ ( 1 + ( C Ω d 0 m a x ) 2 ν + λ θ ) ρ ) u t L 2 ( Ω ) 2 ( ρ ( C Ω 2 θ 2 + ( ν θ α 2 2 π ) ) + ϵ 2 ( ν θ α 2 π ) ) u L 2 ( Ω ) 2 ( 2 R 2 ( λ ) + ϵ ( ( C Ω D ) 2 ν + λ P θ ) ) i = 1 n w i ( . , t ; 1 ) L 2 ( Ω ) 2 ( t ; 1 ) + ( ϵ ρ ) θ 2 u L 2 ( Ω ) 2 ( h α , θ , ϑ 0 ( u L 2 ( Ω ) ) ϑ 0 θ ) λ θ i = 1 n e i ( t ) exp ( e i ( t ) ) ( ϵ ρ θ e ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 ,
where  e = max i ( e x p ( e i m a x ) ) .
By taking  ϵ  and  ρ ϵ  sufficiently small, we can deduce ( ν > θ α 2 π )
( E ( t ) + ϵ ( Φ + Ψ ) ) ( t ) + ρ E ( t ) ( ϵ ρ ) θ 2 u L 2 ( Ω ) 2 ( h α , θ , ϑ 0 ( u L 2 ( Ω ) ) ϑ 0 θ ) .
Since  0 < E ( 0 ) < E a d I , we have that from (46)–(47) that (since  E a d A  or  2 ( ν θ α 2 2 π ) C ˜ Ω 2 ϑ 0 )
( E ( t ) + ϵ ( Φ + Ψ ) ) ( t ) ρ E ( t ) .
Prove now the following energy equivalence (for  ϵ  sufficiently small)
α 1 E ( t ) E ( t ) + ϵ ( Φ + Ψ ) ( t ) α 2 E ( t ) ,
where  α i > 0 , for  i = 1 , 2 . Since
( Φ + Ψ ) ( t ) ( u t , u ) + λ θ i = 1 n e i ( t ) 0 1 exp ( η e i ) w i ( . , t ; η ) L 2 ( Ω ) 2 d η 2 θ ( u t L 2 ( Ω ) 2 + θ 2 u L 2 ( Ω ) 2 + 2 i = 1 n λ e i ( t ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 ) 2 θ ( u t L 2 ( Ω ) 2 + J ( u ) ( t ) + i = 1 n λ e i ( t ) w i ( . , t ; . ) L 2 ( Ω × ( 0 , 1 ) ) 2 ) 2 θ E ( t )
then
α 1 E ( t ) E ( t ) + ϵ ( Φ + Ψ ) ( t ) α 2 E ( t ) ,
where  α 1 = 1 ϵ 2 θ > 0  and  α 2 = 1 + ϵ 2 θ > 0 . From (59) and (58), we can obtain that  E ( t ) + ϵ ( Φ + Ψ ) ( t ) β 1 e x p ( ρ t )  (where  β 1 > 0  is constant depending on  E ( 0 ) + ϵ ( Φ + Ψ ) ( 0 ) α 1  and  α 2 ), and by using again (59), we can deduce the result of the theorem. □

5. Conclusions

In this work, we have studied the existence and exponential stability of global solutions to nonlinear logarithmic Klein–Gordon type equations with multiple-time varying delays and a nonlocal term in a bounded domain. The logarithmic Klein–Gordon equation is the relativistic version of the logarithmic Schrodinger equation. Such logarithmic nonlinearity effects often arise in various areas of physical sciences and engineering. The introduction of retarded arguments is to reflect the different after-effects. Various time-varying delay configurations occur naturally in various areas of physics, biologics, and engineering. This natural phenomenon is due to the fact that the instantaneous rate of change of such systems does not only depend on their current time but rather on their previous history as well. Moreover, in many realistic application fields, introducing time delays into mathematical modeling has long proven to be unavoidable for correctly representing the behavior of real-world systems. It is also known that time delay is a non-negligible constraint in the process and may induce complex behaviors in the dynamical system, e.g., instability, oscillations, chaos, and poor performances. Therefore, these behaviors and aspects, by taking into account different sources of delays, motivate the study of multiple time-varying delay effects on properties of dynamical systems.
In this respect, in order to take into account the influence of different sources of time delays in the velocity of signal transmission, nonlinear Klein–Gordon wave type models with logarithmic nonlinearity are modified by incorporating multiple time delays and a nonlocal operator. The presence of nonlocal term and multiple time delays in a system leads to a more complex analysis.
The proposed strategy consists in controlling the instabilities by imposing some suitable conditions involving different functions and parameters representing the multiple time-varying delays. After obtaining the local existence result of the solutions by using Faedo–Galerkin’s method and logarithmic Sobolev inequality, the global existence is derived. To show the exponential stability result under appropriate conditions, the potential well and perturbed energy methods are applied.
A future objective is to simulate and validate numerically the developed theoretical results. These studies will be the subject of a forthcoming paper. It would also be interesting to investigate the blow-up behavior of solutions of the considered problem. Moreover, the developed analysis in this work can be further applied to the studied model but
  • With dynamical boundary conditions on regular non-cylindrical domains
  • With switching time delays
  • Or with Kirchhoff–Carrier type operators, i.e., by replacing the term  K ( x )  in (1) by the operator  K ( x , u L 2 ( Ω ) 2 ) .

Funding

This research received no external funding.

Data Availability Statement

No data was used for this manuscript.

Acknowledgments

The author thanks the referees for their valuable comments and suggestions to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this annex we give some lemmas used in the document.
Lemma A1.
The potential depth d satisfies
d d I = θ 2 exp ( 2 ( ν θ α 2 2 π ) C ˜ Ω 2 + ϑ 0 + N θ ( 1 + ln α ) θ ) = θ 2 exp ( 2 ( ν θ α 2 2 π ) C ˜ Ω 2 + ϑ 0 θ ) e N α N , d I d A = θ 2 exp ( ϑ 0 θ ) e N α N ,
where  C ˜ Ω = C Ω 1 C Ω  is the Poincaré constant and  α > 0  is defined in assumption (AAA). Moreover, we have
ρ 0 = θ 2 exp ( 2 ( ν θ α 1 2 2 π ) C ˜ Ω 2 + ϑ 0 θ ) e N α 0 N < d I < ρ 1 = θ 2 exp ( 2 ( ν θ α 0 2 2 π ) C ˜ Ω 2 + ϑ 0 θ ) e N α 1 N , ρ 0 a = θ 2 exp ( ϑ 0 θ ) e N α 0 N < d A < ρ 1 a = θ 2 exp ( ϑ 0 θ ) e N α 1 N .
Proof. 
According to assumption (AAA), Lemma 1.1 and relation (3), we can deduce that (for all  v H 0 1 ( Ω ) )
0 = I ( r ¯ ( v ) v ) ( ( ν θ α 2 2 π ) C ˜ Ω 2 + N θ 2 ( 1 + ln α ) + ϑ 0 2 θ 2 ln r ¯ ( v ) v L 2 ( Ω ) 2 ) r ¯ ( v ) v L 2 ( Ω ) 2
and then
r ¯ ( v ) v L 2 ( Ω ) 2 exp ( 2 ( ν θ α 2 2 π ) C ˜ Ω 2 + N θ ( 1 + ln α ) + ϑ 0 θ ) .
Consequently,
d J ( r ¯ ( v ) v ) = θ 2 r ¯ ( v ) v L 2 ( Ω ) 2 d I = θ 2 exp ( 2 ( ν θ α 2 2 π ) C ˜ Ω 2 + N θ ( 1 + ln α ) + ϑ 0 θ ) d A = θ 2 exp ( ϑ 0 θ ) e N α N .
According to (6), we obtain easily the relations (A2). This completes the proof. □
 
Lemma A2.
Let  v 0  be in  H 0 1 ( Ω ) , we have
1
If  0 < v L 2 ( Ω ) 2 R N ( λ m i n , ν , ϑ 0 )  then  I ( v ) 0 ,
where  R N ( λ m i n , ν , ϑ 0 ) = λ m i n e N + ϑ 0 ( 2 π ν ) N 2  and  λ m i n  is the first eigenvalue of  Δ  with homogeneous Dirichlet boundary conditions, that is
λ m i n = inf v H 0 1 ( Ω ) { 0 } v L 2 ( Ω ) 2 v L 2 ( Ω ) 2
(and then  λ m i n 1 = C Ω  is the Poincaré constant)
2
If  I ( v ) < 0 , then  v L 2 ( Ω ) 2 > R N ( λ m i n , ν , ϑ 0 ) .
Proof. 
From Lemma 1.1, for any positive constant  α , we have
I ( v ) = ν v L 2 ( Ω ) 2 + ϑ 0 2 v L 2 ( Ω ) 2 θ Ω v 2 ln v d x ( ν α 2 2 π ) v L 2 ( Ω ) 2 + 1 2 v L 2 ( Ω ) 2 ( θ N ( 1 + ln α ) + ϑ 0 ) θ ln v L 2 ( Ω ) 2 ) .
For  α  such that  ν θ α 2 2 π = 0  (i.e.,  α = 2 π ν θ ), we can deduce that
I ( v ) 1 2 v L 2 ( Ω ) 2 ( θ N 2 ( 2 + ln ( 2 π ν ) ) + ϑ 0 ) θ ln v L 2 ( Ω ) 2 ) .
We can deduce that:
-
if  0 < v L 2 ( Ω ) 2 R N ( λ m i n , ν , ϑ 0 ) , then  I ( v ) 0 ,
where  R N ( λ m i n , ν , ϑ 0 ) = λ m i n e N + ϑ 0 ( 2 π ν ) N 2 ,
-
if  I ( v ) < 0  then  v L 2 ( Ω ) 2 > R N ( λ m i n , ν , ϑ 0 ) .
This completes the proof. □
 
Lemma A3.
Assume that the hypotheses (AAA)(H2) and the initial/initial history conditions  ( ϕ 0 , ϕ 1 , h 0 ) H 0 1 ( Ω ) × L 2 ( Ω ) × L 2 ( Q 0 )  hold. We have, for  E a ] 0 , d ] ,
(i) 
if there exists a real number  t 1 [ 0 , T m a x )  such that  u ( t 1 ) W 2  and  E ( t 1 ) < E a d , then  u ( t ) W 2 , for every  t [ t 1 , T m a x ) . Moreover, (for every  t [ t 1 , T m a x ) )
θ 2 u ( t ) L 2 ( Ω ) 2 > E a ,
(ii) 
if there exists a real number  t 1 [ 0 , T m a x )  such that  u ( t 1 ) W 1  and  E ( t 1 ) < E a d , then  u ( t ) W 1 , for every  t [ t 1 , T m a x ) .  Moreover, (for every  t [ t 1 , T m a x ) )
θ 2 u ( t ) L 2 ( Ω ) 2 < E a ,
T m a x  is the maximal existence time of the weak solution to problem (1).
 
Proof. 
Without loss of generality, we assume that  t 1 = 0 .
(i)
If  u ( 0 ) W 2  and  E ( 0 ) < E a d , then the solution u satisfies  u ( t ) W 2  and  E ( t ) < d , for every  t [ 0 , T m a x ) . In fact, since E is a nonincreasing function, then  E ( t ) E ( 0 ) < d , for every  t [ 0 , T m a x ) . Assume that there exists  t 0  such that  u ( t 0 ) N , then from the definition of d and (11), we can deduce that  d J ( u ( t 0 ) ) E ( t 0 ) E ( 0 ) < d , which is impossible. Consequently,  u ( t ) W 2 , for every  t [ 0 , T m a x ) . Moreover, from the definition of d and the fact that  I ( u ( t ) ) < 0  (since  u ( t ) W 2 ), we can deduce that
E a d J ( r ¯ ( u ) u ) = exp ( 2 I ( u ) θ u L 2 ( Ω ) 2 ) θ 2 u L 2 ( Ω ) 2 < θ 2 u L 2 ( Ω ) 2 .
(ii)
If  u ( 0 ) W 1  and  E ( 0 ) < E a d , then the solution u satisfies  u ( t ) W 1  and  E ( t ) < d , for every  t [ 0 , T m a x ) . In fact, since E is a nonincreasing function, then  E ( t ) E ( 0 ) < E a d , for every  t [ 0 , T m a x ) . Suppose that there exists  t 0  such that  u ( t 0 ) N , then from the definition of d and (11), we can deduce that  d J ( u ( t 0 ) ) E ( t 0 ) E ( 0 ) < d , which is impossible. Consequently,  u ( t ) W 1 , for every  t [ 0 , T m a x ) . Moreover, from the definition of  E ( t )  and the positivity of  I ( u ( t ) ) , we obtain  θ 2 u ( t ) L 2 ( Ω ) 2 E ( t ) E ( 0 ) < E a . This completes the proof.

References

  1. Belmiloudi, A. Robust control problems associated with time-varying delay nonlinear parabolic equations. IMA J. Math. Control Inf. 2003, 20, 305–334. [Google Scholar] [CrossRef]
  2. Bentrcia, T.; Mennouni, A. On the solution behavior of a nonlinear time-fractional Klein–Gordon equation: Theoretical study and numerical validation. Commun. Nonlinear Sci. Numer. Simul. 2023, 125, 107384. [Google Scholar] [CrossRef]
  3. Kolkovska, N.; Dimova, M.; Kutev, N. Nonexistence of global solutions to Klein-Gordon equations with variable coefficients power-type nonlinearities. Open Math. 2023, 21, 20220584. [Google Scholar] [CrossRef]
  4. Killip, R.; Stovall, B.; Visan, M. Blow-up behaviour for the nonlinear Klein–Gordon equation. Math. Ann. 2014, 358, 289–350. [Google Scholar] [CrossRef]
  5. Xu, R.Z. Global existence, blow up and asymptotic behavior of solutions for nonlinear Klein-Gordon equation with dissipative term. Math. Methods Appl. Sci. 2010, 33, 831–844. [Google Scholar]
  6. De Martino, S.; Falanga, M.; Godano, C.; Lauro, G. Logarithmic Schrodinger-like equation as a model for magma transport. Europhys. Lett. 2003, 63, 472–475. [Google Scholar] [CrossRef]
  7. Krolikowski, W.; Edmundson, D.; Bang, O. Unified model for partially coherent solitons in logarithmically nonlinear media. Phys. Rev. E 2000, 61, 3122–3126. [Google Scholar] [CrossRef]
  8. Wazwaz, A.M. Gaussian solitary waves for the logarithmic-KdV and the logarithmic-KP equations. Phys. Scripta 2014, 89, 095206. [Google Scholar] [CrossRef]
  9. Cherfils, L.; Miranville, A.; Zelik, S. The Cahn-Hilliard equation with logarithmic potentials. Milan J. Math. 2011, 79, 561–596. [Google Scholar] [CrossRef]
  10. Barrow, J.D.; Parsons, P. Inflationary models with logarithmic potentials. Phys. Rev. D 1995, 52, 5576–5587. [Google Scholar] [CrossRef]
  11. Bialynicki-Birula, I.; Mycielski, J. Wave equations with logarithmic nonlinearities. Bull. Acad. Polon. Sci. 1975, 3, 461–466. [Google Scholar]
  12. Buljan, H.; Šiber, A.; Soljačić, M.; Schwartz, T.; Segev, M.; Christodoulides, D.N. Incoherent white light solitons in logarithmically saturable noninstantaneous nonlinear media. Phys. Rev. E 2003, 68, 036607. [Google Scholar] [CrossRef]
  13. Cherstvy, A.G.; Safdari, H.; Metzler, R. Anomalous diffusion, nonergodicity, and ageing for exponentially and logarithmically time-dependent diffusivity: Striking differences for massive versus massless particles. J. Phys. D Appl. Phys. 2021, 54, 195401. [Google Scholar] [CrossRef]
  14. Enqvist, K.; McDonald, J. Q-balls and baryogenesis in the MSSM. Phys. Lett. B 1998, 425, 309–321. [Google Scholar] [CrossRef]
  15. Rosen, G. Dilatation covariance and exact solutions in local relativistic field theories. Phys. Rev. 1969, 183, 1186. [Google Scholar] [CrossRef]
  16. Scott, T.C.; Zhang, X.D.; Mann, R.B.; Fee, G.J. Canonical reduction for dilatonic gravity in 3 + 1 dimensions. Phys. Rev. D 2016, 93, 084017. [Google Scholar] [CrossRef]
  17. Al-Gharabli, M.M.; Messaoudi, S. Existence and a general decay result for a plate equation with nonlinear damping and a logarithmic source term. J. Evol. Equ. 2018, 18, 105–125. [Google Scholar] [CrossRef]
  18. Bartkowski, K.; Górka, P. One-dimensional Klein-Gordon equation with logarithmic nonlinearities. J. Phys. A 2008, 41, 355201. [Google Scholar] [CrossRef]
  19. Cazenave, T.; Haraux, A. Équations d’évolution avec non linéarité logarithmique. Ann. Fac. Des Sci. Toulouse Série 1980, 2, 21–51. [Google Scholar] [CrossRef]
  20. Chen, H.; Luo, P.; Liu, G. Global solution and blow-up of a semilinear heat equation with logarithmic nonlinearity. J. Math. Anal. Appl. 2015, 422, 84–98. [Google Scholar] [CrossRef]
  21. Gomes, D.A.; Pimentel, E. Time-Dependent Mean-Field Games with Logarithmic Nonlinearities. SIAM J. Math. Anal. 2015, 47, 3798–3812. [Google Scholar] [CrossRef]
  22. Gorka, P. Logarithmic quantum mechanics: Existence of the ground state. Found. Phys. Lett. 2006, 19, 591–601. [Google Scholar] [CrossRef]
  23. Gorka, P. Logarithmic Klein-Gordon equation. Acta Phys. Polon. B 2009, 40, 59–66. [Google Scholar]
  24. Ha, T.G.; Park, S.H. Blow-up phenomena for a viscoelastic wave equation with strong damping and logarithmic nonlinearity. Adv. Differ. Equ. 2020, 235, 17. [Google Scholar] [CrossRef]
  25. Han, X.S. Global existence of weak solutions for a logarithmic wave equation arising from Q-ball dynamics. Bull. Korean Math. Soc. 2013, 50, 275–283. [Google Scholar] [CrossRef]
  26. Hiramatsu, T.; Kawasaki, M.; Takahashi, F. Numerical study of Q-ball formation in gravity mediation. J. Cos. Astr. Phys. 2010, 6, 008. [Google Scholar] [CrossRef]
  27. Hu, Q.; Zhang, H.; Liu, G. Global exitence and exponential growth of the solution for the logarithmic Boussinesq-type equation. J. Math. Anal. Appl. 2016, 436, 990–1001. [Google Scholar] [CrossRef]
  28. Lian, W.; Ahmed, M.S.; Xu, R. Global existence and blow up of solution for semilinear hyperbolic equation with logarithmic nonlinearity. Nonlinear Anal. 2019, 184, 239–257. [Google Scholar] [CrossRef]
  29. Yan, J.; Qian, X.; Zhang, H.; Song, S. Two regularized energy-preserving finite difference methods for the logarithmic Klein-Gordon equation. J. Comput. Appl. Math. 2021, 393, 113478. [Google Scholar] [CrossRef]
  30. Zhang, H.; Liu, G.; Hu, Q. Exponential decay of energy for a logarithmic wave equation. J. Partial Differ. Equ. 2015, 28, 269–277. [Google Scholar]
  31. Ait Benhassi, E.M.; Ammari, K.S.; Boulite, S.; Maniar, L. Feedback stabilization of a class of evolution equations with delay. J. Evol. Equ. 2009, 9, 103–121. [Google Scholar] [CrossRef]
  32. Belmiloudi, A. Dynamical behavior of nonlinear impulsive abstract partial differential equations on networks with multiple time-varying delays and mixed boundary conditions involving time-varying delays. J. Dyn. Control Syst. 2015, 21, 95–146. [Google Scholar] [CrossRef]
  33. Belmiloudi, A. Stabilization, Optimal and Robust Control. Theory and Applications in Biological and Physical Sciences; Springer: London, UK, 2008. [Google Scholar]
  34. Belmiloudi, A. Nonlinear optimal control problems of degenerate parabolic equations with logistic time varying delays of convolution type. Nonlinear Anal. 2005, 63, 1126–1152. [Google Scholar] [CrossRef]
  35. Chellaoua, H.; Boukhatem, Y. Blow-up result for an abstract evolution problem with infinite memory and time-varying delay. Appl. Anal. 2022, 101, 4574–4597. [Google Scholar] [CrossRef]
  36. Feng, B.W. Well-posedness and exponential stability for a plate equation with time-varying delay and past history. Z. Angew. Math. Physik. 2017, 68, 1–24. [Google Scholar] [CrossRef]
  37. Kafini, M.; Messaoudi, S. Local existence and blow up of solutions to a logarithmic nonlinear wave equation with delay. Appl. Anal. 2020, 99, 530–547. [Google Scholar] [CrossRef]
  38. Kirane, M.; Said-Houari, B. Existence and asymptotic stability of viscoelastic wave equation with a delay. Z. Angew. Math. Phys. 2011, 62, 1065–1082. [Google Scholar] [CrossRef]
  39. Liu, G.; Diao, L. Energy decay of the solution for a weak viscoelastic equation with a time-varying delay. Acta. Appl. Math. 2018, 155, 9–19. [Google Scholar] [CrossRef]
  40. Liu, W. General decay rate estimate for the energy of a weak viscoealstic equation with internal time-varying delay term. Taiwan. J. Math. 2013, 17, 2101–2115. [Google Scholar] [CrossRef]
  41. Nicaise, S.; Pignotti, C. Stability and instability results of the wave equation with a delay term in the boundary or internal feedbacks. SIAM J. Control Opt. 2006, 45, 1561–1585. [Google Scholar] [CrossRef]
  42. Nicaise, S.; Pignotti, C.; Valein, J. Exponential stability of the wave equation with boundary time-varying delay. Discret. Contin. Dyn. Syst. Ser. S 2011, 4, 693–722. [Google Scholar] [CrossRef]
  43. Nieto-Guadarrama, F.; Villavicencio, J. Zitterbewegung and Klein-tunneling phenomena for transient quantum waves. Phys. Rev. A 2020, 101, 042104. [Google Scholar] [CrossRef]
  44. Park, S.H. Global existence, energy decay and blow-up of solutions for wave equations with time delay and logarithmic source. Adv. Differ. Equ. 2020, 631, 17. [Google Scholar] [CrossRef]
  45. Piskin, E.; Yuksekkaya, H. Local existence and blow up of solutions for a logarithmic nonlinear viscoelastic wave equation with delay. Comput. Methods Differ. Equ. 2021, 9, 623–636. [Google Scholar]
  46. Yang, X.G.; Zhang, J.; Wang, S. Stability and dynamics of a weak viscoelastic system with memory and nonlinear time-varying delay. Discret. Contin. Dyn. Syst. 2020, 40, 1493–1515. [Google Scholar] [CrossRef]
  47. Gross, L. Logarithmic Sobolev inequalities. Am. J. Math. 1975, 97, 1061–1083. [Google Scholar] [CrossRef]
  48. Lions, J.L. Quelques Méthodes de Résolution des Problèmes aux Limites non Linéaires; Dunod; Gauthier-Villars: Paris, France, 1969. [Google Scholar]
  49. Ambrosetti, A.; Rabinowitz, P.H. Dual variational methods in critical point theory and applications. J. Funct. Anal. 1973, 14, 349–381. [Google Scholar] [CrossRef]
  50. Payne, L.E.; Sattinger, D.H. Saddle points and instability of nonlinear hyperbolic equations. Isr. Math. J. 1975, 22, 273–303. [Google Scholar] [CrossRef]
  51. Nakao, M. Decay of solutions for some nonlinear evolution equations. J. Math. Anal. Appl. 1997, 60, 542–549. [Google Scholar] [CrossRef]
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

Belmiloudi, A. On Behavior of Solutions for Nonlinear Klein–Gordon Wave Type Models with a Logarithmic Nonlinearity and Multiple Time-Varying Delays. Axioms 2024, 13, 29. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13010029

AMA Style

Belmiloudi A. On Behavior of Solutions for Nonlinear Klein–Gordon Wave Type Models with a Logarithmic Nonlinearity and Multiple Time-Varying Delays. Axioms. 2024; 13(1):29. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13010029

Chicago/Turabian Style

Belmiloudi, Aziz. 2024. "On Behavior of Solutions for Nonlinear Klein–Gordon Wave Type Models with a Logarithmic Nonlinearity and Multiple Time-Varying Delays" Axioms 13, no. 1: 29. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13010029

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