English

Automated dimensional analysis for PDEs

Mathematical Software 2026-01-13 v1 Numerical Analysis Numerical Analysis

Abstract

Physical units are fundamental to scientific computing. However, many finite element frameworks lack built-in support for dimensional analysis. In this work, we present a systematic framework for integrating physical units into the Unified Form Language (UFL). We implement a symbolic Quantity class to track units within variational forms. The implementation exploits the abelian group structure of physical dimensions. We represent them as vectors in Qn\mathbb{Q}^n to simplify operations and improve performance. A graph-based visitor pattern traverses the expression tree to automate consistency checks and factorization. We demonstrate that this automated nondimensionalization functions as the simplest form of Full Operator Preconditioning. It acts as a physics-aware diagonal preconditioner that equilibrates linear systems prior to assembly. Numerical experiments with the Navier--Stokes equations show that this improves the condition number of the saddle-point matrix. Analysis of Neo-Hooke hyperelasticity highlights the detection of floating-point cancellation errors in small deformation regimes. Finally, the Poisson--Nernst--Planck system example illustrates the handling of coupled multiphysics problems with derived scaling parameters. Although the implementation targets the FEniCSx framework, the concepts are general and easily adaptable to other finite element libraries using UFL, such as Firedrake or DUNE.

Keywords

Cite

@article{arxiv.2601.06535,
  title  = {Automated dimensional analysis for PDEs},
  author = {Michal Habera and Andreas Zilian},
  journal= {arXiv preprint arXiv:2601.06535},
  year   = {2026}
}

Comments

29 pages

R2 v1 2026-07-01T08:58:55.255Z