Cette thèse est consacrée à l’étude des modèles microscopiques pour la dérivation de la conduction de la chaleur. Démontrer rigoureusement une équation diffusive macroscopique à partir d’une description microscopique du système est à aujourd’hui encore un problème ouvert. On étudie un système décrit par l’équation de Schrödinger linéaire discrète (DLS) en dim 1, perturbé par une dynamique stochastique conservative. On peut montrer que le système a une limite hydrodynamique donnée par la solution de l’équation de la chaleur. Quand le système est rattaché aux bords à deux réservoirs de Langevin à deux différents potentiels chimiques, on peut montrer que l’état stationnaire, dans la limite vers l'infinie, satisfait la loi de Fourier. On étudie une chaine des oscillateurs anharmonique immergée en un réservoir de chaleur avec un gradient de température. On exerce une tension, variable dans le temps, à une des deux extrémités de la chaine, et l’autre reste fixe. On montre que sous un changement d’échelle diffusive dans l’espace et dans le temps, la distribution d’étirement de la chaine évolue selon un équation diffusive non-linéaire. On développe des estimations qui reposent sur l’hypocoercitivité entropique. La limite macroscopique peut être utilisée pour modéliser les transformations thermodynamique isothermiques entre états stationnaire de non-équilibre. / The object of research of this thesis is the derivation of heat equation from the underlying microscopic dynamics of the system. Two main models have been studied: a microscopic system described by the discrete Schrödinger equation and an anharmonic chain of oscillators in presence of a gradient of temperature. The first model considered is the one-dimensional discrete linear Schrödinger (DLS) equation perturbed by a conservative stochastic dynamics, that changes the phase of each particles, conserving the total norm (or number of particles). The resulting total dynamics is a degenerate hypoelliptic diffusion with a smooth stationary state. It has been shown that the system has a hydrodynamical limit given by the solution of the heat equation. When it is coupled at the boundaries to two Langevin thermostats at two different chemical potentials, it has been proven that the stationary state, in the limit to infinity, satisfies the Fourier’s law. The second model considered is a chain of anharmonic oscillators immersed in a heat bath with a temperature gradient and a time varying tension applied to one end of the chain while the other side is fixed to a point. We prove that under diffusive space-time rescaling the volume strain distribution of the chain evolves following a non-linear diffusive equation. The stationary states of the dynamics are of non-equilibrium and have a positive entropy production, so the classical relative entropy methods cannot be used. We develop new estimates based on entropic hypocoercivity, that allows to control the distribution of the positions configurations of the chain. The macroscopic limit can be used to model isothermal thermodynamic transformations between non-equilibrium stationary states. CEMRACS project on simulating Rayleigh- Taylor and Richtmyer-Meshkov turbulent mixing zones with a probability density function method at last.

