Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces Spectral-like Neural Discretisation (SpeND), a novel mesh-free method that utilizes neural networks to optimize stencil weights for discrete differential operators, ensuring polynomial consistency in a way that is independent of reference solutions. By framing the selection of degrees of freedom as a learning problem, SpeND achieves superior accuracy in approximating the modal response of spectral operators across a wider frequency band compared to traditional methods like LABFM and finite differences. The method demonstrates a significant improvement in dispersion and dissipation error while maintaining the expected convergence rates, showcasing its potential for enhanced computational efficiency in fluid dynamics simulations.
Neural networks can optimize mesh-free discretizations to achieve unprecedented accuracy across a wider frequency band than traditional methods.
Meshfree methods such as smoothed particle hydrodynamics (SPH) with kernel corrections, radial basis function-generated finite differences (RBF-FD), and the local anisotropic basis function method (LABFM) construct discrete differential operators by imposing polynomial consistency on a local stencil. For stencils containing more nodes than there are consistency constraints, the resulting linear system is underdetermined, and the remaining degrees of freedom are fixed implicitly by the choice of kernel, basis preconditioning, or a minimum-norm condition. Polynomial consistency constrains the operator only in the low-wavenumber limit, and no part of the construction selects for accuracy at the wavenumbers where fine-scale content resides. We introduce Spectral-like Neural Discretisation (SpeND), in which the choice of those degrees of freedom is cast as a learning problem: stencil weights are parametrised by a neural network conditioned on the local node geometry, trained to approximate the modal response of a spectral operator over the resolvable band. A hard-constrained projection layer maps the network output onto the affine subspace of consistent weights, so that polynomial consistency holds exactly by construction rather than as a penalty. Training is self-supervised and physics-agnostic, requiring no reference solutions; the objective minimises dispersion and dissipation error over a prescribed band-limited function space. Modal analysis on disordered two-dimensional node distributions shows that the learned fourth-order operator follows the exact response over a substantially wider band than either explicit LABFM at equal stencil size or fourth-order finite differences on a structured grid, whilst recovering the expected fourth-order convergence rate under refinement.