The Dynamics of Force Fluctuations in Shear Granular Flows

Tämän työn tavoitteena oli tutkia rakeisen materiaalin kinematiikkaa ja rakentaa koelaitteisto rakeisen materiaalin leikkausjännitysvirtauksien tutkimiseen. Kokeellisessa osassa on keskitytty sisäisiin voimaheilahteluihin ja niiden ymmärtämiseen. Teoriaosassa on käyty läpi rakeisen materiaalin yleisiä ominaisuuksia ja lisäksi on esitetty kaksi eri tapaa mallintaa fysikaalisien ominaisuuksien heilahteluja rakeisessa materiaalissa. Nämä kaksi esitettyä mallinnusmenetelmää ovat skalaarinen q-malli ja simulointi. Skalaarinen q-malli määrittelee jokaiseen yksittäiseen rakeeseen kohdistuvan jännityksen, rakeen ollessa osa 2- tai 3-dimensionaalista asetelmaa. Tämän mallin perusidea on kuvata jännityksien epähomogeenisuutta, joka johtuu rakeiden satunnaisasettelusta. Simulointimallinnus perustuu event-driven algoritmiin, missä systeemin dynamiikkaa kuvataan yksittäisillä partikkelien törmäyksillä. Törmäyksien vaiheet ratkaistiin käyttämällä liikemääräyhtälöitä ja restituution määritelmää. Teoriaosuudessa käytiin vielä pieniltä osin läpi syitä jännitysheilahteluihin ja rakeisen materiaalin lukkiintumiseen. Tutkimuslaitteistolla tutkittiin rakeisen materiaalin käyttäytymistä rengasmaisessa leikkausjännitysvirtauksessa. Tutkimusosuuden päätavoitteena oli mitata partikkelien kosketuksista ja törmäyksistä johtuvia hetkellisiä voimaheilahteluja rengastilavuuden pohjalta. Rakeisena materiaalina tutkimuksessa käytettiin teräskuulia. Jännityssignaali ajan funktiona osoittaa suurta heilahtelua, joka voi olla jopa kertalukua keskiarvosta suurempaa. Tällainen suuren amplitudin omaava heilahtelu on merkittävä haittapuoli yleisesti rakeisissa materiaaleissa käytettyjen jatkuvuusmallien kanssa. Tällainen heilahtelu tekee käytetyt jatkuvuusmallit epäpäteviksi. Yleisellä tasolla jännityksien todennäköisyysjakauma on yhtäpitävä skalaarisen q-mallin tuloksien kanssa. Molemmissa tapauksissa todennäköisyysjakaumalla on eksponentiaalinen muoto.

Impact of open-ocean convection on particle fluxes and sediment dynamics in the deep margin of the Gulf of Lions.

Abstract. The deep outer margin of the Gulf of Lions and the adjacent basin, in the western Mediterranean Sea, are regularly impacted by open-ocean convection, a major hydrodynamic event responsible for the ventilation of the deep water in the western Mediterranean Basin. However, the impact of open-ocean convection on the flux and transport of particulate matter remains poorly understood. The variability of water mass properties (i.e., temperature and salinity), currents, and particle fluxes were monitored between September 2007 and April 2009 at five instrumented mooring lines deployed between 2050 and 2350-m depth in the deepest continental margin and adjacent basin. Four of the lines followed a NW-SE transect, while the fifth one was located on a sediment wave field to the west. The results of the main, central line SC2350 ("LION") located at 42 02.50 N, 4 410 E, at 2350-m depth, show that open-ocean convection reached midwater depth ( 1000-m depth) during winter 2007-2008, and reached the seabed ( 2350-m depth) during winter 2008-2009. Horizontal currents were unusually strong with speeds up to 39 cm s−1 during winter 2008-2009. The measurements at all 5 different locations indicate that mid-depth and near-bottom currents and particle fluxes gave relatively consistent values of similar magnitude across the study area except during winter 2008-2009, when near-bottom fluxes abruptly increased by one to two orders of magnitude. Particulate organic carbon contents, which generally vary between 3 and 5 %, were abnormally low ( 1 %) during winter 2008-2009 and approached those observed in surface sediments (0.6 %). Turbidity profiles made in the region demonstrated the existence of a bottom nepheloid layer, several hundred meters thick, and related to the resuspension of bottom sediments. These observations support the view that open-ocean deep convection events in the Gulf of Lions can cause significant remobilization of sediments in the deep outer margin and the basin, with a subsequent alteration of the seabed likely impacting the functioning of the deep-sea ecosystem.

