Source author record

Marek Behr

Marek Behr appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

A multiphysics modeling approach for in-stent restenosis: Theoretical aspects and finite element implementation

Development of in silico models are intrinsic in understanding disease progression in soft biological tissues. Within this work, we propose a fully-coupled Lagrangian finite element framework which replicates the process of in-stent restenosis observed post stent implantation in a coronary artery. Coupled advection-reaction-diffusion reactions are set up that track the evolution of the concentrations of the platelet-derived growth factor, the transforming growth factor-$β$, the extracellular matrix, and the density of the smooth muscle cells. A continuum mechanical description of growth incorporating the evolution of arterial wall constituents is developed, and a suitable finite element implementation discussed. Qualitative validation of the computational model are presented by emulating a stented artery. Patient-specific data can be integrated into the model to predict the risk of restenosis and thereby assist in tuning of stent implantation parameters to mitigate the risk.

preprint2022arXiv

Numerical Prediction and Post-Test Numerical Analysis of the ASDMAD Wind Tunnel Tests in ETW

This paper presents numerical results in comparison with experimental data of nominally static aeroelastic polars of the second ASDMAD (Aero-Structural Dynamics Methods for Airplane Design) campaign conducted in the European Transonic Windtunnel. For the simulation the modular solver package SOFIA combined with the DLR solver TAU on the fluid side and the in-house solver FEAFA on the structural side is applied. First the experimental setup and the aeroelastic solver are introduced, followed by a detailed analysis of the results. Comparisons of pressure distributions and the global lift coefficients are presented for variations of angle of attack, Mach number and aerodynamic loading factor. Additionally, the paper contains preliminary investigations of dynamic results for excited vibration applying periodical inner force couples in the wing root area.

preprint2022arXiv

Spline-Based Space-Time Finite Element Approach for Fluid-Structure Interaction Problems With a Focus on Fully Enclosed Domains

Non-Uniform Rational B-Spline (NURBS) surfaces are commonly used within Computer-Aided Design (CAD) tools to represent geometric objects. When using isogeometric analysis (IGA), it is possible to use such NURBS geometries for numerical analysis directly. Analyzing fluid flows, however, requires complex three-dimensional geometries to represent flow domains. Defining a parametrization of such volumetric domains using NURBS can be challenging and is still an ongoing topic in the IGA community. With the recently developed NURBS-enhanced finite element method (NEFEM), the favorable geometric characteristics of NURBS are used within a standard finite element method. This is achieved by enhancing the elements touching the boundary by using the NURBS geometry itself. In the current work, a new variation of NEFEM is introduced, which is suitable for three-dimensional space-time finite element formulations. The proposed method makes use of a new mapping which results in a non-Cartesian formulation suitable for fluid-structure interaction (FSI). This is demonstrated by combining the method with an IGA formulation in a strongly-coupled partitioned framework for solving FSI problems. The framework yields a fully spline-based representation of the fluid-structure interface through a single NURBS. The coupling conditions at the fluid-structure interface are enforced through a Robin-Neumann type coupling scheme. This scheme is particularly useful when considering incompressible fluids in fully Dirichlet-bounded and curved problems, as it satisfies the incompressibility constraint on the fluid for each step within the coupling procedure. The accuracy and performance of the introduced spline-based space-time finite element approach and its use within the proposed coupled FSI framework are demonstrated using a series of two- and three-dimensional benchmark problems.

preprint2021arXiv

A novel approach to fluid-structure interaction simulations involving large translation and contact

In this work, we present a novel method for the mesh update in flow problems with moving boundaries, the phantom domain deformation mesh update method (PD-DMUM). The PD-DMUM is designed to avoid remeshing; even in the event of large, unidirectional displacements of boundaries. The method combines the concept of two mesh adaptation approaches: (1) The virtual ring shear-slip mesh updatemethod (VR-SSMUM); and (2) the elastic mesh update method (EMUM). As in the VR-SSMUM, the PD-DMUMextends the fluid domain by a phantom domain; the PD-DMUM can thus locally adapt the element density. Combined with the EMUM, the PD-DMUMallows the consideration of arbitrary boundary movements. In this work, we apply the PD-DMUM in two test cases. Within the first test case, we validate the PD-DMUM in a 2D Poiseuille flow on a moving background mesh. Subsequently the fluid-structure interaction (FSI) problem in the second test case serves as a proof of concept. More, we stress the advantages of the novel method with regard to conventional mesh update approaches.

preprint2021arXiv

Numerical Design of Distributive Mixing Elements

This paper presents a novel shape-optimization technique for the design of mixing elements in single-screw extruders. Extruders enable the continuous production of constant-cross-section profiles. Equipped with one or several screw-shaped rotors, extruders transport polymer particles towards the outlet. Due to shear heating, melting is induced and a melt stream is created, which can be further processed. While many variants of multi-screw extruders exist, a significant share of all extrusion machines is made up of single-screw extruders due to their comparatively low operating costs and complexity. While the reduced complexity yields economic benefits, single-screw extruders' mixing capabilities, i.e., their ability to produce a melt with a homogeneous material and temperature distribution, suffer compared to multi-screw extruders. To compensate for this shortcoming, so-called mixing elements are added to the screw to enhance mixing by recurring flow reorientations. In view of the largely unintuitive flow characteristics of polymer melts, we present an optimization framework that allows designing these mixing elements numerically based on finite-element simulations of the melt flow. To reduce the computational demand required by shape optimization of a complete mixing section, we only focus on the shape optimization of a single mixing element. This paper presents advances in three aspects of numerical design: (1) A combination of free-form deformation and surface splines is presented, allowing to parameterize the mixing element's shape by very few variables. (2) The combination of this concept with a linear-elasticity-based mesh update method to deform the computational domain without the need for remeshing is demonstrated. (3) A simple yet robust and sensitive objective formulation to assess distributive mixing in laminar flows based on a measure for the interfacial area is proposed.

