Liquid film flow is an important flow type in many applications of process engineering. For supporting experiments, theoretical and numerical investigations are required. The present state of the art is to model the liquid film flow with Navier--Stokes-based methods, whereas the lattice Boltzmann method is employed here. The final model has been developed within this treatise by means of a two-phase flow and a heat transfer model, and boundary and initial conditions. All these sub-models have been applied to simple test cases. It could be found that the two-phase model is capable of solving flow phenomena with a large density ratio which has been shown impressively in conjunction with wall boundary conditions. The heat transfer model was tested against spectral method results with a transient non-uniform flow field. It was possible to find optimal parameters for computation. The final model has been applied to steady-state film flow, and showed very good agreement to OpenFOAM simulations. Tests with transient film flow demonstrated that the model is also able to predict these flow phenomena. / Flüssigkeitsfilmströmungen kommen in vielen verfahrenstechnischen Prozessen zum Einsatz. Zur Unterstützung von Experimenten sind theoretische und numerische Untersuchungen nötig. Stand der Technik ist es, Navier--Stokes-basierte Modelle zu verwenden, wohingegen hier die Lattice-Boltzmann-Methode verwendet wird. Das finale Modell wurde unter Verwendung eines Zweiphasen- und eines Wärmeübertragungsmodell entwickelt und geeignete Rand- und Anfangsbedingungen formuliert. Alle Untermodelle wurden anhand einfacher Testfälle überprüft. Es konnte herausgefunden werden, dass das Zweiphasenmodell Strömungen großer Dichteunterschiede rechnen kann, was eindrucksvoll im Zusammenhang mit Wandrandbedingungen gezeigt wurde. Das Wärmeübertragungsmodell wurde gegen eine Spektrallösung anhand eines transienten und nichtuniformen Strömungsproblemes getestet. Stationäre Filmströmungen zeigten sehr gute Übereinstimmungen mit OpenFOAM-Lösungen und instationäre Berechungen bewiesen, dass das Model auch solche Strömungen abbilden kann.

Cette étude s'intéresse aux risques associés à la formation d'une ATmosphère EXplosive (ATEX) née d'une fuite d'hydrogène et de sa dispersion dans l'air ambiant. La fuite a été modélisée par un jet turbulent à densité variable libre, impactant sur une sphère de diamètre 20mm ou sur une plaque plane. Dans un premier temps, les champs de vitesses et de concentration ont été obtenus expérimentalement en proche sortie grâce à des mesures de Vélocimétrie par Images de Particules (PIV) et de Fluorescence Induite par Plan Laser sur l'acétone (PLIF). La turbulence et le mélange ont été caractérisés pour le cas d'un jet libre ou en présence d'un obstacle. A partir de ces mesures, la structure générale de l'écoulement a été étudiée à partir des champs moyens et fluctuants par comparaison avec les données de la bibliographie. Puis, les données issues des fluctuations ont été analysées statistiquement par l'étude des fonctions de densité de probabilité du scalaire. Ces travaux se sont poursuivis avec la mise en relation des résultats expérimentaux avec ceux obtenues par des simulations numériques DNS (Direct Numerical Simulation) utilisant la méthode Boltzmann sur Réseau (LBM) d'un scalaire passif dans un jet d'air. Cette étude a permis de recueillir et d'analyser des données supplémentaires sur le mélange d'un jet à masse volumique variable libre ou impactant. Ces données sont directement applicables à la maitrise des risques liés aux fuites d'hydrogène. / This study examines the risks associated with the formation of an explosive atmosphere from a hydrogen leak and its dispersion into the air. We considered the leak as a turbulent jet with density variable, free and impinging a $20,mm$ diameter sphere or a flat plate. Firstly, velocity and scalar fields have been measured experimentally in the near field through Particle Image Velocimetry (PIV) and acetone Planar Laser Induced Fluorescence (LIF). Turbulence and mixing have been defined in the case of free jet and impinging jet. From this measurements, the flow structure has been presented from the mean and fluctuating flow measurements by comparison with literature data. Next, the fluctuation scalar fields are studied with the probability density function method. Finally, a comparison has been conducted between the experiments and direct numerical simulation (DNS) of turbulence based on the lattice Boltzmann method (LBM) for passive scalar in air jet. This study is gathering and analyzing data on the mixing of jet with density variable, free and impinging jet. Such data is directly useful to identify and control risks incurred due to hydrogen leak.

