fix unintuitive behavior for scatterd + total fields for diel. sphere

Author: HoBeZwe

src/planeWave/scattered.jl

@@ -46,10 +46,9 @@ function scatteredfield(sphere::Sphere, excitation::PlaneWave, point, quantity::
 
     eps = parameter.relativeAccuracy
 
-    ST = SVector{3,Complex{T}}
-    F = ST(0.0, 0.0, 0.0)
+    F = SVector{3,Complex{T}}(0.0, 0.0, 0.0)
 
-    A₀ = amplitude(sphere, excitation::PlaneWave, quantity, r)
+    A₀ = amplitude(sphere, excitation, quantity, r)
 
     A₀ == 0.0 && return F # Inside of PEC return zero field
 
@@ -84,7 +83,30 @@ function scatteredfield(sphere::Sphere, excitation::PlaneWave, point, quantity::
 
     end
 
-    return convertSpherical2Cartesian(A₀ .* F, point_sph)
+    Fin = inside(sphere, excitation, point, quantity; parameter=parameter)
+
+    return convertSpherical2Cartesian(A₀ .* F, point_sph) + Fin
+end
+
+
+function inside(sphere::Sphere, excitation::PlaneWave{T,R,C}, point, quantity::Field; parameter) where {T,R,C}
+
+    return SVector{3,Complex{R}}(0.0, 0.0, 0.0) # no correction needed
+end
+
+function inside(sphere::DielectricSphere, excitation::PlaneWave{T,R,C}, point, quantity::FarField; parameter) where {T,R,C}
+
+    return SVector{3,Complex{R}}(0.0, 0.0, 0.0) # no correction needed
+end
+
+function inside(sphere::DielectricSphere, excitation::PlaneWave{T,R,C}, point, quantity::Field; parameter) where {T,R,C}
+
+    # inside the sphere the incident field has to be substracted to get only the scattered part
+    if norm(point) < sphere.radius
+        return -field(excitation, point, quantity; parameter=parameter)
+    end
+
+    return SVector{3,Complex{R}}(0.0, 0.0, 0.0)
 end
 
 

src/totalFields.jl

@@ -1,7 +1,6 @@
 
 inside(sphere::Sphere) = 0.0
 inside(sphere::PECSphere) = sphere.radius
-# inside(sphere::DielectricSphere) = sphere.radius
 
 
 """

test/planeWave_dielectric.jl

@@ -77,7 +77,7 @@
 
     # E-Field
     EF₂ = scatteredfield(sp, ex, ElectricField(points_cartNF))
-    EF₁ = scatteredfield(sp, ex, ElectricField(points_cartNF_inside))
+    EF₁ = field(sp, ex, ElectricField(points_cartNF_inside))
 
     diff_EF₂ = norm.(EF₂ - EF₂MoM) ./ maximum(norm.(EF₂))  # worst case error
     diff_EF₁ = norm.(EF₁ - EF₁MoM) ./ maximum(norm.(EF₁))  # worst case error
@@ -90,7 +90,7 @@
     HF₁MoM = hfield(𝓣k1, -(1 / η1)^2 .* m, RT, 𝓚k1, -j, RT, points_cartNF_inside)
 
     HF₂ = scatteredfield(sp, ex, MagneticField(points_cartNF))
-    HF₁ = scatteredfield(sp, ex, MagneticField(points_cartNF_inside))
+    HF₁ = field(sp, ex, MagneticField(points_cartNF_inside))
 
     diff_HF₂ = norm.(HF₂ - HF₂MoM) ./ maximum(norm.(HF₂))  # worst case error
     diff_HF₁ = norm.(HF₁ - HF₁MoM) ./ maximum(norm.(HF₁))  # worst case error