Le problème d'attaque des bétons par sulfates demeure d'actualité dans le domaine du stockage des déchets radioactifs, parce que certains déchets contiennent des sulfates. Le présent rapport contribue à l'étude des mécanismes de dégradation des enrobés de déchets radioactifs par sulfate. Les systèmes expansifs CaO-Al2O3-SO3-H2O sans et avec alcalins ont été étudiés. Dans le système sans alcalins, c'est la formation de l'ettringite secondaire par apport d'eau externe qui conduit à l'expansion importante par effet stérique. Dans le système avec alcalins, les comportements d'expansion sont complètement différents. Un nouveau minéral baptisé "phase U", de type AFm incorporant du sodium, a été mis en évidence. Il semble être le responsable de l'expansion observée. Les conditions de formation, le produit de solubilité, les méthodes de synthèse et les distances interréticulaires de la phase U ont été discutés. Les différentes configurations des enrobés réels ont été ensuite étudiées, le comportement étant lié notamment à la phase U, très peu étudiée jusqu'à présent en milieu de ciment. Nous avons retenu comme mécanismes expansifs possibles soit la formation de la phase U secondaire, soit la transformation de la phase U en ettringite, soit la transformation de la thénardite en mirabilite. La première hypothèse a été validée sur un système simplifié, la modélisation de toutes les configurations des enrobés par les bilans massique et volumique nous permet de déterminer théoriquement la répartition massique des hydrates et d'estimer l'expansion par formation de la phase U secondaire. La deuxième hypothèse a été aussi parfaitement vérifiée : l'observation de l'expansion et l'analyse minéralogique des échantillons contenant la phase U et subissant la lixiviation dans différentes solutions nous permettent de supposer que l'expansion est provoquée par la répulsion électrique interparticulaire d'ettringite colloïdale, issue de la transformation de la phase U initiale et formée en présence de chaux. La cinétique de transformation de la phase U en ettringite a été suivie expérimentalement et analysée avec la théorie de la diffusion. La troisième hypothèse n'a pas été validée. En conclusion, la phase U peut provoquer une expansion par deux principaux mécanismes se manifestant individuellement ou simultanément : la formation de la phase U secondaire et la transformation de îa phase U en ettringite colloïdale en présence de chaux. Le présent travail constitue une première étape importante et nécessaire pour prédire la durabilité d'une structure du stockage de déchets radioactifs sulfatés.

