imported>Karsten11 am 24. April 2023 um 13:48 Uhr

2023-04-24T13:48:40Z

Neue Seite

Die '''Finite-Integral-Methode''' basiert auf der '''Finite Integration Theorie''' (FIT) und ist ein [[Numerik|numerisches]] [[Simulation]]sverfahren zur [[näherungsfrei]]en Lösung der [[Maxwellsche Gleichungen|elektromagnetischen Grundgleichungen]] nach [[James Clerk Maxwell|Maxwell]]. Sie bildet die mathematische Grundlage von Simulationsprogrammen für elektromagnetische Probleme wie z. B. MAFIA und ''CST MICROWAVE STUDIO®''.

== Grundlagen ==
Die Finite-Integral-Methode, zuerst 1976 von [[Thomas Weiland]] vorgestellt, löst die nach Maxwell benannten elektromagnetischen Grundgleichungen

:<math>
\oint_{C(A)} \vec{E} \cdot d \vec{s} \;=\; - \int_A \frac{\partial}{\partial t}\vec{B} \cdot d \vec{A}
</math>

:<math>
\oint_{C(A)} \vec{H} \cdot d \vec{s} \;=\; \int_A \left( \frac{\partial}{\partial t}\vec{D} + \vec{J}\right)\cdot d \vec{A}
</math>

:<math>
\oint_{A(V)} \vec{D} \cdot d \vec{A} \;=\; \int_V \rho \cdot dV
</math>

:<math>
\oint_{A(V)} \vec{B} \cdot d \vec{A} \;=\; 0
</math>

näherungsfrei in deren Integralform und mittels [[Integralapproximation]] die Materialgleichungen

:<math>
\vec{D} \;=\; \varepsilon_0 \varepsilon_r \vec{E}
</math>

:<math>
\vec{B} \;=\; \mu_0 \mu_r \vec{H}
</math>

:<math>
\vec{J} \;=\; \kappa \vec{E} + \vec{J_s}
</math>

in [[Diskretisierung|diskretisierter]] Form.

== Vorgehensweise ==

Das gesamte Problemgebiet wird in ein erstes (oder primäres) dreidimensionales Netz von einzelnen, möglichst kleinen Gitterzellen (engl.: mesh cells) mit den Materialeigenschaften <math>\varepsilon_r, \mu_r, \kappa</math> unterteilt, die jede bzgl. ihrer elektrischen Kantenspannung <math>e_i = E_i \cdot a</math> und ihres magnetischen Flusses durch die Randflächen <math>b_j = B_j \cdot A_j</math> berechnet wird.

Zusätzlich wird ein orthogonal zum ersten Gitternetz angesetztes, zweites (duales) Gitterzellennetz bzgl. der magnetischen Kantenspannung <math>h_i = H_i \cdot a</math> und des elektrischen Flusses durch die Randflächen <math>d_j = D_j \cdot A_j</math> unter Berücksichtigung der Stetigkeitsbedingungen berechnet.
[[Datei:FIT Grid Integral.svg|mini|200px|Kurvenintegral über die elektrische Feldstärke in einer Gitterzelle]]
Durch die quaderförmige Form der Gitterzellen vereinfacht sich das [[Kurvenintegral|Konturintegral]] der elektrischen Feldstärke <math>\oint_{C(A)}{\vec{E}\cdot d\vec{s}}=\int_{{{C}_{1}}}{\vec{E}\cdot d\vec{s}}+\int_{{{C}_{2}}}{\vec{E}\cdot d\vec{s}}-\int_{{{C}_{3}}}{\vec{E}\cdot d\vec{s}}-\int_{{{C}_{4}}}{\vec{E}\cdot d\vec{s}}</math> zur Summe <math>\sum_{i=1}^4 e_i</math> der Kantenspannungen einer Quaderwand der Gitterzelle.

Die zeitliche Ableitung des magnetischen Flusses durch die Randfläche der Gitterzelle wird nun dieser Summe gleichgesetzt, so dass sich folgende Gleichung ergibt:

:<math>
\sum_{i=1}^4 e_i \;=\; e_1 + e_2 - e_3 - e_4 \;=\; - \frac{\partial}{\partial t} b_n
</math>

Diese Berechnung muss für alle sechs Randflächen einer Gitterzelle wiederholt werden. In Matrixschreibweise ergibt sich das Gleichungssystem

