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 .
More precisely, we consider equations of the form
under null initial data
. Here,
, and
P stands for a polynomial with complex coefficients with respect to
, a polynomial with rational powers with respect to
and
, holomorphic with respect to
z on a horizontal strip
for some
, and holomorphic with respect to the perturbation parameter
on a small disc centered at the origin, say
for some
. Throughout the present work,
stands for the dilation operator on
t variable for some fixed
, i.e.,
We adopt the following notation , for any .
The forcing term is constructed under certain growth conditions and turns out to be a holomorphic function in , for .
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],
under null initial data
, and where
, and
P is a polynomial with respect to its first three variables, with holomorphic coefficients on
, and the forcing term is holomorphic on
. 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
(see (
21) and (
23)), for some appropriate
and
(see
Section 3 for their definition). The function
is obtained from a fixed point argument (see Proposition 6), and belongs to
, 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
, and a second on domains of the form
, where
is a finite sector,
is an unbounded sector and where
is a bounded sector that depends on
, and tends to infinity with
approaching the origin. The sets
and
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 -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
, and consider the holomorphic solutions of the problem (
1),
defined on
w.r.t.
for all
. In Theorem 6, we prove that, for some
,
and some adequate domain
, the function
with values in the Banach space of holomorphic and bounded functions on
, say
, admits a formal power series
as
-Gevrey asymptotic expansion on
for some
, for every
.
The proof of this result leans on the application of accurate estimates related to the -branch of Lambert W function (see Lemma 3).
Concerning the
outer solutions of the main problem under study, we consider a good covering
, and the holomorphic solutions of (
1),
defined on
w.r.t.
for all
. In Theorem 7, we prove that
is an outer solution of (
1) with values in the Banach space
of holomorphic and bounded functions on
. Moreover, there exists a formal power series
, which is the common
q-Gevrey asymptotic expansion of some positive order of each solution (
5) on
, for
.
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.
3. Statement of the Main Problem and Auxiliary Equations
Let
and
with
, and
be integer numbers. We also consider a real number
, and choose real numbers
, and
for every
and
. We assume that
In addition to that, there exist natural numbers
with
Let us also fix polynomials
, and
for all
and
, with complex coefficients such that
where
for some
and
. Observe that this condition implies
. These polynomials are chosen in such a way that
for
, and
The main problem under study in this work is
under initial data
.
The function
f is a holomorphic function in
for every
, where
stands for the horizontal strip
. It is constructed as follows. Let
be a continuous function, continuous in
, for some
, entire with respect to the first variable, and holomorphic with respect to the third one in
. We moreover assume that there exists
such that
for every
, some
,
, and where
satisfies (
16).
In view of the definition of
q-Laplace transform and the results described in
Section 2.1, one can define
for some fixed
. Here,
. The function
turns out to be holomorphic on the set
. Observe that, in the previous construction, given
, one can choose
with
and
for all
and with
.
For every
and
, the function
is constructed in the following way:
turns out to be a holomorphic function on
whenever
. In addition to that, we assume that uniform bounds with respect to the perturbation parameter are satisfied, i.e., there exist
with
As mentioned in the title of the work, Equation (
17) has the property to be singularly perturbed at
. Namely, under the above hypotheses, when
equals 0, the nature of the equation changes drastically and is reduced to an ODE with constant coefficients
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
Assuming the solution is of the form (
21), the expression
solves
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
for some appropriate function
and
. We refer to
Section 2.1 for the definitions of the elements involved in the previous expression.
Let us consider a second auxiliary equation:
We define the polynomial
by
whose factorization is given by
with
for all
, where
.
Let
be such that the infinite sector
of bisecting direction
d satisfies the following geometric construction: there exists
such that
, for all
,
and
. The previous is a feasible condition for an appropriate choice of small enough
and large enough
, in view of (
14) and the definition of
. More precisely, one chooses
such that
has positive distance to 1, for all
,
, and
. The previous choice of
d yields
for some
, valid for all
and all
.
Proposition 6. Let be an infinite sector with vertex at the origin satisfying the previous geometric conditions. Then, if (see (18) for its definition) and are small enough, there exists such that Equation (24) admits a unique solution which belongs to with , for all . Proof. Let
. We consider the map
defined by
Given
, let
with
. In view of (
11)–(
13), (
26) and from Lemma 2, Propositions 4 and 5, one has
for every
and
. In addition to this, Lemma 2 and (
26) yield
Let
be small enough satisfying
Then, the estimates (
27) and (
28) yield to
. In other words, the operator
restricted to
is such that
.
On the other hand, let
. Analogously to (
27), one arrives at
by choosing
such that
We conclude that the map
is contractive. The classical fixed point theory in complete metric spaces states the existence of a unique fixed point for
, say
, in
, with
. For every
, the function
is a solution of (
24) in view of the definition of the operator
. Holomorphy of the map
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.
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
and
be an infinite sector with vertex at
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 (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 . Let us first recall the notion of good covering in .
Definition 5. Let be an integer. For every , we choose a finite sector with vertex at the origin such that:
, and if and only if with (under the convention that ).
, for some neighborhood of the origin .
A family of sectors under these assumptions is known as a good covering in .
Let
and
be sectors following the construction in
Section 4.
Definition 6. Let be a good covering in . We also fix a bounded sector, , and an unbounded sector , both with vertex at the origin. For all , let be an infinite sector of bisecting direction . We say that the set is admissible if the following conditions hold:
For every , and we have .
For every , and we have .
Observe that given an admissible set
, the sectors
, the choice of the sectors
,
and
, fixed in
Section 2.1, entail the existence of
with
In addition to this, it holds that
for some
and
, for some
,
,
.
Definition 7. Let be an admissible set.
Then, is known as a family of sectors associated with the good covering .
Theorem 1. Let be a good covering in . For every , we choose such that satisfies the geometric conditions of Proposition 6. Let be a family of sectors associated with the good covering . If there exist small enough such that, ifthen, for every , the problem (17) admits a solution , which defines a bounded and holomorphic function in , for any fixed . Proof. Let
. Proposition 6 guarantees the existence of
such that the Equation (
24) admits a unique solution
which belongs to
for every
, and the map
is holomorphic in
. Taking into account that
is associated with the good covering
and the properties of Laplace transform stated in
Section 2.1, one can construct
in the form (
23), which is well defined, bounded and continuous function on
, where
for
, and
is an infinite sector. The function
is holomorphic w.r.t.
on
. As a matter of fact,
is a solution of (
22) in view of the properties relating both equalities in
Section 9. We finally define
following (
21):
which turns out to be a holomorphic solution of (
17), defined on
, for any fixed
. □
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 . Let be a good covering of . We also consider the admissible set . Let be a natural number satisfying Let be a bounded domain, such that the good covering satisfies the following condition: for all we can select (which depends on ) such that, for every and , the complex number belongs to . We define the set .
In case , , , for , then we say that 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 be a good covering of . Then, there exists such that, for every , , , for any fixed and , one hasfor some , and a real number chosen to be small enough. Proof. Let
and consider consecutive solutions
of (
17), constructed in Theorem 1. We recall that the function
stands for the analytic continuation of a function
, holomorphic w.r.t.
in a neighborhood of the origin
, to the infinite sector
. This entails that the difference
can be written in the following form, after an appropriate path deformation in
, avoiding the roots of
(see
Figure 1).
We write
, where
where
for
, and
stands for the arc of circle from
to
.
Assume that
and
,
for some
and
. Owing to the estimates leading to (
29), displayed in (
49), we derive that
for some large enough
. Moreover, if
is small enough, then there exists
such that
for every
,
,
, and
. An analogous upper bound is attained for
. Concerning
, one can apply the bounds stated in Proposition 6, and (
6) to get the existence of
such that
The result follows taking into account (
34) and (
35), under the condition (
31). □
5.2. Outer Solutions of the Main Problem
Definition 9. Let and be a good covering of , and consider an admissible set . Assume that is such that for some fixed , which is independent of ϵ, then we say that 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 be a good covering of . Then, there exists such that, for every , , with , for any fixed and , one hasfor some and . Proof. Let
and consider consecutive solutions
of (
17), constructed in Theorem 1. We proceed to write the difference of two consecutive solutions in the form
for all
, with
substituted by
in the expressions of
in the proof of Theorem 2. Let
and
, with
,
for some
and
. Owing to analogous bounds as those leading to (
50) and (
51), we arrive at
for some
. Similar estimates hold for
. On the other hand, direct computations yield
for some
. This concludes the proof. □
8. Proof of Proposition 7
In this section, we give proof of the technical Proposition 7. We consider that
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
where
We split into the sum of and , where the first element is associated with the integration in and the second is concerned with the integration restricted to , for some . We study each part of the splitting:
We have
for some
, which, after the change of variable
, equals
which is bounded for every
.
On the other hand, we assume
are such that
. Then, the positive function
defined by
is monotone increasing on
for any choice of
. Therefore,
where
Let
. We make the change of variable
in the last integral and arrive at
Taking into account that
we get that
, where the splitting is done on the integral by cutting the integration path into
and
for
and
, respectively.
for some
which does not depend on
x. Concerning
, we proceed analogously to arrive at
for some
. In the last sequence of inequalities, we have made the change of variable
. The application of Stirling formula
leads us to
for some
. From (
46)–(
48), we conclude
for some
. 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
into
. The quantity
is upper bounded by a constant for every
and
. We now proceed to give upper estimates on
. Let us choose
such that
is monotone increasing on
. It holds that
where
for some
. The conclusion follows from this last upper bound.