En physique de l’atmosphère, la reconstruction d’une source polluante à partir des mesures de capteurs est une question importante. Elle permet en effet d’affiner les paramètres des modèles de dispersion servant à prévoir la propagation d’un panache de polluant, et donne aussi des informations aux primo-intervenants chargés d’assurer la sécurité des populations. Plusieurs méthodes existent pour estimer les paramètres de la source, mais leur application est coûteuse à cause de la complexité des modèles de dispersion. Toutefois, cette complexité est souvent nécessaire, surtout lorsqu’il s’agit de traiter des cas urbains où la présence d’obstacles et la météorologie instationnaire imposent un niveau de précision important. Il est aussi vital de tenir compte des différents facteurs d’incertitude, sur les observations et les estimations. Les travaux menés dans le cadre de cette thèse ont pour objectif de développer une méthodologie basée sur l’inférence bayésienne adaptative couplée aux méthodes de Monte Carlo pour résoudre le problème d’estimation du terme source. Pour cela, nous exposons d’abord le contexte scientifique du problème et établissons un état de l’art. Nous détaillons ensuite les formulations utilisées dans le cadre bayésien, plus particulièrement pour les algorithmes d’échantillonnage d’importance adaptatifs. Le troisième chapitre présente une application de l’algorithme AMIS dans un cadre expérimental, afin d’exposer la chaîne de calcul utilisée pour l’estimation de la source. Enfin, le quatrième chapitre se concentre sur une amélioration du traitement des calculs de dispersion, entraînant un gain important de temps de calcul à la fois en milieu rural et urbain. / In atmospheric physics, reconstructing a pollution source is a challenging but important question : it provides better input parameters to dispersion models, and gives useful information to first-responder teams in case of an accidental toxic release.Various methods already exist, but using them requires an important amount of computational resources, especially as the accuracy of the dispersion model increases. A minimal degree of precision for these models remains necessary, particularly in urban scenarios where the presence of obstacles and the unstationary meteorology have to be taken into account. One has also to account for all factors of uncertainty, from the observations and for the estimation. The topic of this thesis is the construction of a source term estimation method based on adaptive Bayesian inference and Monte Carlo methods. First, we describe the context of the problem and the existing methods. Next, we go into more details on the Bayesian formulation, focusing on adaptive importance sampling methods, especially on the AMIS algorithm. The third chapter presents an application of the AMIS to an experimental case study, and illustrates the mechanisms behind the estimation process that provides the source parameters’ posterior density. Finally, the fourth chapter underlines an improvement of how the dispersion computations can be processed, thus allowing a considerable gain in computation time, and giving room for using a more complex dispersion model on both rural and urban use cases.