Analysis of flow dynamics in right ventricular outflow tract.

BACKGROUND: The mechanism behind early graft failure after right ventricular outflow tract (RVOT) reconstruction is not fully understood. Our aim was to establish a three-dimensional computational fluid dynamics (CFD) model of RVOT to investigate the hemodynamic conditions that may trigger the development of intimal hyperplasia and arteriosclerosis. METHODS: Pressure, flow, and diameter at the RVOT, pulmonary artery (PA), bifurcation of the PA, and left and right PAs were measured in 10 normal pigs with a mean weight of 24.8 ± 0.78 kg. Data obtained from the experimental scenario were used for CFD simulation of pressure, flow, and shear stress profile from the RVOT to the left and right PAs. RESULTS: Using experimental data, a CFD model was obtained for 2.0 and 2.5-L/min pulsatile inflow profiles. In both velocity profiles, time and space averaged in the low-shear stress profile range from 0-6.0 Pa at the pulmonary trunk, its bifurcation, and at the openings of both PAs. These low-shear stress areas were accompanied to high-pressure regions 14.0-20.0 mm Hg (1866.2-2666 Pa). Flow analysis revealed a turbulent flow at the PA bifurcation and ostia of both PAs. CONCLUSIONS: Identified local low-shear stress, high pressure, and turbulent flow correspond to a well-defined trigger pattern for the development of intimal hyperplasia and arteriosclerosis. As such, this real-time three-dimensional CFD model may in the future serve as a tool for the planning of RVOT reconstruction, its analysis, and prediction of outcome.

From caging to rouse dynamics in polymer melts with intramolecular barriers: a critical test of the mode coupling theory

By means of computer simulations and solution of the equations of the mode coupling theory (MCT),we investigate the role of the intramolecular barriers on several dynamic aspects of nonentangled polymers. The investigated dynamic range extends from the caging regime characteristic of glass-formers to the relaxation of the chain Rouse modes. We review our recent work on this question,provide new results, and critically discuss the limitations of the theory. Solutions of the MCT for the structural relaxation reproduce qualitative trends of simulations for weak and moderate barriers. However, a progressive discrepancy is revealed as the limit of stiff chains is approached. This dis-agreement does not seem related with dynamic heterogeneities, which indeed are not enhanced by increasing barrier strength. It is not connected either with the breakdown of the convolution approximation for three-point static correlations, which retains its validity for stiff chains. These findings suggest the need of an improvement of the MCT equations for polymer melts. Concerning the relaxation of the chain degrees of freedom, MCT provides a microscopic basis for time scales from chain reorientation down to the caging regime. It rationalizes, from first principles, the observed deviations from the Rouse model on increasing the barrier strength. These include anomalous scaling of relaxation times, long-time plateaux, and nonmonotonous wavelength dependence of the mode correlators.

Cooperative and Competitive Spreading Dynamics on the Human Connectome.

Increasingly detailed data on the network topology of neural circuits create a need for theoretical principles that explain how these networks shape neural communication. Here we use a model of cascade spreading to reveal architectural features of human brain networks that facilitate spreading. Using an anatomical brain network derived from high-resolution diffusion spectrum imaging (DSI), we investigate scenarios where perturbations initiated at seed nodes result in global cascades that interact either cooperatively or competitively. We find that hub regions and a backbone of pathways facilitate early spreading, while the shortest path structure of the connectome enables cooperative effects, accelerating the spread of cascades. Finally, competing cascades become integrated by converging on polysensory associative areas. These findings show that the organizational principles of brain networks shape global communication and facilitate integrative function.

Analysis of Crankshaft Using Multibody Dynamics

Työn tavoitteena oli löytää tarkka menetelmä kampiakselin vaurioitumisriskin laskentaan vertailemalla eri laskentamenetelmiä. Lopuksi suoritettiin simulointi monikappalejärjestelmälle käyttäen elastisia malleja todellisista rakenteista. Simulointiohjelmana käytettiin AVL:n kehittämää Excite:ia.

Performance Analysis of Two Quantum Reaction Dynamics Codes: Time-Dependent and Time-Independent Strategies

