Improved radial basis function fluid–structure coupling via efficient localized implementation

Previous work by the authors has developed a universal interpolation scheme, using radial basis functions (RBFs), which results in a unified formulation for robust fluid–structure interpolation and high‐quality mesh motion. The method has several significant advantages. Primarily, all volume mesh, structural mesh, and flow‐solver‐type dependence is removed entirely, as all operations are performed on totally arbitrary point clouds of any form. Hence, all connectivity requirements are removed from both the coupling and mesh motion problems. Furthermore, only matrix‐vector multiplications are required during unsteady simulation because dependence relations are computed once prior to any simulation and then remain constant. This property means that the method is both perfectly parallel and totally independent from the flow‐solver. However, the full method is expensive, since the dependence matrix between two sets of points is N × N. The fluid–structure coupling behaviour can also be influenced by parameters used in the interpolation. To alleviate these difficulties a more efficient form of the RBF fluid–structure coupling is presented, which also greatly reduces the interpolation parameter influence. A pointwise form of the partition of unity approach is developed that localizes the interpolation, with results presented for static aeroelastic simulations of the Brite‐Euram multi‐disciplinary optimization wing using a very fine mesh containing 58 000 surface points. It is shown that a 58 × reduction in data size is achieved, and equally importantly the interpolation has a much smaller influence on final aeroelastic results. Copyright © 2009 John Wiley & Sons, Ltd.     

Aeroacoustic source term computation based on radial basis functions

In low Mach number aeroacoustics, the known disparity of length scales makes it possible to apply well-suited simulation models using different meshes for flow and acoustics. The workflow of these hybrid methodologies include performing an unsteady flow simulation, computing the acoustic sources, and simulating the acoustic field. Therefore, hybrid methods seek for robust and flexible procedures, providing a conservative mesh to mesh interpolation of the sources while ensuring high computational efficiency. We propose a highly specialized radial basis function interpolation for the challenges during hybrid simulations. First, the computationally efficient local radial basis function interpolation in conjunction with a connectivity-based neighbor search technique is presented. Second, we discuss the computation of spatial derivatives based on radial basis functions. These derivatives are computed in a local-global approach, using a Gaussian kernel on local point stencils. Third, radial basis function interpolation and derivatives are used to compute complex aeroacoustic source terms. These ingredients are necessary to provide flexible source term calculations that robustly connect flow and acoustics. Finally, the capabilities of the presented approach are shown in a numerical experiment with a co-rotating vortex pair.     

Parallel efficient mesh motion using radial basis functions with application to multi‐bladed rotors

Radial basis functions are used to provide a solution to the problem of mesh motion for unsteady aerodynamic simulation. The method is independent of connectivity and produces high‐quality meshes, but is expensive for large meshes in its full form. Hence, the efficiency of the technique has been greatly improved here by reducing the number of surface points used to define deformations of the surface, and the minor error in position that this implies at other surface points is corrected with a simple decaying perturbation, thus splitting the method into a primary basis function method and a secondary local correction method. This means that the exact surface is retained, but the mesh motion is significantly faster, while splitting the motion into two stages allows both the methods to work on appropriate problems given their relative strengths. An example deformation for a 5×10⁶ cell helicopter rotor mesh with an exaggerated cyclic pitch motion shows excellent mesh quality, thus validating a scheme that is also simple, robust and readily parallelized. Copyright © 2009 John Wiley & Sons, Ltd.     

A radial basis function algorithm for simplified fluid‐structure data transfer

In loosely‐coupled aeroelastic computation, the aerodynamic and elastomechanical models are based on different grids and eventually a simplified structural model, like a shell or stick, is considered. CFD tools are applied over the aerodynamic grid, whereas CSM tools are over the structural one. The meshes are usually non‐conforming; thus, three‐dimensional and non‐intrusive interpolation procedures are necessary to transfer structural deformations to the aerodynamic surface grid and aerodynamic loads to the structure. For that purpose, an interpolator based on radial basis functions has been developed. This tool does not require grid connectivities and can be applied to any three‐dimensional data. We have developed two strategies to deal with two special cases of interest. The first is when the structure is represented by a beam model and there is not a proper structural mesh close to the aerodynamic surface. The second is when managing a full configuration aircraft and the problem has been split into some small blocks in such a way that there are aerodynamic nodes belonging to more than one block. Copyright © 2014 John Wiley & Sons, Ltd.     

