On a parabolic partial differential equation and system modeling a production planning problem

1.   Introduction

Several years ago, Bensoussan, Sethi, Vickson and Derzko ^[1] have been considered the case of a factory producing one type of economic goods and observed that it is necessary to solve the simple partial differential equation

where $\sigma \in \left(0, \infty \right)$ denotes the diffusion coefficient, $\alpha \in \left[ 0, \infty \right)$ represents psychological rate of time discount, $x\in \mathbb{R}^{N}$ is the product vector, $z: = z_{s}^{\alpha }\left(x\right)$ denotes the value function and $\left\vert x\right\vert ^{2}$ is the loss function.

Regime switching refers to the situation when the characteristics of the state process are affected by several regimes (e.g., in finance bull and bear market with higher volatility in the bear market).

It is important to point out that, when dealing with regime switching, we can describe a wide variety of phenomena using partial differential equations. In ^[1], the authors Cadenillas, Lakner and Pinedo ^[2] adapted the model problem in ^[1] to study the optimal production management characterized by the two-state regime switching with limited/unlimited information and corresponding to the system

where $\sigma _{1}, \sigma _{2}\in \left(0, \infty \right)$ denote the diffusion coefficients, $\alpha _{1}, \alpha _{2}\in \left[ 0, \infty \right)$ represent the psychological rates of time discount from what place the exponential discounting, $x\in \mathbb{R}^{N}$ is the product vector, $u_{r}^{s}: = u_{r}^{s}\left(x\right)$ ($r = 1, 2$) denotes the value functions, $\left\vert x\right\vert ^{2}$ is the loss function, $\rho \in \left[ -1, 1 \right]$ is the correlation coefficient and $a_{nm}$ ($n, m = 1, 2$) are the elements of the Markov chain's rate matrix, denoted by $G = [\vartheta _{nm}]_{2\times 2}$ with

the diagonal elements $\vartheta _{nn}$ may be expressed as $\vartheta _{nn} = -\underset{m\neq n}{\Sigma }\vartheta _{nm}$.

Furthermore, in civil engineering, Dong, Malikopoulos, Djouadi and Kuruganti ^[3] applied the model described in ^[2] to the study of the optimal stochastic control problem for home energy systems with solar and energy storage devices; the two regimes switching are the peak and the peak energy demands.

After that, there have been numerous applications of regime switching in many important problems in economics and other fields, see the works of: Capponi and Figueroa-López ^[4], Elliott and Hamada ^[5], Gharbi and Kenne ^[6], Yao, Zhang and Zhou ^[7] and Wang, Chang and Fang ^[8] for more details. Other different research studies that explain the importance of regime switching in the real world are ^[9,10].

In this paper, we focus on the following parabolic partial differential equation and system, corresponding to (1.1)

and (1.2) respectively

where $z_{s}^{\alpha }$ is the solution of (1.1) and $\left(u_{1}^{s}\left(x\right), u_{2}^{s}\left(x\right) \right)$ is the solution of (1.2). The existence and the uniqueness for the case of (1.1) is proved by ^[10] and the existence for the system case of (1.2) by ^[11].

From the mathematical point of view the problem (1.3) has been extensively studied when the space $\mathbb{R}^{N}$ is replaced by a bounded domain and when $\alpha = 0$. In particular, some great results can be found in the old papers of Barles, Porretta ^[12] and Tchamba ^[13]. More recently, but again for the case of a bounded domain, $\alpha = 0$ and in the absence of the gradient term, the problem (1.3) has been also discussed by Alves and Boudjeriou ^[14]. The interest of these authors ^[12,13,14] is to give an asymptotic stable solution at infinity for the considered equation, i.e., a solution which tends to the stationary Dirichlet problem associated with (1.3) when the time go to infinity.

Next, we propose to find a similar result as of ^[12,13,14], for the case of equation (1.3) and system (1.4) that model some real phenomena. More that, our first interest is to provide a closed form solution for (1.3) and (1.4). Our second objective is inspired by the paper of ^[14,15], and it is to solve the parabolic partial differential equation

where $T < \infty$ and $B_{R}$ is a ball of radius $R > 0$ with origin at the center of $\mathbb{R}^{N}$.

