Ce travail de thèse est consacré à la modélisation et simulation numérique d'écoulements diphasiques liquide-gaz mettant en jeu des transferts de chaleur. La simulation de configurations où la prise en compte des effets de compressibilité de la phase gazeuse est indispensable (micropompes, microactionneurs, etc...) a nécessité l'utilisation d'un modèle original, considérant le liquide incompressible et le gaz compressible sous l'hypothèse faible Mach. Lors de cette thèse, ce modèle a été implémenté dans un code diphasique prenant en compte l'interface à l'aide d'une méthode de front-tracking. Des cas tests ont été développés spécifiquement afin de vérifier la conservation de l'énergie pour des configurations de complexité croissante. Les résultats des cas tests ont permis de mettre en évidence la difficulté à assurer la conservation de l'énergie lorsque l'interface n'est pas discontinue mais lissée, comme c'est le cas dans la méthode de front-tracking standard. Une méthode de traitement d'interface hybride a été proposée, rétablissant le caractère discontinu de l'interface avec la reconstruction d'une fonction indicatrice de phase échelon, tandis que le déplacement de l'interface est assuré d'un pas de temps à l'autre à l'aide du front-tracking. Les résultats obtenus avec cette nouvelle méthode hybride sont très satisfaisants, la méthode hybride permettant d'assurer la conservation de l'énergie et de la masse avec précision dans les simulations.

L'application d'un champ électrique en extraction liquide-liquide, permet d'améliorer les coefficients de transfert et de jouer sur les aires interfaciales par des effets de dispersion et de coalescence. La possibilité d'appliquer un champ à un milieu liquide-liquide, dépend des propriétés de ce milieu. Une méthode de mesure de la permittivité et de la conductivité est mise au point et appliquée aux milieux diphasiques semblables à ceux traités dans les colonnes d'extraction liquide-liquide. Un nouveau type de cellule à l'interface plane permet l'étude du transfert sous champ électrique pour les systèmes de l'hydrométallurgie. L'influence favorable du champ est attribuée à un phénomène d'orientation des molécules. Une colonne à plateaux électrifiés, mettant en jeu des phénomènes de coalescence et de dispersion est réalisée et testée. Avec la meilleure géométrie d'électrodes, elle permet de multiplier l'efficacité par quatre. Une petite installation pilote capable de traiter des débits beaucoup plus importants, est réalisée et permet d'approfondir l'étude hydrodynamique sous champ électrique.

Les travaux menés, depuis plusieurs années, au CORIA ont abouti à la construction d'un outil numérique (ARCHER) permettant la simulation numérique directe d'écoulements diphasiques et notamment l'atomisation d'un jet liquide à haute vitesse. Ce type de simulation permet de capturer les phénomènes d'atomisation au voisinage de l'injecteur difficilement caractérisables par les outils expérimentaux actuels. Ces simulations requièrent des conditions d'injection délicates à évaluer a priori car elles dépendent des caractéristiques de l'écoulement au sein de l'injecteur. Or, certains jets présentent une grande sensibilité à ces conditions d'injection. Dès lors, il est nécessaire de simuler l'écoulement au sein de l'injecteur afin d'appréhender la nature de cette sensibilité. L'utilisation d'un maillage cartésien par le code ARCHER conjuguée à la volonté de simuler le système d'atomisation dans son ensemble ont orienté ces travaux vers l'utilisation d'une méthode de frontière immergée. Ces travaux ont ainsi permis de reproduire des écoulements au sein d'injecteurs de forme quelconque tout en conservant le maillage cartésien d'origine, précieux tant pour l'efficacité du solveur que pour sa précision. Dans un premier temps, l'implantation dans le code ARCHER d'une méthode de frontière immergée a été réalisée et testée sur des configurations de canal et de conduite et de l'écoulement autour d'un cylindre. L'application de cette méthode a porté sur la simulation de l'écoulement au sein d'un injecteur triple disque mono-trou et a notamment permis de caractériser l'origine de l'écoulement secondaire formé dans l'orifice de décharge. Afin d'évoluer vers la construction d'un outil numérique capable de simuler le système d'atomisation dans son ensemble, un couplage entre la méthode de frontière immergée et la méthode Ghost fluid a été nécessaire. La version bi-dimensionnelle développée a été testée sur la relaxation d'une goutte posée sur une paroi. Cette version a permis de simuler des écoulements au sein de canaux à différents rapports de longueur sur diamètre et l'écoulement au sein d'une buse convergente. La simulation simultanée de l'écoulement interne et externe a permis de lier les fluctuations de vitesses des écoulements internes à la création de surface engendrée sur les écoulements externes.

