Cassiopee Forum
CFD python modules
S'il vous plaît connectez vous ou Enregistrement'.
Date et heure en ce moment : dimanche, 19 novembre 2017, 00:27

Nouvelle discussion | Répondre
Cassiopee Forum :: Discussions, suggestions, bug report :: Post :: formulation extractPoint  ::
Sylvain Mouton
Junior Member
Image


Messages: 66
formulation extractPoint (mardi, 18 octobre 2011, 15:16) citation  
Je m'interroge sur le fonctionnement de la fonction extractPoint dans le cas où l'on traite un maillage aux noeuds et des valeurs solution aux centres (cas le plus fréquent avec elsA). Il me semble en regardant les sources que extractPoint commence par reconstituer un maillage en centres avant d´appliquer l'extraction sur ce nouveau maillage.
Cela ne me paraît pas correct car sur les bords des blocs, un point initialement dans le volume décrit par le maillage aux noeuds peut se trouver hors du maillage en centres, et donc impossible à interpoler.
Je préfèrerais nettement une interpolation d'ordre zéro, dans laquelle la fonction retourne simplement la valeur de la solution contenue dans la cellule où se trouve le point. Cette approche serait cohérente avec la formulation volume fini retenue.
On pourrait ensuite imaginer des interpolations d'ordre un ou plus qui tiendrait compte des valeurs dans les cellules voisines, au prix d'un certain lissage de la solution.

Ip enregistré Statut: déconnecté Profil | Site Web 
Ordre des réponses: Première réponse en dernier :: Première réponse en premier
Sylvain Mouton
Junior Member
Image


Messages: 66
RE: formulation extractPoint (mercredi, 30 novembre 2011, 15:57) citation  
Evidemment, il est possible pour contourner ce comportement d'interpoler les valeurs aux centres sur les noeuds du maillage (center2Node) puis de faire extractPoint sur ce nouveau résultat ou tout est en centre. Avec cette solution, plus de "trous" entre les blocks du maillage, mais par contre toujours le problème du lissage de la solution.

Ip enregistré Statut: déconnecté Profil | Site Web 
Sylvain Mouton
Junior Member
Image


Messages: 66
RE: formulation extractPoint (lundi, 5 décembre 2011, 11:40) citation  
Et sur un de mes cas, si j'extrait le champ au point en question après avoir fait center2Node, j'obtiens un Mach isentropique de 0.0725 (pour une itération donnée, au cours du calcul). Sans center2Node, j'obtiens 0.0906 (cette dernière valeur est en accord avec tecplot).
Ceci sont les valeurs obtenues lors d'une exécution en parallèle (le proc qui extrait ne dispose que de quelques blocks dumpés par elsAxdt dans l'output-tree, auquel je rajoute à la volée les coordonnées du maillage). Lorsque je traite le champ en entier en séquentiel à la fin du calcul (en lisant le fichier résultat), les deux valeurs sont à nouveau à peu près identiques et correctes...
Bizarre, non ?

Ip enregistré Statut: déconnecté Profil | Site Web 
Nouvelle discussion | Répondre

Software PBLang 4.65 © 2002-2003 by Martin Senftleben
Image