Peut-on modéliser de manière fiable les effets collatéraux (en termes de quantité ou concentration de produits éjectés) et les conséquences d’une explosion en milieu clos, et en présence de produits chimiques liquides ? Pour répondre à cette vaste question, qui soulève spontanément de nombreux sous-problèmes, les travaux de thèse se sont déroulés en trois temps : 1/une étude qualitative et semi-quantitative du scénario général, afin de comprendre le déroulement chronologique des évènements, et d’émettre les hypothèses nécessaires à une modélisation, 2 /L’étude systématique des effets d’une explosion en milieu clos, en présence ou non de produits liquides. Avant de modéliser la dégradation du produit liquide soumis à une température et à une pression élevées, les expérimentations préliminaires ont en effet fait apparaître la nécessité de quantifier dans le temps et dans l’espace, les effets thermiques et mécaniques d’une explosion à volume constant, 3/L’élaboration de modèles 0D(t) à partir des conclusions précédentes afin, d’une part, d’estimer la durée de vie d’une phase liquide dans un environnement thermodynamique contraint et, d’autre part, de démontrer la possibilité de modéliser le problème global de manière réduite. En effet, le terme source d’un tel évènement ne peut être modélisé par une libération ponctuelle de produit : il s’agit d’une libération étendue dans le temps, par le biais d’un écoulement chaud a priori diphasique et de débit variable. Les couplages des phénomènes, observés expérimentalement, rendent nécessaires : 1 - Une modélisation instationnaire de l’évolution de la pression et une estimation du niveau de température atteint dans le volume d’étude, après détonation d’une charge explosive, 2 - Une modélisation de la libération de la phase liquide dans l’enceinte, sous forme de gouttes millimétriques ou de gouttelettes microniques 3 - Une modélisation instationnaire des transferts couplés de masse et d’énergie entre la phase liquide et la phase gazeuse en présence et prise en compte d’éventuelles réactions chimiques à haute pression et haute température 4 - Et enfin, une modélisation instationnaire des rejets à la brèche. L’étude d’une explosion à volume constant a montré qu’il est possible de modéliser de manière simple la montée continue en pression de l’enceinte par une fonction exponentielle croissante. Pour une configuration de référence donnée – explosion d’une sphère d’explosif dans un parallélépipède – la valeur maximale de pression est directement proportionnelle au taux de chargement en explosif, sur l’intervalle [0,01 – 0,6] kg/m3. Le passage à une géométrie différente ou plus complexe demande l’introduction d’un coefficient correctif pour traduire l’amplification (ou l’atténuation) de la combustion turbulente des produits de détonation avec l’air ambiant. En ce qui concerne le champ de température par contre, notre analyse a montré qu’il coexiste des zones chaudes et des zones dites « froides » et que la valeur de température homogène finale calculée à partir d’un code thermochimique ne peut constituer qu’une simple indication. Seule une estimation du volume respectif de ces zones a été proposée ici. Nous avons établi que les propriétés physico-chimiques des produits stockés sont un point clef du problème et on suppose ces données connues pour une gamme de produits chimiques liquides à pression ambiante, communément utilisés dans l’industrie. Seul le phénomène d’évaporation a été développé dans ce mémoire. L’introduction de réactions chimiques entre constituants se traduirait dans les modèles par des termes sources supplémentaires liés à l’apparition ou la disparition d’espèces. / Is it possible to model collateral effects due to an explosion (on a chemical facility for instance) occuring in a closed volume containing liquid chemical products storage units ?This thesis deals with a zerodimensionnal modelisation of such a 3D complex problem to asses the final thermodynamic state of chemical products released in the atmosphere. Developped sub-models take into account:- the unsteady time histories of the internal overpressure and temperature,- the unsteady liquid ejection (droplets sizes)- the unsteady modelisation of the local heat and mass transfers between the gas phase and the liquid phase- the unsteady ejection process of the resulting multiphase mixture in the environment.Models and sub-models are validated thanks to many experimental results.

Cette thèse est consacrée à l'analyse, à l'application et à l'extension bidimensionnelle, d'un nouveau schéma aux volumes finis (SRNH) proposé récemment pour une classe de système non homogène. L'analyse de stabilité du schéma, d'abord dans le cas scalaire ensuite dans le cas de systèmes, mène à une nouvelle formulation où intervient le signe de la matrice Jacobienne du système de lois de bilan considéré. Pour le système de Saint Venant avec terme de pente, on montre formellement que le schéma SRNHS vérifie la C-propriété exacte introduite pour les schémas équilibres par Bermùdez et Vázquez. Les résultats numériques 1D et 2D, en particulier du cas de rupture de barage sur un fond en forme de marche, montrent le degrés d'efficacité du schéma. Pour le système diphasiques des zones de non hyperbolicité peuvent exister, avec apparition de valeurs propres complexes dans la Jacobienne du système. On montre que pour les configurations faiblement non hyperboliques, on peut calculer le signe de la Jacobienne par l'algorithme de Newton-Schultz. Pour les configurations plus raides, où la méthode précédente ne fonctionne plus, on a recours à la méthode de perturbation par densité. Dans les deux cas évoqués, les tests numériques montrent que l'on approche la solution exacte du problème de Ransom avec une grande précision, et que l'on conserve la stabilité des calculs même avec un maillage de finesse relativement élevée.

Cette thèse est dédiée à l'étude de problèmes de couplage en espace entre différents modèles d'écoulements compressibles. Nous considérons des formulations monodimensionnelles où l'interface de couplage est mince, fixe et séparant deux régions de l'espace correspondant aux deux modèles à coupler. L'objectif de notre travail consiste à définir une condition de couplage à l'interface et à résoudre numériquement le problème de couplage muni de cette condition. Après un état de l'art non exhaustif sur le couplage de systèmes hyperboliques de lois de conservation, nous proposons une nouvelle formulation de condition de couplage basée sur l'ajout d'un terme source mesure agissant exactement sur l'interface de couplage. Nous supposons, dans un premier temps, que le poids associé à ce terme source est connu et constant. Deux solveurs de Riemann sont développés dont une approche par relaxation préservant les solutions équilibres du problème de couplage. Cette méthode par relaxation est reprise par la suite dans le cadre d'un problème d'optimisation sous contraintes pour déterminer un poids dynamique en temps selon différentes motivations de transmission à l'interface. Dans une seconde partie, nous développons un solveur de Riemann approché pour un modèle bifluide à deux pressions dans le cas d'un écoulement diphasique isentropique par phase. Le modèle en question a pour particularité de comprendre des termes non conservatifs que l'on réécrit alors sous la forme de termes sources mesures. L'approche par relaxation établie dans la partie précédente est alors étendue au cas du modèle bifluide, moyennant une estimation a priori des contributions non conservatives. Cette méthode nous permet, dans un dernier chapitre, de résoudre numériquement le problème de couplage interfacial entre un modèle bifluide à deux pressions et un modèle de drift-flux grâce à l'approche dite du modèle père.

