build based on b5fb626

dev/.documenter-siteinfo.json

@@ -1 +1 @@
-{"documenter":{"julia_version":"1.10.8","generation_timestamp":"2025-02-05T11:19:52","documenter_version":"1.8.0"}}
+{"documenter":{"julia_version":"1.10.8","generation_timestamp":"2025-02-05T12:01:30","documenter_version":"1.8.0"}}

dev/index.html

@@ -24,15 +24,10 @@
  [LP 0,1,1.4345691826486056,1.6,Nothing]
  [LP 0,2,1.4137165538840926,1.6,Nothing]</code></pre><p>Computation of the guided modes <span>$\ell=0$</span> with the cladding:</p><pre><code class="language-julia hljs">julia&gt; m2=multi_step_fiber_modes(1.6,0,[4,20],[1.44,1.41,1.43])
 1-element Vector{Mode}:
- [LP 0,1,1.4345691826495002,1.6,Nothing]
- ```
-The LP₀₂ becomes leaky when adding the cladding
-
-To compute its losses, we have to use the FEM solver with a PML:</code></pre><p>julia julia&gt; epsilon<em>trench=x-&gt;(1.44-0.03<em>(x[1]&gt;=4)+0.02</em>(x[1]&gt;20))^2; julia&gt; model</em>trench = CartesianDiscreteModel((0,50),5000); julia&gt; ms1D<em>trench=FEM1D(1.6,0,epsilon</em>trench,model<em>trench,field=true,dPML=5,approx</em>neff=m[2].neff) 1-element Vector{Mode{ScalarFieldFEM1D}}:  [Mode LP n°1,1.4137165449018512 + 1.0825770569063048e-8im,1.6,ScalarFieldFEM1D] julia&gt; losses(ms1D_trench[1])*1E6 369.26063441366387</p><pre><code class="nohighlight hljs">Losses of the LP₀₂ mode are 369 dB/km.
-
-### Finding bend losses using the FEM mode solver with a PML
-Let us consider the step_index fiber described below:
-- Core: radius 3.5 µm, refractive index 1.4457
-- Cladding: refractive index 1.4378
-
-Computation of the fundamental mode at $\lambda=1.55$ µm:</code></pre><p>julia julia&gt; ms<em>straight=multi</em>step<em>fiber</em>modes(1.55,0,3.5,[1.4457,1.4378],field=true) 1-element Vector{Mode}:  [LP 0,1,1.4414392004035022,1.55,Nothing]  ```</p><p>To compute the funcamental mode when the fiber is bent, we need to use a 2D mesh generated with GMSH:  <code>julia julia&gt; using GridapGmsh julia&gt; model = GmshDiscreteModel(&quot;./models/example5.msh&quot;) Info    : Reading &#39;./models/example5.msh&#39;... Info    : 10 entities Info    : 14768 nodes Info    : 29536 elements Info    : Done reading &#39;./models/example5.msh&#39; UnstructuredDiscreteModel()</code> and define the relative permittivity for a straight fiber and for fiber bent with a radius of curvature of 10 mm:  <code>julia julia&gt; epsilon2D=x-&gt;(1.4457-(1.4457-1.4378)*(hypot(x[1],x[2])&gt;=3.5))^2; julia&gt; epsilon2D_bent=x-&gt;epsilon2D(x)*(1+x[1]/10E3)^2; julia&gt; ms2D_bent=FEM2D(1.55,epsilon2D_bent,model,field=true,neigs=1,dPML=3,approx_neff=ms_straight[1].neff) 1-element Vector{Mode{ScalarFieldFEM2D}}:  [Mode LP n°1,1.4414582242320717 + 1.5096368959704397e-9im,1.55,ScalarFieldFEM2D] losses(ms2D_bent[1])*1E6</code>  In this case, the bend losses are 53 dB/km</p><h2 id="Credits"><a class="docs-heading-anchor" href="#Credits">Credits</a><a id="Credits-1"></a><a class="docs-heading-anchor-permalink" href="#Credits" title="Permalink"></a></h2><p>OpticalFibers.jl is maintained by Olivier Vanvincq (<a href="https://www.univ-lille.fr/">University of Lille</a>, <a href="https://phlam.univ-lille.fr/">PhLAM laboratory</a>).</p></article><nav class="docs-footer"><a class="docs-footer-nextpage" href="Common/">Common functions »</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="auto">Automatic (OS)</option><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="catppuccin-latte">catppuccin-latte</option><option value="catppuccin-frappe">catppuccin-frappe</option><option value="catppuccin-macchiato">catppuccin-macchiato</option><option value="catppuccin-mocha">catppuccin-mocha</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.8.0 on <span class="colophon-date" title="Wednesday 5 February 2025 11:19">Wednesday 5 February 2025</span>. Using Julia version 1.10.8.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>
+ [LP 0,1,1.4345691826495002,1.6,Nothing]</code></pre><p>The LP₀₂ becomes leaky when adding the cladding</p><p>To compute its losses, we have to use the FEM solver with a PML:</p><pre><code class="language-julia hljs">julia&gt; epsilon_trench=x-&gt;(1.44-0.03*(x[1]&gt;=4)+0.02*(x[1]&gt;20))^2;
+julia&gt; model_trench = CartesianDiscreteModel((0,50),5000);
+julia&gt; ms1D_trench=FEM1D(1.6,0,epsilon_trench,model_trench,field=true,dPML=5,approx_neff=m[2].neff)
+1-element Vector{Mode{ScalarFieldFEM1D}}:
+ [Mode LP n°1,1.4137165449018512 + 1.0825770569063048e-8im,1.6,ScalarFieldFEM1D]
+julia&gt; losses(ms1D_trench[1])*1E6
+369.26063441366387</code></pre><p>Losses of the LP₀₂ mode are 369 dB/km.</p><h3 id="Finding-bend-losses-using-the-FEM-mode-solver-with-a-PML"><a class="docs-heading-anchor" href="#Finding-bend-losses-using-the-FEM-mode-solver-with-a-PML">Finding bend losses using the FEM mode solver with a PML</a><a id="Finding-bend-losses-using-the-FEM-mode-solver-with-a-PML-1"></a><a class="docs-heading-anchor-permalink" href="#Finding-bend-losses-using-the-FEM-mode-solver-with-a-PML" title="Permalink"></a></h3><p>Let us consider the step_index fiber described below:</p><ul><li>Core: radius 3.5 µm, refractive index 1.4457</li><li>Cladding: refractive index 1.4378</li></ul><p>Computation of the fundamental mode at <span>$\lambda=1.55$</span> µm: <code>julia julia&gt; ms_straight=multi_step_fiber_modes(1.55,0,3.5,[1.4457,1.4378],field=true) 1-element Vector{Mode}:  [LP 0,1,1.4414392004035022,1.55,Nothing]</code></p><p>To compute the funcamental mode when the fiber is bent, we need to use a 2D mesh generated with GMSH:  <code>julia julia&gt; using GridapGmsh julia&gt; model = GmshDiscreteModel(&quot;./models/example5.msh&quot;) Info    : Reading &#39;./models/example5.msh&#39;... Info    : 10 entities Info    : 14768 nodes Info    : 29536 elements Info    : Done reading &#39;./models/example5.msh&#39; UnstructuredDiscreteModel()</code> and define the relative permittivity for a straight fiber and for fiber bent with a radius of curvature of 10 mm:  <code>julia julia&gt; epsilon2D=x-&gt;(1.4457-(1.4457-1.4378)*(hypot(x[1],x[2])&gt;=3.5))^2; julia&gt; epsilon2D_bent=x-&gt;epsilon2D(x)*(1+x[1]/10E3)^2; julia&gt; ms2D_bent=FEM2D(1.55,epsilon2D_bent,model,field=true,neigs=1,dPML=3,approx_neff=ms_straight[1].neff) 1-element Vector{Mode{ScalarFieldFEM2D}}:  [Mode LP n°1,1.4414582242320717 + 1.5096368959704397e-9im,1.55,ScalarFieldFEM2D] losses(ms2D_bent[1])*1E6</code>  In this case, the bend losses are 53 dB/km</p><h2 id="Credits"><a class="docs-heading-anchor" href="#Credits">Credits</a><a id="Credits-1"></a><a class="docs-heading-anchor-permalink" href="#Credits" title="Permalink"></a></h2><p>OpticalFibers.jl is maintained by Olivier Vanvincq (<a href="https://www.univ-lille.fr/">University of Lille</a>, <a href="https://phlam.univ-lille.fr/">PhLAM laboratory</a>).</p></article><nav class="docs-footer"><a class="docs-footer-nextpage" href="Common/">Common functions »</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="auto">Automatic (OS)</option><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="catppuccin-latte">catppuccin-latte</option><option value="catppuccin-frappe">catppuccin-frappe</option><option value="catppuccin-macchiato">catppuccin-macchiato</option><option value="catppuccin-mocha">catppuccin-mocha</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.8.0 on <span class="colophon-date" title="Wednesday 5 February 2025 12:01">Wednesday 5 February 2025</span>. Using Julia version 1.10.8.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>