A volume penalization approach to simulate magnetohydrodynamic (MHD) flows in confined domains is presented. Here the incompressible visco-resistive MHD equations are solved using parallel pseudo-spectral solvers in Cartesian geometries. The volume penalization technique is an immersed boundary method which is characterized by a high flexibility for the geometry of the considered flow. In the p...

We present new algorithms for the parallelization of Eulerian–Lagrangian interaction operations in immersed boundary method. Our rely on two well-studied parallel primitives: key-value sort and segmented reduce. The use these primitives allows us to implement our both graphics processing units (GPUs) other shared-memory architectures. strong weak scaling tests problems involving scattered point...

A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immersed interfaces to be discretized with a good accuracy on a compact stencil. Auxiliary unknowns are created at existing grid locations to increase the degrees of freedom of the initia...

In this study, the immersed boundary-thermal lattice Boltzmann method has been used to simulate non-Newtonian fluid flow over a heated circular cylinder. The direct-forcing algorithm has been employed to couple the off-lattice obstacles and on-lattice fluid nodes. To investigate the effect of boundary sharpness, two different diffuse interface schemes are considered to interpolate the velocity ...

The immersed boundary method is a numerical approach for simulating elastic structures which interact with a fluid flow. In many physical systems thermal fluctuations become significant at small scales and play a fundamental role. In this paper stochastic numerical methods are developed which extend the immersed boundary approach to account for thermal fluctuations by including appropriate stoc...