--- title: Quasilinear Diffusion Theory url: https://www.emergentmind.com/topics/quasilinear-diffusion-theory type: topic --- # Quasilinear Diffusion Theory Quasilinear diffusion theory denotes a family of diffusion frameworks in which the effective transport law is not fixed a priori but is generated by the state of the system or by an averaged interaction mechanism. In plasma kinetics, it describes the slow evolution of a distribution function under weak, resonant wave–particle interactions, yielding diffusion equations in velocity, momentum, or action space. In nonlinear PDE theory, it designates equations whose highest-order part depends on the unknown or its derivatives, such as $-\operatorname{div}\!\big(\mu(|\nabla u|)\nabla u\big)=f$, together with nonlocal, stochastic, kinetic, and weighted reaction–diffusion variants [2208.09477] [2101.10137] [2206.11415]. ## 1. Foundational kinetic formulation In the standard kinetic setting, quasilinear diffusion arises from a weakly perturbed Vlasov system after separating a slowly evolving background distribution from fast oscillatory fluctuations. For a uniform magnetized plasma, the quasilinear evolution of a gyrotropic distribution $f_0(v_\parallel,v_\perp)$ can be written in divergence form, \[ \frac{1}{\Omega}\frac{\partial f_0}{\partial \tau} = \frac{\partial}{\partial \mathbf{v}}\cdot \left(\mathbf{D}\cdot \frac{\partial f_0}{\partial \mathbf{v}}\right), \] with diffusion tensor \[ \mathbf{D} = \sum_{\ell}\mathrm{Re}(-i\Delta_\ell)\, \tilde{\mathbf{v}}_\ell^\ast\tilde{\mathbf{v}}_\ell, \qquad \mathrm{Re}(-i\Delta_\ell) = \pi\Omega\,\delta(\omega-k_\parallel v_\parallel-\ell\Omega). \] The resonance condition \[ \omega-k_\parallel v_\parallel-\ell\Omega=0 \] is therefore the organizing principle of the theory: only resonant particles contribute to secular diffusion [2208.09477]. The classical assumptions are weak perturbations, random phases, separation of time scales, and dominance of resonant interactions. In this regime, the cumulative effect of many small wave-induced kicks is diffusive rather than ballistic. In invariant velocity variables for a uniform field, the natural coordinates are $p_\parallel=Mv_\parallel$ and $\mu=Mv_\perp^2/(2B_0)$, and the standard Kennel–Engelmann equation becomes a symmetric $2\times2$ diffusion system in $(p_\parallel,\mu)$ or, equivalently, in $(\mathcal E,\mu)$, where the energy kick depends only on the electric field [2208.09477]. An analogous quasilinear closure appears for long-range Vlasov systems. For the Hamiltonian Mean Field model, one writes \[ f(\theta,p,t)=f_0(p,t)+f_1(\theta,p,t), \] linearizes the Vlasov equation around the slowly varying homogeneous part $f_0$, and derives a diffusion equation \[ \frac{\partial f_0}{\partial t} = \frac{\partial}{\partial p} \left[D(p,t)\frac{\partial f_0}{\partial p}\right]. \] For even single-humped distributions in the HMF model, only the $k=\pm1$ modes matter, the unstable eigenfrequency is purely imaginary, and the diffusion coefficient reduces to \[ D(p,t)=\frac{2\chi(t)\omega_I(t)}{p^2+\omega_I^2(t)}. \] This quasilinear theory works reasonably well for weakly unstable initial conditions and predicts the energy marking the out-of-equilibrium phase transition between unmagnetized and magnetized quasi-stationary states [1607.04165]. ## 2. Hamiltonian and invariant-space reformulations A central development in modern plasma quasilinear theory is its Hamiltonian reformulation. Instead of expressing the wave forcing directly through $(\delta \mathbf E,\delta \mathbf B)$, one introduces the gauge-invariant effective potential \[ \delta\Psi=\delta\Phi-\frac{\mathbf v}{c}\cdot\delta\mathbf A, \qquad \delta H=q\,\delta\Psi, \] and rewrites the Vlasov equation with a noncanonical Poisson bracket. In this representation, the second-order quasilinear evolution is generated by \[ \frac{\partial f_0}{\partial \tau} = \frac12\left\langle\overline{\{\delta H,\delta g\}}\right\rangle, \] which makes gauge invariance, entropy production, and the dyadic structure of the diffusion tensor explicit [2208.09477]. In a uniform magnetic field, the natural invariant space is two-dimensional. In a nonuniform axisymmetric magnetic field, the appropriate invariants become \[ \mathbf J=({\sf J}_{\rm g},\mathcal E,{\sf J}_{\rm d}), \qquad {\sf J}_{\rm d}=q\psi/c, \] so quasilinear transport is intrinsically three-dimensional in invariant space. The corresponding evolution equation is \[ \frac{\partial F_0}{\partial \tau} = \frac{1}{\tau_{\rm b}} \frac{\partial}{\partial {\sf J}^i} \left( \tau_{\rm b} D_{\rm QL}^{ij} \frac{\partial F_0}{\partial {\sf J}^j} \right), \] and the diffusion tensor has the explicit $3\times3$ form \[ \mathbf D_{\rm QL} = \sum_{\ell,k,m} \begin{pmatrix} \ell^2 & \ell\omega_k & \ell m\\ \omega_k\ell & \omega_k^2 & \omega_k m\\ m\ell & m\omega_k & m^2 \end{pmatrix} \Gamma_{\ell km}. \] This tensor couples gyro/pitch-angle transport, energy diffusion, and radial diffusion in a single object [2208.09477]. The same structure is obtained abstractly in canonical action–angle variables: \[ \frac{\partial F_0}{\partial \tau} = \frac{\partial}{\partial \mathbf J}\cdot \left( \mathbf D_{\rm QL}\cdot \frac{\partial F_0}{\partial \mathbf J} \right), \qquad \mathbf D_{\rm QL} = \sum_{\mathbf m,k} \mathbf m\mathbf m\, \pi\delta(\omega_k-\mathbf m\cdot\boldsymbol\Omega)\, |\delta\tilde{\mathcal H}_{\mathbf m,k}(J)|^2. \] This canonical form shows that diffusion is along resonance directions $\mathbf m$ in action space and that the tensor is a sum of dyads weighted by resonant quasilinear potentials [2208.09477]. A closely related conservation-law refinement appears in the extension of the Kennel–Engelmann tensor from two to four dimensions. The standard magnetized tensor respects the energy–parallel-momentum relation but neglects perpendicular momentum absorption. Enforcing four-momentum conservation leads to \[ D^{\mu\nu} = D^{KK}\frac{k^\mu k^\nu}{\omega^2} = D_0\,\frac{k^\mu k^\nu v_\perp^2}{\omega^2}, \] which, after transformation to constants-of-motion variables $(\epsilon,\mu,p_\phi)$, yields the diffusion path \[ \mathbf w = v_\perp \begin{pmatrix} 1\\[2pt] \dfrac{n\Omega}{|B|\omega}\\[4pt] \dfrac{n_\phi}{\omega} \end{pmatrix}. \] This matches the form required by action-angle Hamiltonian theory after bounce averaging [2511.09532]. ## 3. Positivity, conservation, and numerical structure For toroidal plasmas, bounce averaging of the classical Kennel–Engelmann coefficients can destroy positive definiteness because the parallel inhomogeneity makes the resonance kernel asymmetric under the bounce average. A positive-definite alternative is obtained by evaluating the phase integral along the trajectory before averaging. The resulting bounce-averaged diffusion tensor can be written as \[ \mathcal D(\mathbf v_c) = \frac{q^2}{2m^2 t_d}\sum_n \mathbf T_n^\ast \mathbf T_n, \] which is manifestly positive semidefinite. In this form, resonant contributions are expressed through Airy-function factors near resonance points, while nonresonant contributions arise near outer-midplane, inner-midplane, and trapping-tip locations where the phase curvature is small. The construction includes both resonant and non-resonant contributions, and the correlations between the consecutive resonances and in many poloidal periods [1707.08030]. A related finite-Larmor-radius reduction has been developed for ion cyclotron heating. Starting from the kinetic energy change $\dot W$ and matching it to the quasilinear Fokker–Planck operator, one derives reduced coefficients that preserve the diffusion directions, wave polarizations, and H-theorem of the full Kennel–Engelmann model. For the fundamental damping, the lowest-order reduction corresponds to $J_0\to1$ and neglect of higher Bessel terms; for the second harmonic damping, the lowest nonzero contribution is $O((k_\perp\rho_i)^2)$ and corresponds to $J_1\simeq k_\perp\rho_i/2$ [1704.07283]. Conservation can also be enforced at the discretization level. In a three-dimensional momentum / three-dimensional spectral quasilinear model with cylindrical symmetry, the particle pdf and wave sed satisfy coupled bilinear equations that admit an unconditionally conservative weak form. A conservative Galerkin discretization with continuous basis functions for the particle pdf, discontinuous basis functions for the wave sed, and a consistent quadrature rule preserves particle number, momentum, and energy rigorously, independently of the singular transition probability. The resonance manifold is integrated by a marching simplex algorithm, which converts the singular kernel into quadrature over simplices on the manifold [2212.07229]. These developments show that positivity and conservation are not merely numerical conveniences. They encode the structural content of quasilinear diffusion: entropy production, compatibility with wave action, and preservation of the geometric diffusion directions determined by resonance. ## 4. Variational and elliptic-parabolic PDE formulations In nonlinear PDE theory, quasilinear diffusion refers to equations whose principal part depends on the unknown or on its gradient. A prototypical elliptic model is \[ -\operatorname{div}\big(\mu(|\nabla u|)\nabla u\big)=f \quad\text{in }\Omega, \] with mixed boundary conditions on $\Gamma_1\cup\Gamma_2$. In the variational setting, \[ a(u;v,w) = \int_\Omega \mu(|\nabla u|)\,\nabla v\cdot\nabla w\,dx, \qquad H(u)=\int_\Omega \psi(|\nabla u|^2)\,dx-\langle b,u\rangle, \] and the weak problem is equivalent to minimizing $H$ over the admissible affine set $K$ [2101.10137]. The classical Kačanov scheme freezes the diffusion law at the previous iterate: \[ a(u^n;u^{n+1},v)=\langle b,v\rangle \quad \forall v\in X, \] or, equivalently, \[ u^{n+1}=u^n-A(u^n)^{-1}H'(u^n). \] Its traditional convergence theory assumes that $\mu$ is continuously differentiable, decreasing, uniformly bounded above and below, and associated with a uniformly convex energy. The key limitation is the monotonicity requirement $\mu'(t)\le 0$, which excludes shear-thickening and non-monotone laws even though the iteration often converges numerically [2101.10137]. A modified Kačanov iteration introduces damping, \[ u^{n+1} = u^n-\delta(u^n)A(u^n)^{-1}H'(u^n), \] and proves convergence from an energy decrease estimate \[ H(u^n)-H(u^{n+1}) \ge \gamma\|u^{n+1}-u^n\|_Y^2. \] For quasilinear diffusion, this estimate follows from the two-sided structural condition \[ m_\mu (t-s)\le \mu(t^2)t-\mu(s^2)s\le M_\mu (t-s), \qquad t\ge s\ge 0, \] which yields strong monotonicity and Lipschitz continuity of $H'$ without requiring $\mu$ to be decreasing or differentiable. The resulting solver converges for decreasing, non-monotone, and strongly non-monotone viscosity laws, and adaptive damping can both accelerate convergence and recover convergence where the classical Kačanov iteration fails [2101.10137]. A time-dependent quasilinear reaction–diffusion model displays a different aspect of the theory: \[ \partial_t u = \Delta u^m + (1+|x|)^\sigma u^p, \qquad N\ge 3,\quad 1