A finite-element toolbox for the stationary Gross-Pitaevskii equation with rotation

We present a new numerical system using classical finite elements with mesh adaptivity for computing stationary solutions of the Gross-Pitaevskii equation. The programs are written as a toolbox for FreeFem++ (www.freefem.org), a free finite-element software available for all existing operating systems. This offers the advantage to hide all technical issues related to the implementation of the finite element method, allowing to easily implement various numerical algorithms.Two robust and optimised numerical methods were implemented to minimize the Gross-Pitaevskii energy: a steepest descent method based on Sobolev gradients and a minimization algorithm based on the state-of-the-art optimization library Ipopt. For both methods, mesh adaptivity strategies are implemented to reduce the computational time and increase the local spatial accuracy when vortices are present. Different run cases are made available for 2D and 3D configurations of Bose-Einstein condensates in rotation. An optional graphical user interface is also provided, allowing to easily run predefined cases or with user-defined parameter files. We also provide several post-processing tools (like the identification of quantized vortices) that could help in extracting physical features from the simulations. The toolbox is extremely versatile and can be easily adapted to deal with different physical models

Interception and deviation of near Earth objects via solar collector strategy

A solution to the asteroid deviation problem via a low-thrust strategy is proposed. This formulation makes use of the proximal motion equations and a semi-analytical solution of the Gauss planetary equations. The average of the variation of the orbital elements is computed, together with an approximate expression of their periodic evolution. The interception and the deflection phase are optimised together through a global search. The low-thrust transfer is preliminary designed with a shape based method; subsequently the solutions are locally refined through the Differential Dynamic Programming approach. A set of optimal solutions are presented for a deflection mission to Apophis, together with a representative trajectory to Apophis including the Earth escape

Variational assimilation for xenon dynamical forecasts in neutronic using advanced background error covariance matrix modelling

Data assimilation method consists in combining all available pieces of information about a system to obtain optimal estimates of initial states. The different sources of information are weighted according to their accuracy by the means of error covariance matrices. Our purpose here is to evaluate the efficiency of variational data assimilation for the xenon induced oscillations forecasts in nuclear cores. In this paper we focus on the comparison between 3DVAR schemes with optimised background error covariance matrix B and a 4DVAR scheme. Tests were made in twin experiments using a simulation code which implements a mono-dimensional coupled model of xenon dynamics, thermal, and thermal–hydraulic processes. We enlighten the very good efficiency of the 4DVAR scheme as well as good results with the 3DVAR one using a careful multivariate modelling of B

Optimising Spatial and Tonal Data for PDE-based Inpainting

Some recent methods for lossy signal and image compression store only a few selected pixels and fill in the missing structures by inpainting with a partial differential equation (PDE). Suitable operators include the Laplacian, the biharmonic operator, and edge-enhancing anisotropic diffusion (EED). The quality of such approaches depends substantially on the selection of the data that is kept. Optimising this data in the domain and codomain gives rise to challenging mathematical problems that shall be addressed in our work. In the 1D case, we prove results that provide insights into the difficulty of this problem, and we give evidence that a splitting into spatial and tonal (i.e. function value) optimisation does hardly deteriorate the results. In the 2D setting, we present generic algorithms that achieve a high reconstruction quality even if the specified data is very sparse. To optimise the spatial data, we use a probabilistic sparsification, followed by a nonlocal pixel exchange that avoids getting trapped in bad local optima. After this spatial optimisation we perform a tonal optimisation that modifies the function values in order to reduce the global reconstruction error. For homogeneous diffusion inpainting, this comes down to a least squares problem for which we prove that it has a unique solution. We demonstrate that it can be found efficiently with a gradient descent approach that is accelerated with fast explicit diffusion (FED) cycles. Our framework allows to specify the desired density of the inpainting mask a priori. Moreover, is more generic than other data optimisation approaches for the sparse inpainting problem, since it can also be extended to nonlinear inpainting operators such as EED. This is exploited to achieve reconstructions with state-of-the-art quality. We also give an extensive literature survey on PDE-based image compression methods

Stochastic level-set method for shape optimisation

We present a new method for stochastic shape optimisation of engineering structures. The method generalises an existing deterministic scheme, in which the structure is represented and evolved by a level-set method coupled with mathematical programming. The stochastic element of the algorithm is built on the methods of statistical mechanics and is designed so that the system explores a Boltzmann-Gibbs distribution of structures. In non-convex optimisation problems, the deterministic algorithm can get trapped in local optima: the stochastic generalisation enables sampling of multiple local optima, which aids the search for the globally-optimal structure. The method is demonstrated for several simple geometrical problems, and a proof-of-principle calculation is shown for a simple engineering structure.Comment: 17 pages, 10 fig