Detailed description of the projetcs

A. Mesh generation and adaptation.

A1. Mesh adaptation for an Arbitrary Lagrangian Eulerian (ALE) method
R. Abgrall, C. Dobrzynski, P.H. Maire (Bordeaux I)

Context
The physical model that will be considered throughtout this study is based on the equations of fluid dynamics. More specifically, we will be concerned by the simulation of multimaterial compressible fluid flows, such as those encountered in the domain of Inertial Confinement Fusion (ICF). These flows are characterized by strong shock waves and very high aspect ratios. For such numerical simulations, the Lagrangian framework, in which the mesh moves with the local fluid velocity is well adapted since this method deals with the interfaces in a natural manner. Unfortunately, when the vorticity and the shear become too strong this method can lead to non convex and tangled cells. This severe mesh distortion often implies a failure of the simulation.
Subject
One way to overcome this problem is to use an Arbitrary Lagrangian Eulerian (ALE) methodology in which the motion of the grid is chosen arbitrarily, see [4]. The strategy of ALE is to exploit this degree of freedom to improve both the accuracy and the robustness of the simulation. The main elements of most ALE algorithms are an explicit Lagrangian phase, a rezone phase in which the nodes of the Lagrangian grids are moved to improve the geometric quality of the grid and a remap phase in which the Lagrangian solution is transferred to the new grid. We have developed a new cell-centered second order Lagrangian scheme [1,5] which is the basement of our ALE methodology. Our rezoning step utilizes the local minimization of nodally based objective functions [3]. The remapping step is based on an unstructured extension of the simplified face based donor-cell method [6]. This ALE methodology enables to improve the robustness of our calculations. However, the rezoning step which is only based on pure geometrical criteria, tends to produce cells of equal area. This effect can decrease the accuracy brought by the Lagrangian phase. The aim of this study is to introduce a mesh adaptation algorithm [2] into the rezoning procedure in order to maintain the accuracy of the physical solution. The study will be done with quadrilateral meshes.

