Fehmi Cirak

Arnaud Vadeboncoeur, Ömer Deniz Akyildiz, Ieva Kazlauskaite et al.

We introduce a physics-driven deep latent variable model (PDDLVM) to learn simultaneously parameter-to-solution (forward) and solution-to-parameter (inverse) maps of parametric partial differential equations (PDEs). Our formulation leverages conventional PDE discretization techniques, deep neural networks, probabilistic modelling, and variational inference to assemble a fully probabilistic coherent framework. In the posited probabilistic model, both the forward and inverse maps are approximated as Gaussian distributions with a mean and covariance parameterized by deep neural networks. The PDE residual is assumed to be an observed random vector of value zero, hence we model it as a random vector with a zero mean and a user-prescribed covariance. The model is trained by maximizing the probability, that is the evidence or marginal likelihood, of observing a residual of zero by maximizing the evidence lower bound (ELBO). Consequently, the proposed methodology does not require any independent PDE solves and is physics-informed at training time, allowing the real-time solution of PDE forward and inverse problems after training. The proposed framework can be easily extended to seamlessly integrate observed data to solve inverse problems and to build generative models. We demonstrate the efficiency and robustness of our method on finite element discretized parametric PDE problems such as linear and nonlinear Poisson problems, elastic shells with complex 3D geometries, and time-dependent nonlinear and inhomogeneous PDEs using a physics-informed neural network (PINN) discretization. We achieve up to three orders of magnitude speed-up after training compared to traditional finite element method (FEM), while outputting coherent uncertainty estimates.

Arnaud Vadeboncoeur, Ieva Kazlauskaite, Yanni Papandreou et al.

We introduce a new class of spatially stochastic physics and data informed deep latent models for parametric partial differential equations (PDEs) which operate through scalable variational neural processes. We achieve this by assigning probability measures to the spatial domain, which allows us to treat collocation grids probabilistically as random variables to be marginalised out. Adapting this spatial statistics view, we solve forward and inverse problems for parametric PDEs in a way that leads to the construction of Gaussian process models of solution fields. The implementation of these random grids poses a unique set of challenges for inverse physics informed deep learning frameworks and we propose a new architecture called Grid Invariant Convolutional Networks (GICNets) to overcome these challenges. We further show how to incorporate noisy data in a principled manner into our physics informed model to improve predictions for problems where data may be available but whose measurement location does not coincide with any fixed mesh or grid. The proposed method is tested on a nonlinear Poisson problem, Burgers equation, and Navier-Stokes equations, and we provide extensive numerical comparisons. We demonstrate significant computational advantages over current physics informed neural learning methods for parametric PDEs while improving the predictive capabilities and flexibility of these models.

Kosala Bandara, Fehmi Cirak

We introduce the isogeometric shape optimisation of thin shell structures using subdivision surfaces. Both triangular Loop and quadrilateral Catmull-Clark subdivision schemes are considered for geometry modelling and finite element analysis. A gradient-based shape optimisation technique is implemented to minimise compliance, i.e. to maximise stiffness. Different control meshes describing the same surface are used for geometry representation, optimisation and finite element analysis. The finite element analysis is performed with subdivision basis functions corresponding to a sufficiently refined control mesh. During iterative shape optimisation the geometry is updated starting from the coarsest control mesh and proceeding to increasingly finer control meshes. This multiresolution approach provides a means for regularising the optimisation problem and prevents the appearance of sub-optimal jagged geometries with fine-scale oscillations. The finest control mesh for optimisation is chosen in accordance with the desired smallest feature size in the optimised geometry. The proposed approach is applied to three optimisation examples, namely a catenary, a roof over a rectangular domain and a freeform architectural shell roof. The influence of the geometry description and the used subdivision scheme on the obtained optimised curved geometries is investigated in detail.

Musabbir Majeed, Fehmi Cirak

We present an isogeometric analysis technique that builds on manifold-based smooth basis functions for geometric modelling and analysis. Manifold-based surface construction techniques are well known in geometric modelling and a number of variants exist. Common to all is the concept of constructing a smooth surface by blending together overlapping patches (or, charts), as in differential geometry description of manifolds. Each patch on the surface has a corresponding planar patch with a smooth one-to-one mapping onto the surface. In our implementation, manifold techniques are combined with conformal parametrisations and the partition-of-unity method for deriving smooth basis functions on unstructured quadrilateral meshes. Each vertex and its adjacent elements on the surface control mesh have a corresponding planar patch of elements. The star-shaped planar patch with congruent wedge-shaped elements is smoothly parameterised with copies of a conformally mapped unit square. The conformal maps can be easily inverted in order to compute the transition functions between the different planar patches that have an overlap on the surface. On the collection of star-shaped planar patches the partition of unity method is used for approximation. The smooth partition of unity, or blending functions, are assembled from tensor-product b-spline segments defined on a unit square. On each patch a polynomial with a prescribed degree is used as a local approximant. To obtain a mesh-based approximation scheme, the coefficients of the local approximants are expressed in dependence of vertex coefficients. This yields a basis function for each vertex of the mesh which is smooth and non-zero over a vertex and its adjacent elements. Our numerical simulations indicate the optimal convergence of the resulting approximation scheme for Poisson problems and near optimal convergence for thin-plate and thin-shell problems.

