Modeling capillarity and two-phase flow in granular media: from porescale to network scale

Numerical simulations at the pore scale are a way to study the behavior of multiphase flows encountered in many natural processes and industrial applications. In this work, liquid morphology and capillary action are examined at the pore-scale by means of the multicomponent Shan-Chen lattice Boltzmann method (LBM). The accuracy of the numerical model is first contrasted with theoretical solutions. The numerical results are extended to complex microstructures beyond the pendular regime. The LBM has been employed to simulate multiphase flow through idealized granular porous media under quasi-static primary drainage conditions. LBM simulations provide an excellent description of the fluid-fluid interface displacement through the grains. Additionally, the receding phase trapped in the granular media in form of pendular bridges or liquid clusters is well captured. Unfortunately, such simulations require a significant computation time. A 2D model (Throat-Network model) based on analytical solutions is proposed to mimic the multiphase flow with very reduced computation cost, therefore, suitable to replace LBM simulations when the computation resources are limited. The approach emphasizes the importance of simulating at the throat scale rather than the pore body scale in order to obtain the local capillary pressure - liquid content relationships. The Throat-Network model is a starting point for a hybrid model proposed to solve 3D problems. The hybrid model combines the efficiency of the pore-network approach and the accuracy of the LBM at the pore scale to optimize the computational resources. The hybrid model is based on the decomposition of the granular assembly into small subsets, in which LBM simulations are performed to determine the main hydrostatic properties (entry capillary pressure, capillary pressure - liquid content relationship and liquid morphology for each pore throat). Despite the reduction of computation time, it is still not negligible and not affordable for large granular packings. Approximations by the Incircle and the MS-P method, which predict hydrostatic properties, are contrasted with the results provided by LBM and the hybrid model. Relatively accurate predictions are given by the approximations.Per tal d’estudiar els fluxos multifàsics presents a molts processos naturals i industrials és indispensable entendre les propietats físiques dels sistemes multifàsics a escala microscòpica. La morfologia dels fluids i les forces capil·lars s’investiguen a l’escala del porus mitjançant el ”multicomponent Shan-Chen lattice Boltzmann method (LBM)”. La precisió del model numèric ha estat contrastada amb solucions teòriques. Els resultats numèrics s’han estès a microestructures líquides complexes més enllà del règim pendular. El LBM ha estat emprat per simular fluxos multifàsics a través de medis porosos sota condicions quasi-estàtiques de drenatge. Les simulacions dutes a terme mitjançant el LBM proporcionen una descripció excel·lent del moviment de la interfície entre fluids a través de les partícules sòlides. Durant el drenatge, les simulacions numèriques són capaces de reproduir l’efecte del fluid atrapat dins el medi granular en forma de ponts o estructures líquides complexes. Malauradament, aquestes simulacions requereixen un temps de computació molt elevat. Per tal d’optimitzar els recursos de computació, proposem un model 2D (model Throat-Network) basat en solucions analítiques que permet reproduir fluxos multifàsics a través d’un conjunt de discs amb un temps de computació molt reduït. Per tant, aquest mètode és una alternativa que pot substituir les simulacions LBM quan els recursos de computació són escassos. El model Throat-Network destaca la importància de tractar el problema a l’escala de la gola del porus per tal d’obtenir les relacions pressió capil·lar - volum locals. Aquest enfocament és un punt de partida pel model híbrid que es presenta per resoldre els problemes en 3D. El model híbrid combina l’eficàcia del model ”Pore-Network” i la precisió del LBM a l’escala del porus. El model híbrid es basa en la descomposició d’una mostra granular en subdominis més petits, els quals corresponen a les goles dels porus (la gola dels porus és l’espai que connecta dos porus adjacents). Les simulacions LBM s’executen per a cada un dels subdominis per tal de determinar les propietats hidroestàtiques més rellevants (pressió capil·lar d’entrada, la corba de pressió capil·lar - grau de saturació i la morfologia líquida per cada una de les goles del porus). Malgrat la reducció significativa en el cost computacional del model híbrid, els temps de càlcul no són menyspreables i poc realistes per mostres granulars de grans dimensions. Les aproximacions donades pels mètodes de l’”Incircle” i el MS-P, que permeten estimar les propietats hidroestàtiques, han estat contrastades amb els resultats obtinguts amb LBM i el model híbrid.Les simulations numériques à l’échelle du pore sont fréquemment utilisées pour étudier le comportement des écoulements multiphasiques largement rencont des structures liquides et l’actiorés dans phénomènes naturels et applications industrielles. Dans ce travail, la morphologien capillaire sont examinées à l’échelle des pores par la méthode de Boltzmann sur réseau (LBM) à plusieurs composants selon le modèle de Shan-Chen. Les résultats numériques obtenus sont en bon accord avec les solutions théoriques. Les simulations numériques sont étendues à microstructures complexes au-delà du régime pendulaire. La LBM a été utilisée pour modéliser l’écoulement multiphasique à travers un milieu poreux idéalisé dans des conditions de drainage primaire quasi-statique. Les simulations LBM ont fourni une excellente description du déplacement de l’interface fluide-fluide à travers les grains. Pendant le drainage, les simulations LBM sont capables de reproduire la déconnexion d’une phase dans le milieu granulaire sous la forme de ponts pendulaires ou structures liquides complexes. Malheureusement, le temps de calcul nécessaire pour ce type de simulations est assez élevé. Afin d’optimiser les ressources de calcul, nous présentons un modèle 2D (modèle Throat-Network) basé sur des solutions analytiques pour décrire l’écoulement biphasique à travers un ensemble de disques dans un temps de calcul très réduit, donc le modèle 2D est susceptible de remplacer les simulations LBM lorsque les ressources de calcul sont limitées. L’approche souligne l’importance de simuler le problème a l’échelle de la gorge du pore pour obtenir les relations volume - pression capillaire locales. Le modèle Throat-Network est un point de départ pour le modèle hybride proposé pour résoudre les problèmes en 3D. Le modèle hybride combine l’efficacité de l’approche réseau de pores et la précision du LBM à l’échelle des pores. Le modèle hybride est basé sur la décomposition de l’échantillon en petits sous-domaines, dans lesquels des simulations LBM sont effectuées pour déterminer les propriétés hydrostatiques principales (pression capillaire d’entrée, courbe de drainage primaire et morphologie du liquide pour chaque gorge du pore). Malgré la réduction significative des temps de calcul obtenus avec le modèle hybride, le temps n’est pas négligeable et les modélisations numériques d’échantillons de grandes tailles ne sont pas réalistes. Les approximations données par les méthodes Incircle et MS-P, qui prédisent les propriétés hydrostatiques, sont comparées à celles de LBM et du modèle hybride

