Letteratura scientifica selezionata sul tema "Maillage en blocs structurés"

La caractéristique principale des sociétés modernes est l’extraordinaire développement des organisations qui investissent, contrôlent, gèrent l’ensemble des registres de la vie sociale : les rapports sociaux sont essentiellement des rapports organisés, c’est-à-dire des rapports structurés par les organisations. L’auteur propose ici d’analyser les conséquences de cette recomposition. Si la société de classes liée au capitalisme industriel était une société rigide, hiérarchisée et répressive, la société duale qui se profile actuellement est éclatée, fragile et oppressive. D’un côté, elle tend à faire de l’individu son propre réfèrent en le renvoyant constamment à lui-même, à la nécessité de « s’auto-réaliser », de l’autre, elle met en place un maillage organisationnel de l’espace social pour maîtriser les conséquences de son développement.

The Maillard reaction between reducing sugars and amino groups of biomolecules generates complex structures known as AGEs (advanced glycation endproducts). These have been linked to protein modifications found during aging, diabetes and various amyloidoses. To investigate the contribution of alternative routes to the formation of AGEs, we developed a mathematical model that describes the generation of CML [Nε-(carboxymethyl)lysine] in the Maillard reaction between glucose and collagen. Parameter values were obtained by fitting published data from kinetic experiments of Amadori compound decomposition and glycoxidation of collagen by glucose. These raw parameter values were subsequently fine-tuned with adjustment factors that were deduced from dynamic experiments taking into account the glucose and phosphate buffer concentrations. The fine-tuned model was used to assess the relative contributions of the reaction between glyoxal and lysine, the Namiki pathway, and Amadori compound degradation to the generation of CML. The model suggests that the glyoxal route dominates, except at low phosphate and high glucose concentrations. The contribution of Amadori oxidation is generally the least significant at low glucose concentrations. Simulations of the inhibition of CML generation by aminoguanidine show that this compound effectively blocks the glyoxal route at low glucose concentrations (5 mM). Model results are compared with literature estimates of the contributions to CML generation by the three pathways. The significance of the dominance of the glyoxal route is discussed in the context of possible natural defensive mechanisms and pharmacological interventions with the goal of inhibiting the Maillard reaction in vivo.

AbstractA series of alternative structures for housing are proposed through a design-research prototyping process. Nine proto-structures are developed through the conception and realization of columns cast in concrete. The series explores specific structural principles at 1to1 scale, which are further architecturally tested as speculative towers for urban living at 1to100 scale.At 1 to 1 scale, concrete as process rather than just concrete as material sets the main methodology for the design-research. Trial and error experiments, closely related to the properties of the material (liquid to solid formation), seek to put forward new techniques of formwork design and construction procedures that are more flexible and more sustainable than past and current systems.Design analysis, informed by the work of the early ‘structural rationalists’ (F. Candela, P.L. Nervi, H. Isler, R. Maillard, E. Dieste, et al.), considers the transformation of structural languages in an attempt to revive an architecture for vertical living (point-block).

International audience The computation of determinants intervenes in many scientific applications, as for example in the localization of eigenvalues of a given matrix A in a domain of the complex plane. When a procedure based on the application of the residual theorem is used, the integration process leads to the evaluation of the principal argument of the complex logarithm of the function g(z) = det((z + h)I - A)/ det(zI - A), and a large number of determinants is computed to insure that the same branch of the complex logarithm is followed during the integration. In this paper, we present some efficient methods for computing the determinant of a large sparse and block structured matrix. Tests conducted using randomly generated matrices show the efficiency and robustness of our methods. Le calcul de déterminants intervient dans certaines applications scientifiques, comme parexemple dans le comptage du nombre de valeurs propres d’une matrice situées dans un domaineborné du plan complexe. Lorsqu’on utilise une approche fondée sur l’application du théorème desrésidus, l’intégration nous ramène à l’évaluation de l’argument principal du logarithme complexe de lafonction g(z) = det((z + h)I − A)/ det(zI − A), en un grand nombre de points, pour ne pas sauterd’une branche à l’autre du logarithme complexe. Nous proposons dans cet article quelques méthodesefficaces pour le calcul du déterminant d’une matrice grande et creuse, et qui peut être transforméesous forme de blocs structurés. Les résultats numériques, issus de tests sur des matrices généréesde façon aléatoire, confirment l’efficacité et la robustesse des méthodes proposées.

