Discrete Dynamics in Nature and Society

Discrete Dynamics in Nature and Society / 2009 / Article

Research Article | Open Access

Volume 2009 |Article ID 827087 | 12 pages | https://doi.org/10.1155/2009/827087

On the Cauchy Problem of a Quasilinear Degenerate Parabolic Equation

Academic Editor: Guang Zhang
Received21 Apr 2009
Accepted25 Aug 2009
Published11 Nov 2009


By Oleinik's line method, we study the existence and the uniqueness of the classical solution of the Cauchy problem for the following equation in : , provided that is suitable small. Results of numerical experiments are reported to demonstrate that the strong solutions of the above equation may blow up in finite time.

1. Introduction

We consider the following Cauchy problem: This problem arises in financial mathematics recently; more and more mathematicians have been interested in it. In [1], Antonelli et al. introduced a new model for agents' decision under risk, in which the utility function is the solution to (1.1)-(1.2); they also proved, by means of probability methods, the existence of a continuous viscosity solution of (1.1)-(1.2), which satisfies for every , , under the assumption that is uniformly Lipschitz continuous function. In [2], Citti et al. studied the interior regularity properties of this problem; they proved that the viscosity solutions are indeed classical solutions. On the other hand, Antonelli and Pascucci [3] showed that the solution found in [1] can be also considered as a distributional solution.

However, all the above results are obtained when is suitably small; say, the solution is local. The global weak solutions of the Cauchy problem for a more general class of equations, that contains (1.1), are obtained in [47], and so forth. This kind of solutions, however, is few regular and does not satisfy condition (1.3) in general.

In this paper, we will solve the Cauchy problem (1.1)-(1.2) in another simpler way and get the result as [2] again. Moreover, some examples are provided by numerical computation. The results of computation show that the strong solutions of the above equation may blow-up in finite time, though there exist the global weak solutions.

2. Line Method

In order to describe our method, we have to quote the well-known Prandtl system for a nonstationary boundary layer arising in an axially symmetric incompressible flow past a solid body, it has the form in a domain , where and are given functions. If we introduce the Crocco variables: we obtain the following equation for : Oleinik and Samokhin [8] had done excellent work in the boundary theory by the line method. Comparing this equation with (1.1), it is natural to conjecture that we are able to solve problem (1.1)-(1.2) by Qleinik's method.

Consider the following initial boundary problem: where ; its first-order derivatives and are all bounded.

Definition 2.1. A function is said to be a solution of problem (2.4)-(2.5) if has first-order derivatives in (2.4) which is continuous in , and its derivative is continuous; satisfies (2.4) in , together with condition (2.5).

The solution of problem (2.4)-(2.5) will be constructed as the limit of a sequence , , which consists of solutions of the equations As we take a function which is smooth in .

Suppose that for some nonnegative number and when ,

Lemma 2.2. Let be a smooth function such that in , for . Then everywhere . Let be a smooth function such that in , for . Then everywhere in .

Proof. Let us prove the first statement of Lemma 2.2. The difference satisfies the inequality
Let . Then If we choose large enough, by the maximal principle, we know everywhere in .
Let us construct functions satisfying the conditions of Lemma 2.2. To this end, we define a twice continuously differentiable even function such that for , for , where is a function, .
When , if we chose large enough and small enough.
When , by the same reason.
Let , , . Then if we chose large enough and small enough.
Similarly, we are able to prove the second statement of Lemma 2.2.

Thus we have the following.

Lemma 2.3. Suppose that satisfies (2.9) and , then
The smooth functions , can be constructed as in [8], and we omit details here.

Let where . We will show that there exist positive constants and such that the conditions for , , imply that for .

First, we rewrite (2.6) as

Applying the operator to (2.17), then

By (2.9), (2.15), and Cauchy inequality, we are able to get where depends on and its derivatives are up to the second. Let with a positive constant to be chosen later. Then if we choose according to such that . If attains its positive maximum at , then where the constant does not depend on . At the same time, the positive maximum of in cannot be attained by maximal principle. Thus we have

So, if we let small enough such that and set then

In order to estimate the second derivatives of in , consider the function Applying the operator to both sides of (2.17), we find that

At the same time, we can calculate that and so we have By the introduced assumption that the first-order and second-order derivatives of , , and are all bounded and using Cauchy inequality, we can get from (2.29) that By the transformation , if we chose large enough, we are able to show that there exist positive constants and such that the conditions for , , imply that for . Thus we have the following.

Theorem 2.4. Let be the solutions of problems (2.4)-(2.5), then the derivatives of up to the second-order are uniformly bounded with respect to in the domain with a small positive number .

Now let us establish uniform convergence of in . For we obtain the following equation from (2.6):

Let . Then where we have chosen small enough such that and .