Localized fluidization in granular materials: Theoretical and numerical study

We present analytical and numerical results on localized fluidization within a granular layer subjected to a local injection of fluid. As the injection rate increases the three different regimes previously reported in the literature are recovered: homogeneous expansion of the bed, fluidized cavity in which fluidization starts developing above the injection area, and finally the chimney of fluidized grains when the fluidization zone reaches the free surface. The analytical approach is at the continuum scale, based on Darcy's law and Therzaghi's effective stress principle. It provides a good description of the phenomenon as long as the porosity of the granular assembly remains relatively homogeneous, i.e. for small injection rates. The numerical approach is at the particle scale based on the coupled DEM-PFV method. It tackles the more heterogeneous situations which occur at larger injection rates. The results from both methods are in qualitative agreement with data published independently. A more quantitative agreement is achieved by the numerical model. A direct link is evidenced between the occurrence of the different regimes of fluidization and the injection aperture. While narrow apertures let the three different regimes be distinguished clearly, larger apertures tend to produce a single homogeneous fluidization regime. In the former case, it is found that the transition between the cavity regime and the chimney regime for an increasing injection rate coincides with a peak in the evolution of inlet pressure. Finally, the occurrence of the different regimes is defined in terms of the normalized flux and aperture

Correlation functions of 2D Black Holes in JT gravity

Treballs Finals de Grau de Física, Facultat de Física, Universitat de Barcelona, Curs: 2023, Tutor: Cristiano GermaniThis work shows the possibility that the information loss paradox might be solved when the higher order of the steepest descent expansion in the gravity path integral is taken into account. The decaying behaviour of the correlation functions, which are the main signal of this information paradox, is replaced by a linear growing at large Lorentzian times, which is a signal of an information recovery. This has been shown in the 2D Jackiw-Teitelboim gravit

Modeling capillarity and two-phase flow in granular media: from porescale to network scale

