Proper generalized decomposition with synthetic preconditioning for multiscale parametric phonon transport
Abstract
Deterministic solution of the phonon Boltzmann transport equation is costly for two reasons: the distribution depends on high-dimensional phase space, and source iteration converges slowly when scattering is frequent. We address both for the steady linearized gray Callaway model in two dimensions by combining proper generalized decomposition (PGD) with the general synthetic iterative scheme (GSIS). The distribution is a sum of products of a space--azimuth factor and a polar factor. This splitting is intrinsically low-rank: along each characteristic the polar angle only stretches the in-plane optical depth, so the rank is bounded, uniformly in the Knudsen numbers, by the numerical rank of the kernel exp(−τ/sinϕ). PGD--GSIS updates the factors against a macroscopic source frozen within each alternating step and refreshed by the GSIS moment equations. A Fourier analysis shows that pure PGD, although implicit in the polar factor, moves macroscopic information only through lagged spatial sweeps and stalls like source iteration, whereas the synthetic refresh sees only third-order moments of the representation error. On an isothermal cavity, cross-plane conduction and phonon Poiseuille flow, PGD--GSIS reproduces full-order solutions at ranks 3--14, with heat-flux errors below 6×10−3 from the diffusive to the ballistic regime and memory independent of the polar resolution. Treating the inverse Knudsen numbers as coordinates of the polar factor yields one spatial basis over four decades of both; a joint update of all polar factors, which costs only small nodal linear systems, is required for accuracy at large Knudsen numbers.
Related
- Quantum Domain Decomposition for Preconditioning the Finite Element Method
- Mollified-sharp decomposition: a probabilistic regularization of parametric POD for shock-bearing flows
- An implicit multiscale solver for steady kinetic equations based on wave-particle decomposition
- A wave-particle decomposition framework for multiscale kinetic transport: continuous-spectrum equations and coupled wave-particle iteration
- A new parametric study of f(T) teleparallel gravity with generalized Chaplygin gas: Confronting observational data
BibTeX
@misc{liu2026proper,
title = {{Proper generalized decomposition with synthetic preconditioning for multiscale parametric phonon transport}},
author = {Jia Liu and Wei Su},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.27880},
eprint = {2610.27880},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.27880},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Proper generalized decomposition with synthetic preconditioning for multiscale parametric phonon transport
AU - Jia Liu
AU - Wei Su
PY - 2026
DA - 2026/10/20
PB - arXiv
DO - 10.48550/arxiv.2610.27880
UR - https://arxiv.org/abs/2610.27880
AB - Deterministic solution of the phonon Boltzmann transport equation is costly for two reasons: the distribution depends on high-dimensional phase space, and source iteration converges slowly when scattering is frequent. We address both for the steady linearized gray Callaway model in two dimensions by combining proper generalized decomposition (PGD) with the general synthetic iterative scheme (GSIS). The distribution is a sum of products of a space--azimuth factor and a polar factor. This splitting is intrinsically low-rank: along each characteristic the polar angle only stretches the in-plane optical depth, so the rank is bounded, uniformly in the Knudsen numbers, by the numerical rank of the kernel exp(−τ/sinϕ). PGD--GSIS updates the factors against a macroscopic source frozen within each alternating step and refreshed by the GSIS moment equations. A Fourier analysis shows that pure PGD, although implicit in the polar factor, moves macroscopic information only through lagged spatial sweeps and stalls like source iteration, whereas the synthetic refresh sees only third-order moments of the representation error. On an isothermal cavity, cross-plane conduction and phonon Poiseuille flow, PGD--GSIS reproduces full-order solutions at ranks 3--14, with heat-flux errors below 6×10^−3 from the diffusive to the ballistic regime and memory independent of the polar resolution. Treating the inverse Knudsen numbers as coordinates of the polar factor yields one spatial basis over four decades of both; a joint update of all polar factors, which costs only small nodal linear systems, is required for accuracy at large Knudsen numbers.
ER -