Un grand nombre d'écoulements naturels et industriels mettent en jeu des particules (sédimentation,lit fluidisé, sprays...). Les écoulements chargés en particules sont bien décrits numériquement sous l'hypothèse des particules plus petites que toutes les échelles de l'écoulement. Cette thèse consiste à simuler numériquement une turbulence homogène et isotrope soutenue chargée en particules dont la taille est supérieure à l'échelle de Kolmogorov. Pour se faire une méthode de simulation a été développée au sein du code Thétis puis validée. L'originalité de cette méthode consiste en l'utilisation d'une approche de pénalisation associée à la viscosité dans la zone solide. Les particules sont transportées de façon lagrangienne. Les principaux résultats concernent trois simulations faisant varier le rapport de densité entre le fluide et le solide. Chaque simulation simule le mouvement de 512particules avec un diamètre 22 fois plus grand que l'échelle spatiale de Kolmogorov remplissant ainsi3% du volume total. La dispersion des particules est étudiée et montre des comportements comparables à ceux observés pour des particules ponctuelles. Un intérêt particulier est porté sur le régime collisionnel. On observe que la corrélation des vitesses avec le fluide environnant réduit le nombre de chocs frontaux par rapport au cas théorique de particules d'un gaz dense. L'effet de la prise en compte du fluide visqueux entre les particules (couche de lubrification) lors de la collision a été étudiée. L'écoulement moyen à l'échelle des particules est aussi analysé, mettant en évidence l'existence d'une couche de dissipation autour des particules. / Many applications and natural environment flows make use of particles (sedimentation, fluidized bed,sprays...). Particle laden flows are described correctly by numerical methods when the particles are smaller than all other spatial scales of the flow. This thesis involves the numerical simulation of a particle laden sustained homogeneous isotropic turbulence whose particle's size is larger than the Kolmogovov spatial scale. A numerical method has been developed and validated in the numerical code Thetis. The novelty of this method is the viscosity penalization approach. The particles are tracked by a Lagrangian way. The main results obtained are related to three simulations where the density ratio between the solid and the fluid varies. Each simulation reproduces the movement of 512particles whose diameter is 22 times the Kolmogorov spatial scale (3% volumetric solid fraction).The dispersion of particles is studied and has similar behavior than those observed with point particles simulations. The collision regime is also investigated. It is shown that he number of frontal collision is lower than its estimate for kinetic theory of gazes because there is a correlation between the particles velocity and the surrounding fluid. The modification of the collision regime when the lubrication film between particles at collision is taken into account is studied. Finally, the averaged flow around particles is analyzed and shows that there is a dissipation layer around particles.

Ce travail a été développé dans le cadre du projet européen FACETS-ITN, dans le domaine des Neurosciences Computationnelles. Son but est d'améliorer la compréhension des réseaux de neurones stochastiques de taille finie, pour des sources corrélées à caractère aléatoire et pour des matrices de connectivité biologiquement réalistes. Ce résultat est obtenu par l'analyse de la matrice de corrélation du réseau et la quantification de la capacité de codage du système en termes de son information de Fisher. Les méthodes comprennent diverses techniques mathématiques, statistiques et numériques, dont certaines ont été importés d'autres domaines scientifiques, comme la physique et la théorie de l'estimation. Ce travail étend de précédents résultats fondées sur des hypothèses simplifiées qui ne sont pas réaliste d'un point de vue biologique et qui peuvent être pertinents pour la compréhension des principes de travail liés cerveau. De plus, ce travail fournit les outils nécessaires à une analyse complète de la capacité de traitement de l'information des réseaux de neurones, qui sont toujours manquante dans la communauté scientifique.

