Vof2d_compute_radius

Bonjour,
Je souhaite utilser cette routine pour calculer par la suite une surface à partir d’une fraction de phase. Quelle est la définition exacte de radius dans cette subroutine?
Merci par avance
Cédric

Cédric,

vof2d_compute_radius(normalx, normaly, cx, cy, vof, vof_case, radius)

Dans cette routine, radius correspond à la distance (au sens de la projection orthogonale) entre l’interface et coin inférieur gauche de la cellule.

Tous les arguments de la routine sont définis dans la configuration de référence, ce qui signifie qu’il faut appliquer des transformations à la normale et aux dimensions de la cellule (cx,cy) avant de l’utiliser. Cette opération s’effectue en utilisant la routine vof2d_get_reference.

La routine retournera radius dans la configuration de référence. Pour retrouver cette distance dans la configuration réelle, il faut encore appliquer des transformations. Pour la cellule (i,j), ces transformations s’écrivent :

select case (vof_case)
case(2)
   radius = radius + normalx*dx(i)
case(3,6,7)
   radius = -radius
case(4)
   radius = radius + normaly*dy(j)
end select

Tu peux trouver un exemple de cette procédure dans le fichier src/lib/discretization/multiphase/vof/plic/vof2d.f90 entre les lignes 217 et 235.

Antoine

Merci Antoine! C’est parfait pour moi!
Cédric