[Paper Review] On stiff problems via Dirichlet forms

Liping Li, Wenjie Sun|ArXiv.org|Apr 8, 2018

Stochastic processes and statistical mechanics19 references3 citations

TL;DR

This paper establishes a rigorous probabilistic and analytical framework for stiff problems in one-dimensional thermal conduction using Dirichlet forms, identifying three distinct phase transitions—impermeable, semi-permeable, and permeable—based on the total thermal resistance of a singular barrier. It introduces the snapping out Markov process as the probabilistic limit of diffusions with shrinking barriers and derives corresponding flux boundary conditions via Mosco convergence of Dirichlet forms.

ABSTRACT

The stiff problem is concerned with a thermal conduction model with a singular barrier of zero volume. In this paper, we shall build the phase transitions for the stiff problems in one-dimensional space. It turns out that every phase transition definitely depends on the total thermal resistance of the barrier, and the three phases correspond to the so-called impermeable pattern, semi-permeable pattern and permeable pattern of thermal conduction respectively. For each pattern, the related boundary condition of the flux at the barrier is also derived. Mathematically, we shall introduce and explore the so-called snapping out Markov process, which is the probabilistic counterpart of semi-permeable pattern in the stiff problem.

Motivation & Objective

• To rigorously characterize phase transitions in stiff problems arising from singular thermal barriers with zero volume.
• To establish a probabilistic counterpart—snapping out Markov processes—for semi-permeable thermal conduction via Dirichlet form theory.
• To derive explicit flux boundary conditions at the barrier point based on the total thermal resistance.
• To unify the limiting behavior of heat equations with shrinking barriers using Mosco convergence of Dirichlet forms.
• To generalize Lejay’s resolvent-based approach to a more systematic and generalizable framework using symmetric Markovian Dirichlet forms.

Proposed method

• Utilizes symmetric Markovian Dirichlet forms on $L^2(bR, m)$ to represent diffusion processes and their generators.
• Applies Mosco convergence to analyze the limit of Dirichlet forms associated with shrinking barriers of small conductivity.
• Introduces the snapping out Brownian motion (SNOB) as the probabilistic limit of diffusions with shrinking barriers.
• Derives boundary conditions for the flux at the barrier by analyzing the generator and resolvent of the limiting process.
• Uses killing, time change, darning, and piecing out transforms to construct and analyze the limiting Markov process.
• Applies the Hille-Yosida theorem and properties of quasi-regular Dirichlet forms to ensure existence and regularity of solutions.

Experimental results

Research questions

• RQ1 How do different levels of total thermal resistance in a singular barrier lead to distinct phase transitions in stiff problems?
• RQ2 What is the probabilistic limit of diffusion processes with shrinking barriers of vanishing conductivity?
• RQ3 How can the snapping out Markov process be rigorously constructed via Dirichlet form theory?
• RQ4 What are the precise flux boundary conditions at the barrier point in each phase, and how do they depend on the total resistance?
• RQ5 How does Mosco convergence of Dirichlet forms ensure the convergence of solutions to the heat equation with singular barriers?

Key findings

• Three distinct phases emerge based on the total thermal resistance $\bar{\gamma}$: $\bar{\gamma} = \infty$ (impermeable), $0 < \bar{\gamma} < \infty$ (semi-permeable), and $\bar{\gamma} = 0$ (permeable).
• For $\bar{\gamma} = \infty$, the flux satisfies $u'(0+) = u'(0-) = 0$, corresponding to a reflecting (impermeable) boundary condition.
• For $0 < \bar{\gamma} < \infty$, the flux satisfies $a(0+)u'(0+) = a(0-)u'(0-) = \frac{\kappa}{2}(u(0+) - u(0-))$ with $\kappa = 2/\bar{\gamma}$, representing a semi-permeable interface.
• For $\bar{\gamma} = 0$, the solution $u_t$ is continuous at $0$, and the weak solution is unique in $\mathscr{H}(\bbR)$, indicating full permeability.
• The snapping out Markov process arises as the weak limit of diffusions with shrinking barriers and corresponds to the semi-permeable phase.
• The boundary condition in the semi-permeable case reduces to the impermeable one in the limit $\bar{\gamma} \to \infty$, confirming consistency across phases.