Pour la gestion en temps réel des étiages comme pour celle des crues, l'une des principales missions des gestionnaires de bassin consiste à prévenir les situations de crise, généralement définies par des franchissements de débits seuils en différents points stratégiques, afin d'en limiter les dégâts, tant économiques et écologiques que matériels et humains. La thèse propose une méthodologie de synthèse et d'identification d'un modèle de simulation des débits, intégrant les phénomènes de transferts hydrauliques (propagation en cours d'eau) et hydrologiques (transformation de la pluie en débit), ainsi que les transferts relatifs aux prélèvements dans le cas des étiages. Le modèle proposé, construit dans le souci constant des contraintes de l'opérationnel telles que la robustesse face aux aléas du terrain, l'assimilation des données en temps réel ou la commande automatique des vannes de barrage, présente un module hydraulique à base physique simplifiée couplé, selon une approche semi-distribuée, à un module hydrologique global conceptuel. Cette démarche permet la prise en compte de débits latéraux non jaugés dans le transfert en cours d'eau. L'approche est validée à travers des exemples théoriques et divers cas d'application réels : transfert en rivière soumis à l'influence de la marée, synthèse d'un contrôleur boucle ouverte sur un canal d'irrigation, simulation de débits à l'exutoire de divers bassins versants français, prévision des crues, synthèse d'observateurs d'états pour la reconstitution de prélèvements.

La coulée continue est une étape importante dans la chaîne de fabrication de l'acier. Ce procédé agit comme un échangeur thermique qui permet de solidifier l'acier en fusion. De nombreuses études ont révélé qu'une partie importante des défauts de surface du produit solidifié est provoquée par les perturbations du niveau d'acier liquide en lingotière. Chaque ligne de production dispose d'un système de régulation de niveau dont l'objectif est de maintenir un niveau stable d'acier liquide. En présence de certaines perturbations ou lors des phases transitoires, les performances de ces systèmes de régulation se dégradent fortement. L'objectif principal de la thèse est de développer des solutions de contrôle robuste de niveau afin de rejeter les perturbations internes et externes à la boucle de régulation quelles que soient les conditions de coulée. Un outil d'estimation a été développé pour diagnostiquer et rejeter certaines perturbations. L'activité de recherche a porté également sur le recensement et l'évaluation des lois de commande utilisées actuellement dans les usines ou proposées par les fournisseurs dans le but d'établir un référentiel de performances que nous avons ensuite amélioré en mettant en oeuvre de nouvelles stratégies de régulation. Enfin, une commande adaptative multi-modèles combinant à la fois les fonctions de régulation robuste de niveau et d'identification des paramètres a été implantée. Une maquette hydraulique a été utilisée à plusieurs reprises pour reproduire le fonctionnement de la coulée continue sur site afin d'évaluer le gain en performance de la solution retenue. Les solutions proposées présentent un caractère générique et peuvent s'appliquer à d'autres procédés sidérurgiques avec seulement quelques ajustements.

La chronologie des recherches théoriques et numériques exposées dans cet ouvrage est étroitement liée à notre attachement à lexpérimental et aux préoccupations pratiques des grands projets de constructions hydrauliques. Les applications variées qui sont envisagées quotidiennement au HACH (Hydrodynamique Appliquée et Constructions Hydrauliques Secteur Mécanique des Solides, des Fluides et des Structures - Département ArGEnCo Université de Liège), dont plusieurs exemples sont présentés, ont permis didentifier aussi bien des pistes damélioration des modèles numériques existants au sein de la suite WOLF(suite logicielle intégrée développée par le HACH pour la modélisation et loptimisation des écoulements à surface libre), que des besoins en nouvelles potentialités de représentation. Cest de là que sont nés la plupart des développements détaillés ici, avec comme objectif principal laugmentation des capacités des logiciels du HACH pour laide à la gestion et au design douvrages hydrauliques. Le dimensionnement de structures hydrauliques, tout comme la majeure partie des problèmes touchant à lhydrodynamique des écoulements à surface libre rencontrés dans la pratique, nécessite des études à grande et petite échelles, à laide de modèles adaptés localement aux phénomènes représentés, et donc variables au sein dune même étude selon le point particulier examiné. Dautre part, loptimisation de la forme ou du fonctionnement douvrages hydrauliques requiert des outils de calibration robustes permettant une approche objective. Ces derniers doivent de plus être couplés à des logiciels de simulation dont la représentativité et les temps de calcul sont en accord avec les exigences de loptimisation. Ces deux constats principaux posent le cadre des recherches entreprises et sont à la base de la définition des objectifs du travail. Les premiers chapitres de développements sattachent à décrire le cadre numérique de la suite logicielle WOLF, et plus particulièrement des codes 1D et 2D. La rationalisation des temps de calcul des simulations filaires à grande échelle est améliorée à ce stade par limplémentation au sein du logiciel unidimensionnel dun schéma dintégration temporelle implicite, en parallèle du schéma explicite existant. De même, certaines améliorations sont déjà apportées au logiciel WOLF2D. Elles concernent lextension automatique du domaine de calcul en fonction des mailles immergées émergées, le traitement des termes sources de pente de fond et de frottement en accord avec une approche énergétique, la conservation volumique du modèle ainsi que la réorganisation de la filière de résolution des équations de façon à permettre une totale latitude dans le choix du modèle de calcul utilisé. Lenrichissement des potentialités de modélisation des outils de la suite WOLF est ensuite poursuivi par létablissement, limplémentation et lapplication à des cas théoriques et pratiques, dun modèle de turbulence adapté aux caractéristiques spécifiques de la simulation des écoulements intégrée sur la hauteur deau. Afin de dépasser les limitations, tant informatiques que de discrétisation, inhérentes à lutilisation dune taille de maille unique sur un domaine de simulation bidimensionnel, des développements sont entrepris pour permettre de mener les calculs bidimensionnels sur un maillage structuré cartésien constitué de plusieurs zones à la taille de maille différente. Dans le même ordre didée, des développements sont réalisés pour permettre lutilisation de lensemble des outils de modélisation disponibles au sein du logiciel WOLF2D de façon simultanée et automatique au sein dune même simulation. Létape ultime de ces développements est réalisée en couplant les modèles uni- et bidimensionnel au sein dune simulation unique. Une fois toutes ces étapes de développement des modèles hydrodynamiques franchies, un logiciel adapté de calage automatique de paramètres est mis au point. Utilisable avec nimporte quel code de la suite WOLF, il permet la calibration de lensemble des paramètres physiques présents au sein des logiciels hydrodynamiques, ainsi que la résolution de tout problème doptimisation ou de calibration à laide des mêmes modèles. La parallélisation de ce logiciel basé sur la méthode des Algorithmes Génétiques ouvre la porte à des applications denvergure, dont le dimensionnement hydraulique douvrages sur base de simulations bidimensionnelles découlements turbulents à surface libre sur maillage multibloc. A chaque étape, la validation des développements est réalisée de façon rigoureuse par la confrontation des résultats numériques à ceux de cas test, théoriques et expérimentaux, originaux ou proposés dans la littérature. De nombreux exemples dapplications pratiques illustrent l'ensemble des potentialités de modélisation développées et prouvent leur applicabilité et leur efficacité face aux problèmes pratiques rencontrés par les ingénieurs hydrauliciens./ The chronology of the theoretical and numerical researches outlined in this text is closely linked to our attachment to the physical modelling and to the practical concerns related to important hydraulic projects. The various applications carried out daily in the HACH (Applied Hydrodynamics and Hydraulic Constructions - Solid, Fluid and Structures Mechanics ArGEnCo Department University of Liege), several of the being presented here, enabled to identify the ways to improve the existing numerical models of WOLF (integrated software package developed by the HACH for free surface flows modelling), as well as needs in new representation potentialities. The most of the developments described in this text were born from there, with as main objective to increase the capacities of the HACHs models to help to the management and the design of hydraulic constructions. The hydraulic structures design, and the main part of the practical problems related to free surface flows hydrodynamics, requires both large and small scales studies, with reliable models locally suited to the represented phenomena, and thus variable in a single study depending on the particular point examined. On another hand, the shape and working optimization of hydraulic structures requires robust calibration tools allowing an objective approach. The latter have to be linked to simulation software with modelling capacities and computation times in agreement with the optimization requirements. Both these assessments specify the framework of the researches of this thesis and are at the roots of the definition of the objectives of this work. The first chapters of developments aim to describe the numerical framework of the WOLF package, and more particularly the 1D and 2D models. The rationalization of the computation times for large scale one-dimensional simulations is improved at this stage by implementing in the 1D solver an implicit time integration scheme in parallel to the existing explicit one. Similarly, some improvements have been brought to the solver WOLF2D. They concern the automatic extension of the computation domain regarding the wet and dry cells, the treatment of the slope and roughness source terms with an energetic approach, the water volume conservation as well as the reorganization of the equations resolution scheme in order to be totally free in the choice of the mathematical model to solve. The enrichment of the modelling potentialities of the WOLF software is then pursued by setting, implementing and applying an original turbulence model suited to the specific characteristics of the depth integrated flow modelling. In order to get over the computer and modelling limitations inherent in the use of a single mesh size to model a 2D domain, developments have been carried out to allow the realisation of two-dimensional computations on a structured Cartesian grid composed of several areas with different mesh sizes. In the same way, developments have been carried out to use the whole of the mathematical models available in WOLF2D simultaneously and automatically in a single simulation. The ultimate step of these developments consisted in linking the 1D and the 2D models in a single computation. Following all these steps of developments of hydrodynamic models, a suitable tool for parameters automatic calibration has been set up. Usable with any solver of the WOLF package, it allows the calibration of the whole of the physical parameters present in the hydrodynamic solvers, as well as the resolution of any calibration or optimization problem with these models. The final parallelisation of this Genetic Algorithms based tool opens the door to complex applications, such as the hydraulic design of structures on the basis of two-dimensional simulations of free surface turbulent flows on multiblock grid. At each step, the validation of the developments is rigorously performed through the comparison of the numerical results with those from test cases, theoretical or experimental, original or from the literature. Several practical applications examples illustrate the whole of the developed modelling potentialities and prove their applicability and efficiency facing the most of practical problems encountered by hydraulic engineers.

The hydraulic and removal efficiencies of a Horizontal Flow Treatment Wetlands (HF TW) were investigated through an internal three dimensional grid of sampling ports. Tracer tests and regular monitoring of water quality parameters were performed. Results demonstrated that the HF TW has generally good hydraulic and volumetric efficiencies, with relatively low dead zones. The application of models developed by chemical engineering provided the number of tanks in series and the calculated detention times which were input as parameters in the multi flow with dispersion hydraulic model. This second model identified that water fluxes were not homogeneous with depth inside the TW and 60% of the flow was along the bottom layer of the bed. It also indicated the water flow velocities, which were faster on the bottom of the bed, and the axial dispersion, which was higher where flow velocity was lower. The reviewed inflow rate distribution allowed review for all layers of the nominal detention time and of the hydraulic indexes, which are developed by the chemical engineering theory, and based on the incorrect assumption of homogeneous systems. The P-k-C* degradation model was applied in order to define degradation k-rate values of BOD and COD and the frequency distribution profiles were developed. The degradation rate coefficients for BOD ranked from less than 10 m/yr to more than 300 m/yr. Significant higher degradation rates were observed for all the bottom layers and for the closest sampling line from the inlet. The results of COD were similar to those observed for BOD. Finally, the specific pollutants of nitrogen and phosphorus were analysed for total nitrogen (TN) and total phosphorus (TP). The global trend of the TN degradation coefficient values was a slow and regular decrease over length, having systematically higher degradation coefficients for the bottom layers. The saturation of the media sites for sorption capacity of TP was demonstrated being in progress. Les efficacités hydraulique et épuratoire du filtre planté à écoulement horizontal ont été mesurées grâce à un dispositif de prélèvement tridimensionnel dans un filtre en fonctionnement. Des essais de traçage ont été réalisés ainsi que des analyses régulières de la qualité de leau en tous les points déchantillonnage. Les résultats ont démontré que le filtre présente des efficacités hydraulique et volumétrique généralement bonnes, avec relativement peu de zones mortes. Lapplication des modèles mathématiques issus du génie chimique ont permis la détermination du nombre N de réservoirs en série et de calculer les temps de séjour ; ceux-ci ont ensuite été introduit dans le modèle hydraulique de dispersion multi-couches. Ce second modèle a identifié que les flux en eau nétaient pas homogènes avec la profondeur à lintérieur du système et que 60 % de lécoulement se concentraient dans la couche de fond du filtre. Il a également fournit les vitesses découlement par couche, qui se présentent comme plus rapides dans le fond du filtre, et la dispersion axiale, dont la valeur augmente en corrélation avec une diminution de la vitesse. La distribution du débit par horizon a pu être revue comme non homogène avec la profondeur et a permis de recalculer les index hydrauliques et les temps de séjour par horizon en évitant lhypothèse erronée des modèles du génie chimique de lhomogénéité du système. Le modèle de dégradation des polluants P-k-C* a été appliqué pour définir les valeurs des coefficients de dégradation de la DBO, la DCO et lazote total. Les profils de distribution de fréquence ont été dressés. Les différences de coefficients de dégradation de la DBO sont très hautement significatives avec la profondeur et la longueur, en ayant des valeurs plus importantes pour la couche de fond et pour la première ligne de prélèvement la plus proche de la zone de distribution de leffluent dans le filtre. Les résultats obtenus pour la DCO ont été semblables. Les coefficients de dégradation de lazote total ont montré une décroissance lente et régulière avec la longueur, et des valeurs systématiquement plus importantes avec la profondeur. Finalement, il a été démontré que les sites dadsorption du phosphore total sont en cours de saturation.

La présente étude porte sur l'hydrogéologie des formations sédimentaires de couverture d'âge Néoprotérozoïque (Infracambrien) à Tertiaire de la plaine du Gondo. La superficie de la plaine du Gondo atteint 30 000 km2, partagée entre le Mali et le Burkina Faso. Elle renferme une nappe libre dans les dépôts détritiques du Continental Terminal (CT) et une nappe plus profonde qui est en charge. La nappe libre est en dépression piézométrique et le toit de la nappe s'approfondit au fur et à mesure que l'on va vers le centre du bassin du Gondo. La présente étude traite du contexte naturel de la plaine du Gondo en général, et du bassin du Sourou à Lery en particulier, avant de s'appesantir sur la lithologie et la géométrie du système aquifère de la plaine du Gondo. Ensuite, les résultats des diverses investigations utilisant une approche interdisciplinaire (télédétection, géophysique, chimie des majeurs, géochimie isotopique, pompage d'essai, modélisation hydrogéologique) ont permis de mieux comprendre la recharge des nappes et leur fonctionnement hydrodynamique. Des solutions d'exploitation de la nappe profonde sont proposées en vue de satisfaire les besoins d'alimentation en eau potable en milieu urbain et semi-urbain. Sous l'effet conjuguée de l'infiltration dans les mares endoréiques et bas-fonds et en dessous du lac du barrage du Sourou, une hausse piézométrique de près de 0,9 m à 5 m / an est enregistrée. Un modèle de recharge des nappes en régime permanent a été utilisé. Ce modèle prend en compte l'exfiltration par évaporation et l'infiltration supplémentaire dans les cônes d'épandage en bordure du bassin du Gondo.

