From 9579cfafbc21f14afec106696ede9ad21be4e863 Mon Sep 17 00:00:00 2001
From: Guillaume Demesy <guillaume.demesy@fresnel.fr>
Date: Thu, 21 Nov 2019 12:56:58 +0100
Subject: [PATCH] test 2D
---
DiffractionGratings/grating3D.pro | 7 ++-----
.../grating3D_data_2Dlamellar.geo | 16 ++++++++--------
2 files changed, 10 insertions(+), 13 deletions(-)
diff --git a/DiffractionGratings/grating3D.pro b/DiffractionGratings/grating3D.pro
index 19699c0..7525723 100644
--- a/DiffractionGratings/grating3D.pro
+++ b/DiffractionGratings/grating3D.pro
@@ -10,7 +10,6 @@ Include "grating3D_materials.pro"
myDir = "run_results3D/";
Group {
-
// SubDomains
PMLbot = Region[1];
L_6_temp = Region[2];
@@ -322,6 +321,7 @@ PostProcessing {
Quantity {
{ Name u ; Value { Local { [ {u} ]; In Omega; Jacobian JVol; } } }
{ Name Etot ; Value { Local { [ {u}+E1[] ]; In Omega; Jacobian JVol; } } }
+ { Name Htot ; Value { Local { [ H1[]-I[]/(mur[]*mu0*om0)*{Curl u}]; In Omega; Jacobian JVol; } } }
{ Name Edif ; Value { Local { [ {u}+E1d[] ]; In Omega; Jacobian JVol; } } }
{ Name E1 ; Value { Local { [ E1[] ]; In Omega; Jacobian JVol; } } }
{ Name H1y ; Value { Local { [CompY[H1[]] ]; In Omega; Jacobian JVol; } } }
@@ -329,7 +329,7 @@ PostProcessing {
{ Name source ; Value { Local { [ source[] ]; In Omega; Jacobian JVol; } } }
{ Name epsr_xx; Value { Local { [ CompXX[epsr[]] ]; In Omega; Jacobian JVol; } } }
{ Name Damp_pml_top; Value { Local { [Damp_pml_top[] ]; In Omega; Jacobian JVol; } } }
- { Name Poy_tot; Value { Local { [ 0.5*Re[Cross[{u}+E1[] , Conj[ H1[]-I[]/(mur[]*mu0*om0)*{Curl u}]]] ]; In Omega; Jacobian JVol; } } }
+ { Name Poy_tot; Value { Local { [ 0.5*Re[Cross[{u}+E1[] , Conj[ H1[]-I[]/(mur[]*mu0*om0)*{Curl u}]]] ]; In Omega; Jacobian JVol; } } }
For k In {2:6}
{ Name Abs_L~{k} ; Value { Integral { [ ep0*om0 * 0.5*Fabs[eps_im_L~{k}]*(SquNorm[{u}+E1[]]) / (Pinc*period_x*period_y) ] ; In L~{k} ; Integration I1 ; Jacobian JVol ; } } }
@@ -364,9 +364,6 @@ PostProcessing {
}
}
-
-
-
PostOperation {
{ Name postop_helmholtz_vector; NameOfPostProcessing postpro_helmholtz_vector ;
Operation {
diff --git a/DiffractionGratings/grating3D_data_2Dlamellar.geo b/DiffractionGratings/grating3D_data_2Dlamellar.geo
index bf7a3b8..0366ab6 100644
--- a/DiffractionGratings/grating3D_data_2Dlamellar.geo
+++ b/DiffractionGratings/grating3D_data_2Dlamellar.geo
@@ -46,10 +46,10 @@ DefineConstant[
eps_re_L_6 = {-44.9757 , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
eps_im_L_6 = {2.9524 , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
- paramaille = {15 , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
+ paramaille = {10 , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
lc_scat = {30 , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
PML_top = {lambda0 , Name StrCat[pp5,"/4PML top thickness [nm]"]},
- PML_bot = {lambda0 , Name StrCat[pp5,"/5PML bot thickness [nm]"]},
+ PML_bot = {lambda0/5 , Name StrCat[pp5,"/5PML bot thickness [nm]"]},
Nmax = {1 , Name StrCat[pp5,"/6Number of non specular order to output [-]"]},
refine_mesh_L_1= {1 , Name StrCat[pp5,"/7refine layers/1refine mesh layer 1 [-]"]},
refine_mesh_L_2= {1 , Name StrCat[pp5,"/7refine layers/2refine mesh layer 2 [-]"]},
@@ -61,12 +61,12 @@ DefineConstant[
InterpSampling = { 30 , Name StrCat[pp6,"/0Interpolation grid step [nm]"]},
Flag_interp_cubic = { 0 , Name StrCat[pp6,"/1Interpolate on cubic grid?"], Choices {0,1} },
FlagOutEtotCuts = { 0 , Name StrCat[pp6,"/2Output Total Electric Field cuts?"] , Choices {0,1} },
- FlagOutHtotCuts = { 1 , Name StrCat[pp6,"/2Output Total Magnetic Field cuts?"] , Choices {0,1} },
- FlagOutEscaCuts = { 0 , Name StrCat[pp6,"/3Output Scattered Electric Field cuts?"] , Choices {0,1} },
- FlagOutPoyCut = { 0 , Name StrCat[pp6,"/4Output Poynting cuts?"] , Choices {0,1} },
- FlagOutEtotFull = { 0 , Name StrCat[pp6,"/5Total Electric Field Full Output?"] , Choices {0,1} },
- FlagOutEscaFull = { 0 , Name StrCat[pp6,"/6Scattered Electric Field Full Output?"] , Choices {0,1} },
- FlagOutPoyFull = { 0 , Name StrCat[pp6,"/7Poynting Full Output?"] , Choices {0,1} }
+ FlagOutHtotCuts = { 1 , Name StrCat[pp6,"/3Output Total Magnetic Field cuts?"] , Choices {0,1} },
+ FlagOutEscaCuts = { 0 , Name StrCat[pp6,"/4Output Scattered Electric Field cuts?"] , Choices {0,1} },
+ FlagOutPoyCut = { 0 , Name StrCat[pp6,"/5Output Poynting cuts?"] , Choices {0,1} },
+ FlagOutEtotFull = { 1 , Name StrCat[pp6,"/6Total Electric Field Full Output?"] , Choices {0,1} },
+ FlagOutEscaFull = { 0 , Name StrCat[pp6,"/7Scattered Electric Field Full Output?"] , Choices {0,1} },
+ FlagOutPoyFull = { 0 , Name StrCat[pp6,"/8Poynting Full Output?"] , Choices {0,1} }
];
lambda_m = lambda0;
--
GitLab