Les codes de simulation numérique reposant sur des méthodes de type éléments et volumes finis requièrent de discrétiser le domaine étudié – par exemple une pièce mécanique telle qu’un moteur, une aile d’avion, une turbine, etc. – à l’aide d’un maillage. En dimension 3, un maillage est un ensemble composé d’éléments volumique simples, le plus souvent des tétraèdres ou des hexaèdres, qui partitionnent le domaine d’étude. Le choix de tétraèdres ou d’hexaèdres est principalement dicté par l’application (interaction fluide-structure, hydrodynamique, etc.). Si la génération automatique de maillages tétraédriques est un processus relativement maîtrisé aujourd’hui, générer des maillages hexaédriques est toujours un problème ouvert. Ceci est problématique pour les applications qui justement nécessitent impérativement des maillages hexaédriques puisque leur génération se fait de façon semi-automatique, ce qui peut prendre plusieurs semaines à plusieurs mois de temps ingénieur ! Alors que le temps consacré au processus de simulation numérique à proprement parler tend à diminuer du fait de la puissance des machines utilisées, le goulot d’étranglement est désormais dans la préparation des données, à savoir obtenir un modèle de CAO adapté au calcul, puis en générer un maillage.C’est dans ce contexte que s’inscrit la thèse proposée en suivant une approche hybride mêlant :1. Le développement d’algorithmes (semi)-automatiques pour générer et modifier des maillages hexaédriques structurés par blocs ;2. La mise en place d’un logiciel graphique interactif dédié à la manipulation de structures de blocs. Les mécanismes d’interaction seront en outre utilisés pour guider les algorithmes dans leurs prises de décision, que ce soit à l’initialisation (critères à apposer sur des entités particulières de CAO) ou en cours d’algorithme (décision entre plusieurs options sur lesquelles l’algorithme ne peut se prononcer seul).L’objectif de cette thèse n’est donc pas de fournir une solution automatique universelle, ce qui semble inatteignable actuellement, mais plutôt de réduire le temps ingénieur consacré à la génération du maillage en fournissant des outils plus adaptés. Dans cette optique, nous proposons de placer l’étude dans le prolongement de [LED10, KOW12, GAO15, GAO17], où est considéré le problème de simplification et d’enrichissement de maillages hexaédriques par insertion et suppression de couches de mailles. Dans tous ces travaux, les algorithmes proposés sont des algorithmes simples de type « glouton » où le maillage est modifié pas à pas pour converger vers une solution finale Ef : A chaque étape Ei, on fait l’hypothèse que la « meilleure » solution Ef sera obtenue en faisant le choix « optimal » pour Ei. Or en recherche opérationnelle, une telle approche est connue comme perfectible dès lors que le problème d’optimisation traité est non linéaire. L’idée est donc d’utiliser des approches usuelles en recherche opérationnelle et plus spécifiquement des systèmes multi-agents, couplées à des outils interactifs, pour permettre la génération de structures de blocs sur des CA0 complexes
The numerical simulation codes based on element and finite volume methods require discretizing the studied domain - for example a mechanical part such as a motor, an airplane wing, a turbine, etc. - using a mesh. In 3 dimensions, a mesh is a set composed of simple volumic elements, most often tetrahedrons or hexahedra, which partition the field of study. The choice of tetrahedrons or hexahedra is mainly dictated by the application (fluid-structure interaction, hydrodynamics, etc.). If the automatic generation of tetrahedral meshes is a relatively controlled process today, generating hexahedral meshes is still an open problem. This is problematic for applications that just imperatively require hexahedral meshes since their generation is done semi-automatically, which can take several weeks to several months of engineer time! While the time devoted to the digital simulation process itself tends to decrease due to the power of the machines used, the bottleneck is now in the preparation of the data, namely to obtain a CAD model adapted to the computation, then generate a mesh.It is in this context that we propose to develop a hybrid approach combining:1. The development of (semi) -automatic algorithms for generating and modifying block-structured hexahedral meshes;2. The implementation of an interactive graphic software dedicated to the manipulation of block structures. The interaction mechanisms will also be used to guide the algorithms in their decision-making, whether at initialization (criteria to be affixed to particular CAD entities) or in the course of an algorithm (decision between several options on which the algorithm cannot pronounce itself).The objective of this thesis is therefore not to provide a universal automatic solution, which seems unattainable at present, but rather to reduce the engineering time devoted to the generation of the mesh by providing more adapted tools. In this context, we propose to place the study as an extension of [LED10, KOW12, GAO15, GAO17], which considered the problem of simplification and enrichment of hexahedral meshes by insertion and removal of mesh layers. In all these works, the proposed algorithms are simple 'greedy' algorithms where the mesh is modified step by step to converge towards a final solution Ef: At each step Ei, one makes the hypothesis that the 'best' solution Ef will be obtained by making the 'optimal' choice for Ei. However, in operational research, such an approach is known to be perfectible since the problem of optimization treated is nonlinear. The idea is to use usual approaches in operational research and more specifically multi-agent systems, coupled to interactive tools, to allow the generation of block structures on complex CA0s