Eky Febrianto, Yiren Wang, Burigede Liu et al.

Quantum computing holds the promise of solving computational mechanics problems in polylogarithmic time, meaning computational time scales as $\mathscr{O}((\log N)^c)$, where $N$ is the problem size and $c$ a constant. We propose a quantum spectral method with polylogarithmic complexity for solving non-periodic boundary value problems with arbitrary Dirichlet boundary conditions. Our method extends the recently proposed approach by Liu et al. (2025), in which periodic problems are discretised using truncated Fourier series. In such spectral methods, the discretisation of boundary value problems with constant coefficients leads to a set of algebraic equations in the Fourier space. We implement the respective diagonal solution operator by first approximating it with a polynomial and then quantum encoding the polynomial. The mapping between the physical and Fourier spaces is accomplished using the quantum Fourier transform (QFT). To impose zero Dirichlet boundary conditions, we double the domain size and reflect all physical fields antisymmetrically. The respective reflection matrix defines the quantum sine transform (QST) by pre- and post-multiplying with the QFT. For non-zero Dirichlet boundary conditions, the solution is decomposed into a boundary-conforming and a homogeneous part. The homogenous part is determined by solving a problem with a suitably modified forcing vector. We illustrate the basic approach with a Dirichlet-Poisson problem and demonstrate its generality by applying it to a fractional stochastic PDE for modelling spatial random fields. We discuss the circuit implementation of the proposed approach and provide numerical evidence confirming its polylogarithmic complexity.

Igor Kavrakov, Yaswanth Sai Jetti, Ahmet Oguzhan Yuksel et al.

We present an approach for synthesising observational data with elastodynamic finite element models by extending the statistical finite element method (statFEM) framework. The proposed formulation adopts a Bayesian filtering approach to account for uncertainties in the data, the finite element model, and the discrepancies between the model and the physical system. Observational data are assimilated while the state of the spatially discretised finite element problem is advanced in time using the stochastic variant of the explicit Newmark scheme. The prior probability density of the state is obtained by solving an incremental probabilistic forward problem, and the corresponding posterior density is obtained by conditioning on the data available at each time step. In the probabilistic forward problem, spatio-temporal Gaussian random fields representing the forcing, model misspecification, and material parameters are specified via their stochastic PDE formulation. The resulting non-Gaussian prior probability density is approximated using a perturbation approach, yielding a Gaussian posterior with closed-form mean and covariance. The hyperparameters of the random field representing model misspecification are calibrated by maximising the marginal likelihood of the data. The proposed approach is illustrated on one- and two-dimensional elastodynamic examples with synthetic data.

Xiao Xiao, Fehmi Cirak

Shape optimisation of thin-shell structures requires a flexible, differentiable geometric representation suitable for gradient-based optimisation. We propose a neural parametric representation (NRep) for the shell mid-surface based on a neural network with periodic activation functions. The NRep is defined using a multi-layer perceptron (MLP), which maps the parametric coordinates of mid-surface vertices to their physical coordinates. A structural compliance optimisation problem is posed to optimise the shape of a thin-shell parameterised by the NRep subject to a volume constraint, with the network parameters as design variables. The resulting shape optimisation problem is solved using a gradient-based optimisation algorithm. Benchmark examples with classical solutions demonstrate the effectiveness of the proposed NRep. The approach exhibits potential for complex lattice-skin structures, owing to the compact and expressive geometry representation afforded by the NRep.

Yiren Wang, Michael Ortiz, Fehmi Cirak

The computational cost of concurrent multiscale finite element methods is dominated by the repeated solution of microscopic representative volume element (RVE) problems at macroscopic quadrature points. In this work, we introduce a quantum-classical framework for multiscale finite element analysis (QAFE$^2$) that leverages quantum parallelism to fundamentally alter the scaling of RVE-based homogenisation. At the single-RVE level, the proposed quantum solver attains polylogarithmic complexity with respect to the microscopic discretisation size, yielding an exponential asymptotic speedup over the best available classical solvers. More importantly, QAFE$^2$ exploits quantum superposition and entanglement to evaluate, in a single quantum execution, the entire ensemble of RVE problems associated with all macroscopic quadrature points. This capability is a form of intrinsic quantum concurrency with no classical analogue. Numerical experiments on one- and two-dimensional model problems with known analytical solutions confirm the accuracy of the proposed formulation and verify the theoretical computational scaling and parallel performance.