The fixed‐mesh ALE approach applied to solid mechanics and fluid–structure interaction problems

In this paper we propose a method to solve Solid Mechanics and fluid–structure interaction problems using always a fixed background mesh for the spatial discretization. The main feature of the method is that it properly accounts for the advection of information as the domain boundary evolves. To achieve this, we use an Arbitrary Lagrangian–Eulerian (ALE) framework, the distinctive characteristic being that at each time step results are projected onto a fixed, background mesh. For solid mechanics problems subject to large strains, the fixed‐mesh (FM)‐ALE method avoids the element stretching found in fully Lagrangian approaches. For FSI problems, FM‐ALE allows for the use of a single background mesh to solve both the fluid and the structure. Copyright © 2009 John Wiley & Sons, Ltd.     

Immersed stress method for fluid–structure interaction using anisotropic mesh adaptation

This paper presents advancements toward a monolithic solution procedure and anisotropic mesh adaptation for the numerical solution of fluid–structure interaction with complex geometry. First, a new stabilized three‐field stress, velocity, and pressure finite element formulation is presented for modeling the interaction between the fluid (laminar or turbulent) and the rigid body. The presence of the structure will be taken into account by means of an extra stress in the Navier–Stokes equations. The system is solved using a finite element variational multiscale method. We combine this method with anisotropic mesh adaptation to ensure an accurate capturing of the discontinuities at the fluid–solid interface. We assess the behavior and accuracy of the proposed formulation in the simulation of 2D and 3D time‐dependent numerical examples such as the flow past a circular cylinder and turbulent flows behind an immersed helicopter in a forward flight. Copyright © 2013 John Wiley & Sons, Ltd.     

A monolithic strategy for fluid–structure interaction problems

In this work, we present a new monolithic strategy for solving fluid–structure interaction problems involving incompressible fluids, within the context of the finite element method. This strategy, similar to the continuum dynamics, conserves certain properties, and thus provides a rational basis for the design of the time‐stepping strategy; detailed proofs of the conservation of these properties are provided. The proposed algorithm works with displacement and velocity variables for the structure and fluid, respectively, and introduces no new variables to enforce velocity or traction continuity. Any existing structural dynamics algorithm can be used without change in the proposed method. Use of the exact tangent stiffness matrix ensures that the algorithm converges quadratically within each time step. An analytical solution is presented for one of the benchmark problems used in the literature, namely, the piston problem. A number of benchmark problems including problems involving free surfaces such as sloshing and the breaking dam problem are used to demonstrate the good performance of the proposed method. Copyright © 2010 John Wiley & Sons, Ltd.     

Immersed particle method for fluid–structure interaction

A method for treating fluid–structure interaction of fracturing structures under impulsive loads is described. The coupling method is simple and does not require any modifications when the structure fails and allows fluid to flow through openings between crack surfaces. Both the fluid and the structure are treated by meshfree methods. For the structure, a Kirchhoff–Love shell theory is adopted and the cracks are treated by introducing either discrete (cracking particle method) or continuous (partition of unity‐based method) discontinuities into the approximation. Coupling is realized by a master–slave scheme where the structure is slave to the fluid. The method is aimed at problems with high‐pressure and low‐velocity fluids, and is illustrated by the simulation of three problems involving fracturing cylindrical shells coupled with fluids. Copyright © 2009 John Wiley & Sons, Ltd.     

A numerical method for heat transfer problems using collocation and radial basis functions