Cette thèse est consacrée à la simulation numérique de problèmes d'écoulements de fluides raides régis par des systèmes de lois de bilan non homogènes, dans des configurations monodimensionnelles et bidimensionnelles. La méthode numérique utilisée est une extension d'un schéma à deux pas (SRNH), comportant un paramètre \alpha^n_(j+\frac(1)(2)) ajustable, proposé par le professeur F.Benkhaldoun dans un cadre monodimensionnel. Ainsi, en un premier temps on a introduit une variante SRNHR, obtenue en remplaçant la vitesse numérique (\frac(\Delta x)(\Delta t)) par la vitesse de Rusanov locale, en vue de l'extension du schéma au cas bidimensionnel. Par la suite, une analyse de stabilité du schéma, révèle que celui-ci peut être d'ordre 1 ou 2 selon la valeur du paramètre \alpha^n_(j+\frac(1)(2)). Une stratégie de variation de ce paramètre, basée sur la théorie des limiteurs a alors été adoptée. Le schéma peut ainsi être rendu d'ordre 1 dans les zones à forte variation de l'écoulement, et d'ordre 2, là où l'écoulement est régulier. Ensuite on a établi les conditions pour que ce schéma respecte la C-propriété exacte introduite par Bermùdez et Vazquez. Une étude d'implémentation des conditions aux limites, adaptée à ce schéma, a également été menée en se basant sur les invariants de Riemann. Dans la deuxième partie de la thèse, on a appliqué ce schéma à des systèmes monophasiques homogènes et non homogènes. Par exemple on a réalisé la simulation du problème de rupture de barrage sur une marche, pour des configurations 1D et 2D, en menant en particulier une étude de convergence numérique via la détermination des courbes d'erreurs. Enfin, on a utilisé le schéma pour la simulation numérique de systèmes diphasiques (Ransom 1D et 2D).

L'étalonnage des radiomètres embarqués sur satellite, utilisés pour l'observation de la Terre, est basé sur l'utilisation de sources de référence, capables de générer des températures de brillance ou des températures de bruit bien distinctes : une source chaude et une source froide. Dans le cas de la source froide, l'utilisation de sources externes (comme par exemple le ciel froid) est la méthode la plus communément utilisée. L'utilisation d'une charge froide active (ACL pour Active Cold Load) réalisée à base de transistors présente de nombreux avantages en termes d'intégration et de simplification du système. La première partie de cette thèse est consacrée au principe de la radiométrie. Après avoir rappelé l'importance des techniques de télédétection utilisées pour l'observation de la Terre, un état de l'art, non exhaustif, des différents types de radiomètre est présenté. L'un des éléments importants de l'instrument, quel que soit l'architecture choisie, est la phase de calibrage. De ce fait, un état de l'art sur les sources de référence, et plus particulièrement sur les charges froides actives est exposé. On constate notamment, que de part leur faible niveau de bruit et de leur couverture fréquentielle, les transistors à effet de champ (FET) sont de bons candidats. Nous présentons dans le second chapitre, un nouveau type de charge froide active basée sur l'utilisation d'un transistor bipolaire à hétérojonction (HBT) SiGe en bande L. Nous détaillons la technologie et la topologie du circuit et nous montrons les résultats obtenus concernant l'adaptation, la température de bruit et la sensibilité aux variations de température ambiante. En vue de comparer les différentes technologies en termes de stabilité à long terme, nous présentons également la conception et la caractérisation d'une ACL avec un PHEMT sur GaAs. La troisième partie concerne la réalisation d'un radiomètre spécifique, dans le but d'évaluer la stabilité à long terme des charges froides qui ont été déve loppées. La conception de ce radiomètre a été articulée autour de deux points critiques : limiter les variations de gain du système pour accroitre sa stabilité et obtenir une très bonne sensibilité. Pour cela, les choix technologiques et techniques ainsi que les performances du système sont détaillées. Enfin, la dernière partie de cette thèse est consacrée à l'étude de stabilité à long terme des ACLs qui ont été réalisées. Celles-ci sont très stables sur le court-terme, et présentent des caractéristiques intéressantes sur le long terme. Les différences concernant les comportements et les performances des deux circuits sont abordées.

