Existence of Solutions of the third term of the Connaughton-Newell Model with a source term

1. Problem Statement

The paper addresses the mathematical well-posedness of the Connaughton-Newell equation, a model derived as an approximation of three-wave kinetic equations using a non-linear coagulation-fragmentation framework. This model describes the energy distribution among weakly interacting waves in wave turbulence.

Specifically, the author focuses on the third operator ($S_3$) of the Connaughton-Newell equation, which represents specific energy transfer mechanisms (associated with back-scattering). The paper investigates the existence of solutions for the evolution equation:
$$\frac{\partial N_\omega}{\partial t} = S_3(N_\omega) + g_\omega(t)$$
where:

• $N_\omega(t)$ is the wave frequency spectrum (wave action density).
• $g_\omega(t)$ is a non-trivial source term.
• The interaction kernel $K_3(\omega, \mu)$ is assumed to be constant (normalized to 1).

The core mathematical challenge is to prove that, under specific physical constraints (non-negativity and integrability), a solution exists for this integro-differential equation over a time interval $[0, T]$.

2. Methodology

The proof employs a combination of integral equation analysis, Riccati equation theory, and fixed-point arguments via auxiliary functions and generating functions. The methodology proceeds in the following steps:

A. Formulation and Assumptions

The author simplifies the operator $S_3$ with a constant kernel to derive the specific equation:
$$\frac{\partial N_\omega}{\partial t} = N_\omega \int_0^\infty N_\mu d\mu + \int_0^\infty N_{\omega+\mu}N_\mu d\mu + g_\omega(t)$$
Four key assumptions are imposed to ensure physical realism and mathematical tractability:

1. (A1) The source term $g_\omega$ is locally integrable in time.
2. (A2) Both the initial condition $N_\omega(0)$ and the source $g_\omega(t)$ are non-negative.
3. (A3) The initial total mass (integral of $N_\omega$) is finite.
4. (A4) The total source term is locally integrable in time.

B. Reduction to a Riccati Equation

The author integrates the partial differential equation (PDE) with respect to $\omega$ to define the total mass $N(t) = \int_0^\infty N_\omega d\omega$. This reduces the system to a scalar Riccati equation:
$$N'(t) = 2N(t)^2 + g(t)$$
where $g(t) = \int_0^\infty g_\omega(t) d\omega$. The non-negativity of $g(t)$ and $N(0)$ ensures that $N(t)$ remains non-negative.

C. Integral Representation and Auxiliary Bounds

Using the Riccati solution $N(t)$, the author rewrites the original PDE as an integral equation:
$$N_\omega(t) = e^{\psi(t)} N_\omega(0) + e^{\psi(t)} \int_0^t e^{-\psi(\tau)} \left( \int_0^\infty N_{\omega+\mu}N_\mu d\mu + g_\omega(\tau) \right) d\tau$$
where $\psi(t) = \int_0^t N(s) ds$.

To prove convergence, the author constructs auxiliary functions:

• $p_\omega(t)$: An upper bound for the solution derived from the initial data and source term.
• $d_\omega(t)$: The supremum of the solution up to time $t$.
• $e_\omega(t)$: A sequence of functions defined recursively to bound the non-linear integral term.

D. Generating Functions and Uniform Convergence

The author introduces generating functions (Laplace-like transforms) $E(t, z)$ and $P(t, z)$ corresponding to the auxiliary functions. This transforms the recursive integral inequality into an algebraic quadratic equation:
$$E(t, z) = P(t, z) + m(t)E(t, z)^2$$
By analyzing the discriminant of this quadratic equation, the author demonstrates that for a sufficiently small time interval $[a, b]$, the series converges, implying that $\int_0^\infty N_\omega d\omega$ is finite and converges uniformly.

E. Extension to Global Existence

Finally, the author uses a continuation argument (Lemma 5). By defining a set $B$ of times where the integral of the solution equals the Riccati solution $N(t)$, and showing that the supremum of this set must be the entire interval $[0, T]$ (via contradiction), the local existence is extended to global existence on $[0, T]$.

3. Key Contributions

1. Existence Proof for $S_3$ Term: The paper provides the first rigorous proof of the existence of solutions specifically for the third operator of the Connaughton-Newell model in the presence of a source term.
2. Handling Non-Linear Source Terms: Unlike previous works that may have focused on pure coagulation or fragmentation, this work integrates a non-trivial source term $g_\omega(t)$, which is crucial for modeling driven turbulent systems.
3. Analytical Technique: The paper successfully adapts the method of generating functions and Riccati reduction to handle the specific non-linear structure of the back-scattering term ($S_3$), which differs from the standard Smoluchowski coagulation equation.
4. Physical Constraints: The proof strictly adheres to physical constraints (non-negativity and mass conservation/integrability), ensuring the mathematical solution is physically relevant.

4. Main Results

• Theorem 1: Under assumptions (A1)–(A4), the equation $\frac{\partial N_\omega}{\partial t} = S_3(N_\omega) + g_\omega(t)$ admits at least one solution on the interval $[0, T]$.
• Consistency: The proof establishes that the total mass of the solution, $M(t) = \int_0^\infty N_\omega d\omega$, is identical to the solution $N(t)$ of the associated Riccati equation $N' = 2N^2 + g(t)$.
• Uniform Convergence: The integral of the solution spectrum converges uniformly on compact time intervals, validating the interchange of integration and differentiation used in the derivation.

5. Significance

This work is significant for the field of wave turbulence theory and kinetic equations:

• Bridging Physics and Analysis: It bridges the gap between the physical modeling of wave energy transfer (Connaughton-Newell) and rigorous mathematical analysis, confirming that the model is well-posed under realistic conditions.
• Foundation for Numerical Schemes: Proving existence is a prerequisite for developing stable numerical schemes. This result supports the validity of numerical methods (like those proposed by Das & Tran) used to simulate these complex wave interactions.
• Generalization: By including a source term, the model becomes applicable to a wider range of physical scenarios, such as driven plasma fusion devices or solar wind turbulence, where external energy injection is constant or time-dependent.

In summary, the paper successfully establishes the mathematical foundation for a specific, complex component of wave turbulence models, ensuring that the theoretical predictions made by the Connaughton-Newell equation are grounded in rigorous existence theory.