Next Article in Journal
On Mann Viscosity Subgradient Extragradient Algorithms for Fixed Point Problems of Finitely Many Strict Pseudocontractions and Variational Inequalities
Previous Article in Journal
Approximate Solutions of Time Fractional Diffusion Wave Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On a q—Analog of a Singularly Perturbed Problem of Irregular Type with Two Complex Time Variables

by
Alberto Lastra
1,*,† and
Stéphane Malek
2,†
1
Departamento de Física y Matemáticas, University of Alcalá, Ap. de Correos 20, E-28871 Alcalá de Henares (Madrid), Spain
2
Laboratoire Paul Painlevé, University of Lille 1, 59655 Villeneuve d’Ascq CEDEX, France
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2019, 7(10), 924; https://doi.org/10.3390/math7100924
Submission received: 9 July 2019 / Revised: 27 September 2019 / Accepted: 30 September 2019 / Published: 3 October 2019

Abstract

:
The analytic solutions of a family of singularly perturbed q-difference-differential equations in the complex domain are constructed and studied from an asymptotic point of view with respect to the perturbation parameter. Two types of holomorphic solutions, the so-called inner and outer solutions, are considered. Each of them holds a particular asymptotic relation with the formal ones in terms of asymptotic expansions in the perturbation parameter. The growth rate in the asymptotics leans on the 1 -branch of Lambert W function, which turns out to be crucial.

1. Introduction

The intent of this study is to provide analytic solutions and their parametric asymptotic expansions to a family of singularly perturbed q-difference-differential equations in the complex domain, in which two time variables act, for some fixed q > 1 .
More precisely, we consider equations of the form
Q ( z ) u ( t , z , ϵ ) = P ( t , z , ϵ , t 2 , σ q ; t 1 , σ q ; t 2 ) u ( t , z , ϵ ) + f ( t , z , ϵ ) ,
under null initial data u ( t 1 , 0 , z , ϵ ) u ( 0 , t 2 , z , ϵ ) 0 . Here, Q ( X ) C [ X ] , and P stands for a polynomial with complex coefficients with respect to t : = ( t 1 , t 2 ) , t 2 , a polynomial with rational powers with respect to σ q ; t 1 and σ q ; t 2 , holomorphic with respect to z on a horizontal strip H β : = { z C : | Im ( z ) | < β } for some β > 0 , and holomorphic with respect to the perturbation parameter ϵ on a small disc centered at the origin, say D ( 0 , ϵ 0 ) for some ϵ 0 > 0 . Throughout the present work, σ q ; t stands for the dilation operator on t variable for some fixed q > 1 , i.e.,
σ q ; t ( f ( t ) ) : = f ( q t ) .
We adopt the following notation σ q ; t δ ( f ( t ) ) = f ( q δ t ) , for any δ Q .
The forcing term f ( t , z , ϵ ) is constructed under certain growth conditions and turns out to be a holomorphic function in C 🟉 × C 🟉 × H β × ( D ( 0 , ϵ 0 ) { 0 } ) , for 0 < β < β .
The precise assumptions on the elements involved in the main equation under study are detailed in Section 3.
The problem under study (1) turns out to be a q-analog of the main problem studied in [1],
Q ( z ) u ( t , z , ϵ ) = P ( t 1 k 1 + 1 t 1 , t 2 k 2 + 1 t 2 , z , z , ϵ ) u ( t , z , ϵ ) + f ( t , z , ϵ ) ,
under null initial data u ( t 1 , 0 , z , ϵ ) u ( 0 , t 2 , z , ϵ ) 0 , and where Q ( X ) C [ X ] , and P is a polynomial with respect to its first three variables, with holomorphic coefficients on H β × D ( 0 , ϵ 0 ) , and the forcing term is holomorphic on C 🟉 × C × H β × ( D ( 0 , ϵ 0 ) { 0 } ) . The analytic solutions and their asymptotic expansions of those singularly perturbed partial differential equations are obtained in [1]. More precisely, the so-called inner solutions are holomorphic solutions of (2), holomorphic on domains in time that depend on the perturbation parameter and approach infinity, admits Gevrey asymptotic expansion of certain positive order, with respect to ϵ , whereas the so-called outer solutions are holomorphic solutions of (2), holomorphic on a product of finite sectors with vertex at the origin with respect to the time variables, admit Gevrey asymptotic expansion of a different positive order, with respect to ϵ . In this regard, it is worth mentioning recent works on the action of two complex time variables in different singular perturbed problems on partial differential equations such as [1,2,3].
The previous phenomena of the existence of different asymptotic expansions regarding different domains of the actual solution of the main problem are enhanced in the present work, in the sense that a different nature on the asymptotic expansions is observed: Gevrey and q-Gevrey asymptotic expansions.
In this work, we search for the analytic solutions of the main problem as the inverse Fourier transform and q-Laplace transform of a positive order in the form
u ( t , z , ϵ ) = 1 ( 2 π ) 1 / 2 π q 1 / k 1 L γ ω ( u , m , ϵ ) 1 Θ q 1 / k 1 u ϵ λ 1 t 1 exp ( 1 ϵ λ 2 t 2 k 2 u ) d u u e i z m d m
(see (21) and (23)), for some appropriate γ R and k 1 , k 2 , λ 1 , λ 2 > 0 (see Section 3 for their definition). The function ω ( τ , m , ϵ ) is obtained from a fixed point argument (see Proposition 6), and belongs to q Exp ( k 1 , β , μ , α ) d , a Banach space of holomorphic functions with q-exponential growth and exponential decay with respect to τ and m, respectively (see Definition 4).
The form (3) of the analytic solutions is motivated on the shape of those of the main problem in [1], mixing both time variables in a common Laplace operator. However, the present work is fully self-contained and independent from the constructions performed in [1] and can be read without any knowledge of the results obtained therein.
A first family of analytic solutions of (1) is constructed on domains of the form T 1 × T 2 , ϵ × H β × E h 1 , and a second on domains of the form T 1 × ( T 2 D ( 0 , ρ 2 ) ) × H β × E h 2 0 , where T 1 is a finite sector, T 2 is an unbounded sector and where T 2 , ϵ T 2 is a bounded sector that depends on ϵ E h 1 , and tends to infinity with ϵ approaching the origin. The sets E ̲ 0 = ( E h 2 0 ) 0 h 2 ι 2 1 and E ¯ = ( E h 1 ) 0 h 1 ι 1 1 represent good coverings (see Definition 5).
Different path deformations performed on the analytic solutions give rise to Theorems 2 and 3, where upper bounds on the difference of two consecutive solutions are attained (consecutive solutions in the sense that they are related to consecutive sectors in a good covering). Such bounds are related to null Gevrey and q-Gevrey asymptotic expansions of some positive order. As a matter of fact, the previous differences allow for applying a novel ( ( q , k ) ; s ) -version of the cohomological criteria known as a Ramis–Sibuya theorem. Such result is related to functions admitting q-Gevrey asymptotic expansions of order k and a Gevrey sub-level of order s; see Theorem 4. We also apply a q-analog of Ramis–Sibuya Theorem; see Theorem 5.
The main two results of the present work are Theorems 6 and 7 relating the analytic solutions of (1) to their formal power series expansions obtaining asymptotic results of different nature. Such solutions are known as inner and outer solutions (see Definitions 8 and 9, resp.). Such asymptotic solutions have also been observed in the previous study [1], in the framework of singularly perturbed PDEs. However, the different nature of the asymptotic expansions regarding the outer and inner solutions is a novel phenomenon that has firstly been observed in the present study.
The inner and outer expansions appear in the study of matched asymptotic expansions (see [4,5], among others for the classical theory). In the work [6] by Fruchard and Schäfke, the method of matching is developed, studying the nature of such asymptotic expansions, under Gevrey settings.
We fix a good covering ( E h 1 ) 0 h 1 ι 1 1 , and consider the holomorphic solutions of the problem (1), u h 1 ( t , z , ϵ ) defined on E h 1 w.r.t. ϵ for all 0 h 1 ι 1 1 . In Theorem 6, we prove that, for some μ 2 > λ 2 , θ h 1 and some adequate domain χ 2 T 2 , the function
ϵ u h 1 ( t 1 , x 2 ϵ μ 2 e 1 θ h 1 , z , ϵ )
with values in the Banach space of holomorphic and bounded functions on T 1 × χ 2 × H β , say F 1 , admits a formal power series u ^ ( ϵ ) F 1 [ [ ϵ ] ] as ( ( q , K ) ; S ) -Gevrey asymptotic expansion on E h 1 for some K , S > 0 , for every 0 h 1 ι 1 1 .
The proof of this result leans on the application of accurate estimates related to the 1 -branch of Lambert W function (see Lemma 3).
Concerning the outer solutions of the main problem under study, we consider a good covering ( E h 2 0 ) 0 h 2 ι 2 1 , and the holomorphic solutions of (1), u h 2 ( t , z , ϵ ) defined on E h 2 0 w.r.t. ϵ for all 0 h 2 ι 2 1 . In Theorem 7, we prove that
ϵ u h 2 ( t , z , ϵ ) , 0 h 2 ι 2 1
is an outer solution of (1) with values in the Banach space F 2 of holomorphic and bounded functions on T 1 × ( T 2 D ( 0 , ρ 2 ) ) × H β . Moreover, there exists a formal power series u ^ 0 ( ϵ ) F 2 [ [ ϵ ] ] , which is the common q-Gevrey asymptotic expansion of some positive order of each solution (5) on E h 2 0 , for 0 h 2 ι 2 1 .
In recent years, an increasing interest in the study of the asymptotic behavior of solutions to q-difference-differential equations in the complex domain has been observed. New theories giving rise to q-analogs of the classical theory of Borel–Laplace summability have been discussed and studied, as in the case of the work [7], by Tahara, where the author also provides information about q-analogs of Borel and Laplace transforms and related properties on convolution or Watson-type results. The use of procedures based on the Newton polygon is also exploited in recent studies, such as the work [8] by Tahara and Yamazawa. In addition, it is worth mentioning the study of q-analogs of Briot–Bouquet type partial differential equations by Yamazawa in [9]. Integral transforms involving special functions have also been considered in the study of q-difference-differential equations in [10,11]. Other references in this context by the authors and collaborators are listed in the references.
Different kinds of advanced/delayed partial differential equations are the cornerstone of mathematical models that have been recently studied. Examples of such studies have been applied to tsunamis and rogue waves that can be found in [12]. We also refer to other studies such as [13,14], and the references therein.
The outline of the work is as follows. In Section 2, we recall some known facts about formal q-Borel transform, analytic q-Laplace transform and inverse Fourier transform together with some properties that are applied to transform the main equation under study into auxiliary problems. Afterwards, we provide the definition and related properties of some Banach spaces involved in the construction of the solution. In Section 3, we state the main problem under study (17) and two auxiliary equations. The elements involved in them in addition to the domains of existence and upper bounds of the solutions of such equations are detailed. Section 4 is devoted to the existence and description of the domain of existence for the auxiliary Equation (22) and associated estimates. In the following section, Section 5, we provide analytic solutions of (17) (see Theorem 1) and estimates on the difference of two of them (see Theorems 2 and 3). After a brief summary on q-asymptotic expansions in the first part of Section 6, and the description of Ramis–Sibuya type theorems (see Theorems 4 and 5), we provide formal power series expansions in the perturbation parameter of the analytic solutions and relate them asymptotically in adequate domains. These results are attained in Theorems 6 and 7. The work concludes with a brief section of conclusions and two technical sections, Section 8 and Section 9, left to the end of the work in order to not interfere with our reasonings.