Les travaux présentés s’intéressent à la réponse vibratoire et au champ acoustique émis par une structure immergée et excitée par une couche limite turbulente, dans le domaine des bas nombres d’ondes et pour un nombre de Mach faible. Ce travail s’inscrit dans la problématique d’amélioration de laprédiction du bruit rayonné dans ce type de configurations, et peut trouver son application à la discrétion acoustique des navires, ou à la caractérisation du bruit rayonné par des structures externes excitées par un écoulement.Numériquement, une analyse modale de la réponse de la structure en formulation (u,p,φ) est réalisée à l’aide du code élément finis Code_Aster. L’excitation est modélisée par une somme d’ondes planes de pression dont la densité spectrale est obtenue à partir des modèles d’excitation pariétale disponibles dans la littérature. Une analyse harmonique sur base modale est réalisée pour chaque cas de chargement.Cette approche permet la prise en compte du couplage fluide-structure dans le cas d’un fluide lourd et présente l’avantage de s’affranchir des hypothèses généralement faites, de fluide léger et d’orthogonalité des déformées modales.Les résultats issus de la modélisation numérique sont comparés à des données expérimentales, concernant le comportement vibratoire d’un dispositif constitué d’une plaque plane excitée par un écoulement généré en tunnel hydrodynamique. Les résultats numériques et expérimentaux observés sont proches,qu’il s’agisse du comportement global, du niveau spectral moyen en déplacement ou du niveau de pression acoustique mesuré. En complément, l’influence de défauts, constitués de marche montantes et descendantes de hauteur inférieure à l’épaisseur de la couche limite, sur l’excitation et la réponse de la structure est explorée expérimentalement. / The following work consist in the study of the vibroacoustic response of a structure submerged in fluid, under a turbulent boundary layer flow, the response of the structure is driven by the low wavenumber behaviour, for a small Mach number. This work aims at providing better means of predicting the noise radiated in such setups, mainly regarding stealthiness of ships and submarines and noise radiated by outer structures.A numerical modal analysis based on the (u,p,φ) formulation available in the finite element software Code_Aster is performed. The pressure induced by the boudary layer is then described as a sum of plane waves and several harmonical analysis are performed on the reduced problem, projected on the (u,p,φ) modal basis, one for each term of the sum. This allows us to account for the fluid-structure interaction (inertial and acoustic) in confined and infinite fluid domains. Most numerical models found in scientific papers are making the assumption of a light fluid, or a fluid loaded plate, thus not taking clearly into account the fluid-strucure interaction or only the inertialpart. Here the interaction due to the acoustic field radiated by the plate is fully accounted for.The validity of the proposed numerical method is assesed and numerical results are compared to data obtained from an experimental setup used within a hydrodynamic tunnel. Numerically, a good reproduction of the behaviour of the plate is obtained, both in terms of displacement and spectral levels. The acoustic levels are also compared to their numerical counterparts at the position of the transducer. Moreover, an experimantal analysis is performed, for backward and forward steps of height smaller than the thickness of the boundary layer, in order to investigate the influence of such configurations on the boundary layer excitation and on the vibroacoustic response.

La non-reproductibilité numérique apparait dans divers domaines d'application de la simulation HPC. En effet, les différentes distributions d'un calcul parallèle peuvent mener à des résultats numériques différents, à cause des particularités de l'arithmétique flottante. Le besoin de reproductibilité numérique est motivé pour le débogage, le test et la validation des codes de calcul scientifique. Nous nous intéressons aux simulations par éléments finis en hydrodynamique implémentées dans le logiciel openTelemac qui est largement utilisé pour des applications industrielles et scientifiques. Nous identifions et analysons les sources de cette non-reproductibilité. Nous définissons et implementons comment récupérer la reproductibilité numérique de deux modules d'openTelemac. Nous mesurons que le sur-coût en terme de temps de calcul de la version reproductible est tout à fait raisonnable en pratique. / HPC simulations in various scientific domains suffer from failures of numerical reproducibility because of floating-point arithmetic peculiarities. Different distributions of a parallel computation may yield different numerical results. Numerical reproducibility is a requested feature to facilitate the debug, the validation and the test of industrial or large software. In this thesis, we focus on the openTelemac software that implements finite element simulation for industrial and scientific hydrodynamics. We identify and analyze the sources of this reproducibility failure. We define and implement how to recover numerical reproducibility in two openTelemac modules. We also measure that the running time extra-cost of the reproducible version is reasonable enough in practice.

Non disponible / Non disponible

