A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis
Abstract
A stable partitioned algorithm is developed for fluidstructure interaction (FSI) problems involving viscous incompressible flow and rigid bodies. This addedmass partitioned (AMP) algorithm remains stable, without subiterations, for light and even zero mass rigid bodies when addedmass and viscous addeddamping effects are large. The scheme is based on a generalized Robin interface condition for the fluid pressure that includes terms involving the linear acceleration and angular acceleration of the rigid body. Added mass effects are handled in the Robin condition by inclusion of a boundary integral term that depends on the pressure. Addeddamping effects due to the viscous shear forces on the body are treated by inclusion of addeddamping tensors that are derived through a linearization of the integrals defining the force and torque. Addeddamping effects may be important at low Reynolds number, or, for example, in the case of a rotating cylinder or rotating sphere when the rotational moments of inertia are small. In this second part of a twopart series, the general formulation of the AMP scheme is presented including the form of the AMP interface conditions and addeddamping tensors for general geometries. A fully secondorder accurate implementation of the AMP scheme is developed in two dimensions basedmore »
 Authors:
 Rensselaer Polytechnic Inst., Troy, NY (United States)
 Publication Date:
 Research Org.:
 Rensselaer Polytechnic Inst., Troy, NY (United States)
 Sponsoring Org.:
 USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC21)
 OSTI Identifier:
 1353365
 Grant/Contract Number:
 AC5207NA27344
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Journal of Computational Physics
 Additional Journal Information:
 Journal Volume: 343; Journal Issue: C; Journal ID: ISSN 00219991
 Publisher:
 Elsevier
 Country of Publication:
 United States
 Language:
 English
 Subject:
 97 MATHEMATICS AND COMPUTING; fluidstructure interaction; moving overlapping; grids; incompressible NavierStokes; partitioned schemes; addedmass; addeddamping; rigid bodies
Citation Formats
Banks, J. W., Henshaw, W. D., Schwendeman, D. W., and Tang, Qi. A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis. United States: N. p., 2017.
Web. doi:10.1016/j.jcp.2017.01.015.
Banks, J. W., Henshaw, W. D., Schwendeman, D. W., & Tang, Qi. A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis. United States. doi:10.1016/j.jcp.2017.01.015.
Banks, J. W., Henshaw, W. D., Schwendeman, D. W., and Tang, Qi. Fri .
"A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis". United States.
doi:10.1016/j.jcp.2017.01.015. https://www.osti.gov/servlets/purl/1353365.
@article{osti_1353365,
title = {A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis},
author = {Banks, J. W. and Henshaw, W. D. and Schwendeman, D. W. and Tang, Qi},
abstractNote = {A stable partitioned algorithm is developed for fluidstructure interaction (FSI) problems involving viscous incompressible flow and rigid bodies. This addedmass partitioned (AMP) algorithm remains stable, without subiterations, for light and even zero mass rigid bodies when addedmass and viscous addeddamping effects are large. The scheme is based on a generalized Robin interface condition for the fluid pressure that includes terms involving the linear acceleration and angular acceleration of the rigid body. Added mass effects are handled in the Robin condition by inclusion of a boundary integral term that depends on the pressure. Addeddamping effects due to the viscous shear forces on the body are treated by inclusion of addeddamping tensors that are derived through a linearization of the integrals defining the force and torque. Addeddamping effects may be important at low Reynolds number, or, for example, in the case of a rotating cylinder or rotating sphere when the rotational moments of inertia are small. In this second part of a twopart series, the general formulation of the AMP scheme is presented including the form of the AMP interface conditions and addeddamping tensors for general geometries. A fully secondorder accurate implementation of the AMP scheme is developed in two dimensions based on a fractionalstep method for the incompressible NavierStokes equations using finite difference methods and overlapping grids to handle the moving geometry. Here, the numerical scheme is verified on a number of difficult benchmark problems.},
doi = {10.1016/j.jcp.2017.01.015},
journal = {Journal of Computational Physics},
number = C,
volume = 343,
place = {United States},
year = {Fri Jan 20 00:00:00 EST 2017},
month = {Fri Jan 20 00:00:00 EST 2017}
}
Web of Science

An addedmass partitioned (AMP) algorithm is described for solving fluid–structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully secondorder accurate and stable, without subtimestep iterations, even for very light structures when addedmass effects are strong. The fluid, governed by the incompressible Navier–Stokes equations, is solved in velocitypressure form using a fractionalstep method; large deformations are treated with a mixed EulerianLagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler–Bernoulli beam model, and these equations are solved in a Lagrangian frame usingmore »

A quadrature of the twodimensional equations of the hydrodynamics of an incompressible ideal liquid and a model of the flow around convex symmertic bodies with suction or injection
A class of solutions of the steadystate hydrodynamic equations of an incompressible liquid has been considered in the Euler approximation when the pressure field is described by an additive function of the form P(x;y) = F(x) + G(y). It turns out that this class of solutions includes potential flows such as wave motion of a liquid with a finite depth, a flow inside a right angle, and a hyersonic gas flow around a sphere in the approximation of constant density near the stagnation point. Among nonpotential flows, this class includes, in particular, a hyersonic flow around a cylinder. The resultsmore » 
A stable and accurate partitioned algorithm for conjugate heat transfer
We describe a new partitioned approach for solving conjugate heat transfer (CHT) problems where the governing temperature equations in different material domains are timestepped in a implicit manner, but where the interface coupling is explicit. The new approach, called the CHAMP scheme (Conjugate Heat transfer Advanced Multidomain Partitioned), is based on a discretization of the interface coupling conditions using a generalized Robin (mixed) condition. The weights in the Robin condition are determined from the optimization of a condition derived from a local stability analysis of the coupling scheme. The interface treatment combines ideas from optimizedSchwarz methods for domaindecomposition problems togethermore »Cited by 3