Cette étude concerne la caractérisation des écoulements diphasiques bouillants convectifs à haute pression rencontrés dans les réacteurs à eau sous pression (REP). La simulation de ces écoulements est aujourd'hui identifiée comme une voie possible d'amélioration pouvant conduire à la compréhension des mécanismes physiques menant à la crise d'ébullition en réacteur. La première partie de ce travail présente un modèle bi-dimensionnel quasi-établi capable de prédire de façon indépendante les profils de température et de vitesse dans un écoulement diphasique. Le découplage des équations implique de disposer de paramètres d'entrée (taux de vide, vitesse). Ce modèle est basé sur une approche de type modèle de mélange et sur la fermeture des termes de transport turbulent avec le concept de viscosité turbulente. La seconde partie généralise le modèle au cas bi-dimensionnel non-établi en proposant un outil qui résout de façon couplée toutes les équations de bilan et qui est basé sur l'utilisation d'un modèle original de type modèle homogène local avec relaxation thermodynamique. Une confrontation des résultats du modèle à des résultats expérimentaux fournis par la banque de données DEBORA a révélé que notre approche semblait suffisante pour rendre compte d'une bonne partie des données expérimentales en conditions REP. Mais néanmoins qu'elle présentait quelques limites dans des conditions poches du flux critique. Ce travail a permis de mettre en évidence les paramètres sensibles du modèle qui sont aujourd'hui bien identifiés à savoir les mécanismes de transport turbulent d'énergie et le choix du temps de relaxation. / This study is a contribution to the modelling of multidimentional high pressure boiling flows relative to PWR. Numerical simulation of such two-phase flows is considered to be an interesting way for the DNB understanding. The first part of this study exposes a two-dimentional steady state twophase flows model abble to predict velocity and temperature profiles in tube. The mixture balanced equations are used with the eddy diffusivity concept to close the turbulent transport terms. The second part is devoted to the development of the model in the general two dimentional case. Contrary to the steady state model, this model is indenpendant of experimental data and implies the use of an original local homogeneous relaxation modèle (HRM). The results obtained from the comparison with the data bank DEBORA reveals that in a mixture approch two submodel are sufficients to obtain a physial good description of turbulent boiling flows. Some limitations appear at conditions close to DNB conditions. The turbulent closures and the relaxation time in the HRM model have been clearly identified as the most important and sensitive parameters in the model.

Cette thèse porte sur l'étude des ballotements dans les réservoirs de satellites à poste, lors des phases de manoeuvre à faible accélération. En effet la bulle de gaz d'hélium servant à pressuriser le réservoir se met en mouvement, générant ainsi des perturbations sur la stabilité globale du satellite. Afin de mener à bien cette étude, des méthodes numériques ont été développées, avec une méthode de frontières immergées pour prendre en compte les parois du réservoir.Le code est utilise la méthode Level-Set pour capturer l'interface, et gère les sauts à l'aide de la méthode Ghost-Fluid. Un solveur BlackBox Multigrid est également développé pour améliorer lesperformances de calcul. Une étude est présentée dans le dernier chapitre pour définir quelques lois de comportements en fonction des vitesses et accélérations générées lors des manoeuvres. / The core study of this PhD thesis is the sloshing in satellite tanks, during low acceleration maneuvers. Indeed the helium bubble used to pressurize the tank moves, thus generating perturbations on the global stability of the satellite. In order to understand this problem, numerical schemes have been developed, such as an immersed boundary method to model the tank wall. The numerical tool uses a Level-Set function coupled to a Ghost Fluid Method to track the interface and to account for the jump conditions.A BlackBox Multigrid Solver have been developed to improve computational cost. Finally a study is presented in the last chapter to predict the behaviour of the fluids with a varying rotational speed generated during some classical maneuvers.

