Aug 2026· Nature Communications· Vol 17· 0 citations· 51 references
Medicine
TL;DR
A compact nonlinear-manifold representation, based on Gaussian splatting, for solving partial differential equations with discontinuities and thin shear layers, enabling sharp features to be captured with far fewer parameters than conventional numerical or neural-network-based solvers.
Abstract
Conventional linear discretizations, including high-order schemes, often require prohibitively many degrees of freedom to resolve sharp, localized features, and the practical advantages of nonlinear, physics-informed representations over fixed linear spaces remain unclear. Here we introduce a compact nonlinear-manifold representation, based on Gaussian splatting, for solving partial differential equations with discontinuities and thin shear layers. Mobile, anisotropic Gaussian kernels are evolved or optimized to satisfy the governing equations: time-dependent problems through a projection-based scheme on the manifold, and steady forward and inverse problems through direct residual minimization. The representation offers closed-form spatial derivatives and a highly compact description of sharp features. Across smooth, shock-dominated and shear-layer benchmarks, it attains accuracy comparable to classical solvers for smooth flows while remaining markedly more compact in strongly nonlinear regimes, and it displays an approximate error invariance as the Reynolds number increases—departing from the error accumulation inherent to linear representations. Simulating shocks and thin shear layers typically requires extremely fine computational grids. Here, the authors adapt mobile, stretchable Gaussian kernels originally developed in computer graphics to directly satisfy the governing equations, enabling sharp features to be captured with far fewer parameters than conventional numerical or neural-network-based solvers.
We use a large database of direct numerical simulations to investigate the transition of the Rayleigh--Taylor instability to turbulence and its evolution toward a late-time self-similar regime. In addition to tracking the growth of the mixing layer through the mean heavy-fluid concentration profile, we analyze one-dime...
Téo Granger, B. Nadiga, B. Gréa et al.· 0 citations
Reduced-order models (ROMs) are widely used to accelerate the solution of parametrized partial differential equations, with projection-based methods relying on a linear representation of the governing equations in a low-dimensional space. Their efficiency, however, is fundamentally limited by the Kolmogorov N-width of...
Neural operators provide efficient surrogates for spatiotemporal PDE systems, but purely data-driven formulations often accumulate substantial errors during long-horizon autoregressive prediction and may fail to exploit available governing-equation structure. Existing approaches incorporate physics primarily through re...
In this work, we develop a second-order optimization framework for physics-informed neural networks (PINNs) applied to high-dimensional parametric partial differential equations (PDEs). The framework is built on the Gauss--Newton pullback metric, which provides an operator-informed notion of curvature in parameter spac...
This work proposes the Fourier-enhanced alternating Levenberg--Marquardt PINN (FALM-PINN), an optimization framework that decouples representation learning from coefficient fitting within a single nonconvex optimization objective.
Yu-Lun Wu, Matthieu Barreau, Miguel Aguiar et al.· 0 citations
This work introduces a common framework based on the discretization of functional Gauss--Newton problems by finite families of linear measurements and shows that, through an appropriate duality pairing, the linear measurements can be represented by test functions.
Nilo Schwencke, Roland Maier· 0 citations
We use cookies to run the site and, with your consent, for analytics and to show ads.
See our Cookie Policy.