Mass-conserving solutions to coagulation-fragmentation equations with balanced growth
Philippe Lauren\c{c}ot (IMT)

TL;DR
This paper studies a special class of coagulation-fragmentation equations where coagulation and fragmentation are balanced, proving existence and uniqueness of mass-conserving solutions for small initial mass.
Contribution
It introduces a specific class of coefficients ensuring mass conservation and establishes existence and uniqueness results under small initial mass conditions.
Findings
Existence of mass-conserving solutions for small initial mass
Uniqueness of solutions in the considered class
Balance between coagulation and fragmentation coefficients
Abstract
A specific class of coagulation and fragmentation coefficients is considered for which the strength of the coagulation is balanced by that of the multiple fragmentation. Existence and uniqueness of mass-conserving solutions are proved when the initial total mass is sufficiently small.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Mass-conserving solutions to coagulation-fragmentation equations with balanced growth
Philippe Laurençot
Institut de Mathématiques de Toulouse, UMR 5219, Université de Toulouse, CNRS
F–31062 Toulouse Cedex 9, France
Abstract.
A specific class of coagulation and fragmentation coefficients is considered for which the strength of the coagulation is balanced by that of the multiple fragmentation. Existence and uniqueness of mass-conserving solutions are proved when the initial total mass is sufficiently small.
Key words and phrases:
coagulation, fragmentation, mass conservation
1991 Mathematics Subject Classification:
45K05
1. Introduction
The time evolution of the size distribution function of a population of particles growing by successive mergers and spontaneously breaking apart into smaller pieces is described by the coagulation-fragmentation equation
[TABLE]
respectively. In (1.1), denotes the size distribution function of particles of size at time and the sizes of the particles vary as a consequence of coagulation and fragmentation reactions. In (1.1c), denotes the coagulation kernel which is a non-negative and symmetric function defined on and is the rate at which two particles with respective sizes and merge. The first term in (1.1c) accounts for the formation of particles of size resulting from the coalescence of particles with respective sizes and while the second term describes the disappearance of particles of size after coagulation with other particles. In (1.1d), is the overall fragmentation rate of particles with size and the first term in (1.1d) describes the loss of particles of size due to spontaneous breakage. The second term in (1.1d) accounts for the appearance of new particles of size resulting from the break-up of a particle of size and is the daughter distribution function providing the distribution of particles with size resulting from the breakage of a particle with size . We only consider here the case where there is no loss of matter during splitting; that is,
[TABLE]
In fact, there is also no loss of matter during coalescence events, so that conservation of matter is expected throughout time evolution. From a mathematical point of view, this means that the first moment of the distribution function at time , which somehow corresponds to the total mass of the system of particles, remains unchanged as time goes by; that is,
[TABLE]
It is however by now well-known that the conservation of matter may fail to be true and there are basically two mechanisms responsible for the infringement of (1.3). On the one hand, for coagulation kernels increasing rapidly with size, a runaway growth might occur, thereby producing particles with infinite size. Since the distribution function does not include such particles, matter escapes from the system of particles with finite size to form giant clusters and the first moment of decays after some time. This phenomenon is known as gelation and takes place when increases faster than as and the fragmentation is weak enough [6, 13, 12, 20, 21, 28, 29]. On the other hand, if there is no coagulation and the overall fragmentation rate is unbounded as , small particles break apart at a high rate and thus reduce to particles with size zero or dust which are also not accounted for in the distribution function . Another loss of matter occurs in that case and is known as shattering [2, 6, 15, 30]. To be more specific, when the coagulation and fragmentation coefficients are given by
[TABLE]
with , , and , the following results are available: on the one hand, shattering is observed when in (1.4b) and there is no coagulation [2, 15, 30]. On the other hand, the gelation phenomenon takes place when in (1.4a), in (1.4b), and the total mass of the initial condition is sufficiently large compared to the parameter in (1.4b) [12, 13, 21, 39]. Mass-conserving solutions exist when and, either , or and [1, 9, 12, 14, 24, 37, 41], which is in accordance with the numerical simulations and formal computations performed in [35, 39]. Let us point out that the above mentioned results do not depend on the daughter distribution function when only a finite number of fragments are formed during breakup; that is, when for all . This corresponds to the choice in (1.4b). When , the fact that an infinite number of fragments are produced after breakup alters the situation [26, 39].
According to the previous discussion, the value appears to be a borderline case with respect to the occurrence of the gelation phenomenon. The purpose of this paper is to contribute to a better understanding of this particular case for which nothing much seems to be known, except for the particular choice of parameters
[TABLE]
In that case, explicit computations performed in [35, 39] reveal that gelation is likely to take place when is suitably small and and a proof that it indeed occurs when may be found in [6]. On the opposite, numerical simulations performed in [35] indicate that there are mass-conserving solutions when and the ratio is large enough. Still for , it follows from [10, Remark 5.2] that there is a stationary solution with total mass . Finally, when , the initial value problem (1.1) has at least one mass-conserving solution whatever the values of , , and in (1.3) and (1.4) and the value of [26].
We shall here extend these results to a larger class of coagulation and fragmentation coefficients, showing that mass-conserving solutions exist for a suitable range of the values of the parameters and , provided the ratio is large enough. More precisely, consider
[TABLE]
When
[TABLE]
the assumptions (1.5e) and (1.5f) are clearly satisfied for all , while (1.5h) requires . Let us point out that (1.5h) is implied by (1.5a) when .
Remark 1.1**.**
Clearly for all , so that is non-empty. We also observe that, if , then also belongs to for all when and for all when .
An easy consequence of (1.5f) and Remark 1.1 is that
[TABLE]
Indeed,
[TABLE]
and the right-hand side of the previous inequality is finite since . We then set
[TABLE]
and define the weighted -space and the moment of order of for by
[TABLE]
We also denote the positive cone of by while denotes the space endowed with its weak topology.
With the above notation, the main contribution of this paper is the existence and uniqueness of mass-conserving weak solutions to (1.1) on when .
Theorem 1.2**.**
Let , , and be coagulation and fragmentation coefficients satisfying (1.5) and fix
[TABLE]
Consider a non-negative initial condition such that
[TABLE]
and assume that
[TABLE]
- (a)
There is a mass-conserving weak solution to (1.1) on in the sense of Definition 2.1 below, which additionally belongs to for all .
- (b)
If for some , then for all ,
- (c)
If belongs to , then there is a unique mass-conserving weak solution to (1.1) on such that for all .
The proof of Theorem 1.2 (a)-(b) relies on a weak compactness method in a weighted -space, an approach pioneered in [37] and further developed in [12, 14, 17, 18, 21, 22, 24]. The starting point of the analysis is to construct a sequence of approximations of (1.1) for which the existence of a solution follows by an application of Banach fixed point theorem. To guarantee the -weak compactness of the sequence of approximating solutions thus obtained, several estimates are to be derived. In particular, the behaviour of the sequence of the approximating solutions for small and large sizes has to be controlled and the balance between the homogeneity of the coagulation kernel and the overall fragmentation rate is the main difficulty to be overcome here. In fact, special care has to be paid to the control for large sizes and this is where the smallness condition on the total mass comes into play. The key idea is here to obtain joint estimates in and in when the total mass of the initial condition does not exceed . Uniform integrability and time equicontinuity are next derived, the former relying on the monotonicity properties of the coagulation kernel [7, 11, 23, 33], and complete the proof of the weak compactness of the sequence of approximating solutions. Passing to the limit as the approximating parameter diverges to infinity is then done along the lines of the above mentioned references. The uniqueness statement in Theorem 1.2 (c) is proved along the lines of [14, Theorem 2.9].
Remark 1.3**.**
Note that the value defined in (1.8) is finite whatever the value of . On the one hand, that mass-conserving solutions cannot exist when exceeds a positive threshold value is supported by the occurrence of gelation for large values of when , , and with , see (1.6), [6]. However, it is likely that : for instance, it is expected that when , , and (), and in that case. On the other hand, the analysis performed in [26] reveals that gelation does not take place when , , and (), so that the outcome of Theorem 1.2 (a) is far from being optimal in that case. However, being valid for all values of in , the proof of Theorem 1.2 does not exploit the specific non-integrability features of when and a more refined analysis seems to be required. We hope to return to that problem in the near future.
2. Weak Solutions: Existence
In this section, we assume that , , and are coagulation and fragmentation coefficients satisfying (1.5) and that the initial condition belongs to and possesses the integrability properties (1.10) and (1.11) with satisfying (1.9). We also fix
[TABLE]
and a positive real number satisfying
[TABLE]
Throughout this section, and are positive constants depending only on , , , , , , , , , and . Dependence upon additional parameters is indicated explicitly.
2.1. Weak Solutions
We begin this section with the definition of a weak solution to (1.1) which is slightly different from the usual one in order to account for the possible non-integrability of the daughter distribution [6]. It involves the spaces of test functions , , which are defined by and
[TABLE]
Definition 2.1**.**
Let , , and be coagulation and fragmentation coefficients satisfying (1.5) and consider satisfying (1.9) and (2.1) and a non-negative function . For , a weak solution to (1.1) on is a non-negative function
[TABLE]
such that, for all and ,
[TABLE]
where
[TABLE]
In addition, is a mass-conserving weak solution to (1.1) on if and for all .
The existence of weak solutions and mass-conserving weak solutions to coagulation-fragmentation equations has been investigated for various classes of coagulation and fragmentation coefficients, see [1, 3, 4, 5, 6, 9, 11, 12, 14, 17, 18, 21, 22, 24, 31, 32, 36, 37, 41] and the references therein. As for the uniqueness issue, we refer to [1, 3, 4, 5, 6, 9, 14, 16, 19, 25, 34, 38] and the references therein.
2.2. Approximation
We fix an initial condition satisfying (1.10) and (1.11) and recall that . Let be an integer and define
[TABLE]
and
[TABLE]
Denoting the coagulation and fragmentation operators with and instead of and by and , respectively, it follows from (1.5), (2.6), and (2.7) by a classical Banach fixed point argument that there is a unique non-negative function
[TABLE]
which solves
[TABLE]
We extend to by setting for and . It readily follows from (2.9a) that, for and satisfying ,
[TABLE]
Choosing , , in (2.10) gives
[TABLE]
hence, thanks to the symmetry of ,
[TABLE]
After integration with respect to time, we find
[TABLE]
Recalling that for and , the previous identity becomes
[TABLE]
and we infer from (2.8), (2.11), and the non-negativity of and that
[TABLE]
2.3. Moment Estimates
Let us first draw a couple of consequences of (1.5).
Lemma 2.2**.**
Let , , and be coagulation and fragmentation coefficients satisfying (1.5). Then
[TABLE]
[TABLE]
Proof.
The first assertion (2.13) readily follows from (1.5a). Consider next . On the one hand, it follows from the Cauchy-Schwarz inequality that
[TABLE]
On the other hand, using Young’s inequality gives
[TABLE]
and the proof of (2.14) is complete. ∎
The next result is at the heart of the control of solutions to (1.1) for large sizes.
Lemma 2.3**.**
For ,
[TABLE]
Proof.
Fix and define
[TABLE]
Then for with for . Then and studying the variation of reveals that there is such that in and in . Consequently, reaches its maximum twice, at and , from which the non-positivity of on follows. ∎
After this preparation, we are in a position to state and prove the first estimates for for small and large sizes.
Lemma 2.4**.**
Let and set . There is depending on such that, for ,
[TABLE]
Proof.
Since , there holds
[TABLE]
[TABLE]
Consequently, we infer from (1.5b), (2.10) (with , ), and the non-negativity of and that, for ,
[TABLE]
As the choice of ensures that , it follows from (2.12) and Hölder’s inequality that
[TABLE]
which gives, together with (2.15) and Young’s inequality,
[TABLE]
We next set , , and notice that
[TABLE]
by Lemma 2.2 and Lemma 2.3 and
[TABLE]
Also, owing to (1.5e) and (1.7),
[TABLE]
Collecting the previous estimates, it follows from (1.5b), (2.10) with , , and (2.12) that, for ,
[TABLE]
Combining (2.16) and (2.17) leads us to
[TABLE]
Hence, after integration with respect to time,
[TABLE]
for . We finally use the inequality
[TABLE]
to complete the proof. ∎
We next turn to estimates for moments of order .
Lemma 2.5**.**
Consider and assume additionally that . There is depending on such that
[TABLE]
Proof.
Since , we argue as at the beginning of the proof of Lemma 2.4 to derive from (2.10) that
[TABLE]
The range of and (1.5a) next imply that and we infer from Hölder’s inequality that
[TABLE]
Combining the previous inequalities leads us to
[TABLE]
hence, after integration with respect to time,
[TABLE]
for . It then follows from (2.2), (2.8), Lemma 2.4, and Hölder’s inequality that
[TABLE]
where
[TABLE]
Lemma 2.5 then readily follows. ∎
We end up the study of the evolution of moments with moments of higher order.
Lemma 2.6**.**
Let and assume additionally that . For every , there is depending on and such that
[TABLE]
Proof.
We first recall that, since , there is depending only on such that
[TABLE]
see [6] or [8, p. 216] for instance. Let and and define
[TABLE]
A straightforward consequence of (2.2) and Lemma 2.4 is that
[TABLE]
It follows from (1.5c) and the previous inequality that
[TABLE]
On the one hand, owing to (2.13) and the range of , both and belong to and we infer from (2.12) and Hölder’s inequality that
[TABLE]
and
[TABLE]
On the other hand, since and due to (2.13) and the choice of , it follows from (2.12), (2.13), and Hölder’s inequality
[TABLE]
and
[TABLE]
Collecting the above estimates leads us to
[TABLE]
Now, since
[TABLE]
and by (1.5d), (1.5e), and (1.5f), it follows from (1.5b), (2.10), and (2.20) that
[TABLE]
By (2.12) and Hölder’s inequality,
[TABLE]
so that
[TABLE]
where
[TABLE]
Observe that and , so that the first term on the right-hand side of (2.21) is strictly dominated by the last one.
Now, let and consider . It readily follows from (2.19) that there is large enough such that
[TABLE]
Taking in (2.21) and using the above inequality as well as Young’s inequality, we end up with
[TABLE]
We then infer from the comparison principle that
[TABLE]
for , and we use the obvious bound due to (2.8) to complete the proof. ∎
2.4. Uniform Integrability
The outcome of the previous section guarantees that there is no leak of matter for small and large sizes. The purpose of the next result is to prevent concentration of matter at a finite size by showing the uniform integrability of the sequence in . To this end, we first recall that, since , a refined version of the de la Vallée-Poussin theorem [6, 27] ensures that there is a non-negative and convex function with concave first derivative such that
[TABLE]
As a consequence of (2.22), the convexity of , and the concavity of , there holds
[TABLE]
and
[TABLE]
It is furthermore possible to construct such that
[TABLE]
see [6]. We next fix two additional parameters and which only depend on and and satisfy
[TABLE]
Lemma 2.7**.**
Let . There is depending on and such that
[TABLE]
In the proof of Lemma 2.7, we first exploit the monotonicity of the coagulation kernel, the subadditivity of the weight function , and the specific choice (2.6) of the truncation of to show that the contribution of the coagulation term is non-positive, a property which has been uncovered in [7] and subsequently used in [6, 11, 23, 33]. The analysis of the contribution of the fragmentation term is more delicate and the specific choice of the additional parameters and come into play there.
Proof.
Let . On the one hand, it follows from Fubini’s theorem, the symmetry of , and the subadditivity of that
[TABLE]
Owing to (2.24) (with and ), we further obtain
[TABLE]
We finally infer from the monotonicity of and the non-negativity (2.23) of that
[TABLE]
On the other hand, using again Fubini’s theorem and (2.24) (with and ),
[TABLE]
where
[TABLE]
Since , we infer from (1.5b), (1.5d), and (2.25) that
[TABLE]
At this point, we notice that, due to (2.26) and the non-positivity of ,
[TABLE]
so that and is finite. Consequently, we deduce from (1.5b), (2.23), and the estimate on that
[TABLE]
Moreover, using once more (2.26), we realize that
[TABLE]
and we infer from (2.12), Lemma 2.4 (with ), and Hölder’s and Young’s inequalities that
[TABLE]
and
[TABLE]
Consequently,
[TABLE]
It now follows from (2.9a), (2.27), and (2.28) that, for ,
[TABLE]
hence, after integration with respect to time,
[TABLE]
We then use (2.8), (2.22b), Lemma 2.4 (with )), and the monotonicity of to complete the proof. ∎
2.5. Time Equicontinuity
According to the Dunford-Pettis theorem, the estimates derived in Sections 2.3 and 2.4 along with the superlinearity (2.22a) of at infinity imply the weak compactness in of the sequence for all . To show the convergence of the sequence , we are left with the compactness with respect to the time variable which we establish now.
Lemma 2.8**.**
For , there exists depending on such that
[TABLE]
with , .
Proof.
Let and . On the one hand, we infer from Fubini’s theorem that
[TABLE]
Using the subadditivity of and (1.5c), we further obtain
[TABLE]
Since by (1.5a) and (1.9), we deduce from (2.2), (2.12), Hölder’s inequality, and, either Lemma 2.4 (with ) if , or Lemma 2.5 (with ) if , that
[TABLE]
On the other hand, using (1.5b), (1.5d), (1.5f), Fubini’s theorem, and (see (2.1) and Remark 1.1), we find
[TABLE]
Hence, thanks to (1.5a), (2.1), and (2.12),
[TABLE]
We now combine (2.9a), (2.29), and (2.30) to conclude that, for ,
[TABLE]
Therefore, for , we deduce from (2.2), Lemma 2.4 (with ), and Fubini’s theorem that
[TABLE]
. It then follows from (2.1), (2.12), and (2.31) that, for and ,
[TABLE]
The previous inequality being valid for any , we may choose when and otherwise and thereby complete the proof of Lemma 2.8. ∎
Corollary 2.9**.**
For , there exists depending on such that
[TABLE]
where
[TABLE]
Proof.
Consider and . We infer from (2.1), (2.2), Lemma 2.4 (with ), and Lemma 2.8 that
[TABLE]
and Corollary 2.9 readily follows since the previous inequality is valid for all . ∎
2.6. Convergence
We are now in a position to complete the proof of the first statements in Theorem 1.2.
Proof of Theorem 1.2 (a): Existence.
Let . Since and are superlinear at infinity, it follows from Lemma 2.4 (with ), Lemma 2.7, and the Dunford-Pettis theorem that there is a weakly compact subset of such that for all and . In addition, the function defined in Corollary 2.9 satisfies as , so that the sequence is equicontinuous in at any by Lemma 2.8 and Corollary 2.9. According to a variant of the Arzelà-Ascoli theorem [40, Theorem A.3.1], these properties imply the compactness of the sequence in and in , recalling that denotes the space endowed with its weak topology. A diagonal process then allows us to construct a subsequence of (not relabeled) and a non-negative function such that
[TABLE]
A first consequence of, either (2.32) if , or (2.32), Lemma 2.5 (with ), and Fatou’s lemma if , is that
[TABLE]
Consider next . Thanks to (2.2) and Lemma 2.4 (with ),
[TABLE]
On the one hand, we may pass to the limit as in (2.11) with the help of (2.8), (2.32), and (2.34) to conclude that
[TABLE]
On the other hand, we infer from (2.32), (2.35), and Fatou’s lemma that
[TABLE]
We are left with showing that satisfies the weak formulation (2.3) of (1.1) for all . Consider thus and observe that, sice , it follows from (2.10) after integration with respect to time that
[TABLE]
First, since , we infer from (2.1) and (2.34) that
[TABLE]
Therefore,
[TABLE]
Next, since and by (2.1) and Remark 1.1, it follows from (1.5d) and (1.5f) that, for ,
[TABLE]
Then, using (1.5b), (2.35), (2.37), and (2.40), we obtain, for ,
[TABLE]
Now, belongs to by (1.5b) and (2.40) which, together with (2.32), implies that the first term in the right-hand side of the previous inequality converges to zero as . Consequently, for all ,
[TABLE]
Letting in the above inequality gives, since ,
[TABLE]
Finally, since by (1.5a) and (1.9), the by now classical argument designed in [37], and further developed in [6, 12, 24], can be applied to deduce from Lemma 2.5 (with ), (2.32), and (2.33) that
[TABLE]
Collecting (2.39), (2.41), and (2.42), we may take the limit in (2.38) to obtain (2.3), after handling the left-hand side of (2.38) with (2.8) and (2.32). ∎
Proof of Theorem 1.2 (b): Local boundedness of moments.
Consider such that and let be the mass-conserving weak solution to (1.1) on given by (2.32). For , it readily follows from (2.32), Lemma 2.5 (if ) or Lemma 2.6 (if ), and Fatou’s lemma that as claimed. ∎
3. Weak Solutions: Uniqueness
The final step of the proof of Theorem 1.2 is the uniqueness of mass-conserving weak solutions to (1.1) on having a finite moment of sufficiently high order and the proof given below follows quite closely the approach developed in [14, 25, 34].
Proof of Theorem 1.2 (c): Uniqueness.
Let be an initial condition in and two mass-conserving weak solutions and to (1.1) on which satisfy
[TABLE]
Observe that the existence of at least such a solution is guaranteed by Theorem 1.2 (a)-(b) as . According to Definition 2.1 and the choice (1.9) of , there also holds
[TABLE]
We set for , , and . Let . We infer from (1.1) that, for ,
[TABLE]
On the one hand, since , it follows from the convexity of and the subadditivity of and that, for ,
[TABLE]
Therefore, by (1.5c),
[TABLE]
Observing that all the exponents in the powers of range in and all the exponents in the powers of range in by (1.5a), it further follows from several applications of Young’s inequality that
[TABLE]
On the other hand, we infer from (1.5d), (1.5e), (1.5f), and (1.5h) that, for ,
[TABLE]
Inserting the previous estimates in (3.3) and using (1.5b), (3.1), (3.2), and the inequality , , lead us to
[TABLE]
Integrating with respect to time gives
[TABLE]
for . Hence, for all , and the claimed uniqueness follows as is arbitrary in . ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] J. M. Ball and J. Carr , The discrete coagulation-fragmentation equations: existence, uniqueness, and density conservation , J. Statist. Phys., 61 (1990), pp. 203–234.
- 2[2] J. Banasiak , Shattering and non-uniqueness in fragmentation models—an analytic approach , Phys. D, 222 (2006), pp. 63–72.
- 3[3] J. Banasiak and W. Lamb , Global strict solutions to continuous coagulation-fragmentation equations with strong fragmentation , Proc. Roy. Soc. Edinburgh Sect. A, 141 (2011), pp. 465–480.
- 4[4] , Analytic fragmentation semigroups and continuous coagulation-fragmentation equations with unbounded rates , J. Math. Anal. Appl., 391 (2012), pp. 312–322.
- 5[5] J. Banasiak, W. Lamb, and M. Langer , Strong fragmentation and coagulation with power-law rates , J. Engrg. Math., 82 (2013), pp. 199–215.
- 6[6] J. Banasiak, W. Lamb, and Ph. Laurençot , Analytic methods for coagulation-fragmentation models . Book in preparation.
- 7[7] A. V. Burobin , Existence and uniqueness of the solution of the Cauchy problem for a spatially nonhomogeneous coagulation equation , Differentsial ′ nye Uravneniya, 19 (1983), pp. 1568–1579.
- 8[8] J. M. C. Clark and V. Katsouros , Stably coalescent stochastic froths , Adv. in Appl. Probab., 31 (1999), pp. 199–219.
