Abstract and Applied Analysis

Volume 2014, Article ID 463409, 7 pages

http://dx.doi.org/10.1155/2014/463409

## Existence of Solutions of a Partial Integrodifferential Equation with Thermostat and Time Delay

^{1}Laboratoire de Physique Théorique de la Matière Condensée, Sorbonne Universités, UPMC Univ Paris 06, UMR 7600, Place Jussieu, Case Courrier 121, 75252 Paris Cedex 05, France^{2}CNRS, UMR 7600 LPTMC, Paris, France^{3}Department of Management, Polytechnic University of Marche, 60121 Ancona, Italy

Received 1 July 2014; Accepted 11 August 2014; Published 14 October 2014

Academic Editor: Xiao-Jun Yang

Copyright © 2014 Carlo Bianca et al. 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.

#### Abstract

This paper deals with the mathematical analysis of a retarded partial integrodifferential equation that belongs to the class of thermostatted kinetic equations with time delay. Specifically, the paper is devoted to the proof of the existence and uniqueness of mild solutions of the related Cauchy problem. The main result is obtained by employing integration along the characteristic curves and successive approximations sequence arguments. Applications and perspective are also discussed within the paper.

#### 1. Introduction

The derivation and analysis of mathematical frameworks have recently gained much attention in the applied sciences and specifically for the modeling of complex phenomena occurring in biology, chemistry, vehicular traffic, crowd/swarm dynamics, economics, and social systems. At the origin of complex dynamics, there are interactions that occur in nonlinear fashion and randomly among the elements (particles, cells, and pedestrians) composing the complex system [1]. Complex systems are also characterized by emergent properties, which are properties of the system as a whole which do not exist at the individual element level.

Different ordinary differential equations and partial differential equations have been derived with the aim of obtaining an accurate description of complex phenomena. Moreover, thermostatted kinetic frameworks have been developed for the mathematical modeling of complex systems in physics and life sciences [2–5]. The goal of the thermostatted kinetic models is the possibility of modeling the interactions among the particles at the microscopic scale. In particular, these frameworks allow modeling the ability of the particles to express strategies.

The present paper deals with further developments of the thermostatted kinetic theory proposed in [6]. Specifically, this paper is devoted to the mathematical analysis of a partial integrodifferential equation with thermostat and time delay that belongs to the class of the thermostatted kinetic theory frameworks. The paper focuses on the proof of the existence and uniqueness of mild solutions of the related Cauchy problem. The main result is obtained by employing integration along the characteristic curves and successive approximations sequence arguments. Applications and perspective are also discussed within the paper.

The time delay is introduced in order to take into account the fact that most of the emerging behaviours occurring in complex systems at a certain time are strictly related to the interactions among the particles of the system at a previous time. In the pertinent literature, mathematical models with time delays have been proposed only in ODE-based models; see, among others, [7–16]. The introduction of the time delay has provoked the onset of fluctuations and Hopf bifurcation; see [17–20].

It is worth stressing that, to the best of our knowledge, this is the first time that time delay is introduced into a thermostatted partial integrodifferential equation (kinetic).

The contents of the present paper are developed through four more sections, which follow this introduction. In detail, Section 2 deals with the retarded partial integrodifferential equation and the related Cauchy problem; Section 3 is concerned with some preliminary results that are needed for the proof of the main result that is outlined in Section 4. Finally, Section 5 is devoted to a critical analysis of the proposed mathematical equation including research perspective and applications.

#### 2. The Retarded Integrodifferential Equation

This paper is devoted to the result about the existence and uniqueness of a solution of the following Cauchy problem: where is the unknown function, , , and is the following integral operator: with , the time delay, and the initial datum.

In the partial integrodifferential equation with time delay
is the conservative interaction operator, which is splitted into the gain (of particles into state ) operator and the loss (of particles into state ) operator :
Bearing all of the above in mind, and under suitable integrability assumptions on , the th-order moment of is defined as follows:
In particular, the zero-order, first-order, and second-order moments represent the* density* (mass),* mean activation* (linear momentum), and* activation energy* (kinetic energy), respectively. The term
is a damping operator that allows the control of the activation energy. This term is based on the Gaussian isokinetic thermostat (the interested reader is referred, among others, to [21–23]).

In what follows, we assume that for (initial function).

#### 3. Preliminary Results

This section is concerned with some preliminary results that are at the basis of the main result of the present paper.

Lemma 1. *The gain operator (4) satisfies, for all functions and , the following identity:
*

*Proof. *It is obtained by straightforward calculations.

*The main result is based on the following assumptions on the probability density function .(A _{1})The probability density function satisfies, for all , the following identity:
which models the conservation of particles.(A_{2}) The probability density function is an even function with respect to ; then, in particular,
(A_{3}) The probability density function satisfies, for all , the following identity:
*

