L'ingénierie tissulaire intervient pour restaurer le tissu osseux. Parmi les traitements possibles, l'utilisation d'un bioréacteur à perfusion permet l'amplification in vitro de cellules souches ou osseuses prélevées chez le patient avant réimplantation. La contrainte de cisaillement générée par l'écoulement stimule mécaniquement les cellules et amplifie la production tissulaire. Cette technique souffre cependant de sa conception empirique et doit à présent être optimisée. L'objectif de cette thèse est l'étude et la modélisation de la croissance tissulaire et de la prolifération cellulaire à l'échelle du bioréacteur. En particulier, il s'agit de comprendre l'impact de l'écoulement sur la formation du tissu. Pour cela, une double approche de modélisation et d'expérimentation a été adoptée. Des expériences de culture cellulaire ont permis de mettre au jour la prolifération préférentielle des cellules près des parois du bioréacteur comme conséquence de l'hétérogénéité du support, et l'évolution de la morphologie du tissu. Un modèle prédisant le devenir des cellules ainsi que la croissance tissulaire à l'échelle du bioréacteur est proposé. L'aspect multi-échelles du problème est pris en considération et les procédures d'homogénéisation sont menées à bien grâce à la méthode de prise de moyenne volumique. / Bone tissue engineering aims at restoring bone tissues. Among the possible treatments, the use of a perfusion bioreactor allows the amplification in vitro of the patient bone or stem cells prior to implantation. The advantage of using such bioreactors is two-fold: in addition to greatly improving species transport, tissue production is enhanced. Although promising, this technique suffers from its empirical conception and now needs to be optimized. The purpose of this thesis is to study and model tissue growth and cell proliferation under a fluid flow of culture medium at the scale of the bioreactor. In particular, we wish to understand the impact of fluid flow on tissue formation. To this end, a double approach of experimentation and modeling has been adopted. Cell culture experiments in a perfusion bioreactor highlighted the preferential cell proliferation in the parietal region as a consequence of the heterogeneity of the scaffold, and the evolution of the tissue morphology. A model for predicting the cell's fate along with tissue growth at the scale of the bioreactor is proposed. The hierarchy of the system is considered and the upscaling procedures are carried out with the Volume Averaging Method.

