FE/BE coupling for an acoustic fluid-structure interaction problem. Residual a posteriori error estimates

This is the author's accepted manuscript. The final published article is available from the link below. Copyright © 2011 John Wiley & Sons, Ltd.In this paper, we developed an a posteriori error analysis of a coupling of finite elements and boundary elements for a fluid–structure interaction problem in two and three dimensions. This problem is governed by the acoustic and the elastodynamic equations in time-harmonic vibration. Our methods combined integral equations for the exterior fluid and FEMs for the elastic structure. It is well-known that because of the reduction of the boundary value problem to boundary integral equations, the solution is not unique in general. However, because of superposition of various potentials, we consider a boundary integral equation that is uniquely solvable and avoids the irregular frequencies of the negative Laplacian operator of the interior domain. In this paper, two stable procedures were considered; one is based on the nonsymmetric formulation and the other is based on a symmetric formulation. For both formulations, we derived reliable residual a posteriori error estimates. From the estimators we computed local error indicators that allowed us to develop an adaptive mesh refinement strategy. For the two-dimensional case we performed an adaptive algorithm on triangles, and for the three-dimensional case we used hanging nodes on hexahedrons. Numerical experiments underline our theoretical results.DFG German Research Foundatio

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial value problem. Progress in characteristic evolution is traced from the early stage of 1D feasibility studies to 2D axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to extend such simulations to null infinity where the waveform from the binary inspiral and merger can be unambiguously computed. This has now been accomplished by Cauchy-characteristic extraction, where data for the characteristic evolution is supplied by Cauchy data on an extraction worldtube inside the artificial outer boundary. The ultimate application of characteristic evolution is to eliminate the role of this outer boundary by constructing a global solution via Cauchy-characteristic matching. Progress in this direction is discussed.Comment: New version to appear in Living Reviews 2012. arXiv admin note: updated version of arXiv:gr-qc/050809

Transparent boundary conditions based on the Pole Condition for time-dependent, two-dimensional problems

The pole condition approach for deriving transparent boundary conditions is extended to the time-dependent, two-dimensional case. Non-physical modes of the solution are identified by the position of poles of the solution's spatial Laplace transform in the complex plane. By requiring the Laplace transform to be analytic on some problem dependent complex half-plane, these modes can be suppressed. The resulting algorithm computes a finite number of coefficients of a series expansion of the Laplace transform, thereby providing an approximation to the exact boundary condition. The resulting error decays super-algebraically with the number of coefficients, so relatively few additional degrees of freedom are sufficient to reduce the error to the level of the discretization error in the interior of the computational domain. The approach shows good results for the Schr\"odinger and the drift-diffusion equation but, in contrast to the one-dimensional case, exhibits instabilities for the wave and Klein-Gordon equation. Numerical examples are shown that demonstrate the good performance in the former and the instabilities in the latter case

A simple preconditioned domain decomposition method for electromagnetic scattering problems

We present a domain decomposition method (DDM) devoted to the iterative solution of time-harmonic electromagnetic scattering problems, involving large and resonant cavities. This DDM uses the electric field integral equation (EFIE) for the solution of Maxwell problems in both interior and exterior subdomains, and we propose a simple preconditioner for the global method, based on the single layer operator restricted to the fictitious interface between the two subdomains.Comment: 23 page