La représentation des phénomènes météorologiques est un enjeu essentiel en Informatique Graphique dans l’optique d’obtenir des scènes extérieures visuellement réalistes. Depuis quelques années, les chercheurs se sont penchés sur la simulation de la pluie en synthèse d’images. Ce mémoire a pour but de présenter deux méthodes de rendu : une première méthode pour le rendu de la pluie, et une seconde pour le rendu des égouttements provenant des arbres. Notre modèle de rendu de pluie est une méthode multi-échelle permettant de tenir compte des propriétés locales et globales de ce phénomène. Nous avons réussi à mettre en corrélation la densité des gouttes proches de l’observateur (échelle mésoscopique) avec l’atténuation de la visibilité (échelle macroscopique), en fonction d’un seul paramètre global cohérent : l’intensité des précipitations. Cette méthode a fait l’objet d’une publication [56] dans la revue Computers & Graphics en 2015. D’autre part, nous basons notre méthode pour le rendu des égouttements sur une approche phénoménologique. Une telle approche s’avère plus adaptée pour gérer les grands espaces, et présente l’avantage d’être indépendante de la complexité de la scène. L’approche choisie s’appuie sur des expérimentations et des mesures hydrologiques effectuées par des chercheurs spécialisés dans l’hydrologie des forêts. Nous proposons ainsi un modèle d’égouttement cohérent, tenant compte des propriétés intrinsèques à chaque type d’arbres. Cette méthode a fait elle aussi l’objet d’une présentation à la conférence Eurographics Symposium on Rendering (EGSR) ainsi qu’une publication [57] au journal Computer Graphics Forum (CGF) en 2016. / This dissertation aims to present a coherent multiscale model for real-time rain rendering which takes into account local and global properties of rainy scenes. Our goal is to simulate visible rain streaks close to the camera as well as the progressive loss of visibility induced by atmospheric phenomena. Our model proposes to correlate the attenuation of visibility, which is due in part to the extinction phenomenon, and the distribution of raindrops in terms of rainfall intensity and camera's parameters. Furthermore, this method proposes an original rain streaks generation based on spectral analysis and sparse convolution theory. This allows an accurate control of rainfall intensity and streaks appearance, improving the global realism of rainy scenes. We also aim at rendering interactive visual effects inherent to complex interactions between trees and rain in real-time in order to increase the realism of natural rainy scenes. Such a complex phenomenon involves a great number of physical processes influenced by various interlinked factors and its rendering represents a thorough challenge in Computer Graphics. We approach this problem by introducing an original method to render drops dripping from leaves after interception of raindrops by foliage. Our method introduces a new hydrological model representing interactions between rain and foliage through a phenomenological approach. Our model reduces the complexity of the phenomenon by representing multiple dripping drops with a new fully functional form evaluated per-pixel on-the-fly and providing improved control over density and physical properties. Furthermore, an efficient real-time rendering scheme, taking full advantage of latest GPU hardware capabilities, allows the rendering of a large number of dripping drops even for complex scenes.

