paramrf 0.37.2__tar.gz → 0.37.3__tar.gz

{paramrf-0.37.2/paramrf.egg-info → paramrf-0.37.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: paramrf
-Version: 0.37.2
+Version: 0.37.3
 Summary: Parametric radio frequency modeling
 Author-email: Gary Allen <gvcallen@gmail.com>
 Project-URL: homepage, https://github.com/gvcallen/paramrf

{paramrf-0.37.2 → paramrf-0.37.3/paramrf.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: paramrf
-Version: 0.37.2
+Version: 0.37.3
 Summary: Parametric radio frequency modeling
 Author-email: Gary Allen <gvcallen@gmail.com>
 Project-URL: homepage, https://github.com/gvcallen/paramrf

{paramrf-0.37.2 → paramrf-0.37.3}/pmrf/models/components/lines/empirical.py

@@ -20,10 +20,14 @@ class PhysicalLine(AbstractImmittanceLine):
 
     The per-unit-length parameters are
     $$R = 2 z_n \alpha_c$$
-    $$L = \frac{z_n \sqrt{\varepsilon_r}}{c}$$
+    $$L = \frac{z_n \sqrt{\varepsilon_r}}{c} + \frac{R}{\omega}$$
     $$G = \frac{2 \alpha_d}{z_n}$$
     $$C = \frac{\sqrt{\varepsilon_r}}{z_n c}$$
 
+    The skin-effect term adds conductor reactance to the series impedance:
+    $$Z = R + j\omega z_n\sqrt{\varepsilon_r}/c + jR.$$
+    This contribution vanishes at DC.
+
     Example
     --------
     .. code-block:: python
@@ -93,7 +97,11 @@ class PhysicalLine(AbstractImmittanceLine):
         G = G_val * ones
         C = C_val * ones
         
-        return ImmittanceResult.from_rlgc(R, L, G, C, freq.w)
+        return ImmittanceResult(
+            Z=R + 1j * (freq.w * L + R),
+            Y=G + 1j * freq.w * C,
+            omega=freq.w,
+        )
     
 
 class DatasheetLine(AbstractImmittanceLine):
@@ -108,10 +116,14 @@ class DatasheetLine(AbstractImmittanceLine):
 
     and the per-unit-length parameters are
     $$R = 2 z_n \alpha_c$$
-    $$L = \frac{z_n}{v_f c}$$
+    $$L = \frac{z_n}{v_f c} + \frac{R}{\omega}$$
     $$G = \frac{2 \alpha_d}{z_n}$$
     $$C = \frac{1}{z_n v_f c}$$
 
+    The skin-effect term adds conductor reactance to the series impedance:
+    $$Z = R + j\omega z_n/(v_f c) + jR.$$
+    This contribution vanishes at DC.
+
     Example
     --------
     .. code-block:: python
@@ -186,5 +198,9 @@ class DatasheetLine(AbstractImmittanceLine):
         L = (zn / (vf * c)) * jnp.ones_like(w)
         C = (1.0 / (zn * vf * c)) * jnp.ones_like(w)
         
-        return ImmittanceResult.from_rlgc(R, L, G, C, w)
+        return ImmittanceResult(
+            Z=R + 1j * (w * L + R),
+            Y=G + 1j * w * C,
+            omega=w,
+        )
 

{paramrf-0.37.2 → paramrf-0.37.3}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "paramrf"
-version = "0.37.2"
+version = "0.37.3"
 description = "Parametric radio frequency modeling"
 readme = "README.rst"
 authors = [

{paramrf-0.37.2 → paramrf-0.37.3}/tests/test_models/test_lines.py

@@ -25,8 +25,8 @@ from pmrf.materials import (
     DjordjevicSarkarDielectric,
     RoughConductor,
 )
-from pmrf.materials.surface_impedance import TescheTubeSurfaceImpedance
-from pmrf.models.components.lines.coaxial import TescheCoaxialFormulation
+from pmrf.materials.surface_impedance import HalfSpaceSurfaceImpedance, TescheTubeSurfaceImpedance
+from pmrf.models.components.lines.coaxial import SchelkunoffCoaxialFormulation, TescheCoaxialFormulation
 from pmrf.models.components.lines.microstrip import KirschningJansenMicrostripDispersion
 
 @pytest.fixture
@@ -118,6 +118,38 @@ def test_datasheet_line_attenuation(basic_freq):
     
     assert jnp.allclose(jnp.abs(s[:, 1, 0]), expected_mag, atol=1e-4)
 
+
+@pytest.mark.parametrize("line_type", [DatasheetLine, PhysicalLine])
+@pytest.mark.parametrize("length", [1.0, 10.0])
+def test_empirical_line_matches_half_space_coax(line_type, length):
+    """Skin-effect loss has the same resistance and reactance as a coax half-space."""
+    d_in, d_out = 0.92e-3, 3.0e-3
+    ep_r, tand, sigma = 2.1, 1e-4, 5.8e7
+    a, b = d_in / 2, d_out / 2
+    zn = np.sqrt(mu_0 / epsilon_0) * np.log(b / a) / (2 * np.pi * np.sqrt(ep_r))
+    k1 = (2000 / np.log(10)) * np.sqrt(np.pi * 1e6 * mu_0 / sigma) * (1 / a + 1 / b) / (4 * np.pi * zn)
+    k2 = (2000 / np.log(10)) * np.pi * 1e6 * np.sqrt(ep_r) * tand / c
+
+    with jax.enable_x64():
+        freq = Frequency(start=10, stop=500, npoints=491, unit='MHz')
+        if line_type is DatasheetLine:
+            line = DatasheetLine(length=length, zn=zn, vf=1 / np.sqrt(ep_r), k1=k1, k2=k2)
+        else:
+            line = PhysicalLine(length=length, zn=zn, ep_r=ep_r, A=k1 / 1e5, f_A=1.0, tand=tand)
+        reference = CoaxialLine(
+            d_in=d_in, d_out=d_out, length=length,
+            dielectric=ConstantDielectric(ep_r=ep_r, tand=tand),
+            conductor=BulkConductor(sigma=sigma),
+            formulation=SchelkunoffCoaxialFormulation(
+                inner_impedance=HalfSpaceSurfaceImpedance(),
+                shield_impedance=HalfSpaceSurfaceImpedance(),
+            ),
+        )
+
+        difference = np.abs(np.asarray(line.s(freq) - reference.s(freq)))
+
+    assert np.max(difference) < 1e-8
+
 def test_coaxial_line_impedance(basic_freq):
     """Verify CoaxialLine matches the analytical Zc formula for a lossless coax."""
     # RG-58 dimensions roughly: ep_r=2.25, d_in=0.9mm, d_out=2.95mm