Efficient Bernstein Quasi-Trefftz Discretization for Heterogeneous High-Frequency Helmholtz Problems
arXiv:2609.05849
Abstract
We develop an efficient Bernstein quasi-Trefftz discretization for heterogeneous high-frequency Helmholtz problems on unstructured triangular meshes. The method first imposes strong continuity through Bernstein--Bézier smoothness relations on local macro-patches and then reduces the resulting conforming polynomial space by enforcing projected Helmholtz residual constraints. On two-triangle patches, this compresses the local dimension from quadratic growth in the polynomial degree to a trace-sized space of dimension under the natural rank condition. Variable matrix-valued coefficients are handled through local polynomial projection, allowing coefficient-approximation effects to be separated from discretization error. A graph-residual formulation couples the reduced patch spaces and produces a sparse global system in compressed coordinates. The implementation combines explicit continuation, batched coefficient projection, batched residual construction, and stable local kernel extraction by QR factorization. Numerical experiments on heterogeneous unstructured meshes demonstrate high-order accuracy, roundoff-level conformity, robust high-frequency resolution, and substantial reductions in local and global computational cost. Turning-point and penetrable-scattering examples further illustrate the flexibility of the approach.