La propreté inclusionnaire reste un enjeu majeur en élaboration des métaux par voie liquide. La flottation, principal procédé retenu en métallurgie secondaire pour éliminer les particules d’inclusions, consiste à injecter des bulles de gaz au sein du réacteur. Lors de leur ascension, les bulles vont capter les plus grosses inclusions et favoriser la collision et l’agrégation des particules. Dans le but de quantifier les phénomènes influents à l’échelle des inclusions sur la dynamique d’agrégation entre deux inclusions à proximité des bulles, un modèle numérique 3D a été développé. L’écoulement local est modélisé par un cisaillement plan permanent et résolu par une méthode de Boltzmann sur réseau. Le couplage entre les particules et le fluide a été assuré par une méthode de frontière immergée permettant de calculer la perturbation hydrodynamique engendrée par la présence des particules et de mettre à jour les interactions entre particules pour leur suivi lagrangien. Les simulations numériques réalisées ont mis en évidence que les effets hydrodynamiques ont une influence non négligeable sur le comportement des inclusions. Des sections efficaces de collision ont pu être extraites, à partir desquelles ont été calculés des noyaux d’agrégation, données macroscopiques rendant compte des effets à petite échelle. Une première application de ce travail a été menée avec le calcul des fréquences d’agrégation d’un train de bulle dans un réacteur canal à partir de résultats de simulations DNS. Les noyaux d’agrégation ont également été exploités en vue de déterminer, à partir de résultats RANS de l’hydrodynamique d’une poche d’acier, l’évolution de la concentration d’inclusions par un bilan de population global. / Inclusion cleanliness remains a major challenge faced in process metallurgy in liquid phase. Flotation, the main process used in secondary metallurgy to remove inclusions, consists in injecting gas bubbles into the reactor. Rising gas bubbles entrap the biggest inclusions at their surface or in their wake. Besides, they promote collision and aggregation among particles. A 3D numerical model has been developed in order to quantify the roles of the prevailing phenomena on aggregation dynamics between inclusions in the vicinity of bubbles. At inclusion (mesoscopic) scale, the turbulent flow is locally modeled by a steady plane shear flow which is solved using a lattice-Boltzmann method. The coupling between both liquid and solid phases is ensured using an immersed boundary method. This method resolves the hydrodynamic perturbation induced by particles, and hence their interactions that are, in turn, used to update their Lagrangian tracking. The conducted numerical simulations bring out the influence of hydrodynamic effects on inclusion behavior. Collision cross sections have been determined from which ensuing aggregation kernels have been calculated. Such cross sections could provide macroscopic models to represent local particle dynamics. A first application of these results is presented to calculate aggregation frequencies in bubble swarms in a channel flow reactor that was simulated using DNS. Additionally, evolution of inclusion populations in molten steel has been determined from RANS simulation of a liquid steel ladle by means of a global population balance implementing the aggregation kernels determined in the present work.