Les phénomènes d'érosion sont une cause majeure dans la rupture des ouvrages hydrauliques tels que les digues ou les barrages. Le territoire compte 10000 km de digues, d'où un enjeu important en terme de risque d'inondation. Une approche plus générale est essentielle à une compréhension fine des phénomènes. La modélisation des mécanismes de l'érosion doit permettre l'élaboration de solutions pour sécuriser les ouvrages et les rendre moins sensibles à ces phénomènes. Ce travail est consacré à la modélisation d'une interface entre un milieu de sol et un milieu d'écoulement. Le mouvement de cette interface est régi par l'érosion du sol sous l'effet de l'écoulement. L'écoulement est calculé par les équations de Stokes munies d'un terme de pénalisation permettant la représentation unifiée des comportements des deux milieux. Dans une première partie, nous aborderons la problématique de la représentation de l'érosion. La seconde partie traite plus précisément de la représentation unifiée des comportements des milieux, notamment grâce à l'utilisation des méthodes de pénalisation. La troisième partie constitue une revue des principales méthodes permettant le suivit d'une interface. Une description précise de la méthode Level Set est donnée. Dans la quatrième partie, nous abordons la résolution numérique des équations d'écoulement. La cinquième partie rassemble l'ensemble des tests de la méthode. Enfin, la dernière partie est consacrée à l'application du modèle sur trois cas : le " Hole Erosion Test ", le phénomène de suffusion et l'érosion d'un lit de rivière autour d'une pile de pont.