*Lemma 2. If the function satisfies assumptions (9), (10), and (11), then
*

*Proof. *Condition (9) implies
Bearing in mind condition (10), we have
Finally, condition (11) implies
Therefore, the proof of the lemma is concluded.

*Bearing the previous lemma in mind, the evolution equation for the 1st-order moment of is stated in the following results.*

*Theorem 3. Assume that assumptions (9), (10), and (11) hold. If there exists a nonnegative solution of the partial integrodifferential equation (3) such that as , then the 1st-order moment is solution of the following delayed first-order nonlinear ordinary differential equation:
*

*Proof. *The integral operator can be written as follows:
Multiplying both sides of by , integrating over , and considering (10), we have
Since
then we have the proof.

*Corollary 4. Assume that assumptions (9), (10), and (11) hold. If there exists a nonnegative solution of the Cauchy problem (1) such that (i),(ii) as ,*

then the 1st-order moment reads as follows: where

*Proof. *The proof is obtained by coupling the delayed differential equation (18) with the initial condition . If or , then we have for all . Otherwise, the unique solution is function (22).

*Theorem 5. Let be an odd number and . Assume that assumptions (9), (10), and (11) hold. Then, the th-order moment of satisfies the following delayed ordinary differential equation:
*

*Proof. *The proof follows by multiplying both sides of (3) by , taking into account assumptions (9), (10), and (11), and performing integration by parts on the thermostat term.

*According to Corollary 4, (3) can be rewritten as follows:
and, after integrating along the characteristics, (25) reads as follows:
where
being
Bearing all of the above in mind, the integral form of (26) is
where .*

*4. Existence of Mild Solutions*

*4. Existence of Mild Solutions**Definition 6. *A function is said to be a mild solution to the Cauchy problem (1) on the time interval if and is solution to the following integral equation:
where

*Lemma 7. Let be the following successive approximations sequence:
where is a nonnegative function such that . Then, admits, as goes to infinity, a nonnegative limit such that .*

*Proof. *Since is a nonnegative function, then , . Moreover,
and .

We now prove by induction that the sequence is monotone, and specifically
and .

Assume as the induction hypothesis that, for some , we have and . Then,
Taking into account property (8), equation (36) thus reads as follows:
and by using the induction hypothesis we conclude that the sequence is monotone. Moreover,
Taking into account property (12), we have
and by using the induction hypothesis we have
Bearing all of the above in mind, we conclude that the sequence has a nonnegative limit such that , as . Then, the Levi theorem implies that satisfies the following equation:
whose unique solution is . Therefore, the lemma is completely proved.

*The main result of this paper is the following.*

*Theorem 8. Let be a given nonnegative function such that . Then, there exists a unique nonnegative mild solution
to the Cauchy problem (1).*

*Proof. *Since Lemma 7 states that solves (31), in order to prove that is a mild solution of (1), it is enough to show that . In order to prove that , we consider the following successive approximations sequence:
or the following equivalent form of (43):
The assumption on implies that the zero-order moment . Assume now as induction hypothesis that , for some . Integrating both sides of (44) over with respect to , and using (28), we obtain
Taking into account property (12) and by using the induction hypothesis, the right-hand side of (45) thus reads as follows:
Therefore, for all , we have .

Multiplying both sides of (44) by , and integrating over with respect to , we have
Taking into account (12), (13) and repeating the computations developed in [6], it is easy to prove by induction on that
Multiplying both sides of (44) by and integrating over with respect to , we have
and according to [6] it is easy to prove by induction on that
Let . Then,
Since and by construction , we have
Therefore, in , and since is bounded, then .

We now prove the uniqueness of the solution. Let be any solution to (31). The positivity of the operators and implies that, for all , , and then . Since , we thus have .

*Corollary 9. Let be a given nonnegative function such that (i);(ii).*

Then, the mild solution of the Cauchy problem (1) belongs to where

*Proof. *It is sufficient to remember (24) and note that assumption (ii) implies that
which allows us to conclude that .

*5. Applications and Research Perspectives*

*5. Applications and Research Perspectives**The main goal of the present paper refers to the proof of the global existence of mild solutions of a thermostatted partial integrodifferential equation with time delay. As already mentioned in Section 2, this equation can be proposed for the modeling of complex systems in nature and society where only the interactions at the microscopic scales are affected by time delay. Specifically, the partial integrodifferential equation with time delay
can be proposed as a general mathematical framework for the modeling of complex systems composed by a large number of interacting particles and subjected to the external force field . The overall state of the system is described by the distribution function (statistical description). Particles are able to express a specific function; this ability of the particles is modeled by the variable . Moreover, is the encounter rate between particles with states (or ) and . Finally, is the probability density that a particle with state ends up into the state after the interaction with the particle with state . The action of the external force field is controlled by the thermostat term that, as shown in Section 4, maintains constant the first-order and the second-order moments (number density and activation energy of the system).*