Ce travail de thèse de doctorat, effectué au sein du Département de Chimie Appliquée (au Laboratoire de Génie Chimique, LGC), vise à réaliser létude de lécoulement dun fluide dans des géométries complexes rencontrées en Génie Chimique (les empilages aléatoires et structurés, souvent utilisés dans des colonnes dadsorption et de distillation) par la méthode de Boltzmann sur réseau (Lattice Boltzmann Method, LBM) et à développer un outil numérique de simulation sous la forme dun code de calcul tridimensionnel. Dans le but daméliorer les performances de ces éléments complexes, qui conduisent à une efficacité et une sélectivité accrue des procédés se produisant dans les colonnes de distillation et dadsorption, il est donc intéressant dexaminer avec plus de précision la description des phénomènes de transport se produisant à léchelle locale dans ces milieux. Dans un premier temps, les concepts et les équations macroscopiques qui dérivent de lécoulement dun fluide ont été succinctement passés en revue, ainsi que le modèle de turbulence qui a été utilisé. Nous présentons ensuite, les éléments théoriques de la méthode de Boltzmann sur réseau, après quelques détails sur le modèle Gaz sur réseau, dont elle est dérivée. La méthode de Boltzmann sur réseau se démarque des schémas traditionnels en résolvant numériquement une équation basée sur la physique statistique. Les codes des modèles SRTLBE (équation de Boltzmann sur réseau à un seul temps de relaxation) et MRTLBE (équation de Boltzmann sur réseau à plusieurs temps de relaxation) à trois dimensions, de même que des conditions aux limites convenables ont été écrits dans un langage de programmation, puis validées sur des cas types de la littérature : un écoulement laminaire dans un conduit de section carrée et un écoulement de jet turbulent rectangulaire à 3D. Les codes implémentés dans une grille de calcul ont permis de simuler des écoulements de fluides dans des structures très complexes qui ont été obtenues par une technique expérimentale avancée : la microtomographie à rayons X. Létude dune mousse métallique, dun filtre à charbons actifs et dun filtre de capture daérosols ont été présentés. Les résultats numériques ont été comparés avec succès aux mesures expérimentales et à des corrélations de la littérature. Une autre étude a concerné létude de lécoulement à travers deux feuilles de Mellapak 250 Y. Les champs et contours de vitesse ainsi obtenus ont été soigneusement analysés et commentés à la fois en régime laminaire et en régime turbulent. Les résultats hydrodynamiques ont été comparés avec des mesures expérimentales ainsi quavec des valeurs calculées par des corrélations disponibles dans la littérature. Les modèles existant ont été analysés. Enfin, ces simulations (LBM) ont été comparées à celles effectuées dans les mêmes conditions opératoires, avec le code de calcul classique FLUENT.

No description available.

We propose a hypothesis and a numerical model for the formation of branching plagioclase textures visible at both macroscopic (∼cm to ∼m) and microscopic scale within melagabbro of the Isle of Rum, Scotland, based on macroscopic, microscopic observations and relevant geological history. The plagioclase crystals are typically linked as twins and form meshes of planar stellate structures (m-scale) with a large range in geometrical organization from patchy to radiating. Evidence of macroscopic crystal aggregation and alignment is attributed to interfacial free energy minimization at the microscopic scale during growth. Accordingly, a binary immiscible Lattice Boltzmann model was developed to simulate diffusion of simplified plagioclase in the melt phase. Isothermal phase transitions modelled via first order chemical reactions are subsequently coupled with stochastic dynamics at the crystal growth front to simulate energy minimization processes including twinning during crystallization in an igneous environment. The solid phase and the liquid phase are coupled with a temporal flexibility that sets the overall ratio between the rate of diffusion and chemical enrichment in the liquid state and the rate of crystallization. The parameter space of the model is explored extensively, followed by a reasonable transcription of physical parameters and an estimation of other parameters to construct realistic simulation scenarios yielding synthetic plagioclase stellates. The results are presented, analyzed and discussed. They appear to be in reasonable qualitative agreement with observations, and several aspects of the natural stellates such as the stellate spacing and long continuous stretches of plagioclase with epitaxial junctions seem to be in reasonable quantitative agreement with observations.

