Related papers: GPU accelerated fast multipole boundary element me…
Three algorithm are proposed to evaluate volume potentials that arise in boundary element methods for elliptic PDEs. The approach is to apply a modified fast multipole method for a boundary concentrated volume mesh. If $h$ is the meshwidth…
We apply a novel optimization scheme from the image processing and machine learning areas, a fast Primal-Dual method, to achieve controllable and realistic fluid simulations. While our method is generally applicable to many problems in…
A discrete unified gas kinetic scheme (DUGKS) coupled with the immersed boundary (IB) method is developed to perform interface-resolved simulation of particle-laden flows. The present method (IB-DUGKS) preserves the respective advantages of…
Our objective is to stabilise and accelerate the time-domain boundary element method (TDBEM) for the three-dimensional wave equation. To overcome the potential time instability, we considered using the Burton--Miller-type boundary integral…
High-fidelity simulations of unsteady fluid flow are now possible with advancements in high-performance computing hardware and software frameworks. Since computational fluid dynamics (CFD) computations are dominated by linear algebraic…
The optimal control of a globally unstable two-dimensional separated boundary layer over a bump is considered using augmented Lagrangian optimization procedures. The present strategy allows of controlling the flow from a fully developed…
An energy stable finite element scheme within arbitrary Lagrangian Eulerian (ALE) framework is derived for simulating the dynamics of millimetric droplets in contact with solid surfaces. Supporting surfaces considered may exhibit…
Optimally convergent (with respect to the regularity) quadratic finite element method for two dimensional obstacle problem on simplicial meshes is studied in (Brezzi, Hager, Raviart, Numer. Math, 28:431--443, 1977). There was no analogue of…
We use existing 3D Discrete Element simulations of simple shear flows of spheres to evaluate the radial distribution function at contact that enables kinetic theory to correctly predict the pressure and the shear stress, for different…
We present teraflop-scale calculations of biomolecular electrostatics enabled by the combination of algorithmic and hardware acceleration. The algorithmic acceleration is achieved with the fast multipole method (FMM) in conjunction with a…
Deformable elastic bodies in viscous and viscoelastic media constitute a large portion of synthetic and biological complex fluids. We present a parallelized 3D-simulation methodology which fully resolves the momentum balance in the solid…
We present a new method to estimate terminal bubble velocities during first-order phase transitions in a plasma in local equilibrium. The method relies on calculating the extrema of a modified potential function for the scalar field…
We present a highly general implementation of fast multipole methods on graphics processing units (GPUs). Our two-dimensional double precision code features an asymmetric type of adaptive space discretization leading to a particularly…
We present a scalable dissipative particle dynamics simulation code, fully implemented on the Graphics Processing Units (GPUs) using a hybrid CUDA/MPI programming model, which achieves 10-30 times speedup on a single GPU over 16 CPU cores…
The multi-direct-forcing immersed boundary method allows for a small velocity error of the no-slip condition in moving-particle problems but suffers from numerical instability if simulation parameters are not carefully chosen. This study…
This paper introduces open-source computational fluid dynamics software named open computational fluid dynamic code for scientific computation with graphics processing unit (GPU) system (OpenCFD-SCU), developed by the authors for direct…
This work focuses on model preparation for electrostatic simulations of CAD designs to realize a rapid virtual prototyping concept. We present a boundary element method (BEM) allowing discontinuous fields between surfaces. The corresponding…
Mesoscopic simulations of hydrocarbon flow in source shales are challenging, in part due to the heterogeneous shale pores with sizes ranging from a few nanometers to a few micrometers. Additionally, the sub-continuum fluid-fluid and…
The Lagrangian vortex method offers an alternative numerical approach for direct numerical simulation of turbulence. The fact that it uses the fast multipole method (FMM)--a hierarchical algorithm for N-body problems with highly scalable…
Chaotic free surface flows are challenging problems to simulate numerically, mainly due to the significant changes in geometry and frequent topological changes. Methods that track the evolution of the fluid in a Lagrangian formulation are a…