Documents published in Scipedia

• PSYDAC: a high-performance IGA library in Python

  Y. Güçlü, S. Hadjout, A. Ratnani

  eccomas2022.

  
  

  Abstract

  Psydac is a Python 3 library for the solution of partial differential equations, which combines the convenience of a domain specific language with the speed of a high-performance parallel engine. Its main focus is on isogeometric analysis using tensor-product B-spline finite elements; to this end it uses an optimized sparse format called 'stencil matrix', which drastically reduces memory storage compared to the popular CSR/CSC formats. It supports multi-patch mapped geometries, and finite element exterior calculus. It can distribute each domain patch across many MPI processes, with multiple OpenMP threads operating in each block. The users of Psydac define a weak form of the model equations through SymPDE, an extension of Sympy that provides the mathematical expressions and checks their semantic validity. Simple mappings can be defined analytically, and multi-patch NURBS geometries can be imported from file. Once a finite element discretization is chosen, Psydac maps abstract concepts onto concrete objects, the basic building blocks being MPI-distributed vectors and matrices. Python code is automatically generated for the model-specific operations, namely matrix and vector assembly, and user-defined diagnostics. Finally, Psydac accelerates all computationally intensive operations using Pyccel, a transpiler which converts Python code to either C or Fortran. We present the library design, the typical usage workflow, the user interface for a simple 2D example, and the parallel scaling results in a large 3D simulation. In addition we show a few complex applications in fluid dynamics and electromagnetism, where the accuracy of the solver is verified against manufactured and reference solutions.

  Abstract
  Psydac is a Python 3 library for the solution of partial differential equations, which combines the convenience of a domain specific language with the speed of a high-performance [...]

  • 21
  •  read

• Advanced Experiments on Gaussian Process-based Multi-fidelity Methods over Diverse Mathematical Characteristics

  S. Yildiz, H. Pehlivan-Solak, M. Diez, O. Goren, M. Nikbay

  eccomas2022.

  
  

  Abstract

  Advanced applications of multi-fidelity surrogate modelling techniques provide significant improvements in optimization and uncertainty quantification studies in many engineering fields. Multi-fidelity surrogate modelling can efficiently save the design process from the computational time burden caused by the need for numerous computationally expensive simulations. However, no consensus exists about which multi-fidelity surrogate modelling technique usually exhibits superiority over the other methods given for certain conditions. Therefore, the present paper focuses on assessing the performances of the Gaussian Process-based multi-fidelity methods across selected benchmark problems, especially chosen to capture diverse mathematical characteristics, by experimenting with their learning processes concerning different performance criteria. In this study, a comparison of Linear-Autoregressive Gaussian Process and NonlinearAutoregressive Gaussian Process methods is presented by using benchmark problems that mimic the behaviour of real engineering problems such as localized behaviours, multi-modality, noise, discontinuous response, and different discrepancy types. Our results indicate that the considered methodologies were able to capture the behaviour of the actual function sufficiently within the limited amount of budget for 1-D cases. As the problem dimension increases, the required number of training data increases exponentially to construct an acceptable surrogate model. Especially in higher dimensions, i.e. more than 5-D, local error metrics reveal that more training data is needed to attain an efficient surrogate for Gaussian Process based strategies.

  Abstract
  Advanced applications of multi-fidelity surrogate modelling techniques provide significant improvements in optimization and uncertainty quantification studies in many engineering [...]

  • 8
  •  read

• Simulation of wind induced excitation of a membrane structure with ponding water

  N. Kodunthirappully Narayanan, R. Wüchner, J. Degroote

  eccomas2022.

  
  

  Abstract

  This paper proposes a new partitioned coupling approach to simulate the wind induced excitation of a membrane structure with ponding water. This approach uses three different solvers to simulate wind, water and membrane structure. The main assumption here is that the interaction between the wind and water can be neglected due to the small depth and small fetch of the water, relative to the size of the membrane structure. This assumption results in a coupling strategy where the structural solver independently interacts with the wind and water solver. The results from this method is compared with a straightforward approach, where a two-phase solver, modeling the wind and water, is coupled to a structural solver. The obtained results agreed very well with the reference modeling approach, where all the interactions are taken into account. Furthermore, the proposed method was found to be computationally more efficient.

  Abstract
  This paper proposes a new partitioned coupling approach to simulate the wind induced excitation of a membrane structure with ponding water. This approach uses three different [...]

  • 5
  •  read