Contexte
Dans le cadre du projet Laser Méga joule (LMJ) le laboratoire CELIA (Université Bordeaux I) développe un programme de recherche concernant la fusion par confinement inertiel (FCI) en attaque directe. La simulation numérique ainsi que la restitution d'expérience de physique des plasmas est un des points clés de ce programme de recherche. Dans ce cadre, le CELIA a mis au point des outils numériques dédiés à  la conception et à  l'interprétation d'expériences de FCI. L'hydrodynamique de l'implosion d'une cible FCI est un processus très complexe. L'écoulement instationnaire multi-matériaux qui en résulte est le siège d'ondes de choc intenses au sein d'une géométrie caractérisée par de forts rapports d'aspect. La simulation de ces phénomènes consiste à  résoudre les équations de la dynamique des gaz en coordonnées de Lagrange. Le formalisme Lagrangien utilise un maillage lié à  la matière mais les opérateurs aux dérivées partielles sont exprimés en coordonnées Euleriennes. Cette écriture est bien adaptée à  la description de ce type d'écoulements. On s'affranchit ainsi du problème de la diffusion numérique inhérente aux termes de convection dans le formalisme Eulerien.
Sujet
Nous avons développé un nouveau schéma Lagrangien bidimensionnel pour la résolution des équations de la dynamique des gaz [1]. Ce schéma centré aux mailles est conservatif et entropique, les flux sont calculés à  partir d'un solveur aux noeuds qui peut s'interpréter comme une extension bidimensionnelle du solveur acoustique de Godunov. Comme toutes les méthodes Lagrangiennes, ce schéma souffre d'un manque de robustesse lorsque la vorticité et le cisaillement au sein de l'écoulement deviennent trop importants. La conséquence de ce manque de robustesse est l'apparition de mailles non-convexes voire croisées. Une stratégie possible pour remédier à  ce type de problème consiste à  régulariser le maillage afin d'améliorer la qualité géométrique des mailles. Les variables de l'écoulement sont ensuite projetées de l'ancien maillage vers le maillage régularisé. Cette méthode dénommée Arbitrary Lagrangian Eulerian (ALE) permet d'améliorer la robustesse des codes de calcul Lagrangiens. Toutefois, l'utilisation de cette méthode reste délicate car l'étape de régularisation, qui est basée sur la minimisation de fonctions de coà»t purement géométriques [3], ne prend pas suffisamment en compte la physique de l'écoulement.
Le but de ce projet est d'étudier une autre méthode de régularisation pour ces problèmes. Celle-ci sera basée sur une métrique issue du calcul de la solution physique (on utilisera, par exemple, le calcul de métrique décrit dans [2]) et tentera de bouger les noeuds du maillage en fonction de cette métrique. L'étude se fera principalement sur des maillages quadrangulaires.
References
[1] R. Abgrall, J. Breil, P.-H. Maire, J. Ovadia, A cel l centered Lagrangian scheme for two dimensional compressible flow problems http://hal.inria.fr/inria-00113542
[2] F. Alauzet, P.J. Frey, Estimateur d'erreur géométrique et métriques anisotropes pour l'adaptation de maillage. Partie 1 : aspects théoriques Rapport de recherche RR-4759, INRIA, 2003.
[3] P. M. Knupp, Achieving finite element mesh quality via optimization of the Jacobian matrix norm and associated quantities. Part I-a framework for the surface mesh optimization. Int. J. Numer. Meth. Engng 48,401-420 (2000).
[4] C.W. Hirt, A.A Amsden and J.L. Cook, An arbitrary Lagrangian-Eulerian computing method for all flow speeds. Journal of Computional Physics, (14), 227-253, 1974.
[5] P. H. Maire and J. Breil, A second order cell-centered Lagrangian scheme for two-dimensional compressible flows problems, Proceedinds of ICFD'07, Reading, UK, 2007
[6] L. G. Margolin, M. J. Shashkov, Second-order sign-preserving remapping on general grids, Journal of Computional Physics, (184) 266-298, 2003.

Back...
A2. Local Topological Modification of Hexahedral Meshes Using Dual-Based Operations.
T. Tautges (Argonne Natl. Lab.), Ph. Pébay (Sandia Natl. Lab.)
Local topological modification is widely used to improve mesh quality after automatic generation of tetrahedral and quadrilateral meshes. These same techniques are also used to support adaptive refinement of these meshes. In contrast, few methods are known for locally modifying the topology of hexahedral meshes. Most efforts to do this have been based on fixed transition templates or global refinement (ref-bern&eppstien). However, a dual-based approach to looking at this problem has also been studied, which so far has described known local hex mesh modifications using a set of only three fundamental operations (ref-usnccm). In this project, we will study the relation between the three dual-based operations described by Tautges and the insertion of a Boy surface described in Ref. 1. We will also investigate the algebra of these operations and attempt to construct a group structure for the set. Depending on progress, we may also pursue the application of these operations to real hex mesh examples.
1. Marshall Bern, David Eppstein, Jeff Erickson, "Flipping Cubical Meshes", ACM Computer Science Archive, http://www.arXiv.org/cs.CG/0108020, June 20, 2002.
2. Timothy J. Tautges, "Local Topological Modification of Hexahedral Meshes Using Dual-Based Operations", 8th U.S. National Conference on Computational Mechanics, Austin, TX, July 2005.
Back...
A3. Anisotropic mesh adaptation and levelsets for Stokes simulations.
P. Frey (LJLL, UPMC), F. ALauzet (INRIA-Rocquencourt)
Curvature-based adaptation seems a promising way of improving the accuracy of numerical simulations in the vicinity of the interface between two viscous fluids. I this project, we will attempt to perform numerical simulations on viscous test cases modelled by Stokes equations using a coupling between mesh adaptation and levelsets techniques. This work is supported and funded by an INRIA ARC grant. Back...
A4. DEFORMA: High deformation and mesh reconnection.
B. Després, S. Del Pino (CEA-DAM DIF)
High deformation flows may generate pathological meshes. This is particularly true in the Lagrangian context. So one needs specific algorithms, to reconnect the mesh lines for example. This project is twofolds: developp original algorithms; test them for a model problem with a 2D lagragian test code that will be provided. The test problem is the fall of a water bubble in a glass filled with air.
DEFORMA: Grandes déformations et reconnection de maillage
Les calculs d'écoulements fluides en grande déformation produisent des configurations de maillage facilement pathologiques. On se place dans le contexte lagrangien ou le maillage suit la matière. Cela nécessite des algorithmes spécifiques en particulier pour gerer la reconnection des lignes de maillage. Le projet se fixe comme objectif d'étudier un problème modèle de chute d'une bulle d'eau dans un verre rempli d'air, et d'apporter des solutions algorithmiques originales. Un code d'étude lagrangien 2D servira comme point de départ des études numériques.
Back...
A5. Geometry correction
G. Alleon (EADS, Toulouse)
Industrial CAD models often contains patterns which make them difficult to use for numerical computation (i.e. difficult to mesh). The goal of this project is to provide jCAE based tools to remove some of such elements. The available software include:
  • Open CASCADE: A C++ geometry engine which is able to work with real industrial CAD.
  • occjava: The part of jCAE which allow access to Open CASCADE in Java The geometrical elements which have to be removed are the following (by priority order):
    1. Through hole and non-through hole removal
      • Detection algorithm
      • suppression algorithm
    2. Bevel and fillet removal
      • Detection algorithm
      • suppression algorithm to be announced
    The create tool will be tested on models available into the OpenCASCADE shape gallery.

    Back...
