Search

This work utilizes the Immersed Boundary Conformal Method (IBCM) to analyze Kirchhoff-Love and Reissner-Mindlin shell structures within an immersed domain framework. Immersed boundary methods involve embedding complex geometries within a background grid, which allows for great flexibility in modeling intricate shapes and features despite the simplicity of the approach. The IBCM method introduces additional layers conformal to the boundaries, allowing for the strong imposition of Dirichlet boundary conditions and facilitating local refinement. In this study, the construction of boundary layers is combined with high-degree spline-based approximation spaces to further increase efficiency. The Nitsche method, employing non-symmetric average operators, is used to couple the boundary layers with the inner patch, while stabilizing the formulation with minimal penalty parameters. High-order quadrature rules are applied for integration over cut elements and patch interfaces. Numerical experiments demonstrate the efficiency and accuracy of the proposed formulation, highlighting its potential for complex shell structures modeled through Kirchhoff-Love and Reissner-Mindlin theories. These tests include the generation of conformal interfaces, the coupling of Kirchhoff-Love and Reissner-Mindlin theories, and the simulation of a damaged shell.

This work focuses on the coupling of trimmed shell patches using Isogeometric Analysis, based on higher continuity splines that seamlessly meet the $C^1$ requirement of Kirchhoff-Love-based discretizations. Weak enforcement of coupling conditions is achieved through the symmetric interior penalty method, where the fluxes are computed using their correct variationally consistent expression that was only recently proposed and is unprecedentedly adopted herein in the context of coupling conditions. The constitutive relationships account for generically laminated materials, although the proposed tests are conducted under the assumption of uniform thickness and lamination sequence. Numerical experiments assess the method for an isotropic and a laminated plate, as well as an isotropic hyperbolic paraboloid shell from the new shell obstacle course. The boundary conditions and domain force are chosen to reproduce manufactured analytical solutions, which are taken as reference to compute rigorous convergence curves in the $L^2$, $H^1$, and $H^2$ norms, that closely approach optimal ones predicted by theory. Additionally, we conduct a final test on a complex structure comprising five intersecting laminated cylindrical shells, whose geometry is directly imported from a STEP file. The results exhibit excellent agreement with those obtained through commercial software, showcasing the method's potential for real-world industrial applications., Comment: Engineering with Computers. 2024

Mass lumping techniques are commonly employed in explicit time integration schemes for problems in structural dynamics and both avoid solving costly linear systems with the consistent mass matrix and increase the critical time step. In isogeometric analysis, the critical time step is constrained by so-called "outlier" frequencies, representing the inaccurate high frequency part of the spectrum. Removing or dampening these high frequencies is paramount for fast explicit solution techniques. In this work, we propose robust mass lumping and outlier removal techniques for nontrivial geometries, including multipatch and trimmed geometries. Our lumping strategies provably do not deteriorate (and often improve) the CFL condition of the original problem and are combined with deflation techniques to remove persistent outlier frequencies. Numerical experiments reveal the advantages of the method, especially for simulations covering large time spans where they may halve the number of iterations with little or no effect on the numerical solution., Comment: 31 pages, 16 figures. Submitted manuscript

One of the most important objectives of solar physics is the physical understanding of the solar atmosphere, the structure of which is also described in terms of the density (N) and temperature (T) distributions of the atmospheric matter. Several multi-frequency analyses show that the characteristics of these distributions are still debated, especially for the outer coronal emission. We aim to constrain the T and N distributions of the solar atmosphere through observations in the centimetric radio domain. We employ single-dish observations from two of the INAF radio telescopes at the K-band frequencies (18 - 26 GHz). We investigate the origin of the significant brightness temperature ($T_B$) level that we detected up to the upper corona ($\sim 800$ Mm of altitude with respect to the photospheric solar surface). To probe the physical origin of the atmospheric emission and to constrain instrumental biases, we reproduced the solar signal by convolving specific 2D antenna beam models. The analysis of the solar atmosphere is performed by adopting a physical model that assumes the thermal bremsstrahlung as the emission mechanism, with specific T and N distributions. The modelled $T_B$ profiles are compared with those observed by averaging solar maps obtained during the minimum of solar activity (2018 - 2020). The T and N distributions are compatible (within $25\%$ of uncertainty) with the model up to $\sim 60$ Mm and $\sim 100$ Mm of altitude, respectively. The analysis of the role of the antenna beam pattern on our solar maps proves the physical nature of the atmospheric emission in our images up to the coronal tails seen in our $T_B$ profiles. The challenging analysis of the coronal radio emission at higher altitudes, together with the data from satellite instruments will require further multi-frequency measurements., Comment: 12 pages, 8 figures, 5 tables, accepted by A&A; v1