2. Review on Certain Integral Operators: Study of Banach Spaces Involved in the Problem

We briefly describe the foundations of the analytic and formal operators that will allow the transformation of the main problem under study in terms of auxiliary equations. The solutions of such equations belong to certain Banach spaces that are constructed subsequently.

2.1. Review of Some Formal and Analytic Operators

This section is devoted to recalling some of the basic facts on formal and analytic transformations corresponding to the q-analogs of those appearing in the classical Borel–Laplace summability theory. These tools were developed in [15,16].
Through the whole section, q > 1 stands for a real number, and k 1 is a positive integer. E stands for a complex Banach space.
Definition 1.
Given f ^ ( T ) = n 0 a n T n E [ [ T ] ] , the formal q-Borel transform of order k of f ^ ( T ) is defined by
B ^ q ; 1 / k ( f ^ ( T ) ) ( τ ) = n 0 a n τ n ( q 1 / k ) n ( n 1 ) / 2 E [ [ τ ] ] .
The next result, whose proof can be found in Proposition 5 [17], is crucial in order to transform the main equation into an auxiliary one lying in the q-Borel plane.
Proposition 1.
Let m N and j Q . For every f ^ ( T ) E [ [ T ] ] , it holds that
B ^ q ; 1 / k ( T m σ q ; T j f ^ ( T ) ) ( τ ) = τ m ( q 1 / k ) m ( m 1 ) / 2 σ q ; τ j m k B ^ q ; 1 / k ( f ^ ( T ) ) ( τ ) .
The q-analog of Laplace transform used in this work was introduced in [18], and makes use of a kernel given by the Jacobi Theta function of order k, given by
Θ q 1 / k ( x ) = n Z q n ( n 1 ) 2 k x n , x C 🟉 .
We recall that Jacobi Theta function satisfies
Θ q 1 / k q m k x = q m ( m + 1 ) 2 k x m Θ q 1 / k ( x ) , m Z , x C 🟉 .
From Lemma 4.1 [19], given any δ ˜ > 0 , there exists C k , q > 0 (independent of δ ˜ ) such that
Θ q 1 / k ( x ) C q , k δ ˜ exp k 2 log 2 | x | log ( q ) | x | 1 / 2 ,
for all x C 🟉 such that | 1 + x q m k | > δ ˜ , for every m Z . This last property allows for defining a q-analog of Laplace transform with appropriate properties.
Definition 2.
Let ρ > 0 and S d be an unbounded sector with vertex at 0, and bisecting direction d R . Let f : S d D ( 0 , ρ ) E be a holomorphic function, continuous up to the boundary, such that there exist K > 0 and α R with
f ( τ ) E K exp k 2 log 2 | τ | log ( q ) + α log | τ | , τ S d , | τ | ρ ,
f ( τ ) E K , τ D ¯ ( 0 , ρ ) .
We choose an argument γ R within the set of arguments in S d and define the q-Laplace transform of order k of f in direction γ by
L q ; 1 / k γ ( f ( τ ) ) ( T ) = 1 π q 1 / k L γ f ( u ) Θ q 1 / k u T d u u ,
where L γ = { r e γ 1 : r ( 0 , ) } , and π q 1 / k : = log ( q ) k n 0 ( 1 q n + 1 k ) 1 .
The proof of the next results can be found in detail in Lemma 4 and Proposition 6 [17].
Lemma 1.
Let δ ˜ > 0 . In the situation of Definition 2, the integral transform L q ; 1 / k γ ( f ( τ ) ) ( T ) defines a bounded holomorphic function on the domain R γ , δ ˜ D ( 0 , r 1 ) for every 0 < r 1 q ( 1 / 2 α ) / k / 2 , where
R γ , δ ˜ = T C 🟉 : 1 + r e γ 1 T > δ ˜ , f o r a l l r 0 .
The value of L q ; 1 / k γ ( f ( τ ) ) ( T ) does not depend on the choice of γ under e γ 1 S d .
Proposition 2.
Let f be as in Definition 2, and let δ ˜ > 0 . Then, for all σ 0 , one has
T σ σ q j ( L q ; 1 / k γ f ( τ ) ) ( T ) = L q ; 1 / k γ τ σ ( q 1 / k ) σ ( σ 1 ) / 2 σ q j σ k f ( τ ) ( T )
for all T R γ , δ ˜ D ( 0 , r 1 ) , where 0 < r 1 < q ( 1 2 α ) / k / 2 .
We conclude with the definition and properties regarding the inverse Fourier transform.
Definition 3.
Let μ , β R . We write E ( β , μ ) for the vector space of continuous functions h : R C such that
h ( m ) ( β , μ ) = sup m R ( 1 + | m | ) μ exp ( β | m | ) | h ( m ) | < .
The pair ( E ( β , μ ) , · ( β , μ ) ) turns out to be a Banach space.
Proposition 3.
Let μ > 1 , β > 0 and f E ( β , μ ) . The inverse Fourier transform of f is defined by
F 1 ( f ) ( x ) = 1 ( 2 π ) 1 / 2 f ( m ) exp ( i x m ) d m , x R .
This function can be extended to an analytic function on the horizontal strip H β = { z C : | I m ( z ) | < β } . Moreover, it holds that the function ϕ ( m ) = i m f ( m ) belongs to E ( β , μ 1 ) and
z F 1 ( f ) ( z ) = F 1 ( ϕ ) ( z ) , z H β .
In addition to this, it holds that the convolution product of f E ( β , μ ) and g E ( β , μ ) , defined by
ψ ( m ) = 1 ( 2 π ) 1 / 2 f ( m m 1 ) g ( m 1 ) d m 1 ,
is such that ψ E ( β , m ) , and F 1 ( f ) ( z ) F 1 ( g ) ( z ) = F 1 ( ψ ) ( z ) , for every z H β .

2.2. Banach Spaces of Functions of Q-Exponential Growth and Exponential Decay

In this section, we state the definition of the complex Banach space q E x p ( k , β , μ , α ) d . The analytic solution of the main equation under study is built departing from one element in such Banach space via q-Laplace transform and inverse Fourier transform. Similar versions of this Banach space have already appeared in previous works by the authors such as [19,20], and the contribution [17] of the second author.
In the whole section, we fix real numbers β , μ , k > 0 , q , δ > 1 and α . We also set d R and choose an infinite sector S d of bisecting direction d, with vertex at the origin, and the closed disc D ¯ ( 0 , ρ ) , for some ρ > 0 .
Definition 4.
Let q E x p ( k , β , μ , α ) d be the vector space of all continuous complex valued functions ( τ , m ) f ( τ , m ) on ( S d D ¯ ( 0 , ρ ) ) × R which are holomorphic with respect to the first variable on S d D ( 0 , ρ ) and such that
f ( τ , m ) ( k , β , μ , α ) : = sup τ S d D ¯ ( 0 , ρ ) m R ( 1 + | m | ) μ e β | m | exp k log 2 ( | τ | + δ ) 2 log ( q ) α log ( | τ | + δ ) | f ( τ , m ) |
is finite. The pair ( q E x p ( k , β , μ , α ) d , · ( k , β , μ , α ) ) is a Banach space.
In the forthcoming results, we describe some properties on the elements of the previous Banach space when combined with certain operators. The first result is a direct consequence of the previous definition.
Lemma 2.
Let ( τ , m ) a ( τ , m ) be a bounded continuous function defined on ( S d D ¯ ( 0 , ρ ) ) × R , holomorphic with respect to τ on the set S d D ( 0 , ρ ) . For every f ( τ , m ) q E x p ( k , β , μ , α ) d , one has that a ( τ , m ) f ( τ , m ) belongs to q E x p ( k , β , μ , α ) d , and it holds that
a ( τ , m ) f ( τ , m ) ( k , β , μ , α ) C 1 f ( τ , m ) ( k , β , μ , α ) ,
with C 1 : = sup τ S d D ¯ ( 0 , ρ ) m R | a ( τ , m ) | .
Proposition 4.
Let γ 1 , γ 2 , γ 3 0 such that
γ 2 k γ 3 + γ 1 .
Let a γ 1 ( τ ) be a bounded holomorphic function defined on S d D ( 0 , ρ ) with ( 1 + | τ | ) γ 1 | a γ 1 ( τ ) | 1 for every τ S d D ( 0 , ρ ) . Then, for all f q E x p ( k , β , μ , α ) d , it holds that
a γ 1 ( τ ) τ γ 2 σ q ; τ γ 3 f ( τ , m ) ( k , β , μ , α ) C 2 f ( τ , m ) ( k , β , μ , α )
for some C 2 > 0 .
Proof. 
Let f q E x p ( k , β , μ , α ) d . Then, it holds that
a γ 1 ( τ ) τ γ 2 σ q ; τ γ 3 f ( τ , m ) ( k , β , μ , α ) = sup τ S d D ¯ ( 0 , ρ ) m R C 21 ( | τ | ) ( 1 + | m | ) μ e β | m | exp k log 2 ( | τ | q γ 3 + δ ) 2 log ( q ) α log | τ | q γ 3 + δ f τ q γ 3 , m f ( τ , m ) ( k , β , μ , α ) sup τ S d D ¯ ( 0 , ρ ) C 21 ( | τ | ) ,
with
C 21 ( | τ | ) = | τ | γ 2 ( 1 + | τ | ) γ 1 exp k log 2 ( | τ | + δ ) log 2 ( | τ | q γ 3 + δ ) 2 log ( q ) × exp α log ( | τ | + δ ) + α log | τ | q γ 3 + δ .
Observe that σ q ; τ γ 3 ( S d D ¯ ( 0 , ρ ) ) S d D ¯ ( 0 , ρ ) . In order to give upper bounds for C 21 ( | τ | ) , we have
log 2 ( | τ | q γ 3 + δ ) log 2 ( | τ | + δ ) = log 2 ( | τ | + δ q γ 3 ) log 2 ( | τ | + δ ) + ( log ( q ) γ 3 ) 2 2 γ 3 log ( q ) log ( | τ | + δ q γ 3 ) .
The mean value theorem guarantees that
log 2 ( | τ | + δ q γ 3 ) log 2 ( | τ | + δ ) 2 δ ( q γ 3 1 ) max x [ | τ | + δ , | τ | + δ q γ 3 ] log ( x ) x C 22 ,
for some positive constant C 22 that does not depend on τ S d D ¯ ( 0 , ρ ) . This entails that
C 21 ( | τ | ) exp k C 22 2 log ( q ) + ( log ( q ) γ 3 ) 2 | τ | γ 2 ( 1 + | τ | ) γ 1 | τ | / q γ 3 + δ | τ | + δ α 1 ( | τ | + δ q γ 3 ) k γ 3 ,
for every τ S d D ¯ ( 0 , ρ ) . We conclude the result by taking into account (7). □
The next result is stated without proof, heavily rests on Proposition 2 [21], and it is based on bounds stated in Lemma 2.2 [22].
Proposition 5.
Let R 1 , R 2 C [ X ] with d e g ( R 1 ) d e g ( R 2 ) , and such that R 1 ( i m ) 0 for all m R . Let μ > d e g ( R 2 ) + 1 . Given f E ( β , μ ) and g q E x p ( k , β , μ , α ) d , then the function
Φ ( τ , m ) : = 1 R 1 ( i m ) f ( m m 1 ) R 2 ( i m 1 ) g ( τ , m 1 ) d m 1
belongs to q E x p ( k , β , μ , α ) d , and there exists C 2 > 0 such that
Φ ( τ , m ) ( k , β , μ , α ) C 2 f ( m ) ( β , μ ) g ( τ , m ) ( k , β , μ , α ) .