Les renforts à base de fibres de lin sont aujourd’hui une alternative capable de concurrencer les fibres synthétiques conventionnelles puisqu’ils sont écologiques, économiques et présentent des propriétés mécaniques intéressantes. Cependant, leur inconvénient majeur est leur absorption d’eau potentiellement importante qui affecte leurs propriétés mécaniques. Ce projet de recherche propose d’étudier le comportement hygroscopique et le couplage hygro-mécanique dans les composites lin / époxy. Cette étude repose sur une approche expérimentale multi-échelle et une modélisation du comportement visco-élasto-plastique avec prise en compte du couplage hygro-mécanique des composites renforcés par des fibres de lin. Les cinétiques de diffusion dans l’époxy et dans le composite ont été modélisées par une loi de type Langmuir et Fick respectivement. Les coefficients d’hygro-expansion des composites et des fils élémentaires qui constituent le renfort tissu ont été déterminés expérimentalement. Une étude de l’influence du conditionnement jusqu’à saturation à différentes humidités relatives sur le comportement mécanique dans les trois directions du stratifié a également été menée. Cette étude a montré l’existence d’une teneur en eau optimale pour laquelle les propriétés mécaniques sont optimales. L’émergence d’un comportement à deux régions linéaires a été mise en évidence et attribuée à la présence d’hétérogénéités locales au sein du renfort tissu. Des essais de fluage / recouvrance et de relaxation / effacement ont permis de mettre en place un modèle visco-élasto-plastique avec prise en compte du couplage hygro-mécanique. Ce modèle offre de bonnes capacités de prédiction et permettra de prévoir le comportement des structures composites renforcés par des fibres de lin en atmosphère humide. / Flax fibre reinforcements are nowadays an alternative able to compete with conventional synthetic fibres since they are ecological, economic and have interesting mechanical properties. However, their major drawback is their potentially significant water absorption which affects their mechanical properties. This research project proposes to study the hygroscopic behaviour and hygro-mechanical coupling in flax / epoxy composites. This study is based on a multi-scale experimental approach. A modelling of visco-elasto-plastic behaviour taking into account the hygro-mechanical coupling within flax /epoxy composites is established. The diffusion kinetics in composites were modelled by a Fick law. However, the diffusion kinetics in epoxy were modelled by a Langmuir law. The hygro-expansion coefficients of the composites and the elementary yarns that constitute the fabric reinforcement were determined experimentally. A study of the influence of conditioning until saturation at different relative humidities on the mechanical behaviour in the three main directions of the laminates was conducted. This study showed the existence of an optimal water content for which the mechanical properties are maximum. The emergence of a two-linear-region behaviour was pointed out and attributed to the presence of local heterogeneities within the fabric reinforcement. Creep / recovery and stress relaxation tests were exploited in order to develop a visco-elastoplastic model with consideration of the hygro-mechanical coupling. This model offers good predictive capabilities and could be used to predict the behaviour of flax fibres reinforced composite structures in humid atmospheres.

Ein metrischer Graph mit einem selbstadjungierten, negativen Laplace-Operator wird Quantengraph genannt. In dieser Arbeit werden Transporteigenschaften zufälliger Laplace-Operatoren betrachtet. Dazu wird die Multiskalenanalyse (MSA) von euklidischen Räumen auf metrische Graphen angepasst. Eine Überdeckung der metrischen Graphen wird aus gleichmäßig polynomiellem Wachstum und der gleichmäßigen Beschränkung der Kantenlängen gewonnen. Als Hilfsmittel für die MSA werden eine Combes-Thomas-Abschätzung und eine Geometrische Resolventenungleichung bewiesen. Zusammen mit einer Wegner-Abschätzung und der Existenz von verallgemeinerten Eigenfunktionen wird mittels der modifizierten MSA spektrale Lokalisierung (d.h. reines Punktspektrum) mit polynomiell fallenden Eigenfunktionen am unteren Rand des Spektrums für negative Laplace-Operatoren mit zufälligem Potential geschlossen. Dabei sind alle Randbedingungen, die eine nach unten beschränkten Operator liefern, wählbar. / We prove spectral localization for infinite metric graphs with a self-adjoint Laplace operator and a random potential. Therefor we adapt the multiscale analysis (MSA) from the euclidean case to metric graphs. In the MSA a covering of the graph is needed which is obtained from a uniform polynomial growth of the graph. The geometric restrictions of the graph contain a uniform bound on the edge lengths. As boundary conditions we allow all settings which give a lower bounded self-adjoint operator with an associated quadratic form. The result is spectral localization (i.e. pure point spectrum) with polynomially decaying eigenfunctions in a small interval at the ground state energy.

<p>Intra-Arterial Chemotherapy (IAC) is a preferred treatment for the primary liver cancer, despite its adverse side-effects. During IAC, a mixture of chemotherapeutic drugs, e.g. Doxorubicin, is injected into an artery supplying the tumor. A fraction of Doxorubicin is absorbed by the tumor, but the remaining drug passes into systemic circulation, causing irreversible heart failure. The efficiency and safety of the IAC can be improved by chemical filtration of the excessive drugs with a catheter-based Chemofilter device, as proposed by a team of neuroradilogists. </p> <p>The objective of my work was to optimize the hemodynamic and drug binding performance of the Chemofilter device, using Computational Fluid Dynamics (CFD) modeling. For this, I investigated the performance of two distinct Chemofilter configurations: 1) a porous “Chemofilter basket” formed by a lattice of micro-cells and 2) a non-porous “honeycomb Chemofilter” consisting of parallel hexagonal channels. A multiscale modeling approach was developed to resolve the flow through a representative section of the porous membrane and subsequently characterize the overall performance of the device. A heat and mass transfer analogy was utilized to facilitate the comparison of alternative honeycomb configurations. </p> A multiphysics approach was developed for modeling the electrochemical binding of Doxorubicin to the anionic surface of the Chemofilter. An effective diffusion coefficient was derived based on dilute and concentrated solution theory, to account for the induced migration of ions. Computational predictions were supported by results of <i>in-vivo</i> studies performed by collaborators. CFD models showed that the honeycomb Chemofilter is the most advantageous configuration with 66.8% drug elimination and 2.9 mm-Hg pressure drop across the device. Another facet of the Chemofilter project was its surface design with shark-skin inspired texturing, which improves the binding performance by up to 3.5%. Computational modeling enables optimization of the chemofiltration device, thus allowing the increase of drug dose while reducing systemic toxicity of IAC.