The Sun is an extraordinary workbench, from which several fundamental astronomical parameters can be measured with high precision. Among these parameters, the solar radius $R_{\odot}$ plays an important role in several aspects, such as in evolutionary models. Despite the efforts in obtaining accurate measurements of $R_{\odot}$, the subject is still debated and measurements are puzzling and/or lacking in many frequency ranges. We aimed to determine the mean, equatorial, and polar radii of the Sun ($R_c$, $R_{eq}$, and $R_{pol}$) in the frequency range 18.1 - 26.1 GHz. We employed single-dish observations from the newly-appointed Medicina "Gavril Grueff" Radio Telescope and the Sardinia Radio Telescope (SRT) throughout 5 years, from 2018 to mid-2023, in the framework of the SunDish project for solar monitoring. Two methods to calculate the radius at radio frequencies are considered and compared. To assess the quality of our radius determinations, we also analysed the possible degrading effects of the antenna beam pattern on our solar maps, using two 2D-models. We carried out a correlation analysis with the evolution of the solar cycle through the calculation of Pearson's correlation coefficient $\rho$. We obtained several values for the solar radius - ranging between 959 and 994 arcsec - and $\rho$, with typical errors of a few arcsec. Our $R_{\odot}$ measurements, consistent with values reported in literature, suggest a weak prolatness of the solar limb ($R_{eq}$ > $R_{pol}$), although $R_{eq}$ and $R_{pol}$ are statistically compatible within 3$\sigma$ errors. The correlation analysis using the solar images from Grueff shows (1) a positive correlation between the solar activity and the temporal variation of $R_c$ (and $R_{eq}$) at all observing frequencies, and (2) a weak anti-correlation between the temporal variation of $R_{pol}$ and the solar activity at 25.8 GHz., Comment: 18 pages, 12 figures, 6 tables, accepted by A&A; v1

An a posteriori error estimator based on an equilibrated flux reconstruction is proposed for defeaturing problems in the context of finite element discretizations. Defeaturing consists in the simplification of a geometry by removing features that are considered not relevant for the approximation of the solution of a given PDE. In this work, the focus is on Poisson equation with Neumann boundary conditions on the feature boundary. The estimator accounts both for the so-called defeaturing error and for the numerical error committed by approximating the solution on the defeatured domain. Unlike other estimators that were previously proposed for defeaturing problems, the use of the equilibrated flux reconstruction allows to obtain a sharp bound for the numerical component of the error. Furthermore, it does not require the evaluation of the normal trace of the numerical flux on the feature boundary: this makes the estimator well-suited for finite element discretizations, in which the normal trace of the numerical flux is typically discontinuous across elements. The reliability of the estimator is proven and verified on several numerical examples. Its capability to identify the most relevant features is also shown, in anticipation of a future application to an adaptive strategy.

