<HTML>
<HEAD>
<TITLE>Refinement in a box region</TITLE>
</HEAD>
<BODY>
<FONT FACE="Arial"><SPAN STYLE='font-size:11pt'>Good morning!<BR>
<BR>
I am meshing a canoe hull inside a box to use in an OpenFoam simulation and want to refine the water-air interface region. I first import an msh file created by importing an STL file, then reclassifying and re-saving the resulting new surface. These surfaces are then subtracted from the simulation box with the goal of meshing the region outside the hull. I looked at various posts and decided on using the Field functionality, but my approach seems to not be working. Can anyone point to an obvious problem?<BR>
<BR>
Jan Theron<BR>
<BR>
</SPAN><FONT SIZE="2"><SPAN STYLE='font-size:10pt'>-- <BR>
Dr. Johannes N. Theron<BR>
<BR>
Borsteler Chaussee 85-99A, Haus 9<BR>
22453 Hamburg<BR>
Germany<BR>
Tel: +49 176-877-57828<BR>
</SPAN></FONT><SPAN STYLE='font-size:11pt'><BR>
<HR ALIGN=CENTER SIZE="3" WIDTH="95%">Mesh.Algorithm3D = 1; //1=delaunay and 4=frontal<BR>
Mesh.Optimize=0;<BR>
Mesh.RemeshAlgorithm=1;<BR>
Geometry.Tolerance = 1.e-12;<BR>
Merge "4.8m-double-ender-rev4-decked-double-flipped-origin-in-center-refined1-reclassified.msh";<BR>
//RefineMesh;<BR>
CreateTopology;<BR>
<BR>
Surface Loop(300)={2,3,4};<BR>
<BR>
minX = -2.6;<BR>
maxX = 2.6;<BR>
minY = -1;<BR>
maxY = 1;<BR>
minZ = -.5;<BR>
maxZ = 0.6;<BR>
<BR>
IFc = 0.125; // Interface<BR>
WLt = 0.05; // Half thickness of interface layer<BR>
<BR>
pt = 0.05;<BR>
<BR>
Field[1] = Box;<BR>
Field[1].VIn = pt/10;<BR>
Field[1].VOut = pt*2;<BR>
Field[1].XMin = minX;<BR>
Field[1].XMax = maxX;<BR>
Field[1].YMin = minY;<BR>
Field[1].YMax = maxY;<BR>
Field[1].ZMin = IFc-WLt;<BR>
Field[1].ZMax = IFc+WLt;<BR>
<BR>
BackgroundField = pt;<BR>
<BR>
p1=newp; Point(p1) = {minX, minY, minZ, pt};<BR>
p2=newp; Point(p2) = {maxX, minY, minZ, pt};<BR>
p3=newp; Point(p3) = {maxX, maxY, minZ, pt};<BR>
p4=newp; Point(p4) = {minX, maxY, minZ, pt};<BR>
p5=newp; Point(p5) = {minX, maxY, maxZ, pt};<BR>
p6=newp; Point(p6) = {maxX, maxY, maxZ, pt};<BR>
p7=newp; Point(p7) = {maxX, minY, maxZ, pt};<BR>
p8=newp; Point(p8) = {minX, minY, maxZ, pt};<BR>
k2=newreg; Line(k2) = {p5, p8};<BR>
k3=newreg; Line(k3) = {p8, p1};<BR>
k4=newreg; Line(k4) = {p1, p4};<BR>
k5=newreg; Line(k5) = {p4, p5};<BR>
k6=newreg; Line(k6) = {p6, p7};<BR>
k7=newreg; Line(k7) = {p7, p2};<BR>
k8=newreg; Line(k8) = {p2, p3};<BR>
k9=newreg; Line(k9) = {p3, p6};<BR>
k10=newreg; Line(k10) = {p5, p6};<BR>
k11=newreg; Line(k11) = {p4, p3};<BR>
k12=newreg; Line(k12) = {p8, p7};<BR>
k13=newreg; Line(k13) = {p1, p2};<BR>
ll14=newreg; Line Loop(ll14) = {k2, k3, k4, k5};<BR>
Plane Surface(15) = {ll14};<BR>
ll16=newreg; Line Loop(ll16) = {k10, -k9, -k11, k5};<BR>
Plane Surface(17) = {ll16};<BR>
ll18=newreg; Line Loop(ll18) = {k12, k7, -k13, -k3};<BR>
Plane Surface(19) = {ll18};<BR>
ll20=newreg; Line Loop(ll20) = {k6, k7, k8, k9};<BR>
Plane Surface(21) = {ll20};<BR>
ll22=newreg; Line Loop(ll22) = {k11, -k8, -k13, k4};<BR>
Plane Surface(23) = {ll22};<BR>
ll24=newreg; Line Loop(ll24) = {k2, k12, -k6, -k10};<BR>
Plane Surface(25) = {ll24};<BR>
Surface Loop(26) = {15, 25, 19, 21, 23, 17};<BR>
Volume(600)={26,300};<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
</SPAN></FONT>
</BODY>
</HTML>