Phase antenna center (FCA) and its search in Ansys HFSS

When calculating the high-frequency technique with the use of mirror reflecting systems (parabolic mirrors), there is always the problem of searching for the phase center of the antenna (FCA), since proper operation of the mirror is possible only if the focus is on an antenna (called the feed, feeder, feedhorn) which has a phase wave front in the form of a sphere, and the center of this sphere is in the focus of the mirror. With any deviations, both the shape of the phase front from the sphere, and the displacement of the FCA from the focus of the mirror - the efficiency of the mirror system falls, because its radiation pattern is distorted. 3r33333.

3r33333.

Although the topic of the FCA search is quite relevant even in everyday life, because apart from traditional satellite television antennas, parabolic antennas for WiFi, WiMAX and cellular (UMTS /3G, LTE /4G) have become common - nevertheless, this topic is poorly lit in the literature and users often confuse phase DN with the usual radiation pattern. 3r33333.

3r33333.

In videos about computer simulation programs, you can sometimes find practical instructions on how to search for FCA, but usually there is not even a minimal explanation of what we are looking for and what we get. 3r33333.

3r33333.

Therefore, in order to fill the gap, we will write a small article with practical examples.

3r33333.

3r33333. ** Phase radiation pattern ** - is the dependence of the phase of the electromagnetic field emitted by the antenna on the angular coordinates. 3r33333.

(A.P. Pudovkin, Yu.N. Panasyuk, A.A. Ivankov - 3—3–350. Fundamentals of the theory of antennas 3–3–351.) 3–3–3369.

3r33333.

Since in the far-field zone of the antenna, the field vectors E and H are in phase, then the phase DN also equally relates to the electrical and magnetic components of the electromagnetic field emitted by the antenna. 3r33333.

Phase DN is denoted by the following letter Psi:

** Ψ = Ψ (θ, φ) ** , with r = const. 3r33333.

3r33333.

If Ψ (θ, φ) = const with r = const, then this means that the antenna forms a phase wave front in the form of a sphere. 3r33333.

The center of this sphere, in which the origin of the coordinate system is located, is called the phase center of the antenna (FCA). 3r33333.

3r33333.

The phase center of the antenna is the point at which a single emitter of a spherical wave can be placed, equivalent to the antenna system under consideration with respect to the phase of the field created. 3r33333.

(Drabkin A.L., Zuzenko V.L. 3-3-350. Antenna-feeder devices 3-3-351.) 3-33369.

3r33333.

FCs do not have all antennas. For antennas with a phase center and a multilobe amplitude DN with distinct zeros between them, the field phase in adjacent lobes differs by π (180 °). 3r33333.

3r33333.

The relationship between the amplitude and phase patterns of the same antenna is illustrated in 3r-3369.

3r362. 3r33333.

3r33333.

In real antennas, the phase center is usually considered within the limited angles of the main lobe of the radiation pattern. The position of the phase center depends on the frequency of the signal used, the direction of radiation /reception of the antenna, its polarization and other factors. Some antennas do not have a phase center in the conventional sense. 3r33333.

3r33333.

In the simplest cases, for example, with a parabolic antenna, the phase center coincides with the focus of the paraboloid and can be determined from geometrical considerations. In more complex cases, such as horn antennas, the position of the phase center is not obvious and requires appropriate measurements. 3r33333.

3r33333.

Field measurements of the phase center are very laborious (especially in a wide frequency band). 3r33333.

In CAD-simulators of electromagnetic fields, the calculation of the FCA is a very simple task, but it still requires several manual manipulations, since it is executed by brute force and requires a small initial setup of the function that we are going to brute force. 3r33333.

3r33333.

For practical calculations, let's take a real parabolic feed for the Ku-band - LNB range from Inverto, the Black Ultra series. 3r33333.

This feeder has the following form (in section)

3r33333.

A pea-sized ball will be the FCA, but we still don’t know it and our task is to find its position. 3r33333.

In the example we will use such input:

- The calculation frequency is 11538.5 MHz (wavelength is ??? mm)

- linear horizontal polarization (in the Y axis)

- the antenna itself is directed along the X axis, i.e. the main radiation direction is θ = 9? φ = 0

3r33333.

The calculation of the traditional parameters of the far field (Far Field) in Ansys HFSS gives such a pattern in 3D and 2D

3r3102. 3r33333.

3r3105. 3r33333.

3r33333.

The instantaneous values of the voltage (Volt /meter) of the electric field (E-field) depending on the phase

3r33112. 3r33333.

3r33333.

Integral E-field strength (for> 1 wave rotation)

3r3119. 3r33333.

3r33333.

All such parameters of the far field (Far-Field) both in field measurements and in CAD-simulations are calculated on an infinite sphere - Infinite Sphere. The test antenna or its computer model is placed in the center of such a sphere, and the measuring probe moves along the perimeter of such a sphere and measures the amplitude, polarization (amplitude of one of the components) and the phase of the EM wave. The probe can be fixed permanently and rotate the test antenna. The main thing to:

- the distance was always the same (that is, it was the measuring sphere)

- the radius of the sphere was large enough (at least 5-10 lambd), so that measurements were carried out only in that region of space where the electric field vectors E and magnetic H are in phase, i.e. none of the components prevails and is not shifted in phase (has no reactivity) at the expense of charge carriers that are in the metal conductors of the antenna or at the expense of charged dielectric molecules. 3r33333.

3r33333.

In ** Ansys HFSS ** to measure the far field, you must create at least one infinite sphere: Radiation -> Insert Far Field Setup -> Infinite Sphere

3r33333.

3r33333.

φ and θ can always be specified from 0 to 36? but in order to save time on calculations it is sometimes rational to limit the angle of interest to a certain sector. When specifying a step of 1 degree, the full sphere will occupy 360 * 360 = ??? design points, and at a step of 0.1 degrees almost 13 million. To create 3D /2D reports of the pattern, a step of 2-3 degrees is usually sufficient (1?400 design points at a step 3 degrees). Step 1 degree or less, it makes sense to use only for the analysis of the slice

3r3143. 3r33333.

3r33333.

In the “Coordinate System” tab, each sphere necessarily has its own center of coordinates. By default, there is always a global coordinate center of the project[0, 0, 0]. If desired, you can add any number of other relative coordinates. As elements of the geometry of the model and the user sphere "Infinite Sphere" can be assigned relative to the global center of coordinates or relative to the user. We will use this below. 3r33333.

3r33333.

The diverging phase wave front was visible on the E-field animation above. The EM wave forms concentric circles, similar to circles on the water from an abandoned stone. The phase center is the point at which such a stone was thrown. It can be seen that his position is somewhere in the mouth of the horn, but its exact position is not obvious. 3r33333.

3r33333.

The FCA search method is based on the fact that we are looking at the direction of the E-field vector (its phase) along the surface of an infinitely distant sphere. 3r33333.

3r33333.

For the demonstration, we will create 2 animations with vectors of the E-field on a sphere with a radius of 4 lambda (this is not an infinite sphere, but for a better scale of drawing such a radius is enough). 3r33333.

3r33333.

In the first animation, the center of the sphere is placed exactly in the FCA

3r33333.

3r3168. 3r33333.

3r33333.

In the second animation, the center is located at the point of project ? ? 0 (looking ahead to say that it is ??? mm behind the FCA) 3r33369.

3r33175. 3r33333.

3r33333.

On the surface of the first sphere (it is a curve, it is not a plane), it can be seen that the vectors move synchronously. Their amplitude (magnitude) is different, because the antenna's antenna has a maximum in the center (up to 14.4 dBi) which gradually fades by 2 times (-3 dB) at angles of ± 20 °. 3r33333.

We are not interested in color /length, but in the direction of the vector. So that they all move synchronously (in phase). 3r33333.

3r33333.

In the first animation, all the vectors move synchronously, as if the ball rotates right and left. 3r33333.

In the second animation, the vectors are asynchronous, some have already changed the direction of movement, others are not yet. The surface of this sphere constantly undergoes surface tension /deformation. 3r33333.

3r33333.

The first sphere is in the FCA, the second is not in the FCA. 3r33333.

3r33333.

The task of finding the FCA using this method is to move (brutally) the Infinite Sphere with a small step until the phase spread in the area of interest to us (we are only interested in the main radiation lobe) becomes minimal (ideally zero). 3r33333.

3r33333.

But before moving on to brute force, we first understand how it is possible to display phase NAM in HFSS. 3r33333.

3r33333.

In the “Results -> Create Far Field Report” far-field reports, we can derive either a traditional rectangular plot (Rectangular plot) or a 2D pie chart (Radiation pattern) where along a single axis (for example, X) we derive the dependence of the angular coordinate (for example, θ), and Y axis - phase values at these angles θ. 3r33333.

3r33333.

The report we need is the rE - “radiated (radiated) E field” 3r-3369.

For each angle[φ, θ]on the infinite sphere, the complex number (vector) of the electric field is calculated. 3r33333.

When building ordinary amplitude graphs (radiation pattern, distribution of radiation power in a direction), we are interested in the amplitude (mag) of this field, which can be obtained either as mag (rE) or immediately using the more convenient variable Gain (power is relative to the power at the excitation port and isotropic radiator). 3r33333.

When constructing the phase DN, we are interested in the imaginary part of the complex number (vector phase) in polar notation (in degrees). To do this, use the mathematical function ang_deg (angle____ degrees) or cang_deg (accumulated_-angle____ degrees)

3r33333.

For the LNA Inverto Black Ultra antenna, the phase DN in the XZ plane (φ = 0) with horizontal excitation polarization (rEY) has the following form 3r3369.

3r33333.

3r33333.

3r33333.

Theta = 90 angle is radiation forward, Theta = 0 up, Theta = 180 down. 3r33333.

The values of ** ang_deg ** vary from -180 to +18? the angle of 181 ° is the angle of -179 °, so the graph has the shape of a saw when passing through points ± 180 °. 3r33333.

