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.

We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range, and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.

It is shown that the single relaxation time (SRT) version of the Lattice Boltzmann (LB) equation permits to compute the permeability of Darcy's flows in porous media within a few percent accuracy. This stands in contrast with previous claims of inaccuracy, which we relate to the lack of recognition of the physical dependence of the permeability on the Knudsen number.

In the current study, a direct-forcing immersed boundary-non-Newtonian lattice Boltzmann method (IB-NLBM) is developed to investigate the sedimentation and interaction of particles in shear-thinning and shear-thickening fluids. In the proposed IB-NLBM, the non-linear mechanics of non-Newtonian particulate flows is detected by combination of the most desirable features of immersed boundary and lattice Boltzmann methods. The noticeable roles of non-Newtonian behavior on particle motion, settling velocity and generalized Reynolds number are investigated by simulating benchmark problem of one-particle sedimentation under the same generalized Archimedes number. The effects of extra force due to added accelerated mass are analyzed on the particle motion which have a significant impact on shear-thinning fluids. For the first time, the phenomena of interaction among the particles, such as Drafting, Kissing, and Tumbling in non-Newtonian fluids are investigated by simulation of two-particle sedimentation and twelve-particle sedimentation. The results show that increasing the shear-thickening behavior of fluid leads to a significant increase in the kissing time. Moreover, the transverse position of particles for shear-thinning fluids during the tumbling interval is different from Newtonian and the shear-thickening fluids. The present non-Newtonian particulate study can be applied in several industrial and scientific applications, like the non-Newtonian sedimentation behavior of particles in food industrial and biological fluids.

Questo volume, che esce a ridosso del centenario della morte di Eugenio Elia Levi (1883-1917), completa la ricerca svolta dai curatori per rintracciare e mettere a disposizione degli storici tutta la documentazione relativa a questo grande matematico. Dopo un primo volume, pubblicato da Pristem-Bocconi e basato sull'epistolario del Levi, qui i curatori aggiungono una grande quantità di documenti ed informazioni che, anche nella parte già edita, sono difficili da identificare e reperire. Di particolare interesse sono quelli relativi all'impegno didattico profuso da E.E Levi nell'ambito della Mathesis e in collaborazione con Giuseppe Lombardo Radice, alla sua vita accademica, alla sua intensa attività politica a favore dell'interventismo e alle circostanze che hanno portato alla sua tragica morte e che per la prima volta vengono ricostruite in modo dettagliato.