*From the applications viewpoint, we consider a simple model for the evolution of malignancy in tumor cells. Specifically, we assume that the variable models the magnitude of the malignancy of tumor cells and is a delta Dirac function (deterministic output of a pair interaction) depending on the microscopic state of the interacting particles:
and finally we assume that the malignancy of tumor cells increases when cells interact with each other with rate . Accordingly, we have
where is a positive parameter. Bearing all of the above in mind, the integral operator reads as follows:
and then delayed equation (55) now reads as follows:
*

*However, thermostatted equation (55) does not include the role of the space and velocity variables; then, applications refer to the modeling of complex phenomena that are homogeneous in space and velocity. The mathematical analysis performed in the present paper has to be thus generalized for taking also into account the dynamics described by these variables. Moreover, (55) refers to complex systems where the mutual interactions do not produce modification in the number density (conservative interactions).*

*Research perspectives include the possibility of performing an asymptotic analysis by parabolic (low-field) and hyperbolic (high-field) scalings (see [24–34]) with the aim of obtaining the dynamics of the system at the macroscopic scale. This is a work in progress and results will be presented in due course.*

*Conflict of Interests*

*Conflict of Interests**The authors declare that there is no conflict of interests regarding the publication of this paper.*

*Acknowledgment*

*Acknowledgment**Carlo Bianca and Annie Lemarchand were partially supported by L’Agence Nationale de la Recherche (ANR T-KiNeT Project).*

*References*

*References*

