In this paper we analyze least recently used (LRU) caches operating under the shot noise requests model (SNM). The SNM was recently proposed in [S. Traverso et al., ACM Comput. Comm. Rev., 43 (2013), pp. 5-12] to better capture the main characteristics of today's video on demand traffic. We investigate the validity of Che's approximation [H. Che, Y. Tung, and Z. Wang, IEEE J. Selected Areas Commun., 20 (2002), pp. 1305-1314] through an asymptotic analysis of the cache eviction time. In particular, we provide a law of large numbers, a large deviation principle, and a central limit theorem for the cache eviction time, as the cache size grows large. Finally, we derive upper and lower bounds for the "hit" probability in tandem networks of caches under Che's approximation.
caching systems
performance evaluation
asymptotic analysis
We provide a Clark-Ocone formula for square-integrable functionals of a general temporal point process satisfying only a mild moment condition, generalizing known results on the Poisson space. Some classical applications are given, namely a deviation bound and the construction of a hedging portfolio in a pure-jump market model. As a more modern application, we provide a bound on the total variation distance between two temporal point processes, improving in some sense a recent result in this direction.
Clark-Ocone formula
point processes
Malliavin calculus
conditional intensity
deviation inequalities
option hedging
In this work a two level heuristic algorithm is described for a nearly real-time multi-vehicle many-to-many Dial-A-Ride Problem (DARP). This algorithm is ready to support a Demand Responsive Transportation System in which we face the problem of quickly evaluate a good-quality schedule for the vehicles and provide fast response to the users. The insertion heuristic is double dynamic nearly real-time and the objective function is to minimize the variance between the requested and scheduled time of pickup and delivery. In the first level, after a customer web-request, the heuristic returns an answer about the possibility to insert the request into the accepted reservations, and therefore in a vehicle schedule, or reject the request. In the second level, during the time elapsed between a request and the following, and after a reshuffling of the order of the incoming accepted requests, the same heuristic works for the whole set of accepted requests, trying to optimize the solution. We intensively tested the algorithm with a requests-generating software that has allowed us to show the competitive advantage of this web-based architecture.
Vehicle routing; Heuristics; Dial-A-Ride; Transportation planning; Public transport
The geometry of passive tracer trajectories is studied in two different types of rotating turbulent flows; rotating Rayleigh-Bénard convection (RBC; experiments and direct numerical simulations) and rotating electromagnetically forced turbulence (EFT; experiments). This geometry is fully described by the curvature and torsion of trajectories, and from these geometrical quantities we can subtract information on the typical flow structures at different rotation rates. Previous studies, focusing on nonrotating homogeneous isotropic turbulence (HIT), show that the probability density functions (PDFs) of curvature and torsion reveal pronounced power laws. However, the set-ups studied here involve inhomogeneous turbulence, and in RBC the flow near the horizontal plates is definitely anisotropic. We investigate how the typical shapes of the curvature and torsion PDFs, including the pronounced scaling laws, are influenced by this level of anisotropy and inhomogeneity and how this effect changes with rotation. A first effect of rotation is observed as a shift of the curvature and torsion PDFs towards higher values in the case of RBC and towards lower values in the case of EFT. This shift is related to the length scale of typical vortical structures that decreases with rotation in RBC, but increases with rotation in EFT, explaining the opposite shifts of the curvature (and torsion) PDFs. A second remarkable observation is that in RBC the HIT scaling laws are always recovered, as long as the boundary layer (BL) is excluded. This suggests that these scaling laws are very robust and hold as long as we measure in the turbulent bulk. In the BL of the RBC cell, however, the scaling deviates from the HIT prediction for lower rotation rates. This scaling behavior is found to be consistent with the coupling between the boundary layer dynamics and the bulk flow, which changes under rotation. In particular, it is found that the active coupling of the Ekman-type boundary layer with the bulk flow suppresses anisotropy in the BL region for increasing rotation rates.
This study shows the land cover mapping accuracy retrievable by the TASI-600 thermal airborne multispectral sensor and describes some of the classification results tested on the thermal preprocessed data for a rural area. In the paper is provided an overview of the principal TASI-600 characteristics, i.e. 32 spectral bands in the 8.0-11.5 ?m spectral range, and land cover classification performances. A full assessment of the TASI-600 spectral bands has been also obtained by ranking them in order to understanding their role in land cover classification. Results accuracies have been validated using available ground truth. The study highlights that the new generation of multi/hyperspectral thermal sensors opens up interesting opportunities for accurate land cover classification.
Classification accuracies
Land cover mapping
Multispectral thermal data
TASI-600
We propose a mathematical model for the transport of DNA plasmids from the extracellular matrix up to the cell nucleus. The model couples two phenomena: the electroporation process, describing the cell membrane permeabilization to plasmids and the intracellular transport enhanced by the presence of microtubules. Numerical simulations of cells with arbitrary geometry, in 2D and 3D, and a network of microtubules show numerically the importance of the microtubules and the electroporation on the effectiveness of the DNA transfection, as observed by previous biological data. The paper proposes efficient numerical tools for forthcoming optimized procedures of cell transfection.
Finite difference method in irregular domains
Intracellular transport
Mathematical biology
Microtubules
Plasmids transport
1/c expansion of nonminimally coupled curvature-matter gravity models and constraints from planetary precession
Riccardo March
;
Jorge Paramos
;
Orfeu Bertolami
;
Simone Dell'Agnello
The effects of a nonminimally coupled curvature-matter model of gravity on a perturbed Minkowski metric are presented. The action functional of the model involves two functions $f^1(R)$ and $f^2(R)$ of the Ricci scalar curvature $R$.
This work expands upon previous results, extending the framework developed there to compute corrections up to order $O\left(1\slash c^4\right)$ of the 00 component of the metric tensor.
It is shown that additional contributions arise due to both the non-linear form $f^1(R)$ and the nonminimal coupling $f^2(R)$, including exponential contributions that cannot be expressed as an expansion in powers of $1/r$. Some possible experimental implications are assessed with application to perihelion precession.
The theoretical understanding of Gartner's "hype curve" is an interesting open question in deciding the strategic actions to adopt in presence of an incoming technology. In order to describe the hype behaviour quantitatively, we propose a mathematical approach based on a rate equation, similar to that used to describe quantum level transitions. The model is able to describe the hype curve evolution in many relevant conditions, which can be associated to various market parameters. Different hype curves, describing the time evolution of a new technology market penetration, are then obtained within a single coherent mathematical approach. We have also used our theoretical model to describe the time evolution of the number of scientific publications in different fields of scientific research. Data are well described by our model, so we present a statistical analysis and forecasting potentiality of our approach. We note that the hype peak of inflated expectations is very smooth in the case of scientific publications, probably due to the high level of awareness and the deep preliminary understanding which is necessary to carry on a research project. Our model is anyway flexible enough to describe many patterns of increasing interest on a new idea, leading to a hype behaviour or other time evolution.
Hype cycle; Hype-type mathematical approach; Rate equation; Technological development analysis
Hypomorphic mutations in DNA-methyltransferase DNMT3B cause majority of the rare disorder Immunodeficiency, Centromere instability and Facial anomalies syndrome cases (ICF1). By unspecified mechanisms, mutant-DNMT3B interferes with lymphoid-specific pathways resulting in immune response defects. Interestingly, recent findings report that DNMT3B shapes intragenic CpG-methylation of highly-transcribed genes. However, how the DNMT3B-dependent epigenetic network modulates transcription and whether ICF1-specific mutations impair this process remains unknown. We performed a transcriptomic and epigenomic study in patient-derived B-cell lines to investigate the genome-scale effects of DNMT3B dysfunction. We highlighted that altered intragenic CpG-methylation impairs multiple aspects of transcriptional regulation, like alternative TSS usage, antisense transcription and exon splicing. These defects preferentially associate with changes of intragenic H3K4me3 and at lesser extent of H3K27me3 and H3K36me3. In addition, we highlighted a novel DNMT3B activity in modulating the self-regulatory circuit of sense-antisense pairs and the exon skipping during alternative splicing, through interacting with RNA molecules. Strikingly, altered transcription affects disease relevant genes, as for instance the memory-B cell marker CD27 and PTPRC genes, providing us with biological insights into the ICF1-syndrome pathogenesis. Our genome-scale approach sheds light on the mechanisms still poorly understood of the intragenic function of DNMT3B and DNA methylation in gene expression regulation.
Active research in nanotechnology contemplates the use of nanomaterials for environmental engineering applications. However, a primary challenge is understanding the effects of nanomaterial properties on industrial device performance and translating unique nanoscale properties to the macroscale. One emerging example consists of graphene oxide (GO) membranes for separation processes. Thus, here we investigate how individual GO properties can impact GO membrane characteristics and water permeability. GO chemistry and morphology were controlled with easy-to-implement photoreduction and sonication techniques and were quantitatively correlated, offering a valuable tool for accelerating characterization. Chemical GO modification allows for fine control of GO oxidation state, allowing control of GO architectural laminate (GOAL) spacing and permeability. Water permeability was measured for eight GOALs characterized by different GOAL chemistry and morphology and indicates that GOAL nanochannel height dictates water transport. The experimental outputs were corroborated with mesoscale water transport simulations of relatively large domains (thousands of square nanometers) and indicate a no-slip Darcy-like behavior inside the GOAL nanochannels. The experimental and simulation evidence presented in this study helps create a clearer picture of water transport in GOAL, and can be used to rationally design more effective and efficient GO membranes.
Enhanced computational performance of the lattice Boltzmann model for simulating micron- and submicron-size particle flows and non-Newtonian fluid flows
Basagaoglu Hakan
;
Harwell John R
;
Hoa Nguyen
;
Succi Sauro
Significant improvements in the computational performance of the lattice-Boltzmann (LB) model, coded in FORTRAN90, were achieved through application of enhancement techniques. Applied techniques include optimization of array memory layouts, data structure simplification, random number generation outside the simulation thread(s), code parallelization via OpenMP, and intra- and inter-timestep task pipelining. Effectiveness of these optimization techniques was measured on three benchmark problems: (i) transient flow of multiple particles in a Newtonian fluid in a heterogeneous fractured porous domain, (ii) thermal fluctuation of the fluid at the sub-micron scale and the resultant Brownian motion of a particle, and (iii) non-Newtonian fluid flow in a smooth-walled channel. Application of the aforementioned optimization techniques resulted in an average 21x performance improvement, which could significantly enhance practical uses of the LB models in diverse applications, focusing on the fate and transport of nano-size or micron-size particles in non-Newtonian fluids. (C) 2016 Elsevier B.V. All rights reserved.
Computational methods in fluid dynamics
Hydrodynamics
Lattice-Boltzmann
FORTRAN
Optimization
Computational performance of SequenceL coding of the lattice Boltzmann method for multi-particle flow simulations
Basagaoglu Hakan
;
Blount Justin
;
Blount Jarred
;
Nelson Bryant
;
Succi Sauro
;
Westhart Phil M
;
Harwell John R
This paper reports, for the first time, the computational performance of SequenceL for mesoscale simulations of large numbers of particles in a microfluidic device via the lattice-Boltzmann method. The performance of SequenceL simulations was assessed against the optimized serial and parallelized (via OpenMP directives) FORTRAN90 simulations. At present, OpenMP directives were not included in interparticle and particle-wall repulsive (steric) interaction calculations due to difficulties that arose from inter-iteration dependencies between consecutive iterations of the do-loops. SequenceL simulations, on the other hand, relied on built-in automatic parallelism. Under these conditions, numerical simulations revealed that the parallelized FORTRAN90 outran the performance of SequenceL by a factor of 2.5 or more when the number of particles was 100 or less. SequenceL, however, outran the performance of the parallelized FORTRAN90 by a factor of 1.3 when the number of particles was 300. Our results show that when the number of particles increased by 30-fold, the computational time of SequenceL simulations increased linearly by a factor of 1.5, as compared to a 3.2-fold increase in serial and a 7.7-fold increase in parallelized FORTRAN90 simulations. Considering SequenceL's efficient built-in parallelism that led to a relatively small increase in computational time with increased number of particles, it could be a promising programming language for computationally-efficient mesoscale simulations of large numbers of particles in microfluidic experiments. (C) 2016 Elsevier B.V. All rights reserved.
Computational methods in fluid dynamics
Hydrodynamics
Lattice-Boltzmann
SequenceL
In this paper, an alternative approach to implement initial and boundary conditions in the lattice Boltzmann method is presented. The main idea is to approximate the nonequilibrium component of distribution functions as a third-order power series in the lattice velocities and formulate a procedure to determine boundary node distributions by using fluid variables, consistent with such an expansion. The velocity shift associated with the body force effects is included in this scheme, along with an approximation to determine the mass density in complex geometries. Different strategies based on the present scheme are developed to implement velocity and pressure conditions for arbitrarily shaped boundaries, using the D2Q9, D3Q15, D3Q19 and D3Q27 lattices, in two and three space dimensions, respectively. The proposed treatment is tested against several well-established problems, showing second-order spatial accuracy and often improved behavior as compared to various existing methods, with no appreciable computational overhead.
We revisit the integer lattice (IL) method to numerically solve the Vlasov-Poisson equations, and show that a slight variant of the method is a very easy, viable, and efficient numerical approach to study the dynamics of self-gravitating, collisionless systems. The distribution function lives in a discretized lattice phase-space, and each time-step in the simulation corresponds to a simple permutation of the lattice sites. Hence, the method is Lagrangian, conservative, and fully time-reversible. IL complements other existing methods, such as N-body/ particle mesh (computationally efficient, but affected by Monte Carlo sampling noise and two-body relaxation) and finite volume (FV) direct integration schemes (expensive, accurate but diffusive). We also present improvements to the FV scheme, using a moving-mesh approach inspired by IL, to reduce numerical diffusion and the time-step criterion. Being a direct integration scheme like FV, IL is memory limited (memory requirement for a full 3D problem scales as N-6, where N is the resolution per linear phase- space dimension). However, we describe a new technique for achieving N-4 scaling. The method offers promise for investigating the full 6D phase- space of collisionless systems of stars and dark matter.
gravitation
methods: numerical
stars: kinematics and dynamics
galaxies: kinematics and dynamics
dark matter
Roughness as a Route to the Ultimate Regime of Thermal Convection
Toppaladoddi Srikanth
;
Succi Sauro
;
Wettlaufer John S
We use highly resolved numerical simulations to study turbulent Rayleigh-Benard convection in a cell with sinusoidally rough upper and lower surfaces in two dimensions for Pr = 1 and Ra = [4 x 10(6), 3 x 10(9)]. By varying the wavelength. at a fixed amplitude, we find an optimal wavelength lambda(opt) for which the Nusselt-Rayleigh scaling relation is (Nu - 1 proportional to Ra-0.483), maximizing the heat flux. This is consistent with the upper bound of Goluskin and Doering [J. Fluid Mech. 804, 370 (2016)] who prove that Nu can grow no faster than O(Ra-1/2) as Ra -> infinity, and thus with the concept that roughness facilitates the attainment of the so-called ultimate regime. Our data nearly achieve the largest growth rate permitted by the bound. When lambda << lambda(opt) and lambda >> lambda(opt), the planar case is recovered, demonstrating how controlling the wall geometry manipulates the interaction between the boundary layers and the core flow. Finally, for each Ra, we choose the maximum Nu among all., thus optimizing over all lambda, to find Nu(opt) - 1 = 0.01xRa(0.444).
Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved spaces: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature-induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.
Hybrid lattice Boltzmann method on overlapping grids
Di Ilio G
;
Chiappini D
;
Ubertini S
;
Bella G
;
Succi S
In this work, a hybrid lattice Boltzmann method (HLBM) is proposed, where the standard lattice Boltzmann implementation based on the Bhatnagar-Gross-Krook (LBGK) approximation is combined together with an unstructured finite-volume lattice Boltzmann model. The method is constructed on an overlapping grid system, which allows the coexistence of a uniform lattice nodes spacing and a coordinate-free lattice structure. The natural adaptivity of the hybrid grid system makes the method particularly suitable to handle problems involving complex geometries. Moreover, the provided scheme ensures a high-accuracy solution near walls, given the capability of the unstructured submodel of achieving the desired level of refinement in a very flexible way. For these reasons, the HLBM represents a prospective tool for solving multiscale problems. The proposed method is here applied to the benchmark problem of a two-dimensional flow past a circular cylinder for a wide range of Reynolds numbers and its numerical performances are measured and compared with the standard LBGK ones.
Growth in static and controlled environments such as a Petri dish can be used to study the spatial population dynamics of microorganisms. However, natural populations such as marine microbes experience fluid advection and often grow up in heterogeneous environments. We investigate a generalized Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation describing single species population subject to a constant flow field and quenched random spatially inhomogeneous growth rates with a fertile overall growth condition. We analytically and numerically demonstrate that the non-equilibrium steady-state population density develops a flow-driven striation pattern. The striations are highly asymmetric with a longitudinal correlation length that diverges linearly with the flow speed and a transverse correlation length that approaches a finite velocity-independent value. Linear response theory is developed to study the statistics of the steady states. Theoretical predictions show excellent agreement with the numerical steady states of the generalized FKPP equation obtained from Lattice Boltzmann simulations. These findings suggest that, although the growth disorder can be spatially uncorrelated, correlated population structures with striations emerge naturally at sufficiently strong advection. (C) 2016 Elsevier B.V. All rights reserved.
Population dynamics
Spatially quenched disorder
Generalized FKPP equation
Pattern formation
A challenging task in the management of Protected Areas is the conservation of natural habitats and native endangered species through the optimization of control strategies for invasive plant or animal species, typically competing for the use of resources in a fragmented habitat [1]. We review two cases of control strategies on the wolf-wild boar populations in a Southern Italy Protected Area belonging to the Natura 2000 network [2,3]. The challenge for the regional authorities is to plan conservation policies able to maintain the population of wolves while limiting the presence of wild boars, here considered invasive because of their harmfulness on cultivated areas. The first strategy [2] considers the impact of control policies on predator-prey dynamics in fragmented habitats by simulating different dynamical scenarios theoretically analysed with the aggregation method. The key warning from the model is that a very careful combination of control - through proper planning programs - and migration processes among patches of habitats - through the existing suitable ecological corridors - must be used in order to limit the wild-boar population while preserving wolves from extinction. The second strategy has been developed to apply the Z-control approach to a generalized predator-prey system [3]. It considers the specific case of indirect control of the prey (invasive) population. The key role of the model design parameter is stressed and the critical values of the design parameter are found, delimiting the parameter range for the successful application of the Z-method. A further development is the optimization of a control strategy by taking into account the spatio-temporal data related to the control problem of an invasive species over a wide natural protected area. That approach will be applied to the Alta Murgia National Park, where a EU LIFE+ project is underway to eradicate Ailanthus altissima, included in the list of the most invasive alien plant species in Europe causing serious damages both in protected and urban areas [4]. The Alta Murgia National Park is one of the study site of an on-going H2020 project (ECOPOTENTIAL).
This work has been carried out within the H2020 project `ECOPOTENTIAL: Improving Future Ecosystem Benefits Through Earth Observations', coordinated by CNR-IGG (http://www.ecopotential-project.eu). The project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 641762.
Bibliography
1)Baker, C. M., Target the Source: Optimal Spatiotemporal Resource Allocation for Invasive Species Control, CONS. LETTERS, pp 1-8, 2016, doi: 10.1111/conl.12236
2)Lacitignola, D.; Diele, F.; Marangi, C., Dynamical scenarios from a two-patch predator-prey system with human control - Implications for the conservation of the wolf in the Alta Murgia National Park ECOLOGICAL MODELLING, Vol. 316, pp 28-40, 2015,
doi: 10.1016/j.ecolmodel.2015.07.027
3)Lacitignola, D.; Diele, F.; Marangi, C.; Provenzale A., On the dynamics of a generalized predator-prey system with Z-type control, MATHEMATICAL BIOSCIENCES, vol. 280, pp 10-23, 2016, doi: 10.1016/j.mbs.2016.07.011
4)Casella F., Vurro M. , Ailanthus altissima (tree of heaven): Spread and harmfulness in a case-study urban area, Arboricultural Journal: The International Journal of Urban Forestry, 35(3), pp 172-181, 2013, doi: 10.1080/03071375.2013.852352
optimal control
invasive species
mathematical model
protected areas