Le monoxyde de carbone déshydrogénase (CODH) catalyse la réduction réversible du dioxyde de carbone par son site actif. En utilisant une méthode électrochimique appelée voltammétrie de film protéique, nous étudions le mécanisme enzymatique de CODH en immobilisant l'enzyme à une surface d'électrode de graphite où le transfert direct d'électrons est possible. Traditionnellement, pour empêcher la déplétion du substrat à l'électrode, les électrochimistes utilisent des électrodes tournantes (RDE). Cependant, comme la CODH est très active, même la RDE ne peut pas empêcher l'épuisement, ce qui masque les caractéristiques cinétiques importantes de l’enzyme. Nous ne pouvons pas résoudre le problème avec RDE, puisque nous l’utilisons déjà à la vitesse maximum. Par conséquent, nous devons concevoir une nouvelle cellule électrochimique. Pour cela, nous avons utilisé des simulations de dynamique des fluides computationnelles pour explorer diverses géométries afin d'en trouver une appropriée. Nous avons commencé par valider notre méthode numérique avec la solution théorique bien définie de la cellule réelle de RDE. Après la bonne validation, nous avons déterminé les vitesses de transport de masse au sein de plusieurs géométries et à basé sur l'optimisation des paramètres géométriques, nous avons atteint notre conception appropriée. Ce nouveau prototype a une électrode graphite uniformément accessible avec un taux de transport trois fois plus rapide que le RDE à des vitesses de solution acceptables. Nous avons construit, mis en place avec succès le système pour caractériser ses performances de transport, et trouvé un excellent accord entre les résultats numériques et expérimentaux / Carbon monoxide dehydrogenase (CODH) catalyzes the reversible reduction of carbon dioxide by its active site. Thus, CODH participates in the first step of fuel production. Using an electrochemical method called protein film voltammetry, we study the enzymatic mechanism of CODH by immobilizing the enzyme at a graphite electrode surface where direct electron transfer is possible. Traditionally, to prevent depletion of the substrate at the electrode, electrochemists use rotating electrodes (RDE). However, since CODH is very active, even RDE cannot prevent depletion, which masks the important kinetic characteristics of the enzyme and complicates the analysis of the enzymatic response.We cannot solve the problem with RDE, since we already use it at maximum speed. Therefore, we must completely change our approach and design a new electrochemical cell. For this, we used computational fluid dynamics (CFD) simulations to explore various geometries to find a suitable one. We began by validating our numerical method with the well-defined theoretical solution of the real cell of RDE. After good validation, we determined the mass transport velocities within several proposed geometries of the flow cell of hydrodynamic channel and jet electrodes. Based on the optimization of geometric parameters, we have achieved our proper design of jet electrode. This new prototype has a uniformly accessible graphite electrode with a transport rate three times faster than the RDE at acceptable solution speeds. We have successfully built and implemented the system to characterize its transport performance. We found an excellent agreement between the numerical and experimental results

L'étude des dynamiques collectives, observables chez de nombreuses espèces animales, a motivé dans les dernières décennies un champ de recherche actif et transdisciplinaire. De tels comportements sont souvent modélisés par de la matière active, c'est-à-dire par des modèles dans lesquels chaque individu est caractérisé par une vitesse propre qui tend à s'aligner avec celle de ses voisins.Dans un modèle fondateur proposé par Vicsek et al., ainsi que dans de nombreux modèles de matière active liés à ce dernier, une transition de phase entre un comportement chaotique à forte température, et un comportement global et cohérent à faible température, a été observée. De nombreuses preuves numériques de telles transitions de phase ont été obtenues dans le cadre des dynamiques collectives. D'un point de vue mathématique, toutefois, ces systèmes actifs sont encore mal compris. Plusieurs résultats ont été obtenus récemment sous une approximation de champ moyen, mais il n'y a encore à ce jour que peu d'études mathématiques de modèles actifs faisant intervenir des interactions purement microscopiques.Dans ce manuscrit, nous décrivons un système de particules actives sur réseau interagissant localement pour aligner leurs vitesses. L'objet de cette thèse est l'obtention rigoureuse, à l'aide du formalisme des limites hydrodynamiques pour les gaz sur réseau, de la limite macroscopique de ce système hors-équilibre, qui pose de nombreuses difficultés techniques et théoriques. / Collective dynamics can be observed among many animal species, and have given rise in the last decades to an active and interdisciplinary field of study. Such behaviors are usually modeled by active matter, in which each individual is self-driven and tends to align its velocity with that of its neighbors.In a classical model introduced by Vicsek & al., as well as in numerous related active matter models, a phase transition between chaotic behavior at high temperature and global order at low temperature can be observed. Even though ample evidence of these phase transitions has been obtained for collective dynamics, from a mathematical standpoint, such active systems are not fully understood yet. Some progress has been achieved in the recent years under an assumption of mean-field interactions, however to this day, few rigorous results have been obtained for models involving purely local interactions.In this manuscript, we describe a lattice active particle system interacting locally to align their velocities. This thesis aims at rigorously obtaining, using the formalism developed for hydrodynamic limits of lattice gases, the scaling limit of this out-of-equilibrium system, for which numerous technical and theoretical difficulties arise.

