The stochastic heat inclusion with fractional time driven by time-space Brownian white noise

Basic Information

• Paper ID: 2506.08275
• Title: The stochastic heat inclusion with fractional time driven by time-space Brownian white noise
• Authors: Olfa Draouil (University of Tunis), Rahma Yasmina Moulay Hachemi (University of Saïda, Algeria), Bernt Øksendal (University of Oslo, Norway)
• Classification: math.PR (Probability Theory)
• Publication Date: October 9, 2025
• Paper Link: https://arxiv.org/abs/2506.08275v2

Abstract

This paper investigates the fractional-time stochastic heat inclusion problem driven by additive time-space Brownian white noise. The fractional time derivative is interpreted as the Caputo derivative of order α ∈ (0,2). The main results include: (a) if a solution exists, it is a fixed point of a specific set-valued mapping; (b) conversely, any fixed point of this mapping is a solution to the heat inclusion. A solution Y(t,x) is called a mild solution if and only if EY²(t,x) < ∞ for all t,x. The study demonstrates that when α=1 and d=1, or α≥1 and d∈{1,2}, the solution is mild; whereas when α<1, the solution is not mild for any spatial dimension d.

Research Background and Motivation

Problem Background

1. Physical Significance: Fractional-order stochastic heat inclusions can model substance transport processes in porous media affected by noise and uncertainty
2. Mathematical Challenges: Combining fractional derivatives, set-valued mappings, and Brownian noise to capture complex physical phenomena
3. Application Value: The theoretical framework directly addresses real-world environmental challenges

Research Motivation

• Limitations of Classical Cases: When α=1, normal diffusion is recovered, but reality exhibits superdiffusion (α>1) and subdiffusion (α<1) phenomena
• Uncertainty Modeling: The drift term G(Y(t,x)) adopts a set-valued form to simulate uncertainty in sinks/sources
• Noise Effects: The time-space white noise W(t,x) represents disturbances in stochastic or noisy media

Core Contributions

1. Existence Theory: Proves the existence of solutions to the fractional-order stochastic heat inclusion by establishing equivalence between solutions and fixed points of set-valued mappings
2. Mild Solution Criteria: Completely characterizes necessary and sufficient conditions for the existence of mild solutions, extending results from the classical α=1 case
3. Mathematical Framework: Establishes a systematic approach for handling fractional-order stochastic partial differential inclusions
4. Application Examples: Provides concrete applications to nitrate transport in sandy porous aquifers, including numerical computations

Methodology Details

Problem Formulation

The fractional-order stochastic heat inclusion problem is formulated as:

∂^α/∂t^α Y(t,x) ∈ λΔY(t,x) + G(Y(t,x)) + σW(t,x)

where:

• Y(t,x) is an unknown S′-valued process
• ∂^α/∂t^α is the Caputo derivative of order α
• Δ is the spatial Laplacian operator
• G is a set-valued mapping
• W(t,x) is time-space white noise

Boundary conditions:

• Y(0,x) = δ(x) (Dirac measure)
• lim_{x→±∞} Y(t,x) = 0

Core Mathematical Tools

1. Mittag-Leffler Functions

• One-parameter: E_α(z) = Σ_^∞ z^k/Γ(αk+1)
• Two-parameter: E_{α,β}(z) = Σ_^∞ z^k/Γ(αk+β)

2. Caputo Fractional Derivative

D^α f(x) = 1/Γ(n-α) ∫_0^x f^(n)(u)du/(x-u)^{α+1-n}

where n-1 < α < n.

3. Set-Valued Function Integration

Using the Aumann definition:

∫_a^b G(t)dt = {∫_a^b g(t)dt : g∈G}

Main Theoretical Results

Theorem 3.1 (Equivalence Theorem)

A solution Y∈Y is a solution to the fractional heat inclusion if and only if Y is a fixed point of the integral inclusion:

Y(t,x) ∈ I₁(t,x) + I₂(t,x,Y(·)) + I₃(t,x)

where the terms correspond respectively to:

