Experimental mixing parameterization due to multiphase fluid–structure interactions

Experimental estimates of the rate at which energy is extracted from the barotropic tide at corrugated sloping topography are presented. To this end, a series of experimental simulations of the process of generation of internal tides (i.e., internal waves of the tidal frequency) over a corrugated slope in stratified fluid were performed. An oceanic interior is modeled as four-layer stratified fluid. The main focus in these studies was to obtain the relation for the potential energy available for mixing as a function of a slope of a corrugated slope. The available potential energy available for partial mixing to the topographic slope was compared with the available potential energy sufficient for complete mixing of the four layers. The experimental data were compared with the analytic results of a linear theory and found in agreement with a recent theoretically predicted scenario showing that the dominant contribution of the energy distribution in internal wave field resides in waves of the lowest allowed frequency.     

An ALE approach to mechano‐chemical processes in fluid–structure interactions

Mathematical modeling and simulation of fluid–structure interaction problems are in the focus of research already for a longer period. However, taking into account also chemical reactions, leading to structural changes, including changes of mechanical properties of the solid phase, is rather new but for many applications is highly important area. This paper formulates a model system for reactive flow and transport in a vessel system, the penetration of chemical substances into the solid wall. Inside the wall, reactions take place that lead to changes of volume and of the mechanical properties of the wall. Numerical algorithms are developed and used to simulate the dynamics of such a mechano‐chemical fluid–structure interaction problem. As a proof of concept scenario, plaque formation in blood vessels is chosen. The arbitrary Lagrangian Eulerian approach (ALE) is chosen to solve the systems numerically. Temporal discretization of the fully coupled monolithic model is accomplished by backward Euler scheme and spatial discretization by stabilized finite elements. The numerical approach is verified by numerical tests, and effective methods to maintain mesh qualities under large deformations are described. For realistic system parameters, the simulations show that the plaque formation in blood vessel is a long‐time effect. The time scale of the formation is in the simulation of comparable order as in reality. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order implicit Runge–Kutta time integrators for fluid‐structure interactions

This paper presents an approach to develop high‐order, temporally accurate, finite element approximations of fluid‐structure interaction (FSI) problems. The proposed numerical method uses an implicit monolithic formulation in which the same implicit Runge–Kutta (IRK) temporal integrator is used for the incompressible flow, the structural equations undergoing large displacements, and the coupling terms at the fluid‐solid interface. In this context of stiff interaction problems, the fully implicit one‐step approach presented is an original alternative to traditional multistep or explicit one‐step finite element approaches. The numerical scheme takes advantage of an arbitrary Lagrangian–Eulerian formulation of the equations designed to satisfy the geometric conservation law and to guarantee that the high‐order temporal accuracy of the IRK time integrators observed on fixed meshes is preserved on arbitrary Lagrangian–Eulerian deforming meshes. A thorough review of the literature reveals that in most previous works, high‐order time accuracy (higher than second order) is seldom achieved for FSI problems. We present thorough time‐step refinement studies for a rigid oscillating‐airfoil on deforming meshes to confirm the time accuracy on the extracted aerodynamics reactions of IRK time integrators up to fifth order. Efficiency of the proposed approach is then tested on a stiff FSI problem of flow‐induced vibrations of a flexible strip. The time‐step refinement studies indicate the following: stability of the proposed approach is always observed even with large time step and spurious oscillations on the structure are avoided without added damping. While higher order IRK schemes require more memory than classical schemes (implicit Euler), they are faster for a given level of temporal accuracy in two dimensions. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of flow problems with moving mechanical components,fluid–structure interactions and two-fluid interfaces

The application of a stabilized space–time finite element formulation to problems involving fluid–structure interactions and two-fluid interfaces is discussed. Two sample problems are presented and the method is validated by comparison with a test problem. © 1997 John Wiley & Sons, Ltd.     

Stability analysis of numerical interface conditions in fluid–structure thermal analysis

This paper analyses the numerical stability of coupling procedures in modelling the thermal diffusion in a solid and a fluid with continuity of temperature and heat flux at the interface. A simple one-dimensional model is employed with uniform material properties and grid density in each domain. A number of different explicit and implicit algorithms are considered for both the interior equations and the boundary conditions. The analysis shows that in general these are stable provided that Dirichlet boundary conditions are imposed on the fluid and Neumann boundary conditions are imposed on the solid; in each case the imposed values are obtained from the other domains. © 1997 John Wiley & Sons, Ltd.     

Strongly coupled simulation of fluid–structure interaction in a Francis hydroturbine

