-
Implicit-Explicit Time Integration for the Immersed Wave Equation
Authors:
Christian Faßbender,
Tim Bürchner,
Philipp Kopp,
Ernst Rank,
Stefan Kollmannsberger
Abstract:
Immersed boundary methods simplify mesh generation by embedding the domain of interest into an extended domain that is easy to mesh, introducing the challenge of dealing with cells that intersect the domain boundary. Combined with explicit time integration schemes, the finite cell method introduces a lower bound for the critical time step size. Explicit transient analyses commonly use the spectral…
▽ More
Immersed boundary methods simplify mesh generation by embedding the domain of interest into an extended domain that is easy to mesh, introducing the challenge of dealing with cells that intersect the domain boundary. Combined with explicit time integration schemes, the finite cell method introduces a lower bound for the critical time step size. Explicit transient analyses commonly use the spectral element method due to its natural way of obtaining diagonal mass matrices through nodal lumping. Its combination with the finite cell method is called the spectral cell method. Unfortunately, a direct application of nodal lumping in the spectral cell method is impossible due to the special quadrature necessary to treat the discontinuous integrand inside the cut cells. We analyze an implicit-explicit (IMEX) time integration method to exploit the advantages of the nodal lumping scheme for uncut cells on one side and the unconditional stability of implicit time integration schemes for cut cells on the other. In this hybrid, immersed Newmark IMEX approach, we use explicit second-order central differences to integrate the uncut degrees of freedom that lead to a diagonal block in the mass matrix and an implicit trapezoidal Newmark method to integrate the remaining degrees of freedom (those supported by at least one cut cell). The immersed Newmark IMEX approach preserves the high-order convergence rates and the geometric flexibility of the finite cell method. We analyze a simple system of spring-coupled masses to highlight some of the essential characteristics of Newmark IMEX time integration. We then solve the scalar wave equation on two- and three-dimensional examples with significant geometric complexity to show that our approach is more efficient than state-of-the-art time integration schemes when comparing accuracy and runtime.
△ Less
Submitted 23 October, 2023;
originally announced October 2023.
-
Isogeometric Multi-Resolution Full Waveform Inversion based on the Finite Cell Method
Authors:
Tim Bürchner,
Philipp Kopp,
Stefan Kollmannsberger,
Ernst Rank
Abstract:
Full waveform inversion (FWI) is an iterative identification process that serves to minimize the misfit of model-based simulated and experimentally measured wave field data, with the goal of identifying a field of parameters for a given physical object. The inverse optimization process of FWI is based on forward and backward solutions of the (elastic or acoustic) eave equation. In a previous paper…
▽ More
Full waveform inversion (FWI) is an iterative identification process that serves to minimize the misfit of model-based simulated and experimentally measured wave field data, with the goal of identifying a field of parameters for a given physical object. The inverse optimization process of FWI is based on forward and backward solutions of the (elastic or acoustic) eave equation. In a previous paper [1], we explored opportunities of using the finite cell method (FCM) as the wave field solver to incorporate highly complex geometric models. Furthermore, we demonstrated that the identification of the model's density outperforms that of the velocity -- particularly in cases where unknown voids characterized by homogeneous Neumann boundary conditions need to be detected. The paper at hand extends this previous study: The isogeometric finite cell analysis (IGA-FCM) -- a combination of isogeometric analysis (IGA) and FCM -- is applied for the wave field solver, with the advantage that the polynomial degree and subsequently also the sampling frequency of the wave field can be increased quite easily. Since the inversion efficiency strongly depends on the accuracy of the forward and backward wave field solution and of the gradient of the functional, consistent and lumped mass matrix discretization are compared. The resolution of the grid describing the unknown material density is the decouple from the knot span grid. Finally, we propose an adaptive multi-resolution algorithm that refines the material grid only locally using an image processing-based refinement indicator. The developed inversion framework allows fast and memory-efficient wave simulation and object identification. While we study the general behavior of the proposed approach on 2D benchmark problems, a final 3D problem shows that it can also be used to identify voids in geometrically complex spatial structures.
△ Less
Submitted 31 May, 2023;
originally announced May 2023.
-
Two-scale analysis and design of spaceframes with complex additive manufactured nodes
Authors:
Oguz Oztoprak,
Alexander Paolini,
Pierluigi D'Acunto,
Ernst Rank,
Stefan Kollmannsberger
Abstract:
The advancements in additive manufacturing (AM) technology have allowed for the production of geometrically complex parts with customizable designs. This versatility benefits large-scale space-frame structures, as the individual design of each structural node can be tailored to meet specific mechanical and other functional requirements. To this end, however, the design and analysis of such space-f…
▽ More
The advancements in additive manufacturing (AM) technology have allowed for the production of geometrically complex parts with customizable designs. This versatility benefits large-scale space-frame structures, as the individual design of each structural node can be tailored to meet specific mechanical and other functional requirements. To this end, however, the design and analysis of such space-frames with distinct structural nodes needs to be highly automated. A critical aspect in this context is automated integration of the local 3D features into the 1D large-scale models. In the present work, a two-scale modeling approach is developed to improve the design and linear-elastic analysis of space frames with complex additively manufactured nodes. The mechanical characteristics of the 3D nodes are numerically reduced through an automated dimensional reduction process based on the Finite Cell Method (FCM) and substructuring. The reduced stiffness quantities are assembled in the large-scale 1D model which, in turn, enables efficient structural analysis. The response of the 1D model is passed on to the local model, enabling fully resolved 3D linear-elastic analysis. The proposed approach is numerically verified on a simplified beam example. Furthermore, the workflow is demonstrated on a tree canopy structure with additively manufactured nodes with bolted connections. The form of the large-scale structure is found based on the Combinatorial Equilibrium Modeling framework, and the different designs of the local structural nodes are based on generative exploration of the design space. It is demonstrated that the proposed methodology effectively automates the design and analysis of space-frame structures with complex, distinct structural nodes.
△ Less
Submitted 17 January, 2023;
originally announced January 2023.
-
Immersed boundary parametrizations for full waveform inversion
Authors:
Tim Bürchner,
Philipp Kopp,
Stefan Kollmannsberger,
Ernst Rank
Abstract:
Full Waveform Inversion (FWI) is a successful and well-established inverse method for reconstructing material models from measured wave signals. In the field of seismic exploration, FWI has proven particularly successful in the reconstruction of smoothly varying material deviations. In contrast, non-destructive testing (NDT) often requires the detection and specification of sharp defects in a spec…
▽ More
Full Waveform Inversion (FWI) is a successful and well-established inverse method for reconstructing material models from measured wave signals. In the field of seismic exploration, FWI has proven particularly successful in the reconstruction of smoothly varying material deviations. In contrast, non-destructive testing (NDT) often requires the detection and specification of sharp defects in a specimen. If the contrast between materials is low, FWI can be successfully applied to these problems as well. However, so far the method is not fully suitable to image defects such as voids, which are characterized by a high contrast in the material parameters. In this paper, we introduce a dimensionless scaling function $γ$ to model voids in the forward and inverse scalar wave equation problem. Depending on which material parameters this function $γ$ scales, different modeling approaches are presented, leading to three formulations of mono-parameter FWI and one formulation of two-parameter FWI. The resulting problems are solved by first-order optimization, where the gradient is computed by an ajdoint state method. The corresponding Fréchet kernels are derived for each approach and the associated minimization is performed using an L-BFGS algorithm. A comparison between the different approaches shows that scaling the density with $γ$ is most promising for parameterizing voids in the forward and inverse problem. Finally, in order to consider arbitrary complex geometries known a priori, this approach is combined with an immersed boundary method, the finite cell method (FCM).
△ Less
Submitted 16 September, 2022;
originally announced September 2022.
-
Quantitative reconstruction of defects in multi-layered bonded composites using fully convolutional network-based ultrasonic inversion
Authors:
Jing Rao,
Fangshu Yang,
Huadong Mo,
Stefan Kollmannsberger,
Ernst Rank
Abstract:
Ultrasonic methods have great potential applications to detect and characterize defects in multi-layered bonded composites. However, it remains challenging to quantitatively reconstruct defects, such as disbonds and kissing bonds, that influence the integrity of adhesive bonds and seriously reduce the strength of assemblies. In this work, an ultrasonic method based on the supervised fully convolut…
▽ More
Ultrasonic methods have great potential applications to detect and characterize defects in multi-layered bonded composites. However, it remains challenging to quantitatively reconstruct defects, such as disbonds and kissing bonds, that influence the integrity of adhesive bonds and seriously reduce the strength of assemblies. In this work, an ultrasonic method based on the supervised fully convolutional network (FCN) is proposed to quantitatively reconstruct defects hidden in multi-layered bonded composites. In the training process of this method, an FCN establishes a non-linear mapping from measured ultrasonic data to the corresponding velocity models of multi-layered bonded composites. In the predicting process, the trained network obtained from the training process is used to directly reconstruct the velocity models from the new measured ultrasonic data of adhesively bonded composites. The presented FCN-based inversion method can automatically extract useful features in multi-layered composites. Although this method is computationally expensive in the training process, the prediction itself in the online phase takes only seconds. The numerical results show that the FCN-based ultrasonic inversion method is capable to accurately reconstruct ultrasonic velocity models of the high contrast defects, which has great potential for online detection of adhesively bonded composites.
△ Less
Submitted 11 September, 2021;
originally announced September 2021.
-
Uncertainty quantification of microstructure variability and mechanical behaviour of additively manufactured lattice structures
Authors:
Nina Korshunova,
Iason Papaioannou,
Stefan Kollmannsberger,
Daninel Straub,
Ernst Rank
Abstract:
Process-induced defects are the leading cause of discrepancies between as-designed and as-manufactured additive manufacturing (AM) product behavior. Especially for metal lattices, the variations in the printed geometry cannot be neglected. Therefore, the evaluation of the influence of microstructural variability on their mechanical behavior is crucial for the quality assessment of the produced str…
▽ More
Process-induced defects are the leading cause of discrepancies between as-designed and as-manufactured additive manufacturing (AM) product behavior. Especially for metal lattices, the variations in the printed geometry cannot be neglected. Therefore, the evaluation of the influence of microstructural variability on their mechanical behavior is crucial for the quality assessment of the produced structures. Commonly, the as-manufactured geometry can be obtained by computed tomography (CT). However, to incorporate all process-induced defects into the numerical analysis is often computationally demanding. Thus, commonly this task is limited to a predefined set of considered variations, such as strut size or strut diameter. In this work, a CT-based binary random field is proposed to generate statistically equivalent geometries of periodic metal lattices. The proposed random field model in combination with the Finite Cell Method (FCM), an immersed boundary method, allows to efficiently evaluate the influence of the underlying microstructure on the variability of the mechanical behavior of AM products. Numerical analysis of two lattices manufactured at different scales shows an excellent agreement with experimental data. Furthermore, it provides a unique insight into the effects of the process on the occurring geometrical variations and final mechanical behavior.
△ Less
Submitted 17 March, 2021;
originally announced March 2021.
-
Bending behavior of additively manufactured lattice structures: numerical characterization and experimental validation
Authors:
Nina Korshunova,
Gianluca Alaimo,
Seyyed Bahram Hosseini,
Massimo Carraturo,
Alessandro Reali,
Jarkko Niiranen,
Ferdinando Auricchio,
Ernst Rank,
Stefan Kollmannsberger
Abstract:
Selective Laser Melting (SLM) technology has undergone significant development in the past years providing unique flexibility for the fabrication of complex metamaterials such as octet-truss lattices. However, the microstructure of the final parts can exhibit significant variations due to the high complexity of the manufacturing process. Consequently, the mechanical behavior of these lattices is s…
▽ More
Selective Laser Melting (SLM) technology has undergone significant development in the past years providing unique flexibility for the fabrication of complex metamaterials such as octet-truss lattices. However, the microstructure of the final parts can exhibit significant variations due to the high complexity of the manufacturing process. Consequently, the mechanical behavior of these lattices is strongly dependent on the process-induced defects, raising the importance on the incorporation of as-manufactured geometries into the computational structural analysis. This, in turn, challenges the traditional mesh-conforming methods making the computational costs prohibitively large. In the present work, an immersed image-to-analysis framework is applied to efficiently evaluate the bending behavior of AM lattices. To this end, we employ the Finite Cell Method (FCM) to perform a three-dimensional numerical analysis of the three-point bending test of a lattice structure and compare the as-designed to as-manufactured effective properties. Furthermore, we undertake a comprehensive study on the applicability of dimensionally reduced beam models to the prediction of the bending behavior of lattice beams and validate classical and strain gradient beam theories applied in combination with the FCM. The numerical findings suggest that the SLM octet-truss lattices exhibit size effects, thus, requiring a flexible framework to incorporate high-order continuum theories.
△ Less
Submitted 22 January, 2021;
originally announced January 2021.
-
Image-based numerical characterization and experimental validation of tensile behavior of octet-truss lattice structures
Authors:
Nina Korshunova,
Gianluca Alaimo,
Seyyed Bahram Hosseini,
Massimo Carraturo,
Alessandro Reali,
Jarkko Niiranen,
Ferdinando Auricchio,
Ernst Rank,
Stefan Kollmannsberger
Abstract:
The production of lightweight metal lattice structures has received much attention due to the recent developments in additive manufacturing (AM). The design flexibility comes, however, with the complexity of the underlying physics. In fact, metal additive manufacturing introduces process-induced geometrical defects that mainly result in deviations of the effective geometry from the nominal one. Th…
▽ More
The production of lightweight metal lattice structures has received much attention due to the recent developments in additive manufacturing (AM). The design flexibility comes, however, with the complexity of the underlying physics. In fact, metal additive manufacturing introduces process-induced geometrical defects that mainly result in deviations of the effective geometry from the nominal one. This change in the final printed shape is the primary cause of the gap between the as-designed and as-manufactured mechanical behavior of AM products. Thus, the possibility to incorporate the precise manufactured geometries into the computational analysis is crucial for the quality and performance assessment of the final parts. Computed tomography (CT) is an accurate method for the acquisition of the manufactured shape. However, it is often not feasible to integrate the CT-based geometrical information into the traditional computational analysis due to the complexity of the meshing procedure for such high-resolution geometrical models and the prohibitive numerical costs. In this work, an embedded numerical framework is applied to efficiently simulate and compare the mechanical behavior of as-designed to as-manufactured octet-truss lattice structures. The parts are produced using laser powder bed fusion (LPBF). Employing an immersed boundary method, namely the Finite Cell Method (FCM), we perform direct numerical simulations (DNS) and first-order numerical homogenization analysis of a tensile test for a 3D printed octet-truss structure. Numerical results based on CT scan (as-manufactured geometry) show an excellent agreement with experimental measurements, whereas both DNS and first-order numerical homogenization performed directly on the 3D virtual model (as-designed geometry) of the structure show a significant deviation from experimental data.
△ Less
Submitted 11 March, 2021; v1 submitted 14 December, 2020;
originally announced December 2020.
-
Finite Cell Method for functionally graded materials based on V-models and homogenized microstructures
Authors:
Benjamin Wassermann,
Nina Korshunova,
Stefan Kollmannsberger,
Ernst Rank,
Gershon Elber
Abstract:
This paper proposes an extension of the finite cell method (FCM) to V-rep models, a novel geometric framework for volumetric representations. This combination of an embedded domain approach (FCM) and a new modeling framework (V-rep) forms the basis for an efficient and accurate simulation of mechanical artifacts, which are not only characterized by complex shapes but also by their non-standard int…
▽ More
This paper proposes an extension of the finite cell method (FCM) to V-rep models, a novel geometric framework for volumetric representations. This combination of an embedded domain approach (FCM) and a new modeling framework (V-rep) forms the basis for an efficient and accurate simulation of mechanical artifacts, which are not only characterized by complex shapes but also by their non-standard interior structure. These objects gain more and more interest in the context of the new design opportunities opened by additive manufacturing, mainly when graded or micro-structured material is applied. Two different types of functionally graded materials (FGM) are considered: The first one, multi-material FGM, is described using the V-rep models' inherent property to assign different properties throughout the interior of a domain. The second, single-material FGM -- which is heterogeneously micro-structured -- characterizes the effective material behavior of representative volume elements by homogenization and performs large-scale simulations using the embedded domain approach.
△ Less
Submitted 20 July, 2020;
originally announced July 2020.
-
A 3D benchmark problem for crack propagation in brittle fracture
Authors:
Lisa Hug,
Stefan Kollmannsberger,
Zohar Yosibash,
Ernst Rank
Abstract:
We propose a full 3D benchmark problem for brittle fracture based on experiments as well as a validation in the context of phase-field models. The example consists of a series of four-point bending tests on graphite specimens with sharp V-notches at different inclination angles. This simple setup leads to a mixed mode (I + II + III) loading which results in complex yet stably reproducible crack su…
▽ More
We propose a full 3D benchmark problem for brittle fracture based on experiments as well as a validation in the context of phase-field models. The example consists of a series of four-point bending tests on graphite specimens with sharp V-notches at different inclination angles. This simple setup leads to a mixed mode (I + II + III) loading which results in complex yet stably reproducible crack surfaces. The proposed problem is well suited for benchmarking numerical methods for brittle fracture and allows for a quantitative comparison of failure loads and propagation paths as well as initiation angles and the fracture surface. For evaluation of the crack surfaces image-based 3D models of the fractured specimen are provided along with experimental and numerical results. In addition, measured failure loads and computed load-displacement curves are given. To demonstrate the applicability of the benchmark problem, we show that for a phase-field model based on the Finite Cell Method and multi-level hp-refinement the complex crack surface as well as the failure loads can be well reproduced.
△ Less
Submitted 15 October, 2019;
originally announced October 2019.
-
Direct structural analysis of domains defined by point clouds
Authors:
László Kudela,
Stefan Kollmannsberger,
Umut Almac,
Ernst Rank
Abstract:
This contribution presents a method that aims at the numerical analysis of solids represented by oriented point clouds. The proposed approach is based on the Finite Cell Method, a high-order immersed boundary technique that computes on a regular background grid of finite elements and requires only inside-outside information from the geometric model. It is shown that oriented point clouds provide s…
▽ More
This contribution presents a method that aims at the numerical analysis of solids represented by oriented point clouds. The proposed approach is based on the Finite Cell Method, a high-order immersed boundary technique that computes on a regular background grid of finite elements and requires only inside-outside information from the geometric model. It is shown that oriented point clouds provide sufficient information for these point-membership classifications. Further, we address a tessellation-free formulation of contour integrals that allows to apply Neumann boundary conditions on point clouds without having to recover the underlying surface. Two-dimensional linear elastic benchmark examples demonstrate that the method is able to provide the same accuracy as those computed with conventional, continuous surface descriptions, because the associated error can be controlled by the density of the cloud. Three-dimensional examples computed on point clouds of historical structures show how the method can be employed to establish seamless connections between digital shape measurement techniques and numerical analyses.
△ Less
Submitted 5 June, 2019;
originally announced June 2019.
-
Integrating CAD and Numerical Analysis: 'Dirty Geometry' handling using the Finite Cell Method
Authors:
Benjamin Wassermann,
Stefan Kollmannsberger,
Shuohui Yin,
László Kudela,
Ernst Rank
Abstract:
This paper proposes a computational methodology for the integration of Computer Aided Design (CAD) and the Finite Cell Method (FCM) for models with "dirty geometries". FCM, being a fictitious domain approach based on higher order finite elements, embeds the physical model into a fictitious domain, which can be discretized without having to take into account the boundary of the physical domain. The…
▽ More
This paper proposes a computational methodology for the integration of Computer Aided Design (CAD) and the Finite Cell Method (FCM) for models with "dirty geometries". FCM, being a fictitious domain approach based on higher order finite elements, embeds the physical model into a fictitious domain, which can be discretized without having to take into account the boundary of the physical domain. The true geometry is captured by a precise numerical integration of elements cut by the boundary. Thus, an effective Point Membership Classification algorithm that determines the inside-outside state of an integration point with respect to the physical domain is a core operation in FCM. To treat also "dirty geometries", i.e. imprecise or flawed geometric models, a combination of a segment-triangle intersection algorithm and a flood fill algorithm being insensitive to most CAD model flaws is proposed to identify the affiliation of the integration points. The present method thus allows direct computations on geometrically and topologically flawed models. The potential and merit for practical applications of the proposed method is demonstrated by several numerical examples.
△ Less
Submitted 30 April, 2019; v1 submitted 5 October, 2018;
originally announced October 2018.
-
From geometric design to numerical analysis: A direct approach using the Finite Cell Method on Constructive Solid Geometry
Authors:
Benjamin Wassermann,
Stefan Kollmannsberger,
Tino Bog,
Ernst Rank
Abstract:
During the last ten years, increasing efforts were made to improve and simplify the process from Computer Aided Design (CAD) modeling to a numerical simulation. It has been shown that the transition from one model to another, i.e. the meshing, is a bottle-neck. Several approaches have been developed to overcome this time-consuming step, e.g. Isogeometric Analysis (IGA), which applies the shape fun…
▽ More
During the last ten years, increasing efforts were made to improve and simplify the process from Computer Aided Design (CAD) modeling to a numerical simulation. It has been shown that the transition from one model to another, i.e. the meshing, is a bottle-neck. Several approaches have been developed to overcome this time-consuming step, e.g. Isogeometric Analysis (IGA), which applies the shape functions used for the geometry description (typically B-Splines and NURBS) directly to the numerical analysis. In contrast to IGA, which deals with boundary represented models (B-Rep), our approach focuses on parametric volumetric models such as Constructive Solid Geometries (CSG). These models have several advantages, as their geometry description is inherently watertight and they provide a description of the model's interior. To be able to use the explicit mathematical description of these models, we employ the Finite Cell Method (FCM). Herein, the only necessary input is a reliable statement whether an (integration-) point lies inside or outside of the geometric model. This paper mainly discusses such point-in-membership tests on various geometric objects like sweeps and lofts, as well as several geometric operations such as filleting or chamfering. We demonstrate that, based on the information of the construction method of these objects, the point-in-membership-test can be carried out efficiently and robustly
△ Less
Submitted 28 September, 2018;
originally announced September 2018.
-
A design-through-analysis approach using the Finite Cell Method (ECCOMAS Congress 2016)
Authors:
Benjamin Wassermann,
Tino Bog,
Stefan Kollmannsberger,
Ernst Rank
Abstract:
Modern 3D Computer-Aided-Design (CAD) systems use mainly two types of geometric models. Classically, objects are defined by a Boundary Representation (B-Rep), where only the objects' surfaces with their corresponding edges and nodes are stored. One disadvantage concerning a numerical simulation is that B-Rep models are not necessarily watertight. These 'dirty geometries' cause major difficulties i…
▽ More
Modern 3D Computer-Aided-Design (CAD) systems use mainly two types of geometric models. Classically, objects are defined by a Boundary Representation (B-Rep), where only the objects' surfaces with their corresponding edges and nodes are stored. One disadvantage concerning a numerical simulation is that B-Rep models are not necessarily watertight. These 'dirty geometries' cause major difficulties in computational analysis because even basic geometric operations such as point-in-membership tests fail, not to mention meshing as required by classical boundary conforming finite element methods. Alternatively, objects may be represented by Constructive Solid Geometry (CSG), which is strongly related to Procedural Modeling (PM). In this context, the model is created using Boolean operations on primitives. The modeling process is then either stored as a sequence (PM), or as a construction tree (CSG). In contrast to B-Rep models, CSG models are intrinsically water-tight. To run a finite element simulation on a water-tight CSG model, two alternatives are possible: (i) it can either be converted to a B-Rep-model to obtain a finite element mesh or (ii) its implicit description can be used directly by applying an embedded domain approach, like the Finite Cell Method (FCM). In this contribution, we present a design-through analysis methodology using CSG and FCM. A crucial point in FCM is a fast and reliable point-in-membership test which can be directly derived from the CSG model. We present the outline of the modeling approach, the realization of the point-in-membership test as a sequence of CSG-operations, and discuss advantages and limitations on complex models of relevance in mechanical engineering.
△ Less
Submitted 26 September, 2018;
originally announced September 2018.
-
Residual stresses in metal deposition modeling: discretizations of higher order
Authors:
Ali Özcan,
Stefan Kollmannsberger,
John N. Jomo,
Ernst Rank
Abstract:
This article addresses the research question if and how the finite cell method, an embedded domain finite element method of high order, may be used in the simulation of metal deposition to harvest its computational efficiency. This application demands for the solution of a coupled thermo-elasto-plastic problem on transient meshes within which history variables need to be managed dynamically on non…
▽ More
This article addresses the research question if and how the finite cell method, an embedded domain finite element method of high order, may be used in the simulation of metal deposition to harvest its computational efficiency. This application demands for the solution of a coupled thermo-elasto-plastic problem on transient meshes within which history variables need to be managed dynamically on non-boundary conforming discretizations. To this end, we propose to combine the multi-level hp-method and the finite cell method. The former was specifically designed to treat high-order finite element discretizations on transient meshes, while the latter offers a remedy to retain high-order convergence rates also in cases where the physical boundary does not coincide with the boundary of the discretization. We investigate the performance of the method at two analytical and one experimental benchmark.
△ Less
Submitted 25 September, 2018; v1 submitted 20 September, 2018;
originally announced September 2018.
-
Robust and parallel scalable iterative solutions for large-scale finite cell analyses
Authors:
John N. Jomo,
Frits de Prenter,
Mohamed Elhaddad,
Davide D'Angella,
Clemens V. Verhoosel,
Stefan Kollmannsberger,
Jan S. Kirschke,
Vera Nübel,
Harald van Brummelen,
Ernst Rank
Abstract:
The finite cell method is a highly flexible discretization technique for numerical analysis on domains with complex geometries. By using a non-boundary conforming computational domain that can be easily meshed, automatized computations on a wide range of geometrical models can be performed. Application of the finite cell method, and other immersed methods, to large real-life and industrial problem…
▽ More
The finite cell method is a highly flexible discretization technique for numerical analysis on domains with complex geometries. By using a non-boundary conforming computational domain that can be easily meshed, automatized computations on a wide range of geometrical models can be performed. Application of the finite cell method, and other immersed methods, to large real-life and industrial problems is often limited due to the conditioning problems associated with these methods. These conditioning problems have caused researchers to resort to direct solution methods, which signifi- cantly limit the maximum size of solvable systems. Iterative solvers are better suited for large-scale computations than their direct counterparts due to their lower memory requirements and suitability for parallel computing. These benefits can, however, only be exploited when systems are properly conditioned. In this contribution we present an Additive-Schwarz type preconditioner that enables efficient and parallel scalable iterative solutions of large-scale multi-level hp-refined finite cell analyses.
△ Less
Submitted 5 September, 2018; v1 submitted 4 September, 2018;
originally announced September 2018.
-
Parallelization of the multi-level hp-adaptive finite cell method
Authors:
John N. Jomo,
Nils Zander,
Mohamed Elhaddad,
Ali Özcan,
Stefan Kollmannsberger,
Ralf-Peter Mundani,
Ernst Rank
Abstract:
The multi-level hp-refinement scheme is a powerful extension of the finite element method that allows local mesh adaptation without the trouble of constraining hanging nodes. This is achieved through hierarchical high-order overlay meshes, a hp-scheme based on spatial refinement by superposition. An efficient parallelization of this method using standard domain decomposition approaches in combinat…
▽ More
The multi-level hp-refinement scheme is a powerful extension of the finite element method that allows local mesh adaptation without the trouble of constraining hanging nodes. This is achieved through hierarchical high-order overlay meshes, a hp-scheme based on spatial refinement by superposition. An efficient parallelization of this method using standard domain decomposition approaches in combination with ghost elements faces the challenge of a large basis function support resulting from the overlay structure and is in many cases not feasible. In this contribution, a parallelization strategy for the multi-level hp-scheme is presented that is adapted to the scheme's simple hierarchical structure. By distributing the computational domain among processes on the granularity of the active leaf elements and utilizing shared mesh data structures, good parallel performance is achieved, as redundant computations on ghost elements are avoided. We show the scheme's parallel scalability for problems with a few hundred elements per process. Furthermore, the scheme is used in conjunction with the finite cell method to perform numerical simulations on domains of complex shape.
△ Less
Submitted 3 July, 2018;
originally announced July 2018.
-
A high-performance interactive computing framework for engineering applications
Authors:
Jovana Knežević,
Ralf-Peter Mundani,
Ernst Rank
Abstract:
To harness the potential of advanced computing technologies, efficient (real time) analysis of large amounts of data is as essential as are front-line simulations. In order to optimise this process, experts need to be supported by appropriate tools that allow to interactively guide both the computation and data exploration of the underlying simulation code. The main challenge is to seamlessly feed…
▽ More
To harness the potential of advanced computing technologies, efficient (real time) analysis of large amounts of data is as essential as are front-line simulations. In order to optimise this process, experts need to be supported by appropriate tools that allow to interactively guide both the computation and data exploration of the underlying simulation code. The main challenge is to seamlessly feed the user requirements back into the simulation. State-of-the-art attempts to achieve this, have resulted in the insertion of so-called check- and break-points at fixed places in the code. Depending on the size of the problem, this can still compromise the benefits of such an attempt, thus, preventing the experience of real interactive computing. To leverage the concept for a broader scope of applications, it is essential that a user receives an immediate response from the simulation to his or her changes. Our generic integration framework, targeted to the needs of the computational engineering domain, supports distributed computations as well as on-the-fly visualisation in order to reduce latency and enable a high degree of interactivity with only minor code modifications. Namely, the regular course of the simulation coupled to our framework is interrupted in small, cyclic intervals followed by a check for updates. When new data is received, the simulation restarts automatically with the updated settings (boundary conditions, simulation parameters, etc.). To obtain rapid, albeit approximate feedback from the simulation in case of perpetual user interaction, a multi-hierarchical approach is advantageous. Within several different engineering test cases, we will demonstrate the flexibility and the effectiveness of our approach.
△ Less
Submitted 2 July, 2018;
originally announced July 2018.
-
Dynamic load balancing strategies for hierarchical p-FEM solvers
Authors:
Ralf-Peter Mundani,
Alexander Düster,
Jovana Knežević,
Andreas Niggl,
Ernst Rank
Abstract:
Equation systems resulting from a p-version FEM discretisation typically require a special treatment as iterative solvers are not very efficient here. Applying hierarchical concepts based on a nested dissection approach allow for both the design of sophisticated solvers as well as for advanced parallelisation strategies. To fully exploit the underlying computing power of parallel systems, dynamic…
▽ More
Equation systems resulting from a p-version FEM discretisation typically require a special treatment as iterative solvers are not very efficient here. Applying hierarchical concepts based on a nested dissection approach allow for both the design of sophisticated solvers as well as for advanced parallelisation strategies. To fully exploit the underlying computing power of parallel systems, dynamic load balancing strategies become an essential component.
△ Less
Submitted 2 July, 2018;
originally announced July 2018.
-
Computational steering of complex flow simulations
Authors:
Atanas Atanasov,
Hans-Joachim Bungartz,
Jérôme Frisch,
Miriam Mehl,
Ralf-Peter Mundani,
Ernst Rank,
Christoph van Treeck
Abstract:
Computational Steering, the combination of a simulation back-end with a visualisation front-end, offers great possibilities to exploit and optimise scenarios in engineering applications. Due to its interactivity, it requires fast grid generation, simulation, and visualisation and, therefore, mostly has to rely on coarse and inaccurate simulations typically performed on rather small interactive com…
▽ More
Computational Steering, the combination of a simulation back-end with a visualisation front-end, offers great possibilities to exploit and optimise scenarios in engineering applications. Due to its interactivity, it requires fast grid generation, simulation, and visualisation and, therefore, mostly has to rely on coarse and inaccurate simulations typically performed on rather small interactive computing facilities and not on much more powerful high-performance computing architectures operated in batch-mode. This paper presents a steering environment that intends to bring these two worlds - the interactive and the classical HPC world - together in an integrated way. The environment consists of efficient fluid dynamics simulation codes and a steering and visualisation framework providing a user interface, communication methods for distributed steering, and parallel visualisation tools. The gap between steering and HPC is bridged by a hierarchical approach that performs fast interactive simulations for many scenario variants increasing the accuracy via hierarchical refinements in dependence of the time the user wants to wait. Finally, the user can trigger large simulations for selected setups on an HPC architecture exploiting the pre-computations already done on the interactive system.
△ Less
Submitted 2 July, 2018;
originally announced July 2018.
-
Framework for the hybrid parallelisation of simulation codes
Authors:
Ralf-Peter Mundani,
Marko Ljucović,
Ernst Rank
Abstract:
Writing efficient hybrid parallel code is tedious, error-prone, and requires good knowledge of both parallel programming and multithreading such as MPI and OpenMP, resp. Therefore, we present a framework which is based on a job model that allows the user to incorporate his sequential code with manageable effort and code modifications in order to be executed in parallel on clusters or supercomputer…
▽ More
Writing efficient hybrid parallel code is tedious, error-prone, and requires good knowledge of both parallel programming and multithreading such as MPI and OpenMP, resp. Therefore, we present a framework which is based on a job model that allows the user to incorporate his sequential code with manageable effort and code modifications in order to be executed in parallel on clusters or supercomputers built from modern multi-core CPUs. The primary application domain of this framework are simulation codes from engineering disciplines as those are in many cases still sequential and due to their memory and runtime demands prominent candidates for parallelisation.
△ Less
Submitted 1 July, 2018;
originally announced July 2018.
-
Ensuring domain consistency in an adaptive framework with distributed topology for fluid flow simulations
Authors:
Christoph Ertl,
Ralf-Peter Mundani,
Ernst Rank
Abstract:
Top-tier parallel computing clusters continue to accumulate more and more computational power with more and better CPUs and Networks. This allows, especially for environmental simulations, computations with larger domain sizes and better resolution. One of the challenges becoming increasingly important is the decomposition and distribution of the overall work load. State-of-the-art parallel codes…
▽ More
Top-tier parallel computing clusters continue to accumulate more and more computational power with more and better CPUs and Networks. This allows, especially for environmental simulations, computations with larger domain sizes and better resolution. One of the challenges becoming increasingly important is the decomposition and distribution of the overall work load. State-of-the-art parallel codes usually use solutions that involve complete knowledge of the domain topology, which will lead to communication and memory bottlenecks when computing very large domains. To meet this challenge, the authors propose a new strategy for decentralised domain management, based on a proven hierarchic data structure. On the way of developing a framework where individual sub-domains only have local knowledge of their surroundings, this contribution describes the communication patterns used in ensuring a consistent domain, without the need for expensive global broadcast messages. Furthermore, the routines necessary to deal with adaptive changes in domain topology, due to refinement, coarsening and migration of sub-domains to different computational resources, are discussed in detail.
△ Less
Submitted 1 July, 2018;
originally announced July 2018.
-
Interactive data exploration for high-performance fluid flow computations through porous media
Authors:
Nevena Perović,
Jérôme Frisch,
Ralf-Peter Mundani,
Ernst Rank
Abstract:
Huge data advent in high-performance computing (HPC) applications such as fluid flow simulations usually hinders the interactive processing and exploration of simulation results. Such an interactive data exploration not only allows scientiest to 'play' with their data but also to visualise huge (distributed) data sets in both an efficient and easy way. Therefore, we propose an HPC data exploration…
▽ More
Huge data advent in high-performance computing (HPC) applications such as fluid flow simulations usually hinders the interactive processing and exploration of simulation results. Such an interactive data exploration not only allows scientiest to 'play' with their data but also to visualise huge (distributed) data sets in both an efficient and easy way. Therefore, we propose an HPC data exploration service based on a sliding window concept, that enables researches to access remote data (available on a supercomputer or cluster) during simulation runtime without exceeding any bandwidth limitations between the HPC back-end and the user front-end.
△ Less
Submitted 30 June, 2018;
originally announced July 2018.
-
A sliding window technique for interactive high-performance computing scenarios
Authors:
Ralf-Peter Mundani,
Jérôme Frisch,
Vasco Varduhn,
Ernst Rank
Abstract:
Interactive high-performance computing is doubtlessly beneficial for many computational science and engineering applications whenever simulation results should be visually processed in real time, i.e. during the computation process. Nevertheless, interactive HPC entails a lot of new challenges that have to be solved - one of them addressing the fast and efficient data transfer between a simulation…
▽ More
Interactive high-performance computing is doubtlessly beneficial for many computational science and engineering applications whenever simulation results should be visually processed in real time, i.e. during the computation process. Nevertheless, interactive HPC entails a lot of new challenges that have to be solved - one of them addressing the fast and efficient data transfer between a simulation back end and visualisation front end, as several gigabytes of data per second are nothing unusual for a simulation running on some (hundred) thousand cores. Here, a new approach based on a sliding window technique is introduced that copes with any bandwidth limitations and allows users to study both large and small scale effects of the simulation results in an interactive fashion.
△ Less
Submitted 29 June, 2018;
originally announced July 2018.
-
Phase-field modeling of brittle fracture with multi-level hp-FEM and the finite cell method
Authors:
Sindhu Nagaraja,
Mohamed Elhaddad,
Marreddy Ambati,
Stefan Kollmannsberger,
Laura De Lorenzis,
Ernst Rank
Abstract:
The difficulties in dealing with discontinuities related to a sharp crack are overcome in the phase-field approach for fracture by modeling the crack as a diffusive object being described by a continuous field having high gradients. The discrete crack limit case is approached for a small length-scale parameter that controls the width of the transition region between the fully broken and the undama…
▽ More
The difficulties in dealing with discontinuities related to a sharp crack are overcome in the phase-field approach for fracture by modeling the crack as a diffusive object being described by a continuous field having high gradients. The discrete crack limit case is approached for a small length-scale parameter that controls the width of the transition region between the fully broken and the undamaged phases. From a computational standpoint, this necessitates fine meshes, at least locally, in order to accurately resolve the phase-field profile. In the classical approach, phase-field models are computed on a fixed mesh that is a priori refined in the areas where the crack is expected to propagate. This on the other hand curbs the convenience of using phase-field models for unknown crack paths and its ability to handle complex crack propagation patterns. In this work, we overcome this issue by employing the multi-level hp-refinement technique that enables a dynamically changing mesh which in turn allows the refinement to remain local at singularities and high gradients without problems of hanging nodes. Yet, in case of complex geometries, mesh generation and in particular local refinement becomes non-trivial. We address this issue by integrating a two-dimensional phase-field framework for brittle fracture with the finite cell method (FCM). The FCM based on high-order finite elements is a non-geometry-conforming discretization technique wherein the physical domain is embedded into a larger fictitious domain of simple geometry that can be easily discretized. This facilitates mesh generation for complex geometries and supports local refinement. Numerical examples including a comparison to a validation experiment illustrate the applicability of the multi-level hp-refinement and the FCM in the context of phase-field simulations.
△ Less
Submitted 23 April, 2018;
originally announced April 2018.