Now showing1 - 9 of 9
  • Some of the metrics are blocked by your 
    Item type:Publication, Access status: with files ,
    Adaptive tolerances for staggered solution algorithms in fluid-structure interaction
    Staggered solution algorithms are a well-known alternative to monolithic approaches for solving strongly coupled multi-field problems. A coupling between the subproblems is achieved by exchanging coupling quantities between the solvers. Implicit coupling schemes are realized by letting the solvers solve each load or time step repeatedly until convergence up to a given coupling tolerance is achieved. In fluid-structure interaction as well as a variety of other problem classes, the equations governing the individual fields are nonlinear. Accordingly, each solver performs an inner iterative solution procedure that terminates once an inner tolerance is reached. The basic idea of this work builds on the possibility to adaptively adjust these inner tolerances based on carefully designed rules while preserving the black-box nature of the solvers. The resulting coupling scheme yields significant improvements in computational efficiency compared to classical schemes where the inner solver tolerances are held fixed. This is demonstrated in several numerical examples. The idea is tested in combination with state-of-the-art convergence acceleration schemes and can be realized within any staggered solution approach by only minor modifications to the participating solvers.
    Publicationtype: Journal Article
      56  29
  • Some of the metrics are blocked by your 
    Item type:Publication, Access status: with files ,
    Adjoint-assisted robust shape optimization of an idealized arterial bypass graft using the FOSM method
    (Springer, 2025-06-01) ; ;
    This paper presents the shape optimization of an idealized arterial bypass graft under uncertainties. The underlying blood flow problem is numerically solved using computational fluid dynamics (CFD) simulations, able to account for mechanical hemolysis and non-Newtonian viscosity properties of blood. The employed hemolysis and non-Newtonian models utilize a set of parameters that are considered uncertain in this work. To this end, the optimization problem is rendered robust and targets the minimization of the first two statistical moments of a hemolysis index. The propagation of the uncertainties to the hemolysis index is done through the first-order second-moment method (FOSM). The necessary derivatives are computed by the adjoint method and finite differences. Several steady-state robust shape optimization simulations of the idealized bypass graft are presented. Selected optimized shapes are further assessed by means of additional fluid–structure interaction (FSI) simulations under unsteady conditions.
    Publicationtype: Journal Article
      30  22
  • Some of the metrics are blocked by your 
    Item type:Publication,
    On the efficiency of explicit and semi-explicit immersed boundary finite element methods for wave propagation problems
    (American Institute of Mathematical Sciences, 2025-03-01) ; ; ; ;
    Immersed boundary methods have attracted substantial interest in the last decades due to their potential for computations involving complex geometries. Often these cannot be efficiently discretized using boundary-fitted finite elements. Immersed boundary methods provide a simple and fully automatic discretization based on Cartesian grids and tailored quadrature schemes that account for the geometric model. It can thus be described independently of the grid, e.g., by image data obtained from computed tomography scans. The drawback of such a discretization lies in the potentially small overlap between certain elements in the grid and the geometry. These badly cut elements with small physical support pose a particular challenge for nonlinear and/or dynamic simulations. In this work, we focus on problems in structural dynamics and acoustics and concentrate on solving them with explicit time-marching schemes. In this context, badly cut elements can lead to unfeasibly small critical time step sizes. We investigate the performance of implicit-explicit time marching schemes and two stabilization methods developed in previous works as potential remedies. While these have been studied before with regard to their effectiveness in in-creasing the critical time step size, their numerical efficiency has only been considered in terms of accuracy per degree of freedom. In this paper, we eval-uate the computation time required for a given accuracy, which depends not only on the number of degrees of freedom but also on the selected spatial dis-cretization, the sparsity patterns of the system matrices, and the employed time-marching scheme.
    Publicationtype: Journal Article
      15
  • Some of the metrics are blocked by your 
    Item type:Publication,
    Hemodynamics in arterial bypass graft anastomoses with varying cuff sizes and proximal flow paths: a fluid–structure interaction study
    (Taylor & Francis, 2025) ; ;
    This article investigates the effect of the cuff size of arterial bypass grafts and the flow conditions on the hemodynamics in the anastomosis (connection) to the artery, using numerical simulations. We consider a fluid-structure interaction problem which is solved based on a partitioned scheme. Additionally, we employ computational fluid dynamics to investigate the effect of a rigid wall assumption. The work focuses on clinically relevant hemodynamic quantities associated with the development of intimal hyperplasia. We also include a model for the prediction of hemolysis into the simulation. The results show that even minor changes of the cuff size can result into significant differences in the corresponding quantities of interest. The importance of the inflow path is shown to be lower than that of the cuff size. The usually employed rigid wall assumption is found to be adequate to address wall shear stress oscillations but falls short on predicting maximum and minimum wall shear stress-related quantities of interest.
    Publicationtype: Journal Article
      30
  • Some of the metrics are blocked by your 
    Item type:Publication, Access status: with files ,
    Experimental characterization of acoustic damping materials
    Due to their ease of use and low cost, passive damping methods are a preferred mean for the reduction of noise in many engineering applications. This applies in particular to foam materials, which exhibit good acoustic and mechanical damping properties. The selection of suitable materials is usually carried out experimentally and can be very labor‐intensive and time‐consuming. For this reason, it is helpful to develop qualified numerical methodsthat can be used for material design and selection. Hence, foam materials must be characterized experimentally in order to enable a later comparison with vibroacoustic simulations. This contribution, therefore, aims at providing suitable parameters for the use in numerical analyses and their validation. Firstly, the microstructure of a foam specimen is captured by means of a CT scan. In addition to providing the geometry for the multi‐physics simulations, this measurement is also used to determine characteristic foam features, for example, strut thickness and pore size distribution. In the second step, the frequency‐dependent stiffness and damping properties of the material are determined by a special experimental setup utilizing an electrodynamic shaker. Here, the dynamic system is approximated as a single‐mass oscillator, which is sufficiently accurate for low frequencies. These properties will later be used in the numerical model to evaluate different parameter identification approaches. In the third and last step of the experimental campaign, measurements with an impedance tube are conducted to obtain the coefficient of absorption. This material parameter is particularly suitable for comparing experiments and simulations. Finally, the correlation between the experimental results is examined to provide a deeper understanding of the foam materials.
    Publicationtype: Journal Article
      18  19
  • Some of the metrics are blocked by your 
    Item type:Publication, Access status: with files ,
    An automatic simulation pipeline for coupled simulations of acoustic damping materials
    Foamed materials are widely used to reduce noise due to their comparably good acoustic damping behavior. However, out of a large variety of these materials a suitable candidate has to be identified for each application. This is a challenging process that is typically guided by experiments and experience. While numerical simulations could support these experiments and reduce the effort to a great extent, no suitable discretization approach has yet been established that can fully capture the complex geometry of the foam. A fully resolved model is desirable in order to yield reliable predictions that can then be used to establish homogenized models. We established a monolithic coupling approach based on the finite cell method (FCM) that realizes a vibroacoustic simulation in this sense. The fluid and the structure domain are discretized by Cartesian grids and the geometry defined based on computed tomography scans is accounted for during the quadrature of the weak form. Our simulation in the time domain makes use of explicit time marching schemes and is therefore limited by a critical time step size. This is known to be arbitrarily low for discretizations with the FCM containing cells with arbitrarily small support. As a remedy against this we use the classical ‐stabilization technique and investigate its potentials and limitations.
    Publicationtype: Journal Article
      33  14
  • Some of the metrics are blocked by your 
    Item type:Publication, Access status: with files ,
    An analysis of high order FEM and IGA for explicit dynamics : mass lumping and immersed boundaries
    (Wiley, 2024-08-30) ; ;
    Hughes, Thomas J.R.
    ;
    ;
    Sangalli, Giancarlo
    We investigate the behavior of different shape functions for the discretization of hyperbolic problems. In particular, we consider classical Lagrange polynomials and B-splines. The studies focus on the performance of the these functions as a spatial discretization approach combined with an explicit time marching scheme. In this regard, a major concern is the maximum eigenvalue that imposes restrictions on the critical time step size and suitable lumping techniques that yield a diagonal mass matrix. The accuracy of the discretization methods is assessed in an asymptotic manner in terms of the convergence of eigenvalues and eigenvectors. Further, the global accuracy is investigated in terms of the full spectrum. The results show that B-spline discretization with a consistent mass matrix are more accurate than those based on Lagrange shape functions, which holds true in the boundary-fitted as well as in the immersed setting. On the other hand, Lagrange shape functions are more robust with respect to standard lumping techniques, which cannot be directly applied for B-splines without loss of accuracy. In general, we observe that none of the standard lumping schemes yields optimal results for B-splines, even in the boundary-fitted setting. For the immersed setting, also Lagrange shape functions show a drop in accuracy which depends on the position of the boundary that cuts the element. Several remedies are considered in order to overcome these issues, including interpolatory B-spline bases as well as eigenvalue stabilization methods. While accuracy and stability can be improved using these remedies, we conclude from our study that it is still an open question, how to design a discretization method that achieves large critical time step sizes in combination with a diagonal mass matrix and high accuracy in the immersed setting. We note that these considerations primarily relate to linear structural dynamics applications, such as for example, structural acoustics. In nonlinear problems, such as automotive crash dynamics, other considerations predominate. An example of a one-dimensional elastic-plastic bar impacting a rigid wall is illustrative.
    Publicationtype: Journal Article
      26  16
  • Some of the metrics are blocked by your 
    Item type:Publication, Access status: with files ,
    An eigenvalue stabilization technique for immersed boundary finite element methods in explicit dynamics
    (Elsevier, 2024-07-15) ; ; ;
    Löhnert, S.
    ;
    The application of immersed boundary methods in static analyses is often impeded by poorly cut elements (small cut elements problem), leading to ill-conditioned linear systems of equations and stability problems. While these concerns may not be paramount in explicit dynamics, a substantial reduction in the critical time step size based on the smallest volume fraction χ of a cut element is observed. This reduction can be so drastic that it renders explicit time integration schemes impractical. To tackle this challenge, we propose the use of a dedicated eigenvalue stabilization (EVS) technique. The EVS-technique serves a dual purpose. Beyond merely improving the condition number of system matrices, it plays a pivotal role in extending the critical time increment, effectively broadening the stability region in explicit dynamics. As a result, our approach enables robust and efficient analyses of high-frequency transient problems using immersed boundary methods. A key advantage of the stabilization method lies in the fact that only element-level operations are required. This is accomplished by computing all eigenvalues of the element matrices and subsequently introducing a stabilization term that mitigates the adverse effects of cutting. Notably, the stabilization of the mass matrix Mc of cut elements – especially for high polynomial orders p of the shape functions – leads to a significant raise in the critical time step size Δtcr. To demonstrate the efficiency of our technique, we present two specifically selected dynamic benchmark examples related to wave propagation analysis, where an explicit time integration scheme must be employed to leverage the increase in the critical time step size.
    Publicationtype: Journal Article
      23  195
  • Some of the metrics are blocked by your 
    Item type:Publication,
    Vibroacoustic simulations of acoustic damping materials using a fictitious domain approach
    (Elsevier, 2024-01-06) ;
    Marter, Paul
    ;
    ; ;
    The numerical investigation of acoustic damping materials, such as foams, constitutes a valuable enhancement to experimental testing. Typically, such materials are modeled in a homogenized way in order to reduce the computational effort and to circumvent the need for a computational mesh that resolves the complex micro-structure. However, to gain detailed insight into the acoustic behavior, e.g., the transmittance of noise, such fully resolved models are mandatory. The meshing process can still be drastically simplified by using a fictitious domain approach. We propose the finite cell method, which combines the fictitious domain approach with high-order finite elements and resolves the complex geometry using special quadrature rules. In order to take into account the fluid-filled pores of a typical damping material, a coupled vibroacoustic problem needs to be solved. To this end, we construct two separate finite cell discretizations and prescribe coupling conditions at the interface in the usual manner. The only difference to a classical boundary fitted approach to vibroacoustics is that the fluid-solid interface is immersed into the respective discretization and does not correspond to the element boundaries. The proposed enhancement of the finite cell method for vibroacoustics is verified based on a comparison with commercial software and used within an exemplary application.
    Publicationtype: Journal Article
      13