Difference between revisions of "Lecture 4."
(→Method of Moments (MoM)CVEL - Electromagnetic Modeling (MoM)) |
(→Basis and Weighting Function) |
||
| Line 68: | Line 68: | ||
==== Basis and Weighting Function ==== | ==== Basis and Weighting Function ==== | ||
<blockquote> | <blockquote> | ||
| − | The most common equation form solved by CEM modeling codes based on the Method of Moments is the Electric Field Integral Equation (EFIE). This is an equation of the form, | + | The most common equation form solved by CEM modeling codes based on the Method of Moments is the ''Electric Field Integral Equation'' (EFIE). This is an equation of the form, |
<math>\vec{E}=f_e(\vec{J},\vec{M})</math>, | <math>\vec{E}=f_e(\vec{J},\vec{M})</math>, | ||
| − | where E is the impressed (i.e. source) electric field and J and M are the induced electric and magnetic current densities, respectively. The EFIE will be discussed further in the section describing the Boundary Element Method. Generally, codes that solve a form of the EFIE excel at modeling open (unbounded) geometries in which the electric field dominates in the near-field region of the source. | + | where <math>\vec{E}</math> is the impressed (i.e. source) electric field and <math>\vec{J}</math> and <math>\vec{M}</math> are the induced electric and magnetic current densities, respectively. The EFIE will be discussed further in the section describing the Boundary Element Method. Generally, codes that solve a form of the EFIE excel at modeling open (unbounded) geometries in which the electric field dominates in the near-field region of the source. |
| − | Another equation solved by Moment Method codes is the Magnetic Field Integral Equation (MFIE), which has the general form, | + | Another equation solved by Moment Method codes is the ''Magnetic Field Integral Equation'' (MFIE), which has the general form, |
| − | <math></math> | + | <math>\vec{H}=f_m(\vec{J},\vec{M})</math>, |
| − | where H is the impressed (i.e. source) magnetic field intensity. Codes that solve a form of the MFIE are best suited for modeling geometries with circulating currents, where the magnetic near field is dominant. | + | where <math>\vec{H}</math> is the impressed (i.e. source) magnetic field intensity. Codes that solve a form of the MFIE are best suited for modeling geometries with circulating currents, where the magnetic near field is dominant. |
Moment Method codes based on the EFIE or MFIE alone, may exhibit unstable behavior when the modeling surfaces form a resonant cavity at a particular frequency. To avoid this, many moment method codes solve a linear combination of the EFIE and MFIE known as a Combined Field Integral Equation (CFIE). This requires more calculations to fill the matrix, but results in a more stable solution when the modeling surface is large enough to support an interior resonance. | Moment Method codes based on the EFIE or MFIE alone, may exhibit unstable behavior when the modeling surfaces form a resonant cavity at a particular frequency. To avoid this, many moment method codes solve a linear combination of the EFIE and MFIE known as a Combined Field Integral Equation (CFIE). This requires more calculations to fill the matrix, but results in a more stable solution when the modeling surface is large enough to support an interior resonance. | ||
| Line 84: | Line 84: | ||
Some CEM modeling codes employ the Method of Moments to solve other equations. For example, static modeling codes often solve a form of Laplace's equation relating electric field strengths to charge densities or magnetic field strengths to current densities. The Generalized Multiple Technique (GMT), which is described in another section of this report, employs a moment method to solve equations for the electric field generated by multipole sources. | Some CEM modeling codes employ the Method of Moments to solve other equations. For example, static modeling codes often solve a form of Laplace's equation relating electric field strengths to charge densities or magnetic field strengths to current densities. The Generalized Multiple Technique (GMT), which is described in another section of this report, employs a moment method to solve equations for the electric field generated by multipole sources. | ||
</blockquote> | </blockquote> | ||
| − | |||
==== Basis and Weighting Function ==== | ==== Basis and Weighting Function ==== | ||
Revision as of 08:24, 23 February 2019
|
Lecture 4. | |
|
|
|
| Animated cut through diagram of a typical fuel injector [1] | Animated cut through diagram of a typical fuel injector [Click to see animation]. |
|
Oktató
|
További oktatók:
|
Contents
A feladat célja
A hallgató megismerje a végeselem-módszer főbb lépéseit, mint a modell előkészítése (geometria elkészítése vagy importálása), anyagparaméterek, peremfeltételek és gerjesztés megadása az üzemanyag befecskendező szolenoidjának szimulációján keresztül.
A feladat megoldásához szükséges ismeretek
- A végeselem-módszer lépései;
- A sztatikus mágneses térre vonatkozó elméleti ismeretek (anyagok definiálásához, gerjesztés megadásához);
- A geometria elkészítéséhez CAD rendszer ismerete.
Boundary Element Method (BEM)[2]
Method of Moments (MoM)[3]
In the 1960s, R.F. Harrington and others applied a technique called the Method of Moments to the solution of electromagnetic field problems. The Method of Moments (also called the Method of Weighted Residuals) is a technique for solving linear equations of the form,
[math]\mathcal{L}(\phi)=f[/math]
where the functional [math]\mathcal{L}(\bullet)[/math] is a linear operator, [math]f[/math] is a known excitation or forcing function, and [math]\phi[/math] is an unknown quantity. To solve this problem on a digital computer, we start by expressing the unknown solution as a series of basis or expansion functions, [math]v_n[/math],
[math]\phi = \sum_{n=1}^{N}a_n v_n[/math]
where [math]a_n[/math] are unknown coefficients describing the amplitude of each term in the series.
Now instead of one equation with a continuous unknown quantity, [math]f[/math], we have an equation with [math]N[/math] scalar unknowns,
[math]\mathcal{L}(a_1v_1+a_2v_2+\dots+a_Nv_N)=f[/math].
To solve for the values of [math]a_n[/math], we need [math]N[/math] linearly independent equations; so [math]a_n[/math] different weighting or testing functions, [math]w_n[/math], are applied. This yields the following system of [math]N[/math] equations in [math]N[/math] unknowns:
This linear system of equations has the form,
[math]\textbf{A}\textbf{x}=\textbf{b}[/math]
where the elements of [math]\textbf{A}[/math] are known quantities that can be calculated from the linear operator, the functional [math]\mathcal{L}(\bullet)[/math], and the chosen basis and weighting functions. The elements of \textbf{b} are determined by applying the weighting functions to the known forcing function. The unknown elements of [math]\textbf{x}[/math] can be found by solving the matrix equation. After solving for [math]\textbf{x}[/math] (i.e. the unknown coefficients, [math]a_n[/math]), the value of [math]\phi[/math] is determined using [math]\phi = \sum_{n=1}^{N}a_n v_n[/math] equation.
The Method of Moments (MoM) can be used to solve a wide range of equations involving linear operations including integral and differential equations. This numerical technique has many applications other than electromagnetic modeling; however the MoM is widely used to solve equations derived from Maxwell's equations. In general, moment method codes generate and solve large, dense matrix equations and most of the computational resources required are devoted to filling and solving this matrix equation. The particular form of the equations that is solved and the choice of basis and weighting functions have a great impact on the size of this matrix and ultimately the suitability of a given moment method code to model a given geometry.
Basis and Weighting Function
The most common equation form solved by CEM modeling codes based on the Method of Moments is the Electric Field Integral Equation (EFIE). This is an equation of the form,
[math]\vec{E}=f_e(\vec{J},\vec{M})[/math],
where [math]\vec{E}[/math] is the impressed (i.e. source) electric field and [math]\vec{J}[/math] and [math]\vec{M}[/math] are the induced electric and magnetic current densities, respectively. The EFIE will be discussed further in the section describing the Boundary Element Method. Generally, codes that solve a form of the EFIE excel at modeling open (unbounded) geometries in which the electric field dominates in the near-field region of the source.
Another equation solved by Moment Method codes is the Magnetic Field Integral Equation (MFIE), which has the general form,
[math]\vec{H}=f_m(\vec{J},\vec{M})[/math],
where [math]\vec{H}[/math] is the impressed (i.e. source) magnetic field intensity. Codes that solve a form of the MFIE are best suited for modeling geometries with circulating currents, where the magnetic near field is dominant.
Moment Method codes based on the EFIE or MFIE alone, may exhibit unstable behavior when the modeling surfaces form a resonant cavity at a particular frequency. To avoid this, many moment method codes solve a linear combination of the EFIE and MFIE known as a Combined Field Integral Equation (CFIE). This requires more calculations to fill the matrix, but results in a more stable solution when the modeling surface is large enough to support an interior resonance.
Some CEM modeling codes employ the Method of Moments to solve other equations. For example, static modeling codes often solve a form of Laplace's equation relating electric field strengths to charge densities or magnetic field strengths to current densities. The Generalized Multiple Technique (GMT), which is described in another section of this report, employs a moment method to solve equations for the electric field generated by multipole sources.
Basis and Weighting Function
An appropriate choice of basis and weighting functions can make a tremendous difference in the number of elements, [math]N[/math], required to obtain an accurate solution. Since the solution is represented as a summation of basis functions, it is important to choose basis functions that accurately represent the solution with a small number of terms. For example, when solving for the current distribution on a surface, the basis functions should be current distribution elements that can be summed together in a way that is able to efficiently approximate any overall current distribution that might result from the analysis.
Weighting functions should be chosen that maximize the linear independence of the various weighted forms of the equation. Often, the best choices of weighting functions are functions that are identical to the basis functions. Moment method techniques that employ identical basis and weighting functions are called Galerkin techniques.
