<div dir="ltr"><div>Hi Sebastian, <br></div><div><br></div><div>You're looking for another way to achieve the same result, right? <br></div><div>Can you tell about what you've tried already? </div><div><br></div><div>Any OpenFOAM issues are sort of out of scope for Gmsh but knowing what error messages you're getting, etc. would be very helpful.   <br></div><div><br></div><div>Sincerely, <br></div><div>Max <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Sep 25, 2019 at 12:31 PM Sebastian Tietjen <<a href="mailto:Sebastian_Tietjen@web.de">Sebastian_Tietjen@web.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear all,</div><div><br></div><div>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. <br></div><div><br></div><div>I would be pleasd about any problem to fix my problem.</div><div><br></div><div>Regards</div><div><br></div><div>Sebastian Tietjen  <br></div><div><br></div><div><br></div><div><br></div><div>SetFactory("OpenCASCADE");<br>//+<br>Mesh.RecombineAll=1;<br>g2=1.5;<br>g3=1.5;<br>A3=4;<br>A1=2;<br>L2=1/A3;<br>L3=A3*L2;<br>L1=A1*L2;<br>N1=10;//128;<br>N2=16;//156;<br>N3=20;//312<br><br>//vertical direction - x3 (y)<br>one[0] = 1; <br>layer[0] = 0; <br>For i In {2:N3+1}<br>    one[i] = 1;   <br>       X3=g3*(2*(i-1)/N3-1);<br> layer[i]=(L3/2)*(1+Tanh (X3)/Tanh (g3));<br>EndFor <br>//Extrusion part in vertical direction<br>Point(1) = {0, 0, 0, 1.0};<br>Extrude {0,L3,0} {    <br>Point{1}; <br>Layers{one[], layer[]}; }<br>//horizontal direction - x2 (x)<br>oneX[0] = 1; <br>layerX[0] = 0; <br>For j In {2:N2+1}   <br>        oneX[j] = 1;   <br>      X2=g2*(2*(j-1)/N2-1);<br> layerX[j]=(L2*4/2)*(1+Tanh (X2)/Tanh (g2));<br>EndFor <br>//Extrusion part in 1. vertical direction<br>Extrude {L2,0,0} {    <br>Line{1}; <br>Layers{oneX[], layerX[]}; }<br>//+<br><br>//+Extrude surface<br>Extrude {0, 0, L1} {<br>  Surface{1}; Layers{N1}; Recombine;<br>}<br>//+<br>Transfinite Surface {1};<br>//+<br>Physical Surface("Back") = {7};<br>//+<br>Physical Surface("Front") = {1};<br>//+<br>Physical Surface("HotWall") = {5};<br>//+<br>Physical Surface("ColdWall") = {3};<br>//+<br>Physical Surface("Top") = {4};<br>//+<br>Physical Surface("Bottom") = {6};<br>//+<br>Physical Volume("Fluid") = {1};<br>//+<br><br></div></div>
_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
<a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">http://onelab.info/mailman/listinfo/gmsh</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Max Orok</div><div>Contractor<br></div></div><div dir="ltr"><a href="http://www.mevex.com" target="_blank">www.mevex.com</a><br><br><img src="https://docs.google.com/uc?export=download&id=1fHTIiW4OMUjQr1iOkspQ7wiEsxunoOs0&revid=0B6x5w-5zVaEjSkpwbm5oY29jbG1XMzJoYldXTmJpNGFtb3dVPQ" width="164" height="42"><br></div></div></div></div></div></div>