Les écoulements de convection naturelle en espace confiné sont généralement turbulents (grandes dimensions et/ou écarts de température importants). L'approfondissement de la connaissance de ce régime et des transferts qui y sont associés semble donc essentielle. Cette étude expérimentale concerne les mécanismes de transferts thermiques turbulents en s'appuyant sur des mesures couplées vitesse/température dans un écoulement de convection naturelle turbulente à haut nombre de Rayleigh (Ra1-1= 1,2x1011) au sein d'une cavité différentiellement chauffée de rapport de forme vertical (hauteur/largeur) égale à 4. Une chaîne de mesure PIV permet les acquisitions de vitesse tandis que la température est mesurée par micro-thermocouple de type K (0=12,7 μm). Les deux mesures étant synchronisées via un générateur de pulses type BNC®. Une attention particulière a été portée à la détermination des conditions aux limites en température des parois adiabatiques et au post traitement des données de PIV à l'aide d'une décomposition orthogonale basée sur l'énergie cinétique des champs de vitesse (POD). Cela a rendu possible, d'une part le calcul expérimental de quantités liées à la turbulence comme les flux d'enthalpie et des nombres adimensionnés pour la turbulence (Prandtl, diffusivité ... ) et, d 'autre part, l'évaluation de la contribution des autres termes des équations de Navier-Stokes (gradient de pression et terme source volumique). Enfin, une discussion sur les échelles caractéristiques de l'écoulement est menée et une comparaison avec des résultats de simulations numériques est apportée. / Natural convection flows in confined spaces are often turbulent (large dimensions and/or temperature difference ). Improving the knowledge of this type of regime and of the associated heat and mass transfers seems, therefore, to be essential. This experirnental work studies heat transfer mechanisms using coupled and synchronized measurernents of temperature and velocity in a turbulent natural convection flow at high Rayleigh number (Ra1-1=I.2 x1011) in a differentially heated cavity with a vertical aspect ratio (height/width) of 4. Velocity measurernents are acquired by PIV and temperature is measured with K-type micro-thermocouple (0=12.7 μm). Both measurementsystems are synchronized using a BNC® pulse generator. Thermal boundary conditions of the adiabatic walls and post treatment of PIV data using an orthogonal decomposition based on the kinetic energy of the flow - POD - have been considered with caution. It makes possible, on the one hand, the experimental assessment of turbulent quantities like enthalpy fluxes and dimensionless nurnber (Prandtl, diffusivity ... ) and, on the other hand the estimation of the contribution of the other terms of Navier-Stokes equations (pressure gradient and source term). Finally a discussion on characteristic scales of the flow and a comparison with numerical simulations are provided.