Ce travail de thèse porte sur la représentation et la génération de maillages hexaédriques structurés par blocs. Il n'existe pas à ce jour de méthode permettant de générer des structures de blocs satisfaisantes pour n'importe quel domaine géométrique. En pratique, des ingénieurs experts génèrent ces maillages avec des logiciels interactifs, ce qui nécessite parfois plusieurs semaines de travail. De plus, l'ajout d'opérations de modification dans ces logiciels interactifs est un travail délicat pour maintenir la cohérence de la structure de blocs et sa relation avec le domaine géométrique à discrétiser. Afin d'améliorer ce processus, nous proposons tout d'abord de définir des opérations de manipulation de maillages hexaédriques se basant sur l'utilisation du modèle des cartes généralisées. Ensuite, en considérant des structures de blocs obtenues à l'aide de la méthode des Polycubes, nous fournissons des méthodes optimisant la topologie de ces structures pour satisfaire des contraintes de nature géométrique. Nous proposons ainsi une première méthode en dimension 2, qui considère une approche locale du problème en s'appuyant sur l'expérience des ingénieurs manipulant des logiciels interactifs. Puis nous proposons une seconde méthode utilisant cette fois la méta-heuristique d'optimisation par colonie de fourmis pour la sélection de feuillets en dimension 3
This thesis deals with the representation and generation of block-structured hexahedral meshes. To date, there is no method for generating satisfactory block structures for any geometric domain. In practice, expert engineers generate these meshes using interactive software, which can take several weeks to complete. Moreover, adding modification operations in these interactive softwares is a delicate task to maintain the coherence of the block structure and its relationship with the geometric domain to be discretized. In order to improve this process, we first propose to define hexahedral mesh manipulation operations based on the use of the generalized map model. Then, by considering block structures obtained using the Polycube method, we provide methods for optimizing the topology of these structures to satisfy constraints of a geometric nature. We propose a first method in dimension 2, which considers a local approach to the problem based on the experience of engineers working with interactive software. We then propose a second method, this time using ant colony optimization meta-heuristics for leaf selection in dimension 3