Les conséquences liées à la présence de gradients thermiques sur le transfert de matière en milieu poreux sont encore aujourd’hui mal appréhendées, essentiellement en raison de la complexité induite par la présence de phénomènes couplés (thermodiffusion ou effet Soret). Le but de cette thèse est d’étudier et de comprendre l’influence que peut avoir un gradient thermique sur l’écoulement d’un mélange. L’objectif principal est de déterminer les coefficients effectifs modélisant les transferts de chaleur et de matière en milieux poreux, et en particulier le coefficient de thermodiffusion effectif. En utilisant la technique de changement d’échelle par prise de moyenne volumique nous avons développé un modèle macroscopique de dispersion incluant la thermodiffusion. Nous avons étudié en particulier l'influence du nombre de Péclet et de la conductivité thermique sur la thermodiffusion. Les résultats ont montré que pour de faibles nombres de Péclet, le nombre de Soret effectif en milieu poreux est le même que dans un milieu libre, et ne dépend pas du ratio de la conductivité thermique (solide/liquide). À l'inverse, en régime convectif, le nombre de Soret effectif diminue. Dans ce cas, un changement du ratio de conductivité changera le coefficient de thermodiffusion effectif. Les résultats théoriques ont montré également que, lors de la diffusion pure, même si la conductivité thermique effective dépend de la connectivité de la phase solide, le coefficient effectif de thermodiffusion est toujours constant et indépendant de la connectivité de la phase solide. Le modèle macroscopique obtenu par cette méthode est validé par comparaison avec des simulations numériques directes à l'échelle des pores. Un bon accord est observé entre les prédictions théoriques provenant de l'étude à l’échelle macroscopique et des simulations numériques au niveau de l’échelle de pores. Ceci démontre la validité du modèle théorique proposé. Pour vérifier et consolider ces résultats, un dispositif expérimental a été réalisé pour mesurer les coefficients de transfert en milieu libre et en milieu poreux. Dans cette partie, les nouveaux résultats expérimentaux sont obtenus avec un système du type « Two-Bulb apparatus ». La diffusion et la thermodiffusion des systèmes binaire hélium-azote et hélium-dioxide de carbone, à travers des échantillons cylindriques remplis de billes de différents diamètres et propriétés thermiques, sont mesurées à la pression atmosphérique. La porosité de chaque milieu a été déterminée par la construction d'une image 3D de l'échantillon par tomographie. Les concentrations sont déterminées par l'analyse en continu de la composition du mélange de gaz dans les ampoules à l’aide d’un catharomètre. La détermination des coefficients de diffusion et de thermodiffusion est réalisée par confrontation des relevés temporels des concentrations avec une solution analytique modélisant le transfert de matière entre deux ampoules. Les résultats sont en accord avec les résultats théoriques. Cela permet de conforter l’influence de la porosité des milieux poreux sur les mécanismes de diffusion et de thermodiffusion. / A multicomponent system, under nonisothermal condition, shows mass transfer with cross effects described by the thermodynamics of irreversible processes. The flow dynamics and convective patterns in mixtures are more complex than those of one-component fluids due to interplay between advection and mixing, solute diffusion, and thermal diffusion (or Soret effect). This can modify species concentrations of fluids crossing through a porous medium and leads to local accumulations. There are many important processes in nature and industry where thermal diffusion plays a crucial role. Thermal diffusion has various technical applications, such as isotope separation in liquid and gaseous mixtures, identification and separation of crude oil components, coating of metallic parts, etc. In porous media, the direct resolution of the convection-diffusion equations are practically impossible due to the complexity of the geometry; therefore the equations describing average concentrations, temperatures and velocities must be developed. They might be obtained using an up-scaling method, in which the complicated local situation (transport of energy by convection and diffusion at pore scale) is described at the macroscopic scale. At this level, heat and mass transfers can be characterized by effective tensors. The aim of this thesis is to study and understand the influence that can have a temperature gradient on the flow of a mixture. The main objective is to determine the effective coefficients modelling the heat and mass transfer in porous media, in particular the effective coefficient of thermodiffusion. To achieve this objective, we have used the volume averaging method to obtain the modelling equations that describes diffusion and thermodiffusion processes in a homogeneous porous medium. These results allow characterising the modifications induced by the thermodiffusion on mass transfer and the influence of the porous matrix properties on the thermodiffusion process. The obtained results show that the values of these coefficients in porous media are completely different from the one of the fluid mixture, and should be measured in realistic conditions, or evaluated with the theoretical technique developed in this study. Particularly, for low Péclet number (diffusive regime) the ratios of effective diffusion and thermodiffusion to their molecular coefficients are almost constant and equal to the inverse of the tortuosity coefficient of the porous matrix, while the effective thermal conductivity is varying by changing the solid conductivity. In the opposite, for high Péclet numbers (convective regime), the above mentioned ratios increase following a power law trend, and the effective thermodiffusion coefficient decreases. In this case, changing the solid thermal conductivity also changes the value of the effective thermodiffusion and thermal conductivity coefficients. Theoretical results showed also that, for pure diffusion, even if the effective thermal conductivity depends on the particle-particle contact, the effective thermal diffusion coefficient is always constant and independent of the connectivity of the solid phase. In order to validate the theory developed by the up-scaling technique, we have compared the results obtained from the homogenised model with a direct numerical simulation at the microscopic scale. These two problems have been solved using COMSOL Multiphysics, a commercial finite elements code. The results of comparison for different parameters show an excellent agreement between theoretical and numerical models. In all cases, the structure of the porous medium and the dynamics of the fluid have to be taken into account for the characterization of the mass transfer due to thermodiffusion. This is of great importance in the concentration evaluation in the porous medium, like in oil reservoirs, problems of pollution storages and soil pollution transport. Then to consolidate these theoretical results, new experimental results have been obtained with a two-bulb apparatus are presented. The diffusion and thermal diffusion of a helium-nitrogen and helium-carbon dioxide systems through cylindrical samples filled with spheres of different diameters and thermal properties have been measured at the atmospheric pressure. The porosity of each medium has been determined by construction of a 3D image of the sample made with an X-ray tomograph device. Concentrations are determined by a continuous analysing the gas mixture composition in the bulbs with a katharometer device. A transient-state method for coupled evaluation of thermal diffusion and Fick coefficients in two bulbs system has been proposed. The determination of diffusion and thermal diffusion coefficients is done by comparing the temporal experimental results with an analytical solution modelling the mass transfer between two bulbs. The results are in good agreement with theoretical results and emphasize the porosity of the medium influence on both diffusion and thermal diffusion process. The results also showed that the effective thermal diffusion coefficients are independent from thermal conductivity ratio and particle-particle touching.