3. Statement of the Main Problem and Auxiliary Equations

Let β > 0 and k 1 , k 1 , k 2 1 with k 1 > k 1 , and D 1 , D 2 2 be integer numbers. We also consider a real number q > 1 , and choose real numbers Δ D 1 D 2 , d D 1 , δ ˜ D 2 0 , and Δ 1 2 , δ 1 , d 1 , δ ˜ 2 0 for every 1 1 D 1 1 and 1 2 D 2 1 . We assume that
δ ˜ 2 k 1 k 1 1 d 1 + d D 1 + δ ˜ D 2 k 1 δ 1 , 1 1 D 1 1 , 1 2 D 2 1 .
In addition to that, there exist natural numbers λ 1 , λ 2 with
Δ D 1 D 2 = λ 1 d D 1 + λ 2 k 2 δ ˜ D 2 ,
k 1 δ 1 < d 1 , λ 1 d 1 + λ 2 k 2 δ ˜ 2 < Δ 1 2 , 1 1 D 1 1 , 1 2 D 2 1 .
Let us also fix polynomials Q , R D 1 D 2 , and R 1 2 for all 1 1 D 1 1 and 1 2 D 2 1 , with complex coefficients such that
Q ( i m ) R D 1 D 2 ( i m ) S Q , R D 1 D 2 , m R ,
where
S Q , R D 1 D 2 : = z C : | z | r Q , R D 1 D 2 , | arg ( z ) d Q , R D 1 D 2 | η Q , R D 1 D 2 ,
for some r Q , R D 1 D 2 , η Q , R D 1 D 2 > 0 and d Q , R D 1 D 2 R . Observe that this condition implies deg ( R D 1 D 2 ) deg ( Q ) . These polynomials are chosen in such a way that R D 1 D 2 ( i m ) 0 for m R , and
deg ( R D 1 D 2 ) deg ( R 1 2 ) , 1 1 D 1 1 , 1 2 D 2 1 .
Let μ > 1 such that
μ > deg ( R 1 2 ) + 1 , 1 1 D 1 1 , 1 2 D 2 1 .
The main problem under study in this work is
Q ( z ) u ( t , z , ϵ ) = ϵ Δ D 1 D 2 ( t 2 k 2 + 1 t 2 ) δ ˜ D 2 t 1 d D 1 σ q ; t 1 d D 1 k 1 R D 1 D 2 ( z ) u ( t , z , ϵ ) + 1 1 D 1 1 1 2 D 2 1 ϵ Δ 1 2 t 1 d 1 σ q ; t 1 δ 1 σ q ; t 2 1 k 2 δ 1 d 1 k 1 ( t 2 k 2 + 1 t 2 ) δ ˜ 2 c 1 2 ( z , ϵ ) R 1 2 ( z ) u ( t , z , ϵ ) + f ( t , z , ϵ ) ,
under initial data u ( t 1 , 0 , z , ϵ ) u ( 0 , t 2 , z , ϵ ) 0 .
The function f is a holomorphic function in C 🟉 × C 🟉 × H β × ( D ( 0 , ϵ 0 ) { 0 } ) for every 0 < β < β , where H β stands for the horizontal strip H β = { z C : | Im ( z ) | < β } . It is constructed as follows. Let ψ ( τ , m , ϵ ) be a continuous function, continuous in C × R × D ( 0 , ϵ 0 ) , for some ϵ 0 > 0 , entire with respect to the first variable, and holomorphic with respect to the third one in D ( 0 , ϵ 0 ) . We moreover assume that there exists C ψ > 0 such that
| ψ ( τ , m , ϵ ) | C ψ ( 1 + | m | ) μ e β | m | exp k 1 2 log ( q ) log 2 ( | τ | + δ ) + α log ( | τ | + δ ) ,
for every ( τ , m , ϵ ) C × R × D ( 0 , ϵ 0 ) , some α R , δ > 0 , and where μ satisfies (16).
In view of the definition of q-Laplace transform and the results described in Section 2.1, one can define
F ( T , m , ϵ ) = 1 π q 1 / k 1 L γ ψ ( u , m , ϵ ) 1 Θ q 1 / k 1 u T 1 e 1 T 2 k 2 u d u u ,
for some fixed γ R . Here, T : = ( T 1 , T 2 ) . The function
f ( t , z , ϵ ) = F 1 m F ( ϵ λ 1 t 1 , ϵ λ 2 t 2 , m , ϵ ) ( z )
turns out to be holomorphic on the set C 🟉 × C 🟉 × H β × ( D ( 0 , ϵ 0 ) { 0 } ) . Observe that, in the previous construction, given T 1 , T 2 C 🟉 , one can choose γ R with cos ( γ k 2 arg ( T 2 ) ) > 0 and | 1 + ( u q m / k ) / T 1 | > 0 for all m Z and with u L γ .
For every 0 1 D 1 1 and 0 2 D 2 1 , the function c 1 2 ( z , ϵ ) is constructed in the following way:
c 1 2 ( z , ϵ ) = 1 ( 2 π ) 1 / 2 C 1 2 ( m , ϵ ) e i z m d m , ( z , ϵ ) H β × D ( 0 , ϵ 0 ) .
c 1 2 ( z , ϵ ) turns out to be a holomorphic function on H β × D ( 0 , ϵ 0 ) whenever m C 1 2 ( m , ϵ ) E ( β , μ ) . In addition to that, we assume that uniform bounds with respect to the perturbation parameter are satisfied, i.e., there exist C 1 2 > 0 with
sup ϵ D ( 0 , ϵ 0 ) C 1 2 ( m , ϵ ) ( β , μ ) C 1 2 .
As mentioned in the title of the work, Equation (17) has the property to be singularly perturbed at ϵ = 0 . Namely, under the above hypotheses, when ϵ equals 0, the nature of the equation changes drastically and is reduced to an ODE with constant coefficients
Q ( z ) u ( t , z , 0 ) = f ( t , z , 0 ) .

3.1. Study of Auxiliary Equations

