Electromagnetic Waves in Lossy Medium

Uniform Plane Wave in Lossy Medium

Power Dissipation

  • A uniform plane wave propagates through a lossy medium.
  • Objective: Find the power dissipated in a specific volume.
  • Given parameters:
    • Power density: σ=2S/m\,\sigma = 2 \,\text{S/m}
    • Electric field amplitude: Eo=275V/mE_o = 275 \,\text{V/m}
    • Frequency: f=75GHzf = 75 \,\text{GHz}
    • Relative permittivity: εr=72\varepsilon_r = 72
    • Length: l=3δl = 3\delta
    • Area: A=1m2A = 1 \,\text{m}^2
    • Where δ\delta is skin depth

Lossy Medium Characteristics

  • Distinction between lossy and low-loss mediums.
  • Power density calculation:
    • σωε=22π×75×109×72×8.854×1012=0.00667\frac{\sigma}{\omega \varepsilon} = \frac{2}{2 \pi \times 75 \times 10^9 \times 72 \times 8.854 \times 10^{-12}} = 0.00667
    • εε=0.00667\frac{\varepsilon''}{\varepsilon'} = 0.00667
  • With f=75GHzf = 75 \,\text{GHz}, E<em>o=275V/mE<em>o = 275 \,\text{V/m}, σ=2S/m\sigma = 2 \,\text{S/m}, ε</em>r=72\varepsilon</em>r = 72, l=3δl = 3\delta, and A=1m2A = 1 \,\text{m}^2

Further Calculations

  • Attenuation constant:
    • α=ωμε21+(σωε)21σ2με\alpha = \omega \sqrt{\frac{\mu \varepsilon}{2}} \sqrt{\sqrt{1 + (\frac{\sigma}{\omega \varepsilon})^2} - 1} \approx \frac{\sigma}{2} \sqrt{\frac{\mu}{\varepsilon}}
  • Intrinsic impedance:
    • η<em>c=με=μ</em>0ε<em>0ε</em>r=120πεr=120π7244.43Ω\eta<em>c = \sqrt{\frac{\mu}{\varepsilon}} = \sqrt{\frac{\mu</em>0}{\varepsilon<em>0 \varepsilon</em>r}} = \frac{120\pi}{\sqrt{\varepsilon_r}} = \frac{120\pi}{\sqrt{72}} \approx 44.43 \Omega
  • With f=75GHzf = 75 \,\text{GHz}, E<em>o=275V/mE<em>o = 275 \,\text{V/m}, σ=2S/m\sigma = 2 \,\text{S/m}, ε</em>r=72\varepsilon</em>r = 72, l=3δl = 3\delta, and A=1m2A = 1 \,\text{m}^2

Power Dissipated

  • Power dissipated is the difference between input and output power: P<em>diss=P</em>inPoutP<em>{diss} = P</em>{in} - P_{out}
  • Where P<em>inP<em>{in} is the input power and P</em>outP</em>{out} is the output power
  • The Poynting vector: Sav=12Re[E×H]\vec{S}_{av} = \frac{1}{2} \text{Re}[\vec{E} \times \vec{H}^*]
  • P=ASP = A S
  • With l=3δl = 3\delta and A=1m2A = 1 \,\text{m}^2

Power Density Calculation

  • S<em>av=E</em>o22η=27522×44.43=851W/m2S<em>{av} = \frac{|E</em>o|^2}{2 \eta} = \frac{275^2}{2 \times 44.43} = 851 \,\text{W/m}^2
  • P<em>diss=A(S</em>av,inSav,out)=8512.11848.9WP<em>{diss} = A (S</em>{av,in} - S_{av,out}) = 851 - 2.11 \approx 848.9 \,\text{W}
  • S<em>out=S</em>ine2αz=851e2αzS<em>{out} = S</em>{in} e^{-2\alpha z} = 851 e^{-2\alpha z}, where α=0.00248\alpha = 0.00248
  • The power decays exponentially with distance z.
  • With l=3δl = 3\delta and A=1m2A = 1 \,\text{m}^2

2018 Exam Q4

Problem Statement

  • A linearly polarized uniform plane wave in a lossy medium (μμ<em>0,εε</em>0\mu \neq \mu<em>0, \varepsilon \neq \varepsilon</em>0) is propagating in the +z direction.
  • At f=1MHzf = 1 \,\text{MHz}, the intrinsic impedance is ηc=7845Ω\eta_c = 78 \angle 45^\circ \,\Omega and the skin depth is δ=0.25m\delta = 0.25 \,\text{m}.
  • The electric field magnitude is E<em>1=100V/mE<em>1 = 100 \,\text{V/m}, and its phase at z=0z = 0 and t=0t = 0 is ϕ</em>1=60\phi</em>1 = 60^\circ.
  • The field vector at z=0z = 0 and t=0t = 0 points in the +y direction.

Objectives

  1. Determine the wavelength in the medium (λ\lambda).
  2. Determine the conductivity of the medium (σ\sigma).
  3. Determine the intrinsic impedance of the medium (ηc\eta_c).
  4. Obtain expressions for the instantaneous electric field E(z,t)\vec{E}(z, t).
  5. Obtain expressions for the instantaneous magnetic field H(z,t)\vec{H}(z, t).
  6. Find the time-averaged power density vector Sav\vec{S}_{av}.

Solution

Good Conductor Approximation
  • Since it's a good conductor, α=β=1δ\alpha = \beta = \frac{1}{\delta}.
  • And ηc=(1+j)ασ=jωμσ\eta_c = (1 + j) \frac{\alpha}{\sigma} = \sqrt{\frac{j \omega \mu}{\sigma}}.
Wavelength Calculation
  • λ=2πβ=2πδ=2π(0.25)=1.5708m\lambda = \frac{2\pi}{\beta} = 2 \pi \delta = 2 \pi (0.25) = 1.5708 \,\text{m}
Conductivity Calculation
  • δ=2ωμσσ=2ωμδ2=22πfμδ2=1πfμδ2\delta = \sqrt{\frac{2}{\omega \mu \sigma}} \Rightarrow \sigma = \frac{2}{\omega \mu \delta^2} = \frac{2}{2 \pi f \mu \delta^2} = \frac{1}{\pi f \mu \delta^2}
  • σ=1π×1×106×4π×107×0.252=4.05S/m\sigma = \frac{1}{\pi \times 1 \times 10^6 \times 4\pi \times 10^{-7} \times 0.25^2} = 4.05 \,\text{S/m}
  • Alternative:
    • α=πfμσσ=α2πfμ=4πfμ\alpha = \sqrt{\pi f \mu \sigma} \Rightarrow \sigma = \frac{\alpha^2}{\pi f \mu} = \frac{4}{\pi f \mu}
Intrinsic Impedance Calculation
  • ηc=(1+j)ασ=(1+j)πfμσ\eta_c = (1 + j) \frac{\alpha}{\sigma} = (1 + j) \sqrt{\frac{\pi f \mu}{\sigma}}
Electric Field Expression
  • E(z,t)=Re[y^E<em>0eαzej(ωtβz+ϕ</em>0)]\vec{E}(z, t) = \text{Re} [\hat{y} E<em>0 e^{-\alpha z} e^{j(\omega t - \beta z + \phi</em>0)}]
  • E(z,t)=y^100e4zcos(2π×106t4z+60)V/m\vec{E}(z, t) = \hat{y} 100 e^{-4z} \cos(2\pi \times 10^6 t - 4z + 60^\circ) \,\text{V/m}
Magnetic Field Expression
  • H=1ηa<em>z^×E=x^E</em>0ηeαzej(ωtβz+ϕ045)\vec{H} = -\frac{1}{\eta} \hat{a<em>z} \times \vec{E} = -\hat{x} \frac{E</em>0}{\eta} e^{-\alpha z} e^{j(\omega t - \beta z + \phi_0 - 45^\circ)}
  • H(z,t)=x^71.6e4zcos(2π×106t4z+15)A/m\vec{H}(z, t) = -\hat{x} 71.6 e^{-4z} \cos(2\pi \times 10^6 t - 4z + 15^\circ) \,\text{A/m}
Time-Averaged Power Density Vector
  • S<em>av=12Re[E×H]=z^12E</em>02ηe2αzcos(θη)\vec{S}<em>{av} = \frac{1}{2} \text{Re} [\vec{E} \times \vec{H}^*] = \hat{z} \frac{1}{2} \frac{|E</em>0|^2}{|\eta|} e^{-2\alpha z} \cos(\theta_{\eta})
  • Sav=z^12100278e8zcos(45)=z^45.31e8zW/m2\vec{S}_{av} = \hat{z} \frac{1}{2} \frac{100^2}{78} e^{-8z} \cos(45^\circ) = \hat{z} 45.31 e^{-8z} \,\text{W/m}^2

2013 Test 1 Q2

Problem Statement

  • A uniform plane, time-harmonic electromagnetic wave radiated by a short dipole antenna is given by:
    • E(R)=(θ^+φ^)E0ReαRejβRV/mE(R) = (\hat{\theta} + \hat{\varphi}) \frac{E_0}{R} e^{-\alpha R} e^{-j\beta R} \,\text{V/m}
  • The antenna is used in a communication system with a submarine in seawater.
    • εr=81\varepsilon_r = 81
    • μr=1\mu_r = 1
    • ε=8×102\varepsilon'' = 8 \times 10^{-2}
  • The system must operate to a depth of at least 1.2 km.
  • Acceptable signal-to-noise ratio to a depth of 3 skin depths.

Solution Approach

  1. Determine the skin depth in seawater.
  2. Determine the attenuation constant.
  3. Calculate the required frequency for the communication system.

Calculations

  • Given: ε=81\varepsilon' = 81, ε=8×102\varepsilon'' = 8 \times 10^{-2}, μr=1\mu_r = 1
    • \frac{\varepsilon''}{\varepsilon'} = \frac{8 \times 10^{-2}}{81} = 0.000987 < < 1
Skin Depth
  • Operating depth: d=1.2km=1200md = 1.2 \,\text{km} = 1200 \,\text{m}

  • d=3δδ=d3=12003=400mδ=1/αd = 3\delta \quad \Rightarrow \quad \delta = \frac{d}{3} = \frac{1200}{3} = 400 \,\text{m} \quad \delta = 1/\alpha

  • Frequency: f=41.6MHzf = 41.6 \,\text{MHz}

  • Intrinsic impedance:

    • η<em>c=με=μ</em>0μ<em>rε</em>0ε<em>r=120πε</em>r=120π81=40π341.89Ω\eta<em>c = \sqrt{\frac{\mu}{\varepsilon}} = \sqrt{\frac{\mu</em>0 \mu<em>r}{\varepsilon</em>0 \varepsilon<em>r}} = \frac{120\pi}{\sqrt{\varepsilon</em>r}} = \frac{120\pi}{\sqrt{81}} = \frac{40\pi}{3} \approx 41.89 \Omega

Magnetic Field Expression

  • H(R)=θ^E<em>0ηeαRR+φ^E</em>0ηeαRR=E(R)η\vec{H}(R) = -\hat{\theta} \frac{E<em>0}{\eta} \frac{e^{-\alpha R}}{R} + \hat{\varphi} \frac{E</em>0}{\eta} \frac{e^{-\alpha R}}{R} = \frac{E(R)}{\eta}

Time-Averaged Power Density Vector

  • S<em>av=12Re[E×H]=R^E</em>022ηe2αRR2\vec{S}<em>{av} = \frac{1}{2} \text{Re} [\vec{E} \times \vec{H}^*] = \hat{R} \frac{|E</em>0|^2}{2 \eta} \frac{e^{-2 \alpha R}}{R^2}
  • Attenuation in dB:
    • Attenuation=10log<em>10S(R)S(0)=10log</em>10e2αR=10log10e2×3=26.01dB\text{Attenuation} = 10 \log<em>{10} \frac{S(R)}{S(0)} = 10 \log</em>{10} e^{-2\alpha R} = 10 \log_{10} e^{-2 \times 3} = -26.01 \,\text{dB}

2019 Test 2 Q4

Problem Statement

  • Given the following fields:
    • E=[z^]E0ejkx\vec{E} = [\hat{z}] E_0 e^{-jkx}
    • H=[y^]H0ejkx\vec{H} = [\hat{y}] H_0 e^{-jkx}

Calculations

  • Need to apply these general formulas:
    • η=με\eta = \sqrt{\frac{\mu}{\varepsilon}}
    • k=ωμεk = \omega \sqrt{\mu \varepsilon}

Numerical Values

  • ωπ=108ω=2π×108rad/s\frac{\omega}{\pi} = 10^8 \Rightarrow \omega = 2 \pi \times 10^8 \,\text{rad/s}
  • f=ω2π=25MHzf = \frac{\omega}{2 \pi} = 25 \,\text{MHz}
  • k=ωμε=2πfμ<em>rμ</em>0ε<em>rε</em>0k = \omega \sqrt{\mu \varepsilon} = 2 \pi f \sqrt{\mu<em>r \mu</em>0 \varepsilon<em>r \varepsilon</em>0}

Rewrite magnetic Field

  • Examine the rewritten form:

    • H=[x^]H<em>0cos(2π108t)\vec{H} = [\hat{x}] H<em>0 \cos(2 \pi 10^{8} t) [y^]H</em>0cos(2π108t)[\hat{y}] H</em>0 \cos(2 \pi 10^{8} t)
    Conclusion
  • Wave polarization: Left-hand circular polarization (LHC).