A6. Co-refinement of fault surfaces
P. Havé, A. Benali, L. Agelas(IFP)

Context:
Modern geologic simulations on complex geometries require specific modeling tools for the accurate representations of heterogeneous components. The associated meshes evolve according to a displacement field provided by a geologist, which must preserve the initial heterogeneity distribution. The displacements are defined a priori at a few time steps, and then continuously interpolated between two given positions. Moreover, topological modifications may appear such as fractures, defining geologic faults. Computing transport of information through a fault is one of the goal of this project. The co-refinement (common refinement) of two non-coincident mesh surfaces aims at providing a mesh which likns the refined surfaces by mean of a continuous bijection. This co-refinement must allow various numerical strategies through the faults, either relying on volumic mesh adaptation or specically designed numerical schemes in the context of extended stencil finite volume methods. The first part of this project consists in defining a co-refinement algorithm which satisfies some properties ensuring the stability of the scheme. The second part consists in developping a prototype to test the algorithms on realistic non-coincident meshes (with large gaps, overlaps, non planar surfaces...). The results of this project are meant to be used in industrial geological codes.

Contexte:
La simulation géologique moderne en géométries complexes nécessite des outils spécifiques pour une représentation précise des composants hétérogènes. Les maillages associés évoluent suivant un champ de déplacement fourni par un géologue, qui doit préserver la distribution initial des hérogénéités. Le déplacement est défini a priori en quelques pas de temps, et est ainsi interpolé continuement entre deux positions connues. De plus, des modifications topologiques peuvent apparaître comme des fractures, définissant ainsi des failles géologiques. Le calcul du transport d'information à  travers les failles est un des buts de ce projet. Le co-raffinement de deux maillages de surfaces non-coïncidents vise à fournir un maillage associant les surfaces par le biais d'une bijection continue. Le co-raffinement doit permettre l'expression de stratégies numériques variées à  travers les failles, soit en s'appuyant sur une adaptation du maillage volumique, soit par des schémas numériques spécifiquement adaptatés dans le contexte des méthodes de volumes finis à  support élargi. La première partie de ce projet consiste en la définition d'un algorithme de co-raffinement respectant certaines propriétés assurant la stabilité du schéma numérique. La seconde partie consiste au développement d'un prototype afin de tester les algorithmes sur des maillages non-coïncidents réalistes (avec des trous, des superpositions, des surfaces non planes...). Les résultats de ce projet seront suceptibles d'être utilisés dans des codes géologiques industriels. Back...
A7. Quad mesh
G. Alleon (EADS)
This projects aims at creating a 3D surfacic quad mesher using jCAE technologies. The available software include several tools.
  • Open CASCADE: A C++ geometry engine which is able to work with real industrial CAD
  • occjava: The part of jCAE which allow access to Open CASCADE in Java
  • amibe: The 3D surface mesher included in jCAE