Let us finish our introduction and start with the main results.

2.   The main results

We use the change of variable

in

to rewrite (1.3) and (1.5) in an equivalent form

where

Our first result is the following.

Theorem 2.1. Assume $\Omega = B_{R}$, $N\geq 3$, $T < \infty$ and $\alpha = 0$.There exists a unique radially symmetric positive solution

of (2.2) increasing in the time variable and such that

where $u_{s}\in C^{2}\left(B_{R}\right) \cap C\left(\overline{B}_{R}\right)$ is the unique positive radially symmetric solution of theDirichlet problem

which will be proved. In addition,

is the unique radially symmetric solution of the problem (1.5).

Instead of the existence results discussed in the papers of ^[12,13,14], in our proof of the Theorem 2.1 we give the numerical approximation of solution $u\left(x, t\right)$.

The next results refer to the entire Euclidean space $\mathbb{R}^{N}$ and present closed-form solutions.

Theorem 2.2. Assume $\Omega = \mathbb{R}^{N}$, $N\geq 1$, $T = \infty$, $\alpha > 0$ and $c\in \left(0, \infty \right)$ is fixed. There exists aunique radially symmetric solution

of (2.2), increasing in the time variable and such that

where $u_{s}^{\alpha }\in C^{2}\left(\mathbb{R}^{N}\right)$ is the uniqueradially symmetric solution of the stationary Dirichlet problem associatedwith (2.2)

Moreover, the closed-form radially symmetric solution of the problem (1.3) is

where

The following theorem is our main result regarding the system (1.4).

Theorem 2.3. Suppose that $N\geq 1$, $\alpha _{1}, \alpha _{2}\in \left(0, \infty \right)$ and\ $a_{11}, a_{22}\in \left[ 0, \infty \right)$ with $a_{11}^{2}+a_{22}^{2}\neq 0$. Then, the system (1.4) has a uniqueradially symmetric convex solution

of quadratic form in the $x$ variable and such that

where

is the radially symmetric convex solution of quadratic form in the $x$variable of the stationary system (1.2) which exists from the resultof ^[11].

Our results complete the following four main works: Bensoussan, Sethi, Vickson and Derzko ^[1], Cadenillas, Lakner and Pinedo ^[2], Canepa, Covei and Pirvu ^[15] and Covei ^[10], which deal with a stochastic control model problem with the corresponding impact for the parabolic case (see ^[13,16] for details).

2.1.   An auxiliary result

To prove our Theorem 2.1, we use a lower and upper solution method and the comparison principle that can be found in ^[17].

Lemma 2.1. If, there exist $\overline{u}\left(x\right)$, $\underline{u}\left(x\right) \in C^{2}\left(B_{R}\right) \cap C\left(\overline{B}_{R}\right)$ two positive functions satisfying

then

and there exists

a solution of (2.4) such that

where $\underline{u}\left(x\right)$ and $\overline{u}\left(x\right)$ arerespectively, called a lower solution and an upper solution of (2.4).

The corresponding result of Lemma 2.1 for the parabolic equations can be found in the work of Pao ^[18] and Amann ^[19]. To achieve our goal, complementary to the works ^{[12,13,14,15]} it can be used the well known books of Gilbarg and Trudinger ^[20], Sattinger ^[17], Pao ^[18] and a paper of Amann ^[19]. Further on, we can proceed to prove Theorem 2.1.

3.   Proof of Theorem 2.1

By a direct calculation, if there exists and is unique, $u_{s}\in C^{2}\left(B_{R}\right) \cap C\left(\overline{B}_{R}\right)$, a positive solution of the stationary Dirichlet problem (2.4) then

is the solution of the problem (2.2) and

is the solution of the problem (1.5) belonging to

We prove that (2.4) has a unique radially symmetric solution. The existence of solution for (2.4) is obtained by a standard monotone iteration and the lower and the upper solution method, Lemma 2.1. Hence, starting from the initial iteration

we construct a sequence $\{u_{s}^{k}\left(x\right) \}_{k\geq 1}$ successively by

and this sequence will be pointwise convergent to a solution $u_{s}\left(x\right)$ of (2.4).