In this section, we preserve the statements and constructions concerning the main problem under (17), and the geometric and algebraic conditions held on the elements involved in the main equation.
We search for solutions of (17) of the form
u ( t , z , ϵ ) = F 1 ( m U ( ϵ λ 1 t 1 , ϵ λ 2 t 2 , m , ϵ ) ) ( z ) .
Assuming the solution is of the form (21), the expression U ( T 1 , T 2 , m , ϵ ) solves
Q ( i m ) U ( T , m , ϵ ) = T 1 d D 1 σ q ; T 1 d D 1 k 1 ( T 2 k 2 + 1 T 2 ) δ ˜ D 2 R D 1 D 2 ( i m ) U ( T , m , ϵ ) + 1 1 D 1 1 1 2 D 2 1 ϵ Δ 1 2 λ 1 d 1 λ 2 k 2 δ ˜ 2 T 1 d 1 σ q ; T 1 δ 1 σ q ; T 2 1 k 2 ( δ 1 d 1 k 1 ) ( T 2 k 2 + 1 T 2 ) δ ˜ 2 × 1 ( 2 π ) 1 / 2 C 1 2 ( m m 1 , ϵ ) R 1 2 ( i m 1 ) U ( T , m 1 , ϵ ) d m 1 + F ( T , m , ϵ ) .
We reduce the study of solutions of (17) to those of (22), which are linked through (21). In order to solve (22), we adapt a recent approach developed in [1] to a new situation involving both partial differential and q-difference operators. We seek for solutions of (22) of the special form
U γ ( T , m , ϵ ) = 1 π q 1 / k 1 L γ ω ( u , m , ϵ ) 1 Θ q 1 / k 1 u T 1 e 1 T 2 k 2 u d u u ,
for some appropriate function ω ( τ , m , ϵ ) and γ R . We refer to Section 2.1 for the definitions of the elements involved in the previous expression.
Let us consider a second auxiliary equation:
Q ( i m ) ω ( τ , m , ϵ ) = k 2 δ ˜ D 2 ( q 1 / k 1 ) d D 1 ( d D 1 1 ) / 2 τ d D 1 + δ ˜ D 2 R D 1 D 2 ( i m ) ω ( τ , m , ϵ ) + 1 1 D 1 1 1 2 D 2 1 ϵ Δ 1 2 λ 1 d 1 λ 2 k 2 δ ˜ 2 k 2 δ ˜ 2 q ( δ 1 d 1 / k 1 ) δ ˜ 2 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 τ δ ˜ 2 + d 1 × 1 ( 2 π ) 1 / 2 C 1 2 ( m m 1 , ϵ ) R 1 2 ( i m 1 ) ω ( q δ 1 d 1 k 1 τ , m 1 , ϵ ) d m 1 + ψ ( τ , m , ϵ ) .
We define the polynomial P m ( τ ) by
P m ( τ ) = Q ( i m ) k 2 δ ˜ D 2 ( q 1 / k 1 ) d D 1 ( d D 1 1 ) / 2 τ d D 1 + δ ˜ D 2 R D 1 D 2 ( i m ) ,
whose factorization is given by
P m ( τ ) = k 2 δ ˜ D 2 R D 1 D 2 ( i m ) ( q 1 / k 1 ) d D 1 ( d D 1 1 ) / 2 = 0 d D 1 + δ ˜ D 2 1 ( τ q ( m ) ) ,
with
q ( m ) = Δ m exp 1 1 d D 1 + δ ˜ D 2 arg Q ( i m ) ( q 1 / k 1 ) d D 1 ( d D 1 1 ) / 2 k 2 δ ˜ D 2 R D 1 D 2 ( i m ) + 2 π d D 1 + δ ˜ D 2 ,
for all 0 d D 1 + δ ˜ D 2 1 , where Δ m = | Q ( i m ) | k 2 δ ˜ D 2 | R D 1 D 2 ( i m ) | ( q 1 / k 1 ) d D 1 ( d D 1 1 ) 2 1 / ( d D 1 + δ ˜ D 2 ) .
Let d R be such that the infinite sector S d of bisecting direction d satisfies the following geometric construction: there exists M 1 > 0 such that | τ q ( m ) | M 1 ( 1 + | τ | ) , for all 0 d D 1 + δ ˜ D 2 1 , m R and τ S d D ¯ ( 0 , ρ ) . The previous is a feasible condition for an appropriate choice of small enough η Q , R D 1 D 2 > 0 and large enough r Q , R D 1 D 2 > 0 , in view of (14) and the definition of q ( m ) . More precisely, one chooses S d such that q ( m ) / τ has positive distance to 1, for all 0 d D 1 + δ ˜ D 2 1 , m R , and τ S d D ¯ ( 0 , ρ ) . The previous choice of d yields
| P m ( τ ) | C P | R D 1 D 2 ( i m ) | ( 1 + | τ | ) d D 1 + δ ˜ D 2 ,
for some C P > 0 , valid for all τ S d D ¯ ( 0 , ρ ) and all m R .
Proposition 6.
Let S d be an infinite sector with vertex at the origin satisfying the previous geometric conditions. Then, if C ψ > 0 (see (18) for its definition) and ϵ 0 > 0 are small enough, there exists ϖ > 0 such that Equation (24) admits a unique solution ω d ( τ , m , ϵ ) which belongs to q E x p ( k 1 , β , μ , α ) d with ω d ( τ , m , ϵ ) ( k 1 , β , μ , α ) ϖ , for all ϵ D ( 0 , ϵ 0 ) .
Proof. 
Let ϵ D ( 0 , ϵ 0 ) . We consider the map H ϵ defined by
H ϵ ( ω ( τ , m ) ) : = 1 P m ( τ ) 1 1 D 1 1 1 2 D 2 1 ϵ Δ 1 2 λ 1 d 1 λ 2 k 2 δ ˜ 2 k 2 δ ˜ 2 q ( δ 1 d 1 / k 1 ) δ ˜ 2 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 τ δ ˜ 2 + d 1 × 1 ( 2 π ) 1 / 2 C 1 2 ( m m 1 , ϵ ) R 1 2 ( i m 1 ) ω ( q δ 1 d 1 k 1 τ , m 1 ) d m 1 + 1 P m ( τ ) ψ ( τ , m , ϵ ) .
Given ϖ > 0 , let ω ( τ , m ) q E x p ( k 1 , β , μ , α ) d with ω d ( τ , m , ϵ ) ( k 1 , β , μ , α ) ϖ . In view of (11)–(13), (26) and from Lemma 2, Propositions 4 and 5, one has
ϵ Δ 1 2 λ 1 d 1 λ 2 k 2 δ ˜ 2 k 2 δ ˜ 2 q ( δ 1 d 1 / k 1 ) δ ˜ 2 ( 2 π ) 1 / 2 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 τ δ ˜ 2 + d 1 P m ( τ ) × C 1 2 ( m m 1 , ϵ ) R 1 2 ( i m 1 ) ω ( q δ 1 d 1 k 1 τ , m 1 ) d m 1 ( k 1 , β , μ , α ) ϵ 0 Δ 1 2 λ 1 d 1 λ 2 k 2 δ ˜ 2 k 2 δ ˜ 2 q ( δ 1 d 1 / k 1 ) δ ˜ 2 ( 2 π ) 1 / 2 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 C P C 2 C 2 C 1 2 ω ( τ , m , ϵ ) ( k 1 , β , μ , α )
for every 0 1 D 1 1 and 0 2 D 2 1 . In addition to this, Lemma 2 and (26) yield
1 P m ( τ ) ψ ( τ , m , ϵ ) ( k 1 , β , μ , α ) 1 C P sup m R 1 | R D 1 D 2 ( i m ) | ψ ( τ , m , ϵ ) ( k 1 , β , μ , α ) 1 C P sup m R 1 | R D 1 D 2 ( i m ) | C ψ .
Let ϵ 0 , C ψ > 0 be small enough satisfying
1 1 D 1 1 1 2 D 2 1 ϵ 0 Δ 1 2 λ 1 d 1 λ 2 k 2 δ ˜ 2 k 2 δ ˜ 2 q ( δ 1 d 1 / k 1 ) δ ˜ 2 ( 2 π ) 1 / 2 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 C P C 2 C 2 C 1 2 ϖ + 1 C P sup m R 1 | R D 1 D 2 ( i m ) | C ψ ϖ .
Then, the estimates (27) and (28) yield to H ( ω ( τ , m ) ( k 1 , β , μ , α ) ϖ . In other words, the operator H ϵ restricted to D ¯ ( 0 , ϖ ) q E x p ( k 1 , β , μ , α ) d is such that H ϵ ( D ¯ ( 0 , ϖ ) ) D ¯ ( 0 , ϖ ) .
On the other hand, let ω 1 ( τ , m ) , ω 2 ( τ , m ) D ¯ ( 0 , ρ ) q E x p ( k 1 , β , μ , α ) d . Analogously to (27), one arrives at
H ϵ ( ω 1 ( τ , m ) ) H ϵ ( ω 2 ( τ , m ) ) ( k 1 , β , μ , α ) 1 2 ω 1 ( τ , m ) ω 2 ( τ , m ) ( k 1 , β , μ , α )
by choosing ϵ 0 > 0 such that
1 1 D 1 1 1 2 D 2 1 ϵ 0 Δ 1 2 λ 1 d 1 λ 2 k 2 δ ˜ 2 k 2 δ ˜ 2 q ( δ 1 d 1 / k 1 ) δ ˜ 2 ( 2 π ) 1 / 2 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 C P C 2 C 2 C 1 2 1 2 .
We conclude that the map H ϵ : D ¯ ( 0 , ϖ ) q E x p ( k 1 , β , μ , α ) d q E x p ( k 1 , β , μ , α ) d is contractive. The classical fixed point theory in complete metric spaces states the existence of a unique fixed point for H ϵ , say ω d ( τ , m , ϵ ) , in q E x p ( k 1 , β , μ , α ) d , with ω d ( τ , m , ϵ ) ( k 1 , β , μ , α ) ϖ . For every ϵ D ( 0 , ϵ 0 ) , the function ω d ( τ , m , ϵ ) is a solution of (24) in view of the definition of the operator H ϵ . Holomorphy of the map D ( 0 , ϵ 0 ) ϵ ω d ( τ , m , ϵ ) is derived from the construction of the fixed point. □
In order to prove that the solutions of (22) and (24) are related via (23), we need to clarify how operators involved in (22) are transformed into the corresponding ones in (24). This is left until the end of the work in Section 9 not to interfere with our line of reasoning.

4. Domains of Existence for the Solutions of (22) and Associated Estimates

In this section, we describe appropriate domains on the time variables in which the solution of the main problem under consideration is well defined, within appropriate geometric conditions. Let T ˜ 1 be a bounded sector with a vertex at the origin such that there exists δ 1 > 0 with
1 + r e γ 1 T 1 δ 1 ,
for all r 0 and all T 1 T ˜ 1 , and γ being an argument in S d . We also fix an unbounded sector T ˜ 2 , with vertex at the origin, such that
γ k 2 arg ( T 2 ) π 2 + δ 2 , π 2 δ 2 ,
for some δ 2 > 0 and well chosen γ arg ( S d ) . Observe that, in particular, there exists δ 3 > 0 with cos ( γ k 2 arg ( T 2 ) ) > δ 3 .
The next technical result describes accurate bounds for the solutions of (22) in different domains. The proof is left to Section 8.
Proposition 7.
Let U γ ( T , m , ϵ ) be defined in (23), with ω ( τ , m , ϵ ) = ω d ( τ , m , ϵ ) being the function obtained in Proposition 6. The following statements hold:
1
There exist small enough ρ 1 > 0 and large enough ρ 2 > 0 such that
| U γ ( T , m , ϵ ) | C ˜ 1 ( 1 + | m | ) μ e β | m | 1 + | T 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ | T 1 | | T 2 | k 2 δ 3 e k 1 2 log ( q ) log 2 ( | T 2 | k 2 δ 3 ) + α log ( | T 2 | k 2 δ 3 ) 1 + log ( | T 2 | k 2 δ 3 ) 1 / 2 e k 1 log ( q ) log 2 | T 2 | k 2 δ 3 log k 1 log ( q ) log | T 2 | k 2 δ 3 ,
for every T 1 T ˜ 1 with | T 1 | < ρ 1 and T 2 T ˜ 2 with | T 2 | > ρ 2 , and some C ˜ 1 > 0 .
2
There exist small enough ρ 1 , ρ 2 > 0 such that
| U γ ( T , m , ϵ ) | C ˜ 2 ( 1 + | m | ) μ e β | m | 1 + | T 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ | T 1 | × exp ρ δ 3 2 | T 2 | k 2 ,
for every T 1 T ˜ 1 with | T 1 | < ρ 1 and T 2 T ˜ 2 with | T 2 | < ρ 2 , and some C ˜ 2 > 0 .

5. Analytic Solutions of the Main Problem: Inner and Outer Solutions

In this section, we preserve the values of the elements involved in the main problem (17) stated in Section 3. More precisely, we assume (11)–(16), and also the hypotheses on the forcing term (19) in (18) and the coefficients in (20). Let d R and S d be an infinite sector with vertex at 0 C under the geometric condition imposed in Proposition 6. Our main aim is to construct analytic solutions of (17) and their asymptotic behavior in different domains. For this purpose, we consider the analytic solutions as stated in Section 3.1.
Such solutions are defined in families of sectors with respect to the perturbation parameter, conforming good coverings of C 🟉 (see Definition 5). We also provide information about the difference of two solutions in consecutive sectors of the good covering, which will be crucial to determine the asymptotic behavior of the analytic solutions. We refer to consecutive solutions to solutions that are associated with consecutive elements in a fixed good covering of C 🟉 . Let us first recall the notion of good covering in C 🟉 .
Definition 5.
Let ι 2 be an integer. For every 0 h ι 1 , we choose a finite sector with vertex at the origin E h such that:
  • E h D ( 0 , ϵ 0 ) , and E j 1 E j 2 = if and only if 0 j 1 , j 2 ι 1 with | j 1 j 2 | 2 (under the convention that E ι : = E 0 ).
  • h = 0 ι 1 E h = U { 0 } , for some neighborhood of the origin U .
    A family of sectors ( E h ) 0 h ι 1 under these assumptions is known as a good covering in C 🟉 .
Let T ˜ 1 and T ˜ 2 be sectors following the construction in Section 4.
Definition 6.
Let ( E h ) 0 h ι 1 be a good covering in C 🟉 . We also fix a bounded sector, T 1 , and an unbounded sector T 2 , both with vertex at the origin. For all 0 h ι 1 , let S d h be an infinite sector of bisecting direction d h R . We say that the set { T 1 , T 2 , ( E h ) 0 h ι 1 , ( S d h ) 0 h ι 1 } is admissible if the following conditions hold:
  • For every 0 h ι 1 , ϵ E h and t 1 T 1 we have ϵ λ 1 t 1 T ˜ 1 .
  • For every 0 h ι 1 , ϵ E h and t 2 T 2 we have ϵ λ 2 t 2 T ˜ 2 .
Observe that given an admissible set { T 1 , T 2 , ( E h ) 0 h ι 1 , ( S d h ) 0 h ι 1 } , the sectors T ˜ 1 , the choice of the sectors T ˜ 1 , T ˜ 2 and γ : = γ d S d h , fixed in Section 2.1, entail the existence of δ 1 > 0 with
1 + r e γ d h 1 ϵ λ 1 t 1 δ 1 , r 0 , t 1 T 1 , ϵ E h , 0 h ι 1 .
In addition to this, it holds that
γ d h k 2 arg ( ϵ λ 2 t 2 ) π 2 + δ 2 , π 2 δ 2 , ϵ E h , 0 h ι 1 ,
for some δ 2 > 0 and cos ( γ d h k 2 arg ( ϵ λ 2 t 2 ) ) > δ 3 , for some δ 3 > 0 , ϵ E h , 0 h ι 1 .
Definition 7.
Let { T 1 , T 2 , ( E h ) 0 h ι 1 , ( S d h ) 0 h ι 1 } be an admissible set.
Then, { ( S d h ) 0 h ι 1 , T 1 × T 2 } is known as a family of sectors associated with the good covering ( E h ) 0 h ι 1 .
Theorem 1.
Let ( E h ) 0 h ι 1 be a good covering in C 🟉 . For every 0 h ι 1 , we choose d h R such that S d h satisfies the geometric conditions of Proposition 6. Let { ( S d h ) 0 h ι 1 , T 1 × T 2 } be a family of sectors associated with the good covering ( E h ) 0 h ι 1 . If there exist small enough ϵ 0 , C ψ > 0 such that, if
sup ϵ D ( 0 , ϵ 0 ) ψ ( τ , m , ϵ ) ( k 1 , β , μ , α ) C ψ ,
then, for every 0 h ι 1 , the problem (17) admits a solution u h ( t , z , ϵ ) , which defines a bounded and holomorphic function in T 1 × T 2 × H β × E h , for any fixed 0 < β < β .
Proof. 
Let 0 h ι 1 . Proposition 6 guarantees the existence of ϵ 0 , C ψ > 0 such that the Equation (24) admits a unique solution ω d h ( τ , m , ϵ ) which belongs to q E x p ( k 1 , β , μ , α ) d h for every ϵ D ( 0 , ϵ 0 ) , and the map ϵ ω d h ( τ , m , ϵ ) is holomorphic in D ( 0 , ϵ 0 ) . Taking into account that { ( S d h ) 0 h ι 1 , T 1 × T 2 } is associated with the good covering ( E h ) 0 h ι 1 and the properties of Laplace transform stated in Section 2.1, one can construct U γ h ( T , m , ϵ ) in the form (23), which is well defined, bounded and continuous function on T ˜ 1 × T ˜ 2 × R × D ( 0 , ϵ 0 ) , where T ˜ 1 R γ h , δ ˜ D ( 0 , r 1 ) for 0 < r 1 q ( 1 / 2 α ) / k 1 / 2 , and T ˜ 2 is an infinite sector. The function U γ h ( T , m , ϵ ) is holomorphic w.r.t. ( T , ϵ ) on T ˜ 1 × T ˜ 2 × D ( 0 , ϵ 0 ) . As a matter of fact, U γ h is a solution of (22) in view of the properties relating both equalities in Section 9. We finally define u h ( t , z , ϵ ) following (21):
u h ( t , z , ϵ ) = F 1 ( m U γ h ( ϵ λ 1 t 1 , ϵ λ 2 t 2 , m , ϵ ) ) ( z ) ,
which turns out to be a holomorphic solution of (17), defined on T 1 × T 2 × H β × E h , for any fixed 0 < β < β . □
In order to provide the asymptotic behavior of the analytic solutions of (17) in different domains, with respect to the perturbation parameter ϵ , we state the definition of inner and outer solutions of the problem (17).

5.1. Inner Solutions of the Main Problem

Definition 8.
Let ι 1 2 . Let ( E h 1 ) 0 h 1 ι 1 1 be a good covering of C 🟉 . We also consider the admissible set { T 1 , T 2 , ( E h 1 ) 0 h 1 ι 1 1 , ( S d h 1 ) 0 h 1 ι 1 1 } . Let μ 2 > 0 be a natural number satisfying
μ 2 > λ 2 , k 1 λ 1 2 > k 1 ( ( μ 2 λ 2 ) k 2 ) 2 .
Let χ 2 be a bounded domain, such that the good covering ( E h 1 ) 0 h 1 ι 1 1 satisfies the following condition: for all 0 h 1 ι 1 1 we can select θ h 1 R (which depends on E h 1 ) such that, for every x 2 χ 2 and ϵ E h 1 , the complex number t 2 = x 2 ϵ μ 2 e θ h 1 1 belongs to T 2 . We define the set T 2 , ϵ , μ 2 : = { x 2 ϵ μ 2 e θ h 1 1 : x 2 χ 2 } .
In case ϵ E h 1 , t 2 T 2 , ϵ , μ 2 , t 1 T 1 , z H β for 0 < β < β , then we say that u h 1 ( t , z , ϵ ) represents an inner solution of (17).
Theorem 2.
Under the assumptions of Theorem 1 and the constraints on the inner solutions of the main problem of Definition 8, let ( E h 1 ) 0 h 1 ι 1 be a good covering of C 🟉 . Then, there exists C i n n > 0 such that, for every 0 h 1 ι 1 1 , ϵ E h 1 E h 1 + 1 , t 2 T 2 , ϵ , μ 2 , z H β for any fixed 0 < β < β and t 1 T 1 , one has
| u h 1 + 1 ( t , z , ϵ ) u h 1 ( t , z , ϵ ) | C i n n | ϵ | Δ 1 exp k 1 2 log ( q ) log 2 C ^ 1 | ϵ | λ 1 + k 1 log ( q ) log C ^ 2 | ϵ | ( μ 2 λ 2 ) k 2 log log ( C ^ 3 | ϵ | ( μ 2 λ 2 ) k 2 ) ,
for some C ^ 1 , C ^ 2 , C ^ 3 > 0 , Δ 1 R and a real number k 1 ( 0 , k 1 ) chosen to be small enough.
Proof. 
Let 0 h 1 ι 1 1 and consider consecutive solutions u h 1 , u h 1 + 1 of (17), constructed in Theorem 1. We recall that the function ω d h 1 ( τ , m , ϵ ) stands for the analytic continuation of a function ω ( τ , m , ϵ ) , holomorphic w.r.t. τ in a neighborhood of the origin D ( 0 , ρ ) , to the infinite sector S d h 1 . This entails that the difference u h 1 + 1 ( t , z , ϵ ) u h 1 ( t , z , ϵ ) can be written in the following form, after an appropriate path deformation in τ , avoiding the roots of P m ( τ ) (see Figure 1).
We write u h 1 + 1 ( t , z , ϵ ) u h 1 ( t , z , ϵ ) = J 1 J 2 + J 3 , where
J 1 = 1 ( 2 π ) 1 / 2 1 π q 1 / k 1 L γ h 1 + 1 , ρ / 2 ω d h 1 + 1 ( u , m , ϵ ) 1 Θ q 1 / k 1 u ϵ λ 1 t 1 e 1 ϵ λ 2 t 2 u d u u e i z m d m , J 2 = 1 ( 2 π ) 1 / 2 1 π q 1 / k 1 L γ h 1 , ρ / 2 ω d h 1 ( u , m , ϵ ) 1 Θ q 1 / k 1 u ϵ λ 1 t 1 e 1 ϵ λ 2 t 2 u d u u e i z m d m , J 3 = 1 ( 2 π ) 1 / 2 1 π q 1 / k 1 C h 1 , h 1 + 1 , ρ / 2 ω ( u , m , ϵ ) 1 Θ q 1 / k 1 u ϵ λ 1 t 1 e 1 ϵ λ 2 t 2 u d u u e i z m d m ,
where L γ j , ρ / 2 = [ ρ / 2 e γ j 1 , ) for j = h 1 , h 1 + 1 , and C h 1 , h 1 + 1 , ρ / 2 stands for the arc of circle from ρ / 2 e γ h 1 1 to ρ / 2 e γ h 1 + 1 1 .
Assume that ϵ E h 1 E h 1 + 1 and t 2 T 2 , ϵ , μ 2 , z H β for some 0 < β < β and t 1 T 1 . Owing to the estimates leading to (29), displayed in (49), we derive that dist ( χ 2 ϵ λ 2 μ 2 , 0 ) ρ 2 for some large enough ρ 2 . Moreover, if ρ 0 > 0 is small enough, then there exists C ˜ 8 > 0 such that
| J 1 | C ˜ 8 | ϵ λ 1 t 1 | 1 / 2 e k 1 2 log ( q ) log 2 ρ / 2 | ϵ λ 1 t 1 | | ϵ λ 2 μ 2 x 2 | k 2 δ 3 × e k 1 2 log ( q ) log 2 ( | ϵ λ 2 μ 2 x 2 | k 2 δ 3 ) + α log ( | ϵ λ 2 μ 2 x 2 | k 2 δ 3 ) 1 + log ( | ϵ λ 2 μ 2 x 2 | k 2 δ 3 ) 1 / 2 × e k 1 log ( q ) log | ϵ λ 2 μ 2 x 2 | k 2 δ 3 log k 1 log ( q ) log | ϵ λ 2 μ 2 x 2 | k 2 δ 3 ( 1 + | m | ) μ e | m | ( β β ) d | m |
for every t 1 T 1 , z H β , ϵ E h 1 E h 1 + 1 , and x 2 χ 2 . An analogous upper bound is attained for | J 2 | . Concerning | J 3 | , one can apply the bounds stated in Proposition 6, and (6) to get the existence of C ˜ 9 > 0 such that
| J 3 | C ˜ 9 2 ρ 3 / 2 | ϵ λ 1 t 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ / 2 | ϵ λ 1 t 1 | exp ρ / 2 | ϵ λ 2 μ 2 x 2 | k 2 δ 3 C ˜ 9 2 ρ 3 / 2 | ϵ λ 1 t 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ / 2 | ϵ λ 1 t 1 | .
The result follows taking into account (34) and (35), under the condition (31). □

5.2. Outer Solutions of the Main Problem

Definition 9.
Let ι 2 2 and ( E h 2 0 ) 0 h 2 ι 2 be a good covering of C 🟉 , and consider an admissible set { T 1 , T 2 , ( E h 2 0 ) 0 h 2 ι 2 1 , ( S d h 2 0 ) 0 h 2 ι 2 1 } . Assume that t 2 T 2 is such that | t 2 | < ρ 2 for some fixed ρ 2 > 0 , which is independent of ϵ, then we say that u h 2 ( t , z , ϵ ) represents an outer solution of (17).
Theorem 3.
Under the assumptions of Theorem 1 and the constraints on the outer solutions of the main problem of Definition 9, let ( E h 2 0 ) 0 h 2 ι 2 be a good covering of C 🟉 . Then, there exists C o u t > 0 such that, for every 0 h 2 ι 2 1 , ϵ E h 2 0 E h 2 + 1 0 , t 2 T 2 with | t 2 | < ρ 2 , z H β for any fixed 0 < β < β and t 1 T 1 , one has
| u h 2 + 1 ( t , z , ϵ ) u h 2 ( t , z , ϵ ) | C o u t | ϵ | Δ 2 exp k 1 2 log ( q ) log 2 C ^ 4 | ϵ | λ 1 C ^ 5 | ϵ | λ 2 k 2 ,
for some C ^ 4 , C ^ 5 > 0 and Δ 2 R .
Proof. 
Let 0 h 2 ι 2 1 and consider consecutive solutions u h 2 , u h 2 + 1 of (17), constructed in Theorem 1. We proceed to write the difference of two consecutive solutions in the form
u h 2 + 1 ( t , z , ϵ ) u h 2 ( t , z , ϵ ) = J 1 J 2 + J 3 ,
for all 0 h 2 ι 2 1 , with h 1 substituted by h 2 in the expressions of J 1 , J 2 , J 3 in the proof of Theorem 2. Let ϵ E h 2 0 E h 2 + 1 0 and t 2 T 2 , with | t 2 | < ρ 2 , z H β for some 0 < β < β and t 1 T 1 . Owing to analogous bounds as those leading to (50) and (51), we arrive at
| J 1 | C ˜ 10 | ϵ λ 1 t 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ | ϵ λ 1 t 1 | exp ρ / 2 2 | ϵ λ 2 t 2 | k 2 δ 3 ,
for some C ˜ 10 > 0 . Similar estimates hold for | J 2 | . On the other hand, direct computations yield
| J 3 | C ˜ 11 exp k 1 2 log ( q ) log 2 ρ / 2 | ϵ λ 1 t 1 | | ϵ λ 1 t 1 | ρ / 2 1 / 2 exp ρ / 2 | ϵ λ 2 t 2 | k 2 δ 1 2 ρ ,
for some C ˜ 11 > 0 . This concludes the proof. □

6. Asymptotic Expansions of Mixed Order

This section is divided in two parts. The first part recalls some facts about q-asymptotic expansions and also describes q-asymptotic expansions which show a sub-Gevrey growth in their estimates, and related results.
In the second part of this section, we provide the existence of a formal solution of the main problem under study, written as a formal power series in the perturbation parameter, and explain the asymptotic relationship between this and the analytic solutions.

6.1. Review on q-Asymptotic Expansions

In the whole subsection, ( F , · ) stands for a complex Banach space.
We first recall the notion of q-Gevrey asymptotic expansions, which can be also found in [20] in more detail.
Definition 10.
Let V be a bounded open sector with vertex at 0 in C , q R with q > 1 . We also fix a positive integer k. We say that a holomorphic function f : V F admits the formal power series f ^ ( ϵ ) = n 0 f n ϵ n F [ [ ϵ ] ] as its q-Gevrey asymptotic expansion of order 1 / k if for every open subsector U of V, i.e., U ¯ { 0 } V , there exist A , C > 0 such that
f ( ϵ ) n = 0 N f n ϵ n F C A N + 1 q N ( N + 1 ) 2 k | ϵ | N + 1 ,
for every ϵ U and N 0 .
Such pure q-asymptotic expansions have been recently studied when dealing with the asymptotic behavior of the solutions of q-difference-differential equations in the complex domain [17,19,20,23].
We recall the definition of asymptotic expansion of mixed order as introduced in the work [24].
Definition 11.
Let V be a bounded open sector of the complex plane with vertex at the origin. Let q > 1 and s 0 be real numbers, and let k be a positive integer. The function f : V F is said to admit the formal series F ^ ( ϵ ) = n 0 f n ϵ n F [ [ ϵ ] ] as its ( ( q , k ) ; s ) Gevrey asymptotic expansion if for every open subsector U of V (i.e., ( U ¯ { 0 } ) V ) there exist A , C > 0 such that
f ( ϵ ) n = 0 N f n ϵ n F C A N N ! s q N ( N + 1 ) 2 k | ϵ | N + 1 ,
for all ϵ U and N 0 .
Observe that the set of functions admitting q-Gevrey asymptotic expansions of order 1 / k coincide with the functions admitting ( ( q , k ) ; 0 ) Gevrey asymptotic expansion.
We now proceed to state a ( ( q , k ) ; s ) -version of the Ramis–Sibuya Theorem. The classical statement of this result involves Gevrey asymptotic expansions, and guarantees s-summability of some power series. This cohomological criterion can be found in Proposition 2 [25] and Lemma XI-2-6 [26].
This version has already been stated in the work [24] where a complete proof can be found therein.
Theorem 4
( ( ( q , k ) ; s ) –Ramis–Sibuya Theorem). Let ( E h ) 0 h ι 1 be a good covering in C 🟉 . For every 0 h ι 1 , we consider a holomorphic function G h : E h F and define Δ h : = G h + 1 G h , which turns out to be a holomorphic function in Z h : = E h E h + 1 (here, E ι and G ι stand for E 0 and G 0 , respectively). Assume that the following statements hold:
  • G h is a bounded function in a vicinity of 0, for every 0 h ι 1 .
  • The function Δ h ( ϵ ) admits null ( ( q , k ) ; s ) Gevrey asymptotic expansion in Z h for every 0 h ι 1 , i.e., there exist C 1 , C 2 > 0 such that
    Δ h ( ϵ ) F C 1 C 2 N N ! s q N ( N + 1 ) 2 k | ϵ | N + 1 ,
    for every ϵ Z h and all N 0 .
Then, there exists G ^ F [ [ ϵ ] ] , which is the common ( ( q , k ) ; s ) -Gevrey asymptotic expansion of G h ( ϵ ) on E h , for all 0 h ι 1 .
In [17], a q-Gevrey version of Ramis–Sibuya is obtained. That version is related to q-Gevrey asymptotic expansions of some positive order k, which coincide with ( q , k ) ; 0 ) -Gevrey asymptotic expansions in our framework. The Ramis–Sibuya theorem in that framework reads as follows:
Theorem 5 (q-RS)).
Let ( F , · F ) be a Banach space and ( E h ) 0 h ι 1 be a good covering in C 🟉 . For every 0 h ι 1 , let G h ( ϵ ) be a holomorphic function from E h into F and let the cocycle Δ h ( ϵ ) = G h + 1 ( ϵ ) G h ( ϵ ) be a holomorphic function from Z h = E h + 1 E h into F (with the convention that E ι = E 0 and G ι = G 0 ). We make further assumptions:
1
The functions G h ( ϵ ) are bounded as ϵ tends to 0 on E h , for all 0 h ι 1 .
2
The function Δ h ( ϵ ) is q-exponentially flat of order k on Z h for all 0 h ι 1 , meaning that there exist two constants C h 1 R and C h 2 > 0 with
Δ h ( ϵ ) F C h 2 | ϵ | C h 1 exp k 2 log ( q ) log 2 | ϵ | ,
for all ϵ Z h , and all 0 h ι 1 .
Then, there exists a formal power series G ^ ( ϵ ) F [ [ ϵ ] ] which is the common q-Gevrey asymptotic expansion of order 1 / k of the function G h ( ϵ ) on E h , for all 0 h ι 1 .