For the quad mesher two approaches could be considered:
  • Advancing-front: This algorithm will start with a 1D discretisation of the edges of the surface to mesh. It will use occjava to query the geometry description
  • Pairing triangles: This algorithm will start with a triangle mesh created by amibe. It will pair triangles to create quad. It will use occjava to query the curvature if the geometry. The mesher will try to create best quality elements, according to common criterias (taper, skew angle and warping angle).
Back...
A7. Maillage volumique - couches limites
L. Francez, J.F. Lague (Distène)
Un besoin industriel est la réalisation automatique de connexions de maillages provenant de différentes méthodes (Octree, tétraédrisation automatique). L'objectif de ce projet est de proposer des algorithmes efficaces afin de générer la zone de raccord entre maillages volumiques de natures différentes. Sujet: le premier travail consistera en la génération de la zone de raccord (ou couche tampon). Les difficultés résident dans la gestion des frontières à forte courbure ou de forte concavité Plusieurs étapes sont donc suggérées.
consultez la description complète. du projet (PDF)
Back...

B. Innovative numerical methods.

B1. Guaranteed, locally efficient, asymptotically exact, and asymptotically robust a posteriori error estimation in the finite element method.
M. Vohralík (LJLL, UPMC)
The theory of a posteriori error estimation in numerical approximation of partial differential equations aims at controlling the error between the numerical approximation and the (unknown) exact solution. The ultimate goal is quadruple:
  1. give an upper bound on the error in the numerical solution, which only uses the approximate solution and can thus be evaluated (guaranteed upper bound);
  2. give an expression for the estimated error locally, for example in each element of the computational mesh, and assure that this estimate on the error represents a lower bound for the actual error, up to a multiplicative constant (local efficiency);
  3. assure that the ratio of the estimated error and the actual error goes to one (asymptotic exactness);
  4. guarantee the three previous properties independently of the data or mesh properties (robustness).
It is really amazing that one may be able to give a bound on the error without the knowledge of the exact solution and it is indeed possible. Three main branches of a posteriori estimates in the finite element method have evolved during the last decades, see the survey books by Ainsworth and Oden and Verfà¼rth. Explicit residual-based estimators are build upon the fact that approximate solution does not satisfy the given partial differential equation, are computationally inexpensive, and usually fulfill the first and second desirable property. However, up to very rare exceptions, the first property is not really satisfied in the strict sense, since one has a computable upper bound up to an unknown multiplicative constant. As such, it is completely sufficient for the usual practice, where one only uses it in order to refine mesh elements with increased error, but not for the actual control of the error. Implicit residual-based estimators like the equilibrated residual method remove much of the above drawbacks, but are computationally expensive and only give the guaranteed upper bound asymptotically. Finally, averaging-based estimators as the celebrated Zienkiewicz-Zhu one are easy to compute, often fulfill the third property, but systematically fail the first one in the strict (guaranteed) sense.
We have recently introduced computationally inexpensive estimators which fulfill all the four desirable properties (the robustness may however only be asymptotic), in the framework of mixed finite element and finite volume methods. The purpose of this project is to adapt these techniques to the finite element case. At first, a model second-order elliptic problem will be considered, featuring however discontinuous and anisotropic coefficients. Extensions to (singularly-perturbed) convection-diffusion-reaction problems are intended later. In a theoretical part, rigorous analysis of the estimator properties will be done. Implementation and extensive numerical simulations will form an important second part the project, with the purpose to verify the properties of the developed estimator and compare it to some well-known ones, like the Zienkiewicz-Zhu, explicit and implicit residual, or some estimators introduced recently.
see also a PDF description of this project.
Back...
B2. Compréhension et mise en place d'un algorithme d'optimisation
A. Rousseau (INRIA-Rhône Alpes)
In a collaboration with the french national agency for development of ecology and energy control (ADEME), we intend to build a new numerical method to compute small scale phenomena without refining our mesh ion atmospheric models. In an existing mesh, we virtually drop some particles that are moved thanks to a system of SDEs (stochastic differential equations) written by S. Pope. We then compute local values of the required fields thanks to averaging them on a cloud of particles. We are thus using Monte-Carlo methods, and aim to study their convergence rate, and finally compare them to classical refinement methods. One (contastable) constraint of Pope's model is the uniform value of the density (rho=cst). That is, the particles have to be uniformly distributed at every time step. This particular problem is the framework of our CEMRACS project. We aim to use D. Bertsekas "Auction Algorithm" in order to move a given cloud of particles to a new position, which is also given in advance, and that realizes the constraint rho=cst. Naturally, the transport cost will have to be minimum. This is a problem of 3D optimal transport, which is known to be difficult.