Indeed, since for each $k$ the right-hand side of (3.1) is known, the existence theory for linear elliptic boundary-value problems implies that $\{u_{s}^{k}\left(x\right) \}_{k\geq 1}$ is well defined, see ^[20].

Let us prove that $\{u_{s}^{k}\left(x\right) \}_{k\geq 1}$ is a pointwise convergent sequence to a solution of (2.4) in $\overline{B}_{R}$. To do this, first we prove that $\{u_{s}^{k}\left(x\right) \}_{k\geq 1}$ is monotone nondecreasing of $k$. We apply the mathematical induction by verifying the first step, $k = 1$.

Now, by the standard comparison principle, Lemma 2.1, we have

Moreover, the induction argument yields the following

i.e., $\{u_{s}^{k}\left(x\right) \}_{k\geq 1}$ is a monotone nondecreasing sequence.

Next, using again Lemma 2.1, we find

where we have used

i.e., Lemma 2.1 confirm.Thus, in view of the monotone and bounded property in (3.3) the sequence $\left\{ { u_{s}^{k}\left(x\right)} \right\}$$_{k\geq 1}$ converges. We may pass to the limit in (3.3) to get the existence of a solution

associated to (2.4), which satisfies

Furthermore, the convergence of $\left\{ {u_{s}^{k}\left(x\right)} \right\}$ is uniformly to $u_{s}\left(x\right)$ in $\overline{B}_{R}$ and $u_{s}\left(x\right)$ has a radial symmetry, see ^[15] for arguments of the proof. The regularity of solution $u_{s}\left(x\right)$ is a consequence of classical results from the theory of elliptic equations, see Gilbarg and Trudinger ^[20]. The uniqueness of $u_{s}\left(x\right)$ follows from a standard argument with the use of Lemma 2.1 and we omit the details.

Clearly, $u\left(x, t\right)$ is increasing in the time variable. The regularity of $u\left(x, t\right)$ follows from the regularity of $u_{s}\left(x\right)$. Letting $t\rightarrow T$ we see that (2.3) holds. The solution of the initial problem (1.5) is saved from (2.1).

Finally, we prove the uniqueness for (2.2). Let

be two solutions of the problem (2.2), i.e., its hold

and

Setting

and subtracting the two equations corresponding to $u$ and $v$ we find

Let us prove that $u\left(x, t\right) -v\left(x, t\right) \leq 0$ in $\overline{B}_{R}\times \left[ 0, T\right]$. If the conclusion were false, then the maximum of

is positive. Assume that the maximum of $w$ in $\overline{B}_{R}\times \left[ 0, T\right]$ is achieved at $\left(x_{0}, t_{0}\right)$. Then, at the point $\left(x_{0}, t_{0}\right) \in B_{R}\times \left[ 0, T\right)$, where the maximum is attained, we have

and

which is a contradiction. Reversing the role of $u$ and $v$ we obtain that $u\left(x, t\right) -v\left(x, t\right) \geq 0$ in $\overline{B}_{R}\times \left[ 0, T\right]$. Hence $u\left(x, t\right) = v\left(x, t\right)$ in $\overline{B}_{R}\times \left[ 0, T\right]$. The proof of Theorem 2.1 is completed.

Finally, our main result, Theorem 2.2 will be obtained by a direct computation.

4.   Proof of Theorem 2.2

In view of the arguments used in the proof of Theorem 2.1 and the real world phenomena, we use a purely intuitive strategy in order to prove Theorem 2.2.

Indeed, for the verification result in the production planning problem, we need $z\left(x, t\right)$ to be almost quadratic with respect to the variable $x$.

More exactly, we observe that there exists and is unique

that solve (2.2), where

and $B$, $D$ are given in (2.8). The condition (4.1) is used to obtain the asymptotic behaviour of solution to the stationary Dirichlet problem associated with (2.2). Then our strategy is reduced to find $B, D\in \left(0, \infty \right)$ and the function $h$ which depends of time and $c\in \left(0, \infty \right)$ such that

or, after rearranging the terms

where (4.1) holds. Now, by a direct calculation we see that the system of equations

