Sono state esposte le sintesi delle conoscenze scientifiche acquisite sino ad oggi sulle Antiche Vie dell'Ambra, dal Mesolitico fino all'Antica Roma e dal Mar Baltico all'Europa Meridionale, al Mar Nero ed al Medio Oriente, con riferimenti ai ritrovamenti più recenti, ai giacimenti, con particolare riguardo all'ambra baltica, e sono state avanzate contemporaneamente ipotesi o riferimenti sui tragitti relativi. Una particolare evidenza è stata data alle vie che attraverso lo scambio commerciale hanno determinato la diffusione delle diverse culture fra i popoli e quindi maggior valore assumono in questo contesto le antichissime vie di epoche in cui i contatti e la comunicazione fra i popoli erano in dipendenza soprattutto delle vie del commercio. Grande importanza hanno assunto, alla luce delle superiori considerazioni, le vie che in epoca protostorica, ed in particolare nelle età del Bronzo e del Ferro, hanno consentito e caratterizzato i commerci dell'Ambra, oltreché di altri materiali, dal Mar Baltico, all'Adriatico, all'Egeo, al Mar Nero ed al Medio Oriente, mettendo in contatto popoli e culture completamente diversi fra loro, contribuendo senza dubbio allo sviluppo delle civiltà e delle conoscenze. Le vie di pace, come quelle dell'ambra, hanno sempre generato cultura.