Le travail de cette thèse aborde le sujet de l’influence des particules non-pesantes et de taille macroscopique sur les écoulements en canal dans des conditions proches du seuil de la transition laminaire-turbulent. Les suspensions sont faiblement concentrées (fraction volumique φ = 5%). Le couplage hydrodynamique existant entre la phase dispersée et la phase continue est résolu numériquement par la Force-Coupling Method, et les particules sont suivies d’une façon lagrangienne. Dans un écoulement laminaire de Couette ou de Poiseuille plan, nous montrons que les contraintes induites par la phase solide augmentent avec l’inertie, et l’influence de la concentration est plus faible qu’en régime de Stokes. Les particules avancent avec un retard dans la direction de l’écoulement et migrent à travers les lignes de courant (effet Segré-Silberberg en Poiseuille). Les vitesses de migration et de glissement s’amplifient avec l’inertie et sont du même ordre de grandeur quand Rep = O(1). Quand les particules sont lâchées librement dans un écoulement de Poiseuille plan en-deça du seuil critique de transition à la turbulence, la suspension initiale- ment homogène (φ = 5%) devient stratifiée, après un temps d’écoulement de plusieurs dizaines d’unités de temps (rapport de la hauteur du canal sur la vitesse moyenne de l’écoulement). Après une centaine d’unités de temps, nous observons le développement d’une instabilité à l’interface entre la zone chargée en particules et la zone de fluide pur. Des motifs dunaires prennent place dans la direction de la vorticité. Ces motifs sont soutenus par des écoulements secondaires d’intensités faibles mais non-nulles. Dans un écoulement au-dessus du seuil de transition, nous avons étudié les profils des phases continues et dispersées et réalisé des visualisations 3D afin de comprendre pourquoi les particules macroscopiques diminuent le nombre de Reynolds critique de relaminarisation de l’écoulement. Nous observons que les particules provoquent une augmentation significative des fluctuations de vitesses dans les directions transverses et qu’elles modifient les structures rotationnelles de l’écoulement, qui deviennent plus petites, plus nombreuses et plus énergétiques (plus grandes vitesses de rotation). Le coefficient de frottement pariétal de l’écoulement de suspension en régime de transition est supérieur à celui de l’écoulement monophasique. Quand le nombre de Reynolds est diminué et que l’écoulement devient finalement laminaire, le coefficient de frottement pariétal rejoint la loi laminaire d’un écoulement monophasique, à condition de substituer la viscosité effective de la suspension à la viscosité du fluide dans l’expression du nombre de Reynolds. D’après nos résultats, la turbulence de l’écoulement de suspension est conservée jusqu’à des nombres de Reynolds bien inférieurs à celui de l’écoulement monophasique en canal, en accord avec les observations ex- périmentales de Matas, Morris et Guazzelli (PRL, 2003) pour une géométrie cylindrique. Par ailleurs, nous montrons que selon le sens de la transition, laminaire → turbulent ou turbulent → laminaire, le nombre de Reynolds critique de transition d’un régime à l’autre n’est pas le même. Nous n’avons pas observé d’influence significative de la concentration en ce qui concerne la valeur du nombre de Reynolds critique de relaminarisation pour les deux concentrations étudiées (φ = 2.5% et 5%). / This PhD addresses the influence of macroscopic and neutrally buoyant particles on the channel flows close to the laminar-turbulent transition regime. The suspension flow is moderately concentrated (solid volumetric concentration _ = 5%). The hydrodynamic coupling between the dispersed and carrier fluid is numerically resolved using the Force-Coupling Method approach. Particle trajectories are obtained by lagrangian tracking. In laminar wall-bounded flows, we show that the stress induced by the solid phase increases with inertia, and that the effect of the concentration is weaker than in the Stokes regime. The particles lag the flow and they migrate across the streamlines (Segré-Silberberg effect in Poiseuille flow). The migration and slip velocities are of the same order of magnitude for Rep = O(1). When the particles are freely suspended in a Poiseuille flow below the transition threshold, the initially homogeneous suspension (_ = 5%) becomes stratified after several ten time units (channel height/average flow velocity). After a hundred time units, the different rheological properties of the suspension segregated parts induce an instability yielding the formation of dune-like patterns, sustained by weak but finite secondary flows. In the fluctuating flow regime, we studied the profiles of the continuous and dispersed phase and realized 3D visualizations in order to understand why finite size particles delay the relaminarization threshold. The particles induce a significant increase of the velocity fluctuations in the transverse directions and they modify the rotational flow structures, which become smaller, more numerous and more energetic (larger rotation velocity). The wall-friction coefficient of the suspension flow in the transition regime is larger than the single-phase flow case. When the Reynolds number is decreased and the flow becomes laminar, the friction coefficient recovers the laminar law of a single phase flow provided that the fluid viscosity is replaced by the effective suspension viscosity in the Reynolds number definition. Our results clearly show that the two-phase channel flow turbulence is conserved down to a threshold well below the single phase flow limit, in agreement with the observations of Matas, Morris et Guazzelli (PRL, 2003) for a cylindrical geometry. In addition, we show that according to the transition direction, i.e. laminar 7! turbulent or turbulent 7! laminar, the switch from a regime to another does not occur at the same critical Reynolds number. Finally, in the limit of moderately concentrated (_ = 2.5−5%) suspension flow in a channel, the concentration has no significant influence on the critical Reynolds number.