Dans le cadre d'un projet avec l'ADEME, nous tâchons de mettre en place une méthode numérique qui permette de calculer des phénomènes de sous-maille sans avoir besoin de raffiner le maillage. Typiquement, dans un maillage existant, on "jette" tout un tas de particules, soumises à  un système d'équations différentielles stochastiques (Pope & al.). Ensuite, on récupère des informations locales en faisant des moyennes sur les particules. Ce sont donc des méthodes de Monte-Carlo qui, on l'espère, seront comparables ou même meilleures, à  des méthodes plus classiques de raffinement de maillage... Une des contraintes (discutable) du modèle de Pope est d'avoir un champ uniforme (rho=cst), et donc des particules uniformément réparties : c'est là -dessus que porte le projet. Le travail consiste à  mettre en place l'Auction Algorithm développé par Bertsekas (US), qui permettra de déplacer les particules dans notre domaine de calcul afin que leur position finale soit uniforme, et en ayant minimisé le coût de transport. Autrement dit, c'est un problème numérique de transport optimal 3D.
Back...
B3. Engineering and flow in coronary arteries
J.F.Gerbeau, I. Vignon-Clémentel (INRIA-Rocquencourt)
The objective is to develop a mathematical modeling and numerical tools to simulate cardiac activity. The main purpose is to manage to model blood flow in the largest coronary arteries which perfuses the myocardium, taking into account its interaction with the myocardium movement. This project will first consist in gathering morphometric coronary data and generating one-dimensional coronary trees on the external surface of the heart, which will lead to a three-dimensional geometry of such trees. Simulations of blood flow in these coronary arteries will then be based on the mesh of this three-dimensional geometry, its displacement due to the myocardium movement and reduced models for boundary conditions, using the ALE code developed at INRIA in the REO project.

