Skip to content

Commit

Permalink
Minor fixes to doc file
Browse files Browse the repository at this point in the history
  • Loading branch information
hughcars committed Sep 4, 2024
1 parent 880c1ff commit 4c2fb95
Show file tree
Hide file tree
Showing 2 changed files with 21 additions and 17 deletions.
Binary file modified docs/src/assets/examples/cavity-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
38 changes: 21 additions & 17 deletions docs/src/examples/cylinder.md
Original file line number Diff line number Diff line change
Expand Up @@ -115,7 +115,7 @@ the above formula for this problem are listed in the table below.
| ``(2,1,2)`` | ``5.290341\text{ GHz}`` | ``7.269033\text{ GHz}`` |

First, we examine the output of the `cavity_pec.json` simulation. The file
`postpro/pec/eig.csv` contains information about the computed eigenfrequencies and
`postpro/cavity_pec/eig.csv` contains information about the computed eigenfrequencies and
associated quality factors:

```
Expand All @@ -142,7 +142,7 @@ obtained by *Palace*. Since the only source of loss in the simulation is the non
dielectric loss tangent, we have ``Q = Q_d = 1/0.0004 = 2.50\times 10^3`` in all cases.

Next, we run `cavity_impedance.json`, which adds the surface impedance boundary condition.
Examining `postpro/impedance/eig.csv` we see that the mode frequencies are roughly
Examining `postpro/cavity_impedance/eig.csv` we see that the mode frequencies are roughly
unchanged but the quality factors have fallen due to the addition of imperfectly conducting
walls to the model:

Expand All @@ -166,7 +166,7 @@ walls to the model:
```

However, the bulk dielectric loss postprocessing results, computed from the energies written
to `postpro/impedance/domain-E.csv`, still give ``Q_d = 1/0.004 = 2.50\times 10^3`` for
to `postpro/cavity_impedance/domain-E.csv`, still give ``Q_d = 1/0.004 = 2.50\times 10^3`` for
every mode as expected.

Focusing on the ``\text{TE}_{011}`` mode with ``f_{\text{TE},010} = 5.00\text{ GHz}``, we
Expand Down Expand Up @@ -221,23 +221,24 @@ may only become significant on sufficiently refined meshes.

This example demonstrates the eigenmode simulation type in *Palace* to solve for the
cutoff-frequencies of a circular waveguide. As with the cavity the interior material to be
Silicon (``\varepsilon_r = 2.08``,
``\tan\delta = 4\times 10^{-4}``), with cylindrical domain radius ``a = 2.74\text{ cm}``,
and length ``d=2*a = 5.48\text{ cm}``.
Periodic boundary conditions (BCs) are applied in the $z$-direction.
According to [[1]](#References), the cutoff frequencies for the transverse electric and
magnetic modes are given by the formulae:
Silicon (``\varepsilon_r = 2.08``, ``\tan\delta = 4\times 10^{-4}``), with cylindrical
domain radius ``a = 2.74\text{ cm}``, and length ``d=2a = 5.48\text{ cm}``, however now
periodic boundary conditions are applied in the $z$-direction. According to
[[1]](#References), the cutoff frequencies for the transverse electric and magnetic modes
are given by the formulae:

```math
f_{\text{TE},nm} = \frac{1}{2\pi\sqrt{\mu\varepsilon}} \frac{p'_{nm}}{a} \qquad
f_{\text{TM},nm} = \frac{1}{2\pi\sqrt{\mu\varepsilon}} \frac{p_{nm}}{a}
\begin{aligned}
f_{\text{TE},nm} &= \frac{1}{2\pi\sqrt{\mu\varepsilon}} \frac{p'_{nm}}{a}\\
f_{\text{TM},nm} &= \frac{1}{2\pi\sqrt{\mu\varepsilon}} \frac{p_{nm}}{a}
\end{aligned}
```

which are identical to those for the cavity modes, in the special case of ``l=0``.

In addition to these pure waveguide modes, there are aliasing cavity
modes corresponding to a full wavelength in the computational domain (``l==2``). In a
practical problem these are suppressed by choosing a smaller value of ``d`` which shifts
modes corresponding to a full wavelength in the computational domain (``l=2``). In a
practical problem these can be suppressed by choosing a smaller value of ``d`` which shifts
such modes to higher frequencies. The relevant modes are tabulated as

| ``(n,m,l)`` | ``f_{\text{TE}}`` | ``f_{\text{TM}}`` |
Expand All @@ -263,6 +264,9 @@ in the `"Boundaries"` object: `waveguide.json` specifies a perfect electric cond
attribute pairs are defined by `"DonorAttributes"` and `"ReceiverAttributes"`, and the
distance between them is given by the `"Translation"` vector in mesh units.

The file `postpro/waveguide/eig.csv` contains information about the computed eigenfrequencies and
associated quality factors:

```
m, Re{f} (GHz), Im{f} (GHz), Q,
1.000000000e+00, +2.223255722e+00, +4.446511256e-04, +2.500000155e+03,
Expand All @@ -285,10 +289,10 @@ distance between them is given by the `"Translation"` vector in mesh units.
```

In common with the PEC cavity ``Q = Q_d = 1/0.0004 = 2.50\times 10^3`` in all cases, and all
the anticipated waveguide modes are recovered with ``TE_{1,1}`` having the lowest cutoff
frequency followed by ``TM_{0,1}`` and ``TE_{2,1}``, while the aliasing mode ``TE_{1,1,2}`` has
marginally lower frequency than the waveguide modes ``TE_{0,1}`` and ``TM_{1,1}``
(``4.397\text{ GHz}`` compared to ``4.627\text{ GHz}``).
the anticipated waveguide modes are recovered with ``\text{TE}_{1,1}`` having the lowest
cutoff frequency followed by ``\text{TM}_{0,1}`` and ``\text{TE}_{2,1}``, while the aliasing
mode ``\text{TE}_{1,1,2}`` has marginally lower frequency than the waveguide modes
``\text{TE}_{0,1}`` and ``\text{TM}_{1,1}`` (``4.397\text{ GHz}`` compared to ``4.627\text{ GHz}``) and is thus found first.

## References

Expand Down

0 comments on commit 4c2fb95

Please sign in to comment.