Dans cette thèse, nous étudions la robustesse dans le contexte de la modélisation de systèmes complexes par les automates cellulaires. En effet, si l'on cherche à reproduire un comportement émergent à partir d'un modèle d'automate cellulaire, il nous semble nécessaire de se demander si les comportements observés sont bien le résultat d'interactions entre entités constituantes, ou bien s'ils dépendent d'une définition particulière du modèle. Nous allons ainsi être amenés à considérer la robustesse du modèle, à savoir la résistance de son comportement à de petites variations sur les attributs de sa définition. Dans un premier temps, nous montrons la pertinence de cette approche en considérant plusieurs définitions possibles d'une perturbation de la mise à jour globale et en les appliquant à une classe simple et représentative de modèles d'automates cellulaires, les Automates Cellulaires Elémentaires. Nous observons que, malgré le fait que nos perturbations soient proches et qu'une majorité des modèles considérés ne change pas de comportement, quelques cas particuliers montrent des changements qualitatifs du comportement que nous étudions plus en détail. Dans un second temps, nous appliquons cette approche en nous penchant sur un modèle particulier d'automate cellulaire, qui simule le phénomène de formation d'essaim à partir d'un modèle évolué d'automate cellulaire, le gaz sur réseau. Nous explorons la robustesse du comportement du modèle en considérant la perturbation de deux attributs du modèle, la forme de la grille cellulaire et la mise à jour globale, et en tirons les conclusions sur la relation entre l'observation du comportement et la définition précise du modèle / In this thesis, we study the role of robustness in the context of the modelling of complex systems by cellular automata. Indeed, if we consider a cellular automaton which aims at reproducing an emergent behaviour from a similar structure, we want to determine whether its observed dynamics are the result of the interaction of entities, or whether it depends a precise definition of the model. We thus consider the model's robustness, that is, the resistance of the behaviour to small perturbations on the model features. First, we show the relevance of this approach by considering several definitions of a perturbation of the global updating and by applying them to a simple and representative class of cellular automata, the Elementary Cellular Automata. We observe that, despite the fact that most models show little or no change between the different perturbations, some particular cases show qualitative changes that we study in detail. Second, we apply this approach to a particular model of cellular automata, which simulates a swarming behaviour based on a lattice-gas model. We then explore the model robustness by considering the pertubations of two of the model's attributes, the lattice shape and the global updating, and discuss the relationship between the observation of the behaviour and the precise definitions of the model

Ce travail a été développé dans le cadre du projet européen FACETS-ITN, dans le domaine des Neurosciences Computationnelles. Son but est d’améliorer la compréhension des réseaux de neurones stochastiques de taille finie, pour des sources corrélées à caractère aléatoire et pour des matrices de connectivité biologiquement réalistes. Ce résultat est obtenu par l’analyse de la matrice de corrélation du réseau et la quantification de la capacité de codage du système en termes de son information de Fisher. Les méthodes comprennent diverses techniques mathématiques, statistiques et numériques, dont certaines ont été importés d’autres domaines scientifiques, comme la physique et la théorie de l’estimation. Ce travail étend de précédents résultats fondées sur des hypothèses simplifiées qui ne sont pas réaliste d’un point de vue biologique et qui peuvent être pertinents pour la compréhension des principes de travail liés cerveau. De plus, ce travail fournit les outils nécessaires à une analyse complète de la capacité de traitement de l’information des réseaux de neurones, qui sont toujours manquante dans la communauté scientifique. / The brain is the most complex system in the known universe. Its nested structure with small-world properties determines its function and behavior. The analysis of its structure requires sophisticated mathematical and statistical techniques. In this thesis we shed new light on neural networks, attacking the problem from different points of view, in the spirit of the Theory of Complexity and in terms of their information processing capabilities. In particular, we quantify the Fisher information of the system, which is a measure of its encoding capability. The first technique developed in this work is the mean-field theory of rate and FitzHugh-Nagumo networks without correlations in the thermodynamic limit, through both mathematical and numerical analysis. The second technique, the Mayer’s cluster expansion, is taken from the physics of plasma, and allows us to determine numerically the finite size effects of rate neurons, as well as the relationship of the Fisher information to the size of the network for independent Brownian motions. The third technique is a perturbative expansion, which allows us to determine the correlation structure of the rate network for a variety of different types of connectivity matrices and for different values of the correlation between the sources of randomness in the system. With this method we can also quantify numerically the Fisher information not only as a function of the network size, but also for different correlation structures of the system. The fourth technique is a slightly different type of perturbative expansion, with which we can study the behavior of completely generic connectivity matrices with random topologies. Moreover this method provides an analytic formula for the Fisher information, which is in qualitative agreement with the other results in this thesis. Finally, the fifth technique is purely numerical, and uses an Expectation-Maximization algorithm and Monte Carlo integration in order to evaluate the Fisher information of the FitzHugh-Nagumo network. In summary, this thesis provides an analysis of the dynamics and the correlation structure of the neural networks, confirms this through numerical simulation and makes two key counterintuitive predictions. The first is the formation of a perfect correlation between the neurons for particular values of the parameters of the system, a phenomenon that we term stochastic synchronization. The second, which is somewhat contrary to received opinion, is the explosion of the Fisher information and therefore of the encoding capability of the network for highly correlated neurons. The techniques developed in this thesis can be used also for a complete quantification of the information processing capabilities of the network in terms of information storage, transmission and modification, but this would need to be performed in the future.