- Y. Bar-Yam,
*Dynamics of Complex Systems*, Studies in Nonlinearity, Westview Press, 2003. - C. Bianca, “Thermostatted kinetic equations as models for complex systems in physics and life sciences,”
*Physics of Life Reviews*, vol. 9, no. 4, pp. 359–399, 2012. View at Publisher · View at Google Scholar · View at Scopus - C. Bianca, “How do mutative events modify moments evolution in thermostatted kinetic models?”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 19, no. 7, pp. 2155–2159, 2014. View at Publisher · View at Google Scholar · View at MathSciNet - C. Bianca, C. Dogbe, and L. Guerrini, “A thermostatted kinetic framework with particle refuge for the modeling of tumors hiding,”
*Applied Mathematics & Information*, vol. 8, no. 2, pp. 469–473, 2014. View at Google Scholar - C. Bianca, M. Ferrara, and L. Guerrini, “High-order moments conservation in thermostatted kinetic models,”
*Journal of Global Optimization*, vol. 58, no. 2, pp. 389–404, 2014. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - C. Bianca, “Modeling complex systems by functional subsystems representation and thermostatted-KTAP methods,”
*Applied Mathematics and Information Sciences*, vol. 6, no. 3, pp. 495–499, 2012. View at Google Scholar · View at Scopus - P. J. Cunningham and W. J. Wangersky,
*Time Lag in Population Models*, Yale, 1958. - T. Erneux,
*Applied Delay Differential Equations*, Springer, New York, NY, USA, 2009. View at MathSciNet - C. T. H. Baker, G. A. Bocharov, and C. A. H. Paul, “Mathematical modelling of the interleukin-2 T-cell system: a comparative study of approaches based on ordinary and delay differential equations,”
*Journal of Theoretical Medicine*, vol. 2, pp. 117–128, 1997. View at Google Scholar - J. Bélair, M. C. Mackey, and J. M. Mahay, “Age-structured and two-delay models for erythropoiesis,”
*Mathematical Biosciences*, vol. 128, no. 1-2, pp. 317–346, 1995. View at Publisher · View at Google Scholar · View at Scopus - Y. Kuang,
*Delay Differential Equations with Applications in Popula tion Dynamics*, Mathematics in Science and Engineering, Academic Press, Boston, Mass, USA, 1993. View at MathSciNet - H. M. Byrne, “The effect of time delays on the dynamics of avascular tumor growth,”
*Mathematical Biosciences*, vol. 144, no. 2, pp. 83–117, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - M. J. Piotrowska, “Hopf bifurcation in a solid avascular tumour growth model with two discrete delays,”
*Mathematical and Computer Modelling*, vol. 47, no. 5-6, pp. 597–603, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - C. Bianca, M. Ferrara, and L. Guerrini, “Hopf bifurcations in a delayed-energy-based model of capital accumulation,”
*Applied Mathematics and Information Sciences*, vol. 7, no. 1, pp. 139–143, 2013. View at Google Scholar · View at Scopus - C. Bianca and L. Guerrini, “Existence of limit cycles in the Solow model with delayed-logistic population growth,”
*The Scientific World Journal*, vol. 2014, Article ID 207806, 8 pages, 2014. View at Publisher · View at Google Scholar - C. Bianca, M. Ferrara, and L. Guerrini, “Qualitative analysis of a retarded mathematical framework with applications to living systems,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 736058, 7 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Ma, “Global Hopf bifurcation in the Leslie-Gower predator-prey model with two delays,”
*Nonlinear Analysis: Real World Applications*, vol. 13, no. 1, pp. 370–375, 2012. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - Y. Song, J. Wei, and M. Han, “Local and global Hopf bifurcation in a delayed hematopoiesis model,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 14, no. 11, pp. 3909–3919, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. Wei, “Bifurcation analysis in a scalar delay differential equation,”
*Nonlinearity*, vol. 20, no. 11, pp. 2483–2498, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - T. Chen, X. Li, and J. He, “Complex dynamics of credit risk contagion with time-delay and correlated noises,”
*Abstract and Applied Analysis*, vol. 2014, Article ID 456764, 10 pages, 2014. View at Publisher · View at Google Scholar · View at MathSciNet - K. F. Gauss, “Uber ein Neues Allgemeines Grundgesatz der Mechanik (On a New Fundamental Law of Mechanics),”
*Journal für die Reine und Angewandte Mathematik*, vol. 4, pp. 232–235, 1829. View at Google Scholar - G. P. Morriss and C. P. Dettmann, “Thermostats: analysis and application,”
*Chaos*, vol. 8, no. 2, pp. 321–336, 1998. View at Publisher · View at Google Scholar · View at Scopus - D. Ruelle, “Smooth dynamics and new theoretical ideas in nonequilibrium statistical mechanics,”
*Journal of Statistical Physics*, vol. 95, no. 1-2, pp. 393–468, 1999. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - C. Bardos, F. Golse, and D. Levermore, “Fluid dynamic limits of kinetic equations. I. Formal derivations,”
*Journal of Statistical Physics*, vol. 63, no. 1-2, pp. 323–344, 1991. View at Publisher · View at Google Scholar - P. Lions and N. Masmoudi, “From the Boltzmann equations to the equations of incompressible fluid mechanics—I,”
*Archive for Rational Mechanics and Analysis*, vol. 158, no. 3, pp. 173–193, 2001. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - C. Dogbe, “Fluid dynamic limits for gas mixture. I. Formal derivations,”
*Mathematical Models and Methods in Applied Sciences*, vol. 18, no. 9, pp. 1633–1672, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - T. Hillen and H. G. Othmer, “The diffusion limit of transport equations derived from velocity-jump processes,”
*SIAM Journal on Applied Mathematics*, vol. 61, no. 3, pp. 751–775, 2000. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J. Soler, L. L. Bonilla, O. Sánchez, and T. Goudon, “Low-field limit for a nonlinear discrete drift-diffusion model arising in semiconductor superlattices theory,”
*SIAM Journal on Applied Mathematics*, vol. 64, no. 5, pp. 1526–1549, 2004. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - F. Poupaud, “Runaway phenomena and fluid approximation under high fields in semiconductor kinetic theory,”
*Zeitschrift für Angewandte Mathematik und Mechanik*, vol. 72, no. 8, pp. 359–372, 1992. View at Publisher · View at Google Scholar · View at MathSciNet - P. Degond and B. Wennberg, “Mass and energy balance laws derived from high-field limits of thermostatted Boltzmann equations,”
*Communications in Mathematical Sciences*, vol. 5, no. 2, pp. 355–382, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - H. Grad, “Asymptotic theory of the Boltzmann equation,” in
*Rarefied Gas Dynamics (Proceedings of the 3rd International Symposium)*, vol. 1, pp. 26–59, Palais de l'UNESCO, Paris, France, 1962. View at Google Scholar - C. Bianca and A. Lemarchand, “Density evolution by the low-field limit of kinetic frameworks with thermostat and mutations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 20, no. 1, pp. 14–23, 2015. View at Publisher · View at Google Scholar · View at MathSciNet - C. Cercignani,
*The Boltzmann Equation and Its Applications*, vol. 67 of*Applied Mathematical Sciences*, Springer, New York, NY. USA, 1988. View at Publisher · View at Google Scholar · View at MathSciNet - G. C. Papanicolaou, “Asymptotic analysis of transport processes,”
*Bulletin of the American Mathematical Society*, vol. 81, pp. 330–392, 1975. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet

*
*