6.2. Asymptotic Expansions for the Analytic Solutions of the Main Problem

In this section, we preserve all the assumptions made on the elements involved in the main problem (17), detailed in Section 3. Moreover, we depart from the geometric construction of the elements used to construct the analytic solutions of (17), described in Section 4 and Section 5, and under Assumption (31).
The next result shows that the difference of two consecutive solutions of the main problem allow the application of the ( ( q , k ) ; s ) -version of Ramis–Sibuya Theorem, obtained in Theorem 4, in adequate domains.
Lemma 3.
Let k 1 , k 1 , k 2 , λ 1 , λ 2 , μ 2 be positive constants with μ 2 > λ 2 . For every n N , we consider the function
Ψ ( x ) = x n exp k 1 λ 1 2 2 log ( q ) log 2 ( x ) + k 1 ( μ 2 λ 2 ) k 2 log ( q ) log ( x ) log ( log ( x ) ) .
Then, it holds that
Ψ ( x ) C ^ A ^ n n ! S q n ( n + 1 ) K , x > 0 ,
for n n 0 , where n 0 1 is a large enough integer depending on k 1 , k 1 , k 2 , λ 1 , λ 2 , μ 2 with
S = k 1 ( μ 2 λ 2 ) k 2 k 1 λ 1 2 , K = 9 k 1 λ 1 2 7 .
Proof. 
We make the change of variable τ = log ( x ) and consider the auxiliary function Ψ ¯ ( τ ) : = Ψ ( e τ ) . We search for the maximum of Ψ ¯ for τ R . It holds that
Ψ ¯ ( τ ) = Ψ ( e τ ) n k 1 λ 1 2 log ( q ) τ + k 1 ( μ 2 λ 2 ) k 2 log ( q ) ( log ( τ ) + 1 ) .
It is known that the solution of the equation log ( t ) + log ( b ) t = log ( a ) , for some fixed a , b > 0 , is given by t = W ( a log ( b ) ) / log ( b ) , where W is the Lambert W function. The maximum of Ψ ¯ ( τ ) is then attained at
τ = τ n = k 1 ( μ 2 λ 2 ) k 2 k 1 λ 1 2 W 1 exp n log ( q ) + k 1 ( μ 2 λ 2 ) k 2 k 1 ( μ 2 λ 2 ) k 2 k 1 λ 1 2 k 1 ( μ 2 λ 2 ) k 2 ,
where W 1 stands for the 1 -branch of Lambert W function. Let
A : = k 1 ( μ 2 λ 2 ) k 2 , C : = k 1 λ 1 2 k 1 ( μ 2 λ 2 ) k 2 .
Then, τ n can be written in the form τ n = 1 C W 1 ( exp ( n log ( q ) A 1 ) C ) . Regarding Theorem 1 [27], we have
1 2 u u < W 1 ( e u 1 ) < 1 2 u 2 3 u , u > 0 ,
which can be applied to estimate τ n :
1 + 2 u 0 + 2 / 3 u 0 C < 1 C W 1 ( e u 0 1 ) < 1 + 2 u 0 + u 0 C , u 0 = n log ( q ) A log ( C ) ,
leading to
1 C 1 + 2 ( n log ( q ) A log ( C ) ) + 2 3 ( n log ( q ) A log ( C ) ) < τ n τ n < 1 C 1 + 2 n log ( q ) A log ( C ) + n log ( q ) A log ( C ) .
We conclude that Ψ ¯ ( τ ) Ψ ¯ ( τ n ) . Taking into account (41), we arrive at
Ψ ¯ ( τ n ) = exp n τ n k 1 λ 1 2 2 log ( q ) τ n 2 + k 1 ( μ 2 λ 2 ) k 2 log ( q ) τ n log ( τ n ) < C ^ 1 A ^ 1 n exp n 2 log ( q ) 1 A C 2 k 1 λ 1 2 9 C 2 A 2 ξ + n log ( n log ( q ) C A ) k 1 ( μ 2 λ 2 ) k 2 C A ,
for every 0 < ξ < 1 , and some C ^ 1 , A ^ 1 > 0 . The definition of A , C in (40) yields
Ψ ¯ ( τ n ) C ^ 2 A ^ 2 n n ! k 1 ( μ 2 λ 2 ) k 2 k 1 λ 1 2 q n 2 7 9 k 1 λ 1 2 .
Statement (38) follows directly from the last inequality. □
Theorem 6.
Let the assumptions of Theorem 2 hold. Let F 1 be the Banach space of holomorphic and bounded functions on T 1 × χ 2 × H β . For all 0 h 1 ι 1 1 , we consider the function
ϵ u h 1 ( t 1 , x 2 ϵ μ 2 e 1 θ h 1 , z , ϵ ) , ϵ E h 1 ,
which defines a holomorphic and bounded function on E h 1 , with values in F 1 . Then, there exist a formal series u ^ ( ϵ ) F 1 [ [ ϵ ] ] such that for all 0 h 1 ι 1 1 , the function (43) admits u ^ ( ϵ ) as its ( ( q , K ) ; S ) -Gevrey asymptotic expansion on E h 1 , for K and S stated in (39).
Proof. 
Let 0 h 1 ι 1 1 . Regarding (32) in Theorem 2, Lemma 3 and from usual estimates, we guarantee that, for every ϵ E h 1 E h 1 + 1 ,
sup t 1 T 1 , x 2 χ 2 , z H β | u h 1 + 1 ( t , z , ϵ ) u h 1 ( t , z , ϵ ) | Ψ 1 | ϵ | | ϵ | n , n 0 .
Taking into account Lemma 3, we get that
u h 1 + 1 ( t , z , ϵ ) u h 1 ( t , z , ϵ ) F 1 C ^ A ^ n n ! S q n ( n + 1 ) K | ϵ | n , n 0 , ϵ E h 1 E h 1 + 1 .
Theorem 4 states the existence of a formal power series u ^ ( ϵ ) F 1 [ [ ϵ ] ] , which is the common ( ( q , K ) ; S ) -Gevrey asymptotic expansion of the function (43), as a function on E h 1 with values in F 1 , for every 0 h 1 ι 1 1 . □
Theorem 7.
Let the assumptions of Theorem 3 hold. Let F 2 be the Banach space of holomorphic and bounded functions on T 1 × ( T 2 D ( 0 , ρ 2 ) ) × H β . For every 0 h 2 ι 2 1 , we consider the function
ϵ u h 2 ( t , z , ϵ ) , ϵ E h 2 0 ,
which is an outer solution of (17), holomorphic and bounded on E h 2 0 , with values in F 2 . Then, there exists a formal power series u ^ 0 ( ϵ ) F 2 [ [ ϵ ] ] such that the function (44) admits u ^ 0 ( ϵ ) as its q-Gevrey asymptotic expansion of order 1 / ( k 1 λ 1 2 ) on E h 2 0 , for all 0 h 2 ι 2 1 .
Proof. 
Taking into account (36) in Theorem 3, and Theorem 5, we get the existence of a formal power series u ^ 0 ( ϵ ) F 2 [ [ ϵ ] ] which is the common q-Gevrey asymptotic expansion of the function (44), as a function on E h 2 0 with values in F 2 . This holds for all 0 h 2 ι 2 1 . □