L'objectif de cette thèse est d'apporter de nouvelles solutions pour améliorer les performances lors de la résolution numérique d'équations aux dérivées partielles en étudiant une méthode de génération de maillages quadrilatéraux basée sur des champs de croix.Notre approche repose sur l'alignement d'un champ de croix donné par rapport au bord du domaine de calcul, suivi du partitionnement de ce dernier en régions à quatre côtés que l'on maillera ensuite en quadrilatères. Cela nous permet de générer un maillage structuré par bloc tout en préservant la position des singularités du champ de croix initial. Pour ce faire, nous mettons en place une étude théorique des champs de croix nous permettant d'identifier les conditions sous lesquelles un champ de croix donné permet effectivement de partitionner et de mailler le domaine sur lequel il est défini.Cette méthode est d'abord présentée pour les domaines simplement connexes, puis généralisée aux domaines non-simplement connexes. Elle nous permet de gérer les domaines constitués de plusieurs matériaux ainsi que les points singuliers de bord qui, en pratique, permettent de délimiter des portions de la frontière du domaine pour prendre en compte des conditions aux limites mixtes dans le cadre des simulations numériques. Enfin, nous mettons en place la discrétisation de la méthode sur des maillages triangulaires, puis nous proposons une généralisation à des surfaces courbes dans l'espace
The objective of this thesis is to provide new solutions to improve performance in the numerical solution of partial differential equations by studying a method for generating quadrilateral meshes based on cross-field techniques.Our approach relies on aligning a given cross field with the boundary of the computational domain, followed by partitioning the domain into four-sided regions which are then meshed into quadrilaterals. This allows us to generate a structured block mesh while preserving the positions of the singularities of the initial cross field. To achieve this, we conduct a theoretical study of cross fields to identify the conditions under which a given cross field effectively partitions and meshes the domain on which it is defined.This method is initially presented for simply connected domains and then extended to non-simply connected domains. It allows us to handle domains consisting of multiple materials as well as boundary singular points, which in practice delineate portions of the domain boundary to account for mixed boundary conditions in numerical simulations. Finally, we discretize the method on triangular meshes and propose a generalization to curved surfaces in space

Le traitement numérique de problèmes de mécanique des fluides s'appui sur des géométries souvent particulières et complexes. Afin de les représenter au mieux, nous avons développé une méthode de raccordement de maillages multiblocs non-conformes. Elle repose sur l'interpolation non-conservative et implicite des variables sur des interfaces entre les blocs. Le traitement des conditions de raccord a été intégré au code de calcul Aquilon en 2D et en 3D, dans le contexte de la méthode du Lagrangien Augmenté servant à la résolution des équations de Navier-Stokes. La validation numérique 2D a été effectuée sur un ensemble de cas test dont les solutions analytiques sont bien connues ainsi que sur les cas de la marche descendante, de la cavité entraînée et de l'écoulement autour d'un cylindre. Les résultats obtenus montrent la faisabilité et le bien-fondé de l'approche proposée.

Le traitement numérique de problème de mécanique des fluides s'appuit sur des géométries souvent particulières et complexes. Afin de les représenter au mieux, nous avons développé une méthode de raccordement de maillages multiblocs non-conformes. Elle repose sur l'interpolation, non-conservative et implicite des variables des interfaces entre les blocs. Le traitement des conditions de raccord a été intégré au code de calcul Aquilon en 2D et en 3D, dans le contexte de la méthode du Lagrangien Augmentée servant à la résolution des équations de Navier-Stokes. La validation numérique 2D a été effectuée sur un ensemble de cas tests dont les solutions analytiques sont bien connues ainsi que sur les cas de la marche descendante, de la cavité entraînée et de l'écoulement autour d'un cylindre. Les résultats
obtenus montrent la faisabilité et le bien-fondé de l'approche proposée.