The computer simulation of reaction dynamics has nowadays reached a remarkable degree of accuracy. Triatomic elementary reactions are rigorously studied with great detail on a straightforward basis using a considerable variety of Quantum Dynamics computational tools available to the scientific community. In our contribution we compare the performance of two quantum scattering codes in the computation of reaction cross sections of a triatomic benchmark reaction such as the gas phase reaction Ne + H2+ %12. NeH++ H. The computational codes are selected as representative of time-dependent (Real Wave Packet [ ]) and time-independent (ABC [ ]) methodologies. The main conclusion to be drawn from our study is that both strategies are, to a great extent, not competing but rather complementary. While time-dependent calculations advantages with respect to the energy range that can be covered in a single simulation, time-independent approaches offer much more detailed information from each single energy calculation. Further details such as the calculation of reactivity at very low collision energies or the computational effort related to account for the Coriolis couplings are analyzed in this paper.

Self-accelerating dolomite-for-calcite replacement: Self-organized dynamics of burial dolomitization and associated mineralization

A new dynamic model of dolomitization predicts a multitude of textural, paragenetic, geochemical and other properties of burial dolomites. The model is based on two postulates, (1) that the dolomitizing brine is Mg-rich but under saturated with both calcite and dolomite, and (2) that the dolomite-for-calcite replacement happens not by dissolution-precipitation as usually assumed, but by dolomite-growth-driven pressure solution of the calcite host. Crucially, the dolomite-for-calcite replacement turns out to be self-accelerating via Ca2 : the Ca2 released by each replacement increment accelerates the rate of the next, and so on. As a result, both pore-fluid Ca2 and replacement rate grow exponentially.

Multiple processes regulate long-term population dynamics of sea urchins on Mediterranean rocky reefs

We annually monitored the abundance and size structure of herbivorous sea urchin populations (Paracentrotus lividus and Arbacia lixula) inside and outside a marine reserve in the Northwestern Mediterranean on two distinct habitats (boulders and vertical walls) over a period of 20 years, with the aim of analyzing changes at different temporal scales in relation to biotic and abiotic drivers. P. lividus exhibited significant variability in density over time on boulder bottoms but not on vertical walls, and temporal trends were not significantly different between the protection levels. Differences in densities were caused primarily by variance in recruitment, which was less pronounced inside the MPA and was correlated with adult density, indicating density-dependent recruitment under high predation pressure, as well as some positive feedback mechanisms that may facilitate higher urchin abundances despite higher predator abundance. Populations within the reserve were less variable in abundance and did not exhibit the hyper-abundances observed outside the reserve, suggesting that predation effects maybe more subtle than simply lowering the numbers of urchins in reserves. A. lixula densities were an order of magnitude lower than P. lividus densities and varied within sites and over time on boulder bottoms but did not differ between protection levels. In December 2008, an exceptionally violent storm reduced sea urchin densities drastically (by 50% to 80%) on boulder substrates, resulting in the lowest values observed over the entire study period, which remained at that level for at least two years (up to the present). Our results also showed great variability in the biological and physical processes acting at different temporal scales. This study highlights the need for appropriate temporal scales for studies to fully understand ecosystem functioning, the concepts of which are fundamental to successful conservation and management.

Recombination dynamics in ZnO nanowires: Surfaces states versus mode quality factor

In this work, we investigate the influence of finite size on the recombinations dynamics of ZnO nanowires. We demonstrate that diameter as well as lenght of nanowires determine the lifetime of the neutral donor bound excitons. Our findings suggest that while the length is mainly responsible for different mode quality factors of the cavity-like nanowires, the diameter determines the influence of surface states as alternative recombinations channels for the optical modes trapped in the nanocavity. In addition, comparing nanowires grown using different catalyst we show that the surfaces states strongly depend on each precursor characteristics.

The dynamics of coexistence: habitat sharing versus segregation patterns among three sympatric montane vipers

Contact zones of closely related and ecologically similar species constitute rare opportunities to study the evolutionary consequences of past speciation processes. They represent natural laboratories in which strong competition could lead to the exclusion of one species, or the various species may switch into distinct ecological niches. Alternatively, if reproductive isolation has not yet been achieved, they may hybridize. We elucidate the degree of taxon integrity by comparing genetics and habitat use of three similar-sized congeneric viper species, Vipera ammodytes, Viperaaspis, and Viperaberus, of Nadiza Valley in western Slovenia. No hybridization was detected for either mitochondrial or nuclear genomes. Similarly, external intermediacy by a single prestudy viper (probably V.ammodytesxV. aspis) indicates that hybridization occasionally occurs, but should be very rare. Populations of the three related viperids are partially allopatric in Nadiza Valley, but they also coexist in a narrow contact zone in the montane grassland along the south-exposed slope of Mount Stol (1673m a.s.l.). Here, the three species that occupy areas in or near patches of rocky microhabitats (e.g. stone piles, slides, and walls) live in syntopy. However, fine-scale measurements of structural components show partial habitat segregation, in which V.berus becomes more dominant at elevations above 1400m and occupies mostly the mountain ridge and north-exposed slopes of Mount Stol, V.aspis occurs below 1300m and is the only species to inhabit stoneless patches of grass and bushes around 1000m and lower, and V.ammodytes occurs at all elevations up to 1500m, but is restricted to a rocky microhabitat. We suggest that a high degree of microstructure divergence, slightly different environmental niches, and a generally favourable habitat for all three viper species, keep the pressure for mis-mating and hybridization low, although mechanisms such as reduced hybrid inferiority and temporal mating segregation cannot yet be excluded.