7. Conclusions

The present study deals with the asymptotic behavior of the solutions of a family of singularly perturbed q-difference-differential equations in the complex domain. The main results (Theorems 6 and 7) describe different q-Gevrey levels which approximate the analytic solution (constructed in Theorem1) with respect to the perturbation parameter near the origin. This variation is obtained by means of the consideration of different domains in time variables where the analytic solutions are defined.
The appearance of the secondary branch of Lambert W function in the procedure puts into light the asymptotic behavior of the solutions of the equations considered, and the growth rate of its solutions near the singularity. In this sense, the study of such special function has been crucial in our reasoning. From the authors’ point of view, this outstanding behavior will give rise to further asymptotic results in other related q-difference problems that frequently appear in different applications.

8. Proof of Proposition 7

In this section, we give proof of the technical Proposition 7. We consider that U γ ( T , m , ϵ ) is constructed in the form (23).
For the first part of the proof, we take into account the property (6) on Jacobi Theta function, and Proposition 6, to arrive at
| U γ ( T , m , ϵ ) | 1 π q 1 / k 1 0 | ω ( r e γ 1 ) | 1 Θ q 1 / k 1 r e γ 1 T 1 e r | T 2 | k 2 cos ( γ k 2 arg ( T 2 ) ) d r r ϖ π q 1 / k 1 C q , k 1 δ ˜ ( 1 + | m | ) μ e β | m | I ( | T 1 | , | T 2 | ) ,
where
I ( | T 1 | , | T 2 | ) = 0 exp k 1 2 log ( q ) log 2 ( r + δ ) + α log ( r + δ ) exp k 1 2 log ( q ) log 2 r | T 1 | × | T 1 | r 1 / 2 exp r | T 2 | k 2 δ 3 d r r .
We split I ( | T 1 | , | T 2 | ) into the sum of I 1 ( | T 1 | , | T 2 | ) and I 2 ( | T 1 | , | T 2 | ) , where the first element is associated with the integration in ( 0 , ρ ) and the second is concerned with the integration restricted to ( ρ , ) , for some ρ > 0 . We study each part of the splitting:
We have
I 1 ( | T 1 | , | T 2 | ) C 3 | T 1 | 1 / 2 0 ρ exp k 1 2 log ( q ) log 2 r | T 1 | d r r 3 / 2 ,
for some C 3 > 0 , which, after the change of variable r = | T 1 | r , equals
C 3 0 ρ / | T 1 | exp k 1 2 log ( q ) log 2 ( r ) d r ( r ) 3 / 2 ,
which is bounded for every | T 1 | R + .
On the other hand, we assume ρ , ρ 1 > 0 are such that ρ ρ 1 . Then, the positive function ϕ : [ ρ , ) R defined by ϕ ( r ) = log 2 ( r / | T 1 | ) is monotone increasing on [ ρ , ) for any choice of | T 1 | < ρ 1 . Therefore,
I 2 ( | T 1 | , | T 2 | ) 1 ρ 3 / 2 | T 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ | T 1 | I 2.1 ( | T 2 | k 2 / δ 3 ) ,
where
I 2.1 ( | T 2 | k 2 / δ 3 ) = 0 exp k 1 2 log ( q ) log 2 ( r + δ ) + α log ( r + δ ) exp r δ 3 | T 2 | k 2 d r .
Let x = | T 2 | k 2 / δ 3 . We make the change of variable r = x r ˜ in the last integral and arrive at
I 2.1 ( x ) = 0 exp k 1 2 log ( q ) log 2 ( x r ˜ + δ ) + α log ( x r ˜ + δ ) exp r ˜ x d r ˜ .
Taking into account that
log 2 ( x r ˜ + δ ) = log 2 ( x ) + 2 log ( x ) log ( r ˜ + δ x ) + log 2 ( r ˜ + δ x ) ,
log ( x r ˜ + δ ) = log ( x ) + log ( r ˜ + δ x ) ,
we get that I 2.1 ( x ) = I 2.2 ( x ) + I 2.3 ( x ) , where the splitting is done on the integral by cutting the integration path into ( 0 , 1 δ / x ) and ( 1 δ / x , ) for I 2.2 ( x ) and I 2.3 ( x ) , respectively.
I 2.2 ( x ) exp k 1 2 log ( q ) log 2 ( x ) + α log ( x ) x 0 1 δ / x exp k 1 2 log ( q ) log 2 ( r ˜ + δ x ) exp r ˜ d r ˜ C 4 exp k 1 2 log ( q ) log 2 ( x ) + α log ( x ) x
for some C 4 > 0 which does not depend on x. Concerning I 2.3 ( x ) , we proceed analogously to arrive at
I 2.3 e k 1 2 log ( q ) log 2 ( x ) + α log ( x ) x 1 δ x r ˜ + δ x k 1 log ( x ) log ( q ) + α exp k 1 2 log ( q ) log 2 ( r ˜ + δ x ) e r ˜ d r ˜ C 5 e k 1 2 log ( q ) log 2 ( x ) + α log ( x ) x 1 δ x r ˜ + δ x k 1 log ( x ) log ( q ) e r ˜ 2 d r ˜ , C 5 e δ 2 x e k 1 2 log ( q ) log 2 ( x ) + α log ( x ) x 2 k 1 log ( x ) log ( q ) + 1 1 / 2 h k 1 log ( x ) log ( q ) e h d h C 5 e δ 2 x e k 1 2 log ( q ) log 2 ( x ) + α log ( x ) x 2 k 1 log ( x ) log ( q ) + 1 Γ k 1 log ( x ) log ( q ) + 1
for some C 5 > 0 . In the last sequence of inequalities, we have made the change of variable r ˜ + δ x = 2 h . The application of Stirling formula Γ ( x + 1 ) x 2 π x 1 / 2 x x e x leads us to
I 2.3 C 6 ( log ( x ) ) 1 / 2 k 1 log ( x ) log ( q ) k 1 log ( x ) log ( q ) e k 1 2 log ( q ) log 2 ( x ) + α log ( x ) x ,
for some C 6 > 0 . From (46)–(48), we conclude
I 2 ( | T 1 | , | T 2 | ) C 7 | T 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ | T 1 | | T 2 | k 2 δ 3 e k 1 2 log ( q ) log 2 ( | T 2 | k 2 δ 3 ) + α log ( | T 2 | k 2 δ 3 ) + | T 2 | k 2 δ 3 e k 1 2 log ( q ) log 2 | T 2 | k 2 δ 3 + α log | T 2 | k 2 δ 3 log ( | T 2 | k 2 δ 3 ) 1 / 2 e k 1 log ( q ) log | T 2 | k 2 δ 3 log k 1 log ( q ) log | T 2 | k 2 δ 3
for some C 7 > 0 . The first statement of Proposition 7 holds. We give the proof for the second statement.
The first arguments in the proof of the first statement can be followed word by word up to the splitting of I ( | T 1 | , | T 2 | ) into I 1 ( | T 1 | , | T 2 | ) + I 2 ( | T 1 | , | T 2 | ) . The quantity I 1 ( | T 1 | , | T 2 | ) is upper bounded by a constant for every | T 1 | > 0 and | T 2 | > 0 . We now proceed to give upper estimates on I 2 ( | T 1 | , | T 2 | ) . Let us choose ρ 1 , ρ > 0 such that ϕ ( r ) is monotone increasing on [ ρ , ) . It holds that
I 2 ( | T 1 | , | T 2 | ) 1 ρ 3 / 2 | T 1 | 1 / 2 exp k 1 2 log ( q ) log 2 ρ | T 1 | I 2.1 ( | T 2 | ) ,
where
I 2.1 ( | T 2 | ) = ρ exp k 1 2 log ( q ) log 2 ( r + δ ) + α log ( r + δ ) exp r δ 3 2 | T 2 | k 2 exp r δ 3 2 | T 2 | k 2 d r . C 8 exp ρ δ 3 2 | T 2 | k 2 ,
for some C 8 > 0 . The conclusion follows from this last upper bound.