Simple, mesh/grid free, numerical schemes for the solution of heat transfer problems are developed and validated. Unlike the mesh or grid-based methods, these schemes use well-distributed quasi-random collocation points and approximate the solution using radial basis functions. The schemes work in a similar fashion as finite differences but with random points instead of a regular grid system. This allows the computation of problems with complex-shaped boundaries in higher dimensions with no extra difficulty. © 1998 John Wiley & Sons, Ltd.     

Topology optimization for stationary fluid–structure interaction problems using a new monolithic formulation

This paper outlines a new procedure for topology optimization in the steady‐state fluid–structure interaction (FSI) problem. A review of current topology optimization methods highlights the difficulties in alternating between the two distinct sets of governing equations for fluid and structure dynamics (hereafter, the fluid and structural equations, respectively) and in imposing coupling boundary conditions between the separated fluid and solid domains. To overcome these difficulties, we propose an alternative monolithic procedure employing a unified domain rather than separated domains, which is not computationally efficient. In the proposed analysis procedure, the spatial differential operator of the fluid and structural equations for a deformed configuration is transformed into that for an undeformed configuration with the help of the deformation gradient tensor. For the coupling boundary conditions, the divergence of the pressure and the Darcy damping force are inserted into the solid and fluid equations, respectively. The proposed method is validated in several benchmark analysis problems. Topology optimization in the FSI problem is then made possible by interpolating Young's modulus, the fluid pressure of the modified solid equation, and the inverse permeability from the damping force with respect to the design variables. Copyright © 2009 John Wiley & Sons, Ltd.     

Truly monolithic algebraic multigrid for fluid–structure interaction

The coupling of flexible structures to incompressible fluids draws a lot of attention during the last decade. Many different solution schemes have been proposed. In this contribution, we concentrate on the strong coupling fluid–structure interaction by means of monolithic solution schemes. Therein, a Newton–Krylov method is applied to the monolithic set of nonlinear equations. Such schemes require good preconditioning to be efficient. We propose two preconditioners that apply algebraic multigrid techniques to the entire fluid–structure interaction system of equations. The first is based on a standard block Gauss–Seidel approach, where approximate inverses of the individual field blocks are based on a algebraic multigrid hierarchy tailored for the type of the underlying physical problem. The second is based on a monolithic coarsening scheme for the coupled system that makes use of prolongation and restriction projections constructed for the individual fields. The resulting nonsymmetric monolithic algebraic multigrid method therefore involves coupling of the fields on coarse approximations to the problem yielding significantly enhanced performance. Copyright © 2010 John Wiley & Sons, Ltd.     

A new staggered scheme for fluid–structure interaction

Staggered solution procedures represent the most elementary computational strategy for the simulation of fluid–structure interaction problems. They usually consist of a predictor followed by the separate execution of each subdomain solver. Although it is generally possible to maintain the desired order of accuracy of the time integration, it is difficult to guarantee the stability of the overall computation. In the context of large solid over fluid mass ratios, compressible flows and explicit subsolvers, substantial development has been carried out by Felippa, Park, Farhat, Löhner and others. In this work, a new staggered scheme is presented. It is shown that, for a linear model problem, the scheme is second‐order accurate and unconditionally stable. The dependency of the leading truncation error on the solid over fluid mass ratio is investigated. The strategy is applied to two‐dimensional and three‐dimensional fluid–structure interaction problems. It is shown that the conclusions derived from the investigation of the model problem apply. The new strategy extends the applicability of staggered schemes to problems involving relatively small solid over fluid mass ratios and incompressible fluid flow. It is suggested that the proposed scheme has the same range of applicability as the Dirichlet–Neumann or block Gauß–Seidel type strategies. Copyright © 2012 John Wiley & Sons, Ltd.     

A mesh free approach using radial basis functions and parallel domain decomposition for solving three‐dimensional diffusion equations