The values of ** cang_deg ** accumulate if the direction of the phase change is constant. If the phase has made up to 3 full turns (6 times crossed 180 °), then the accumulated value reaches 1070 °. 3r33333.

3r33333.

As it was written in the beginning of the article, the phase and amplitude DN at antennas are usually associated with one another. In the adjacent amplitude lobes (beam), the phases differ by 180 °. 3r33333.

3r33333.

We put one on top of another the graphs of the phase (red /green) and amplitude (purple) DN

3r33333.

3r33333.

Humps on the amplitude NAM clearly follow the phase fractures, as it is written in the books. 3r33333.

3r33333.

We are interested in the phase front only in a certain sector of space, within the main radiation lobe (the remaining petals still shine past the parabolic mirror). 3r33333.

3r33333.

Therefore, we restrict the schedule only to the sector ± 45 °. 3r33333.

3r33333.

3r33333.

Add on the chart the markers MIN (m1) and MAX (m2) which show the greatest phase variation in the studied sector. 3r33333.

In addition, we add the mathematical function pk2pk () which automatically searches for the entire graphics minimum and maximum and shows the difference. 3r33333.

The graph above shows the difference m2-m1 = pk2pk = ??? °

3r33333.

The task of the FCA search is to move with the small step of the Infinite Sphere until the value of the pk2pk function (cang_deg (rE)) is minimized. 3r33333.

3r33333.

To move the Infinite Sphere, you need to create another additional coordinate system: Modeler -> Coordinate System -> Create -> Relative CS -> Offset 3r33333.

3r33333.

since we obviously know that the symmetrical horn of the FCA will be on the X axis (Z = Y = 0), then for Z and Y we set ? and it will only move along the X axis, for which we assign the variable Pos (with an initial value of 0 mm)

3r33333.

3r33333.

To automate the brute-force process, let's create an optimization task. 3r33333.

Optimetrics -> Add -> Parametric, and set the step of the variable Pos 1 mm, in the range from 0 to 100 mm

3r33333.

3r33333.

In the tab “Calculations -> Setup Calculation” we select the report type “Far Field” and the function pk2pk (cang_deg (rEY)). In the “Range Functions” button, we indicate the range from -45 to +45 degrees (or any other interesting)

3r33232. 3r33333.

3r33333.

Run ParametricSetup1 -> Analyze. 3r33333.

The calculation is performed fairly quickly, because All far-field calculations are related to post-processing and do not need to be re-solved by the model. 3r33333.

3r33333.

After the calculation is completed, click ParametricSetup1 -> View analysis results. 3r33333.

3r33312. 3r33333.

3r33333.

We see a clear minimum at a distance of X = 25mm

3r33333.

For higher accuracy, we edit the parametric analysis in the range of 25.0-25.1 mm in increments of ??? mm

3r33333.

3r33333.

We get a clear minimum at X = ??? mm 3r33333.

3r33333.

To visually assess where in the model the FCA turned out, you can draw spheres (Non-model) or points. 3r33333.

Here at the point X = ??? mm 2 spheres (with radius 2 and 4 lambda)

are placed.

3r33333.

3r33338. 3r33333.

3r33333.

Here is the same, in the animation

3r33333. 3r33333.

3r33333.

Here is a larger plan of the plane and the pea at the point X = ???r3r3369.

3r33352. 3r33333.

3r33333.

It is widely believed that in HFSS (and other programs, such as CST), when “3D Plot” graphics are superimposed on the antenna geometryWhich schedule is automatically placed in the FCA. 3r33333.

3r33333.

Unfortunately, this is not the case. The 3D graphic is always superimposed on the center of the coordinate system that was used when setting the “Infinite Sphere” for this graphic. If the global coordinate system was used by default[0, 0, 0], the 3D Plot will be placed at ??? (even if the antenna itself is far away). 3r33333.

In order to combine graphics, in the settings of 3D Plot you need to select such an “Infinite Sphere” (create another one), for which “Relative CS” is set at the FCA point which we found manually. 3r33333.

3r33333.

It should be noted that such an overlay will be truthful only for the sector under study (for example, the main beam of the DN); in the side and rear petals of the FC, it may be in a different place or be non-spherical. 3r33333.

3r33333.

Also note that the “Infinite Sphere” settings have nothing to do with the “Radiation Boundary” boundary condition. The Rad layer can be specified as a rectangle, a cone, a cylinder, a ball, an ellipsoid of rotation, and somehow move its position, shape, and rotation. The position and shape of “Infinite Sphere” will not change. This will always be a sphere (ball) with an infinite (rather large) radius and with a center in a given coordinate system.

It may be interesting

This is my first time i visit here. I found such a substantial number of interesting stuff in your blog especially its examination. Really its inconceivable article. Keep it up.Gulf Coast Western Reviews

Thankyou for sharing the data which is beneficial for me and others likewise to see.Gulf Coast Western Reviews

Positive site, where did u come up with the information on this posting?I have read a few of the articles on your website now, and I really like your style. Thanks a million and please keep up the effective work.GSM Solutions