9. Connection of the Solutions of (22) and (24)

Let ϵ D ( 0 , ϵ 0 ) . We consider ω ( τ , m , ϵ ) q Exp ( k 1 , β , μ , α ) d . Let γ R be an argument of S d .
We depart from
U γ ( T , z , ϵ ) = 1 π q 1 / k 1 L γ ω ( u , m , ϵ ) 1 Θ q 1 / k 1 u T 1 e 1 T 2 k 2 u d u u .
Let T ˜ 1 be a bounded sector and T ˜ 2 be an unbounded sector, both with a vertex at the origin, satisfying the assumptions in Section 5. U γ ( T , m , ϵ ) turns out to be a continuous function on T ˜ 1 × T ˜ 2 × R , holomorphic with respect to T on T ˜ 1 × T ˜ 2 .
Proposition 8.
In the previous framework, the following identities hold for all T 1 T ˜ 1 , T 2 T ˜ 2 , m R :
( i )
For every integer δ 0 , one has
( T 2 k 2 + 1 T 2 ) δ U γ ( T , m , ϵ ) = 1 π q 1 / k 1 L γ ( k 2 u ) δ ω ( u , m , ϵ ) 1 Θ q 1 / k 1 u T 1 e 1 T 2 k 2 u d u u .
( i i )
For every 1 1 D 1 1 , one has
T 1 d 1 σ q ; T 1 δ 1 U γ ( T , m , ϵ ) = 1 π q 1 / k 1 L γ u d 1 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 ω ( q δ 1 d 1 k 1 u , m , ϵ ) e 1 T 2 k 2 q δ 1 d 1 k 1 u 1 Θ q 1 / k 1 u T 1 d u u .
( i i i )
It holds that
T 1 d D 1 σ q ; T 1 d D 1 k 1 U γ ( T , m , ϵ ) = 1 π q 1 / k 1 L γ u d D 1 ( q 1 / k 1 ) d D 1 ( d D 1 1 ) / 2 ω ( u , m , ϵ ) e 1 T 2 k 2 u 1 Θ q 1 / k 1 u T 1 d u u .
( i v )
For every 1 1 D 1 1 and 1 2 D 2 1 , we have
T 1 d 1 σ q ; T 1 δ 1 σ q ; T 2 1 k 2 δ 1 d 1 k 1 U γ ( T , m , ϵ ) = 1 π q 1 / k 1 L γ u d 1 ( q 1 / k 1 ) d 1 ( d 1 1 ) / 2 ω ( q δ 1 d 1 k 1 u , m , ϵ ) e 1 T 2 k 2 u 1 Θ q 1 / k 1 u T 1 d u u .
Proof. 
Statement ( i ) is a consequence of the derivation under the integral sign. Statement ( i i ) is a direct consequence of Proposition 6 in [17] with k : = k 1 , σ : = d 1 , j : = δ 1 and f ( x ) = ω ( x , m , ϵ ) exp ( 1 T 2 k 2 x ) . The third statement is a consequence of Proposition 6 in [17]. The last statement is derived from the application of ( i i ) and the fact that
exp T 2 q 1 k 2 ( δ 1 d 1 k 1 ) k 2 q δ 1 d 1 k 1 u = exp T 2 k 2 u .