Un dispositif expérimental a été développé pour l'étude du couplage de la croissance d'un biofilm de Shewanella oneidensis MR-1 et du transport conservatif de l'Erioglaucine. La croissance du biofilm a été suivie par mesure de conductivité hydraulique et par acquisition d'images à l'aide d'une caméra digitale. La fraction volumique du biofilm a été caractérisée par des essais d'élution d'une macromolécule (i.e. : le Bleu Dextran) par analogie avec les méthodes de chromatographie d'exclusion ou de filtration sur gel. Ainsi au bout de 29 jours, un biofilm quasi-homogène sur l'ensemble de la cellule d'écoulement (0,1×0,1×0,05m3) et équivalent à 50% du volume poral a été formé. L'influence de la croissance du biofilm sur les propriétés de transport du milieu a été évaluée. Les essais de transport conservatif de l'Erioglaucine effectués pour deux vitesses d'injection et à deux stades de croissance du biofilm (17 et 29 jours) ont montré l'influence d'hétérogénéités locales sur les paramètres de transport (i.e. : la porosité, la perméabilité et la dispersion hydrodynamique). Ainsi après 17 jours de culture quand le biofilm occupe partiellement le milieu poreux (moitié inférieure) un modèle à deux équations ou double milieu permet de caractériser le transport conservatif. A contrario après 29 jours de culture où le biofilm occupe tout le milieu poreux, un comportement fickien classique caractérise le transport. Les valeurs théoriques du coefficient de dispersion longitudinale prédites par la méthode de prise de moyenne volumique ont permis de reproduire de manière satisfaisante le comportement observé expérimentalement / An experimental device was performed for the study of coupling the growth of a Shewanella oneidensis MR-1 bacterial biofilm and the non reactive transport of Brilliant Blue FCF. The biofilm growth was monitored by hydraulic conductivity measurements and by image acquisition with a digital camera. The biofilm volume fraction was estimated through tracer experiments with a macromolecular tracer (i.e., Dextran Blue) as in size-exclusion chromatography or gel filtration chromatography. Then after 29 days of bacterial culture a quasi-homogenous biofilm was grown in the whole flow cell (0,1×0,1×0,05m3) occupying about 50% of void space volume. The influence of biofilm growth on porous media transport properties was evaluated. Conservative tracer experiments with Brilliant Blue FCF run at two hydrodynamic conditions and at two growth steps of biofilm (17 and 29 days) showed the influence of local heterogeneities on transport parameters (i.e., porosity, permeability and hydrodynamic dispersion). Then at 17 days of growth when the biofilm partially covers the porous medium (bottom half of the flow cell) a two-equation model or double-layer model was suitable to characterize the conservative transport. A contrario after 29 days of growth, when the biofilm covers the whole porous medium, a classical fickian model was convenient. Numerical values of longitudinal dispersion coefficient from volume averaging well fitted experimental results

Mes différentes activités de recherche qui s'insèrent dans une problématique axée autour du transport multi-échelle en milieu poreux sont détaillées ici au travers de 3 volets: (1) Influence des hétérogénéités sur les processus de transfert (2) Prise en compte des effets de couplage en transport réactif (3) Impact des processus biologiques sur l'hydrodynamique et le transport Un survol des défis qui restent à surmonter et des perspectives à venir de ce travail de recherche concluent le document.