• I₁: Initial condition term
• I₂: Integral of the drift term
• I₃: Integral of the noise term

Theorem 4.2 (Mild Solution Criteria)

• Mild Cases: α=1 and d=1, or α≥1 and d∈{1,2}
• Non-mild Cases: α<1 for any spatial dimension d

Experimental Setup

Numerical Example Parameters

• Fractional order: α = 0.8
• Diffusion coefficient: λ = 0.1
• Noise intensity: σ = 0.5
• Nonlinear term: G(Y) = 0.1Y, 0.5Y
• Spatial domain: x ∈ -50, 50 meters
• Time range: t ∈ 0, 30 days

Physical Interpretation

The solution Y(t,x) represents the logarithm of nitrate concentration, with negative values corresponding to concentrations between 0-1 units. The threshold Y(t,x)=0 corresponds to a concentration of 1 unit, calibrated to the WHO safety limit of 50 mg/L.

Experimental Results

Numerical Simulation Findings

1. Spatial Propagation Pattern: Displays non-uniform diffusion of nitrate from a point source
2. Temporal Evolution: Captures memory effects characteristic of subdiffusion (α<1)
3. Risk Zone Identification: Regions where Y(t,x)>0 correspond to high-risk contamination areas

Consistency Between Theory and Numerics

Although Theorem 4.2 indicates the absence of mild solutions when α<1, numerical results demonstrate stable approximations due to:

• Finite domain truncation effects
• Bounded approximations of singular kernels
• Convergence in the weak solution sense

Related Work

Main Research Lines

1. Fractional-Order Stochastic PDEs: Hu (2019), Kochubel et al. (2021)
2. Set-Valued Differential Inclusions: Aumann integral theory
3. Time-Space White Noise Theory: Dalang & Walsh (1992)
4. Mittag-Leffler Function Applications: Kilbas et al. (2006)

Innovations of This Work

• First systematic study of fractional-order stochastic heat inclusions
• Establishes complete mild solution theory
• Provides applications to practical environmental problems

Conclusions and Discussion

Main Conclusions

1. Existence: Proves solution existence through fixed-point theory
2. Mildness Characterization: Completely determines conditions for mild solution existence
3. Application Value: Successfully models environmental pollution transport problems

Limitations

1. Fixed-Point Uniqueness: Uniqueness for Volterra-type inclusions remains unresolved
2. Numerical Methods: Lacks rigorous numerical analysis theory
3. High-Dimensional Cases: The mild solution problem for d≥3 remains open

Future Directions

1. Volterra Inclusion Theory: Develop new fixed-point theorems
2. Numerical Algorithms: Design efficient solvers for fractional-order stochastic PDEs
3. Application Extensions: Generalize to more complex environmental and biological systems

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Complete mathematical derivations with clear proof strategies
2. Problem Importance: Addresses gaps in fractional-order stochastic PDE theory
3. Practical Value: Provides application examples in environmental science
4. Technical Innovation: Skillfully combines Laplace-Fourier transforms for handling fractional derivatives

Weaknesses

1. Theoretical Completeness: Existence and uniqueness theory for Volterra inclusions needs refinement
2. Numerical Verification: Discrepancies between theoretical and numerical results require deeper analysis
3. Application Limitations: Insufficient parameter calibration and validation for practical use

Impact and Significance

1. Academic Contribution: Provides new tools for fractional-order stochastic analysis
2. Application Prospects: Broad applications in environmental modeling, biological diffusion, and related fields
3. Methodological Value: The approach of using set-valued mappings for uncertainty handling has generalization potential

Applicable Scenarios

• Contaminant transport in porous media
• Anomalous diffusion in biological systems
• Fractional Brownian motion modeling in financial markets
• Memory effect research in materials science

References

The paper primarily references Aumann's set-valued integral theory, Mittag-Leffler function theory, Caputo fractional derivative theory, and work by Hu et al. on stochastic heat equations. The paper cites 25 relevant references covering fractional calculus, stochastic analysis, and partial differential equations.