Ce sujet de thèse s’inscrit dans le cadre des études actuellement menées au CEA/DAM concernant des schémas numériques ALE (Arbitrary-Lagrangian-Eulerian)de type « Lagrange + Projection », dans le contexte des simulations hydrodynamiques mutli-matériaux en grandes déformations. Ces schémas doivent respecter les équations de conservation de la masse, de la quantité de mouvement et de l’énergie totale.Les schémas décalés en temps et en espace sont très utilisés dans les codes industriels. Ils sont robustes et permettent une bonne approximation des comportements complexes, mais sont connus pour ne pas conserver exactement l’énergie totale. Cela pose un problème dans le traitement des chocs, sur maillages raffinés ou dans la simulation des milieux réactifs.En 2016, des travaux originaux on été proposés par A. Llor et. al. pour rendre conservatif ce type de schéma dans un contexte lagrangien (sans projection), notamment en proposant une correction pour retrouver la conservation de l’énergie totale.Le travail de cette thèse a été d’étendre ce schéma lagrangien dans un contexte ALE multi-matériaux (avec interface), en garantissant la conservation de toutes les quantités, le respect du second principe de la thermodynamique et la robustesse. De nombreux cas tests ont été menés (en 2D plan et en 2D axisymétrique) et comparés aux méthodes existantes afin de montrer la pertinence de cette approche. / This PhD subject comes within actual studies managed by CEA/DAM about ALE (Arbitrary-Lagrangian-Eulerian) schemes (with a splitting of Lagrangian and Remapping steps) in the context of hydrodynamic simulations. These numerical schemes have to respect mass, momentum and total energy conservation, which are the fundamental equations of the studied systems.Space- and Time-Staggered are widely used in industrial codes for their simplicity androbustness despite their known lack of exact energy conservation. This is a major drawbackin presence of strong shocks. Among all existing schemes, none of them meet the expectations of robustness, conservation,thermodynamic consistency (both shocks and relaxations capture), accuracy andadaptibility to complex behaviors. Recently, some novel works have been proposed by A.Llor et. al. in order to make conservative this type of scheme in a Lagrangian context (without remapping step). Current remap methods, necessary in large deformations, donot guarantee simultaneously total energy conservation and thermodynamic consistency.This work aims at extending this conservative Lagrangian space- and time-staggeredscheme to a multi-material ALE methodology, keeping its good properties (conservation,accuracy, thermodynamic consistency, robustness) intact. Classical, but demanding, test cases have been performed (both in plane and axisymmetric 2D geometries) and have been compared to existing numerical methods in order to assess the relevance of our approach.

