Charles R. Doering,Khachik V. Sargsyan,Leonard M. Sander
Abstract: We consider extinction times for a class of birth-death processes commonly found in applications, where there is a control parameter which determines whether the population quickly becomes extinct, or rather persists for a long time. We give an exact expression for the discrete case and its asymptotic expansion for large values of the population. We have results below the threshold, at the threshold, and above the threshold (where there is a quasi-stationary state and the extinction time is very long.) We show that the Fokker-Planck approximation is valid only quite near the threshold. We compare our analytical results to numerical simulations for the SIS epidemic model, which is in the class that we treat. This is an interesting example of the delicate relationship between discrete and continuum treatments of the same problem.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the behavior of the probability density function (PDF) of a passive scalar that decays in a periodic domain under the influence of a smooth large - scale velocity field. Specifically, the author focuses on the situation where the normalized passive scalar PDF reaches a self - similar profile within a finite time, namely the so - called "strange or statistical eigenmode regime". The following is a detailed description of this problem:
### Research Background and Problems
1. **Evolution Equation of Passive Scalar**
The evolution of the passive scalar \(\phi(x, y, t)\) is described by the advection - diffusion (AD) equation:
\[
\frac{\partial \phi}{\partial t}+(\vec{u}\cdot\nabla)\phi = \kappa\nabla^{2}\phi
\]
where \(\kappa\) is the molecular diffusion coefficient and \(\vec{u}(x, y, t)\) is the advection velocity field.
2. **Research Objectives**
- To study how the PDF of the normalized passive scalar evolves in a periodic domain, especially in the strange eigenmode.
- To explore the influence of different physical parameters (such as the Péclet number) on the shape of the PDF.
- To verify the PDF characteristics predicted by the Sinai and Yakhot theory, including the power - law distribution and the Gaussian distribution forms.
3. **Strange Eigenmode**
In the strange eigenmode, the PDF of the normalized passive scalar reaches a self - similar profile within a finite time. This self - similarity indicates that the system has reached a non - equilibrium steady state, which is characterized by the combined effect of diffusion limitation and the "folding and filling" effect within a finite domain.
4. **Main Problems**
- **Shape of the PDF**: How does the shape of the PDF change under different Péclet number conditions? In particular, for small Péclet numbers, the PDF exhibits stretching or a pure exponential function; while for large Péclet numbers, the PDF tends to be Gaussian.
- **Conditional Statistics**: In the strange eigenmode, how do the characteristics of conditional diffusion and conditional dissipation affect the PDF?
- **Numerical Verification**: Verify the theoretical predictions through numerical simulations, especially the PDF behavior under different diffusion coefficients and flow intensity conditions.
### Numerical Simulation and Results
The author verified the above - mentioned theoretical predictions through numerical simulations, using two - dimensional lattice mapping and the diffusion approximation in Fourier space to solve the AD equation. The results show that:
- For small diffusion coefficients (\(\kappa\)), the core part of the PDF tends to be Gaussian.
- For strong flow fields, the core part of the PDF also tends to be Gaussian.
- These results support the Sinai and Yakhot theory's prediction that the PDF has a universal Gaussian core under large Péclet number conditions.
In conclusion, this paper aims to gain an in - depth understanding of the PDF behavior of passive scalars under specific conditions and verify these behaviors through theoretical analysis and numerical simulations, providing important insights for fluid mechanics and turbulence research.