L'objectif est de développer des modèles mathématiques et des outils de simulation numérique de l'activité cardiaque. Le but est ici de parvenir à  une simulation de l'écoulement du sang dans les plus grosses artères coronaires pour assurer la perfusion du myocarde en prenant en compte son interaction avec la paroi du myocarde en mouvement.
Après avoir rassemblé des données sur la morphologie des coronaires, il s'agira de générer des arbres coronaires unidimensionnels sur la surface externe d'un coeur. A partir de cette ossature unidimensionnelle, sera construite une géométrie tridimensionnelle. Le maillage de cette géométrie, son déplacement provenant du mouvement du coeur et l'utilisation de modèles réduits en entrée et sortie permettront la simulation de l'écoulement dans ces artères coronaires avec le code de calcul ALE développé dans le projet REO (INRIA).
Back...
B4. Direct simulation of blood flow through a stent
M. Fernandez, J.F.Gerbeau (INRIA-Rocquencourt)
The objective is to model a new technique to clinically treat cerebral aneurysms which consists in placing a stent at the entrance of the aneurysm. The mesh of the stent is sufficiently fine to significantly reduce blood flow in the aneurysm and trigger clotting which naturally acts as a plug. A homogenized mathematical model of such a stent has already been proposed and tested in three-dimensional numerical simulations of blood flow in cerebral aneurysms. Questions that could be answered were at the "macroscopic" level: influence of the stent on blood flow in the aneurysm and flow rates in the lateral branches. Animal experiments that were conducted at INRIA confirmed qualitatively the results. They moreover revealed an interesting phenomenon: the fine structure of the stent could "beneficially" influence blood flow. The scope of the project is to work on this phenomenon. More precisely, the geometry of the stent will be modeled, taking into account the threads, and direct simulations of flow through such a fine geometry will be performed. If possible, industrial definitions of the stent will be used. These simulations will be compared to the above mentioned homogenized model.

On cherche à  modéliser une technique nouvelle de traitement d'anévrismes cérébraux développée récemment. La méthode consiste à  placer un dispositif, appelé stent, à  l'entrée de l'anévrisme. La finesse de la grille du stent entraà®ne une dimininution suffisante du flux dans l'anévrisme pour permettre un colmatage naturel.
Un modèle de stent a déjà  été proposé et implémenté dans des modèles numériques 3D d'écoulement dans des anévrismes cérébraux. Les questions auxquelles ce modèle a permis de répondre étaient d'ordre "macroscopique": influence du stent sur l'écoulement dans l'anévrisme et sur le débit dans les branches latérales. Des expériences animales à  l'INRA ont confirmé qualitativement les résultats obtenus et ont mis en évidence un phénomène inattendu: la structure fine du stent aurait une influence «bénéfique» sur l'écoulement. L'objectif de ce stage est de travailler sur cette question. Plus précisément, nous modéliserons la géométrie du stent, en prenant en compte ses fils, et étudierons un écoulement transverse par simulation directe. Si possible, nous tenterons de travailler directement sur les données du programme de Cardiatis définissant le tressage du stent. Ces simulations seront confrontées aux modèles homogénéisés déjà  utilisés.
Back...
B5. Image-based meshing of aneurims for flow simulations. Predictive factors of the rupture risk.
M. Thiriet (INRIA-Rocquencourt)
This projet has been CANCELLED.
Back...
B6. Adaptation non conforme et régularisation de maillage en dynamique des fluides Lagrangienne.
Ph. Hoch (CEA-DAM DIF)
Le contexte est celui de la dynamique des fluides compressibles bidimensionnelles écrite en variables Lagrangiennes. Les équations sur les vitesses et l'énergie seront résolues par une méthode Lagrangienne donnée (centrée Després/Mazeran ou Maire/Breil ou bien décalée Wilkins).
On considère un écoulement laminaire en géométrie plane avec un maillage initial cartésien. Le bord gauche comprime le domaine complet avec une pression linéaire qui décroit de P2 en bas à  P1 en haut tandis qu'une condition de mur est imposée sur le bord droit. L'approche classique pour obtenir un maillage avec des rapports d'aspects raisonnables passe par une régularisation. On parle dans ce cas d'une méthode ALE non-directe qui comprend une première phase de lissage de la grille puis d'une seconde phase de recalcul des quantités conservatives sur ce nouveau maillage. Malheureusement, ce point de vue ne permet d'avoir des éléments d'aire raisonnable que jusqu'à  un certain moment.
Le but de cette étude est d'abord d'étudier une stratégie de déraffinement local non conforme simplifiée (dite par soudure) de manière à conserver des éléments de type quadrangulaires après le raffinement. Plus précisemment, la ligne de soudure introduite devra remplacer les arêtes des éléments pour lesquels les rapports d'aspect sont mauvais. Dans l'implémentation on ne devra donc calculer que l'intersection d'une ligne (la ligne de soudure) avec quelques éléments (pour avoir une reconstruction conservative). Dans un deuxième temps, on étudiera également le problème symétrique du raffinement local en considérant un problème en détente (le calcul géométrique pour avoir des reconstructions conservatives étant plus simple pour le raffinement).

