In this paper, we demonstrate that vertically vibrating a plate in a cornstarch suspension causes the suspension to vigorously ratchet up the plate. We show that this is a necessary consequence of the fact that cornstarch in water is shear thickening: when the plate moves up it opposes gravity and so the fluid stiffens; when it moves down it works with gravity and so the fluid flows. This produces asymmetric ratcheting that opposes gravity. We find several unusual states that result in this simple experimental system, and we reproduce the essential effect in two different numerical simulations. (C) 2015 AIP Publishing LLC.
It is shown that the combination of generalized Van der Waals equations of state with high-order discrete velocity lattices, permits to simulate the dynamics of liquid droplets at air-water density ratios, with very moderate levels of spurious currents near the droplet interface. Satisfactory agreement with experimental data on droplet collisions at density ratios of order thousand is reported.
Novel risk predictor for thrombus deposition in abdominal aortic aneurysms
Nestola M G C
;
Gizzi A
;
Cherubini C
;
Filippi S
;
Succi S
The identification of the basic mechanisms responsible for cardiovascular diseases stands as one of the most challenging problems in modern medical research including various mechanisms which encompass a broad spectrum of space and time scales. Major implications for clinical practice and pre-emptive medicine rely on the onset and development of intraluminal thrombus in which effective clinical therapies require synthetic risk predictors/indicators capable of informing real-time decision-making protocols. In the present contribution, two novel hemodynamics synthetic indicators, based on a three-band decomposition (TBD) of the shear stress signal, are introduced. Extensive fluid-structure computer simulations of patient-specific scenarios confirm the enhanced risk-prediction capabilities of the TBD indicators. In particular, they permit a quantitative and accurate localization of the most likely thrombus deposition in realistic aortic geometries, where previous indicators would predict healthy operation. The proposed methodology is also shown to provide additional information and discrimination criteria on other factors of major clinical relevance, such as the size of the aneurysm. Copyright (C) EPLA, 2015
We study the Couette flow of a quasi-2d soft-glassy material in a Hele-Shaw geometry. The material is chosen to be above the jamming point, where a yield stress sigma(Upsilon) emerges, below which the material deforms elastically and above which it flows like a complex fluid according to a Herschel-Bulkley (HB) rheology. Simultaneously, the effect of the confining plates is modelled as an effective linear friction law, while the walls aside the Hele-Shaw cell are sufficiently close to each other to allow visible cooperativity effects in the velocity profiles (Goyon et al., 2008). The effects of cooperativity are parametrized with a steady-state diffusion-relaxation equation for the fluidity field f = gamma over dot/sigma , defined as the ratio between shear rate.j/and shear stress cr. For particular rheological flow-curves (Bingham fluids), the problem is tackled analytically: we explore the two regimes sigma >> sigma(Upsilon) and alpha approximate to sigma(Upsilon) and quantify the effect of the extra localisation induced by the wall friction. Other rheo-thinning fluids are explored with the help of numerical simulations based on lattice Boltzmann models, revealing a robustness of the analytical findings. Synergies and comparisons with other existing works in the literature (Barry et al., 2011) are also discussed. (C) 2015 Elsevier B.V. All rights reserved.
We propose numerical simulations of viscoelastic fluids based on a hybrid algorithm combining Lattice-Boltzmann models (LBM) and Finite Differences (FD) schemes, the former used to model the macroscopic hydrodynamic equations, and the latter used to model the polymer dynamics. The kinetics of the polymers is introduced using constitutive equations for viscoelastic fluids with finitely extensible non-linear elastic dumbbells with Peterlin's closure (FENE-P). The numerical model is first benchmarked by characterizing the rheological behavior of dilute homogeneous solutions in various configurations, including steady shear, elongational flows, transient shear and oscillatory flows. As an upgrade of complexity, we study the model in presence of non-ideal multicomponent interfaces, where immiscibility is introduced in the LBM description using the "Shan-Chen" interaction model. The problem of a confined viscoelastic (Newtonian) droplet in a Newtonian (viscoelastic) matrix under simple shear is investigated and numerical results are compared with the predictions of various theoretical models. The proposed numerical simulations explore problems where the capabilities of LBM were never quantified before. (C) 2015 Elsevier Inc. All rights reserved.
Polymers
Lattice Boltzmann models
Binary liquids
Droplet deformation and orientation
Mathematical modelling of viscoelastic flows
The ECOPOTENTIAL project focuses its activities and pilot actions on a targeted set of internationally recognised protected areas (PA) in Europe, European Territories and beyond, including mountain, arid and semi-arid, and coastal and marine ecosystems. Building on the knowledge gained in individual PAs, the ECOPOTENTIAL project addresses cross-scale ecological interactions and landscape-ecosystem dynamics at regional to continental scales, using geostatistical methods and the emerging novel approaches in Macrosystems Ecology, which is addressing long-term and large-scale ecological challenges. ECOPOTENTIAL addresses the entire chain of ecosystem-related services, by (a) developing ecosystem data services, with special emphasis on Copernicus services; (b) implementing model output services to distribute the results of the modelling activities; and (c) estimating current and future ecosystem services and benefits, combining ecosystem functions (supply) with beneficiaries needs (demand). In ECOPOTENTIAL all data, model results and acquired knowledge are made available on common and open platforms, coherent with the Global Earth Observation System of Systems (GEOSS) data sharing principles and fully interoperable with the GEOSS Common Infrastructure (GCI).
Protected Area (PA)
Global Earth Observation System of Systems (GEOSS)
GEOSS Common Infrastructure (GCI)
Sap transport in trees has long fascinated scientists, and a vast literature exists on experimental and modelling studies of trees during the growing season when large negative stem pressures are generated by transpiration from leaves. Much less attention has been paid to winter months when trees are largely dormant but nonetheless continue to exhibit interesting flow behaviour. A prime example is sap exudation, which refers to the peculiar ability of sugar maple (Acer saccharum) and related species to generate positive stem pressure while in a leafless state. Experiments demonstrate that ambient temperatures must oscillate about the freezing point before significantly heightened stem pressures are observed, but the precise causes of exudation remain unresolved. The prevailing hypothesis attributes exudation to a physical process combining freeze-thaw and osmosis, which has some support from experimental studies but remains a subject of active debate. We address this knowledge gap by developing the first mathematical model for exudation, while also introducing several essential modifications to this hypothesis. We derive a multiscale model consisting of a nonlinear system of differential equations governing phase change and transport within wood cells, coupled to a suitably homogenized equation for temperature on the macroscale. Numerical simulations yield stem pressures that are consistent with experiments and provide convincing evidence that a purely physical mechanism is capable of capturing exudation.
tree sap exudation
sugar maple
multiphase flow and transport
phase change
differential equations
periodic homogenization
Mechanotransduction map: simulation model, molecular pathway, gene set
Dent Jennifer E
;
Devescovi Valentina
;
Li Han
;
Di Lena Pietro
;
Lu Youtao
;
Liu Yuanhua
;
Nardini Christine
Results: We here present a molecular map of mechanotransduction, built in CellDesigner to warrant that maximum information is embedded in a compact network format. To validate the map's necessity we tested its redundancy in comparison with existing pathways, and to estimate its sufficiency, we quantified its ability to reproduce biological events with dynamic simulations, using Signaling Petri Networks.
Motivation: Mechanotransduction-the ability to output a biochemical signal from a mechanical input-is related to the initiation and progression of a broad spectrum of molecular events. Yet, the characterization of mechanotransduction lacks some of the most basic tools as, for instance, it can hardly be recognized by enrichment analysis tools, nor could we find any pathway representation. This greatly limits computational testing and hypothesis generation on mechanotransduction biological relevance and involvement in disease or physiological mechanisms.
The dynamics of colloidal intrusions in liquid crystals: A simulation perspective
Foffano G
;
Lintuvuori J S
;
Tiribocchi A
;
Marenduzzo D
Dispersing colloidal particles into liquid crystals provides a promising avenue to build a novel class of materials, with potential applications, among others, as photonic crystals, biosensors, metamaterials and new generation liquid crystal devices. Understanding the physics and dynamical properties of such composite materials is then of high-technological relevance; it also provides a remarkable challenge from a fundamental science point of view due to the intricacies of the hydrodynamic equations governing their dynamical evolution. Here, we provide an overview of our current theoretical understanding of the dynamical and hydrodynamic properties of colloid-liquid crystal composites, focussing on the results obtained from computer simulations; no or very limited previous knowledge of the field of liquid crystals is assumed. While our main emphasis is on the dynamics, we also review a selection of equilibrium results and simulations to provide the necessary background. We start by describing what we know about the simplest possible problem: that of a single particle in a nematic, or cholesteric, liquid crystal. We then consider two particles, and review the conditions which lead to the formation of a dimer; we then again focus on dynamical problems. Finally, we turn to the more complicated case of a dispersion, reviewing here simulations motivated by optical tweezer and rheological experiments. We close by making a list of some of the many open problems in this rapidly developing research field.
Switching hydrodynamics in liquid crystal devices: A simulation perspective
Tiribocchi A
;
Henrich O
;
Lintuvuori J S
;
Marenduzzo D
In liquid crystal devices it is important to understand the physics underlying their switching between different states, which is usually achieved by applying or removing an electric field. Flow is known to be a key determinant of the timescales and pathways of the switching kinetics. Incorporating hydrodynamic effects into theories for liquid crystal devices is therefore important; however this is also highly non-trivial, and typically requires the use of accurate numerical methods. Here, we review some recent advances in our theoretical understanding of the dynamics of switching in liquid crystal devices, mainly gained through computer simulations. These results, as we shall show, uncover interesting new physics, and may be important for future applications.
Scalar ? 4 field theory for active-particle phase separation
R Wittkowski
;
A Tiribocchi
;
J Stenhammar
;
R Allen
;
D Marenduzzo
;
M E Cates
Recent theories predict phase separation among orientationally disordered active particles whose propulsion speed decreases rapidly enough with density. Coarse-grained models of this process show time-reversal symmetry (detailed balance) to be restored for uniform states, but broken by gradient terms; hence, detailed-balance violation is strongly coupled to interfacial phenomena. To explore the subtle generic physics resulting from such coupling, we here introduce 'Active Model B'. This is a scalar ?4 field theory (or phase-field model) that minimally violates detailed balance via a leading-order square-gradient term. We find that this additional term has modest effects on coarsening dynamics, but alters the static phase diagram by creating a jump in (thermodynamic) pressure across flat interfaces. Both results are surprising, since interfacial phenomena are always strongly implicated in coarsening dynamics but are, in detailed-balance systems, irrelevant for phase equilibria.
active phase separation
time reversal symmetry
active matter
Spontaneous motility of passive emulsion droplets in polar active gels
G de Magistris
;
A Tiribocchi
;
C A Whitfield
;
R J Hawkins
;
M E Cates
;
D Marenduzzo
We study by computer simulations the dynamics of a droplet of passive, isotropic fluid, embedded in a polar active gel. The latter represents a fluid of active force dipoles, which exert either contractile or extensile stresses on their surroundings, modelling for instance a suspension of cytoskeletal filaments and molecular motors. When the polarisation of the active gel is anchored normal to the droplet at its surface, the nematic elasticity of the active gel drives the formation of a hedgehog defect; this defect then drives an active flow which propels the droplet forward. In an extensile gel, motility can occur even with tangential anchoring, which is compatible with a defect-free polarisation pattern. In this case, upon increasing activity the droplet first rotates uniformly, and then undergoes a discontinuous nonequilibrium transition into a translationally motile state, powered by bending deformations in the surrounding active medium.
Active emulsion
Liquid crystals
Lattice Boltzmann simulations
We focus on three different convexity principles for local and nonlocal variational integrals. We prove various generalizations of them, as well as their equivalences. Some applications to nonlinear eigenvalue problems and Hardy-type inequalities are given. We also prove a measure-theoretic minimum principle for nonlocal and non- linear positive eigenfunctions.
Nonlinear eigenvalue problems
Maximum Principle
Uniqueness of eigenfunctions
Hardy inequalities
nonlocal equations
Di Pietro R
;
Guarino S
;
Verde N V
;
DomingoFerrer J
In this paper, we survey emerging and established wireless ad-hoc technologies and we highlight their security/privacy features and deficiencies. We also identify open research issues and technology challenges for each surveyed technology. (c) 2014 Elsevier B.V. All rights reserved.
Pervasive mobile and low-end wireless technologies, such as radio-frequency identification (RFID), wireless sensor networks and the impending vehicular ad-hoc networks (VANETs), make the wireless scenario exciting and in full transformation. For all the above (and similar) technologies to fully unleash their potential in the industry and society, there are two pillars that cannot be overlooked: security and privacy. Both properties are especially relevant if we focus on ad-hoc wireless networks, where devices are required to cooperate - e.g. from routing to the application layer - to attain their goals.
In this article, we face the problem of ensuring reliability of a wireless sensor network which is monitoring a given set of points of interest while maximizing its lifetime (i.e., the amount of time over which the monitoring activity can be performed). The two objectives are contrasting. Indeed, the traditional approach to achieve reliability involves providing redundant coverage, which, however, drastically reduces the network lifetime. We propose an alternative strategy where sensors adapt their sensing radii in response to failures to restore feasibility only when needed. We provide Column Generation exact algorithms for both the traditional approach and our variant, as well as a heuristic procedure for the coverage restoration phase. The advantages of our approach are shown by means of computational tests on a set of instances and failure simulations.
In the k-labeled Spanning Forest Problem (kLSF), given a graph G with a label (color) assigned to each edge, and an integer positive value kmax we look for the minimum number of connected components that can be obtained by using at most kmax different labels. The problem is strictly related to the Minimum Labelling Spanning Tree Problem (MLST), since a spanning tree of the graph (i.e. a single connected component) would obviously be an optimal solution for the kLSF, if it can be obtained without violating the bound on kmax. In this work we present heuristic and exact approaches to solve this new problem.