Source author record

Francesco Biscani

Francesco Biscani 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

11works
14topics
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

11 published item(s)

preprint2022arXiv

Reliable event detection for Taylor methods in astrodynamics

We present a novel approach for the detection of events in systems of ordinary differential equations. The new method combines the unique features of Taylor integrators with state-of-the-art polynomial root finding techniques to yield a novel algorithm ensuring strong event detection guarantees at a modest computational overhead. Detailed tests and benchmarks focused on problems in astrodynamics and celestial mechanics (such as collisional N-body systems, spacecraft dynamics around irregular bodies accounting for eclipses, computation of Poincare' sections, etc.) show how our approach is superior in both performance and detection accuracy to strategies commonly employed in modern numerical integration works. The new algorithm is available in our open source Taylor integration package heyoka.

preprint2016arXiv

Differentiable Genetic Programming

We introduce the use of high order automatic differentiation, implemented via the algebra of truncated Taylor polynomials, in genetic programming. Using the Cartesian Genetic Programming encoding we obtain a high-order Taylor representation of the program output that is then used to back-propagate errors during learning. The resulting machine learning framework is called differentiable Cartesian Genetic Programming (dCGP). In the context of symbolic regression, dCGP offers a new approach to the long unsolved problem of constant representation in GP expressions. On several problems of increasing complexity we find that dCGP is able to find the exact form of the symbolic expression as well as the constants values. We also demonstrate the use of dCGP to solve a large class of differential equations and to find prime integrals of dynamical systems, presenting, in both cases, results that confirm the efficacy of our approach.

preprint2016arXiv

On the astrodynamics applications of Weierstrass elliptic and related functions

Weierstrass elliptic and related functions have been recently shown to enable analytical explicit solutions to classical problems in astrodynamics. These include the constant radial acceleration problem, the Stark problem and the two-fixed center (or Euler's) problem. In this paper we review the basic technique that allows for these results and we discuss the limits and merits of the approach. Applications to interplanetary trajectory design are then discussed including low-thrust planetary fly-bys and the motion of an artificial satellite under the influence of an oblate primary including $J_2$ and $J_3$ harmonics.

preprint2015arXiv

A complete and explicit solution to the three-dimensional problem of two fixed centres

We present for the first time an explicit, complete and closed-form solution to the three-dimensional problem of two fixed centres, based on Weierstrass elliptic and related functions. With respect to previous treatments of the problem, our solution is exact, valid for all initial conditions and physical parameters of the system (including unbounded orbits and repulsive forces), and expressed via a unique set of formulae. Various properties of the three-dimensional problem of two fixed centres are investigated and analysed, with a particular emphasis on quasi-periodic and periodic orbits, regions of motion and equilibrium points.

preprint2015arXiv

Explicit Solution to the Constant Radial Acceleration Problem

While the constant radial acceleration problem is known to be integrable and has received some recent attention in an orbital mechanics context, a closed form explicit solution, relating the state variables to a time parameter, has eluded all researchers so far. It is here shown how such a solution exists and is elegantly expressed in terms of the Weierstrass elliptic and related functions. Previously known facts can be derived from the new explicit solution and new insights are revealed.

preprint2014arXiv

The Stark problem in the Weierstrassian formalism

We present a new general, complete closed-form solution of the Stark problem in terms of Weierstrass elliptic and related functions. With respect to previous treatments of the problem, our analysis is exact and valid for all values of the external force field, and it is expressed via unique formulæ valid for all initial conditions and parameters of the system. The simple form of the solution allows us to perform a thorough investigation of the properties of the dynamical system, including the identification of quasi-periodic and periodic orbits, the formulation of a simple analytical criterion to determine the boundness of the trajectory, and the characterisation of the equilibrium points.

preprint2013arXiv

A first-order secular theory for the post-Newtonian two-body problem with spin -- I: The restricted case

We revisit the relativistic restricted two-body problem with spin employing a perturbation scheme based on Lie series. Starting from a post-Newtonian expansion of the field equations, we develop a first-order secular theory that reproduces well-known relativistic effects such as the precession of the pericentre and the Lense-Thirring and geodetic effects. Additionally, our theory takes into full account the complex interplay between the various relativistic effects, and provides a new explicit solution of the averaged equations of motion in terms of elliptic functions. Our analysis reveals the presence of particular configurations for which non-periodical behaviour can arise. The application of our results to real astrodynamical systems (such as Mercury-like and pulsar planets) highlights the contribution of relativistic effects to the long-term evolution of the spin and orbit of the secondary body.

preprint2010arXiv

A Global Optimisation Toolbox for Massively Parallel Engineering Optimisation

A software platform for global optimisation, called PaGMO, has been developed within the Advanced Concepts Team (ACT) at the European Space Agency, and was recently released as an open-source project. PaGMO is built to tackle high-dimensional global optimisation problems, and it has been successfully used to find solutions to real-life engineering problems among which the preliminary design of interplanetary spacecraft trajectories - both chemical (including multiple flybys and deep-space maneuvers) and low-thrust (limited, at the moment, to single phase trajectories), the inverse design of nano-structured radiators and the design of non-reactive controllers for planetary rovers. Featuring an arsenal of global and local optimisation algorithms (including genetic algorithms, differential evolution, simulated annealing, particle swarm optimisation, compass search, improved harmony search, and various interfaces to libraries for local optimisation such as SNOPT, IPOPT, GSL and NLopt), PaGMO is at its core a C++ library which employs an object-oriented architecture providing a clean and easily-extensible optimisation framework. Adoption of multi-threaded programming ensures the efficient exploitation of modern multi-core architectures and allows for a straightforward implementation of the island model paradigm, in which multiple populations of candidate solutions asynchronously exchange information in order to speed-up and improve the optimisation process. In addition to the C++ interface, PaGMO's capabilities are exposed to the high-level language Python, so that it is possible to easily use PaGMO in an interactive session and take advantage of the numerous scientific Python libraries available.

preprint2010arXiv

Multiplication of sparse Laurent polynomials and Poisson series on modern hardware architectures

In this paper we present two algorithms for the multiplication of sparse Laurent polynomials and Poisson series (the latter being algebraic structures commonly arising in Celestial Mechanics from the application of perturbation theories). Both algorithms first employ the Kronecker substitution technique to reduce multivariate multiplication to univariate multiplication, and then use the schoolbook method to perform the univariate multiplication. The first algorithm, suitable for moderately-sparse multiplication, uses the exponents of the monomials resulting from the univariate multiplication as trivial hash values in a one dimensional lookup array of coefficients. The second algorithm, suitable for highly-sparse multiplication, uses a cache-optimised hash table which stores the coefficient-exponent pairs resulting from the multiplication using the exponents as keys. Both algorithms have been implemented with attention to modern computer hardware architectures. Particular care has been devoted to the efficient exploitation of contemporary memory hierarchies through cache-blocking techniques and cache-friendly term ordering. The first algorithm has been parallelised for shared-memory multicore architectures, whereas the second algorithm is in the process of being parallelised. We present benchmarks comparing our algorithms to the routines of other computer algebra systems, both in sequential and parallel mode.

preprint2010arXiv

On the Impact of the Migration Topology on the Island Model

Parallel Global Optimization Algorithms (PGOA) provide an efficient way of dealing with hard optimization problems. One method of parallelization of GOAs that is frequently applied and commonly found in the contemporary literature is the so-called Island Model (IM). In this paper we analyze the impact of the migration topology on the performance of a PGOA which uses the Island Model. In particular we consider parallel Differential Evolution and Simulated Annealing with Adaptive Neighborhood and draw first conclusions that emerge from the conducted experiments.

preprint2010arXiv

Towards a High Fidelity Direct Transcription Method for Optimisation of Low-Thrust Trajectories

We build upon some new ideas in direct transcription methods developed within the Advanced Concepts Team to introduce two improvements to the Sims-Flanagan transcription for low-thrust trajectories. The obtained new algorithm is able to produce an operational trajectory accounting for the real spacecraft dynamics and adapting the segment duration on-line improving the final trajectory optimality.