Modified Kalman Filtering Derived from Non-Maxwellian Distribution Functions in Open Systems
Abstract
Kalman filtering (KF) recursively infers plasma quantities, represented by a state, from noisy diagnostics while propagating uncertainty in the inferred state separately from diagnostic noise. For linear dynamical and measurement models with Gaussian probability density functions (PDFs), the state mean and covariance provide the KF description. We modify this KF for open plasmas with particle and energy sources by extending the physical state from Maxwellian variables to retained NMDF coordinates, whose evolution follows a projection of the nonlinear Landau-Fokker-Planck equation. The same NMDF state is propagated through a fixed diagnostic response to obtain the corresponding measurement PDF. Because a non-Gaussian likelihood can drive the posterior outside the Gaussian family, the mean-covariance representation is extended to a finite moment closure whose coordinates are related to the retained moments through their Jacobian. The recursive prediction-correction structure is preserved, recovering conventional KF results for linear models with Gaussian PDFs. As a proof of concept, we use seven published non-Gaussian Alcator C-Mod Langmuir-probe current PDFs. For each target PDF, NMDF response parameters are calibrated from the other probe PDFs and frozen; this excluded target tests independent prediction. This test shows that different plasma regions require different NMDF structures, instead of the standard assumption of Maxwellian velocity distributions across all regimes and locations. Because the published measurement PDFs do not retain time ordering, this proof of concept does not yet test recursive uncertainty dynamics. Implementing the projected Landau-Fokker-Planck prediction with time-resolved diagnostics is the next step toward full validation of the modified KF.
Related
- Bayesian inference of non-Maxwellian distribution functions from collective Thomson scattering spectra
- Recovering Non-Maxwellian Velocity Distribution Functions from Collective Thomson Scattering Spectra Using Generalized Basis Representations
- High-precision fault inception time detection in active distribution networks using Euclidean distance and Gaussian filtering
- Unconditional authentication in quantum key distribution via hybrid entangled physical unclonable functions
- Final-state effects on the transverse-momentum spectra of charged hadrons in p+Pb and Pb+Pb collisions at √(sNN) = 5.02~TeV using the modified Tsallis distribution
BibTeX
@misc{izacard2026modified,
title = {{Modified Kalman Filtering Derived from Non-Maxwellian Distribution Functions in Open Systems}},
author = {Olivier Izacard},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2608.22050},
eprint = {2608.22050},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.22050},
}
RIS
TY - UNPB
TI - Modified Kalman Filtering Derived from Non-Maxwellian Distribution Functions in Open Systems
AU - Olivier Izacard
PY - 2026
DA - 2026/09/10
PB - arXiv
DO - 10.48550/arxiv.2608.22050
UR - https://arxiv.org/abs/2608.22050
AB - Kalman filtering (KF) recursively infers plasma quantities, represented by a state, from noisy diagnostics while propagating uncertainty in the inferred state separately from diagnostic noise. For linear dynamical and measurement models with Gaussian probability density functions (PDFs), the state mean and covariance provide the KF description. We modify this KF for open plasmas with particle and energy sources by extending the physical state from Maxwellian variables to retained NMDF coordinates, whose evolution follows a projection of the nonlinear Landau-Fokker-Planck equation. The same NMDF state is propagated through a fixed diagnostic response to obtain the corresponding measurement PDF. Because a non-Gaussian likelihood can drive the posterior outside the Gaussian family, the mean-covariance representation is extended to a finite moment closure whose coordinates are related to the retained moments through their Jacobian. The recursive prediction-correction structure is preserved, recovering conventional KF results for linear models with Gaussian PDFs. As a proof of concept, we use seven published non-Gaussian Alcator C-Mod Langmuir-probe current PDFs. For each target PDF, NMDF response parameters are calibrated from the other probe PDFs and frozen; this excluded target tests independent prediction. This test shows that different plasma regions require different NMDF structures, instead of the standard assumption of Maxwellian velocity distributions across all regimes and locations. Because the published measurement PDFs do not retain time ordering, this proof of concept does not yet test recursive uncertainty dynamics. Implementing the projected Landau-Fokker-Planck prediction with time-resolved diagnostics is the next step toward full validation of the modified KF.
ER -