Voir aussi une description PDF du projet.
Back...

C. Scientific Visualisation.

C1. Telemesh : a "Web 2.0" GUI for 2D scientific visualisation and mesh edition.
J.D. Benamou (INRIA-Rocquencourt, Rice Univ.)
.
Modern standardized WEB technologies, broadband network and faster computers are today's keys to fast and painless performant WEB graphic GUI's development. SVG (Scalable Vector Graphic) is a scriptable dynamic markup langage for the WEB. Soon a standard (it is natively supported by Firefox, Opera and Safari, with a pugin in IE) it allows to programm browsable interactive GUI. This new techniques have yet had few or no impact in the scientific computing area. The first part of the proposed project consists in programming a client-server application which enables the user to
  1. Upload a 2-D mesh (and possibly an associated EDP solution) to the server where it will be converted on the fly to SVG and served back to the client (brwowser) for visualization.
  2. Edit the SVG mesh (move/remove/add nodes) inside the browser box.
  3. Convert and download the new mesh in its original format after a new upload/conversion on the server.
The second part of the project would naturally consist in the definition and coding of a "WEB service" API to automate the call of Telemesh functions from other servers. It will then be possible for anyone to give access to scientific simulations codes on a WEB server and take advantage of Telemesh to provide visualization and edition simply through the automatic generation of links to Telemesh server.

Les technologies modernes du WEB et leur standardisation, la puissance des ordinateurs et l'avènement du haut débit permettent aujourd'hui de développer rapidement et facilement des interfaces graphiques performantes sous forme de pages WEB chargées dans un navigateur.
Le langage SVG (Scalable Vector Grahics) en particulier est en passe de devenir un standard, il est supporté de façon native par FireFox, Opera et Safari, avec un plugin par Internet Explorer). Il permet de réaliser des pages Internet graphiques dynamiques et interactives: des interfaces graphiques fonctionnant dans un navigateur. Ces nouvelles techniques ont peu ou pas été utilisées dans le domaine de la simulation numérique d'EDP. Le projet proposé consiste dans un premier temps à  programmer sous une architecture client-serveur un logiciel permettant de
  1. Télécharger (upload) un maillage 2-D (et éventuellement une solution d'EDP associée) vers le serveur afin de le (les) visualiser en retour au format SVG dans un navigateur.
  2. Manipuler coté client le maillage (déplacement/supression/rajout) de noeuds (toujours dans un navigateur).
  3. Récupérer le nouveau maillage dans son format d'origine et/ou lui appliquer un traitement après un échange automatisé avec le serveur.
La suite logique de ce projet consisterait à  définir et programmer une API "Web service" permettant d'automatiser l'appel des fonctions de visualisation à  partir d'autres sites WEB. Il deviendrait alors possible de mettre à  disposition un code de simulation sur un serveur Internet à  travers une interface WEB et d'assurer la visualisation des résultats via la génération (automatique) d'un simple lien vers le serveur de TELEMESH.
Back...
C2. Large mesh decimation and visualisation
G. Alleon (EADS, Toulouse)
CAE include a still very prototypical large mesh viewer. This viewer is based on an out-of-core octree [1]. One limitation is that it currently rely on the geometry for far viewing. In many case the geometry used to do the mesh is not available. The goal of this project is to create a multi-level representation of the mesh for far and mid distance view, using decimation. jCAE already contains a decimation algorithm based on Michael Garland work, but it's currently not design for large mesh (in core implementation). The create mesh viewer should be implemented using VTK Java interface.
[1] P. Cignoni, C. Montani, C. Rocchini, R. Scopigno, External Memory Management and Simplification of Huge Meshes
[2] Michael Garland, Quadric Error Metrics

Back...