Dans le cadre de l'amélioration des méthodes de modélisation du comportement hydraulique des moteurs à ergols liquides, il est important de disposer des méthodes d'identification expérimentale adaptées au besoin de caractérisation des fonctions de transfert des pompes en régime cavitant, notamment pour prédire les risques d'occurrence des phénomènes de couplage fluide-structure comme le Pogo. La matrice de transfert dynamique d'une pompe relie dans le domaine fréquentiel les fluctuations complexes (module + phase) de débit et de pression à l'entrée avec celles de sortie de la pompe. Par ailleurs, les mécanismes physiques qui régissent la réponse d'une pompe soumise à des fluctuations de pression et/ou de débit restent méconnus, ce qui rend aujourd'hui très difficile les prévisions théoriques. Dans le cadre des travaux développés, une nouvelle boucle d'essais a été mise en place permettant la mesure de la pression fluctuante et l'évaluation du débit fluctuant par une approche d'intensimétrie hydroacoustique. Une méthodologie d'essais a été mise en place sur le banc de test pour caractériser les fonctions de transfert de différents organes hydrauliques, notamment en régime cavitant, comme une pompe centrifuge. Ces travaux expérimentaux ont permis également d'alimenter et de valider les modèles de type " système " RLC utilisés pour la simulation du comportement acoustique du circuit. Les méthodologies d'essais, de post-traitement et de modélisation développés dans le cadre de cette thèse seront appliquées par l'industriel lors des essais dédiés sur le banc de test à des composants moteur.