La modélisation des écoulements en réservoir fracturé s’appuie généralement sur une représentation dite double-milieu. La fiabilité du calcul des transferts matrice-fissure est alors essentielle pour les prévisions de production. A l’échelle de la maille d’un simulateur double-milieu, ces transferts sont formulés en supposant un écoulement pseudo-permanent entre un bloc matriciel représentatif, non discrétisé, d’une part, et le réseau de fissures limitrophes homogénéisé d’autre part. Si elle peut suffire à estimer les échanges monophasiques, une telle formulation se révèle par contre inadaptée à la prévision des mécanismes de production polyphasiques. Cette thèse revisite une méthode de sous-maillage des blocs matriciels et l’optimise en vue de la rendre opérationnelle et intégrable aux simulateurs double-milieu tels qu’utilisés à ce jour par l’industrie pétrolière. Sur le plan théorique, un modèle mixte a été construit en couplant le modèle d’écoulement de fracture à grande échelle à un modèle d’écoulement à l’échelle locale du bloc matriciel désormais discrétisé et non plus comme un simple terme-source. Sur le plan pratique, la mise en œuvre de ce modèle mixte a comporté trois étapes principales, (a) la définition d’une méthode générale de discrétisation du bloc matrice réalisant le meilleur compromis entre caractère prédictif et coût de calcul, (b) la transcription numérique du couplage entre d’une part, le modèle double-milieu tel qu’utilisé actuellement et d’autre part, un modèle d’écoulement diphasique au sein du bloc matriciel maillé, (c) la validation des formules de couplage et des méthodes de résolution pour divers cas d’application qui mettent également en évidence l’impact de certaines hypothèses de travail. Le modèle mixte proposé offre de réelles et décisives perspectives d’amélioration des simulateurs double milieu de l’industrie.

