
<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi all,<div class=""><br class=""></div><div class="">I’m trying to create a disk mesh with a given background and, within a rim around the center, a local decrease of the characteristic length that follows a Gaussian (see attached image 1).</div><div class=""><br class=""></div><div class="">Code:</div><div class=""><br class=""></div><div class=""><div class="">SetFactory("OpenCASCADE");</div><div class="">Disk(1) = {0, 0, 0, 4};</div><div class="">Physical Surface("surfacedomain") = {1};</div><div class=""><br class=""></div><div class="">Field[3] = MathEval;</div><div class="">Field[3].F = ".5 - (.5-0.01)*exp(-(1./0.1) *((sqrt(x*x+y*y)-1))*((sqrt(x*x+y*y)-1)))";</div><div class=""><br class=""></div><div class="">Background Field = 3;</div></div><div class=""><br class=""></div><div class="">However, for the parameters above the rim is not closed, what am I getting wrong?</div><div class="">If I change parameters slightly, e.g. 0.1 to 0.2 or 0.01 to 0.005 I get the correct representation. My problem is, that I eventually would like to do this in 3D within a sphere instead of a disk. However, this is visually harder to inspect (or is there a good way to do this?), so I would like to make sure that things work in principle in 2D.</div><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Lars</div><div class=""><br class=""></div><div class=""><img apple-inline="yes" id="E33AAF69-83F6-4B60-8752-FB58884C4160" width="662" height="583" src="cid:55F7C697-63BB-466C-8778-8321B96825EA" class=""></div></body></html>