Ce travail de thèse est consacré à l'étude des écoulements diphasiques dispersés turbulents gaz/gouttes et plus particulièrement à la modélisation du phénomène de dépôt de gouttes en canal horizontal, dont la compréhension et la prédiction sont essentielles pour de nombreuses applications industrielles. Les gouttes sont supposées de taille plus petite que les échelles de longueur caractéristiques de l'écoulement de gaz turbulent, avec une masse volumique grande devant celle de la phase continue, les forces qui agissent sur les gouttes se limitent ainsi à la traînée, à la poussée d'Archimède et à la gravité. Le taux de présence de la phase dispersée est suffisamment important pour tenir compte de l'influence des gouttes sur la turbulence du gaz (couplage à deux sens), mais suffisamment faible pour pouvoir négliger les collisions entre les gouttes. En écoulement horizontal, le dépôt des gouttes en paroi est piloté par deux mécanismes principaux qui agissent en parallèle : la gravité et la diffusion turbulente/vol libre. Cette physique du dépôt est déclinée en deux volets, avec une première étude à l'échelle 3D locale et une seconde étude à l'échelle système 1D. Dans chacune de ces approches, un modèle pour la vitesse de dépôt de gouttes en paroi est développé, puis validé par comparaison à des données expérimentales. Le modèle de dépôt local, établi sous l'hypothèse d'un film liquide infiniment mince et absorbant, est implanté dans le code de simulation numérique NEPTUNE_CFD, puis validé par comparaison aux données expérimentales de Namie & Ueda, qui étudient le dépôt des gouttes en canal horizontal. Une analyse des équations de transport des principales grandeurs moyennes de l'écoulement, ainsi que des transferts d'énergies entre phases, est menée afin de mettre en évidence les phénomènes de couplage et leurs influences sur la turbulence de la phase continue. Le modèle unidimensionnel, développé dans le cadre d'un besoin industriel, est implanté dans le code CATHARE-3 et est confronté aux données de l'expérience REGARD du CEA Grenoble. / Droplets dispersion and deposition in turbulent duct flows are important processes, occurring in numerous environmental and industrial applications. This work is devoted to the study of gas-droplets flows and, more particularly, the objective is to improve the droplets deposition modeling in horizontal flows. Droplets are supposed to be smaller than the Kolmogorov scale, with a density large compared to the density of the gas phase. Under these assumptions, the motion of a droplet is considered to be governed by the drag force, the buoyancy force, and the gravity. Dilute incompressible and isothermal gas-droplets flows are studied, so inter-particle collisions are neglected but two-way coupling is retained, which means that modulation of turbulence by the particles is accounted for. In horizontal flow, droplets reach the wall under the actions of the gravitational settling and the turbulent diffusion. Two approaches will be used in developing this deposition physics with a first study at the 3D local scale and a second one at the 1D scale, realized for an industrial need. For each case, a model is developed for the mean deposition velocity of the droplets, with is implemented in a numerical simulation tool and then validated by comparison to experimental data. The local deposition model is established under the assumptions that the liquid film is extremely thin and perfectly absorbing and is implemented in the Neptune_CFD code. The experience carried out by Namie & Ueda, which consist in small droplets deposition from a turbulent dispersed flow in a horizontal rectangular duct, is simulated. An analysis of the interphase transfer terms in the kinetic energy equations shows the interactions between the dispersed phase and the continuous one and the impact of these phenomena on the turbulence of the gas phase is pointed out. The 1D deposition model is developed for the CATHARE-3 code and experimental data from the REGARD facility of the CEA Grenoble are used for validation.