La simulation numérique adaptative d'écoulements présentant des phénomènes multi-échelles s'effectue généralement au moyen d'une hiérarchie de modèles différents selon l'échelle mise en jeu et le niveau de précision requis. Ce type de modélisation numérique entraîne des problèmes de couplage multi-échelles complexes. Cette thèse est ainsi dédiée au développement, à l'analyse et à la mise en œuvre de méthodes performantes permettant de résoudre des problèmes de couplage en espace de modèles décrits par des systèmes d'équations aux dérivées partielles hyperboliques.Dans une première partie, nous développons et analysons une méthode numérique dédiée au couplage interfacial des équations d'Euler mono-dimensionnelles. Chacun des systèmes de lois de conservation est muni d'une loi de pression distincte et l'interface de couplage séparant ces modèles est supposée fixe et infiniment mince. Les conditions de transmission sont modélisées par un terme source mesure localisé à l'interface de couplage. Le poids associé à cette mesure modélise les pertes de conservation à l'interface (typiquement des pertes de charge) et sa définition permet l'application de plusieurs stratégies de couplage. Notre méthode d'approximation repose sur les techniques d'approximation par relaxation de type Suliciu. La résolution exacte du problème de Riemann pour le système relaxé nous permet de définir un schéma numérique équilibre pour le modèle de couplage. Ce schéma préserve certaines solutions stationnaires du modèle de couplage et est applicable pour des lois de pression générales. L'implémentation de notre méthode permet de mener des expériences numériques illustrant les propriétés de notre schéma. Par exemple, nous montrons qu'il est possible de contrôler l'écoulement à l'interface de couplage en calculant des poids solutions de problèmes d'optimisation sous contraintes.La deuxième partie de cette thèse est dédiée au développement de deux schémas numériques d'ordre arbitrairement élevé en espace pour l'approximation des solutions stationnaires du problème mixte associé au modèle de Jin et Xin. Nos schémas d'approximation reposent sur la méthode de Galerkin discontinue. L’approximation des solutions du problème mixte par notre premier schéma fait intervenir uniquement des erreurs de discrétisation tandis que notre deuxième schéma est constitué à la fois d'erreurs de modélisation et de discrétisation. L'erreur de modélisation provient du remplacement, dans certaines régions spatiales, de la résolution du modèle de relaxation par celle de l'équation scalaire équilibre associée. Sous l'hypothèse d'une interface de couplage éventuellement caractéristique, la résolution du problème de Riemann associé au modèle couplé nous permet de construire un schéma numérique d'ordre arbitrairement élevé prenant en compte l'éventuelle existence de couches limites à l'interface de couplage. Enfin, la mise en œuvre de ces méthodes nous permet d'analyser quantitativement et qualitativement les erreurs de modélisation et de discrétisation commises lors de l'utilisation du schéma couplé. Ces erreurs sont fonction du niveau de raffinement de maillage utilisé, du degré de polynôme choisi et de la position de l'interface de couplage. / The adaptive numerical simulation of multiscale flows is generally carried out by means of a hierarchy of different models according to the specific scale into play and the level of precision required. This kind of numerical modeling involves complex multiscale coupling problems. This thesis is thus devoted to the development, analysis and implementation of efficient methods for solving coupling problems involving hyperbolic models.In a first part, we develop and analyze a coupling algorithm for one-dimensional Euler systems. Each system of conservation laws is closed with a different pressure law and the coupling interface separating these models is assumed fix and thin. The transmission conditions linking the systems are modelled thanks to a measure source term concentrated at the coupling interface. The weight associated to this measure models the losses of conservation and its definition allows the application of several coupling strategies. Our method is based on Suliciu's relaxation approach. The exact resolution of the Riemann problem associated to the relaxed system allows us to design an extremely accurate scheme for the coupling model. This scheme preserves equilibrium solutions of the coupled problem and can be used for general pressure laws. Several numerical experiments assess the performances of our scheme. For instance, we show that it is possible to control the flow at the coupling interface when solving constrained optimization problems for the weights.In the second part of this manuscript we design two high order numerical schemes based on the discontinuous Galerkin method for the approximation of the initial-boundary value problem associated to Jin and Xin's model. Our first scheme involves only discretization errors whereas the second approximation involves both modeling and discretization errors. Indeed in the second approximation, we replace in some regions the resolution of the relaxation model by the resolution of its associated scalar equilibrium equation. Under the assumption of a possible characteristic coupling interface, we exactly solve the Riemann problem associated to the coupled model. This resolution allows us to design a high order numerical scheme which captures the possible boundary layers at the coupling interface. Finally, the implementation of our methods enables us to analyze quantitatively and qualitatively the modeling and discretization errors involved in the coupled scheme. These errors are functions of the mesh size, the degree of the polynomial approximation and the position of the coupling interface.

L'objectif de ce travail est le développement d'un modèle et d'une méthode numérique adaptés à la simulation du ruissellement d'eau de pluie sur des surfaces agricoles. Pour cela, nous utilisons un système d'équations aux dérivées partielles pour les eaux peu profondes : le système de Saint Venant. La pluie et l'infiltration y sont pris en compte par l'ajout de termes source et le couplage avec un modèle d'infiltration. Le système étant hyperbolique et conservatif nous choisissons d'utiliser un schéma aux volumes finis avec reconstruction hydrostatique. Avec cette reconstruction nous obtenons un schéma permettant de traiter les équilibres stationnaires et les interfaces sec/mouillé dues aux événements pluvieux. Tout d'abord, nous effectuons une comparaison systématique de différents flux numériques, de différentes reconstructions d'ordre deux (MUSCL et ENO modifiées) et de différentes méthodes de traitement des frottements sur des solutions analytiques unidimensionnelles. Ces comparaisons nous permettent de choisir une méthode numérique adaptée à la simulation du ruissellement qui en outre capture les phénomènes de type roll-waves. Ensuite nous nous intéressons au traitement de la pluie et de l'infiltration à l'aide du modèle de Green-Ampt. Cette méthode généralisée en dimension deux est validée sur des résultats expérimentaux obtenus à l'INRA d'Orléans et à l'IRD. La méthode numérique choisie et validée a été implémentée dans FullSWOF_2D un logiciel libre écrit en C++ pour la simulation de ruissellement d'eau de pluie sur des surfaces agricoles.