Author Contributions

Both authors have contributed equally for the work, concerning conceptualization, methodology, validation, formal analysis, investigation, resources, writing–original draft preparation, writing– review and editing, visualization, supervision, and project administration.

Funding

Both authors are are partially supported by the Spanish Ministry of Economy, Industry and Competitiveness under the project MTM2016-77642-C2-1-P.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lastra, A.; Malek, S. Boundary layer expansions for initial value problems with two complex time variables. arXiv 2019, arXiv:1904.04886. [Google Scholar]
  2. Lastra, A.; Malek, S. On parametric Gevrey asymptotics for some nonlinear initial value problems in two complex time variables. arXiv 2018, arXiv:1806.04615. [Google Scholar]
  3. Lastra, A.; Malek, S. On parametric Gevrey asymptotics for some initial value problems in two asymmetric complex time variables. Results Math. 2018, 73, 155. [Google Scholar] [CrossRef]
  4. O’Malley, R. Singular Perturbation Methods for Ordinary Differential Equations; Applied Mathematical Sciences; Springer: New York, NY, USA, 1991; Volume 89, viii+225p. [Google Scholar]
  5. Skinner, L. Singular Perturbation Theory; Springer: New York, NY, USA, 2011. [Google Scholar]
  6. Fruchard, A.; Schäfke, R. Composite Asymptotic Expansions; Lecture Notes in Mathematics; Springer: Heidelberg, Germany, 2013; Volume 2066, x+161p. [Google Scholar]
  7. Tahara, H. q-analogues of Laplace and Borel transforms by means of q-exponentials. Ann. Inst. Fourier 2017, 67, 1865–1903. [Google Scholar] [CrossRef]
  8. Tahara, H.; Yamazawa, H. q-analogue of summability of formal solutions of some linear q-difference-differential equations. Opusc. Math. 2015, 35, 713–738. [Google Scholar] [CrossRef]
  9. Yamazawa, H. Holomorphic and singular solutions of a q-analog of the Briot–Bouquet type difference-differential equations. Funkc. Ekvacioj 2016, 59, 185–197. [Google Scholar] [CrossRef]
  10. Pravica, D.W.; Spurr, M.J. Unique summing of formal power series s olutions to advanced and delayed differential equations. Discret. Contin. Dyn. Syst. 2005, 2005, 730–737. [Google Scholar]
  11. Ho, C.-L. On the use of Mellin transform to a class of q-difference-differential equations. Phys. Lett. A 2000, 268, 217–223. [Google Scholar] [CrossRef]
  12. Pravica, D.W.; Randriampiry, N.; Spurr, M.J. q-advanced models for tsunami and rogue waves. Abstr. Appl. Anal. 2012, 2012, 414060. [Google Scholar] [CrossRef]
  13. Pravica, D.W.; Randriampiry, N.; Spurr, M.J. Solutions of a class of multiplicatively advanced differential equations. C. R. Math. Acad. Sci. 2018, 356, 776–817. [Google Scholar] [CrossRef]
  14. Pravica, D.W.; Randriampiry, N.; Spurr, M.J. On q-advanced spherical Bessel functions of the first kind and perturbations of the Haar wavelet. Appl. Comput. Harmon. Anal. 2018, 44, 350–413. [Google Scholar] [CrossRef]
  15. Dreyfus, T. Building meromorphic solutions of q-difference equations using a Borel–Laplace summation. Int. Math. Res. Not. IMRN 2015, 2015, 6562–6587. [Google Scholar] [CrossRef]
  16. Ramis, J.-P. About the growth of entire functions solutions of linear algebraic q-difference equations. Ann. Fac. Sci. Toulouse Math. 1992, 1, 53–94. [Google Scholar] [CrossRef]
  17. Malek, S. On parametric Gevrey asymptotics for a q-analog of some linear initial value problem. Funkc. Ekvacioj 2017, 60, 21–63. [Google Scholar] [CrossRef]
  18. Vizio, L.D.; Zhang, C. On q-summation and confluence. Ann. Inst. Fourier 2009, 59, 347–392. [Google Scholar] [CrossRef]
  19. Lastra, A.; Malek, S. On q-Gevrey asymptotics for singularly perturbed q-difference-differential problems with an irregular singularity. Abstr. Appl. Anal. 2012, 2012, 860716. [Google Scholar] [CrossRef]
  20. Lastra, A.; Malek, S. On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations. Adv. Differ. Equ. 2015, 2015, 344. [Google Scholar] [CrossRef]
  21. Lastra, A.; Malek, S. On singularly perturbed linear initial value problems with mixed irregular and fuchsian time singularities. J. Geom. Anal. 2019. [Google Scholar] [CrossRef]
  22. Costin, O.; Tanveer, S. Short time existence and Borel summability in the Navier–Stokes equation in R 3. Comm. Partial Differ. Equ. 2009, 34, 785–817. [Google Scholar] [CrossRef]
  23. Dreyfus, T.; Lastra, A.; Malek, S. On the multiple scale analysis for some linear partial q-difference and differential equations with holomorphic coefficients. Adv. Differ. Equ. 2019. [Google Scholar] [CrossRef]
  24. Malek, S. On a Partial q-Analog of a Singularly Perturbed Problem with Fuchsian and Irregular Time Singularities. 2019. Available online: https://web.ma.utexas.edu/mp_arc/c/19/19-33.pdf (accessed on 8 July 2019).
  25. Balser, W. From Divergent Power Series to Analytic Functions. Theory and Application of Multisummable Power Series; Lecture Notes in Mathematics; Springer: Berlin, Germany, 1994; Volume 1582, x+108p. [Google Scholar]
  26. Hsieh, P.; Sibuya, Y. Basic Theory of Ordinary Differential Equations; Universitext; Springer: New York, NY, USA, 1999. [Google Scholar]
  27. Chatzigeorgiou, I. Bounds on the Lambert Function and Their Application to the Outage Analysis of User Cooperation. IEEE Commun. Lett. 2013, 17, 1505–1508. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Initial (left) and deformed path (right) for the difference of two consecutive solutions of (17). The symbol “×” represents a root of P m ( τ ) .
Figure 1. Initial (left) and deformed path (right) for the difference of two consecutive solutions of (17). The symbol “×” represents a root of P m ( τ ) .
Mathematics 07 00924 g001

Share and Cite

MDPI and ACS Style

Lastra, A.; Malek, S. On a q—Analog of a Singularly Perturbed Problem of Irregular Type with Two Complex Time Variables. Mathematics 2019, 7, 924. https://doi.org/10.3390/math7100924

AMA Style

Lastra A, Malek S. On a q—Analog of a Singularly Perturbed Problem of Irregular Type with Two Complex Time Variables. Mathematics. 2019; 7(10):924. https://doi.org/10.3390/math7100924

Chicago/Turabian Style

Lastra, Alberto, and Stéphane Malek. 2019. "On a q—Analog of a Singularly Perturbed Problem of Irregular Type with Two Complex Time Variables" Mathematics 7, no. 10: 924. https://doi.org/10.3390/math7100924

APA Style

Lastra, A., & Malek, S. (2019). On a q—Analog of a Singularly Perturbed Problem of Irregular Type with Two Complex Time Variables. Mathematics, 7(10), 924. https://doi.org/10.3390/math7100924

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