Tesi en modalitat de cotutela: Universitat Politècnica de Catalunya i Université Grenoble Alpes. Resums extesos en francès i català a l'apendeix de la tesi.Numerical simulations at the pore scale are a way to study the behavior of multiphase flows encountered in many natural processes and industrial applications. In this work, liquid morphology and capillary action are examined at the pore-scale by means of the multicomponent Shan-Chen lattice Boltzmann method (LBM). The accuracy of the numerical model is first contrasted with theoretical solutions. The numerical results are extended to complex microstructures beyond the pendular regime. The LBM has been employed to simulate multiphase flow through idealized granular porous media under quasi-static primary drainage conditions. LBM simulations provide an excellent description of the fluid-fluid interface displacement through the grains. Additionally, the receding phase trapped in the granular media in form of pendular bridges or liquid clusters is well captured. Unfortunately, such simulations require a significant computation time. A 2D model (Throat-Network model) based on analytical solutions is proposed to mimic the multiphase flow with very reduced computation cost, therefore, suitable to replace LBM simulations when the computation resources are limited. The approach emphasizes the importance of simulating at the throat scale rather than the pore body scale in order to obtain the local capillary pressure - liquid content relationships. The Throat-Network model is a starting point for a hybrid model proposed to solve 3D problems. The hybrid model combines the efficiency of the pore-network approach and the accuracy of the LBM at the pore scale to optimize the computational resources. The hybrid model is based on the decomposition of the granular assembly into small subsets, in which LBM simulations are performed to determine the main hydrostatic properties (entry capillary pressure, capillary pressure - liquid content relationship and liquid morphology for each pore throat). Despite the reduction of computation time, it is still not negligible and not affordable for large granular packings. Approximations by the Incircle and the MS-P method, which predict hydrostatic properties, are contrasted with the results provided by LBM and the hybrid model. Relatively accurate predictions are given by the approximations.Per tal d’estudiar els fluxos multifàsics presents a molts processos naturals i industrials és indispensable entendre les propietats físiques dels sistemes multifàsics a escala microscòpica. La morfologia dels fluids i les forces capil·lars s’investiguen a l’escala del porus mitjançant el ”multicomponent Shan-Chen lattice Boltzmann method (LBM)”. La precisió del model numèric ha estat contrastada amb solucions teòriques. Els resultats numèrics s’han estès a microestructures líquides complexes més enllà del règim pendular. El LBM ha estat emprat per simular fluxos multifàsics a través de medis porosos sota condicions quasi-estàtiques de drenatge. Les simulacions dutes a terme mitjançant el LBM proporcionen una descripció excel·lent del moviment de la interfície entre fluids a través de les partícules sòlides. Durant el drenatge, les simulacions numèriques són capaces de reproduir l’efecte del fluid atrapat dins el medi granular en forma de ponts o estructures líquides complexes. Malauradament, aquestes simulacions requereixen un temps de computació molt elevat. Per tal d’optimitzar els recursos de computació, proposem un model 2D (model Throat-Network) basat en solucions analítiques que permet reproduir fluxos multifàsics a través d’un conjunt de discs amb un temps de computació molt reduït. Per tant, aquest mètode és una alternativa que pot substituir les simulacions LBM quan els recursos de computació són escassos. El model Throat-Network destaca la importància de tractar el problema a l’escala de la gola del porus per tal d’obtenir les relacions pressió capil·lar - volum locals. Aquest enfocament és un punt de partida pel model híbrid que es presenta per resoldre els problemes en 3D. El model híbrid combina l’eficàcia del model ”Pore-Network” i la precisió del LBM a l’escala del porus. El model híbrid es basa en la descomposició d’una mostra granular en subdominis més petits, els quals corresponen a les goles dels porus (la gola dels porus és l’espai que connecta dos porus adjacents). Les simulacions LBM s’executen per a cada un dels subdominis per tal de determinar les propietats hidroestàtiques més rellevants (pressió capil·lar d’entrada, la corba de pressió capil·lar - grau de saturació i la morfologia líquida per cada una de les goles del porus). Malgrat la reducció significativa en el cost computacional del model híbrid, els temps de càlcul no són menyspreables i poc realistes per mostres granulars de grans dimensions. Les aproximacions donades pels mètodes de l’”Incircle” i el MS-P, que permeten estimar les propietats hidroestàtiques, han estat contrastades amb els resultats obtinguts amb LBM i el model híbrid.Les simulations numériques à l’échelle du pore sont fréquemment utilisées pour étudier le comportement des écoulements multiphasiques largement rencont des structures liquides et l’actiorés dans phénomènes naturels et applications industrielles. Dans ce travail, la morphologien capillaire sont examinées à l’échelle des pores par la méthode de Boltzmann sur réseau (LBM) à plusieurs composants selon le modèle de Shan-Chen. Les résultats numériques obtenus sont en bon accord avec les solutions théoriques. Les simulations numériques sont étendues à microstructures complexes au-delà du régime pendulaire. La LBM a été utilisée pour modéliser l’écoulement multiphasique à travers un milieu poreux idéalisé dans des conditions de drainage primaire quasi-statique. Les simulations LBM ont fourni une excellente description du déplacement de l’interface fluide-fluide à travers les grains. Pendant le drainage, les simulations LBM sont capables de reproduire la déconnexion d’une phase dans le milieu granulaire sous la forme de ponts pendulaires ou structures liquides complexes. Malheureusement, le temps de calcul nécessaire pour ce type de simulations est assez élevé. Afin d’optimiser les ressources de calcul, nous présentons un modèle 2D (modèle Throat-Network) basé sur des solutions analytiques pour décrire l’écoulement biphasique à travers un ensemble de disques dans un temps de calcul très réduit, donc le modèle 2D est susceptible de remplacer les simulations LBM lorsque les ressources de calcul sont limitées. L’approche souligne l’importance de simuler le problème a l’échelle de la gorge du pore pour obtenir les relations volume - pression capillaire locales. Le modèle Throat-Network est un point de départ pour le modèle hybride proposé pour résoudre les problèmes en 3D. Le modèle hybride combine l’efficacité de l’approche réseau de pores et la précision du LBM à l’échelle des pores. Le modèle hybride est basé sur la décomposition de l’échantillon en petits sous-domaines, dans lesquels des simulations LBM sont effectuées pour déterminer les propriétés hydrostatiques principales (pression capillaire d’entrée, courbe de drainage primaire et morphologie du liquide pour chaque gorge du pore). Malgré la réduction significative des temps de calcul obtenus avec le modèle hybride, le temps n’est pas négligeable et les modélisations numériques d’échantillons de grandes tailles ne sont pas réalistes. Les approximations données par les méthodes Incircle et MS-P, qui prédisent les propriétés hydrostatiques, sont comparées à celles de LBM et du modèle hybride.Postprint (published version

