TY - JOUR A1 - Colombo, Daniele A1 - Drira, Slah A1 - Frotscher, Ralf A1 - Staat, Manfred T1 - An element-based formulation for ES-FEM and FS-FEM models for implementation in standard solid mechanics finite element codes for 2D and 3D static analysis JF - International Journal for Numerical Methods in Engineering N2 - Edge-based and face-based smoothed finite element methods (ES-FEM and FS-FEM, respectively) are modified versions of the finite element method allowing to achieve more accurate results and to reduce sensitivity to mesh distortion, at least for linear elements. These properties make the two methods very attractive. However, their implementation in a standard finite element code is nontrivial because it requires heavy and extensive modifications to the code architecture. In this article, we present an element-based formulation of ES-FEM and FS-FEM methods allowing to implement the two methods in a standard finite element code with no modifications to its architecture. Moreover, the element-based formulation permits to easily manage any type of element, especially in 3D models where, to the best of the authors' knowledge, only tetrahedral elements are used in FS-FEM applications found in the literature. Shape functions for non-simplex 3D elements are proposed in order to apply FS-FEM to any standard finite element. KW - distorted element KW - ES-FEM KW - FS-FEM KW - non-simplex S-FEM elements KW - S-FEM Y1 - 2022 U6 - https://doi.org/10.1002/nme.7126 SN - 1097-0207 VL - 124 IS - 2 SP - 402 EP - 433 PB - Wiley CY - Chichester ER - TY - JOUR A1 - Morandi, Paolo A1 - Butenweg, Christoph A1 - Breis, Khaled A1 - Beyer, Katrin A1 - Magenes, Guido ED - Ansal, Atilla T1 - Latest findings on the behaviour factor q for the seismic design of URM buildings JF - Bulletin of Earthquake Engineering N2 - Recent earthquakes as the 2012 Emilia earthquake sequence showed that recently built unreinforced masonry (URM) buildings behaved much better than expected and sustained, despite the maximum PGA values ranged between 0.20–0.30 g, either minor damage or structural damage that is deemed repairable. Especially low-rise residential and commercial masonry buildings with a code-conforming seismic design and detailing behaved in general very well without substantial damages. The low damage grades of modern masonry buildings that was observed during this earthquake series highlighted again that codified design procedures based on linear analysis can be rather conservative. Although advances in simulation tools make nonlinear calculation methods more readily accessible to designers, linear analyses will still be the standard design method for years to come. The present paper aims to improve the linear seismic design method by providing a proper definition of the q-factor of URM buildings. These q-factors are derived for low-rise URM buildings with rigid diaphragms which represent recent construction practise in low to moderate seismic areas of Italy and Germany. The behaviour factor components for deformation and energy dissipation capacity and for overstrength due to the redistribution of forces are derived by means of pushover analyses. Furthermore, considerations on the behaviour factor component due to other sources of overstrength in masonry buildings are presented. As a result of the investigations, rationally based values of the behaviour factor q to be used in linear analyses in the range of 2.0–3.0 are proposed. KW - Unreinforced masonry buildings KW - Modern constructions KW - Seismic design KW - Linear elastic analysis KW - Behaviour factor q Y1 - 2022 U6 - https://doi.org/10.1007/s10518-022-01419-7 SN - 1573-1456 SN - 1570-761X VL - 20 IS - 11 SP - 5797 EP - 5848 PB - Springer Nature CY - Cham ER - TY - JOUR A1 - Peere, Wouter A1 - Blanke, Tobias ED - Vernon, Chris T1 - GHEtool: An open-source tool for borefield sizing in Python JF - Journal of Open Source Software N2 - GHEtool is a Python package that contains all the functionalities needed to deal with borefield design. It is developed for both researchers and practitioners. The core of this package is the automated sizing of borefield under different conditions. The sizing of a borefield is typically slow due to the high complexity of the mathematical background. Because this tool has a lot of precalculated data, GHEtool can size a borefield in the order of tenths of milliseconds. This sizing typically takes the order of minutes. Therefore, this tool is suited for being implemented in typical workflows where iterations are required. GHEtool also comes with a graphical user interface (GUI). This GUI is prebuilt as an exe-file because this provides access to all the functionalities without coding. A setup to install the GUI at the user-defined place is also implemented and available at: https://www.mech.kuleuven.be/en/tme/research/thermal_systems/tools/ghetool. KW - geothermal KW - energy KW - borefields KW - sizing Y1 - 2022 U6 - https://doi.org/10.21105/joss.04406 SN - 2475-9066 VL - 7 IS - 76 SP - 1 EP - 4, 4406 ER - TY - JOUR A1 - Schwager, Christian A1 - Flesch, Robert A1 - Schwarzbözl, Peter A1 - Herrmann, Ulf A1 - Teixeira Boura, Cristiano José T1 - Advanced two phase flow model for transient molten salt receiver system simulation JF - Solar Energy N2 - In order to realistically predict and optimize the actual performance of a concentrating solar power (CSP) plant sophisticated simulation models and methods are required. This paper presents a detailed dynamic simulation model for a Molten Salt Solar Tower (MST) system, which is capable of simulating transient operation including detailed startup and shutdown procedures including drainage and refill. For appropriate representation of the transient behavior of the receiver as well as replication of local bulk and surface temperatures a discretized receiver model based on a novel homogeneous two-phase (2P) flow modelling approach is implemented in Modelica Dymola®. This allows for reasonable representation of the very different hydraulic and thermal properties of molten salt versus air as well as the transition between both. This dynamic 2P receiver model is embedded in a comprehensive one-dimensional model of a commercial scale MST system and coupled with a transient receiver flux density distribution from raytracing based heliostat field simulation. This enables for detailed process prediction with reasonable computational effort, while providing data such as local salt film and wall temperatures, realistic control behavior as well as net performance of the overall system. Besides a model description, this paper presents some results of a validation as well as the simulation of a complete startup procedure. Finally, a study on numerical simulation performance and grid dependencies is presented and discussed. KW - Molten salt solar tower KW - Molten salt receiver system KW - Dynamic simulation KW - Two-phase modelling KW - Transient flux distribution Y1 - 2022 U6 - https://doi.org/10.1016/j.solener.2021.12.065 SN - 0038-092X (print) SN - 1471-1257 (online) VL - 232 SP - 362 EP - 375 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Bhattarai, Aroj A1 - May, Charlotte Anabell A1 - Staat, Manfred A1 - Kowalczyk, Wojciech A1 - Tran, Thanh Ngoc T1 - Layer-specific damage modeling of porcine large intestine under biaxial tension JF - Bioengineering N2 - The mechanical behavior of the large intestine beyond the ultimate stress has never been investigated. Stretching beyond the ultimate stress may drastically impair the tissue microstructure, which consequently weakens its healthy state functions of absorption, temporary storage, and transportation for defecation. Due to closely similar microstructure and function with humans, biaxial tensile experiments on the porcine large intestine have been performed in this study. In this paper, we report hyperelastic characterization of the large intestine based on experiments in 102 specimens. We also report the theoretical analysis of the experimental results, including an exponential damage evolution function. The fracture energies and the threshold stresses are set as damage material parameters for the longitudinal muscular, the circumferential muscular and the submucosal collagenous layers. A biaxial tensile simulation of a linear brick element has been performed to validate the applicability of the estimated material parameters. The model successfully simulates the biomechanical response of the large intestine under physiological and non-physiological loads. KW - biaxial tensile experiment KW - anisotropy KW - hyperelastic KW - constitutive modeling KW - damage Y1 - 2022 U6 - https://doi.org/10.3390/bioengineering9100528 SN - 2306-5354 N1 - Der Artikel gehört zum Sonderheft "Computational Biomechanics" VL - 9 IS - 10, Early Access SP - 1 EP - 17 PB - MDPI CY - Basel ER - TY - JOUR A1 - Mandekar, Swati A1 - Holland, Abigail A1 - Thielen, Moritz A1 - Behbahani, Mehdi A1 - Melnykowycz, Mark T1 - Advancing towards Ubiquitous EEG, Correlation of In-Ear EEG with Forehead EEG JF - Sensors N2 - Wearable EEG has gained popularity in recent years driven by promising uses outside of clinics and research. The ubiquitous application of continuous EEG requires unobtrusive form-factors that are easily acceptable by the end-users. In this progression, wearable EEG systems have been moving from full scalp to forehead and recently to the ear. The aim of this study is to demonstrate that emerging ear-EEG provides similar impedance and signal properties as established forehead EEG. EEG data using eyes-open and closed alpha paradigm were acquired from ten healthy subjects using generic earpieces fitted with three custom-made electrodes and a forehead electrode (at Fpx) after impedance analysis. Inter-subject variability in in-ear electrode impedance ranged from 20 kΩ to 25 kΩ at 10 Hz. Signal quality was comparable with an SNR of 6 for in-ear and 8 for forehead electrodes. Alpha attenuation was significant during the eyes-open condition in all in-ear electrodes, and it followed the structure of power spectral density plots of forehead electrodes, with the Pearson correlation coefficient of 0.92 between in-ear locations ELE (Left Ear Superior) and ERE (Right Ear Superior) and forehead locations, Fp1 and Fp2, respectively. The results indicate that in-ear EEG is an unobtrusive alternative in terms of impedance, signal properties and information content to established forehead EEG. KW - in-ear EEG KW - correlation KW - forehead EEG KW - impedance spectroscopy KW - biopotential electrodes Y1 - 2022 U6 - https://doi.org/10.3390/s22041568 SN - 1424-8220 VL - 22 IS - 4 SP - 1 EP - 19 PB - MDPI CY - Basel ER - TY - JOUR A1 - Maurischat, Andreas T1 - Algebraic independence of the Carlitz period and its hyperderivatives JF - Journal of Number Theory KW - Drinfeld modules KW - Periods KW - t-modules KW - Transcendence KW - Higher derivations Y1 - 2022 U6 - https://doi.org/10.1016/j.jnt.2022.01.006 SN - 0022-314X VL - 240 SP - 145 EP - 162 PB - Elsevier CY - Orlando, Fla. ER - TY - JOUR A1 - Schüller-Ruhl, Aaron A1 - Dinstühler, Leonard A1 - Senger, Thorsten A1 - Bergfeld, Stefan A1 - Ingenhag, Christian A1 - Fleischhaker, Robert ED - Mackenzie, Jacob T1 - Direct fabrication of arbitrary phase masks in optical glass via ultra-short pulsed laser writing of refractive index modifications JF - Applied Physics B N2 - We study the possibility to fabricate an arbitrary phase mask in a one-step laser-writing process inside the volume of an optical glass substrate. We derive the phase mask from a Gerchberg–Saxton-type algorithm as an array and create each individual phase shift using a refractive index modification of variable axial length. We realize the variable axial length by superimposing refractive index modifications induced by an ultra-short pulsed laser at different focusing depth. Each single modification is created by applying 1000 pulses with 15 μJ pulse energy at 100 kHz to a fixed spot of 25 μm diameter and the focus is then shifted axially in steps of 10 μm. With several proof-of-principle examples, we show the feasibility of our method. In particular, we identify the induced refractive index change to about a value of Δn=1.5⋅10−3. We also determine our current limitations by calculating the overlap in the form of a scalar product and we discuss possible future improvements. Y1 - 2022 U6 - https://doi.org/10.1007/s00340-022-07928-2 SN - 1432-0649 (Online) SN - 0946-2171 (Print) N1 - Corresponding author: Robert Fleischhaker VL - 128 IS - Article number: 208 SP - 1 EP - 11 PB - Springer CY - Berlin ER - TY - JOUR A1 - Philipp, Mohr A1 - Efthimiou, Nikos A1 - Pagano, Fiammetta A1 - Kratochwil, Nicolaus A1 - Pizzichemi, Marco A1 - Tsoumpas, Charalampos A1 - Auffray, Etiennette A1 - Ziemons, Karl T1 - Image reconstruction analysis for positron emission tomography with heterostructured scintillators JF - IEEE Transactions on Radiation and Plasma Medical Sciences N2 - The concept of structure engineering has been proposed for exploring the next generation of radiation detectors with improved performance. A TOF-PET geometry with heterostructured scintillators with a pixel size of 3.0×3.1×15 mm3 was simulated using Monte Carlo. The heterostructures consisted of alternating layers of BGO as a dense material with high stopping power and plastic (EJ232) as a fast light emitter. The detector time resolution was calculated as a function of the deposited and shared energy in both materials on an event-by-event basis. While sensitivity was reduced to 32% for 100 μm thick plastic layers and 52% for 50 μm, the CTR distribution improved to 204±49 ps and 220±41 ps respectively, compared to 276 ps that we considered for bulk BGO. The complex distribution of timing resolutions was accounted for in the reconstruction. We divided the events into three groups based on their CTR and modeled them with different Gaussian TOF kernels. On a NEMA IQ phantom, the heterostructures had better contrast recovery in early iterations. On the other hand, BGO achieved a better contrast to noise ratio (CNR) after the 15th iteration due to the higher sensitivity. The developed simulation and reconstruction methods constitute new tools for evaluating different detector designs with complex time responses. KW - TOF PET KW - Heterostructure KW - Metascintillator KW - Multiple TOF kernels KW - Image Reconstruction Y1 - 2022 U6 - https://doi.org/10.1109/TRPMS.2022.3208615 SN - 2469-7311 SN - 2469-7303 VL - 7 IS - 1 SP - 41 EP - 51 PB - IEEE CY - New York, NY ER - TY - JOUR A1 - Baringhaus, Ludwig A1 - Gaigall, Daniel T1 - A goodness-of-fit test for the compound Poisson exponential model JF - Journal of Multivariate Analysis N2 - On the basis of bivariate data, assumed to be observations of independent copies of a random vector (S,N), we consider testing the hypothesis that the distribution of (S,N) belongs to the parametric class of distributions that arise with the compound Poisson exponential model. Typically, this model is used in stochastic hydrology, with N as the number of raindays, and S as total rainfall amount during a certain time period, or in actuarial science, with N as the number of losses, and S as total loss expenditure during a certain time period. The compound Poisson exponential model is characterized in the way that a specific transform associated with the distribution of (S,N) satisfies a certain differential equation. Mimicking the function part of this equation by substituting the empirical counterparts of the transform we obtain an expression the weighted integral of the square of which is used as test statistic. We deal with two variants of the latter, one of which being invariant under scale transformations of the S-part by fixed positive constants. Critical values are obtained by using a parametric bootstrap procedure. The asymptotic behavior of the tests is discussed. A simulation study demonstrates the performance of the tests in the finite sample case. The procedure is applied to rainfall data and to an actuarial dataset. A multivariate extension is also discussed. KW - Bootstrapping KW - Collective risk model Y1 - 2022 U6 - https://doi.org/10.1016/j.jmva.2022.105154 SN - 0047-259X SN - 1095-7243 VL - 195 IS - Article 105154 PB - Elsevier CY - Amsterdam ER -