preprint2020arXiv

A Multi-Vector Interface Quasi-Newton Method with Linear Complexity for Partitioned Fluid-Structure Interaction

In recent years, interface quasi-Newton methods have gained growing attention in the fluid-structure interaction community by significantly improving partitioned solution schemes: They not only help to control the inherent added-mass instability, but also prove to substantially speed up the coupling's convergence. In this work, we present a novel variant: The key idea is to build on the multi-vector Jacobian update scheme first presented by Bogaers et al. (2014) and avoid any explicit representation of the (inverse) Jacobian approximation, since it slows down the solution for large systems. Instead, all terms involving a quadratic complexity have been systematically eliminated. The result is a new multi-vector interface quasi-Newton variant whose computational cost scales linearly with the problem size.

preprint2020arXiv

Finite-Element Formulation for Advection-Reaction Equations with Change of Variable and Discontinuity Capturing

We propose a change of variable approach and discontinuity capturing methods to ensure physical constraints for advection-reaction equations discretized by the finite element method. This change of variable confines the concentration below an upper bound in a very natural way. For the non-negativity constraint, we propose to use a discontinuity capturing method defined on the reference element that is combined with an anisotropic crosswind-dissipation operator. This discontinuity capturing cannot completely eliminate negative values but effectively minimizes their occurrence. The proposed methods are applied to different biophysical models and show a good agreement with experimental results for the FDA benchmark blood pump for a physiological red blood cell pore formation model.

preprint2014arXiv

Fully-implicit log-conformation formulation of constitutive laws

Subject of this paper is the derivation of a new constitutive law in terms of the logarithm of the conformation tensor that can be used as a full substitute for the 2D governing equations of the Oldroyd-B, Giesekus and other models. One of the key features of these new equations is that - in contrast to the original log-conf equations given by Fattal and Kupferman (2004) - these constitutive equations combined with the Navier-Stokes equations constitute a self-contained, non-iterative system of partial differential equations. In addition to its potential as a fruitful source for understanding the mathematical subtleties of the models from a new perspective, this analytical description also allows us to fully utilize the Newton-Raphson algorithm in numerical simulations, which by design should lead to reduced computational effort. By means of the confined cylinder benchmark we will show that a finite element discretization of these new equations delivers results of comparable accuracy to known methods.

preprint2014arXiv

Transition Prediction for Scramjet Intakes Using the γ-Re_θ_t Model Coupled to Two Turbulence Models

Due to the thick boundary layers in hypersonic flows, the state of the boundary layer significantly influences the whole flow field as well as surface heat loads. Hence, for engineering applications the efficient numerical prediction of laminar-to-turbulent transition is a challenging and important task. Within the framework of the Reynolds averaged Navier-Stokes equations, Langtry/Menter [1] proposed the $γ$-$Re_{θ_t}$ transition model using two transport equations for the intermittency and $γ$-$Re_{θ_t}$ combined with the Shear Stress Transport turbulence model (SST) [2]. The transition model contains two empirical correlations for onset and length of transition. Langtry/Menter [1] designed and validated the correlations for the subsonic and transonic flow regime. For our applications in the hypersonic flow regime, the development of a new set of correlations proved necessary, even when using the same SST turbulence model [3]. Within this paper, we propose a next step and couple the transition model with the SSG/LRR-$ω$ Reynolds stress turbulence model [4] which we found to be well suited for scramjet intake simulations. First, we illustrate the necessary modifications of the Reynolds stress model and the hypersonic in-house correlations using a hypersonic flat plate test case. Next, the transition model is successfully validated for its use coupled to both turbulence models using a hypersonic double ramp test case. Regardless of the turbulence model, the transition model is able to correctly predict the transition process compared to experimental data. In addition, we apply the transition model combined with both turbulence models to three different fully 3D scramjet intake configurations which are experimentally investigated in wind tunnel facilities. The agreement with the available experimental data is also shown.

preprint2013arXiv

Wavelet-based Adaptive Techniques Applied to Turbulent Hypersonic Scramjet Intake Flows

The simulation of hypersonic flows is computationally demanding due to large gradients of the flow variables caused by strong shock waves and thick boundary or shear layers. The resolution of those gradients imposes the use of extremely small cells in the respective regions. Taking turbulence into account intensives the variation in scales even more. Furthermore, hypersonic flows have been shown to be extremely grid sensitive. For the simulation of three-dimensional configurations of engineering applications, this results in a huge amount of cells and prohibitive computational time. Therefore, modern adaptive techniques can provide a gain with respect to computational costs and accuracy, allowing the generation of locally highly resolved flow regions where they are needed and retaining an otherwise smooth distribution. An h-adaptive technique based on wavelets is employed for the solution of hypersonic flows. The compressible Reynolds averaged Navier-Stokes equations are solved using a differential Reynolds stress turbulence model, well suited to predict shock-wave-boundary-layer interactions in high enthalpy flows. Two test cases are considered: a compression corner and a scramjet intake. The compression corner is a classical test case in hypersonic flow investigations because it poses a shock-wave-turbulent-boundary-layer interaction problem. The adaptive procedure is applied to a two-dimensional confguration as validation. The scramjet intake is firstly computed in two dimensions. Subsequently a three-dimensional geometry is considered. Both test cases are validated with experimental data and compared to non-adaptive computations. The results show that the use of an adaptive technique for hypersonic turbulent flows at high enthalpy conditions can strongly improve the performance in terms of memory and CPU time while at the same time maintaining the required accuracy of the results.