L’augmentation de la puissance informatique disponible permet aux ingénieurs et designers d’aborder par simulation des problèmes de plus en plus complexes (multi-physiques, multi-échelles, géométries intriquées ...). Dans ce contexte, les quadratures à base de discrétisation (FDM, FEM, FVM) montrent leur limite : le besoin d’un grand nombre de sous-domaines qui implique des coûts RAM et CPU prohibitifs. La méthode de Monte-Carlo apparaît plus appropriée, mais son utilisation est verrouillée par la difficulté de générer des modèles probabilistes de systèmes complexes. Pour surpasser ceci, une approche systémique est proposée et implémentée pour créer une preuve de concept appliquée à la simulation des transferts thermiques couplés. Après une étape de validation vis-à-vis de solutions analytiques, l’outil est employé; sur des cas d’illustration (transferts thermiques au sein de bâtiments et dans une centrale solaire) pour étudier ses capacités. L’approche mise en œuvre présente un comportement particulièrement avantageux pour la simulation de systèmes complexes : son temps de calcul ne dépend que des parties influentes du problème. De plus, elles sont automatiquement identifiées, même en présence de géométries étendues ou intriquées, ce qui rend les simulations auto-adaptatives. Par ailleurs, ses performances de calcul ne sont pas corrélées avec le rapport d’échelle caractérisant le système simulé. Ceci en fait une approche douée d’une remarquable capacité à traiter les problèmes à la fois multi-physiques et multi-échelles. En parallèle de l’estimation d’une observable par des chemins d’exploration, l’outil analyse également ces derniers de manière statistique. Ceci lui permet de générer un modèle prédictif réduit de l’observable, procurant ainsi une capacité d’auto-apprentissage à la simulation. Son utilisation peut améliorer les processus d’optimisation et de contrôle-commande, ou simplifier les mesures par méthodes inverses. De plus, elle a aussi permis de mener une analyse par propagation d’incertitudes, affectant les conditions aux frontières, vers l’observable. Enfin, une démonstration d’optimisation, utilisant des modèles réduits générés, a été réalisée. / As computing power increases, engineers and designers tackle increasingly complex problems using simulation (multiphysics, multiscale, intricated geometries ...). In this context, discretization-based quadratures (FDM, FEM, FVM) show their limit: the need of a great number of sub-domains which induces prohibitive consumption of RAM and CPU power. The Monte Carlo method appears to be more appropriate, but the difficulty to build probabilistic models of complex systems forms a bottleneck. A systemic approach is proposed to alleviate it and is implemented to create a proof-of-concept dedicated to the coupled heat transfer simulation. After a successful validation step against analytical solutions, this tool is applied to illustrative cases (emulating heat transfer in buildings and in solar heating systems) in order to study its simulation capabilities.This approach presents a major beneficial behavior for complex systems simulation: the computation time only depends on the influential parts of the problem. These parts are automatically identified, even in intricate or extensive geometries, which makes the simulation self-adaptive. In addition, the computational performance and the system scale ratio are completely uncorrelated. Consequently, this approach shows an exceptional capacity to tackle multiphysics and multiscale problems. Each temperature is estimated using exploration paths. By statistically analyzing these paths during the process, the tool is able to generate a reduced predictive model of this physical quantity, which is bringing a self-learning capacity to the simulation. Its use can significantly improve optimization and control of processes, or simplify inverse measurements. Furthermore, based on this model, an uncertainty propagation analysis has been performed. It quantifies the effect of uncertainties affecting boundary conditions on the temperature. Finally a Particle Swarm Optimization (PSO) process, based on simulations done by the framework, is successfully carried out.

