In order to understand the flow profiles of complex fluids, a crucial issue concerns the emergence of spatial correlations among plastic rearrangements exhibiting cooperativity flow behaviour at the macroscopic level. In this paper, the rate of plastic events in a Poiseuille flow is experimentally measured on a confined foam in a Hele-Shaw geometry. The correlation with independently measured velocity profiles is quantified by looking at the relationship between the localisation length of the velocity profiles and the localisation length of the spatial distribution of plastic events. To complement the cooperativity mechanisms studied in foam with those of other soft glassy systems, we compare the experiments with simulations of dense emulsions based on the lattice Boltzmann method, which are performed both with and without wall friction. Finally, unprecedented results on the distribution of the orientation of plastic events show that there is a non-trivial correlation with the underlying local shear strain. These features, not previously reported for a confined foam, lend further support to the idea that cooperativity mechanisms, originally invoked for concentrated emulsions (Goyon et al., Nature, vol. 454, 2008, pp. 84-87), have parallels in the behaviour of other soft glassy materials.
Background
Cell organization is governed and maintained via specific interactions among its constituent macromolecules. Comparison of the experimentally determined protein interaction networks in different model organisms has revealed little conservation of the specific edges linking ortholog proteins. Nevertheless, some topological characteristics of the graphs representing the networks - namely non-random degree distribution and high clustering coefficient - are shared by networks of distantly related organisms. Here we investigate the role of the topological features of the protein interaction network in promoting cell organization.
Methods
We have used a stochastic model, dubbed ProtNet representing a computer stylized cell to answer questions about the dynamic consequences of the topological properties of the static graphs representing protein interaction networks.
Results
By using a novel metrics of cell organization, we show that natural networks, differently from random networks, can promote cell self-organization. Furthermore the ensemble of protein complexes that forms in pseudocells, which self-organize according to the interaction rules of natural networks, are more robust to perturbations.
Conclusions
The analysis of the dynamic properties of networks with a variety of topological characteristics lead us to conclude that self organization is a consequence of the high clustering coefficient, whereas the scale free degree distribution has little influence on this property.
We consider an initial-boundary value problem for a degenerate pseudoparabolic regularization of a nonlinear forward-backward-forward parabolic equation, with a bounded nonlinearity which is increasing at infinity. We prove existence of suitably defined nonnegative solutions of the problem in a space of Radon measures. Solutions satisfy several monotonicity and regularization properties; in particular, their singular part is nonincreasing and may disappear in finite time. The problem is of intrinsic mathematical interest, but also arises naturally when studying, by time reversal, the spontaneous appearance of singularities in a specific application.
Travelling wave solutions of a parabolic-hyperbolic system for contact inhibition of cell-growth
Bertsch M
;
Hilhorst D
;
Izuhara H
;
Mimura M
;
Wakasa T
We consider a cell growth model involving a nonlinear system of partial differential equations which describes the growth of two types of cell populations with contact inhibition. Numerical experiments show that there is a parameter regime where, for a large class of initial data, the large time behaviour of the solutions is described by a segregated travelling wave solution with positive wave speed c. Here, the word segregated expresses the fact that the different types of cells are spatially segregated, and that the single densities are discontinuous at the moving interface which separates the two populations. In this paper, we show that, for each wave speed c > c, there exists an overlapping travelling wave solution, whose profile is continuous and no longer segregated. We also show that, for a large class of initial functions, the overlapping travelling wave solutions cannot represent the large time profile of the solutions of the system of partial differential equations. The structure of the travelling wave solutions strongly resembles that of the scalar Fisher-KPP equation, for which the special role played by the travelling wave solution with minimal speed has been extensively studied.
a parabolic-hyperbolic system
contact inhibition of cells
Fisher-KPP equation
travelling wave solutions
tumour growth
In [Comm. Appl. Math. Comput. Sci., 4 (2009), pp. 153-175], Barenblatt presents a model for partial laminarization and acceleration of shear flows by the presence of suspended particles of different sizes, and provides a formal asymptotic analysis of the resulting velocity equation. In the present paper we revisit the model. In particular we allow for a continuum of particle sizes, rewrite the velocity equation in a form which involves the Laplace transform of a given function or measure, and provide several rigorous asymptotic expansions for the velocity. The model contributes to a better insight to the extreme velocities in hurricanes, fire storms, and dust storms, and the analysis confirms Barenblatt's conclusion that often the smallest suspended particles are responsible for the extreme flow acceleration at large altitudes.
Entropy-Assisted Computing of Low-Dissipative Systems
Karlin Ilya V
;
Boesch Fabian
;
Chikatamarla Shyam S
;
Succi Sauro
Entropy feedback is reviewed and highlighted as the guiding principle to reach extremely low dissipation. This principle is illustrated through turbulent flow simulations using the entropic lattice Boltzmann scheme.
The Role of Very Low-Reynolds Hydrodynamics on the Transfer of Information Among Active Agents
De Rosis Alessandro
;
Ubertini Francesco
;
Ubertini Stefano
;
Succi Sauro
We investigate the role of hydrodynamic interactions on the decision-making and leader-identification processes within a group of fifty small-size active individuals, immersed in a viscous fluid at very low Reynolds number, . A fraction of the individuals is informed about the spatial location of the target, and moves accordingly along a privileged trajectory. The rest of the group has no access to this information, but may draw indirect benefit by following the trajectory of the informed individuals, through a process of leader-identification. Such process responds to simple behavioral rules ("social" interactions) discussed previously in the literature (Krause and Ruxton in Living in groups, 2002). The above scenario is enriched with two mechanical ingredients: the presence of an obstacle, preventing the informed individuals from following a straight trajectory to the target and hydrodynamic interactions with the surrounding fluid. It is found that hydrodynamic interactions are particularly effective in steering the uninformed individuals towards the target under neutral conditions, i.e. whenever the fraction of informed individuals is around 50 %. At lower fractions, only the informed individuals manage to reach the target, regardless of hydrodynamic interactions. Likewise, at higher fractions, all individuals reach the target, independently of hydrodynamic effects. This shows that, while hydrodynamics is subdominant under most circumstances, it may nonetheless take on a strategic role whenever the informed and uninformed individuals become comparable in number.
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