Nous analysons dans cette thèse différents aspects associés à la modélisation des écoulements à surface libre en eaux peu profondes (Shallow Water). Nous étudions tout d’abord le système d’équations de Saint-Venant à deux dimensions et leur résolution par la méthode numérique des volumes finis, en portant une attention particulière sur les aspects hyperboliques et conservatifs. Ces schémas permettent de traiter les équilibres stationnaires, les interfaces sec/mouillé et aussi de modéliser des écoulements subcritique, transcritique et supercritique. Nous présentons ensuite la théorie de la méthode d’assimilation variationnelle de données adaptée à ce type d’écoulement. Son application au travers des études de sensibilité est longuement discutée dans le cadre de l'hydraulique à surface libre. Après cette partie à caractère théorique, la partie tests commence par une qualification de l’ensemble des méthodes numériques qui sont implémentées dans le code DassFlow, développé à l’Université de Toulouse, principalement à l’IMT mais aussi à l’IMFT. Ce code résout les équations Shallow Water par une méthode de volumes finis et est validé par comparaison avec les solutions analytiques pour des cas tests classiques. Ces mêmes résultats sont comparés avec un autre code d’hydraulique à surface libre aux éléments finis en deux dimensions, Telemac 2D. Une particularité notable du code DassFlow est de permettre l’assimilation variationnelle de données grâce au code adjoint permettant le calcul du gradient de la fonction coût. Ce code adjoint a été obtenu en utilisant l'outil de différentiation automatique Tapenade (Inria). Nous testons ensuite sur un cas réel, hydrauliquement complexe, différentes qualités de Modèles Numériques de Terrain (MNT) et de bathymétrie du lit d’une rivière. Ces informations proviennent soit d’une base de données classique type IGN, soit d’informations LIDAR à très haute résolution. La comparaison des influences respectives de la bathymétrie, du maillage et du type de code utilisé, sur la dynamique d’inondation est menée très finement. Enfin nous réalisons des études cartographiques de sensibilité aux paramètres du modèle sur DassFlow. Ces cartes montrent l’influence respective des différents paramètres ou de la localisation des points de mesure virtuels. Cette localisation optimale de ces points est nécessaire pour une future assimilation de données efficiente. / We analyze in this thesis various aspects associated with the modeling of free surface flows in shallow water approximation. We first study the system of Saint-Venant equations in two dimensions and its resolution with the numerical finite volumes method, focusing in particular on aspects hyperbolic and conservative. These schemes can process stationary equilibria, wetdry interfaces and model subcritical, transcritical and supercritical flows. After, we present the variational data assimilation method theory fitted to this kind of flow. Its application through sensitivity studies is fully discussed in the context of free surface water. After this theoretical part, we test the qualification of numerical methods implemented in the code Dassflow, developed at the University of Toulouse, mainly at l'IMT, but also at IMFT. This code solves the Shallow Water equations by finite volume method and is validated by comparison with analytical solutions for standard test cases. These results are compared with another hydraulic free surface flow code using finite elements in two dimensions: Telemac2D. A significant feature of the Dassflow code is to allow variational data assimilation using the adjoint method for calculating the cost function gradient. The adjoint code was obtained using the automatic differentiation tool Tapenade (INRIA). Then, the test is carried on a real hydraulically complex case using different qualities of Digital Elevation Models (DEM) and bathymetry of the river bed. This information are provided by either a conventional database types IGN or a very high resolution LIDAR information. The comparison of the respective influences of bathymetry, mesh size, kind of code used on the dynamics of flooding is very finely explored. Finally we perform sensitivity mapping studies on parameters of the Dassflow model. These maps show the respective influence of different parameters and of the location of virtual measurement points. This optimal location of these points is necessary for an efficient data assimilation in the future.