A stable partitioned FSI algorithm for rigid bodies and incompressible flow in three dimensions

Banks, Jeffrey W.; Henshaw, William D.; Schwendeman, Donald W.; Tang, Qi

doi:10.1016/j.jcp.2018.06.072

Cited by 31 publications

(20 citation statements)

References 59 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The CBCs along the boundaries with s " r and s " τ are similar to those given in (12) and (13), respectively.…”

Section: The Wave Equation In Two Dimensionssupporting

confidence: 57%

“…For the Neumann conditions, we use second-order accurate approximations of CBC b r1s and CBC τ r1s. For example, the approximation of CBC b r1s in (13) given by…”

Section: The Wave Equation In Two Dimensionsmentioning

confidence: 99%

“…CBCs are also useful for problems involving material interfaces, such as conjugate heat transfer [10] and electromagnetics [5,7]. In recent work, we have developed Added-Mass Partitioned (AMP) schemes for a wide range of fluid-structure interaction (FSI) problems, including schemes for incompressible flows coupled to rigid bodies [11][12][13] and elastic solids [14,15]. These strongly-partitioned schemes incorporate AMP interface conditions derived using CBCs and the physical matching conditions at fluid-solid interfaces in order to overcome added-mass instabilities that can occur for the case of light bodies [16,17].…”

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Local Compatibility Boundary Conditions for High-Order Accurate Finite-Difference Approximations of PDEs

Al¹,

Banks²,

Henshaw³

et al. 2021

Preprint

View full text Add to dashboard Cite

We describe a new approach to derive numerical approximations of boundary conditions for highorder accurate finite-difference approximations. The approach, called the Local Compatibility Boundary Condition (LCBC) method, uses boundary conditions and compatibility boundary conditions derived from the governing equations, as well as interior and boundary grid values, to construct a local polynomial, whose degree matches the order of accuracy of the interior scheme, centered at each boundary point. The local polynomial is then used to derive a discrete formula for each ghost point in terms of the data. This approach leads to centered approximations that are generally more accurate and stable than one-sided approximations. Moreover, the stencil approximations are local since they do not couple to neighboring ghost-point values which can occur with traditional compatibility conditions. The local polynomial is derived using continuous operators and derivatives which enables the automatic construction of stencil approximations at different orders of accuracy. The LCBC method is developed here for problems governed by second-order partial differential equations, and it is verified for a wide range of sample problems, both time-dependent and time-independent, in two space dimensions and for schemes up to sixth-order accuracy.

show abstract

“…The CBCs along the boundaries with s " r and s " τ are similar to those given in (12) and (13), respectively.…”

Section: The Wave Equation In Two Dimensionssupporting

confidence: 57%

“…For the Neumann conditions, we use second-order accurate approximations of CBC b r1s and CBC τ r1s. For example, the approximation of CBC b r1s in (13) given by…”

Section: The Wave Equation In Two Dimensionsmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Local Compatibility Boundary Conditions for High-Order Accurate Finite-Difference Approximations of PDEs

Al¹,

Banks²,

Henshaw³

et al. 2021

Preprint

View full text Add to dashboard Cite

show abstract

“…Given that the cross-sectional area of the reservoir is much larger than that of the pipeline, the water level in the tank can be treated as a constant [30] P = P 0…”

Section: Outlet and Inlet Boundarymentioning

confidence: 99%

Fluid-Structure Interaction Response of a Water Conveyance System with a Surge Chamber during Water Hammer

Guo

Zhou

et al. 2020

Water

View full text Add to dashboard Cite