Memòria Digital de Catalunya

Localització: Barcelona, Biblioteca de l'Ateneu Barcelonès, ms. 606Núm. 2 i nota: "Biblioteca teatral de Pau de Montellà: borrador", a la cobertaAutògrafEsmenes, algunes amb papers encartats, i cancel·lacionsAl f. 6v: PuigcerdàPaper pautatAcotacions de l'autor entre parèntesi

Memòria Digital de Catalunya

Localització: Barcelona, Biblioteca de l'Ateneu Barcelonès, ms. 610Núm. 3 i nota: "Biblioteca teatral de Pau Montellà: 2a copia", a la cobertaAutògrafEsmenes i cancel·lacionsAl f. 19r: BarcelonaPaper pautatAcotacions de l'autor entre parèntesi

Memòria Digital de Catalunya

Localització: Barcelona, Biblioteca de l'Ateneu Barcelonès, ms. 642Còpia de representacióA la coberta: "1a copia"Núm. 2 a la coberta, i nota: "Biblioteca teatral de Pau Montellà"Esmenes i cancel·lacionsAl f. 64v: PuigcerdàPaper pautatAcotacions de l'autor entre parèntesi

Memòria Digital de Catalunya

Localització: Barcelona, Biblioteca de l'Ateneu Barcelonès, ms. 640Còpia de representacióA la coberta, nota: "Biblioteca teatral de Pau Montellà"Esmenes i cancel·lacionsAl f. 1r, informació complementària cancel·lada: "(... ) y nostres detectives de antany o (...)"Al f. 70v: BarcelonaPaper pautatAcotacions de l'autor entre parèntesi

Memòria Digital de Catalunya

Localització: Barcelona, Biblioteca de l'Ateneu Barcelonès, ms. 605Núm. 1 i nota: "Biblioteca teatral de Pau Montellà: borrador", a la cobertaDedicatòria de l'autor a Alejandro SerraclaraAutògrafEsmenes, algunes amb papers encartats, i cancel·lacionsAl f. 16r: PuigcerdàPaper pautatAcotacions de l'autor entre parèntesi