René R. Hiemstra

Maarten M. Hodzelmans, René R. Hiemstra, Joris J. C. Remmers et al.

For higher-order discretizations of explicit dynamics problems, Isogeometric Analysis (IGA) has several favorable properties as compared to classical Finite Element Analysis (FEA). While FEA produces spurious modes at orders beyond linear, this is not the case for IGA. Consequently, fewer degrees of freedom are required for comparable accuracy, larger timesteps can be taken, and the method is more robust for nonlinear problems. If outlier modes are removed, the timestep even becomes virtually independent of the order. In this paper, we investigate how these advantages apply to the poroelastic continuum model. We consider both a primal formulation, wherein our variables are the displacement of the matrix material, the fluid displacement, and the pressure, as well as a reduced form wherein the pressure is eliminated. For our discretizations, we employ divergence-conforming spline spaces. Conforming spline spaces for the fluid displacement ensure inf-sup stability for the primal form, as well as a correct null space in the reduced form. Furthermore, we prove and demonstrate that the two formulations coincide when both displacements are discretized with conforming spline spaces. Through spectral analysis, we find that the aforementioned benefits of IGA do carry over directly to the context of poroelasticity. In 1D, we split the discrete spectrum into fast and slow waves. When normalized against an analytical solution, each of these sub-spectra closely resembles results known in elasticity. Consequently, when poroelasticity is discretized with outlier-free IGA, the timestep is essentially independent of the order. We show this timestep scaling in 2D as well.

René R. Hiemstra, Torsten Keßler, Michael R. A. Abdelmalik

We reduce the eight-dimensional weak form of the bilinear Boltzmann collision operator to a five-dimensional kinematic core by rigidly rotating the laboratory frame to align with the colliding pair and integrating over the $\mathrm{SO}(3)$ rotation group. This reduction yields an exact Wigner--Eckart factorization within a spectral Galerkin framework of associated Laguerre polynomials and spherical harmonics. The decomposition decouples the angular geometry from the scattering physics. The former, represented by Clebsch--Gordan coefficients, is evaluated exactly, while the latter is evaluated to machine precision by a spectrally convergent singular quadrature strategy. By explicitly zeroing specific entries, the macroscopic collision invariants are embedded without approximation. Cache-optimized contractions deliver up to a 37-fold single-core speedup and a 1000-fold memory reduction over standard dense Cartesian formulations. The approach is validated against analytical solutions for Maxwell molecules and infinite-order Chapman--Enskog viscosity coefficients for hard spheres.