Summary: GeneFEAST, implemented in Python, is a gene-centric functional enrichment analysis summarisation and visualisation tool that can be applied to large functional enrichment analysis (FEA) results arising from upstream FEA pipelines. It produces a systematic, navigable HTML report, making it easy to identify sets of genes putatively driving multiple enrichments and to explore gene-level quantitative data first used to identify input genes. Further, GeneFEAST can compare FEA results from multiple studies, making it possible, for example, to highlight patterns of gene expression amongst genes commonly differentially expressed in two sets of conditions, and giving rise to shared enrichments under those conditions. GeneFEAST offers a novel, effective way to address the complexities of linking up many overlapping FEA results to their underlying genes and data, advancing gene-centric hypotheses, and providing pivotal information for downstream validation experiments. Availability: GeneFEAST is available at https://github.com/avigailtaylor/GeneFEAST Contact: avigail.taylor@well.ox.ac.uk, Comment: Main text: 3 pages, 1 figure. Supplementary Information: 16 pages, 3 figures, 2 tables, 4 boxes

This paper addresses the overwhelming computational resources needed with standard numerical approaches to simulate architected materials. Those multiscale heterogeneous lattice structures gain intensive interest in conjunction with the improvement of additive manufacturing as they offer, among many others, excellent stiffness-to-weight ratios. We develop here a dedicated HPC solver that benefits from the specific nature of the underlying problem in order to drastically reduce the computational costs (memory and time) for the full fine-scale analysis of lattice structures. Our purpose is to take advantage of the natural domain decomposition into cells and, even more importantly, of the geometrical and mechanical similarities among cells. Our solver consists in a so-called inexact FETI-DP method where the local, cell-wise operators and solutions are approximated with reduced order modeling techniques. Instead of considering independently every cell, we end up with only few principal local problems to solve and make use of the corresponding principal cell-wise operators to approximate all the others. It results in a scalable algorithm that saves numerous local factorizations. Our solver is applied for the isogeometric analysis of lattices built by spline composition, which offers the opportunity to compute the reduced basis with macro-scale data, thereby making our method also multiscale and matrix-free. The solver is tested against various 2D and 3D analyses. It shows major gains with respect to black-box solvers; in particular, problems of several millions of degrees of freedom can be solved with a simple computer within few minutes., Comment: 30 pages, 12 figures, 2 tables

This contribution presents a model order reduction framework for real-time efficient solution of trimmed, multi-patch isogeometric Kirchhoff-Love shells. In several scenarios, such as design and shape optimization, multiple simulations need to be performed for a given set of physical or geometrical parameters. This step can be computationally expensive in particular for real world, practical applications. We are interested in geometrical parameters and take advantage of the flexibility of splines in representing complex geometries. In this case, the operators are geometry-dependent and generally depend on the parameters in a non-affine way. Moreover, the solutions obtained from trimmed domains may vary highly with respect to different values of the parameters. Therefore, we employ a local reduced basis method based on clustering techniques and the Discrete Empirical Interpolation Method to construct affine approximations and efficient reduced order models. In addition, we discuss the application of the reduction strategy to parametric shape optimization. Finally, we demonstrate the performance of the proposed framework to parameterized Kirchhoff-Love shells through benchmark tests on trimmed, multi-patch meshes including a complex geometry. The proposed approach is accurate and achieves a significant reduction of the online computational cost in comparison to the standard reduced basis method., Comment: 44 pages, 23 figures, 3 tables

This paper presents a PDE-based parameterisation framework for addressing the planar surface-to-volume (StV) problem of finding a valid description of the domain's interior given no more than a spline-based description of its boundary contours. The framework is geared towards isogeometric analysis (IGA) applications wherein the physical domain is comprised of more than four sides, hence requiring more than one patch. We adopt the concept of harmonic maps and propose several PDE-based problem formulations capable of finding a valid map between a convex parametric multipatch domain and the piecewise-smooth physical domain with an equal number of sides. In line with the isoparametric paradigm of IGA, we treat the StV problem using techniques that are characteristic for the analysis step. As such, this study proposes several IGA-based numerical algorithms for the problem's governing equations that can be effortlessly integrated into a well-developed IGA software suite. We augment the framework with mechanisms that enable controlling the parametric properties of the outcome. Parametric control is accomplished by, among other techniques, the introduction of a curvilinear coordinate system in the convex parametric domain that, depending on the application, builds desired features into the computed harmonic map, such as homogeneous cell sizes or boundary layers.