The analysis of transient heat conduction problems in large, complex computational domains is a problem of interest in many technological applications including electronic cooling, encapsulation using functionally graded composite materials, and cryogenics. In many of these applications, the domains may be multiply connected and contain moving boundaries making it desirable to consider meshless methods of analysis. The method of fundamental solutions along with a parallel domain decomposition method is developed for the solution of three‐dimensional parabolic differential equations. In the current approach, time is discretized using the generalized trapezoidal rule transforming the original parabolic partial differential equation into a sequence of non‐homogeneous modified Helmholtz equations. An approximate particular solution is derived using polyharmonic splines. Interfacial conditions between subdomains are satisfied using a Schwarz Neumann–Neumann iteration scheme. Outside of the first time step where zero initial flux is assumed, the initial estimates for the interfacial flux is given from the converged solution obtained during the previous time step. This significantly reduces the number of iterations required to meet the convergence criterion. The accuracy of the method of fundamental solutions approach is demonstrated through two benchmark problems. The parallel efficiency of the domain decomposition method is evaluated by considering cases with 8, 27, and 64 subdomains. Copyright 2004 © John Wiley & Sons, Ltd.     

A parallel multiselection greedy method for the radial basis function–based mesh deformation

Greedy algorithm has been widely adopted for the point selection procedure of radial basis function–based mesh deformation. However, in large deformation simulations with thousands of points selected, the greedy point selection will be too expensive and thus become a performance bottleneck. To improve the efficiency of the point selection procedure, a parallel multiselection greedy method has been developed in this paper. Multiple points are selected at each step to accelerate the convergence speed of the greedy algorithm. In addition, 2 strategies are presented to determine the specific selecting number. The parallelization of the greedy point selection is realized on the basis of a master‐slave model, and a hybrid decomposition algorithm is proposed to address the load imbalance problem. Numerical benchmarks show that both our multiselection method and the parallelization could obviously improve the point selection efficiency. Specifically, total speedups of 20 and 55 are separately obtained for the 3D undulating fish with 10⁶ cell mesh and the 3D rotating hydrofoil with 11 million cell mesh.     

Investigation of an adaptive sampling method for data interpolation using radial basis functions

An adaptive sampling method is presented for optimizing the location of data points in parameter space for multidimensional data interpolation. The method requires a small number of points to begin, and achieves a compromise between space‐filling updates and local refinement in areas where the data are nonlinear, as measured by the Laplacian. A smooth separation function quantifies the sample spacing, and this is blended with the Laplacian to form a criterion on which to assess potential new sample positions. Validation results are presented using two‐dimensional analytic test cases, which demonstrate that the method can recover known optimal designs and gives improvement over data‐independent approaches. In addition, a detailed analysis of the various model parameters is presented. Initial findings are very promising, and it is hoped that further work using the method to generate an aerodynamic database using CFD simulations will lead to a reduction in the number of points required for a given modelling accuracy. Copyright © 2010 John Wiley & Sons, Ltd.     

Educating local radial basis functions using the highest gradient of interest in three dimensional geometries

We present a novel methodology to effectively localize radial basis function approximation methods in three dimensions. The local scheme requires shape parameter‐dependent functions that can be used to approximate gradients of scattered data and to solve partial differential operators. The optimum shape parameter is obtained from the highest gradient of interest, where a known analytical function, when boundary conditions are not present, or a shape parameter‐free global approximation are used to educate the localized scheme. The later option is applicable to problems where the operator needs to be solved multiple times, like in time evolution or stochastic integration. Past shape parameter's optimizations, for two‐dimensional domains, based on the condition number of the interpolant matrix, were unable to provide satisfactory approximations. The applicability of our method is illustrated in the context of an analytical expression interpolation and during a Ginzburg–Landau relaxation of a free energy functional. In general, the optimum shape parameter depends on geometry, node distribution, and density, whereas the approximation errors decrease as the node density and the local stencil size increase. The effective localization of radial basis functions motivates its use in moving boundary problems and accelerates solutions through sparse matrix solvers. Copyright © 2016 John Wiley & Sons, Ltd.     

A Lagrangian finite element approach for the analysis of fluid–structure interaction problems

