PyStratify against NGSolve

← Verification · Cylinders · Point dipole

Under active development. Verification is growing as comparisons are added and rerun: the matrix, its figures and their presentation are not yet in their final form. Spotted an error, or missing a topic? Tell us.

NGSolve computes with finite elements. Each case computes the same quantity with PyStratify and with NGSolve and plots both on the same axes.

Step-index fibre

total decay rates against NGSolve's 3D finite elements

Open in calculator
0.70.80.911.11.20.60.811.21.41.61.822.2λ [µm]Γ / Γ₀
radial, PyStratify azimuthal, PyStratify axial, PyStratify radial, NGSolve azimuthal, NGSolve axial, NGSolve
radialazimuthalaxialPyStratifyNGSolve, 10 points

Source

Source typeElectric point dipole
Offset from origin · r [nm]250
λ min [nm]620
λ max [nm]1200
points60
Cross sectionpR 150 nmr 250 nm⊥ radial⊥∥φ azimuthal∥φ∥z axial∥z
SidepR 150 nmr 250 nm⊥ radial⊥∥φ azimuthal∥φ∥z axial∥z
Cross sectionpR 150 nmr 250 nm⊥ radial⊥∥φ azimuthal∥φ∥z axial∥z
SidepR 150 nmr 250 nm⊥ radial⊥∥φ azimuthal∥φ∥z axial∥z

Cylinder · Electric dipole

LabelMaterialnkThickness [nm]
CFibreconst20150
HVacuumconst10∞