Abstract and Applied Analysis
Volume 2014 (2014), Article ID 463409, 7 pages
http://dx.doi.org/10.1155/2014/463409

Research Article

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

Carlo Bianca,1,2 Luca Guerrini,³ and Annie Lemarchand1,2

¹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
²CNRS, UMR 7600 LPTMC, Paris, France
³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₁)The probability density function satisfies, for all , the following identity:

 which models the conservation of particles.(A₂) The probability density function is an even function with respect to ; then, in particular,

(A₃) 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

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

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

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

Acknowledgment

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

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.

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.

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.

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.

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.

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.

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.

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.

Y. Kuang, Delay Differential Equations with Applications in Popula tion Dynamics, Mathematics in Science and Engineering, Academic Press, Boston, Mass, USA, 1993.

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.

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.

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.

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.

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.

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.

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.

J. Wei, “Bifurcation analysis in a scalar delay differential equation,” Nonlinearity, vol. 20, no. 11, pp. 2483–2498, 2007.

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.

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.

G. P. Morriss and C. P. Dettmann, “Thermostats: analysis and application,” Chaos, vol. 8, no. 2, pp. 321–336, 1998.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

C. Cercignani, The Boltzmann Equation and Its Applications, vol. 67 of Applied Mathematical Sciences, Springer, New York, NY. USA, 1988.

G. C. Papanicolaou, “Asymptotic analysis of transport processes,” Bulletin of the American Mathematical Society, vol. 81, pp. 330–392, 1975.