• Isogeometric analysis of industry applications with Geomiso SEA: a new hybrid software for shell analysis

  P. Karakitsios, V. Tsotoulidi, P. Kolios, G. Mprellas

  eccomas2022.

  
  

  Abstract

  In this paper the new Geomiso SEA software (www.geomiso.com) is proposed for applications on inelastic static isogeometric analysis with shell elements. This hybrid program is applicable to real world structures, while it satisfies the rising need for technical software of dual CAD/CAA nature. It is based on the new isogeometric method, which has attracted a lot of attention for solving boundary value problems. T-spline-based isogeometric shell analysis efficiently handles multi-patch geometries. Geomiso SEA is not just a plug-in, but a complete software solution, used to simulate spline models of structures or machine components, for analyzing their strength and behavior. The utilization of the exact mesh for analysis vanishes geometric errors, while there is no need of repeating the geometry design for refinement purposes. Industry applications on both thick and thin shells are demonstrated with a comparison between Geomiso SEA and FEA programs. This unique solution for seamless integration of the industrial design of shell geometries with its computational realtime testing, appears to be preferable to FEA programs, representing major improvements, such as higher accuracy, and considerably reduced computational cost.

  Abstract
  In this paper the new Geomiso SEA software (www.geomiso.com) is proposed for applications on inelastic static isogeometric analysis with shell elements. This hybrid program [...]

  • 5
  •  read

• Efficient coarse correction for parallel time stepping in plaque growth simulations

  S. Frei, A. Heinlein

  eccomas2022.

  
  

  Abstract

  In order to make the numerical simulation of atherosclerotic plaque growth feasible, a temporal homogenization approach is employed. The resulting macro-scale problem for the plaque growth can be further accelerated by using parallel time integration schemes, such as the parareal algorithm. However, the parallel scalability is dominated by the computational cost of the coarse propagator. Therefore, in this paper, an interpolation-based coarse propagator, which uses growth values from previously computed micro-scale problems, is introduced. For a simple model problem, it is shown that this approach reduces both the computational work for a single parareal iteration as well as the required number of parareal iterations.

  Abstract
  In order to make the numerical simulation of atherosclerotic plaque growth feasible, a temporal homogenization approach is employed. The resulting macro-scale problem for [...]

  • 3
  •  read

• Multi-fidelity active learning for shape optimization problems affected by noise

  J. Wackers, R. Pellegrini, M. Diez, A. Serani, M. Visonneau

  eccomas2022.

  
  

  Abstract

  The efficiency of simulation-driven design optimization based on surrogate models, depends strongly on the suitability of the surrogate model for the simulation data on which it is based. We investigate adaptive surrogate modelling methods that maximize the efficiency and the robustness for any optimization problem. Specific techniques include: adaptive sampling, noise filtering by metamodel tuning, and small initial datasets to give maximum freedom to the adaptation. These methodological advancements are demonstrated for an analytical test problem, as well as the shape optimization of the DTMB 5415 ship model for calm-water resistance.

  Abstract
  The efficiency of simulation-driven design optimization based on surrogate models, depends strongly on the suitability of the surrogate model for the simulation data on which [...]

  • 1
  •  read

• Reduced-order model for large amplitude vibrations of flexible structures coupled with a flow

  T. Flament, J. Deü, A. Placzek, M. Balmaseda, D. Tran

  eccomas2022.

  
  

  Abstract

  This work concerns the numerical modeling of the vibrations of geometrically nonlinear structures coupled with a fluid flow. Firstly, a reduced-order model (ROM) for the geometrically nonlinear structure is proposed. Then, the aforementioned ROM is used to replace a Finite Element solver (FE) in the frame of a fluid-structure partitioned coupling on a two-dimensional example involving vortex induced vibrations.

  Abstract
  This work concerns the numerical modeling of the vibrations of geometrically nonlinear structures coupled with a fluid flow. Firstly, a reduced-order model (ROM) for the geometrically [...]

  • 4
  •  read