Ο στόχος της παρούσας εργασίας είναι η περαιτέρω κατανόηση και ποσοτική σύνδεση φαινομένων μεταφοράς που λαμβάνουν χώρα σε πορώδη μέσα με τα αντίστοιχα φαινόμενα στην κλίμακα λίγων πόρων. Η επέκταση των αποτελεσμάτων από την κλίμακα πόρου στην κλίμακα του πορώδους μέσου δεν είναι προφανής και για το λόγο αυτό η τοπολογία και μορφολογία της πορώδους δομής αντιμετωπίζονται εδώ με δίκτυα πόρων, με έμφαση στα φαινόμενα που λαμβάνουν χώρα στις διασταυρώσεις, αλλά και με ψηφιακές αναπαραστάσεις της δομής με βάση μικροφωτογραφίες δείγματος του υλικού. Συγκεκριμένα, στην εργασία αυτή εξετάζεται η διασπορά μορίων διαλυμένης ουσίας σε δίκτυα πόρων, παρουσιάζεται μία καινούργια τεχνική ανακατασκευής ανομοιογενών πορωδών υλικών και αναπτύσσεται μια μέθοδος προσομοίωσης της ροής αερίων δια μέσου ανακατασκευασμένων πορωδών υλικών στη μεταβατική περιοχή ροής όπου η μέση ελεύθερη διαδρομή των μορίων ενός αερίου είναι συγκρίσιμη με το μέγεθος των πόρων οπότε και παύει να ισχύει η συνήθης παραδοχή του συνεχούς. Η επίδραση της ανάμειξης μέσα σε πόρους ή στις διασταυρώσεις πόρων/ρωγμών στη διασπορά μορίων διαλυμένης ουσίας σε πορώδη μέσα ερευνήθηκε μέσα από την ανάπτυξη και χρήση διαφορετικών τεχνικών προσομοίωσης με έμφαση στις λεπτομέρειες της ροής και της μεταφοράς μάζας στην περιοχή της διασταύρωσης. Βρέθηκε ότι μία νέα μέθοδος τυχαίου περιπάτου αναπαράγει με καλή ακρίβεια το συντελεστή διασποράς σε χαμηλές και μεσαίες τιμές του Peclet, χάρη στο γεγονός ότι λαμβάνει υπ’ όψη την ανάντι της ροής κίνηση των σωματιδίων και τους διαφορετικούς χρόνους παραμονής μέσα σε κάθε κλάδο. Παράλληλα αναπτύχθηκε μία καινοτόμος μέθοδος ανακατασκευής πορωδών μέσων. Η τεχνική στηρίζεται στο διφασικό πρότυπο δικτύου Boltzmann, το οποίο περιγράφει την εξέλιξη συστημάτων υγρού-αερίου υπό την επίδραση της διεπιφανειακής τάσης. Ο μηχανισμός αυτός οδηγεί στη δημιουργία συσχετισμένων δομών, όπου τόσο η μορφολογία του πορώδους μέσου όσο και ο βαθμός συσχέτισής του καθορίζονται από τις λειτουργικές παραμέτρους του προτύπου. Η τεχνική εφαρμόστηκε επιτυχώς σε πραγματικό δείγμα εδάφους με αφετηρία την πληροφορία που δίνεται από μία μικροφωτογραφία μίας στατιστικά χαρακτηριστικής τομής του. Τέλος, μελετήθηκε η ροή αερίων σε πορώδη μέσα, σε πεπερασμένους αριθμούς Knudsen, όπου η μέση διάμετρος των πόρων είναι της ίδιας τάξης με τη μέση ελευθέρα διαδρομή των μορίων του αερίου. Η μελέτη έγινε με τη μεσοσκοπική μέθοδο DSMC. Ο έλεγχος της αξιοπιστίας της μεθόδου και της παρούσας υλοποίησής της έγινε μέσω της μελέτης της ισοθερμοκρασιακής ροής αερίου μεταξύ παραλλήλων πλακών. Παράλληλα υπολογίστηκε το δυναμικό ιξώδες αερίου σε συνθήκες υψηλής αραίωσης και παρουσιάστηκε η εξάρτησή του από τον αριθμό Knudsen. Βρέθηκε ότι τα αποτελέσματα προσεγγίζονται ικανοποιητικά από μία αναλυτική έκφραση τύπου Bosanquet που συσχετίζει το αποτελεσματικό ιξώδες με την τιμή του στο όριο του συνεχούς και με τον αριθμό Knudsen. Επιπρόσθετα μελετήθηκε για πρώτη φορά με τη μέθοδο DMSC η ροή αερίων σε υπολογιστικά ανακατασκευασμένες πορώδεις δομές. Επιβεβαιώθηκε το φαινόμενο του Klinkenberg και η γραμμική εξάρτηση του συντελεστή διαπερατότητας από την αντίστροφη πίεση. Τέλος χρησιμοποιήθηκε μια διαφορετική προσέγγιση στο πρόβλημα υπολογισμού της ροής στη μεταβατική περιοχή μέσω ανάπτυξης προτύπου δικτύου Boltzmann, κατάλληλα τροποποιημένου για ροές σε συνθήκες αραίωσης. Το πρότυπο δοκιμάστηκε τόσο στην περίπτωση ροής μεταξύ παραλλήλων πλακών όσο και σε ροή σε πορώδη μέσα όπου η συμφωνία με τη μέθοδο DSMC βρέθηκε πολύ ικανοποιητική. / The aim of the present study is the further understanding and quantification of transport phenomena in porous media and their connection with the phenomena in the scale of a few pores. The extension of the results from the pore-scale to the scale of the porous medium is not obvious and for this reason the representation of the porous medium is treated both with pore-networks and digital reconstruction. Specifically, in this study it is examined the dispersion of molecules of a solute in porous networks, a new reconstruction technique is presented for heterogeneous granular materials and also a methodology is developed for the study of gas flow in reconstructed porous media in the transient regime, where the mean free path of the gas molecules is comparable with the characteristic length of the pores and thus the continuum description is no longer valid. The effect of the mixing in the pores or the junctions of the pores on the dispersion of molecules of a solute in porous media is examined through various simulation techniques with emphasis on the details of the flow and mass transport in the area of the junction. It was found that a new random-walk technique is reproducing with good accuracy the dispersion coefficient for low and average values of the Peclet number, due to the fact that it takes into account the backwards, with respect to the main direction of the flow, movement of the molecules and the different residence time in each branch. Furthermore, a new reconstruction technique was developed for porous media. The technique is based on 2-phase lattice Boltzmann model, which describes the evolution of a gas-liquid system under the influence of the surface tension. This mechanism leads to the creation of correlated structures, where the morphology of the porous medium and the correlation factor are determined by the operating parameters of the model. The technique was applied successfully for the reconstruction of a real soil sample, starting from the information that is solely given from a microphotograph of a statistically adequate section of the material. Finally, the gas flow through porous media was examined at moderate Knudsen numbers, where the mean diameter of the pores is of the same order of magnitude with the mean free path of the gas molecules. The study was done mainly with the mesoscopic DSMC technique. The credibility of the technique was examined through the study of the isothermal gas flow through parallel plates. Additionally, the dynamic viscosity of a gas under rarefaction conditions was calculated and its dependence on the Knudsen number was shown. It was found that the results are approximated satisfactorily with an analytical Bosanquet-type equation that relates the effective viscosity with its value at the continuum limit and with the Knudsen number. Furthermore, it was studied for the first time with the DSMC method the gas flow through reconstructed porous media. The Klinkenberg effect was confirmed and the linear dependence of the permeability coefficient on the inverse pressure was shown. Finally an alternative approach was used for the calculation of gas flow though porous media in the transient regime through the development of a lattice Boltzmann model suitably modified for rarefied gas flows. The model was tested for the case of flow through parallel plates as well as for the case of flow through porous media and the agreement with the DSMC method was very satisfactory.