has a unique solution that satisfies our expectations, namely,

where $B$ and $D$ are given in (2.8), is a radially symmetric solution of the problem (2.2). The uniqueness of the solution is followed by the arguments in ^[10] combined with the uniqueness proof in Theorem 2.1. The justification of the asymptotic behavior and regularity of the solution can be proved directly, once we have a closed-form solution. Finally, the closed-form solution in (2.7) is due to (2.1)–(4.2) and the proof of Theorem 2.2 is completed.

5.   Proof of Theorem 2.3

One way of solving this system of partial differential equation of parabolic type (1.4) is to show that the system (1.4) is solvable by

for some unique $\beta _{1}, \beta _{2}, \eta _{1}, \eta _{2}\in \left(0, \infty \right)$ and $h_{1}\left(t\right)$, $h_{2}\left(t\right)$ are suitable chosen such that

The main task for the proof of existence of (5.1) is performed by proving that there exist

such that

or equivalently, after grouping the terms

where $h_{1}\left(t\right),$ $h_{2}\left(t\right)$ must satisfy (5.2). Now, we consider the system of equations

To solve (5.3), we can rearrange those equations $1$, $2$ in the following way

We distinguish three cases:

1.in the case $a_{22} = 0$ we have an exact solution for (5.4) of the form

2.in the case $a_{11} = 0$ we have an exact solution for (5.4) of the form

3.in the case $a_{11}\neq 0$ and $a_{22}\neq 0,$ to prove the existence and uniqueness of solution for (5.4) we will proceed as follows. We retain from the first equation of (5.4)

and from the second equation

The existence of $\beta _{1}$, $\beta _{2}\in \left(0, \infty \right)$ for (5.4) can be easily proved by observing that the continuous functions $f_{1}, f_{2}:\left[ 0, \infty \right) \rightarrow \mathbb{R}$ defined by

have the following properties

respectively

The observations (5.5) and (5.6) imply

has at least one solution $\left(\beta _{1}, \beta _{2}\right) \in \left(0, \infty \right) \times \left(0, \infty \right)$ and furthermore it is unique (see also, the references [21,22] for the existence and the uniqueness of solutions).

The discussion from cases 1–3 show that the system (5.4) has a unique positive solution. Next, letting

be the unique positive solution of (5.4), we observe that the equations $3$, $4$ of (5.3) can be written equivalently as a system of linear equations that is solvable and with a unique solution

By defining

we observe that

Using the fact that $G_{a, \alpha }^{-1}$ has all ellements positive and rewriting (5.7) in the following way

we can see that there exist and are unique $\eta _{1}$, $\eta _{2}\in \left(0, \infty \right)$ that solve (5.7). Finally, the equations $5$, $6$, $7$ of (5.3) with initial condition (5.2) can be written equivalently as a solvable Cauchy problem for a first order system of differential equations

with a unique solution and then (5.1) solve (1.4). The rest of the conclusions are easily verified.

6.   Application for the system case

Next, we present an application.

Application 1. Suppose there is one machine producing two products (see ^[23,24], for details). We consider a continuous time Markov chain generator

and the time-dependent production planning problem with diffusion $\sigma _{1} = \sigma _{2} = \frac{1}{\sqrt{2}}$ and let $\alpha _{1} = \alpha _{2} = \frac{1 }{2}$ the discount factor. Under these assumptions, we can write the system (5.4) with our data

which has a unique positive solution

On the other hand, the system (5.7) becomes

which has a unique positive solution

Finally, the system in (5.8) becomes

which has the solution

Next, from

we have

and finally

from where we can write the unique solution of the system (1.4) in the form (5.1).

7.   Discussion

Let us point that in Theorem 2.3 we have proved the existence and the uniqueness of a solution of quadratic form in the $x$ variable and then the existence of other different types of solutions remain an open problem.

8.   Conclusions

Some closed-form solutions for equations and systems of parabolic type are presented. The form of the solutions is unique and tends to the solutions of the corresponding elliptic type problems that were considered.

Acknowledgments

The author is grateful to the anonymous referees for their useful suggestions which improved the contents of this article.

Conflict of interest

The authors declare there is no conflict of interest.