• Fluid-structure interaction simulation of a wire-wrapped tube array using overset grids

  H. Dolfen, D. Van Hauwermeiren, A. Bral, J. Degroote

  eccomas2022.

  
  

  Abstract

  Axial flows in tube bundle geometries are omnipresent in nuclear power plants, including fuel assemblies and heat exchangers. The tubes are often long and slender which makes them susceptible to flow-induced vibrations. Current reactor research investigates the use of wires helicoidally wound around each rod, to preserve their mutual distance. This poses new challenges for numerical simulations. The wire-wrapped geometry is complex, among other things due to small gaps, making the meshing process complicated and cumbersome. This research explores the use of overset grids for such wire-wrapped tube bundles. This meshing method allows grids to overlap and the flow solution is interpolated between the different grids, providing more freedom for meshing. In this case the approach consists of making one grid for a single tube with a wire, termed the component mesh, and placing this grid in a so-called background mesh that takes into account other geometrical features (e.g. the bundle duct). In the case of bundles, the same tube mesh can be repeated several times without additional meshing effort, regardless of the bundle size. The approach is verified using 4 cases found in literature, using the same component mesh either one or multiple times for each case, thus reducing the meshing effort to a minimum thanks to the freedom and re-usability overset offers. The first two cases involve a fluid-structure interaction simulation of a single wire-wrapped tube in a cylindrical domain, the first one simulating a steady deformation and the second one a free vibration of the tube. Good agreement with results in literature was found. The latter two cases are Computational Fluid Dynamics simulations of a 7-pin and 19-pin bundle in a hexagonal duct, and excellent agreement with literature results was obtained. The overset approach was proven beneficial for simulating wire-wrapped bundle geometries: with largely reduced effort the same predictive capabilities can be obtained and potentially even extended.

  Abstract
  Axial flows in tube bundle geometries are omnipresent in nuclear power plants, including fuel assemblies and heat exchangers. The tubes are often long and slender which makes [...]

  • 7
  •  read

• Fuzzy logic based rapid visual screening methodology for structural damage state determination of URM buildings

  N. Bektaş

  eccomas2022.

  
  

  Abstract

  Most of the Unreinforced Masonry (URM) buildings are quite old in Europe based on 'Building stock inventory to assess seismic vulnerability across Europe' [1] report. Following the earthquakes (Albania, Italy, etc.) that occurred in Europe, it was revealed that masonry buildings are extremely vulnerable. While probabilistic and deterministic approaches are important for examining a small number of buildings, they do not offer the opportunity to examine a large building stock in a short period of time. Rapid Visual Screening (RVS) methods are used to identify building preand post-earthquake vulnerability. Several RVS techniques have been presented in literature over last 30 years. Recent earthquakes have highlighted critical necessity of a rapid vulnerability assessment method for pre-earthquake warning, mitigation, preparedness, and post-earthquake damage state assessment of existing buildings. These findings demonstrate the importance of using an accurate RVS technique to inspect buildings. Due to the subjectivity of the screener, these RVS methods contain uncertainty and vagueness. Fuzzy Inference System (FIS) overcomes nonrandom uncertainty and vagueness by considering building characteristics in terms of their degree of truth. This paper introduces a FIS-based SRVS case implementation and compares FIS-based Soft-RVS (S-RVS) to traditional RVS methods for identifying building damage state taking into account rapid visual assessment reports about damage caused by the 2019 Albania earthquake. To determine the damage states of URM buildings, 40 buildings damaged in the 2019 Albania earthquake were analyzed and processed to use in the applied fuzzy logic mathematical model. Initial findings demonstrate that the site-specific FIS-based S-RVS method is capable of accurately determining the damage states of at least half of the buildings.

  Abstract
  Most of the Unreinforced Masonry (URM) buildings are quite old in Europe based on 'Building stock inventory to assess seismic vulnerability across Europe' [1] report. [...]

  • 25
  •  read

• Space-time fluid-structure interaction: formulation and dG(0) time discretization

  D. Khimin, J. Roth, T. Wick

  eccomas2022.

  
  

  Abstract

  This contribution is related to the key-note lecture `Space-time fluid-structure interaction with adjoint-based methods for error estimation and optimization' given in the Minisymposium `Innovative Methods for Fluid-Structure Interaction'. The main objective is twofold. First, we design function spaces and a space-time variational-monolithic formulation of fluid-structure interaction in arbitrary Lagrangian-Eulerian coordinates. Second, we apply a Galerkin-time discretization using discontinuous finite elements of degree r = 0. Therein, the main emphasis is on the correct derivation of the jump terms and the integration of nonlinear time derivatives, as the latter arise due to the arbitrary Lagrangian-Eulerian transformation.

  Abstract
  This contribution is related to the key-note lecture `Space-time fluid-structure interaction with adjoint-based methods for error estimation and optimization' given in [...]

  • 16
  •  read