:<math>
{{\begin{pmatrix}
\ldots & \ldots & \ldots & \ldots & \ldots & \ldots & \ldots \\
1 & \ldots & 1 & \ldots & -1 & \ldots & -1 \\
\ldots & \ldots & \ldots & \ldots & \ldots & \ldots &
\ldots
\end{pmatrix}}}
\cdot {{\begin{pmatrix} e_i\\
\vdots \\ e_j \\ \vdots \\ e_k \\ \vdots \\ e_l\end{pmatrix}}}
\;=\; - \frac{\partial}{\partial t} \; {{\begin{pmatrix} \vdots \\b_n
\\ \vdots
\end{pmatrix}}}
</math>

:<math>\textbf{C} \cdot \vec{e} \;=\; - \frac{\partial}{\partial t} \;\vec{b}</math>

Die beschreibende Matrix <math>\textbf{C}</math> besitzt als Elemente nur die Werte 1, 0, −1.

Analog dazu werden die übrigen Maxwell’schen Gleichungen behandelt. In Matrixschreibweise ergibt sich das Gitter-Maxwell-Gleichungssystem

: <math>
\textbf{C} \cdot \vec{e} \;=\; - \frac{\partial}{\partial t} \;\vec{b}
</math>

: <math>
\textbf{C}_{\text{Dual}} \cdot \vec{h} \;=\; \frac{\partial}{\partial t} \; \vec{d} + \vec{j}
</math>

: <math>
\textbf{S}_{\text{Dual}} \cdot \vec{d} \;=\; \vec{q}
</math>

: <math>
\textbf{S} \cdot \vec{b} \;=\; \vec{0}
</math>

Die Matrix <math>\textbf{C}</math> entspricht dem analytischen [[Rotation eines Vektorfeldes|Rotations-Operator]], die Matrix <math>\textbf{S}</math> entspricht dem analytischen [[Divergenz eines Vektorfeldes|Divergenz-Operator]]. Der Index <math>Dual</math> weist auf die Berechnung der Kantenspannungen und Flüsse im Dualen Gitter hin.

Die Materialgleichungen werden analog zu den Maxwell’schen Gleichungen diskretisiert.

:<math>
\vec{d} \;=\; \textbf{M}_{\varepsilon} \;\cdot\; \vec{e}
</math>

:<math>
\vec{b} \;=\; \textbf{M}_{\mu} \;\cdot\; \vec{h}
</math>

:<math>
\vec{j} \;=\; \textbf{M}_{\kappa} \;\cdot\; \vec{e} \;+\; \vec{j}_s,
</math>
wobei die Materialgrößen orts-, frequenz- und richtungsabhängig sein können.

Die FIT-Methode ist auf alle elektromagnetische Probleme im Zeit- und Frequenzbereich anwendbar, sowohl in der Elektrostatik, als auch in der Elektrodynamik. Durch den speziellen Zuschnitt der FIT-Methode auf die Maxwell’schen Gleichungen und das daraus entstehende diskrete Analogon sind die [[Stetige Funktion|Stetigkeit]]sbedingungen a priori erfüllt und die analytischen Eigenschaften der Vektoroperationen werden beibehalten.

Für elektrodynamische Probleme werden im Frequenzbereich alle zeitlichen Ableitungen durch <math>j\omega</math> ersetzt. Das Ergebnis einer Simulation im Frequenzbereich liefert die Impulsantwort auf ein monofrequentes Eingangssignal.

Im Zeitbereich ist eine breitbandige Anregung mit freien Signalverläufen gestattet. Die Simulationsrechnung beschreibt in diesem Fall das Frequenzverhalten über einen vorab definierten Frequenzbereich.

== Literatur ==

*T. Weiland: ''Eine Methode zur Lösung der Maxwellschen Gleichungen für sechskomponentige Felder auf diskreter Basis'', AEÜ, Band 31, Heft 3, pp. 116–120, 1977
*T. Weiland: ''A Discretization Method for the Solution of Maxwell’s Equations for Six-Component Fields'', Electronics and Communications AEUE, vol. 31, no. 3, pp. 116–120, 1977.

[[Kategorie:Theoretische Elektrotechnik]]
[[Kategorie:Numerische Mathematik]]
[[Kategorie:Elektrodynamik]]

Finite-Integral-Methode - Versionsgeschichte

imported>Karsten11 am 24. April 2023 um 13:48 Uhr