Optimizing feedstock imports with environmental constraints

In addressing the problem of commodity production out of feedstock imports, an eco-environmentally rational agent aims at minimizing the cost of feedstock imports and their transportation, but also the water footprint of the feedstock production process and the water scarcity in the exporting countries. The problem is formulated as a nonlinear program. This study proves the existence of solutions and quantitatively demonstrates that transportation costs and non-uniform feedstock characteristics inhibit feedstock interchangeability. Moreover, it is shown that the interplay between water footprint and water scarcity across countries can inhibit or foster feedstock interchangeability.Comment: 14 page

Turing patterns in a 3D morpho-chemical bulk-surface reaction-diffusion system for battery modeling

In this paper we introduce a bulk-surface reaction-diffusion (BSRD) model in three space dimensions that extends the DIB morphochemical model to account for the electrolyte contribution in the application, in order to study structure formation during discharge-charge processes in batteries. Here we propose to approximate the model by the Bulk-Surface Virtual Element Method on a tailor-made mesh that proves to be competitive with fast bespoke methods for PDEs on Cartesian grids. We present a selection of numerical simulations that accurately match the classical morphologies found in experiments. Finally, we compare the Turing patterns obtained by the coupled 3D BS-DIB model with those obtained with the original 2D version.Comment: 25 pages, 11 figures, 1 tabl

Turing pattern formation on the sphere for a morphochemical reaction-diffusion model for electrodeposition

The present paper deals with the pattern formation properties of a specific morpho- electrochemical reaction-diffusion model on a sphere. The physico-chemical background to this study is the morphological control of material electrodeposited onto spherical parti- cles. The particular experimental case of interest refers to the optimization of novel metal- air flow batteries and addresses the electrodeposition of zinc onto inert spherical supports. Morphological control in this step of the high-energy battery operation is crucial to the energetic efficiency of the recharge process and to the durability of the whole energy- storage device. To rationalise this technological challenge within a mathematical modeling perspective, we consider the reaction-diffusion system for metal electrodeposition intro- duced in [Bozzini et al., J. Solid State Electr.17, 467–479 (2013)] and extend its study to spherical domains. Conditions are derived for the occurrence of the Turing instability phe- nomenon and the steady patterns emerging at the onset of Turing instability are investi- gated. The reaction-diffusion system on spherical domains is solved numerically by means of the Lumped Surface Finite Element Method (LSFEM) in space combined with the IMEX Euler method in time. The effect on pattern formation of variations in the domain size is investigated both qualitatively, by means of systematic numerical simulations, and quan- titatively by introducing suitable indicators that allow to assign each pattern to a given morphological class. An experimental validation of the obtained results is finally presented for the case of zinc electrodeposition from alkaline zincate solutions onto copper spheres

Lumped finite elements for reaction–cross-diffusion systems on stationary surfaces

We consider a lumped surface finite element method (LSFEM) for the spatial approximation of reaction–diffusion equations on closed compact surfaces in R3R3 in the presence of cross-diffusion. We provide a fully-discrete scheme by applying the Implicit–Explicit (IMEX) Euler method. We provide sufficient conditions for the existence of polytopal invariant regions for the numerical solution after spatial and full discretisations. Furthermore, we prove optimal error bounds for the semi- and fully-discrete methods, that is the convergence rates are quadratic in the meshsize and linear in the timestep. To support our theoretical findings, we provide two numerical tests. The first test confirms that in the absence of lumping numerical solutions violate the invariant region leading to blow-up due to the nature of the kinetics. The second experiment is an example of Turing pattern formation in the presence of cross-diffusion on the sphere

Numerical preservation of velocity induced invariant regions for reaction-diffusion systems on evolving surfaces

