在有损耗介质中如何计算电磁波传播速度？

Question

在有损耗介质中，电磁波的传播速度如何计算？传统公式 $ v = c / \sqrt{\mu_r \varepsilon_r} $ 仅适用于理想无损耗介质。当介质存在电导率 $ \sigma $ 或复介电常数 $ \varepsilon = \varepsilon' - j\varepsilon'' $ 时，电磁波会产生衰减和相位延迟。此时，传播常数 $ \gamma = \alpha + j\beta $ 中的相位常数 $ \beta $ 决定了波速，相速度 $ v_p = \omega / \beta $ 需基于复介电常数或损耗角正切重新推导。常见问题是如何从介质参数 $ \varepsilon', \varepsilon'', \mu, \sigma $ 准确求解 $ \beta $，进而得到实际传播速度。这在高频通信、雷达探测和材料表征中尤为关键。

Answer 1

在有损耗介质中电磁波传播速度的精确计算方法

1. 从理想介质到有损耗介质：传播速度公式的演进

在真空或理想无损耗介质中，电磁波的相速度由经典公式给出：

$$ v = \frac{c}{\sqrt{\mu_r \varepsilon_r}} $$

其中 $ c $ 是光速，$ \mu_r $ 和 $ \varepsilon_r $ 分别为相对磁导率和相对介电常数。然而，当介质具有电导率 $ \sigma $ 或复介电常数 $ \varepsilon = \varepsilon' - j\varepsilon'' $ 时，该公式不再适用。

实际工程中，如高频PCB材料（如Rogers）、生物组织、大气层等均表现为有损耗特性，必须引入复传播常数 $ \gamma = \alpha + j\beta $ 来描述波的衰减与相位变化。

2. 复介电常数与电导率的关系建模

在时谐场（$ e^{j\omega t} $）下，电导率 $ \sigma $ 可等效为介电常数的虚部：

$$ \varepsilon_{\text{eff}} = \varepsilon' - j\left( \varepsilon'' + \frac{\sigma}{\omega} \right) $$

因此，总的复介电常数可统一表示为：

$$ \tilde{\varepsilon} = \varepsilon' - j\varepsilon_{\text{total}}'' \quad \text{其中} \quad \varepsilon_{\text{total}}'' = \varepsilon'' + \frac{\sigma}{\omega} $$

这一模型将电导损耗与极化损耗统一纳入分析框架，适用于宽频带电磁仿真与材料表征。

3. 传播常数的推导与相位常数求解

麦克斯韦方程导出的波动方程在均匀各向同性介质中的解依赖于传播常数：

$$ \gamma = j\omega\sqrt{\mu \tilde{\varepsilon}} = \alpha + j\beta $$

将其展开为实部与虚部形式：

$$ \alpha = \omega \sqrt{ \frac{\mu \varepsilon'}{2} \left( \sqrt{1 + \tan^2 \delta} - 1 \right) }, \quad \beta = \omega \sqrt{ \frac{\mu \varepsilon'}{2} \left( \sqrt{1 + \tan^2 \delta} + 1 \right) } $$

其中 $ \tan \delta = \frac{\varepsilon_{\text{total}}''}{\varepsilon'} $ 为损耗角正切，是衡量介质损耗的关键参数。

4. 相速度的准确表达式

一旦获得相位常数 $ \beta $，即可计算实际相速度：

$$ v_p = \frac{\omega}{\beta} = \frac{1}{\sqrt{ \mu \varepsilon' } \cdot \sqrt{ \frac{ \sqrt{1 + \tan^2 \delta} + 1 }{2} }} $$

当 $ \tan \delta \to 0 $（无损耗），上式退化为传统公式；而随着 $ \tan \delta $ 增大，相速度下降且波长缩短。

下表列出不同介质在10 GHz下的典型参数对比：

介质类型	$\varepsilon'$	$\tan\delta$	$\sigma$ (S/m)	$\beta$ (rad/m)	$v_p/c$
Air	1.0	0.0001	0	209.4	0.99995
FEP	2.1	0.0003	0	302.8	0.690
FR-4	4.4	0.02	0.01	458.2	0.457
Silicon	11.7	0.001	0.001	762.3	0.275
Water (20°C)	80	0.1	0.01	1420	0.148
Human Muscle	52	0.6	0.8	1890	0.111
Rogers RO4350B	3.66	0.0037	0	380.1	0.552
GaAs	12.9	0.0002	0	795.6	0.264
Alumina

9.8

0.0004	0	618.4	0.339
Seawater	70	3.5	5.0	3140	0.063

5. 工程应用中的数值求解流程

在高频通信系统设计中，通常需通过以下步骤精确计算传播速度：

1. 获取材料的 $ \varepsilon', \varepsilon'', \sigma, \mu $ 参数（可通过VNA测量或厂商手册）
2. 计算角频率 $ \omega = 2\pi f $
3. 构建有效复介电常数 $ \tilde{\varepsilon} = \varepsilon' - j(\varepsilon'' + \sigma/\omega) $
4. 计算复传播常数 $ \gamma = j\omega\sqrt{\mu \tilde{\varepsilon}} $
5. 提取 $ \beta = \operatorname{Im}(\gamma) $
6. 计算相速度 $ v_p = \omega / \beta $
7. 可进一步求群速度 $ v_g = d\omega/d\beta $ 用于色散分析

6. Python代码实现示例

import numpy as np

def calculate_phase_velocity(f, eps_prime, eps_double_prime, sigma, mu_r=1):
    omega = 2 * np.pi * f
    mu = 4e-7 * np.pi * mu_r
    eps_0 = 8.854e-12
    eps_total_imag = eps_double_prime + sigma / omega
    epsilon_tilde = eps_0 * (eps_prime - 1j * eps_total_imag)
    
    gamma = 1j * omega * np.sqrt(mu * epsilon_tilde)
    beta = np.imag(gamma)
    
    if beta == 0:
        return float('inf')
    vp = omega / beta
    return vp

# 示例：计算FR-4在10GHz下的相速度
vp = calculate_phase_velocity(
    f=10e9,
    eps_prime=4.4,
    eps_double_prime=0.0,
    sigma=0.01
)
print(f"Phase velocity: {vp:.2e} m/s ({vp/3e8:.3f}c)")

7. Mermaid流程图：电磁波在损耗介质中传播分析流程

graph TD
    A[输入介质参数] --> B[ε', ε'', σ, μ, f]
    B --> C{是否考虑电导率?}
    C -- 是 --> D[计算等效复介电常数]
    C -- 否 --> E[直接使用ε' - jε'']
    D --> F[计算γ = jω√(με̃)]
    E --> F
    F --> G[分离α和β]
    G --> H[计算相速度 v_p = ω/β]
    H --> I[输出传播特性]
    I --> J[应用于电路设计/成像算法/信道建模]