We present a coupled simulation approach for fluid–structure–acoustic interactions (FSAI) as an example for strongly surface coupled multi-physics problems. In addition to the multi-physics character, FSAI feature multi-scale properties as a further challenge. In our partitioned approach, the problem is split into spatially separated subdomains interacting via coupling surfaces. Within each subdomain, scalable, single-physics solvers are used to solve the respective equation systems. The surface coupling between them is realized with the scalable open-source coupling tool preCICE described in the “Partitioned Fluid–Structure– Acoustics Interaction on Distributed Data: Coupling via preCICE”. We show how this approach enables the use of existing solvers and present the overall scaling behavior for a three-dimensional test case with a bending tower generating acoustic waves.We run this simulation with different solvers demonstrating the performance of various solvers and the flexibility of the partitioned approach with the coupling tool preCICE. An efficient and scalable in-situ visualization reducing the amount of data in place at the simulation processors before sending them over the network or to a file system completes the simulation environment.@en

Fluid-structure interaction simulations are crucial for many engineering problems. For example, the blood flow around new heart valves or the deployment of airbags during a car crash are often modeled with fluidstructure interaction simulations. Also, to design safe parachutes, simulations are carried out to model the unsteady deformations of the parachute during a jump. Thus, there is an apparent need for multi-physics software codes which can model fluid-structure interaction problems. However, current state-of-the-art solvers cannot be used for design or optimization studies of for example aircraft structures due to long simulation times. This is mainly caused by a large number of coupling iterations needed to reach convergence within each time step for a strongly coupled fluid-structure interaction simulation. Also, a large number of time steps are required to reach an acceptable accuracy in time for unsteady simulations. Hence, there is an urgency for efficiency improvements of fluid-structure interaction solvers. In this thesis, two approaches are investigated to decrease the computational times for a fluid-structure interaction simulation: multi-level acceleration of the coupled problem, and the use of higher order time integration schemes.@en

In this article, we endeavour to find a fast solver for finite volume discretizations for compressible unsteady viscous flows. Thereby, we concentrate on comparing the efficiency of important classes of time integration schemes, namely time adaptive Rosenbrock, singly diagonally implicit (SDIRK) and explicit first stage singly diagonally implicit Runge-Kutta (ESDIRK) methods. To make the comparison fair, efficient equation system solvers need to be chosen and a smart choice of tolerances is needed. This is determined from the tolerance TOL that steers time adaptivity. For implicit Runge-Kutta methods, the solver is given by preconditioned inexact Jacobian-free Newton-Krylov (JFNK) and for Rosenbrock, it is preconditioned Jacobian-free GMRES. To specify the tolerances in there, we suggest a simple strategy of using TOL/100 that is a good compromise between stability and computational effort. Numerical experiments for different test cases show that the fourth order Rosenbrock method RODASP and the fourth order ESDIRK method ESDIRK4 are best for fine tolerances, with RODASP being the most robust scheme.@en

Radial Basis Function (RBF) mesh deformation is one of the most robust mesh deformation methods available. Using the greedy (data reduction) method in combination with an explicit boundary correction, results in an efficient method as shown in literature. However, to ensure the method remains robust, two issues are addressed: 1) how to ensure that the set of control points remains an accurate representation of the geometry in time and 2) how to use/automate the explicit boundary correction, while ensuring a high mesh quality. In this paper, we propose an adaptive RBF mesh deformation method, which ensures the set of control points always represents the geometry/displacement up to a certain (user-specified) criteria, by keeping track of the boundary error throughout the simulation and re-selecting when needed. Opposed to the unit displacement and prescribed displacement selection methods, the adaptive method is more robust, user-independent and efficient, for the cases considered. Secondly, the analysis of a single high aspect ratio cell is used to formulate an equation for the correction radius needed, depending on the characteristics of the correction function used, maximum aspect ratio, minimum first cell height and boundary error. Based on the analysis two new radial basis correction functions are derived and proposed. This proposed automated procedure is verified while varying the correction function, Reynolds number (and thus first cell height and aspect ratio) and boundary error. Finally, the parallel efficiency is studied for the two adaptive methods, unit displacement and prescribed displacement for both the CPU as well as the memory formulation with a 2D oscillating and translating airfoil with oscillating flap, a 3D flexible locally deforming tube and deforming wind turbine blade. Generally, the memory formulation requires less work (due to the large amount of work required for evaluating RBF's), but the parallel efficiency reduces due to the limited bandwidth available between CPU and memory. In terms of parallel efficiency/scaling the different studied methods perform similarly, with the greedy algorithm being the bottleneck. In terms of absolute computational work the adaptive methods are better for the cases studied due to their more efficient selection of the control points. By automating most of the RBF mesh deformation, a robust, efficient and almost user-independent mesh deformation method is presented.@en

The partitioned simulation of fluid-structure interactions offers great flexibility in terms of exchanging flow and structure solver and using existing established codes. However,it often suffers from slow convergence and limited parallel scalability. Quasi-Newton or accelerated fixed-point iterations are a very efficient way to solve the convergence issue. At the same time,they stabilize and speed up not only the standard staggered fluid-structure coupling iterations,but also the variant with simultaneous execution of flow and structure solver that is fairly inefficient if no acceleration methods for the underlying fixed-point iteration are used. In this chapter,we present a review on combinations of iteration patterns (parallel and staggered) and of quasi-Newton methods and compare their suitability in terms of convergence speed,robustness,and parallel scalability. Some of these variants use the so-called manifold mapping that yields an additional speedup by using an approach that can be interpreted as a generalization of the multi-level idea.@en