6,549 research outputs found
cuIBM -- A GPU-accelerated Immersed Boundary Method
A projection-based immersed boundary method is dominated by sparse linear
algebra routines. Using the open-source Cusp library, we observe a speedup
(with respect to a single CPU core) which reflects the constraints of a
bandwidth-dominated problem on the GPU. Nevertheless, GPUs offer the capacity
to solve large problems on commodity hardware. This work includes validation
and a convergence study of the GPU-accelerated IBM, and various optimizations.Comment: Extended paper post-conference, presented at the 23rd International
Conference on Parallel Computational Fluid Dynamics (http://www.parcfd.org),
ParCFD 2011, Barcelona (unpublished
Hybrid finite difference/finite element immersed boundary method
The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian
description of the structural deformations, stresses, and forces along with an Eulerian description of the
momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary
methods described immersed elastic structures using systems of flexible fibers, and even now, most
immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This
work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian
variables that facilitates independent spatial discretizations for the structure and background grid. This
approach employs a finite element discretization of the structure while retaining a finite difference scheme
for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively
contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases
in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors
that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes.
The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse
structural meshes with the immersed boundary method. This work also contrasts two different weak forms
of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations
facilitated by our coupling approach
Removing the Stiffness of Elastic Force from the Immersed Boundary Method for the 2D Stokes Equations
The Immersed Boundary method has evolved into one of the most useful
computational methods in studying fluid structure interaction. On the other
hand, the Immersed Boundary method is also known to suffer from a severe
timestep stability restriction when using an explicit time discretization. In
this paper, we propose several efficient semi-implicit schemes to remove this
stiffness from the Immersed Boundary method for the two-dimensional Stokes
flow. First, we obtain a novel unconditionally stable semi-implicit
discretization for the immersed boundary problem. Using this unconditionally
stable discretization as a building block, we derive several efficient
semi-implicit schemes for the immersed boundary problem by applying the Small
Scale Decomposition to this unconditionally stable discretization. Our
stability analysis and extensive numerical experiments show that our
semi-implicit schemes offer much better stability property than the explicit
scheme. Unlike other implicit or semi-implicit schemes proposed in the
literature, our semi-implicit schemes can be solved explicitly in the spectral
space. Thus the computational cost of our semi-implicit schemes is comparable
to that of an explicit scheme, but with a much better stability property.Comment: 40 pages with 8 figure
A three-dimensional immersed boundary method for non-Newtonian fluids
Fluid-structure-interaction (FSI) phenomenon is common in science and engineering. The fluid involved in an FSI problem may be non-Newtonian such as blood. A popular framework for FSI problems is Peskin's immersed boundary (IB) method. However, most of the IB formulations are based on Newtonian fluids. In this letter, we report an extension of the IB framework to FSI involving Oldroyd-B and FENE-P fluids in three dimensions using the lattice Boltzmann approach. The new method is tested on two FSI model problems. Numerical experiments show that the method is conditionally stable and convergent with the first order of accuracy
Solving Fluid Structure Interaction Problems with an Immersed Boundary Method
An immersed boundary method for the compressible Navier-Stokes equations can be used for moving boundary problems as well as fully coupled fluid-structure interaction is presented. The underlying Cartesian immersed boundary method of the Launch Ascent and Vehicle Aerodynamics (LAVA) framework, based on the locally stabilized immersed boundary method previously presented by the authors, is extended to account for unsteady boundary motion and coupled to linear and geometrically nonlinear structural finite element solvers. The approach is validated for moving boundary problems with prescribed body motion and fully coupled fluid structure interaction problems. Keywords: Immersed Boundary Method, Higher-Order Finite Difference Method, Fluid Structure Interaction
- …