Dans cette thèse, nous explorons le comportement critiques de chaînes de spins gouvernées par des interactions à décroissance algébrique. Dans une première partie, nous étudions le diagramme de phase d'une chaîne de Potts en utilisant un algorithme multicanonique. Nous proposons une nouvelle méthode de détection de l'ordre des transitions de phase exploitant les points spinodaux. A l'aide de cette méthode, nous localisons la ligne séparant les transitions du premier et du second ordre avec une précision sans précédent, et mettons en évidence un effet de taille finie inhabituel. Dans une deuxième partie, nous introduisons une nouvelle méthode multicanonique intégrant un algorithme de mise-à-jour collective des spins. Cette méthode, extrêmement souple, étend considérablement l'intervalle de tailles simulables, et s'avère bien plus précise que les méthodes multicanoniques usuelles. Nous appliquons cette méthode à l'étude d'effets de taille finie dans le cadre de transitions du premier ordre : les résultats suggèrent fortement que les configurations correspondant à des phases en coexistence sont caractérisées par une dimension fractale dépendent du taux de décroissance de l'interaction. Dans un dernier chapitre, nous étudions une chaîne d'Ising régie par des interactions à longue portée en présence de champs aléatoires à distribution bimodale, et prouvons l'existence d'un point tricritique pour des interactions à décroissance lente.

Les nano-objets unidimensionnels alliés présentent des propriétés physiques spécifiques qui résultent à la fois de leur morphologie, de leur taille et de la répartition chimique des atomes. Nous exploitons un modèle d'Ising sur réseau qui rend compte en particulier des effets de ségrégation au sein de nanofils bimétalliques pour obtenir une compréhension fine des effets gouvernant cette répartition à l'équilibre.Dans une première section, nous détaillons l'équilibre d'un nanofil en fonction de sa taille et de sa composition, de manière à mettre en évidence le rôle des effets de taille finie sur la thermodynamique d'équilibre d'objets bimétalliques 1D. Contrairement aux systèmes infinis, l'équilibre dépend de l'ensemble statistique considéré. Ainsi la ségrégation est plus marquée dans l'ensemble canonique, où la concentration du nanofil est imposée, que dans l'ensemble pseudo-Grand Canonique (p-GC) où le nanofil est en équilibre avec un réservoir qui fixe la différence de potentiel chimique entre les espèces. De même, la contrainte de composition dans l'ensemble canonique induit des corrélations chimiques d'occupation des sites qui favorisent davantage les paires hétéroatomiques. Nous montrons que l'écart observé entre les isothermes des deux ensembles croît avec la courbure de l'isotherme canonique et avec l'amplitude des fluctuations de la concentration nominale dans l'ensemble p-GC. Ces fluctuations diminuant avec la taille du nanofil considéré, l'écart entre les ensembles s'annule à la limite thermodynamique. Les effets de taille finie se traduisent par ailleurs par l'apparition, à basse température et pour de petits nanofils, d'une coexistence d'un mode pur en l'espèce ségrégeante et d'un mode de faible concentration nominale constitué principalement de configurations de type cœur-coquille et Janus. Nous développons alors un formalisme permettant de caractériser cette bimodalité.Alors que les résultats évoqués précédemment concernent un nanofil considéré seul, nous étudions dans la deuxième section l'équilibre de l'ensemble des nanofils formant un co-dépôt unidimensionnel inférieur à la mono-couche. Nous montrons que la distribution en taille de ces nanofils varie globalement selon une loi de puissance, quelle que soit la composition du codépôt, de sorte que la ségrégation n'a que peu d'influence sur la microstructure observée. Par contre, en raison du rapport surface/volume et des corrélations chimiques dans ces objets, la composition des nanofils du co-dépôt varie très fortement selon leur taille, les petits nanofils étant plus riches en l'espèce ségrégeante que les plus grands. Enfin, nous étendons le diagramme de bimodalité d'un nanofil seul à l'ensemble des nanofils du co-dépôt et montrons que cette bimodalité est difficilement observable car elle ne concerne que des amas de petite taille qui sont très minoritaires du fait de la cohésion atomique.

