Related papers: A fast, matrix-based method to perform omnidirecti…
The computational cost of the boundary-condition-enforced immersed boundary method (IBM) increases in the order of $\mathcal{O}(N^2)$ as the number of Lagrangian points, $N$, increases. This is due to the time-consuming calculation of the…
A boundary thickening-based direct forcing (BTDF) immersed boundary (IB) method is proposed for fully resolved simulation of incompressible viscous flows laden with finite size particles. By slightly thickening the boundary thickness, the…
We present a flexible discretization technique for computational models of thin tubular networks embedded in a bulk domain, for example a porous medium. These systems occur in the simulation of fluid flow in vascularized biological tissue,…
We describe a new, adaptive solver for the two-dimensional Poisson equation in complicated geometries. Using classical potential theory, we represent the solution as the sum of a volume potential and a double layer potential. Rather than…
Optimal exploitation of supercomputing resources for the evaluation of electrostatic forces remains a challenge in molecular dynamics simulations of very large systems. The most efficient methods are currently based on particle-mesh Ewald…
This work outlines a new multi-physics-compatible immersed rigid body method for Eulerian finite-volume simulations. To achieve this, rigid bodies are represented as a diffuse scalar field and an interface seeding method is employed to…
We present a fast, direct and adaptive Poisson solver for complex two-dimensional geometries based on potential theory and fast multipole acceleration. More precisely, the solver relies on the standard decomposition of the solution as the…
Two methods for solid body representation in flow simulations available in the Pencil Code are the immersed boundary method and overset grids. These methods are quite different in terms of computational cost, flexibility and numerical…
We propose a novel framework for simulating ink as a particle-laden flow using particle flow maps. Our method addresses the limitations of existing flow-map techniques, which struggle with dissipative forces like viscosity and drag, thereby…
We present a methodology for simulating three-dimensional flow of incompressible viscoplastic fluids modelled by generalised Newtonian rheological equations. It is implemented in a highly efficient framework for massively parallelisable…
We provide an alternative Fourier analysis for multigrid applied to the Poisson problem in 1D, based on explicit derivation of spectra of the iteration matrix. The new Fourier analysis has advantages over the existing one. It is easy to…
Measuring stress fields in fluids and soft materials is crucial in various fields such as mechanical engineering, medicine, and bioengineering. However, conventional methods that calculate stress fields from velocity fields struggle to…
Many physical systems of interest involve the close interaction of a flow in a domain with complex, time-varying boundaries. Treatment of boundaries of this nature is cumbersome due to the difficulty in explicitly tracking boundaries that…
The evolution of interfaces is intrinsic to many physical processes ranging from cavitation in fluids to recrystallization in solids. Computational modeling of interface motion entails a number of challenges, many of which are related to…
We extend the shifted boundary method (SBM) to the simulation of incompressible fluid flow using immersed octree meshes. Previous work on SBM for fluid flow primarily utilized two- or three-dimensional unstructured tetrahedral grids.…
Open-vocabulary panoptic reconstruction is essential for advanced robotics perception and simulation. However, existing methods based on 3D Gaussian Splatting (3DGS) often struggle to simultaneously achieve geometric accuracy, coherent…
In order to run Computational Fluid Dynamics (CFD) codes on large scale infrastructures, parallel computing has to be used because of the computational intensive nature of the problems. In this paper we investigate the ADAPT platform where…
We present a method that gives highly accurate electrostatic potentials for systems where we have periodic boundary conditions in two spatial directions but free boundary conditions in the third direction. These boundary conditions are…
Statistically correcting measured cross sections for detector effects is an important step across many applications. In particle physics, this inverse problem is known as unfolding. In cases with complex instruments, the distortions they…
Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems, whose basic concept is to embed physical laws to constrain/inform neural networks, with the need of less data for training…