Für Festigkeitsnachweise hat sich die Methode der Finiten Elemente (FEM) als Goldstandard etabliert. Zwar wird sowohl bei der Modellbildung als auch bei der Auswertung der Ergebnisse nach wie vor eine intellektuelle Eigenleistung gefordert, die Ergebnisse selbst sind aber unter dieser Voraussetzung zuverlässig und tendenziell reproduzierbar. Dank der Leistungsfähigkeit der heutigen Arbeitsplatzrechner werden zunehmend große Produkte betrachtet – Assemblies, die aus einer Vielzahl unterschiedlichster Parts bzw. Baugruppen bestehen. Hier begegnen wir einem neuen Phänomen. Es gibt oft Basisstrukturen, in denen Detailstrukturen enthalten sind, deren geometrische Abmessungen sich um mehrere Größenordnungen von den Gesamtabmessungen unterscheiden können. Eine gemeinsame Elementierung erweist sich dabei als wenig sinnvoll. Ebenso findet man oft eine große Anzahl von Gleichteilen, für die im Prinzip jeweils eine Mustervernetzung genügt. Selbst wenn die FE-Software dies zulassen sollte, bleibt das Problem der extrem unterschiedlichen Elementgrößen innerhalb eines Modells. Das häufig propagierte defeaturing, für das sogar Automatisierungsansätze existieren, ist ebenso wenig zielführend, weil es auf geometrische Details bezogen ist, die nicht notwendig physikalische Funktionselemente darstellen. Gerade die physikalischen Eigenschaften der Parts sollten ja erhalten bleiben und in die Analyse einfließen. In Einzelfällen werden physikalisch motivierte Vereinfachungen praktiziert; so werden Wellen auf Balkenstrukturen reduziert, wenn man sich nur für die mechanischen Eigenschaften von Rotoren interessiert. Eine Verallgemeinerung und Systematisierung solcher Individualansätze auf größere Klassen von Strukturkomponenten ist bisher nicht untersucht worden.

Le travail présenté dans ce mémoire s’intéresse à la prédiction des vibrations d'un noyau de transformateur multicouche, constitué d'un assemblage de tôles ferromagnétiques. Le problème couplé magnéto-mécanique est résolu par une approche séquentielle progressive : la résolution magnétique est suivie d'une résolution mécanique. Un modèle multi-échelle simplifié 3D décrivant les anisotropies magnétiques et magnétostrictives, et considérant les non-linéarités magnétiques et de magnétostriction, est utilisé comme loi de comportement du matériau. La structure du noyau du transformateur est modélisée en 2D. Une technique d'homogénéisation permet de tenir compte du comportement anisotrope de chaque couche afin de définir un comportement moyen pour chaque élément du maillage éléments finis.. Des mesures expérimentales sont ensuite effectuées, permettant d’une part la validation des lois de comportement matériau utilisées, et d’autres part des modèles de comportement structurel statique, du comportement structurel dynamique et de l'estimation du bruit. Différents matériaux et différentes géométries de prototypes de transformateurs sont considérés pour ce travail. Des optimisations structurelles sont finalement proposées grâce à des simulations numériques s’appuyant sur le modèle développé, afin de réduire les vibrations et les émissions de bruit du noyau du transformateur. / This thesis deals with the prediction of the vibration of a multi-layer transformer core made of an assembly of electrical sheets. This magneto-mechanical coupled problem is solved by a stepping finite element method sequential approach: magnetic resolution is followed by mechanical resolution. A 3D Simplified Multi-Scale Model (SMSM) describing both magnetic and magnetostrictive anisotropies is used as the constitutive law of the material. The transformer core structure is modeled in 2D and a homogenization technique is implemented to take the anisotropic behavior of each layer into consideration and define an average behavior at each element of the finite element mesh. Experimental measurements are then carried out, allowing the validation of the material constitutive law, static structural behavior, dynamic structural behavior, and the noise estimation. Different materials geometries are considered for this workStructural optimizations are finally achieved by numerical simulation for lower vibration and noise emission of the transformer cores.

