Bernd Simeon

Alexander Shamanskiy, Bernd Simeon

Isogeometric Analysis (IGA) is a recently introduced computational approach intended to breach the gap between the Finite Element Analysis and the Computer Aided Design worlds. In this work, we apply it to numerically simulate thermal expansion of oil-free twin screw compressors in operation. High global smoothness of IGA leads to a more accurate representation of the compressor geometry. We utilize standard tri-variate B-splines to parametrize the rotors, while the casing is modeled exactly by using NURBS. We employ the Galerkin version of IGA to solve the thermal expansion problem in the stationary case. The results allow to estimate the contraction of the clearance space between the casing and the rotors. The implementation is based on the open source C++ library G+Smo. This work is supported by the European Union within the Project MOTOR: Multi-ObjecTive design Optimization of fluid eneRgy machines.

Clarissa Arioli, Alexander Shamanskiy, Sven Klinkel et al.

This paper deals with a special class of parametrizations for Isogeometric Analysis (IGA). The so-called scaled boundary parametrizations are easy to construct and particularly attractive if only a boundary description of the computational domain is available. The idea goes back to the Scaled Boundary Finite Element Method (SB-FEM), which has recently been extended to IGA. We take here a different viewpoint and study these parametrizations as bivariate or trivariate B-spline functions that are directly suitable for standard Galerkin-based IGA. Our main results are first a general framework for this class of parametrizations, including aspects such as smoothness and regularity as well as generalizations to domains that are not star-shaped. Second, using the Poisson equation as example, we explain the relation between standard Galerkin-based IGA and the Scaled Boundary IGA by means of the Laplace-Beltrami operator. Further results concern the separation of integrals in both approaches and an analysis of the singularity in the scaling center. Among the computational examples we present a planar rotor geometry that stems from a screw compressor machine and compare different parametrization strategies.

Alexander Shamanskiy, Michael Helmut Gfrerer, Jochen Hinz et al.

The construction of volumetric parametrizations for computational domains is a key step in the pipeline of isogeometric analysis. Here, we investigate a solution to this problem based on the mesh deformation approach. The desired domain is modeled as a deformed configuration of an initial simple geometry. Assuming that the parametrization of the initial domain is bijective and that it is possible to find a locally invertible displacement field, the method yields a bijective parametrization of the target domain. We compute the displacement field by solving the equations of nonlinear elasticity with the neo-Hookean material law, and we show an efficient variation of the incremental loading algorithm tuned specifically to this application. In order to construct the initial domain, we simplify the target domain's boundary by means of an L2-projection onto a coarse basis and then apply the Coons patch approach. The proposed methodology is not restricted to a single patch scenario but can be utilized to construct multi-patch parametrizations with naturally looking boundaries between neighboring patches. We illustrate its performance and compare the result to other established parametrization approaches on a range of two-dimensional and three-dimensional examples.

Jan Kleinert, Bernd Simeon

The present article presents a summarizing view at differential-algebraic equations (DAEs) and analyzes how new application fields and corresponding mathematical models lead to innovations both in theory and in numerical analysis for this problem class. Recent numerical methods for nonsmooth dynamical systems subject to unilateral contact and friction illustrate the topicality of this development.

Felix Dietrich, Dennis Merkert, Bernd Simeon

In this paper, we first introduce the reader to the Basic Scheme of Moulinec and Suquet in the setting of quasi-static linear elasticity, which takes advantage of the fast Fourier transform on homogenized microstructures to accelerate otherwise time-consuming computations. By means of an asymptotic expansion, a hierarchy of linear problems is derived, whose solutions are looked at in detail. It is highlighted how these generalized homogenization problems depend on each other. We extend the Basic Scheme to fit this new problem class and give some numerical results for the first two problem orders.