We propose a numerical method to solve the Wigner equation in quantum systems of spinless, non-relativistic particles. The method uses a spectral decomposition into L-2(R-d) basis functions in momentum-space to obtain a system of first-order advection-reaction equations. The resulting equations are solved by splitting the reaction and advection steps so as to allow the combination of numerical techniques from quantum mechanics and computational fluid dynamics by identifying the skew-hermitian reaction matrix as a generator of unitary rotations. The method is validated for the case of particles subject to a one-dimensional (an-)harmonic and Morse potential using finite-differences for the advection part. Thereby, we verify the second order of convergence and observe non-classical behavior in the evolution of the Wigner function. (C) 2015 Elsevier Inc. All rights reserved.

Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin and A. Ajdari, Phys. Rev. Lett., 2009, 103, 036001, we show that cooperativity effects resulting from the nonlocal nature of the fluidity (inverse viscosity) are intimately related to the emergence of shear-banding configurations. This connection materializes through the onset of inhomogeneous compact solutions (compactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. The compacton coexistence with regions of zero fluidity ("non-flowing vacuum") is shown to be stabilized by the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material.

In this article, we study in detail the fluid dynamics system proposed in Clarelli et al. (2013, J. Math. Biol., 66, 1387-1408) to model the formation of cyanobacteria biofilms. After analysing the linear stability of the unique non-trivial equilibrium of the system, we introduce in the model the influence of light and temperature, which are two important factors for the development of a cyanobacteria biofilm. Since the values of the coefficients we use for our simulations are estimated through information found in the literature, some sensitivity and robustness analyses on these parameters are performed. All these elements enable us to control and to validate the model we have already derived and to present some numerical simulations in the 2D and the 3D cases.

We apply the Z-control approach to a generalized predator prey system and consider the specific case of indirect control of the prey population. We derive the associated Z-controlled model and investigate its properties from the point of view of the dynamical systems theory. The key role of the design parameter A. for the successful application of the method is stressed and related to specific dynamical properties of the Z-controlled model. Critical values of the design parameter are also found, delimiting the lambda-range for the effectiveness of the Z-method. Analytical results are then numerically validated by the means of two ecological models: the classical Lotka-Volterra model and a model related to a case study of the wolf wild boar dynamics in the Alta Murgia National Park. Investigations on these models also highlight how the Z-control method acts in respect to different dynamical regimes of the uncontrolled model. (C) 2016 The Authors. Published by Elsevier Inc.

In this paper some recent mathematical models applied to material damage are reviewed. The complexity of damage processes related to Cultural Heritage materials creates the necessity of developing predictive tools in order to monitor and detect surface alterations even before they are visible by naked eyes. The proposed models, elaborated by a research group of the Institute of Applied Mathematics - Research Council of Italy, are based on partial differential equations and well capture the main features of chemical processes (copper corrosion and salt crystallization), which occur on different materials such as stone and copper.

We propose a version of the pure temporal epidemic type aftershock sequences (ETAS) model: the ETAS model with correlated magnitudes. As for the standard case, we assume the Gutenberg-Richter law to be the probability density for the magnitudes of the background events. Instead, the magnitude of the triggered shocks is assumed to be probabilistically dependent on that of the relative mother events. This probabilistic dependence is motivated by some recent works in the literature and by the results of a statistical analysis made on some seismic catalogs [Spassiani and Sebastiani, J. Geophys. Res. 121, 903 (2016)10.1002/2015JB012398]. On the basis of the experimental evidences obtained in the latter paper for the real catalogs, we theoretically derive the probability density function for the magnitudes of the triggered shocks proposed in Spassiani and Sebastiani and there used for the analysis of two simulated catalogs. To this aim, we impose a fundamental condition: averaging over all the magnitudes of the mother events, we must obtain again the Gutenberg-Richter law. This ensures the validity of this law at any event's generation when ignoring past seismicity. The ETAS model with correlated magnitudes is then theoretically analyzed here. In particular, we use the tool of the probability generating function and the Palm theory, in order to derive an approximation of the probability of zero events in a small time interval and to interpret the results in terms of the interevent time between consecutive shocks, the latter being a very useful random variable in the assessment of seismic hazard.

The distribution of the magnitudes of seismic events is generally assumed to be independent on past seismicity. However, by considering events in causal relation, for example, mother-daughter, it seems natural to assume that the magnitude of a daughter event is conditionally dependent on one of the corresponding mother events. In order to find experimental evidence supporting this hypothesis, we analyze different catalogs, both real and simulated, in two different ways. From each catalog, we obtain the law of the magnitude of the triggered events by kernel density. The results obtained show that the distribution density of the magnitude of the triggered events varies with the magnitude of their corresponding mother events. As the intuition suggests, an increase of the magnitude of the mother events induces an increase of the probability of having high values of the magnitude of the triggered events. In addition, we see a statistically significant increasing linear dependence of the magnitude means.

The understanding of water transport in graphene oxide (GO) membranes stands out as a major theoretical problem in graphene research. Notwithstanding the intense efforts devoted to the subject in the recent years, a consolidated picture of water transport in GO membranes is yet to emerge. By performing mesoscale simulations of water transport in ultrathin GO membranes, we show that even small amounts of oxygen functionalities can lead to a dramatic drop of the GO permeability, in line with experimental findings. The coexistence of bulk viscous dissipation and spatially extended molecular friction results in a major decrease of both slip and bulk flow, thereby suppressing the fast water transport regime observed in pristine graphene nanochannels. Inspection of the flow structure reveals an inverted curvature in the near-wall region, which connects smoothly with a parabolic profile in the bulk region. Such inverted curvature is a distinctive signature of the coexistence between single-particle zero-temperature (noiseless) Langevin friction and collective hydrodynamics. The present mesoscopic model with spatially extended friction may offer a computationally efficient tool for future simulations of water transport in nanomaterials. Copyright (C) EPLA, 2016

It is commonly agreed that the most challenging problems in modern science and engineering involve the concurrent and nonlinear interaction of multiple phenomena, acting on a broad and disparate spectrum of scales in space and time. It is also understood that such phenomena lie at the interface between different disciplines, such as physics, chemistry, material science and biology. The multiscale and multi-level nature of these problems commands a paradigm shift in the way they need to be handled, both conceptually and in terms of the corresponding problem-solving computational tools

Hybrid particle-continuum computational frameworks permit the simulation of gas flows by locally adjusting the resolution to the degree of non-equilibrium displayed by the flow in different regions of space and time. In this work, we present a new scheme that couples the direct simulation Monte Carlo (DSMC) with the lattice Boltzmann (LB) method in the limit of isothermal flows. The former handles strong non-equilibrium effects, as they typically occur in the vicinity of solid boundaries, whereas the latter is in charge of the bulk flow, where non-equilibrium can be dealt with perturbatively, i.e. according to Navier-Stokes hydrodynamics. The proposed concurrent multiscale method is applied to the dilute gas Couette flow, showing major computational gains when compared with the full DSMC scenarios. In addition, it is shown that the coupling with LB in the bulk flow can speed up the DSMC treatment of the Knudsen layer with respect to the full DSMC case. In other words, LB acts as a DSMC accelerator.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'. We discuss a unified mesoscale framework (chimaera) for the simulation of complex states of flowing matter across scales of motion. The chimaera framework can deal with each of the three macro-meso-micro levels through suitable 'mutations' of the basic mesoscale formulation. The idea is illustrated through selected simulations of complex micro-and nanoscale flows.

This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'. Hybrid particle-continuum computational frameworks permit the simulation of gas flows by locally adjusting the resolution to the degree of non-equilibrium displayed by the flow in different regions of space and time. In this work, we present a new scheme that couples the direct simulation Monte Carlo (DSMC) with the lattice Boltzmann (LB) method in the limit of isothermal flows. The former handles strong non-equilibrium effects, as they typically occur in the vicinity of solid boundaries, whereas the latter is in charge of the bulk flow, where non-equilibrium can be dealt with perturbatively, i.e. according to Navier-Stokes hydrodynamics. The proposed concurrent multiscale method is applied to the dilute gas Couette flow, showing major computational gains when compared with the full DSMC scenarios. In addition, it is shown that the coupling with LB in the bulk flow can speed up the DSMC treatment of the Knudsen layer with respect to the full DSMC case. In other words, LB acts as a DSMC accelerator.

We present an axisymmetric lattice Boltzmann model based on the Kupershtokh et al. multiphase model that is capable of solving liquid-gas density ratios up to 10(3). Appropriate source terms are added to the lattice Boltzmann evolution equation to fully recover the axisymrnetric multiphase conservation equations. We validate the model by showing that a stationary droplet obeys the Young-Laplace law, comparing the second oscillation mode of a droplet with respect to an analytical solution and showing correct mass conservation of a propagating density wave.

We present the formulation of a kinetic mapping scheme between the Direct Simulation Monte Carlo (DSMC) and the Lattice Boltzmann Method (LBM) which is at the basis of the hybrid model used to couple the two methods in view of efficiently and accurately simulate isothermal flows characterized by variable rarefaction effects. Owing to the kinetic nature of the LBM, the procedure we propose ensures to accurately couple DSMC and LBM at a larger Kn number than usually done in traditional hybrid DSMC--Navier-Stokes equation models. We show the main steps of the mapping algorithm and illustrate details of the implementation. Good agreement is found between the moments of the single particle distribution function as obtained from the mapping scheme and from independent LBM or DSMC simulations at the grid nodes where the coupling is imposed. We also show results on the application of the hybrid scheme based on a simpler mapping scheme for plane Poiseuille flow at finite Kn number. Potential gains in the computational efficiency assured by the application of the coupling scheme are estimated for the same flow.

We introduce a new particle shape which shows preferential rotation in three dimensional homogeneous isotropic turbulence. We call these particles chiral dipoles because they consist of a rod with two helices of opposite handedness, one at each end. 3D printing is used to fabricate these particles with a length in the inertial range and their rotations are tracked in a turbulent flow between oscillating grids. High aspect ratio chiral dipoles preferentially align with their long axis along the extensional eigenvectors of the strain rate tensor, and the helical ends respond to the extensional strain rate with a mean spinning rate that is nonzero. We use Stokesian dynamics simulations of chiral dipoles in pure strain flow to quantify the dependence of spinning on particle shape. Based on the known response to pure strain, we build a model that gives the spinning rate of small chiral dipoles using velocity gradients along Lagrangian trajectories from high resolution direct numerical simulations. The statistics of chiral dipole spinning determined with this model show surprisingly good agreement with the measured spinning of much larger chiral dipoles in the experiments.

We investigate the influence of shear on the gravitational settling of heavy inertial particles in homogeneous shear turbulence (HST). In addition to the well-known enhanced settling velocity, observed for heavy inertial particles in homogeneous isotropic turbulence (HIT), a horizontal drift velocity is also observed in the shearing direction due to the presence of a nonzero mean vorticity (introducing symmetry breaking due to the mean shear). This drift velocity is due to the combination of shear, gravity, and turbulence, and all three of these elements are needed for this effect to occur. We extend the mechanism responsible for the enhanced settling velocity in HIT to the case of HST. Two separate regimes are observed, characterized by positive or negative drift velocity, depending on the particle settling velocity.

Planktonic copepods are small crustaceans that have the ability to swim by quick powerful jumps. Such an aptness is used to escape from high shear regions, which may be caused either by flow perturbations, produced by a large predator (i.e., fish larvae), or by the inherent highly turbulent dynamics of the ocean. Through a combined experimental and numerical study, we investigate the impact of jumping behavior on the small-scale patchiness of copepods in a turbulent environment. Recorded velocity tracks of copepods displaying escape response jumps in still water are here used to define and tune a Lagrangian copepod (LC) model. The model is further employed to simulate the behavior of thousands of copepods in a fully developed hydrodynamic turbulent flow obtained by direct numerical simulation of the Navier-Stokes equations. First, we show that the LC velocity statistics is in qualitative agreement with available experimental observations of copepods in turbulence. Second, we quantify the clustering of LC, via the fractal dimension D2. We show that D2 can be as low as ~2.3 and that it critically depends on the shear-rate sensitivity of the proposed LC model, in particular it exhibits a minimum in a narrow range of shear-rate values. We further investigate the effect of jump intensity, jump orientation, and geometrical aspect ratio of the copepods on the small-scale spatial distribution. At last, possible ecological implications of the observed clustering on encounter rates and mating success are discussed.

We investigate the non-equilibrium hydrodynamic effects on the reactivity of a nanoporous catalytic sample. Numerical simulations using the Lattice Boltzmann Method (LBM) show that non-equilibrium effects enhance the reactivity of the porous sample, in agreement with theoretical predictions [1]. In addition, we provide a quantitative assessment of the reactivity in terms of the thickness of the reactive layer inside the nanoporous catalytic sample. Such an assessment constitutes a first step towards integrated simulations encompassing nanoscale reactivity and transport coefficients within a macroscale description of experimental relevance.

The effects of compressibility on Rayleigh-Taylor instability (RTI) are investigated by inspecting the interplay between thermodynamic and hydrodynamic nonequilibrium phenomena (TNE, HNE, respectively) via a discrete Boltzmann model. Two effective approaches are presented, one tracking the evolution of the local TNE effects and the other focusing on the evolution of the mean temperature of the fluid, to track the complex interfaces separating the bubble and the spike regions of the flow. It is found that both the compressibility effects and the global TNE intensity show opposite trends in the initial and the later stages of the RTI. Compressibility delays the initial stage of RTI and accelerates the later stage. Meanwhile, the TNE characteristics are generally enhanced by the compressibility, especially in the later stage. The global or mean thermodynamic nonequilibrium indicators provide physical criteria to discriminate between the two stages of the RTI.