L'objet de ce travail est de proposer une approche multi-échelle de la fatigue fondée sur l'énergie, et susceptible d'estimer les durées de vie associées à des chargements multidimensionnels variables. Le fondement de la démarche consiste à supposer que l'énergie dissipée à petite échelle régit le comportement à la fatigue. À chaque point matériel, est associée une distribution stochastique de points faibles qui sont susceptibles de plastifier et de contribuer à la dissipation d'énergie sans affecter des contraintes macroscopiques globales. Ceci revient à adopter le paradigme de Dang Van en fatigue polycyclique. La structure est supposée élastique (ou adaptée) à l'échelle macroscopique. De plus, on adopte à l'échelle mésoscopique un comportement élastoplastique avec une dépendance de la fonction de charge plastique non seulement de la partie déviatorique des contraintes, mais aussi de la partie hydrostatique. On considère également un écrouissage cinématique linéaire sous l'hypothèse d'une plasticité associée. Au lieu d'utiliser le nombre de cycles comme variable incrémentale, le concept d'évolution temporelle du chargement est adopté pour un suivi précis de l'historique du chargement réel. L'effet de la contrainte moyenne est pris en compte dans la fonction de charge mésoscopique ; une loi de cumul non linéaire de dommage est également considérée dans le modèle. La durée de vie à la fatigue est ensuite déterminée à l'aide d'une loi de phénoménologique fondée sur la dissipation d'énergie mésoscopique issue du cycle d'accommodation plastique. La première partie du travail a porté sur une proposition d'un modèle de fatigie à gradient de mise en oeuvre plus simple que les précédents modèles. / The aim of this work is to propose a multi-scale approach to energy-based fatigue, which can estimate lifetimes associated with variable multidimensional loading. The foundation of the approach is to assume that the energy dissipated on a small scale governs the fatigue behavior. Each material point is associated to a stochastic distribution of weak points that are likely to plasticize and contribute to the dissipation of energy without affecting global macroscopic stresses. This amounts to adopting Dang Van's paradigm of high cycle fatigue. The structure is supposed to be elastic (or adapted) on a macroscopic scale. In addition, we adopt on the mesoscopic scale an elastoplastic behavior with a dependence of the plastic load function not only of the deviatoric part of the stresses, but also of the hydrostatic part. Linear kinematic hardening is also considered under the assumption of an associated plasticity. Instead of using the number of cycles as an incremental variable, the concept of temporal evolution of the load is adopted for a precise follow-up of the history of the actual loading. The effect of mean stress is taken into account in the mesoscopic yield function; a law of nonlinear accumulation of damage is also considered in the model. Fatigue life is then determined using a phenomenological law based on mesoscopic energy dissipation from the plastic accommodative cycle. The first part of the work focused on a proposal for a fatigue model with a simpler implementation gradient than the previous models.

La friture profonde de type batch a été étudiée dans le projet FUI Fry’In (Réf. AAP17, 2014-2018) dans le but de proposer des innovations de rupture pour des friteuses batch domestiques et professionnelles. La thèse a appuyé le projet sur la maîtrise de deux effets négatifs de la friture : i) la thermo-oxydation de l’huile responsable des mauvaises odeurs et produits de dégradation ainsi que ii) la prise d’huile généralement favorisée au détriment de son égouttage. L’étude a été réalisée en combinant des mesures directes (spectroscopie et imagerie infrarouges en mode ATR, photo-ionisation, mesures DSC, imagerie rapide…) et modélisation multi-échelle (écoulement de l’huile et égouttage lors du retrait, description lagrangienne des réactions en présence d’un écoulement, couplage avec les ciné-tiques de dissolution de l’oxygène). La complexité du processus de thermo-oxydation a été réduite en considérant les hydroperoxydes comme une forme de stockage organique de l’oxygène, qui propage l’oxydation dans des régions en anoxie. Leur décomposition produit de nombreux composés de scission, dont la nature est influencée par les conditions locales de température et de concentration en oxygène. La prise d’huile a été décrite comme le bilan net entre l’huile charriée au moment du retrait et l’huile égouttée. L’égouttage a été étudié sur des barreaux métalliques et des produits réels. Il se conduit à la formation de quatre à huit gouttes en quelques secondes. Les cinétiques de drainage anisothermes ont été prédites par un modèle mécanistique. Le mécanisme spécifique de prise d’huile en cours de friture a été aussi analysé ; il se produit uniquement dans le cas des produits préfrits congelés. / Batch deep-frying has been investigated within the collaborative project FUI Fry’In (ref. AAP17, 2014-2018) with the aim of proposing breakthrough innovations for small and medium size appliances. The PhD thesis was part of the project and focused on two specific adverse effects of deep-frying on food products: oil thermo-oxidation responsible for break-down products and off-flavors, and oil pickup process usually favored relatively to oil dripping. The work was carried out by combing direct measurements (FTIR-ATR spectroscopy and imaging, photoionization, DSC measurements, fast imaging…) and multiscale modeling (oil flow and oil dripping during product re-moval, Lagrangian description of reactions in aniso-thermal flows, coupling with oxygen dissolution kinetics). The complex problem of thermo-oxidation was split into simpler mechanisms by noticing that hydroperoxides are a kind of long-lived form of or-ganic oxygen, which trigger propagation in deep re-gions under anoxia. Their decomposition lead to various scission products, which were shown to be in-fluenced by both local temperature and oxygen con-centration. Oil uptake was described as the net balance between the amount of dragged oil during product removal and oil dripping at the tips of the product. The dripping process studied on both metal-lic sticks and real products occurs in less than few seconds and leads to a formation of four to eight drop-lets. The detailed drainage kinetics in anisothermal conditions were captured and predicted with the pro-posed mechanistic models. The specific mechanism of oil uptake during the immersion stage was eluci-dated and was shown to occur only in parfried frozen products.