Ce travail concerne l'utilisation de la technique de Fracturation Hydraulique (FH) pour exploiter l'énergie géothermique des réservoirs profonds de roches sèches chaudes (HDR). La fracturation hydraulique est réalisée par injection de fluides géothermiques dans des réservoirs partiellement fracturés de faible perméabilité. Les fluides à haute pression sont destinés à faire évoluer les fissures et leur connectivité. Les valeurs de débit/pression auxquelles les fluides géothermiques doivent être pompés, ainsi que le calendrier de pompage pour initier la fracturation hydraulique, dépendent principalement des conditions géostatiques existantes (contraintes géostatiques, pression fluide et température initiales de l'HDR) ainsi que des propriétés des fissures de l'HDR (longueur, épaisseur, densité et distribution directionnelle initiales moyennes de fissures). Tous ces éléments, en sus de leurs effets sur la stabilité des forages, sont analysés dans cette recherche. Des modèles de fracturation, qui sont capables de suivre l'évolution des fissures dans toutes les orientations spatiales possibles, sont utilisés pour obtenir le tenseur anisotrope de perméabilité. Ces modèles sont intégrés dans un code domestique d'éléments finis qui est développé pour résoudre des problèmes aux limites thermo-poroélastiques. Pour supprimer/diminuer les oscillations qui accompagnent les solutions paraboliques et/ou hyperboliques lors de la convection forcée, plusieurs techniques de stabilisation ont dû être implémentées.

Ce travail s’inscrit dans une tendance historiographique récente visant à mieux comprendre la « vocation industrielle (…) beaucoup trop souvent sous-estimée » du Sud-Ouest de la France. Il tente pour cela de reconstituer les trajectoires des activités textiles à l’œuvre dans le département des Basses-Pyrénées au cours de la période contemporaine, à partir de l'étude des entreprises et de leur environnement, des statistiques et discours sur l'industrialisation, et des « cycles de vie » des articles fabriqués et commercialisés. L’étude met d’une part en évidence l’importance, la structuration et la diversité des activités textiles : les Basses-Pyrénées se caractérisent en effet par la grande variété des fabrications, en renouvellement constant (tissage, bonneterie, confection), et s’imposent même comme le principal espace de production de certains articles à la diffusion internationale (bérets, linge basque, lainages des Pyrénées, sandales). Elles présentent en parallèle des filières complètes, parfois uniques dans le Sud-Ouest, pour différentes matières premières (coton, lin, laine, jute).D’autre part, le maintien sur la longue durée des organisations productives diffuses, le recours massif à l’énergie hydraulique et la domination de la petite entreprise familiale s’articulent avec l’adoption souvent rapide des innovations techniques et une forte capacité de reconversion de multiples centres industriels ruraux ou urbains.Enfin, bien que limitée par la concurrence des grandes régions industrielles françaises vis-à-vis desquelles elles peuvent apparaître comme périphériques, les trajectoires textiles des Basses-Pyrénées parviennent à s’inscrire dans différentes échelles, du local à l’international. Elles témoignent d’échanges techniques et l’élargissement du capitalisme familial s’inscrivent dans des logiques endogènes et exogènes. L’exploitation des opportunités régionales joue à ce titre un rôle prépondérant : proximité de l’Espagne, cohérence vestimentaire d’une vaste aire pyrénéenne, extension des débouchés par l’émigration et la demande touristique, exploitation de l’aura régionaliste du Béarn, du Pays Basque et des Pyrénées. La position excentrée des Basses-Pyrénées, tout en étant une source de marginalisation économique, est un facteur de dynamisme industriel certes limité mais solidement inscrit dans la longue durée. / This study comes within the scope of a recent historiographical trend to change perception of the industrialisation of the South-West of France which is often approached from the perspective of underdevelopment . It aims at reconstructing the « paths » of different textile industries in the Basses-Pyrénées / Pyrénées-Atlantiques from late 18th century to late 20th century, based on the study of produce, markets, producting organisations, firms, views and representations. The research points out the importance, structuring and the diversity of textile industries : complete industrial processes can be found in Basses-Pyrénées for several fibres (cotton, linen, wool, jute) and some products are produced and sold by millions (berets, sandals and, to a leaser extent, basque linen, wool of the Pyrénées...). The maintain of “soft” producting organisations, the massive use of hydraulic power and the domination of family businesses are also linked with fast technical changes and strong ability of industrial reconversion inside various industrial countries and small or middle-sized towns.In the end, these industrial “paths” reveal the economical integration of Basses-Pyrénées on different scales, from local to international, as regards raw material, markets or technical, or financial and human movements. The exploitation of regional opportunities plays a major role in the continuity of these textile industries : proximity of Spain, continuing consumption in the pyrenean area, emigration to America, tourism, regionalistic vogue for Bearn, Basque Country and Pyrénées.