Cette thèse porte sur le problème de transport optimal incompressible, un problème introduit par Brenier à la fin des années 80 dans le but de décrire l’évolution d’un fluide incompressible et non-visqueux de façon lagrangienne, ou autrement dit en fixant l’état initial et l’état final de ce fluide, et en minimisant une certaine fonctionnelle sur un ensemble de dynamiques admissibles. Ce manuscrit contient deux parties.Dans la première, on étudie la dépendance du problème de transport optimal incompressible par rapport aux données. Plus précisément, on étudie la dépendance du champ de pression (le multiplicateur de Lagrange associé à la contrainte d’incompressibilité) par rapport à la mesure γ prescrivant l’état initial et l’état final du fluide. On montre au Chapitre 2 par des méthodes variationnelles que le gradient de la pression, en tant qu’élément d’un espace proche du dual des fonctions C^1, dépend de γ de façon hölderienne pour la distance de Monge-Kantorovic. En contrepartie, on montre au Chapitre 4 que pour tout r > 1, le champ de pression dans l'espace de Lebesgue L^r_t L^1_x ne peut pas être une fonction lipschitzienne de γ. Ce résultat est lié au caractère mal-posé de l’équation d’Euler cinétique, une équation cinétique s’interprétant comme l’équation d’optimalité dans le problème de transport optimal incompressible, à laquelle le Chapitre 3 de cette thèse est dédié.Dans une seconde partie, on s’intéresse à la régularisation entropique du problème de transport optimal incompressible, introduit par Arnaudon, Cruzeiro, Léonard et Zambrini en 2017 sous le nom de problème de Brödinger. On prouve au Chapitre 5 que comme dans le cas du transport optimal incompressible, on peut associer à toute solution un champ scalaire de pression agissant comme multiplicateur de Lagrange pour la contrainte d’incompressibilité. On montre ensuite au Chapitre 6 que lorsque le paramètre de régularisation tend vers zéro, le problème de Brödinger converge vers le problème de transport optimal incompressible au sens de la Γ-convergence, et avec convergence des champs de pression. Ce dernier chapitre est issu d'un travail effectué en collaboration avec L. Monsaingeon. / This thesis focuses on Incompressible Optimal Transport, a minimization problem introduced by Brenier in the late 80's, aiming at describing the evolution of an incompressible and inviscid fluid in a Lagrangian way , i.e. by prescribing the state of the fluid at the initial and final times and by minimizing some functional among the set of admissible dynamics. This text is divided into two parts.In the first part, we study the dependence of this optimization problem with respect to the data. More precisely, we analyse the dependence of the pressure field, the Lagrange multiplier corresponding to the incompressibility constraint, with respect to the endpoint conditions, described by a probability measure γ determining the state of the fluid at the initial and final times. We show in Chapter 2 by purely variational methods that the gradient of the pressure field, as an element of a space that is close to the dual of C^1, is a Hölder continuous function of γ for the Monge-Kantorovic distance. On the other hand, we prove in Chapter 4 that for all r>1 the pressure field, as an element of L^r_t L^1_x, cannot be a Lipschitz continuous function of γ for the Monge-Kantorovic distance. This last statement is linked to an ill-posedness result proved in Chapter 3 for the so-called kinetic Euler equation, a kinetic PDE interpreted as the optimality equation of the Incompressible Optimal Transport problem.In the second part, we are interested in the entropic regularization of the Incompressible Optimal Transport problem: the so-called Brödinger problem, introduced by Arnaudon, Cruzeiro, Léonard and Zambrini in 2017. On the one hand, we prove in Chapter 5 that similarly to what happens in the Incompressible Optimal Transport case, to a solution always corresponds a scalar pressure field acting as the Lagrange multiplier for the incompressibility constraint. On the other hand, we prove in Chapter 6 that when the diffusivity coefficient tends to zero, the Brödinger problem converges towards the Incompressible Optimal Transport problem in the sense of Gamma-convergence, and with convergence of the pressure fields. The results of Chapter 6 come from a joint work with L. Monsaingeon.

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.