Evolutionary dynamics of collective action in spatially structured populations.

Many models proposed to study the evolution of collective action rely on a formalism that represents social interactions as n-player games between individuals adopting discrete actions such as cooperate and defect. Despite the importance of spatial structure in biological collective action, the analysis of n-player games games in spatially structured populations has so far proved elusive. We address this problem by considering mixed strategies and by integrating discrete-action n-player games into the direct fitness approach of social evolution theory. This allows to conveniently identify convergence stable strategies and to capture the effect of population structure by a single structure coefficient, namely, the pairwise (scaled) relatedness among interacting individuals. As an application, we use our mathematical framework to investigate collective action problems associated with the provision of three different kinds of collective goods, paradigmatic of a vast array of helping traits in nature: "public goods" (both providers and shirkers can use the good, e.g., alarm calls), "club goods" (only providers can use the good, e.g., participation in collective hunting), and "charity goods" (only shirkers can use the good, e.g., altruistic sacrifice). We show that relatedness promotes the evolution of collective action in different ways depending on the kind of collective good and its economies of scale. Our findings highlight the importance of explicitly accounting for relatedness, the kind of collective good, and the economies of scale in theoretical and empirical studies of the evolution of collective action.

Fluid dynamics simulation of right ventricular outflow tract oversizing.

OBJECTIVES: Repair of the right ventricular outflow tract (RVOT) in paediatric cardiac surgery remains challenging due to the high reoperation rate. Intimal hyperplasia and consequent arteriosclerosis is one of the most important limitation factors for graft durability. Since local shear stress and pressure are predictive elements for intimal hyperplasia and wall degeneration, we sought to determine in an oversized 12-mm RVOT model, with computed fluid dynamics simulation, the local haemodynamical factors that may explain intimal hyperplasia. This was done with the aim of identifying the optimal degree of oversizing for a 12-mm native RVOT. METHODS: Twenty domestic pigs, with a weight of 24.6 ± 0.89 kg and a native RVOT diameter of 12 ± 1.7 mm, had valve conduits of 12, 16, 18 and 20 mm implanted. Pressure and flow were measured at 75, 100 and 125% of normal flow at RVOT at the pulmonary artery, pulmonary artery bifurcation and at the left and right pulmonary arteries. Three-dimensional computed fluid dynamics (CFD) simulation in all four geometries in all flow modalities was performed. Local shear stress and pressure conditions were investigated. RESULTS: Corresponding to 75, 100 and 125% of steady-state flow, three inlet velocity profiles were obtained, 0.2, 0.29 and 0.36 m/s, respectively. At inflow velocity profiles, low shear stress areas, ranged from 0 to 2 Pa, combined with high-pressure areas ranging from 11.5 to 12.1 mmHg that were found at distal anastomosis, at bifurcation and at the ostia of the left and right pulmonary arteries in all geometries. CONCLUSIONS: In all three oversized geometries, the local reparation of shear stress and pressure in the 16-mm model showed a similar local profile as in the native 12 mm RVOT. According to these findings, we suggest oversizing the natural 12-mm RVOT by not more than 4 mm. The elements responsible for wall degeneration and intimal hyperplasia remain very similar to the conditions present in native RVOT.

Nature of the Epidemic Threshold for the Susceptible-Infected-Susceptible Dynamics in Networks

We develop an analytical approach to the susceptible-infected-susceptible epidemic model that allows us to unravel the true origin of the absence of an epidemic threshold in heterogeneous networks. We find that a delicate balance between the number of high degree nodes in the network and the topological distance between them dictates the existence or absence of such a threshold. In particular, small-world random networks with a degree distribution decaying slower than an exponential have a vanishing epidemic threshold in the thermodynamic limit.