This work simulates a complex fluid flow in fluid–structure interaction (FSI). The flow under consideration is governed by Navier–Stokes equations for incompressible viscous fluids and modeled with the finite volume method. Large eddy simulation is used to simulate the unsteady turbulent flow. The structure is represented by a finite element formulation. The present work introduces a strongly coupled partitioned approach that is applied to complex flow in fluid machinery. In this approach, the fluid and structure equations are solved separately using different solvers, but are implicitly coupled into one single module based on sensitivity analysis of the important displacement and stress modes. The applied modes and their responses are used to build up a reduced‐order model. The proposed model is used to predict the unsteady flow fields of a 3D complete passage, involving in stay, guide vanes, and runner blades, for a Francis hydro turbine and FSI is considered. The computational results show that a fairly good convergence solution is achieved by using the reduced‐order model that is based on only a few displacement and stress modes, which largely reduces the computational cost, compared with traditional approaches. At the same time, a comparison of the numerical results of the model with available experimental data validates the methodology and assesses its accuracy. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of fluid–structure interaction by an arbitrary Lagrangian–Eulerian finite element formulation

In this paper, the interaction fluid–rigid body is analysed by a finite element procedure that incorporates the arbitrary Lagrangian–Eulerian (ALE) method into a well‐known two‐step projection scheme. The flow is assumed to be two‐dimensional, incompressible and viscous, with no turbulence models being included. The flow past a circular cylinder at ℛℯ=200 is first analysed, for fixed and oscillating conditions. The dependence of lock‐in upon the shift between the mechanical and the Strouhal frequencies, for a given amplitude of forced vibration, is illustrated. The aerodynamic forces and the wake geometry are compared for locked‐in conditions with different driving frequencies. The behaviour of a rectangular cylinder (B/D=4) at ℛℯ=500 (based on height D) is also analysed. The flutter derivatives associated with aerodynamic damping (H₁* and A₂* in Scanlan's notation) are evaluated by the free oscillation method for several values of reduced flow speed above the Strouhal one (namely for 3≤U*≤8). Torsional flutter was attained at U*≥5, with all the other situations showing stable characteristics. Copyright © 1999 John Wiley & Sons, Ltd.     

A fast method for solving fluid–structure interaction problems numerically

The paper presents a semi‐implicit algorithm for solving an unsteady fluid–structure interaction problem. The algorithm for solving numerically the fluid–structure interaction problems was obtained by combining the backward Euler scheme with a semi‐implicit treatment of the convection term for the Navier–Stokes equations and an implicit centered scheme for the structure equations. The structure is governed either by the linear elasticity or by the non‐linear St Venant–Kirchhoff elasticity models. At each time step, the position of the interface is predicted in an explicit way. Then, an optimization problem must be solved, such that the continuity of the velocity as well as the continuity of the stress hold at the interface. During the Broyden, Fletcher, Goldforb, Shano (BFGS) iterations for solving the optimization problem, the fluid mesh does not move, which reduces the computational effort. The term ‘semi‐implicit’ used for the fully algorithm means that the interface position is computed explicitly, while the displacement of the structure, velocity and the pressure of the fluid are computed implicitly. Numerical results are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Fully coupled fluid–structural interactions using an efficient high resolution upwind scheme

The E-CUSP upwind scheme recently developed by Zha and Hu is extended to the moving grid system and is applied to calculate the flow-induced vibration based on a fully coupled fluid–structure interaction methodology. The scheme is used to calculate the flow-induced vibration of an elastically mounted cylinder, forced pitching airfoil, and an elastically mounted airfoil. The numerical results are compared with experimental data and have demonstrated the accuracy, efficiency and robustness of the new E-CUSP scheme for solving flow-induced vibration problems with moving mesh systems.     

Modelling of pipeline integrity taking into account the fluid–structure interaction

Two mechanical models have been presented in this paper for structural failure prediction of piping systems conveying liquids subjected to pressure transients. One model takes into account the axial fluid–structure interaction (FSI) phenomenon between fluid and pipe motion, whereas the other refers to an extension of the well-known waterhammer formulation. Both models are described by a system of non-linear hyperbolic equations which are solved by using a numerical procedure based upon the operator splitting technique and Glimm's scheme. To implement Glimm's method, it is presented the solution of a 4×4 Riemann problem with discontinuous coefficients. Numerical predictions of both models are presented and compared, so that the influence of the FSI term on the failure analysis is focused on. © 1998 John Wiley & Sons, Ltd.     

Computational algorithms for tracking dynamic fluid–structure interfaces in embedded boundary methods

A robust, accurate, and computationally efficient interface tracking algorithm is a key component of an embedded computational framework for the solution of fluid–structure interaction problems with complex and deformable geometries. To a large extent, the design of such an algorithm has focused on the case of a closed embedded interface and a Cartesian computational fluid dynamics grid. Here, two robust and efficient interface tracking computational algorithms capable of operating on structured as well as unstructured three‐dimensional computational fluid dynamics grids are presented. The first one is based on a projection approach, whereas the second one is based on a collision approach. The first algorithm is faster. However, it is restricted to closed interfaces and resolved enclosed volumes. The second algorithm is therefore slower. However, it can handle open shell surfaces and underresolved enclosed volumes. Both computational algorithms exploit the bounding box hierarchy technique and its parallel distributed implementation to efficiently store and retrieve the elements of the discretized embedded interface. They are illustrated, and their respective performances are assessed and contrasted, with the solution of three‐dimensional, nonlinear, dynamic fluid–structure interaction problems pertaining to aeroelastic and underwater implosion applications. Copyright © 2012 John Wiley & Sons, Ltd.     

A fictitious domain/mortar element method for fluid–structure interaction

A new method for the computational analysis of fluid–structure interaction of a Newtonian fluid with slender bodies is developed. It combines ideas of the fictitious domain and the mortar element method by imposing continuity of the velocity field along an interface by means of Lagrange multipliers. The key advantage of the method is that it circumvents the need for complicated mesh movement strategies common in arbitrary Lagrangian–Eulerian (ALE) methods, usually used for this purpose. Copyright © 2001 John Wiley & Sons, Ltd.     

Injection/suction boundary conditions for fluid–structure interaction simulations in incompressible flow

This paper presents the analysis of injection/suction boundary conditions in the context of the fluid–structure interactions simulation of the incompressible turbulent flow. First, the equations used in the modelling of the fluid and the structure are presented, as well as the numerical methods used in the corresponding solvers. Injection/suction boundary conditions are then presented with details of different implementation alternatives. Arbitrary Lagrangian–Eulerian (ALE) approach was also implemented in order to test the injection/suction boundary conditions. Numerical tests are performed where injection/suction boundary conditions are compared to ALE simulations. These tests include forced movement of the structure and two‐degrees‐of‐freedom structure model simulations. Copyright © 2002 John Wiley & Sons, Ltd.     

Nonlinear acceleration of coupled fluid–structure transient thermal problems by Anderson mixing

Conjugate heat‐transfer problems are typically solved using partitioned methods where fluid and solid subdomains are evaluated separately by dedicated solvers coupled through a common boundary. Strongly coupled schemes for transient analysis require fluid and solid problems to be solved many times each time step until convergence to a steady state. In many practical situations, a fairly simple and frequently employed fixed‐point iteration process is rather ineffective; it leads to a large number of iterations per time step and consequently to long simulation times. In this article, Anderson mixing is proposed as a fixed‐point convergence acceleration technique to reduce computational cost of thermal coupled fluid–solid problems. A number of other recently published methods with applications to similar fluid–structure interaction problems are also reviewed and analyzed. Numerical experiments are presented to illustrate relative performance of these methods on a test problem of rotating pre‐swirl cavity air flow interacting with a turbine disk. It is observed that performance of Anderson mixing method is superior to that of other algorithms in terms of total iteration counts. Additional computational savings are demonstrated by reusing information from previously solved time steps. Copyright © All rights reserved 2012 Rolls‐Royce plc.     

Blood flow dynamics and fluid–structure interaction in patient‐specific bifurcating cerebral aneurysms

Hemodynamics plays an important role in the progression and rupture of cerebral aneurysms. The current work describes the blood flow dynamics and fluid–structure interaction in seven patient‐specific models of bifurcating cerebral aneurysms located in the anterior and posterior circulation regions of the circle of Willis. The models were obtained from 3D rotational angiography image data, and blood flow dynamics and fluid–structure interaction were studied under physiologically representative waveform of inflow. The arterial wall was assumed to be elastic, isotropic and homogeneous. The flow was assumed to be laminar, non‐Newtonian and incompressible. In one case, the effects of different model suppositions and boundary conditions were reported in detail. The fully coupled fluid and structure models were solved with the finite elements package ADINA. The vortex structure, pressure, wall shear stress (WSS), effective stress and displacement of the aneurysm wall showed large variations, depending on the morphology of the artery, aneurysm size and position. The time‐averaged WSS, effective stress and displacement at the aneurysm fundus vary between 0.17 and 4.86 Pa, 4.35 and 170.2 kPa and 0.16 and 0.74 mm, respectively, for the seven patient‐specific models of bifurcating cerebral aneurysms. Copyright © 2008 John Wiley & Sons, Ltd.