A Lagrangian finite element method for the analysis of incompressible Newtonian fluid flows, based on a continuous re‐triangulation of the domain in the spirit of the so‐called Particle Finite Element Method, is here revisited and applied to the analysis of the fluid phase in fluid–structure interaction problems. A new approach for the tracking of the interfaces between fluids and structures is proposed. Special attention is devoted to the mass conservation problem. It is shown that, despite its Lagrangian nature, the proposed combined finite element‐particle method is well suited for large deformation fluid–structure interaction problems with evolving free surfaces and breaking waves. The method is validated against the available analytical and numerical benchmarks. Copyright © 2010 John Wiley & Sons, Ltd.     

FE/FMBE coupling to model fluid–structure interaction

In this paper, a finite element (FE)/fast multipole boundary element (FMBE)‐coupling method is presented for modeling fluid–structure interaction problems numerically. Vibrating structures are assumed to consist of elastic or sound absorbing materials. An FE method (FEM) is used for this part of the solution. This structural sub‐domain is embedded in a homogeneous fluid. The case where the boundary of the structural sub‐domain has a very complex geometry is of special interest. In this case, the BE method (BEM) is a more suitable numerical tool than FEM to account for the sound propagation in the homogeneous fluid. The efficiency of the BEM is increased by using FMBEM. The BE‐surface mesh required is directly generated by the FE‐mesh used to discretize the structural sub‐domain and the absorbing material. This FE/FMBE‐coupling method makes it possible to predict the effects of arbitrarily shaped absorbing materials and vibrating structures on the sound field in the surrounding fluid numerically. The coupling method proposed is used to study the acoustic behavior of the lining of an anechoic chamber and that of an entire anechoic chamber in the low‐frequency range. The numerical results obtained are compared with the experimental data. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of transient flow in pipelines with fluid–structure interaction using method of lines

A mathematical model is presented for transient flow in a pipeline with fluid–structure interaction. Water hammer theory and equations of axial motion for the pipeline are employed and the Poisson, junction and transient shear stress couplings are taken into account, which give rise to four coupled non‐linear, first‐order hyperbolic partial differential equations governing the fluid flow and pipe motion. These equations are discretized in space using the Keller box scheme and the method of lines is employed to reduce the partial differential equations to a system of ordinary differential equations. The resulting system is solved using a backward differentiation formulation method. The effect of transient shear stress on transient flow is investigated and the mechanisms underlying this effect are explored. The results revealed that the influence of transient shear stress can be significant and varies considerably, depending on the boundary conditions, viz, valve closure time. Copyright © 2005 John Wiley & Sons, Ltd.     

A rheological interface model and its space–time finite element formulation for fluid–structure interaction

This contribution discusses extended physical interface models for fluid–structure interaction problems and investigates their phenomenological effects on the behavior of coupled systems by numerical simulation. Besides the various types of friction at the fluid–structure interface the most interesting phenomena are related to effects due to additional interface stiffness and damping. The paper introduces extended models at the fluid–structure interface on the basis of rheological devices (Hooke, Newton, Kelvin, Maxwell, Zener). The interface is decomposed into a Lagrangian layer for the solid‐like part and an Eulerian layer for the fluid‐like part. The mechanical model for fluid–structure interaction is based on the equations of rigid body dynamics for the structural part and the incompressible Navier–Stokes equations for viscous flow. The resulting weighted residual form uses the interface velocity and interface tractions in both layers in addition to the field variables for fluid and structure. The weak formulation of the whole coupled system is discretized using space–time finite elements with a discontinuous Galerkin method for time‐integration leading to a monolithic algebraic system. The deforming fluid domain is taken into account by deformable space–time finite elements and a pseudo‐structure approach for mesh motion. The sensitivity of coupled systems to modification of the interface model and its parameters is investigated by numerical simulation of flow induced vibrations of a spring supported fluid‐immersed cylinder. It is shown that the presented rheological interface model allows to influence flow‐induced vibrations. Copyright © 2010 John Wiley & Sons, Ltd.