We propose and analyse a finite element method with mass lumping (LESFEM) for the numerical approximation of reaction-diffusion systems (RDSs) on surfaces in R3 that evolve under a given velocity field. A fully-discrete method based on the implicit-explicit (IMEX) Euler time-discretisation is formulated and dilation rates which act as indicators of the surface evolution are introduced. Under the assumption that the mesh preserves the Delaunay regularity under evolution, we prove a sufficient condition, that depends on the dilation rates, for the existence of invariant regions (i) at the spatially discrete level with no restriction on the mesh size and (ii) at the fully-discrete level under a timestep restriction that depends on the kinetics, only. In the specific case of the linear heat equation, we prove a semi- and a fully-discrete maximum principle. For the well-known activator-depleted and Thomas reaction-diffusion models we prove the existence of a family of rectangles in the phase space that are invariant only under specific growth laws. Two numerical examples are provided to computationally demonstrate (i) the discrete maximum principle and optimal convergence for the heat equation on a linearly growing sphere and (ii) the existence of an invariant region for the LESFEM-IMEX Euler discretisation of a RDS on a logistically growing surface

A VLBI experiment using a remote atomic clock via a coherent fibre link

We describe a VLBI experiment in which, for the first time, the clock reference is delivered from a National Metrology Institute to a radio telescope using a coherent fibre link 550 km long. The experiment consisted of a 24-hours long geodetic campaign, performed by a network of European telescopes; in one of those (Medicina, Italy) the local clock was alternated with a signal generated from an optical comb slaved to a fibre-disseminated optical signal. The quality of the results obtained with this facility and with the local clock is similar: interferometric fringes were detected throughout the whole 24-hours period and it was possible to obtain a solution whose residuals are comparable to those obtained with the local clock. These results encourage further investigation of the ultimate VLBI performances achievable using fibre dissemination at the highest precision of state-of-the-art atomic clocks

Measuring absolute frequencies beyond the GPS limit via long-haul optical frequency dissemination

Global Positioning System (GPS) dissemination of frequency standards is ubiquitous at present, providing the most widespread time and frequency reference for the majority of industrial and research applications worldwide. On the other hand, the ultimate limits of the GPS presently curb further advances in high-precision, scientific and industrial applications relying on this dissemination scheme. Here, we demonstrate that these limits can be reliably overcome even in laboratories without a local atomic clock by replacing the GPS with a 642-km-long optical fiber link to a remote primary caesium frequency standard. Through this configuration we stably address the 1S0—3P0 clock transition in an ultracold gas of 173Yb, with a precision that exceeds the possibilities of a GPS-based measurement, dismissing the need for a local clock infrastructure to perform beyond-GPS high-precision tasks. We also report an improvement of two orders of magnitude in the accuracy on the transition frequency reported in literature

The First Geodetic VLBI Field Test of LIFT: a 550-km-long Optical Fiber Link for Remote Antenna Synchronization

We present the first field test of the implementation of a coherent optical fiber link for remote antenna synchronization realized in Italy between the Italian Metrological Institute (INRIM) and the Medicina radio observatory of the National Institute for Astrophysics (INAF). The Medicina VLBI antenna participated in the EUR137 experiment carried out in September 2015 using, as reference systems, both the local H-maser and a remote H-maser hosted at the INRIM labs in Turin, separated by about 550 km. In order to assess the quality of the remote clock, the observed radio sources were split into two sets, using either the local or the remote H-maser. A system to switch automatically between the two references was integrated into the antenna field system. The observations were correlated in Bonn and preliminary results are encouraging since fringes were detected with both time references along the full 24 hours of the session. The experimental set-up, the results, and the perspectives for future radio astronomical and geodetic experiments are presented

Preserving invariance properties of reaction–diffusion systems on stationary surfaces

We propose and analyse a lumped surface finite element method for the numerical approximation of reaction–diffusion systems on stationary compact surfaces in R3. The proposed method preserves the invariant regions of the continuous problem under discretization and, in the special case of scalar equations, it preserves the maximum principle. On the application of a fully discrete scheme using the implicit–explicit Euler method in time, we prove that invariant regions of the continuous problem are preserved (i) at the spatially discrete level with no restriction on the meshsize and (ii) at the fully discrete level under a timestep restriction. We further prove optimal error bounds for the semidiscrete and fully discrete methods, that is, the convergence rates are quadratic in the meshsize and linear in the timestep. Numerical experiments are provided to support the theoretical findings. We provide examples in which, in the absence of lumping, the numerical solution violates the invariant region leading to blow-up