Un fait marquant et essentiel de la maladie neurodégénérative d’Alzheimer est la formation de plaques amyloïdes dans le cerveau, résultat de l’agrégation des protéines amyloïde-β (Aβ1-40/1-42). Le développement de nouveaux médicaments requiert la compréhension des mécanismes de formation des fibres amyloïdes et la connaissance de la structure et dynamique des oligomères métastables qui sont les vecteurs principaux de la neurotoxicité. Parce que les simulations atomistiques en solvant explicite ne peuvent pas être réalisées sur de grands systèmes pour des temps très longs, nous avons opté pour un modèle protéique gros grain (CG) avec un solvant implicite. Nous nous sommes intéressés dans ces travaux de thèse à clarifier le rôle d’interactions hydrodynamiques(HI) dans la dynamique de formation des agrégats du peptide Aβ(16-22), connu pour former également des fibres amyloïdes. Ces interactions sont essentielles pour modéliser,dans un solvant implicite, les processus se produisant dans des environnements cellulaires très encombrés. Notre approche est basée sur une méthode multi-échelle et multi-physique qui couple les techniques Lattice Boltzmann et de dynamique moléculaire(LBMD). Dans notre système, les interactions médiées par le solvant aqueux sont incluses naturellement. Pour le système moléculaire, nous avons choisi le modèle gros grain à haute résolution OPEP (Optimized Potential for Efficient Protein structure prediction). Pour la première fois, nous avons effectué des simulations quasi tout-atome pour de très grands systèmes contenant des milliers de peptides Aβ ( 16-22). Après avoir correctement réglé le paramètre clé de notre couplage afin d’obtenir la diffusivité expérimentale des monomères et des oligomères du peptide Aβ ( 16-22), nous avons démontré que les HI accélèrent le processus d’agrégation pour des systèmes de taille moyenne (100 Aβ (16-22) peptides) et grande (1000 Aβ (16-22) peptides). Une caractérisation détaillée de la taille des clusters et de l’organisation structurelle des peptides est présentée. Enfin,nous avons examiné comment la concentration affecte la première phase d’agrégation des peptides et leurs structures. / The self-assembly of misfolded amyloid-β (Aβ 1-40/1-42) proteins into insoluble fibrils is strongly linked to the pathogenesis of Alzheimer’s disease (AD). The development of new drugs requires the understanding of the mechanisms leading to fibril formation, and the knowledge of the dynamics and structures of the early metastable oligomers which are the main neurotoxic species. Because atomistic simulations in explicit solvent cannot be performed on very large systems for a significant time scale, we resort to a coarse grained (CG) protein model with an implicit solvent. Our investigation enlightens the role of hydrodynamic interactions (HI) in the kinetics of β-amyloidogenesis, interactions which are essential, when an implicit solvent is used, to model processes occurring in highly crowded like-cell environments, among others.Our approach is based on a multi-scale and multi-physics method that couples Lattice Boltzmann and Molecular Dynamics (LBMD) techniques. In our scheme the solvent- mediated interactions are included naturally. As a first step, we focus on Aβ (16-22) peptide, known to form amyloid fibril alone, and we adopt the high resolution CG OPEP (Optimized Potential for Efficient Protein structure prediction) model, developed in our laboratory. For the first time, we have performed quasi-all-atom simulations for very large systems containing thousands of Aβ (16-22) peptides. After the correct tuning of the key parameters of our coupling in order to obtain the experimental diffusivity of Aβ (16-22) monomer and small oligomers, we have demonstrated that HI speed up the aggregation process of medium (100 peptides) and large (1000 peptides) systems. A detailed characterization of the fluctuating clusters along the trajectories is presented in terms of their sizes and the structural organization of the peptides. Finally, we have investigated how changes in the concentration affect the early aggregation phase of the peptides and their structures.