Ce mémoire présente des travaux de thèse consacrés à la caractérisation et la modélisation des hétérogénéités spatio-temporelles dans les écoulements granulaires cisaillés entre deux parois rigides. De nombreuses simulations ont permis de révéler le rôle crucial des dimensions de l’écoulement et l’influence des interactions des grains avec les parois (avec leur rugosité) sur le comportement global. Pour des systèmes allant jusqu'à une épaisseur de cent diamètres, des profils de vitesse non homogènes ont été mis en évidence, alors que la contrainte cisaillante est constante dans le volume, mettant en défaut le modèle viscoplastique au sein de l’écoulement. Pour réconcilier ces deux observations, nous avons enrichi le modèle viscoplastique en reliant la viscosité à une variable interne qui porte la perturbation due à la présence des parois. Cette nouvelle formulation de la viscosité permet de rétablir la validité d'une loi de comportement locale prenant en compte simultanément l’épaisseur de l’écoulement, la rugosité des parois et le nombre inertiel.Cette dépendance de la viscosité à une variable interne liée à la connectivité des grains ou à leur agitation à compacité fixée, suggère que, d’une manière générale,les écoulements granulaires doivent être décrits en termes d’au moins trois paramètres en fonction du nombre inertiel : le coefficient de frottement, la compacité et la connectivité. / This PhD work is devoted to the description and modeling of spatiotemporal inhomogeneities in granular flows sheared between two rigid walls. Our extensive simulations reveal the crucial role played by flow dimensions and the interactions of the grains with the walls and their roughness. For granular systems with increasingly larger thickness, non uniform strain profiles are evidenced while the shear stress remains uniform. This observation contradicts the common viscoelastic approach based on inertial number in the bulk of the flow. In order to reconcile these observations, we propose a viscoplastic model by a introducing an internal variable carrying the wall-induced perturbations of the flow. This re-formulation of granular viscosity reaffirms the local rheology by accounting for flow thickness, wall routines and inertialeffects. The well-defined dependence of the viscosity on an internal variable pertaining to grain connectivity or kinematic randomness at constant packing fraction, suggests that granular flows should be described by at least three parameters as a function of the inertial number:friction coefficient, packing fraction and connectivity.

