CN108021757B - Design method for high-power microwave phase inversion - Google Patents

Design method for high-power microwave phase inversion Download PDF

Info

Publication number
CN108021757B
CN108021757B CN201711282811.5A CN201711282811A CN108021757B CN 108021757 B CN108021757 B CN 108021757B CN 201711282811 A CN201711282811 A CN 201711282811A CN 108021757 B CN108021757 B CN 108021757B
Authority
CN
China
Prior art keywords
electric field
phase
calculated
field amplitude
magnetic field
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711282811.5A
Other languages
Chinese (zh)
Other versions
CN108021757A (en
Inventor
刘建卫
孙旭
张莉琴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201711282811.5A priority Critical patent/CN108021757B/en
Publication of CN108021757A publication Critical patent/CN108021757A/en
Application granted granted Critical
Publication of CN108021757B publication Critical patent/CN108021757B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Inorganic Compounds Of Heavy Metals (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention discloses a design method for high-power microwave phase inversion, wherein high-power microwaves sequentially pass through spatial positions X1, X2 and X3 during spatial transmission and are measured according to an infrared camera at the X1 position
Figure DDA0001497792750000011
Assumption at the beginning of iteration
Figure DDA0001497792750000012
Phase phi of1Zero, magnetic field amplitude
Figure DDA0001497792750000013
Calculating the electric field amplitude Ey of the X2 position by using the Stratton-Chu formula2And phase phi2(ii) a The electric field amplitudes and phases of the X2 and X3 positions are sequentially calculated in the above manner, and the X3 position transmits the electric field amplitude Ey calculated to the X1 position4And phase phi4Completing an iteration when the calculated phase phi1，Φ2，Φ3When the phase of the wave beam at the X1, X2 and X3 positions is not changed, the phase of the electromagnetic wave is obtained by inverting the wave beam by using the vector diffraction theory, the method is applicable to the frequency of the electromagnetic wave transversely output by any high-power cyclotron oscillation tube, and is more accurate than the current design method for phase inversion by using the scalar diffraction theory.

Description

Design method for high-power microwave phase inversion
Technical Field
The invention relates to the technical field of electronics, in particular to a design method for high-power microwave phase inversion.
Background
When high-power microwaves (generally, quasi-gaussian beams in a propagation mode) transversely output from an output window of a cyclotron oscillation tube are propagated in space, in order to better correct beam parameters, a subsequent reflecting surface (such as an optical matching unit (MOU)), particularly a phase correction surface, needs to be designed, and electromagnetic field amplitude and phase distribution information of an output beam need to be accurately known. The amplitude and phase information of the electromagnetic field is the basis for designing a subsequent mirror or corrugated waveguide transmission device, and has extremely important function.
The high-power microwave is transmitted in space, and the radiation field of the space is calculated by adopting a scalar diffraction theory method or a vector diffraction theory method. Since the scalar diffraction theory is always adopted abroad to calculate the space radiation field, the phase of a beam at a certain position can be obtained by adopting a phase inversion method based on the scalar diffraction theory. The Russian academy of sciences applied physical research Institute (IAP) published in 1995 "3D wave field retrieval from sensitive research in a raw cross sections" (A.V Chirkov, G.G Dennisov, N.LAleksandov) (Optics Communications, Volume 115, Issues 5-6, 1April, Pages 449-. Up to now, although the radiation field of the space can be calculated by the method of the vector diffraction theory, the phase inversion method based on the vector diffraction theory is not adopted. The vector diffraction theory method is a three-dimensional full-wave electromagnetic wave analysis method, and six components of a radiation field which can be calculated are as follows: expressed as Ex, Ey, Ez, Hx, Hy, Hz in a rectangular coordinate system, the calculated radiation field is more accurate than that calculated by scalar diffraction theory, and more information of the electromagnetic field can be known.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a design method of high-power microwave phase inversion, which opens up a design method of high-power microwave phase inversion adopting a vector diffraction theory, is more accurate than the current design method of phase inversion adopting a scalar diffraction theory, and can effectively solve the problems provided by the background technology.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a design method of high-power microwave phase inversion comprises the following steps:
step 100, the high-power microwave sequentially passes through spatial positions X1, X2 and X3 during spatial transmission, and the electric field amplitude measured by an infrared camera at the X1 position is measured
Figure BDA0001497792730000021
Assumptions at the beginning of the iteration
Figure BDA0001497792730000022
Phase phi of1Zero, magnetic field amplitude
Figure BDA0001497792730000023
Wherein, η0377 omega is free space wave impedance, and the electric field amplitude Ey of the X2 position is calculated by the formula of Stratton-Chu2And phase phi2；
Step 200, Ey to be calculated2Measuring the electric field amplitude at X2 position with infrared camera
Figure BDA0001497792730000027
Instead, the phase phi2The constant, magnetic field amplitude at the X2 position is shown as
Figure BDA0001497792730000024
Wherein, η0377 omega is free space wave impedance, and then the electric field amplitude Ey of the X3 position is calculated by the formula of Stratton-Chu3And phase phi3；
Step 300, Ey to be calculated3Measuring the electric field amplitude at X3 position with infrared camera
Figure BDA0001497792730000028
Instead, the phase is unchanged, the magnetic field amplitude at the X3 position is
Figure BDA0001497792730000025
Wherein, η0377 omega is free space wave impedance, and then the electric field amplitude Ey to the X1 position is calculated through the back transmission of the X3 position4And phase phi4；
Step 400, Ey to be calculated4Measuring the electric field amplitude at X1 position with infrared camera
Figure BDA0001497792730000029
Instead, the phase phi4Constant, magnetic field amplitude
Figure BDA0001497792730000026
Wherein, η0377 omega is free space wave impedance, thus completing one iteration;
step 500, when the calculated phase phi1，Φ2，Φ3When the change is not changed, the phases of the beams at the final X1, X2 and X3 positions are obtained.
Furthermore, thermal paper is arranged at the spatial positions X1, X2 and X3, and burnt spots are formed on the thermal paper by high-power microwaves, so that potential diagrams of energy and the like at the spatial positions X1, X2 and X3 can be measured by an infrared camera.
Further, the Stratton-Chu formula:
Figure BDA0001497792730000031
where μ is the vacuum permeability, ε is the vacuum dielectric constant, ω is the microwave angular frequency, k0Is the wave number in the free space of the electromagnetic wave,
Figure BDA0001497792730000032
and
Figure BDA0001497792730000033
respectively expressed as the electric field intensity and the magnetic field intensity of a certain plane in space, S is the area of a certain plane,
Figure BDA0001497792730000034
is the outer normal direction of a certain plane,
Figure BDA0001497792730000035
is a free space Green function, R is the distance from a point on a certain plane to a space point,
Figure BDA0001497792730000036
respectively, the electric field strength and the magnetic field strength radiated to a point in space in a certain plane.
Further, the Stretton-Chu formula can be used for controlling the electric field intensity and the magnetic field intensity in a rectangular coordinate systemRespectively expressed as:
Figure BDA0001497792730000037
wherein Ex, Ey and Ez are respectively x, y and z components of an electric field, and Hx and Hy Hz are respectively x, y and z components of a magnetic field;
because Ex is 0, Ey is in the infrared camera measurement, Ez is 0, Hx is 0, Ey is 0, Ez is Ey/377, the magnitude of Ey and its phase phi on the next plane can be found out according to the formula straton-Chu.
Compared with the prior art, the invention has the beneficial effects that:
the invention opens up a design method for high-power microwave phase inversion by adopting the vector diffraction theory, obtains the phase of the electromagnetic wave by using the vector diffraction theory, is suitable for the frequency of the electromagnetic wave transversely output by any high-power cyclotron oscillation tube, and is more accurate than the current design method for phase inversion by adopting the scalar diffraction theory.
Drawings
FIG. 1 is a flow chart of a design method of the present invention;
FIG. 2 is a schematic diagram of the propagation of high power microwaves in space according to the present invention;
FIG. 3 is an equipotential diagram at the X1 position of the present invention;
FIG. 4 is an equipotential diagram at the X2 position of the present invention;
FIG. 5 is an equipotential diagram at the X3 position of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the present invention provides a design method of high power microwave phase inversion, which comprises the following steps:
step 100, high-power microwave is in the airThe electric field amplitude measured by an infrared camera at the X1 position passes through the space positions X1, X2 and X3 in turn during time-transfer
Figure BDA0001497792730000041
Assumptions at the beginning of the iteration
Figure BDA0001497792730000042
Phase phi of1Zero, magnetic field amplitude
Figure BDA0001497792730000043
Wherein, η0377 omega is free space wave impedance, and the electric field amplitude Ey of the X2 position is calculated by the formula of Stratton-Chu2And phase phi2；
Step 200, Ey to be calculated2Measuring the electric field amplitude at X2 position with infrared camera
Figure BDA0001497792730000044
Instead, the phase phi2The constant, magnetic field amplitude at the X2 position is shown as
Figure BDA0001497792730000045
Wherein, η0377 omega is free space wave impedance, and then the electric field amplitude Ey of the X3 position is calculated by the formula of Stratton-Chu3And phase phi3；
Step 300, Ey to be calculated3Measuring the electric field amplitude at X3 position with infrared camera
Figure BDA0001497792730000046
Instead, the phase is unchanged, the magnetic field amplitude at the X3 position is
Figure BDA0001497792730000051
Wherein, η0377 omega is free space wave impedance, and then the electric field amplitude Ey to the X1 position is calculated through the back transmission of the X3 position4And phase phi4；
Step 400, Ey to be calculated4Measuring the power at the X1 position with an infrared cameraMagnitude of field
Figure BDA0001497792730000057
Instead, the phase phi4Constant, magnetic field amplitude
Figure BDA0001497792730000052
Wherein, η0377 omega is free space wave impedance, thus completing one iteration;
step 500, when the calculated phase phi1，Φ2，Φ3When the change is not changed, the phases of the beams at the final X1, X2 and X3 positions are obtained.
As shown in fig. 2, which is a schematic diagram of a path of high-power microwaves transmitted in space, information about amplitude and phase distribution of electromagnetic fields at spatial positions X1, X2, and X3 is to be obtained, wherein X1, X2, and X3 may be equally spaced.
For the amplitude of the electromagnetic field at a certain position, the measurement of low power can be carried out, namely, thermal paper is arranged at the corresponding position, burning spots are punched on the thermal paper, energy equipotential diagrams at X1, X2 and X3 can be measured through an infrared camera, the equipotential diagrams of the amplitude of a beam can be known according to the energy equipotential diagrams, and the amplitude at X1, X2 and X3 can be actually measured through the measurement of the amplitude
Figure BDA0001497792730000053
Is shown, according to electromagnetic field theory
Figure BDA0001497792730000054
Figure BDA0001497792730000055
η0The amplitude distribution of the magnetic field of the electromagnetic wave at X1, X2, and X3 is obtained by assuming 377 Ω as free space wave impedance. And calculating the electromagnetic field distribution of quasi-Gaussian beams output by the high-power cyclotron oscillation tube in space by adopting a vector diffraction theory, and only considering the Ey component and the Hz component of the electromagnetic field under the establishment of a rectangular coordinate system.
The vector diffraction theory can be expressed by the Stratton-Chu equation:
Figure BDA0001497792730000056
where μ is the vacuum permeability, ε is the vacuum dielectric constant, ω is the microwave angular frequency, k0Is the wave number in the free space of the electromagnetic wave,
Figure BDA0001497792730000061
and
Figure BDA0001497792730000062
respectively expressed as the electric field intensity and the magnetic field intensity of a certain plane in space, S is the area of a certain plane,
Figure BDA0001497792730000063
is the outer normal direction of a certain plane,
Figure BDA0001497792730000064
is a free space Green function, R is the distance from a point on a certain plane to a space point,
Figure BDA0001497792730000065
respectively, the electric field strength and the magnetic field strength radiated to a point in space in a certain plane.
The electric field strength and the magnetic field strength of the Stratton-Chu formula in a rectangular coordinate system can be respectively expressed as:
Figure BDA0001497792730000066
the Ex, Ey and Ez are x, y and z components of an electric field respectively, and the Hx and Hy Hz are x, y and z components of a magnetic field respectively, because Ex is 0, Ey is in infrared camera measurement, Ez is 0, Hx is 0, Ey is 0, Ez is Ey/377, and the Ey amplitude and the phase phi on the next plane can be obtained according to the Stratton-Chu formula.
In the embodiment, by taking 94GHz Gaussian beams with the beam waist radius of 8mm as an example, the transmission ranges of the electric fields at X1, X2 and X3 positions can be measured by an infrared camera by transmitting 100mm (X1 position), 200mm (X2 position) and 300mm (X3 position) in space
Figure BDA0001497792730000067
The magnetic field amplitude can be obtained by the formula (1):
Figure BDA0001497792730000068
according to the vector diffraction theory, the value of the electromagnetic field is substituted into the Stratton-Chu formula, and the iterative algorithm of the figure 1 is combined, so that phi can be obtained1，Φ2，Φ3The values of (A) are shown in an equipotential chart in FIG. 3, FIG. 4, and FIG. 5.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.

Claims (3)

1. A design method of high-power microwave phase inversion is characterized by comprising the following steps: the method comprises the following steps:
step 100, the high-power microwave sequentially passes through spatial positions X1, X2 and X3 during spatial transmission, and the electric field amplitude measured by an infrared camera at the X1 position is measured
Figure FDA0002292735850000011
Assumptions at the beginning of the iteration
Figure FDA0002292735850000012
Phase phi of1Zero, magnetic field amplitude
Figure FDA0002292735850000013
Wherein, η0377 omega is free space wave impedance, and the X2 position is obtained by calculation through the Stratton-Chu formulaElectric field amplitude Ey of2And phase phi2；
Step 200, Ey to be calculated2Measuring the electric field amplitude at X2 position with infrared camera
Figure FDA0002292735850000014
Instead, the phase phi2The constant, magnetic field amplitude at the X2 position is shown as
Figure FDA0002292735850000015
Wherein, η0377 omega is free space wave impedance, and then the electric field amplitude Ey of the X3 position is calculated by the formula of Stratton-Chu3And phase phi3；
The Stratton-Chu formula:
Figure FDA0002292735850000016
where μ is the vacuum permeability, ε is the vacuum dielectric constant, ω is the microwave angular frequency, k0Is the wave number in the free space of the electromagnetic wave,
Figure FDA0002292735850000017
and
Figure FDA0002292735850000018
respectively expressed as the electric field intensity and the magnetic field intensity of a certain plane in space, S is the area of a certain plane,
Figure FDA0002292735850000019
is the outer normal direction of a certain plane,
Figure FDA00022927358500000110
is a free space Green function, R is the distance from a point on a certain plane to a space point,
Figure FDA00022927358500000111
respectively representThe electric field intensity and the magnetic field intensity radiated to a space point for a certain plane;
step 300, Ey to be calculated3Measuring the electric field amplitude at X3 position with infrared camera
Figure FDA00022927358500000112
Instead, the phase is unchanged, the magnetic field amplitude at the X3 position is
Figure FDA00022927358500000113
Wherein, η0377 omega is free space wave impedance, and then the electric field amplitude Ey of the X1 position is obtained through the inverse transmission calculation of the X3 position4And phase phi4；
Step 400, Ey to be calculated4Measuring the electric field amplitude at X1 position with infrared camera
Figure FDA0002292735850000021
Instead, the phase phi4Constant, magnetic field amplitude
Figure FDA0002292735850000022
Wherein, η0377 omega is free space wave impedance, namely one iteration is completed;
step 500, when the calculated phase phi1，Φ2，Φ3When the change is not changed, the phases of the beams at the final X1, X2 and X3 positions are obtained.
2. The design method of high power microwave phase inversion according to claim 1, characterized in that: thermal paper is arranged at the spatial positions X1, X2 and X3, burning spots are formed on the thermal paper by high-power microwaves, and energy equipotential diagrams at the spatial positions X1, X2 and X3 can be measured by an infrared camera.
3. The design method of high power microwave phase inversion according to claim 1, characterized in that: the Stretton-Chu formula is divided into electric field strength and magnetic field strength in a rectangular coordinate systemExpressed as:
Figure FDA0002292735850000023
wherein Ex, Ey and Ez are respectively x, y and z components of an electric field, and Hx and Hy Hz are respectively x, y and z components of a magnetic field;
because Ex is 0, Ey is in the infrared camera measurement, Ez is 0, Hx is 0, Ey is 0, Ez is Ey/377, and the Ey amplitude and the phase phi on the next plane are calculated according to the Stratton-Chu formula.
CN201711282811.5A 2017-12-07 2017-12-07 Design method for high-power microwave phase inversion Active CN108021757B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711282811.5A CN108021757B (en) 2017-12-07 2017-12-07 Design method for high-power microwave phase inversion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711282811.5A CN108021757B (en) 2017-12-07 2017-12-07 Design method for high-power microwave phase inversion

Publications (2)

Publication Number Publication Date
CN108021757A CN108021757A (en) 2018-05-11
CN108021757B true CN108021757B (en) 2020-03-10

Family

ID=62078914

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711282811.5A Active CN108021757B (en) 2017-12-07 2017-12-07 Design method for high-power microwave phase inversion

Country Status (1)

Country Link
CN (1) CN108021757B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109726448B (en) * 2018-12-11 2023-03-24 中国人民解放军63660部队 Sectional mode converter design method based on iteration method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412983A (en) * 2013-07-22 2013-11-27 电子科技大学 Quasi-optical phase correction surface designing method
CN105093203A (en) * 2015-08-31 2015-11-25 中国科学院电子学研究所 Star subsurface remote sensing detection radar echo simulation and parameter inversion method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412983A (en) * 2013-07-22 2013-11-27 电子科技大学 Quasi-optical phase correction surface designing method
CN105093203A (en) * 2015-08-31 2015-11-25 中国科学院电子学研究所 Star subsurface remote sensing detection radar echo simulation and parameter inversion method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
3D wavebeam field reconstruction from intensity measurements in a few cross sections;A.V. Chirkov 等;《Optics Communications》;19950430;第449-452页 *
Design of Adapted Phase Correcting Mirrors for Gyrotrons;Jianwei Liu1 等;《2013 IEEE 14th International Vacuum Electronic,International Conference》;20131231;第1-2页 *

Also Published As

Publication number Publication date
CN108021757A (en) 2018-05-11

Similar Documents

Publication Publication Date Title
Colas et al. RF current distribution and topology of RF sheath potentials in front of ICRF antennae
Rozhnev et al. Studying characteristics of the slow-wave system of the traveling-wave tube with a sheet electron beam
Liu et al. Vector method for synthesis of adapted phase-correcting mirrors for gyrotron output couplers
Maxworth et al. Multistation observations of the azimuth, polarization, and frequency dependence of ELF/VLF waves generated by electrojet modulation
CN109472066A (en) Reflector antenna random error analysis method based on the displacement of unit center point
Koehn et al. Schemes of microwave heating of overdense plasmas in the stellarator TJ-K
CN108021757B (en) Design method for high-power microwave phase inversion
Kumar et al. Magnetically enhanced THz generation by self-focusing laser in VA-MCNTs
Kamol et al. Three-dimensional standing waves in a microwave oven
Wang et al. Quasi-optical mode converter for a 0.42-THz TE 26 mode pulsed gyrotron oscillator
Shafir et al. Characterization of inductively coupled plasma generated by a quadruple antenna
Kowalski et al. Simple correctors for elimination of high-order modes in corrugated waveguide transmission lines
Hammond et al. Overdense microwave plasma heating in the CNT stellarator
Jianwei et al. Design of quasi-optical mode converter for 94GHz gyrotron
Kochetov et al. Axially symmetric transient electromagnetic fields in a radially inhomogeneous biconical transmission line
Botton et al. Compact 3-D envelope ADI-FDTD algorithm for simulations of coherent radiation sources
Hochman et al. A numerical methodology for efficient evaluation of 2D Sommerfeld integrals in the dielectric half-space problem
Gopal et al. Nonlinear mixing-based terahertz emission in inclined rippled density plasmas
Xu et al. Improvement of the Q-factor measurement in RF cavities
Wait Currents excited on a conducting surface of large radius of curvature
Gayen et al. Resonance and parametric analysis of planar broad-wall longitudinal slot array antennas
Krivets et al. The Smith-Pursell effect amplification of the electromagnetic waves in an open waveguide with a metal-dielectric layer
Sato et al. Near field leaky-wave focusing antenna using inhomogeneous rectangular waveguide
Horvath et al. Investigations of the electron temperature profiles at the WEGA stellarator
Li et al. Design of a Vlasov mode converter of 263 GHz gyrotron oscillator for DNP‐NMR

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant