Mollified-sharp decomposition: a probabilistic regularization of parametric POD for shock-bearing flows
Abstract
This paper introduces the mollified-sharp decomposition, a probabilistic regularization of moving shocks in parametric reduced-order models. Each detected shock location is treated as a random variable with a prescribed probability density. Averaging over this artificial distribution replaces the localized pressure change by a smooth transition whose spatial extent is set by the density width, rather than by direct filtering of the pressure field. Each snapshot is decomposed exactly into a regularized mollified field and a local sharp correction that restores the shock. This probabilistic construction and exact additive split define the general method; the detector, kernel, treatment of multiple shocks, alignment coordinates, and regression are implementation choices. In the present realization, a calibrated indicator detects shocks, a compactly supported Wendland kernel mollifies them, and a peak-normalized and, where necessary, partitioned weight derived from each shock-location probability density defines the centroid, principal axes, and scales of its local alignment domain. Separate POD-GPR models represent the mollified field and aligned corrections, with additional regressions for shock presence and alignment. The method is demonstrated on the transonic airfoil pressure data of Catalani et al. (2023) and compared with a standard POD-GPR model constructed from the same data and POD-energy criterion. The mollified-sharp model reduces the mean test-set relative L2 pressure error by 31.2% and the mean test-set surface-pressure-coefficient error by 33.2%, while also improving the predicted shock locations and pressure changes; the trade-off is a median online evaluation time 2.24 times as long. The construction is applicable in principle to other parameter-dependent fields with moving sharp features, such as moving material interfaces in multiphase flows.
Related
- Proper generalized decomposition with synthetic preconditioning for multiscale parametric phonon transport
- Deflection monitoring of discontinuous built-up beams using FBG measurements with mechanics-guided regularization
- Phase-field modeling of TiO2 nanocarving via reaction with hydrogen-bearing gas
- Functional Dimensional Regularization for the Lee-Yang Model
- Regularization of stochastic fluid dynamics and non-perturbative renormalization of the shear viscosity of incompressible fluids
BibTeX
@misc{schmidt2026mollifiedsharp,
title = {{Mollified-sharp decomposition: a probabilistic regularization of parametric POD for shock-bearing flows}},
author = {Oliver T. Schmidt},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.19532},
eprint = {2609.19532},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.19532},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Mollified-sharp decomposition: a probabilistic regularization of parametric POD for shock-bearing flows
AU - Oliver T. Schmidt
PY - 2026
DA - 2026/09/17
PB - arXiv
DO - 10.48550/arxiv.2609.19532
UR - https://arxiv.org/abs/2609.19532
AB - This paper introduces the mollified-sharp decomposition, a probabilistic regularization of moving shocks in parametric reduced-order models. Each detected shock location is treated as a random variable with a prescribed probability density. Averaging over this artificial distribution replaces the localized pressure change by a smooth transition whose spatial extent is set by the density width, rather than by direct filtering of the pressure field. Each snapshot is decomposed exactly into a regularized mollified field and a local sharp correction that restores the shock. This probabilistic construction and exact additive split define the general method; the detector, kernel, treatment of multiple shocks, alignment coordinates, and regression are implementation choices. In the present realization, a calibrated indicator detects shocks, a compactly supported Wendland kernel mollifies them, and a peak-normalized and, where necessary, partitioned weight derived from each shock-location probability density defines the centroid, principal axes, and scales of its local alignment domain. Separate POD-GPR models represent the mollified field and aligned corrections, with additional regressions for shock presence and alignment. The method is demonstrated on the transonic airfoil pressure data of Catalani et al. (2023) and compared with a standard POD-GPR model constructed from the same data and POD-energy criterion. The mollified-sharp model reduces the mean test-set relative L^2 pressure error by 31.2% and the mean test-set surface-pressure-coefficient error by 33.2%, while also improving the predicted shock locations and pressure changes; the trade-off is a median online evaluation time 2.24 times as long. The construction is applicable in principle to other parameter-dependent fields with moving sharp features, such as moving material interfaces in multiphase flows.
ER -