The objective of this study is to address the difficulty of simplifying the geometric model in which a differential problem is formulated, also called defeaturing, while simultaneously ensuring that the accuracy of the solution is maintained under control. This enables faster and more efficient simulations, without sacrificing accuracy. More precisely, we consider an isogeometric discretisation of an elliptic model problem defined on a two-dimensional hierarchical B-spline computational domain with a complex boundary. Starting with an oversimplification of the geometry, we build a goal-oriented adaptive strategy that adaptively reintroduces continuous geometrical features in regions where the analysis suggests a large impact on the quantity of interest. This strategy is driven by an a posteriori estimator of the defeaturing error based on first-order shape sensitivity analysis, and it profits from the local refinement properties of hierarchical B-splines. The adaptive algorithm is described together with a procedure to generate (partially) simplified hierarchical B-spline geometrical domains. Numerical experiments are presented to illustrate the proposed strategy and its limitations.

We present a novel stabilized isogeometric formulation for the Stokes problem, where the geometry of interest is obtained via overlapping NURBS (non-uniform rational B-spline) patches, i.e., one patch on top of another in an arbitrary but predefined hierarchical order. All the visible regions constitute the computational domain, whereas independent patches are coupled through visible interfaces using Nitsche's formulation. Such a geometric representation inevitably involves trimming, which may yield trimmed elements of extremely small measures (referred to as bad elements) and thus lead to the instability issue. Motivated by the minimal stabilization method that rigorously guarantees stability for trimmed geometries [1], in this work we generalize it to the Stokes problem on overlapping patches. Central to our method is the distinct treatments for the pressure and velocity spaces: Stabilization for velocity is carried out for the flux terms on interfaces, whereas pressure is stabilized in all the bad elements. We provide a priori error estimates with a comprehensive theoretical study. Through a suite of numerical tests, we first show that optimal convergence rates are achieved, which consistently agrees with our theoretical findings. Second, we show that the accuracy of pressure is significantly improved by several orders using the proposed stabilization method, compared to the results without stabilization. Finally, we also demonstrate the flexibility and efficiency of the proposed method in capturing local features in the solution field., Comment: 38 pages, 29 figures. Reviewed version

This work introduces a reduced order modeling (ROM) framework for the solution of parameterized second-order linear elliptic partial differential equations formulated on unfitted geometries. The goal is to construct efficient projection-based ROMs, which rely on techniques such as the reduced basis method and discrete empirical interpolation. The presence of geometrical parameters in unfitted domain discretizations entails challenges for the application of standard ROMs. Therefore, in this work we propose a methodology based on i) extension of snapshots on the background mesh and ii) localization strategies to decrease the number of reduced basis functions. The method we obtain is computationally efficient and accurate, while it is agnostic with respect to the underlying discretization choice. We test the applicability of the proposed framework with numerical experiments on two model problems, namely the Poisson and linear elasticity problems. In particular, we study several benchmarks formulated on two-dimensional, trimmed domains discretized with splines and we observe a significant reduction of the online computational cost compared to standard ROMs for the same level of accuracy. Moreover, we show the applicability of our methodology to a three-dimensional geometry of a linear elastic problem., Comment: 35 pages, 19 figures, 5 tables. Preprint accepted in Computer Methods in Applied Mechanics and Engineering

Removing geometrical details from a complex domain is a classical operation in computer aided design. This procedure simplifies the meshing process, and it enables faster simulations with less memory requirements. However, depending on the partial differential equation that one wants to solve, removing some important geometrical features may greatly impact the solution accuracy. Unfortunately, the effect of geometrical simplification on the accuracy of the problem solution is often neglected or its evaluation is based on engineering expertise, only due to the lack of reliable tools. It is therefore important to have a better understanding of the effect of geometrical model simplification, also called defeaturing, to improve our control on the simulation accuracy along the design and analysis phases. In this work, we consider as a model problem the Poisson equation on a geometry with Neumann features, we consider some finite element discretization of it, and we build an adaptive strategy that is twofold. Firstly, it is able to perform geometrical refinements, that is, to choose at each iteration step which geometrical feature is important to obtain an accurate solution. Secondly, it performs standard mesh refinements; since the geometry changes at each iteration, the algorithm is designed to be used with an immersed method. To drive this adaptive strategy, we introduce an a posteriori estimator of the energy error between the exact solution defined in the exact fully-featured geometry, and the numerical approximation of the solution defined in the defeatured geometry. The reliability of the estimator is proven for very general (potentially trimmed multipatch) geometric configurations, and in particular for IGA with hierarchical B-splines. Finally, numerical experiments are performed to validate the presented theory and to illustrate the capabilities of the proposed adaptive strategy., Comment: 50 pages