If attains its positive maximal value in we can choose large enough such that and then at the maximal point we have If attains its negative minimal value in , we have Notice that at . By (2.34) and (2.35), which means that the series , whose sum has the form , is majorized by a geometrical progression and, therefore, is uniformly convergent. The fact that and its derivatives up to the second-order are bounded implies that the first derivatives of are uniformly convergent as .

It follows from (2.6) that are also uniformly convergent as .

Now, we can take ; then by the above discussion, we have the following theorem.

Theorem 2.5. Suppose that and satisfies (2.9) and is suitable smooth, then there exists a small positive number such that the Cauchy problem (2.4) has a classical solution.

By the way, it is easy to prove the uniqueness of the solution for the Cauchy problem (2.4), and we omit the details here.

3. Computational Examples

In this section, a numerical simulate is made for the equations by differential method. Numerical computation of these examples shows that the strong solutions for the corresponding Cauchy problem of (1.1)-(1.2) will blow-up in finite time.

Let and , , but , . Then instead of studying the Cauchy Problem (1.1)-(1.2), we can study the following initial boundary problem: If , it is clear that if is a classical solution of (3.1), then is a strong solution of the Cauchy problem (1.1)-(1.2).

To dissect domain , suppose that and , stands for the space step-length in the axis and axis , and stands for the time step-length. Let and define . The differential scheme of the original equation is (to ensure numerical stability, here we apply arithmetic averages in order to avoid “oscillation’’ and “shifting’’ of the numerical solution)

So we get the following explicit formula:

Experiment 1. Suppose , , , which itself does not satisfy (1.1); we get the graphs (see Figures 13) where changes according to the changes of when different functions are given to .

Figure 1 shows that when , at , the numerical solutions become oscillatory; at the bifurcation of solutions occurs; when the solutions will blow-up. Similarly Figure 2 shows that when , at the bifurcation of solutions occurs; when the solutions will blow-up. Figure 3 is the spatiotemporal graphs of solutions when at (initial value) and . When , the solutions will blow-up.

Experiment 2. The initial value is unknown in the general situation; so we use random numbers as the initial value and draw graphs (see Figures 4 and 5) where changes as changes when different functions are given to .

Figures 4 and 5 show that even though the initial value is sufficiently small, the blow-up will appear in finite time for the different functions.

The numerical result shows that there is a locality solution of the equation. When becomes larger, the bifurcation of solutions occurs in finite time and blow-up appears. For this problem, it is essential to have a further research.


The research was supported by the Fujian National Science Foundation of China Grant 2008J0198, 2009J1009. The authors would like to thank Professor Zhao Junning for insightful discussions and Professor Xu Chuanju for helpful suggestions and comments.


  1. F. Antonelli, E. Barucci, and M. E. Mancino, “A comparison result for FBSDE with applications to decisions theory,” Mathematical Methods of Operations Research, vol. 54, no. 3, pp. 407–423, 2001. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  2. G. Citti, A. Pascucci, and S. Polidoro, “Regularity properties of viscosity solutions of a non-Hörmander degenerate equation,” Journal de Mathématiques Pures et Appliquées, vol. 80, no. 9, pp. 901–918, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  3. F. Antonelli and A. Pascucci, “On the viscosity solutions of a stochastic differential utility problem,” Journal of Differential Equations, vol. 186, no. 1, pp. 69–87, 2002. View at: Publisher Site | Google Scholar | MathSciNet
  4. A. I. Vol'pert and S. I. Hudjaev, “Cauchy's prblem for degenerate second order quasi-linear parabolic equations,” Mathematics of the USSR-Sbornik, vol. 7, no. 3, pp. 365–387, 1969. View at: Publisher Site | Google Scholar
  5. M. Escobedo, J. L. Vazquez, and E. Zuaua, “Entropy solutions for diffusion-convection equations with partial diffusivity,” Transactions of the American Mathematical Society, vol. 342, no. 2, pp. 829–842, 1994. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  6. H. S. Zhan, The study of the Cauchy problem of a second order quasilinear degenerate parabolic equation and the parallelism of a Riemannian manifold, Ph.D. thesis, Xiamen University's, 2004.
  7. H. S. Zhan and J. N. Zhao, “Uniqueness and stability of the solution to the Cauchy problem for a degenerate quasilinear parabolic equation in multiple space variables,” Chinese Annals of Mathematics. Series A, vol. 26, no. 4, pp. 527–536, 2005. View at: Google Scholar | MathSciNet
  8. O. A. Oleinik and V. N. Samokhin, Mathematical Models in Boundary Layer Theory, vol. 15 of Applied Mathematics and Mathematical Computation, Chapman & Hall/CRC, Boca Raton, Fla, USA, 1999. View at: MathSciNet

Copyright © 2009 Zongqi Liang and Huashui Zhan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

438 Views | 328 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.