Considérant l'effet des parois sur la croissance des bulles, nous avons récemment utilisé des plaques ondulées comme cloisons dans des lits minces fluidisés gaz-solide à multi-compartiments. Des analyses approfondies de la dynamique des bulles et du transfert de chaleur paroi-lit dans des lits fluidisés à cloisons planes (FWBFB) ou ondulées (CWBFB) ont été effectuées pour une variété de déclinaisons pariétales et des conditions d'opération couvrant une large gamme d'angles d'ondulation (θ = 120 °, 90 °), de distances moyennes inter-paroi (C), de hauteurs initiales de repos du lit et des rapports de vitesses superficielles du gaz aux vitesses minimales de barbotage, Ug/Umb. Il a été observé que le débit de gaz nécessaire pour amorcer le barbotage était plus faible dans le cas du CWBFB. Par ailleurs, un réseau de zones cou (distance minimale) et hanche (distance maximale) du CWBFB a également favorisé la rupture des bulles, une haute fréquence des bulles, une faible vitesse de montée des bulles et donc tous convergeant vers une meilleure distribution de gaz. Le CWBFB a offert un fonctionnement de fluidisation gaz-solide stable, une faible hauteur de désengagement de transport (HDT) comparé à FWBFB. CWBFB offert significanlty supérieur coefficient de transfert de mur-à-lit chaleur comapred à FWBFB. Des simulations complètes de type mécanique des fluides numériques (CFD) transitoires Euler-Euler 3-D ont également été menées, ce qui a permis de mieux comprendre les effets de parois ondulées sur l'augmentation de la force de traînée sur les particules dans des zones à hautes pressions convergentes-divergentes des parois ondulées. Ceux-ci ont indiqué des changements notables dans le régime de fluidisation en remplaçant les parois planes par des parois ondulées et ont en outre révélé que les zones cou étaient responsables de la création des instabilités comparativement aux zones hanche. La moyenne dans le temps des contours de la fraction volumique simulée de gaz a corroboré avec les résultats expérimentaux que le CWBFB a offert la meilleure distribution de gaz par rapport à FWBFB. Les profils axiaux de la moyenne dans le temps de la fraction volumique solide simulée ont montré que CWBFB réduit de façon nette la HDT. / With the endeavor of approaching an ideal allothermal gasifier, recently our group proposed a reactor concept of allothermal cyclic multi-compartment bubbling fluidized beds for biomass gasification with steam. The concept consisted of multiple intercalated parallelepipedic slim gasification and combustion compartments to enhance unit heat integration and thermal efficiency while preventing contact between flue gas and syngas to generate a N2-free high-quality biosyngas. However, the efficiency of contacting between gas and particles in bubbling fluidized beds is dictated to a large extent by the bubble dynamics which impacts mixing, heat and mass transfers. Literature showed that the decrease in clearance between flat walls for slim fluidization enclosures or in diameter for cylindrical vessels would make fluidized beds very sensitive to wall effects and prone to operate in slug flow regime. Since the occurrence of slugging in multi-compartment slim beds could reduce their thermal and chemical efficiency, the objective of current work was to devise suitable strategies in treating the incipient bubbles to suppress the slugging behavior of bed. By considering the effect of walls on bubble growth, we recently employed corrugated plates as separating walls in slim multi-compartment gas solid fluidized beds. Thorough analyses of bubble dynamics and wall-to-bed heat transfer in flat- (FWBFB) and corrugated- (CWBFB) wall bubbling fluidized beds were performed for a variety of wall declinations and operating conditions covering a range of corrugation angles (θ=120o, 90o), average inter-wall clearances (C), initial rest bed heights (Hi) and ratios of gas superficial velocity to minimum bubbling velocity, Ug/Umb. It was observed that gas flowrate required to achieve the incipient bubbling condition was lower in case of CWBFB. A network of neck (minimum clearance) and hip (maximum clearance) locations in CWBFB also promoted bubbles breakup, higher bubble frequency, lower bubble rise velocity and thus all converging into a better gas distribution. CWBFB offered stable gas-solid fluidization operation and lower transport disengagement height as compared to FWBFB. During the experimental work, digital image analysis technique and fast response heat flux probes were employed to study the effects of operating and geometrical parameters on bubble dynamics and wall-to-bed heat transfer. Two artificial neural network correlations valid both for FWBFB and CWBFB were recommended for the estimation of bubble frequency and size (equivalent diameter). Full 3-D transient Euler-Euler CFD simulations with kinetic theory of granular flow were also carried out which helped shaping an understanding of the effects of corrugated walls on increasing the drag force on particles in the converging-diverging high-pressure zones in corrugated walls. The dynamic fluctuations in the simulated solid phase volume fraction, granular temperature and granular pressure were monitored to determine their standard deviations. These revealed notable shifts in the fluidization regime by replacing flat walls with corrugated walls and further revealed that necks were responsible for inception of instabilities as compared to hips. Time averaged contours of simulated gas volume fraction corroborated with experimental findings that CWBFB offered better gas distribution as compared to FWBFB. Axial profiles of simulated time averaged solid volume fraction and granular temperature showed that CWBFB significantly reduced the transport disengagement height as compared to FWBFB.