A partitioned solution approach for the fluid-structure interaction of wind and thin-walled structures

In many engineering applications two or more different interacting systems require the numer-ical solution of so-called multifield problems. In civil engineering the interaction of fluid and structures plays an important role, i.e. for fabric tensile structures of light and flexible materials often used for large roof systems, capacious umbrellas or canopies. Whereas powerful numerical simulation techniques have been established in structural engineering as well as in fluid mechan-ics, only relatively few approaches to simulate the interaction of fluids with civil engineering constructions have been presented. To determine the wind loads on complex structures, it is still state-of-the-art to apply semi-empirical, strongly simplifying methods or to perform expensive ex-periments in wind tunnels. In this paper an approach of a coupled fluid-structure simulation will be presented for membrane and thin shell structures. The interaction is described by the struc-tural deformation as response to wind forces, resulting in a modification of the fluid flow domain. Besides a realistic determination of the wind loads, information on the structural stability can be obtained. The so-called partitioned solution is based on an iterative frame algorithm, integrating different codes for Computational Fluid Dynamics (CFD) and for Computational Structural Dy-namics (CSD) in an explicit or an implicit time-stepping procedure. All data exchange between the two different applications is performed via a neutral geometric model provided by a coupling interface. A conservative interpolation method is used for the interpolation of the nodal loads. The time-dependent motion of the structure requires a dynamic modification of the different grids and a redefinition of the Navier-Stokes equations in an Arbitrary Langrangian Eulerian (ALE) formulation. As an example for the present implementation, results of a coupled fluid-structure simulation for a textile membrane canopy will be presented

Inverting elastic dislocations using the Weakly-enforced Slip Method

Earthquakes cause lasting changes in static equilibrium, resulting in global deformation fields that can be observed. Consequently, deformation measurements such as those provided by satellite based InSAR monitoring can be used to infer an earthquake's faulting mechanism. This inverse problem requires a numerical forward model that is both accurate and fast, as typical inverse procedures require many evaluations. The Weakly-enforced Slip Method (WSM) was developed to meet these needs, but it was not before applied in an inverse problem setting. Consequently, it was unknown what effect particular properties of the WSM, notably its inherent continuity, have on the inversion process. Here we show that the WSM is able to accurately recover slip distributions in a Bayesian-inference setting, provided that data points in the vicinity of the fault are removed. In a representative scenario, an element size of 2 km was found to be sufficiently fine to generate a posterior probability distribution that is close to the theoretical optimum. For rupturing faults a masking zone of 20 km sufficed to avoid numerical disturbances that would otherwise be induced by the discretization error. These results demonstrate that the WSM is a viable forward method for earthquake inversion problems. While our synthesized scenario is basic for reasons of validation, our results are expected to generalize to the wider gamut of scenarios that finite element methods are able to capture. This has the potential to bring modeling flexibility to a field that if often forced to impose model restrictions in a concession to computability.Comment: The associated software implementation is openly available in zenodo at https://doi.org/10.5281/zenodo.507179