[Gmsh] boundaries

[ouafa bouloumou]

bonjour,

 

j'essaie de générer un maillage  en vue d'une utilisation sous OpenFoam. pour cela, j'ai défini un ensemble de conditions aux limites grace aux surfaces physiques.

l'ennui c'est que quand je convertie mon maillage sous format OpenFoam, certaines parois ne contiennent aucune faces.

ci joint un aperçu de mon script:

 

Physical Surface("inlet patch") = {fg_canal[21],fg_canal[22],fg_canal[23],    //Amont
                             ns_tube[1],ns_tube[2],ns_tube[3],ns_tube[4],  //Tube
                             ns_tube[5],ns_tube[6],ns_tube[7],ns_tube[8],ns_tube[9]};
                             
Physical Surface("outlet patch") = {fd_canal[24],fd_canal[25],fd_canal[26]};   //Aval

Physical Surface("interface patch") = {s10:s20,  // solide_fluide dans canal
                                 fg_solide[7],fb_solide[8],fd_solide[9],fh_solide[6]}; // solide_fluide dans solide 
                                 
Physical Surface("jonction empty") = {s1:s9, // trou canal
                                      ns_solide[1],ns_solide[2],ns_solide[3],ns_solide[4], // trou solide
                                      ns_solide[5],ns_solide[6],ns_solide[7],ns_solide[8],ns_solide[9]};  

Physical Surface("wall") = {ns_canal[1],ns_canal[2],ns_canal[3],ns_canal[4],ns_canal[5],ns_canal[6], //Top
                      ns_canal[7],ns_canal[8],ns_canal[9],ns_canal[10],ns_canal[11],ns_canal[12],
                      ns_canal[13],ns_canal[14],ns_canal[15],ns_canal[16],ns_canal[17],ns_canal[18],
                      ns_canal[19],ns_canal[20],ns_canal[21],ns_canal[22],ns_canal[23],ns_canal[24],
                      ns_canal[25],ns_canal[26],
                      s21:s26,               //Bottom_canal
                      ns_solide[10],ns_solide[11],ns_solide[12],ns_solide[13],   //Bottom_solide
                      ns_solide[14],ns_solide[15],ns_solide[16],ns_solide[17],
                      ns_solide[18],ns_solide[19],ns_solide[20],
                      fb_canal[26],fb_canal[20],fb_canal[15],fb_canal[11],fb_canal[16],fb_canal[23],   //Front surface
                      fb_solide[20],fb_solide[15],fb_solide[11],fb_solide[16],
                      fh_canal[24],fh_canal[18],fh_canal[14],fh_canal[13],fh_canal[17],fh_canal[21],   //Back surface
                      fh_solide[18],fh_solide[14],fh_solide[13],fh_solide[17],
                      fg_solide[18],fg_solide[19],fg_solide[20],fd_solide[17],fd_solide[12],fd_solide[16],  //bords gauche et droite du solide
                      fg_tube[7],fb_tube[8],fd_tube[9],fh_tube[6]};  // contour cylindre

dans le fichier boundary j'obtiens:

 

 

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

6
(
    jonctionempty
    {
        type            patch;
        nFaces          0;
        startFace       833457;
    }
    interfacepatch
    {
        type            patch;
        nFaces          0;
        startFace       833457;
    }
    wall
    {
        type            patch;
        nFaces          30375;
        startFace       833457;
    }
    inletpatch
    {
        type            patch;
        nFaces          2565;
        startFace       863832;
    }
    outletpatch
    {
        type            patch;
        nFaces          0;
        startFace       866397;
    }
    defaultFaces
    {
        type            patch;
        nFaces          1824;
        startFace       866397;
    }
)

// ************************************************************************* //


 

quelqu'un aurait il une idée? merci à tous,

Ouafa.

 

_________________________________________________________________
Inédit ! Des Emoticônes Déjantées! Installez les dans votre Messenger !
http://www.ilovemessenger.fr/Emoticones/EmoticonesDejantees.aspx
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090708/a5e68078/attachment.html>