Cette thèse porte sur l'étude numérique des fluctuations extrêmes de la force de traînée exercée par un écoulement turbulent sur un corps immergé.Ce type d'évènement, très rare, est difficile à caractériser par le biais d'un échantillonnage direct, puisqu'il est alors nécessaire de simuler l'écoulement sur des durées extrêmement longues. Cette thèse propose une approche différente, basée sur l'application d'algorithmes d'échantillonnage d'événements rares. L'objectif de ces algorithmes, issus de la physique statistique, est de modifier la statistique d'échantillonnage des trajectoires d'un système dynamique, de manière à favoriser l'occurrence d'événements rares. Si ces techniques ont été appliquées avec succès dans le cas de dynamiques relativement simples, l'intérêt de ces algorithmes n'est à ce jour pas clair pour des dynamiques déterministes extrêmement complexes, comme c'est le cas pour les écoulement turbulents.Cette thèse présente tout d'abord une étude de la dynamique et de la statistique associée aux fluctuations extrêmes de la force de traînée sur un obstacle carré fixe immergé dans un écoulement turbulent à deux dimensions. Ce cadre simplifié permet de simuler la dynamique sur des durées très longues, permettant d'échantillonner un grand nombre de fluctuations dont l'amplitude est assez élevée pour être qualifiée d'extrême.Dans un second temps, l'application de deux algorithmes d’échantillonnage est présentée et discutée.Dans un premier cas, il est illustré qu'une réduction significative du temps de calcul d'extrêmes peut être obtenue. En outre, des difficultés liées à la dynamique de l'écoulement sont mises en lumière, ouvrant la voie au développement de nouveaux algorithmes spécifiques aux écoulements turbulents. / This thesis discusses the numerical simulation of extreme fluctuations of the drag force acting on an object immersed in a turbulent medium.Because such fluctuations are rare events, they are particularly difficult to investigate by means of direct sampling. Indeed, such approach requires to simulate the dynamics over extremely long durations.In this work an alternative route is introduced, based on rare events algorithms.The underlying idea of such algorithms is to modify the sampling statistics so as to favour rare trajectories of the dynamical system of interest.These techniques recently led to impressive results for relatively simple dynamics. However, it is not clear yet if such algorithms are useful for complex deterministic dynamics, such as turbulent flows.This thesis focuses on the study of both the dynamics and statistics of extreme fluctuations of the drag experienced by a square cylinder mounted in a two-dimensional channel flow.This simple framework allows for very long simulations of the dynamics, thus leading to the sampling of a large number of events with an amplitude large enough so as they can be considered extreme.Subsequently, the application of two different rare events algorithms is presented and discussed.In the first case, a drastic reduction of the computational cost required to sample configurations resulting in extreme fluctuations is achieved.Furthermore, several difficulties related to the flow dynamics are highlighted, paving the way to novel approaches specifically designed to turbulent flows.

