Source author record

Luca Bonaventura

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

10 published item(s)

preprint2023arXiv

Ensemble Reservoir Computing for Dynamical Systems: Prediction of Phase-Space Stable Region for Hadron Storage Rings

We investigate the ability of an ensemble reservoir computing approach to predict the long-term behaviour of the phase-space region in which the motion of charged particles in hadron storage rings is bounded, the so-called dynamic aperture. Currently, the calculation of the phase-space stability region of hadron storage rings is performed through direct computer simulations, which are resource- and time-intensive processes. Echo State Networks (ESN) are a class of recurrent neural networks that are computationally effective, since they avoid backpropagation and require only cross-validation. Furthermore, they have been proven to be universal approximants of dynamical systems. In this paper, we present the performance reached by ESN based on an ensemble approach for the prediction of the phase-space stability region and compare it with analytical scaling laws based on the stability-time estimate of the Nekhoroshev theorem for Hamiltonian systems. We observe that the proposed ESN approach is capable of effectively predicting the time evolution of the extent of the dynamic aperture, improving the predictions by analytical scaling laws, thus providing an efficient surrogate model.

preprint2020arXiv

A hyperbolic reformulation of the Serre-Green-Naghdi model for general bottom topographies

We present a novel hyperbolic reformulation of the Serre-Green-Naghdi (SGN) model for the description of dispersive water waves. Contrarily to the classical Boussinesq-type models, it contains only first order derivatives, thus allowing to overcome the numerical difficulties and the severe time step restrictions arising from higher order terms. The proposed model reduces to the original SGN model when an artificial sound speed tends to infinity. Moreover, it is endowed with an energy conservation law from which the energy conservation law associated with the original SGN model is retrieved when the artificial sound speed goes to infinity. The governing partial differential equations are then solved at the aid of high order ADER discontinuous Galerkin finite element schemes. The new model has been successfully validated against numerical and experimental results, for both flat and non-flat bottom. For bottom topographies with large variations, the new model proposed in this paper provides more accurate results with respect to the hyperbolic reformulation of the SGN model with the mild bottom approximation recently proposed in "C. Escalante, M. Dumbser and M.J. Castro. An efficient hyperbolic relaxation system for dispersive non-hydrostatic water waves and its solution with high order discontinuous Galerkin schemes, Journal of Computational Physics 2018".

preprint2020arXiv

A novel dowscaling procedure for compositional data in the Aitchison geometry with application to soil texture data

In this work, we present a novel downscaling procedure for compositional quantities based on the Aitchison geometry. The method is able to naturally consider compositional constraints, i.e. unit-sum and positivity. We show that the method can be used in a block sequential Gaussian simulation framework in order to assess the variability of downscaled quantities. Finally, to validate the method, we test it first in an idealized scenario and then apply it for the downscaling of digital soil maps on a more realistic case study. The digital soil maps for the realistic case study are obtained from SoilGrids, a system for automated soil mapping based on state-of-the-art spatial predictions methods.

preprint2020arXiv

Second order fully semi-Lagrangian discretizations of advection-diffusion-reaction systems

We propose a second order, fully semi-Lagrangian method for the numerical solution of systems of advection-diffusion-reaction equations, which employs a semi-Lagrangian approach to approximate in time both the advective and the diffusive terms. Standard interpolation procedures are used for the space discretization on structured and unstructured meshes. The proposed method allows for large time steps, while avoiding the solution of large linear systems, which would be required by an implicit time discretization technique. Numerical experiments demonstrate the effectiveness of the proposed approach and its superior efficiency with respect to more conventional explicit and implicit time discretizations

preprint2019arXiv

Dynamical p-adaptivity for LES of compressible flows in a high order DG framework

We investigate the possibility of reducing the computational burden of LES models by employing locally and dynamically adaptive polynomial degrees in the framework of a high order DG method. A degree adaptation technique especially featured to be effective for LES applications, that was previously developed by the authors and tested in the statically adaptive case, is applied here in a dynamically adaptive fashion. Two significant benchmarks are considered, comparing the results of adaptive and non adaptive simulations. The proposed dynamically adaptive approach allows for a significant reduction of the computational cost of representative LES computation, while allowing to maintain the level of accuracy guaranteed by LES carried out with constant, maximum polynomial degree values.

preprint2019arXiv

IsoGeometric Approximations for Electromagnetic Problems in Axisymmetric Domains

We propose a numerical method for the solution of electromagnetic problems on axisymmetric domains, based on a combination of a spectral Fourier approximation in the azimuthal direction with an IsoGeometric Analysis (IGA) approach in the radial and axial directions. This combination allows to blend the flexibility and accuracy of IGA approaches with the advantages of a Fourier representation on axisymmetric domains. It also allows to reduce significantly the computational cost by decoupling of the computations required for each Fourier mode. We prove that the discrete approximation spaces employed functional space constitute a closed and exact de Rham sequence. Numerical simulations of relevant benchmarks confirm the high order convergence and other computational advantages of the proposed method.

preprint2015arXiv

Flux form Semi-Lagrangian methods for parabolic problems

A semi-Lagrangian method for parabolic problems is proposed, that extends previous work by the authors to achieve a fully conservative, flux-form discretization of linear and nonlinear diffusion equations. A basic consistency and convergence analysis are proposed. Numerical examples validate the proposed method and display its potential for consistent semi-Lagrangian discretization of advection--diffusion and nonlinear parabolic problems.

preprint2015arXiv

Local Exponential Methods: a domain decomposition approach to exponential time integration of PDEs

A local approach to the time integration of PDEs by exponential methods is proposed, motivated by theoretical estimates by A.Iserles on the decay of off-diagonal terms in the exponentials of sparse matrices. An overlapping domain decomposition technique is outlined, that allows to replace the computation of a global exponential matrix by a number of independent and easily parallelizable local problems. Advantages and potential problems of the proposed technique are discussed. Numerical experiments on simple, yet relevant model problems show that the resulting method allows to increase computational efficiency with respect to standard implementations of exponential methods.

preprint2015arXiv

Monotonicity, positivity and strong stability of the TR-BDF2 method and of its SSP extensions

We analyze the one-step method TR-BDF2 from the point of view of monotonicity, strong stability and positivity. All these properties are strongly related and reviewed in the common framework of absolute monotonicity. The radius of absolute monotonicity is computed and it is shown that the parameter value which makes the method L-stable is also the value which maximizes the radius of monotonicity. Two hybrid variants of TR-BDF2 are proposed, that reduce the formal order of accuracy and maximize the absolute monotonicity radius, while keeping the native L-stability useful in stiff problems. Numerical experiments compare these different hybridization strategies with other methods commonly used in the presence of stiff and mildly stiff source terms. The results show that both strategies provide a good compromise between accuracy and robustness at high CFL numbers, without suffering from the limitations of alternative approaches already available in literature.

preprint2014arXiv

Anisotropic dynamic models for Large Eddy Simulation of compressible flows with a high order DG method

The impact of anisotropic dynamic models for applications to LES of compressible flows is assessed in the framework of a numerical model based on high order discontinuous finite elements. The projections onto lower dimensional subspaces associated to lower degree basis function are used as LES filter, along the lines proposed in Variational Multiscale templates. Comparisons with DNS results available in the literature for channel flows at Mach numbers 0.2, 0.7 and 1.5 show clearly that the anisotropic model is able to reproduce well some key features of the flow, especially close to the wall, where the flow anisotropy plays a major role.