Explicit time integration schemes coupled with Galerkin discretizations of time-dependent partial differential equations require solving a linear system with the mass matrix at each time step. For applications in structural dynamics, the solution of the linear system is frequently approximated through so-called mass lumping, which consists in replacing the mass matrix by some diagonal approximation. Mass lumping has been widely used in engineering practice for decades already and has a sound mathematical theory supporting it for finite element methods using the classical Lagrange basis. However, the theory for more general basis functions is still missing. Our paper partly addresses this shortcoming. Some special and practically relevant properties of lumped mass matrices are proved and we discuss how these properties naturally extend to banded and Kronecker product matrices whose structure allows to solve linear systems very efficiently. Our theoretical results are applied to isogeometric discretizations but are not restricted to them., Comment: 29 pages, 26 figures. Accepted manuscript

This contribution explores the combined capabilities of reduced basis methods and IsoGeometric Analysis (IGA) in the context of parameterized partial differential equations. The introduction of IGA enables a unified simulation framework based on a single geometry representation for both design and analysis. The coupling of reduced basis methods with IGA has been motivated in particular by their combined capabilities for geometric design and solution of parameterized geometries. In most IGA applications, the geometry is modelled by multiple patches with different physical or geometrical parameters. In particular, we are interested in nonaffine problems characterized by a high-dimensional parameter space. We consider the Empirical Interpolation Method (EIM) to recover an affine parametric dependence and combine domain decomposition to reduce the dimensionality. We couple spline patches in a parameterized setting, where multiple evaluations are performed for a given set of geometrical parameters, and employ the Static Condensation Reduced Basis Element (SCRBE) method. At the common interface between adjacent patches a static condensation procedure is employed, whereas in the interior a reduced basis approximation enables an efficient offline/online decomposition. The full order model over which we setup the RB formulation is based on NURBS approximation, whereas the reduced basis construction relies on techniques such as the Greedy algorithm or proper orthogonal decomposition (POD). We demonstrate the developed procedure using an illustrative model problem on a three-dimensional geometry featuring a multi-dimensional geometrical parameterization., Comment: 18 pages, 7 figures, 1 table

We present a new solar radio imaging system implemented through the upgrade of the large single-dish telescopes of the Italian National Institute for Astrophysics (INAF), not originally conceived for solar observations. During the development and early science phase of the project (2018-2020), we obtained about 170 maps of the entire solar disk in the 18-26 GHz band, filling the observational gap in the field of solar imaging at these frequencies. These solar images have typical resolutions in the 0.7-2 arcmin range and a brightness temperature sensitivity <10 K. Accurate calibration adopting the Supernova Remnant Cas A as a flux reference, provided typical errors <3% for the estimation of the quiet-Sun level components and for active regions flux measurements. As a first early science result of the project, we present a catalog of radio continuum solar imaging observations with Medicina 32-m and SRT 64-m radio telescopes including the multi-wavelength identification of active regions, their brightness and spectral characterization. The interpretation of the observed emission as thermal bremsstrahlung components combined with gyro-magnetic variable emission pave the way to the use of our system for long-term monitoring of the Sun. We also discuss useful outcomes both for solar physics (e.g. study of the chromospheric network dynamics) and space weather applications (e.g. flare precursors studies)., Comment: 43 pages, 11 figures, 6 tables, accepted for publication in Solar Physics