Surfactants are amphiphilic molecules which self-assemble at the interface in oil-water-surfactant mixtures such that the hydrophobic part, called tail, stays in oil and the remaining part, called head, resides in hydrophilic en-vironment. Depending upon concentration of individual components, these mixtures form several microphases, such as bilayers, micelles, columnar and lamellar phases. A lamellar phase, at equilibrium, is made up of alternat-ing layers of water and oil separated by surfactants, or of alternate layers of water and surfactant bilayers such that the hydrophilic heads are in contact with water. This equilibrium state is rarely achieved in macroscopic samples due to thermodynamic and kinetic constraints; instead, a lamellar fluid is usually disordered with a large number of defects. These defects have significant effect on the flow behaviour of the lamellar mesophase systems. They are known to alter the flow field, resulting stresses and in turn could get distorted or annihilated by the flow. In present work, we analyse this two way coupling between lamellar structure and flow field. The structural and rheological evolution of an initially disordered lamellar phase system under a shear flow is examined using a mesoscale model based on a free energy functional for the concentration field, which is the scaled difference in the concentration between the hydrophilic and hydrophobic components. Two distinct modes of structural evolution are observed depending only on Peclet number, which ratio of inertial forces to mass diffusivity, in-dependent of system size. At low Peclet number, local domains are formed which are then rotated and stretched by shear. A balance between defect creation and annihilation is reached due to which the system never reaches the equilibrium layer configuration. In the opposite limit, partially formed layers break and reform so as to form a nearly aligned lamellar phase con-figuration with residual defects. Viscosity of lamellar phase system increases with layer moduli, differences in viscosity of individual components, fluidity of the lamellae due to shear banding and defect pinning. These factors however, do not have any effect on alignment mechanism.

Turbulence is an open and challenging problem for mathematical approaches, physical modeling and numerical simulations. Numerical solutions contribute significantly to the understand of the nature and effects of turbulence. The focus of this thesis is the development of appropriate numerical methods for the computer simulation of turbulent flows. Many of the existing approaches to turbulence utilize analogies from kinetic theory. Degond & Lemou (J. Math. Fluid Mech., 4, 257-284, 2002) derived a k-✏ type turbulence model completely from kinetic theoretic framework. In the first part of this thesis, a numerical method is developed for the computer simulation based on this model. The Boltzmann equation used in the model has an isotropic, relaxation collision operator. The relaxation time in the collision operator depends on the microscopic turbulent energy, making it difficult to construct an efficient numerical scheme. In order to achieve the desired numerical efficiency, an appropriate change of frame is applied. This introduces a stiff relaxation source term in the equations and the concept of asymptotic preserving schemes is then applied to tackle the stiffness. Some simple numerical tests are introduced to validate the new scheme. In the second part of this thesis, alternative approaches are sought for more efficient numerical techniques. The Lattice Boltzmann Relaxation Scheme (LBRS) is a novel method developed recently by Rohan Deshmukh and S.V. Raghuram Rao for simulating compressible flows. Two different approaches for the construction of implicit sub grid scale -like models as Implicit Large Eddy Simulation (ILES) methods, based on LBRS, are proposed and are tested for Burgers turbulence, or Burgulence. The test cases are solved over a largely varying Reynolds number, demonstrating the efficiency of this new ILES-LBRS approach. In the third part of the thesis, as an approach towards the extension of ILES-LBRS to incompressible flows, an artificial compressibility model of LBRS is proposed. The modified framework, LBRS-ACM is then tested for standard viscous incompressible flow test cases.