From cf14b08c1e25db25a6036649f868524280863c43 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@uliege.be>
Date: Thu, 12 May 2022 13:44:35 +0200
Subject: [PATCH] refine close to feed
---
Antennas/dipole3d.geo | 17 ++++++++++++++---
Antennas/dipole_data.geo | 2 +-
2 files changed, 15 insertions(+), 4 deletions(-)
diff --git a/Antennas/dipole3d.geo b/Antennas/dipole3d.geo
index c2f94f2..d75a7e7 100644
--- a/Antennas/dipole3d.geo
+++ b/Antennas/dipole3d.geo
@@ -1,7 +1,7 @@
-fac = 2 ;
+fac = 5 ;
-Mesh.CharacteristicLengthFactor = fac ;
-Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)
+Mesh.MeshSizeFactor = fac ;
+Mesh.Algorithm3D = 10; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)
Mesh.Optimize = 1 ;
// characteristic lengths & some transfinite number of divisions
@@ -243,3 +243,14 @@ If(CoefGeo==1 && !Flag_PML_Cyl)
Physical Surface(SKINDIPOLEUP) = skindipole1[];
Physical Surface(SKINFEED) = skindipole2[];
EndIf
+
+
+Field[1] = Ball;
+Field[1].VIn = lcd/5;
+Field[1].VOut = lc;
+Field[1].XCenter = 0;
+Field[1].YCenter = 0;
+Field[1].ZCenter = 0;
+Field[1].Radius = 1*delta_gap;
+Field[1].Thickness = 4*delta_gap;
+Background Field = 1;
diff --git a/Antennas/dipole_data.geo b/Antennas/dipole_data.geo
index 78c6a30..3b8a1ef 100644
--- a/Antennas/dipole_data.geo
+++ b/Antennas/dipole_data.geo
@@ -71,7 +71,7 @@ DefineConstant[
rdipole = rdipole*mm; // in [m]
AngleWedge = ((Flag_PML_Cyl==0) ? AngleWedge_deg : AngleWedgeCyl_deg) * deg2rad ;
-CoefGeo = (!Flag_3Dmodel) ? 2*Pi : 2*Pi/AngleWedge; // axisymmetry in 2D, 1/8 or 1/4 of the 3D model
+CoefGeo = (!Flag_3Dmodel) ? 2*Pi : (2*Pi/AngleWedge); // axisymmetry in 2D, 1/8 or 1/4 of the 3D model
Printf("CoefGeo %g", CoefGeo);
--
GitLab