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_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