Dans cette thèse nous présentons des progrès algorithmiques ainsi que plusieurs applications des méthodes de Monte Carlo quantique (QMC) pour simulations à partir des premiers principes. Les améliorations que nous proposons permettent d'étudier par QMC des systèmes de plus grosse taille voire périodiques, avec l'ambition de faire du QMC une alternative valable à la théorie de la fonctionnelle de la densité (DFT). Tous les résultats ont été obtenus par le logiciel TurboRVB. D'abord, nous présentons une implémentation du QMC basée sur la fonction d'onde Jastrow-Geminale qui combine une grande flexibilité avec un traitement précis des corrélations électroniques. On a appliqué une technique originale de plongement pour réduire la taille de la base atomique à la molécule d'eau ainsi qu'à un modèle simplifié du transfert de protons (TP) dans l'eau. Nos résultats ouvrent la voie à l'étude des phénomènes microscopiques tels que le TP directement par QMC. Ensuite, on a amélioré notre méthode afin de simuler les solides cristallins. Grâce à une nouvelle procédure pour choisir de manière appropriée les conditions aux limites, nous avons pu réduire les erreurs de taille finie qui affectent les simulations QMC des solides. Sur la base des techniques développées, nous étudions enfin le supraconducteur FeSe. Le QMC fournit le meilleur résultat concernant sa structure cristalline; via une étude systématique du paysage énergétique à différentes configurations magnétiques, nous montrons un lien fort entre la structure, le magnétisme et les mouvements de charge dans ce matériau, prélude à une compréhension quantitative de la supraconductivité à haute température des premiers principes. / In this thesis we present algorithmic progresses as well as applications of continuum quantum Monte Carlo (QMC) methods for electronic structure calculations by first principles. The improvements we propose allow to tackle much larger molecular as well as extended systems by QMC, with the ultimate goal of making QMC a valid alternative to density functional theory (DFT). All results have been obtained with the TurboRVB software, which we contributed to develop. At first, we present a QMC framework based on the Jastrow-Geminal wavefunction which combines great flexibility with a compact analytical form, while providing at the same time an accurate treatment of electron correlations. We apply an original atomic embedding scheme for reducing the basis set size to the water molecule and to a simple model of proton transfer (PT) in aqueous systems. Our results pave the way to the study of microscopic phenomena such as PT directly by QMC. Afterwards, we extend our QMC framework in order to simulate crystalline solids. We propose a novel procedure to find special values of the boundary conditions which allow to greatly reduce the finite-size errors affecting solid state QMC simulations. Using the techniques previously developed, we study the iron-based superconductor FeSe. We show that QMC provides the best crystal structure predictions on this compound; by means of a systematic study of the energy landscape at different magnetic orderings, we show a strong link between structural, magnetic and charge degrees of freedom in FeSe. Our results represent an important step towards a quantitative understanding of high-temperature superconductivity by first-principles.

Cette thèse est consacrée aux écoulements homogènes de bulles, ainsi qu'à leur couplage avec le transport d'un scalaire et la turbulence. Elle s'intéresse plus spécifiquement aux effets de taille finie, des interactions hydrodynamiques et de la microstructure de la suspension qui sont étudiés à l'aide de simulations numériques directes à l'échelle d'une seule bulle. La dynamique d'une suspension laminaire de bulles induite par la seule gravité est d'abord revisitée. L'influence de la fraction volumique sur la vitesse de dérive des bulles est établie analytiquement et numériquement pour une suspension parfaitement ordonnée, puis des ressemblances entre suspensions ordonnées et suspensions désordonnées sont mises en évidence. Ces résultats sont ensuite mis à profit pour la modélisation du transport d'un scalaire passif au sein d'une suspension laminaire, tel que décrit par une diffusivité effective tensorielle, et des différences essentielles entre systèmes ordonnés et systèmes désordonnés concernant le transport de scalaire sont mises en exergue. Enfin, la turbulence est prise en compte dans les simulations et son interaction avec une bulle de taille finie est caractérisée. Il est montré que le comportement dynamique d'une bulle de taille comparable à la microéchelle de Taylor ressemble qualitativement à celui d'une microbulle, avec, notamment, une préférence pour certaines régions caractéristiques de l'écoulement. Une définition de l'écoulement vu par la bulle compatible avec les modèles standards de masse ajoutée et de portance est finalement proposée / This thesis is devoted to the study of homogeneous bubbly flows and their coupling with scalar transport and turbulence. It focuses on the effects of finite size, hydrodynamic interactions, and suspension microstructure, which are investigated using direct numerical simulations at the bubble scale. The dynamics of laminar buoyancy-driven bubbly suspensions is first revisited. More specifically, the effect of volume fraction on the bubble drift velocity is clarified by connecting numerical results to theory for dilute ordered systems, and similarities between perfectly ordered and free disordered suspensions are evidenced. These results are then used for the modeling of passive scalar transport in laminar suspensions as described by an effective diffusivity tensor, and crucial differences between ordered and disordered systems with respect to scalar transport are highlighted. Lastly, turbulence is included in the simulations, and its interaction with a finite-size bubble is characterized. The behavior of a bubble as large as Taylor microscale is shown to share a number of common features with that of a microbubble, most notably, the flow sampled by the bubble is biased. A definition of the liquid flow seen by the bubble, as it enters in usual models for the added mass and the lift forces, is finally proposed