Cette thèse est consacrée à l'étude des propriétés des transports diphasiques au sein des couches de diffusions (Gas Diffusion Layer = GDL) des piles à combustible PEMFC (Proton Exchange Membrane Fuel Cells). La GDL est faite d'une structure fibreuse (dont l'épaisseur est de quelques centaines de micromètres) traitée généralement avec une matière hydrophobe. Des images numériques de la GDL réelle obtenues par tomographie X sont d'abord analysées afin d'étudier des propriétés telles que la porosité, la perméabilité, ou le tenseur de diffusion. L'écrasement de la GDL est ensuite simulé en utilisant un algorithme comprimant les fibres dans un plan transversal. Les transports diphasiques (invasion quasi statique d'eau liquide) sont modélisés dans des réseaux de pores, milieux représentatifs de l'espace poreux de la GDL, en relation avec le problème de la gestion de l'eau dans les piles PEMFC. Deux algorithmes d'invasion, dénommés algorithmes séquentiel et cinétique, sont développés et comparés pour analyser les distributions de phases au sein des GDL. Un point clé est que l'eau rentre dans la couche poreuse par divers points d'injection indépendants, conduisant à la possibilité de multiples points de percée. Des expériences sur un système microfluidique sont conduites pour valider les algorithmes utilisés. Une étude statistique est menée pour caractériser le nombre de points de percée, les profils de saturation, l'accès au gaz, le transport diffusif, de même que l'influence du piégeage et de la mouillabilité mixte. / This thesis is devoted to the study of transport properties and two-phase flow in the Gas Diffusion Layer (GDL) of Proton Exchange Membrane Fuel Cells (PEMFC). A GDL is a thin fibrous structure (a few hundreds μm thick) treated generally with a hydrophobic agent. Numerical images obtained from X-ray computed tomography X are first exploited to study properties such as the porosity, permeability and diffusion tensors of a real GDL microstructure. The effect of GDL compression is also investigated using an algorithm mimicking the compression in GDL through plane direction. Then two phase flow (quasi-static water invasion) is studied in relation with the water management problem in PEMFC, using a structured pore network representation of the pore space. Two invasion algorithms, referred to as the sequential and the kinetic algorithm respectively, are developed and compared to study the fluid distributions within the GDL. A key point is that water enters the porous layer through multiple independent inlet injection points, leading to the possibility of many breakthrough points. Experiments are conducted on a microfluidic device to validate the algorithms. A numerical statistical study is performed to characterize the breakthrough point statistics, saturation profiles, gas access, diffusion transport as well as the influence of trapping and mixed wettability.

Ce travail a pour thème la simulation numérique d’écoulements diphasiques dans un contexte industriel. En effet, la simulation d’écoulements diphasiques est un domaine qui présente de nombreux défis, en raison de phénomènes complexes qui surviennent, comme la cavitation et autres transferts entre les phases. En outre, ces écoulements se déroulent généralement dans des géométries complexes rendant difficile une résolution efficiente. Les modèles que nous considérons font partie de la catégorie des modèles à interfaces diffuses et permettent de prendre en compte aisément les différents transferts entre les phases. Cette classe de modèles inclut une hiérarchie de sous-modèles pouvant simuler plus ou moins d’interactions entre les phases. Pour mener à bien cette étude nous avons en premier lieu comparé les modèles diphasiques dits à quatre équations et six équations, en incluant les effets de transfert de masse. Nous avons ensuite choisi de nous concentrer sur le modèle à quatre équations. L’objectif majeur de notre travail a alors été d’étendre les méthodes aux résidus distribués à ce modèle. Dans le contexte des méthodes de résolution numérique, il est courant d’utiliser la forme conservative des équations de bilan. En effet, la résolution sous forme non-conservative conduit à une mauvaise résolution du problème. Cependant, résoudre les équations sous forme non-conservative peut s’avérer plus intéressant d’un point de vue industriel. Dans ce but, nous utilisons une approche développée récemment permettant d’assurer la conservation en résolvant un système sous forme non-conservative, à condition que la forme conservative soit connue. Nous validons ensuite notre méthode et l’appliquons à des problèmes en géométries complexes. Finalement, la dernière partie de notre travail est dédiée à étudier la validité des modèles à interfaces diffuses pour des applications à des problèmes industriels réels. On cherche alors, en utilisant des méthodes de quantification d’incertitude, à obtenir les paramètres rendant nos simulations les plus vraisemblables et cibler les éventuels développements pouvant rendre nos simulations plus réalistes. / The topic of this thesis is the numerical simulation of two-phase flows in an industrial framework. Two-phase flows modelling is a challenging domain to explore, mainly because of the complex phenomena involved, such as cavitation and other transfer processes between phases. Furthermore, these flows occur generally in complex geometries, which makes difficult the development of efficient resolution methods. The models that we consider belong to the class of diffuse interface models, and they allow an easy modelling of transfers between phases. The considered class of models includes a hierarchy of sub-models, which take into account different levels of interactions between phases. To pursue our studies, first we have compared the so-called four-equation and six-equation two-phase flow models, including the effects of mass transfer processes. We have then chosen to focus on the four-equation model. One of the main objective of our work has been to extend residual distribution schemes to this model. In the context of numerical solution methods, it is common to use the conservative form of the balance law. In fact, the solution of the equations under a non-conservative form may lead to a wrong solution to the problem. Nonetheless, solving the equations in non-conservative form may be more interesting from an industrial point of view. To this aim, we employ a recent approach, which allows us to ensure conservation while solving a non-conservative system, at the condition of knowing its conservative form. We then validate our method and apply it to problems with complex geometry. Finally, the last part of our work is dedicated to the evaluation of the validity of the considered diffuse interface model for applications to real industrial problems. By using uncertainty quantification methods, the objective is to get parameters that make our simulations the most plausible, and to target the possible extensions that can make our simulations more realistic.

Les besoins énergétiques de la population mondiale ne cessent d’augmenter. Les prévisions indiquent par exemple une forte croissance de la demande du secteur du transport aéronautique. La recherche de systèmes toujours plus performants et moins polluants est nécessaire. Des nouveaux concepts pour la combustion ont été mis au point et appliqués aux turbines à gaz. Parmi eux il existe ceux basés sur la combustion en prémélange pauvre ou en prémélange pauvre pré-vaporisé dans le cas où le carburant utilisé est liquide. Les nouveaux systèmes énergétiques basés sur la combustion en régime pauvre sont prometteurs pour satisfaire les futures normes d’émissions polluantes, mais ils sont plus sensibles aux instabilités de combustion qui limitent leur plage de fonctionnement et peuvent détériorer irréversiblement ces systèmes. Dans ce domaine il reste des questions à aborder. En particulier celle du comportement des flammes tourbillonnaires en combustion diphasique soumises à des perturbations acoustiques. La plupart des moteurs aéronautiques utilisent des flammes de ce type, cependant leur dynamique et leurs interactions mutuelles, quand elles subissent les effets d’une perturbation acoustique, sont loin d’être bien comprises. Ce travail aborde ces questions et apporte des éléments de compréhension sur les mécanismes pilotant la réponse de l’écoulement diphasique et de la flamme, ainsi que des éléments de validation des modèles de prédiction des points de fonctionnement instables. TACC-Spray est le banc expérimental utilisé pour ce travail. Il a été conçu et développé au sein du laboratoire CORIA lors de ce doctorat qui s’inscrit dans le cadre du projet ANR FASMIC. Le système d’injection qui équipe ce banc expérimental reçoit trois injecteurs tourbillonnaires alimentés en combustible liquide (ici n-heptane), développés par le laboratoire EM2C. Ils sont montés en lignes dans le banc, celui-ci représentant ainsi un secteur d’une chambre annulaire. Le montage étant complexe et nouveau, un travail de développement de solutions techniques a été fait pour rendre possible l’équipement du TACC-Spray avec des capteurs de pression, température, photomultiplicateur ainsi que des diagnostiques optiques performants (e.g. LDA, PDA, imagerie à haute cadence). Pour cette étude, le système énergétique, composé par l’écoulement diphasique et la flamme, a été soumis à l’impact d’un mode acoustique transverse excité dans la cavité acoustique. La réponse du système a été étudiée en fonction de son positionnement dans le champ acoustique. Trois bassins d’influence du champ acoustique sur le système énergétique ont été choisis, à savoir: (i) le ventre de pression acoustique caractérisé principalement par des fortes fluctuations de pression, (ii) le ventre d’intensité acoustique présentant de forts gradients de pression et vitesse acoustique, (iii) le ventre de vitesse acoustique avec de fortes fluctuations de vitesse où la fluctuation de pression est résiduelle. L’approche de cette étude a consisté à étudier en premier lieu le système de référence en absence de forçage acoustique, les résultats sont recueillis dans la Partie I de ce manuscrit. En deuxième lieu le système énergétique est placé à chacune des positions d’intérêt dans le champ acoustique et la réponse de l’écoulement d’air sans combustion, la réponse de l’écoulement diphasique avec combustion et finalement celle des flammes, sont étudiées systématiquement. Les résultats de l’étude avec forçage acoustique sont rassemblés dans la Partie II du manuscrit. / The energy needs of population around the word are continuously increasing. For instance, forecasts indicates an important grow of the request of the aeronautic transportation sector. It is necessary to continue the research efforts to get more performants and less contaminating systems. New concepts for combustion have been developed and introduced to the gas turbine industry. Among these concepts it is found technologies based on lean-premixed combustion or lean-premixed prevaporized combustion when liquid fuels are employed. These novel energetic systems, making use of lean combustion, are promising to meet the future norms about pollutant emissions, but this make them more sensitive to combustion instabilities that limit their operating range and can lead to irreversible damage. In this domain, many questions still need to be considered. In particular that of the behavior of two-phase flow swirling flames subjected to acoustic perturbations. Indeed most of aero-engines operate with this type of flames, but the dynamics and mutual interaction of these flames, as they are submitted to acoustic perturbation, are not yet well understood. This work addresses these issues and gives some understanding elements for the mechanisms driving the response of the flow and of the flame to acoustic perturbations and delivers data to validate models predicting unstable operating points.The experimental bench employed for this work is TACC-Spray. It has been designed and developed in the CORIA laboratory during this PhD thesis which is inscribed in the framework of the ANR FASMIC project. The injections system that equips this bench is composed by three swirled injectors fed with a liquid fuel (here n-heptane), developed by the EM2C laboratory. They are linearly arranged in the bench such that this represents an unwrapped sector of an annular chamber. The setup, being new and complex, needed technical solutions developed during this work and applied then in order to equip TACC-Spray with pressure and temperature sensors, a photomultiplier as well as adequate optic diagnostics (LDA, PDA, high speed imaging systems). In this study, the energetic system, composed by the two-phase swirling flow and the spray flame, has been submitted to the impact of a transverse acoustic mode excited within the acoustic cavity. The system response has been studied as a function of its location in the acoustic field. Three basins of influence of the acoustic field on the energetic system have been chosen, namely: (i) the pressure antinode characterized mainly by strong pressure fluctuations, (ii) the intensity antinode where important acoustic pressure and velocity gradients are present, (iii) the velocity antinode with strong velocity fluctuations where the acoustic pressure is residual. The approach of the study presented here is to investigate in first place the energetic system free of acoustic forcing. The results concerning this first study are presented in the Part I of this manuscript. In second place, the energetic system is placed in each of the location of interest within the acoustic field and the response of the air flow without combustion, that of the two-phase flow with combustion and finally that of the spray flames, are systematically investigated. The results of the study under acoustic forcing are shown in Part II of the manuscript.

Au cours de l’utilisation de certains systèmes propulsifs, tels que les moteurs fusées cryotechniques ou les moteurs Diesel, le point de fonctionnement peut varier sur une large plage de pressions. Ces variations de pression peuvent conduire à un changement de régime thermodynamique si la pression critique du fluide est franchie, l’injection initialement diphasique devenant alors transcritique. Ce changement modifie la topologie de l’écoulement, ainsi que la dynamique du mélange, ce qui impacte le comportement de la flamme. L’objectif de cette thèse est de développer une méthodologie originale capable de traiter au sein du même solveur des écoulements sous-critiques ainsi que supercritiques. Pour cela, une extension du solveur AVBP-RG aux écoulements diphasiques sous-critiques est proposée, basée sur des modèles d’interface diffuse. Les développements nécessaires à l’intégration de ces modèles dans le cadre du solveur aux éléments finis sont effectués. Des simulations numériques multidimensionnelles sont ensuite proposées de manière à confronter le modèle à des données exprérimentales, vis-à-vis desquelles un bon accord est observé. Cesrésultats offrent des perspectives encourageantes vers de futures améliorations du modèle et des applications à des configurations industrielles complexes. / In various industrial combustion devices, such as liquid rocket engines at ignition or Diesel engines during the compression stage, the operating point varies over a wide range of pressures. These pressure variations can lead to a change of thermodynamic regime when the critical pressure is exceeded, switching from two-phase injection to transcritical injection. This change modifies the topology of the flow and the mixing, thereby impacting the flame dynamics. The objective of the present Ph.D thesis is to develop an original methodology able to address both subcritical and supercritical flows within the same solver. To achieve this, an extension of the real gas solver AVBP-RG to subcritical two-phase flows is provided, based on diffuse interface models. The required developments for the integration of such models into the finite-element framework of the solver are provided. Multidimensional numerical simulations are led in order to confront the model with experimental data, with which good agreement is observed. These results offer encouraging perspectives regarding further enhancements of the model and applications to complex industrial cases.