L'évolution constante ainsi que la complexification qui s'en suit des architectures matérielles oblige les personnes développant des codes de simulations scientifiques à une mise à jour perpétuelle de leur logiciel et de leurs connaissances afin de maintenir une bonne exploitation des performances de ces plateformes de calcul. L'hétérogénéité récente de ces plateformes et la multiplication du nombres de coeurs par machine posent notamment une nouvelle problématique aux logiciels largement parallélisés via l'utilisation de la bibliothèque d'échange de message MPI, qui se voient perdre en performance au fur et à mesure de l'augmentation du parallélisme local à exploiter. Pour cs,la transition vers l'exploitation d'un parallélisme hybride via l'association de MPI avec Open MP (en parallélisme de boucles) peine à s'imposer et ne résout pas les soucis de passage à l'échelle causé par un manque d'équilibrage de charge. Pour répondre à ce problème (ainsi que ceux de la portabilité d'un code et sa capacité à exploiter les plateformes hétérogènes),l'utilisation de la programmation par tâches a récemment gagné en popularité à travers trois méthodes: parallélisme par boucle, tâches soumises de manière séquentielle (STF, flux de tâche séquentiel) ou paramétrique (PTG: graphe de tâches paramétré). Cette thèse se concentre sur l'effort de transition de cs, un code de simulation en mécanique des fluides sur maillages non structurés (algèbre linéaire creuse) vers l'utilisation d'un support exécutif à base de tâches. Pour se faire, une évaluation des deux derniers paradigmes est proposée via l'utilisation des supports Star PU et PaRSEC. Nous établissons leur intérêt pour l'évolution d'un code tel que cs comparé à l'utilisation de MPI+Open MP et nous les comparons sur des questions de complexité et de performance avec l'existant dans cs. Cette thèse est une collaboration entre le département MFEE de EDF R&D et l'équipe Tadaam qui met en commun Inria Bordeax Sud-Ouest, le LaBRI, l'Université de Bordeaux 1 et Bordeaux INP
The ever-evolving and ultimately complexification of computer architectures enforce scientific developpers to maintain both their software and technical knowledge in order to pertain their ability to fully exploit recent high performance computing platforms. Moreover the birth of heterogeneous plateforms and the large increase in core number per processor lead developpers of software parallelized through distributed paralleslim to see a small but consistent decrease of performance on recent plateforms. For cs, the transission from a full MPI code to a MPI+OpenMP code is not fully satisfying and cannot adress the scalability issue that alack of load balancing causes. This load balancing issue is a complex problemto solve for a code such as cs where partitionning is a NP problem and mustdata being indexed due to the unstructured nature of the meshes used. To solvethis kind of issue and many others (such as portability and heterogeneity), the use of a task based runtime recently won popularity through three different ways: loop parallelism, sequential submission of task(STF) or parameterized task graph (PTG). This thesis is focused on the transition effort of cs, a computational fluid dynamic (CFD) code with unstructured meshes (sparse linear algebra) from a mostly distributed parallelism to the integration of a task based runtime. In order to do this,we investigated the STF and PTG paradigms through the use of StarPU and PaRSECin cs. We demonstrate their advantages over a solution such that MPI + OpenMP as well ascompare their use in terms of complexity and performance with the existing solutions in cs. This thesis was made through the collaboration between the MFEE department of EDF R&D and the Tadaam research team, grouping Inria Bordeaux Sud-Ouest, the LaBRI,University of Bordeaux 1 and Bordeaux INP

La simulation des phénomènes physiques exige souvent l’utilisation d’une discrétisation du milieu sous forme de maillage. Un exemple de ce type de situation est la simulation d’écoulement de fluides et la simulation du stress géomécanique pour les gisements pétroliers. Dans ces cas, le milieu étudié n’est pas homogène et l'hypothèse sur l’homogénéité de ce milieu peut mener à des résultats incorrects. C’est pourquoi la simulation des hétérogénéités est très importante pour ce genre de problèmes.Cette thèse est consacrée à la simulation géostatistique des hétérogénéités sur les maillages non-structurés par les méthodes géostatistiques non-linéaires. Le but de cette thèse est la création d’algorithmes de simulation des hétérogénéités directement sur les maillages non-structurés, sans utiliser les maillages fins réguliers intermédiaires et de l’upscaling. On présente deux modèles théoriques pour les simulations des variables continues sur les maillages non-structurés qui sont les deux versions différentes du modèle Gaussien discret (DGM) - DGM 1 et DGM 2. Le modèle théorique utilisé dans cette thèse permet de convertir le problème de simulation sur un maillage non-structuré en un problème de simulation d’un vecteur Gaussien multivarié et l’application de fonctions de transformation adaptées pour chaque élément du vecteur. La simulation de faciès est aussi envisagée en utilisant une généralisation des modèles pluri-Gaussiens et Gaussien tronqués pour les maillages non-structurés.L’application des méthodes développées est illustrée sur un gisement pétrolier - le cas d’étude X (gisement du gaz offshore)
Simulations of physical phenomenon often require discretizing the medium with a mesh. An example of this type of simulation is the simulation of fluid flow through a porous medium and the evaluation of the geomechanical stress in the petroleum reservoir. The studied medium is often not homogeneous and applying a homogeneity hypothesis can lead to incorrect simulation results. That makes simulation of heterogeneities important for this kind of problems.This thesis is devoted to geostatistical simulations of heterogeneities on unstructured grids using methods of non-linear geostatistics. The objective of this work is the development of algorithms for simulating heterogeneities directly on unstructured grids without using intermediate fine scale regular grids and upscaling. We present two theoretical models for geostatistical simulations of continuous parameters on unstructured grids which are different generalizations of the Discrete Gaussian model (DGM) – DGM 1 and DGM 2. The proposed theoretical models enable converting the problem of geostatistical simulation on an unstructured grid into the well-studied problem of simulating multivariate Gaussian random vectors followed by application of block-dependent transformation functions. The problem of simulating facies is also addressed in this work, for which generalizations of pluri-Gaussian and truncated Gaussian simulation models for unstructured grids are proposed.An application of the proposed methods is demonstrated on a case study X, which is an offshore gas reservoir with a tartan-meshed grid

Le besoin croissant en simulation a conduit à l’élaboration de supercalculateurs complexes et d’un nombre croissant de logiciels hautement parallèles. Ces supercalculateurs requièrent un rendement énergétique et une puissance de calcul de plus en plus importants. Les récentes évolutions matérielles consistent à augmenter le nombre de noeuds de calcul et de coeurs par noeud. Certaines ressources n’évoluent cependant pas à la même vitesse. La multiplication des coeurs de calcul implique une diminution de la mémoire par coeur, plus de trafic de données, un protocole de cohérence plus coûteux et requiert d’avantage de parallélisme. De nombreuses applications et modèles actuels peinent ainsi à s’adapter à ces nouvelles tendances. En particulier, générer du parallélisme massif dans des méthodes d’éléments finis utilisant des maillages non structurés, et ce avec un nombre minimal de synchronisations et des charges de travail équilibrées, s’avèrent particulièrement difficile. Afin d’exploiter efficacement les multiples niveaux de parallélisme des architectures actuelles, différentes approches parallèles doivent être combinées. Cette thèse propose plusieurs contributions destinées à paralléliser les codes et les structures irrégulières de manière efficace. Nous avons développé une approche parallèle hybride par tâches à grain fin combinant les formes de parallélisme distribuée, partagée et vectorielle sur des structures irrégulières. Notre approche a été portée sur plusieurs applications industrielles développées par Dassault Aviation et a permis d’importants gains de performance à la fois sur les multicoeurs classiques ainsi que sur le récent Intel Xeon Phi
The growing need for numerical simulations results in larger and more complex computing centers and more HPC softwares. Actual HPC system architectures have an increasing requirement for energy efficiency and performance. Recent advances in hardware design result in an increasing number of nodes and an increasing number of cores per node. However, some resources do not scale at the same rate. The increasing number of cores and parallel units implies a lower memory per core, higher requirement for concurrency, higher coherency traffic, and higher cost for coherency protocol. Most of the applications and runtimes currently in use struggle to scale with the present trend. In the context of finite element methods, exposing massive parallelism on unstructured mesh computations with efficient load balancing and minimal synchronizations is challenging. To make efficient use of these architectures, several parallelization strategies have to be combined together to exploit the multiple levels of parallelism. This P.h.D. thesis proposes several contributions aimed at overpassing this limitation by addressing irregular codes and data structures in an efficient way. We developed a hybrid parallelization approach combining the distributed, shared, and vectorial forms of parallelism in a fine grain taskbased approach applied to irregular structures. Our approach has been ported to several industrial applications developed by Dassault Aviation and has led to important speedups using standard multicores and the Intel Xeon Phi manycore

La prédiction des concentrations d'espèces polluantes (suies et NOx) rejetées par les systèmes en combustion et l'évaluation de la durée de vie des parois de ce type d'installation passent par une bonne prise en compte des transferts radiatifs dans les outils de simulation de la combustion. Nous avons développé pour cela un code de calculs des transferts radiatifs basé sur la Méthode aux Ordonnées Discrètes et utilisant des maillages non structurés de la CFD. Le rayonnement des gaz de la combustion est pris en compte par un modèle statistique à bandes étroites en k-corrélés. Divers types de quadratures angulaires et trois schémas de dérivation spatiale différents ont été intégrés et sont comparés. Les tests de validations ont permis de montrer les limites aux fortes épaisseurs optiques de l'approximation de type volumes finis faites dans la méthode aux ordonnées discrètes. Les premiers calculs effectués sur des solutions obtenues par LES permettent de déterminer les termes sources radiatifs et les flux incidents aux parois instantanés, ce qui permet d'envisager le couplage avec la combustion, à plus long terme.