La Tomographie à Emission de Positrons (TEP) permet d'obtenir une mesure dynamique et résolue en espace de la concentration d'un traceur radioactif injecté au patient. La quantification du débit sanguin cérébral par TEP repose sur l'utilisation d'un modèle cinétique le reliant à la variation spatio-temporelle de la concentration du traceur dans le cerveau. Différents modèles cinétiques sont proposés dans la littérature. Cependant, la majorité d'entre eux repose sur une modélisation compartimentale de l'organe observé. Dans ce cas, l'organe est subdivisé en un compartiment capillaire échangeant avec un compartiment tissulaire par une cinétique le plus souvent du premier ordre. Les résultats obtenus avec ce type de modèle sous-estiment le débit et ne permettent pas de prédire les premiers instants de la dynamique de répartition du traceur. Ces faiblesses ont été confirmées suite à l'amélioration de la résolution temporelle des tomographes, conduisant à l'élaboration de modèles incorporant plus de réalité physiologique. Cependant, tous ces modèles sont développés pour modéliser les échanges entre la micro-circulation et le tissu environnant à l'échelle d'un capillaire (échelle microscopique). Or la résolution spatiale des tomographes utilisés en clinique ne permet pas de distinguer la micro-circulation et le tissu. L'utilisation de ces modèles cinétiques avec des mesures de concentrations macroscopiques dépasse donc leur cadre théorique de validité et peut introduire des résultats faussés. Dans ce contexte, nous proposons un modèle cinétique basé sur le changement d'échelle (utilisant la méthode de prise de moyenne volumique). Ce changement d'échelle permet de remplacer l'ensemble micro-circulation/tissu par un volume fictif, homogène, dont les propriétés macroscopiques sont calculées à partir des propriétés microscopiques d'un Volume Elémentaire Représentatif (VER) du milieu. Dans un premier temps, afin de pouvoir comparer les résultats de ce modèle avec ceux du modèle compartimental standard, le VER considéré est constitué d'un capillaire unique et de son enveloppe de tissu, puis une complexité géométrique supplémentaire est introduite en considérant un réseau de capillaire isotrope à l'échelle de Darcy. Ces modèles sont utilisés pour identifier le débit à l'aide d'une méthode inverse. Pour cela, l'évolution temporelle du champ de concentration dans notre géométrie de référence, qui ne peut être mesurée par TEP en raison de sa faible résolution spatiale, est déterminée par des simulations numériques ainsi que par des mesures in vitro à l'aide d'un modèle expérimental, également développé au cours de ce travail, permettant de reproduire l'écoulement dans un canal traversant une matrice diffusante (gel d'alginate). / Positron Emission Tomography (PET) provides a dynamic and space-resolved measurement of the concentration field of a radioactive tracer previously injected to the patient. Quantification of cerebral blood flow by PET is based on the use of a kinetic model linking cerebral blood flow to the spatial and temporal variations of tracer concentration in the brain. Various kinetic models have been proposed in the literature. However, most of the mare based on a compartmental approach of the observed organ In this case, the organ is divided in two compartments, the capillary and the tissue, and the exchanges between these two compartments are often described by a first order kinetic model. Results obtained with this kind of model under estimate the flow rate and are notable to predict the first instants of the tracer dynamics distribution. With the continuous improvement of the temporal resolution of PET, these weaknesses have been confirmed, which led to the development of models incorporating more physiological reality. However, all these models have been developed to describe exchanges between micro-circulation and surrounding tissue at the scale of capillary vessels (microscopic scale). Because the spatial resolution of PET inclinical practice is insufficient to allow the distinction between micro-circulation and tissue, using of these models with kinetic measurement of macroscopic concentrations exceeds their theoretical validity and can introduce false results. In this context, we propose a kinetic model based on up-scaling (using the method of volume averaging). This up-scaling technique allows to replace the two previous compartments (tissue and micro-circulation) by an homogeneous fictive volume, whose macroscopic properties are calculated from the microscopic properties of are presentative elementary volume (REV) of the medium. First, in order to compare the results of this model with those of the standard compartmental model, the considered REV consists of a single capillary and its surrounding tissue. Second, additional geometric complexity is introduced by considering an isotropic capillary network at the Darcy scale. These models are used to identify the flow rate using an inverse method. For that purpose, the temporal evolution of concentration field in a geometry of reference, which can't be measured by PET due to its low spatial resolution, is determined by numerical simulations and by in vitro measurements. These measurements are performed using an experimental model developed during this work to reproduce the flow in a channel passing through a diffusive matrix (alginate gel).

Ces travaux de thèse visent à étudier les transferts couplés de chaleur, de masse et d’air au sein des matériaux poreux. Sur le volet de la modélisation, il s’agit de prédire le comportement hygrothermique de ces matériaux, à l’aide d’un modèle macroscopique, qui intègre à la fois l’effet du phénomène de thermodiffusion et celui de la pression totale de l’air s’exerçant sur les parois du bâtiment. Ce modèle, dont les paramètres d’entrée sont évalués expérimentalement, utilise des moteurs de transfert continus, d’où la possibilité de traiter des problèmes de transferts dans les matériaux multicouches. Il présente aussi l’avantage d’admettre, dans certaines configurations, des solutions analytiques d’où la possibilité d’entreprendre des comparaisons avec des solutions numériques. De plus, une justification formelle des équations de bilan de ce modèle a été abordée, moyennant l’utilisation d’une approche à changement d’échelle « micro-macro ». Il s’agit d’affiner la modélisation des transferts hydriques du comportement macroscopique, en utilisant des informations issues de la microstructure. Le passage de l’échelle microscopique à l’échelle macroscopique a été réalisé à l’aide de la méthode d’homogénéisation par prise de moyenne. Une des difficultés de l’utilisation de ce modèle réside dans l’identification des nombreux paramètres caractérisant les propriétés hygrothermiques des matériaux. Une partie du travail a été consacrée à l’évaluation des principales propriétés intrinsèques des matériaux moyennant l’élaboration de différents prototypes expérimentaux au laboratoire. Par ailleurs, une approche expérimentale dédiée à l’évaluation du processus de la thermodiffusion dans les matériaux poreux a été entreprise. Pour cela, une expérimentation relative à la détermination de l’effet du gradient de température et de la dynamique du processus d’échange d’eau à l’intérieur des parois a été mise en place au laboratoire. L’utilisation de la plateforme expérimentale MegaCup du Technical University of Denmark a permis de collecter des données relatives à la sensibilité de l’effet de la thermodiffusion sur les transferts couplés de chaleur, d’air et d’humidité. Une comparaison des résultats expérimentaux et numériques a ensuite été effectuée. Peu d’écarts ont été relevés. Aussi, une investigation expérimentale portant sur la contribution des infiltrations massiques sur les transferts hydriques dans les matériaux de construction a été réalisée. Moyennant le développement d’un banc d’essai, une caractérisation expérimentale du coefficient d’infiltration d’humidité a été entreprise. Ce dernier est utilisé comme paramètre d’entrée des modèles de simulation numérique. / The purpose of this thesis is to study coupled heat air and moisture transfer in porous building materials. Concerning the modeling part, the interest is to predict the hygrothermal behavior, with a macroscopic model, that incorporates simultaneously the effect of thermodiffusion phenomenon and that of total pressure on the building walls. The input parameters are evaluated experimentally using continuous driving potentials, where the ability to deal with problems of transfer in multilayer materials. In some configurations, it presents the advantage to undertake analytical solution that can be confronted with numerical solutions. Furthermore, a formal justification of balance equations of the developed model was addressed through the use of ascaling approach. Then, the modeling of macroscopic moisture transfer behavior, by implementing information from the microstructure can be refined. The transition of the microscopic to macroscopic scale was performed using the mean field homogenization. One of the difficulties with the use of this model lies in the identification of many parameters characterizing the hygrothermal properties of materials. Therefore, a part of the present work was devoted to the evaluation of the main properties of materials through the development of various experimental prototypes in the laboratory. More over, an experimental approach dedicated to the evaluation of the thermodiffusion process in porous materials has been undertaken. In this way, an experimentation concerning the determination of the temperature gradient and dynamics of water exchange process inside walls has been established. Furthermore, the use of the experimental platform MegaCup at theTechnical University of Denmark has collected data on the sensitivity of the thermodiffusion effect. Subsequently, a comparison of the experimental and the numerical results was performed. Few differences were observed. Otherwise, an experimental investigation on the contribution of the mass infiltration of water transfers in building materials was performed. A characterization of the moisture infiltration coefficient was performed through the development of the experimental test. This coefficient was used as an input parameter in the simulation models.

On étudie l'écoulement turbulent d'un fluide sur une canopée, que l'on modélise comme un milieu poreux déformable. Ce milieu poreux est en fait composé d'un tapis de fibres susceptibles de se courber sous la charge hydrodynamique du fluide, et ainsi de créer un couplage fluide-structure à l'échelle d'une hauteur de fibre (honami). L'objectif de la thèse est de développer un modèle macroscopique de cette interaction fluide-structure, afin d'en réaliser des simulations numériques. Une approche numérique de simulation aux grandes échelles est donc mise en place pour capturer les grandes structures de l'écoulement et leur couplage avec les déformations du milieu poreux. Pour cela nous dérivons les équations régissant la grande échelle, au point de vue du fluide ainsi que de la phase solide. À cause du caractère non-local de la phase solide, une approche hybride est proposée. La phase fluide est décrite d'un point de vue Eulerien, tandis que la description de la dynamique de la phase solide nécessite une représentation Lagrangienne. L'interface entre le fluide et le milieu poreux est traitée de manière continue. Cette approche de l'interface fluide/poreux est justifiée par un développement théorique sous forme de bilan de masse et de quantité de mouvement à l'interface. Ce modèle hybride est implémenté dans un solveur écrit en C$++$, à partir d'un solveur fluide disponible dans la librairie CFD \openfoam. Un préalable nécessaire à la réalisation d'un tel modèle macroscopique est la connaissance des phénomènes de la petite échelle en vue de les modéliser. Deux axes sont explorés concernant cet aspect. Le premier consiste à étudier les effets de l'inertie sur la perte de charge en milieu poreux. Un paramètre géométrique est proposé pour caractériser la sensibilité d'une microstructure poreuse à l'inertie de l'écoulement du fluide dans ses pores. L'efficacité de ce paramètre géométrique est validée sur une diversité de microstructures et le caractère général du paramètre est démontré. Une loi asymptotique est ensuite proposée pour modéliser les effets de l'inertie sur la perte de charge, et comprendre comment celle-ci évolue en fonction de la nature de la microstructure du milieu poreux. Le deuxième axe d'étude de la petite échelle consiste à étudier l'effet de l’interaction fluide-structure à l'échelle du pore sur la perte de charge au niveau macroscopique. Comme les cas présentent de grands déplacements de la phase solide, une approche par frontières immergées est proposée. Ainsi deux méthodes numériques sont employées pour appliquer la condition de non-glissement à l'interface fluid/solide: l'une par interface diffuse, l'autre par reconstitution de l'interface. Cela permet une validation croisée des résultats et d'atteindre des temps de calcul acceptables tout en maîtrisant la précision des résultats numériques. Cette étude permet de montrer que l'interaction fluide-structure à l'échelle du pore a un effet considérable sur la perte de charge effective au niveau macroscopique. Des questions fondamentales sont ensuite abordées, telles que la taille d'un élément représentatif ou la forme des équations de transport dans un milieu poreux souple. / We study the turbulent flow of a fluid over a canopy, that we model as a deformable porous medium. This porous medium is more precisely a carpet of fibres that bend under the hydrodynamic load, hence initiating a fluid-structure coupling at the scale of a fibre's height (honami). The objective of the thesis is to develop a macroscopic model of this fluid-structure interaction in order to perform numerical simulations of this process. The volume averaging method is implemented to describe the large scales of the flow and their interaction with the deformable porous medium. An hybrid approach is followed due to the non-local nature of the solid phase; While the large scales of the flow are described within an Eulerian frame by applying the method of volume averaging, a Lagrangian approach is proposed to describe the ensemble of fibres. The interface between the free-flow and the porous medium is handle with a One-Domain- Approach, which we justify with the theoretical development of a mass- and momentum- balance at the fluid/porous interface. This hybrid model is then implemented in a parallel code written in C$++$, based on a fluid- solver available from the \openfoam CFD toolbox. Some preliminary results show the ability of this approach to simulate a honami within a reasonable computational cost. Prior to implementing a macroscopic model, insight into the small-scale is required. Two specific aspects of the small-scale are therefore studied in details; The first development deals with the inertial deviation from Darcy's law. A geometrical parameter is proposed to describe the effect of inertia on Darcy's law, depending on the shape of the microstructure of the porous medium. This topological parameter is shown to efficiently characterize inertia effects on a diversity of tested microstructures. An asymptotic filtration law is then derived from the closure problem arising from the volume averaging method, proposing a new framework to understand the relationship between the effect of inertia on the macroscopic fluid-solid force and the topology of the microstructure of the porous medium. A second research axis is then investigated. As we deal with a deformable porous medium, we study the effect of the pore-scale fluid-structure interaction on the filtration law as the flow within the pores is unsteady, inducing time-dependent fluidstresses on the solid- phase. For that purpose, we implement pore-scale numerical simulations of unsteady flows within deformable pores, focusing for this preliminary study on a model porous medium. Owing to the large displacements of the solid phase, an immersed boundary approach is implemented. Two different numerical methods are compared to apply the no-slip condition at the fluid-solid interface: a diffuse interface approach and a sharp interface approach. The objective is to find the proper method to afford acceptable computational time and a good reliability of the results. The comparison allows a cross-validation of the numerical results, as the two methods compare well for our cases. This numerical campaign shows that the pore-scale deformation has a significant impact on the pressure drop at the macroscopic scale. Some fundamental issues are then discussed, such as the size of a representative computational domain or the form of macroscopic equations to describe the momentum transport within a soft deformable porous medium.