Hi Sebastian, You're looking for another way to achieve the same result, right? Can you tell about what you've tried already?
Any OpenFOAM issues are sort of out of scope for Gmsh but knowing what error messages you're getting, etc. would be very helpful. Sincerely, Max On Wed, Sep 25, 2019 at 12:31 PM Sebastian Tietjen <[email protected]> wrote: > Dear all, > > I try to create a mesh based on a hyperbolic-tangent function. This .geo > file works perfectly in gmsh, however, implementing this in OpenFoam causes > problems. My guess is that OpenFoam does not like how the layer command is > used. I try already couple of stuff to fix it, but it is somehow not > possible. > > I would be pleasd about any problem to fix my problem. > > Regards > > Sebastian Tietjen > > > > SetFactory("OpenCASCADE"); > //+ > Mesh.RecombineAll=1; > g2=1.5; > g3=1.5; > A3=4; > A1=2; > L2=1/A3; > L3=A3*L2; > L1=A1*L2; > N1=10;//128; > N2=16;//156; > N3=20;//312 > > //vertical direction - x3 (y) > one[0] = 1; > layer[0] = 0; > For i In {2:N3+1} > one[i] = 1; > X3=g3*(2*(i-1)/N3-1); > layer[i]=(L3/2)*(1+Tanh (X3)/Tanh (g3)); > EndFor > //Extrusion part in vertical direction > Point(1) = {0, 0, 0, 1.0}; > Extrude {0,L3,0} { > Point{1}; > Layers{one[], layer[]}; } > //horizontal direction - x2 (x) > oneX[0] = 1; > layerX[0] = 0; > For j In {2:N2+1} > oneX[j] = 1; > X2=g2*(2*(j-1)/N2-1); > layerX[j]=(L2*4/2)*(1+Tanh (X2)/Tanh (g2)); > EndFor > //Extrusion part in 1. vertical direction > Extrude {L2,0,0} { > Line{1}; > Layers{oneX[], layerX[]}; } > //+ > > //+Extrude surface > Extrude {0, 0, L1} { > Surface{1}; Layers{N1}; Recombine; > } > //+ > Transfinite Surface {1}; > //+ > Physical Surface("Back") = {7}; > //+ > Physical Surface("Front") = {1}; > //+ > Physical Surface("HotWall") = {5}; > //+ > Physical Surface("ColdWall") = {3}; > //+ > Physical Surface("Top") = {4}; > //+ > Physical Surface("Bottom") = {6}; > //+ > Physical Volume("Fluid") = {1}; > //+ > > _______________________________________________ > gmsh mailing list > [email protected] > http://onelab.info/mailman/listinfo/gmsh > -- Max Orok Contractor www.mevex.com
_______________________________________________ gmsh mailing list [email protected] http://onelab.info/mailman/listinfo/gmsh