Arion Pons, Fehmi Cirak

Modern high-performance combat aircraft exceed conventional flight-envelope limits on maneuverability through the use of thrust vectoring, and so achieve supermaneuverability. With ongoing development of biomimetic unmanned aerial vehicles (UAVs), the potential for supermaneuverability through biomimetic mechanisms becomes apparent. So far, this potential has not been well studied: biomimetic UAVs have not yet been shown to be capable of any of the forms of classical supermaneuverability available to thrust-vectored aircraft. Here we show this capability, by demonstrating how biomimetic morphing-wing UAVs can perform sophisticated multiaxis nose-pointing-and-shooting (NPAS) maneuvers at low morphing complexity. Nonlinear flight-dynamic analysis is used to characterize the extent and stability of the multidimensional space of aircraft trim states that arises from biomimetic morphing. Navigating this trim space provides an effective model-based guidance strategy for generating open-loop NPAS maneuvers in simulation. Our results demonstrate the capability of biomimetic aircraft for air combat-relevant supermaneuverability, and provide strategies for the exploration, characterization, and guidance of further forms of classical and non-classical supermaneuverability in such aircraft.

Sumudu Herath, Xiao Xiao, Fehmi Cirak

Knitting is an effective technique for producing complex three-dimensional surfaces owing to the inherent flexibility of interlooped yarns and recent advances in manufacturing providing better control of local stitch patterns. Fully yarn-level modelling of large-scale knitted membranes is not feasible. Therefore, we use a two-scale homogenisation approach and model the membrane as a Kirchhoff-Love shell on the macroscale and as Euler-Bernoulli rods on the microscale. The governing equations for both the shell and the rod are discretised with cubic B-spline basis functions. For homogenisation we consider only the in-plane response of the membrane. The solution of the nonlinear microscale problem requires a significant amount of time due to the large deformations and the enforcement of contact constraints, rendering conventional online computational homogenisation approaches infeasible. To sidestep this problem, we use a pre-trained statistical Gaussian Process Regression (GPR) model to map the macroscale deformations to macroscale stresses. During the offline learning phase, the GPR model is trained by solving the microscale problem for a sufficiently rich set of deformation states obtained by either uniform or Sobol sampling. The trained GPR model encodes the nonlinearities and anisotropies present in the microscale and serves as a material model for the membrane response of the macroscale shell. The bending response can be chosen in dependence of the mesh size to penalise the fine out-of-plane wrinkling of the membrane. After verifying and validating the different components of the proposed approach, we introduce several examples involving membranes subjected to tension and shear to demonstrate its versatility and good performance.

Zhaowei Liu, Musabbir Majeed, Fehmi Cirak et al.

We introduce a coupled finite and boundary element formulation for acoustic scattering analysis over thin shell structures. A triangular Loop subdivision surface discretisation is used for both geometry and analysis fields. The Kirchhoff-Love shell equation is discretised with the finite element method and the Helmholtz equation for the acoustic field with the boundary element method. The use of the boundary element formulation allows the elegant handling of infinite domains and precludes the need for volumetric meshing. In the present work the subdivision control meshes for the shell displacements and the acoustic pressures have the same resolution. The corresponding smooth subdivision basis functions have the $C^1$ continuity property required for the Kirchhoff-Love formulation and are highly efficient for the acoustic field computations. We validate the proposed isogeometric formulation through a closed-form solution of acoustic scattering over a thin shell sphere. Furthermore, we demonstrate the ability of the proposed approach to handle complex geometries with arbitrary topology that provides an integrated isogeometric design and analysis workflow for coupled structural-acoustic analysis of shells.

Kosala Bandara, Thomas Rüberg, Fehmi Cirak

We develop a new optimisation technique that combines multiresolution subdivision surfaces for boundary description with immersed finite elements for the discretisation of the primal and adjoint problems of optimisation. Similar to wavelets multiresolution surfaces represent the domain boundary using a coarse control mesh and a sequence of detail vectors. Based on the multiresolution decomposition efficient and fast algorithms are available for reconstructing control meshes of varying fineness. During shape optimisation the vertex coordinates of control meshes are updated using the computed shape gradient information. By virtue of the multiresolution editing semantics, updating the coarse control mesh vertex coordinates leads to large-scale geometry changes and, conversely, updating the fine control mesh coordinates leads to small-scale geometry changes. In our computations we start by optimising the coarsest control mesh and refine it each time the cost function reaches a minimum. This approach effectively prevents the appearance of non-physical boundary geometry oscillations and control mesh pathologies, like inverted elements. Independent of the fineness of the control mesh used for optimisation, on the immersed finite element grid the domain boundary is always represented with a relatively fine control mesh of fixed resolution. With the immersed finite element method there is no need to maintain an analysis suitable domain mesh. In some of the presented two- and three-dimensional elasticity examples the topology derivative is used for creating new holes inside the domain.