Fluid-structure interaction (FSI) is a frequent and unstable inherent phenomenon in water conveyance systems. Especially in a system with a surge chamber, valve closing and the subsequent water level oscillation in the surge chamber are the excitation source of the hydraulic transient process. Water-hammer-induced FSI has not been considered in preceding research, and the results without FSI justify further investigations. In this study, an FSI eight-equation model is presented to capture its influence. Both the elbow pipe and surge chamber are treated as boundary conditions, and solved using the finite volume method (FVM). After verifying the feasibility of using FVM to solve FSI, friction, Poisson, and junction couplings are discussed in detail to separately reveal the influence of a surge chamber, tow elbows, and a valve on FSI. Results indicated that the major mechanisms of coupling are junction coupling and Poisson coupling. The former occurs in the surge chamber and elbows. Meanwhile, a stronger pressure pulsation is produced at the valve, resulting in a more complex FSI response in the water conveyance system. Poisson coupling and junction coupling are the main factors contributing to a large amount of local transilience emerging on the dynamic pressure curves. Moreover, frictional coupling leads to the lower amplitudes of transilience. These results indicate that the transilience is induced by the water hammer-structure interaction and plays important roles in the orifice optimization in the surge chamber.Keywords: fluid-structure interaction; pipe flow; finite volume method; water hammer; transient flow IntroductionPipes conveying fluid are prevalent in many fields, including marine, civil engineering, nuclear power industries, petroleum, and water conservancy systems in daily life [1]. As an inherent phenomenon, fluid-structure interaction (FSI) always occurs in water conveyance pipelines [2][3][4]. Because of the existence of FSI, those pipes with few supports or thin walls show poor robustness, where the FSI of pipes is enhanced [5]. Therefore, the FSI responses must be considered when analyzing the characteristics of water conveyance systems [6]. The types of coupling that occur between the pipeline and fluid mainly include friction coupling, Poisson coupling, and junction coupling, among which the former two coupling take place throughout the whole pipeline. However, the last one only happens locally in pipes, including in elbows, branches, valves, boundaries, and variable cross sections [7], where the system coupling is much stronger [8]. Amongst these three coupling forms, friction coupling has the weakest response and shows changes in its magnitude over a long period [9]. Meanwhile, as the most important coupling form, junction coupling produces a pressure head larger Water 2020, 12, 1025 2 of 18 than the classical water hammer, and greatly depends on the robustness of the system [10]. Poisson coupling and friction coupling can greatly influence junction coupling, and, in turn, junction coupling ...

show abstract

“…38,39 To gain a fully stable solution at arbitrary density ratios, much more costly implicit schemes were proposed. [40][41][42] The reason for the numerical instabilities in explicit coupling schemes is attributed to the particles' mass being exceeded by its so-called added mass, that is, the fluid attached to the moving particle. In an explicit coupling scenario, this mass cannot be properly accounted for, which results in excessive accelerations on startup that lead to oscillations.…”

Section: Introductionmentioning

confidence: 99%

Coupling fully resolved light particles with the lattice Boltzmann method on adaptively refined grids

Werner

Rettinger

Rüde

2021

Numerical Methods in Fluids

View full text Add to dashboard Cite

The simulation of geometrically resolved rigid particles in a fluid relies on coupling algorithms to transfer momentum both ways between the particles and the fluid. In this article, the fluid flow is modeled with a parallel lattice Boltzmann method using adaptive grid refinement to improve numerical efficiency. The coupling with the particles is realized with the momentum exchange method.When implemented in plain form, instabilities may arise in the coupling when the particles are lighter than the fluid. The algorithm can be stabilized with a virtual mass correction specifically developed for the lattice Boltzmann method.The method is analyzed for a wide set of physically relevant regimes, varying independently the body-to-fluid density ratio and the relative magnitude of inertial and viscous effects. These studies of a single rising particle exhibit periodic regimes of particle motion as well as chaotic behavior, as previously reported in the literature. The new coupled lattice Boltzmann method is compared with available experimental and numerical results. This serves to validate the presented method and additionally it leads to new physical insight for some of the parameter settings.

show abstract

A stable partitioned FSI algorithm for rigid bodies and incompressible flow in three dimensions

Cited by 31 publications

References 59 publications

Local Compatibility Boundary Conditions for High-Order Accurate Finite-Difference Approximations of PDEs

Local Compatibility Boundary Conditions for High-Order Accurate Finite-Difference Approximations of PDEs

Fluid-Structure Interaction Response of a Water Conveyance System with a Surge Chamber during Water Hammer

Coupling fully resolved light particles with the lattice Boltzmann method on adaptively refined grids

Contact Info

Product

Resources

About