A design method and system for a relativistic klystron based on multi-objective optimization

Abstract

The present invention provides a design method and system for a relativistic klystron based on multi-objective optimization. The method includes: obtaining the simulation design parameters and optimization objectives of the relativistic klystron input by the user; according to the simulation design parameters, using the self-developed klystron large-signal simulation software KlyH based on the one-dimensional disk model to simulate and analyze the relativistic klystron, and using the hybrid multi-objective optimization algorithm to optimize the optimization objectives. Through the integration of the hybrid multi-objective optimization algorithm and the klystron large-signal simulation software KlyH based on the one-dimensional disk model, it is possible to efficiently explore different design schemes of relativistic klystrons in the multi-dimensional parameter space, not only capable of simultaneously optimizing multiple key indicators such as efficiency, interaction length, and power, but also significantly improving the global search ability, quickly generating high-quality designs that meet complex requirements, and greatly improving the reliability and engineering feasibility of the relativistic klystron design in the present invention.

Description

Technical Field

[0001] The present invention relates to the technical field of relativistic klystron design, and specifically relates to a design method and system for a relativistic klystron based on multi-objective optimization. Background Art

[0002] The relativistic klystron amplifier (RKA) has the advantages of high peak power, high efficiency, high gain, and stable output microwave frequency and phase, and is of great significance for promoting the practical application of high-power microwave systems. The RKA is an electro-vacuum O-type microwave device that uses beam-wave interaction to achieve radio frequency (RF) amplification, and has broad application prospects in technical fields such as particle accelerators, pulsed radars, and power synthesis. During the process of the RKA moving towards practical application, the number of RKAs required is often large. If the efficiency of a single RKA can be improved, the total power consumption will decrease exponentially, and the value is quite significant. In the simulation and optimization process of the klystron, one-dimensional simulation has many remarkable advantages. Usually, one-dimensional simulation has a high calculation efficiency, especially in the large-scale parameter scanning and preliminary design stages. And one-dimensional simulation has the advantage of simplifying the model. Users can quickly define parameters when setting up the simulation without having to deal with complex geometric structures. It is very meaningful to conduct in-depth research and improvement on one-dimensional large-signal calculation programs.

[0003] With the significant increase in the demand for bandwidth performance, multi-objective optimization, and design flexibility in modern klystron applications, the limitations of traditional simulation methods have gradually emerged. For example, some one-dimensional simulation software uses fixed parameters and single-objective optimization methods, resulting in an overly long design cycle and difficulty in quickly approaching the target in terms of efficiency and accuracy.

[0004] Therefore, the existing technology still needs to be further developed. Summary of the Invention

[0005] The purpose of the present invention is to overcome the above technical deficiencies, and provide a design method and system for a relativistic klystron based on multi-objective optimization to solve the problems existing in the existing technology.

[0006] To achieve the above technical purpose, according to the first aspect of the present invention, the present invention provides a design method for a relativistic klystron based on multi-objective optimization, including:

[0007] S100. Obtain the simulation design parameters and optimization objectives of the relativistic klystron input by the user;

[0008] S200. According to the simulation design parameters, use the klystron large-signal simulation software KlyH based on a one-dimensional disk model to simulate and simulate the relativistic klystron, and use a hybrid multi-objective optimization algorithm to optimize the optimization objectives to obtain an optimized relativistic klystron.

[0009] Specifically, the simulation design parameters include:

[0010] The microwave parameters, calculation parameters, electron beam parameters, and resonator modulation parameters of the relativistic klystron.

[0011] Specifically, the hybrid multi-objective optimization algorithm includes a first-stage rough search and a second-stage fine optimization.

[0012] Specifically, the method for optimizing the optimization objective by using the hybrid multi-objective optimization algorithm includes:

[0013] S210. Determine the optimization parameters;

[0014] S220. Conduct the first-stage rough search, and use the non-dominated sorting genetic algorithm to generate multiple design schemes for the optimization parameters of the relativistic klystron;

[0015] S230. Conduct the second-stage fine optimization, finely optimize the multiple design schemes for the optimization parameters of the relativistic klystron generated in the first-stage rough search, then generate multiple design schemes for the optimization parameters of the relativistic klystron, and select the optimal scheme from the multiple design schemes.

[0016] Specifically, S220 includes:

[0017] Use the non-dominated sorting genetic algorithm to generate multiple design schemes for the optimization parameters of the relativistic klystron, and screen the Pareto optimal solution set of the non-dominated sorting genetic algorithm according to the physical characteristics of the relativistic klystron;

[0018] The optimization parameters meet the requirements of the preset physical constraint mechanism.

[0019] Specifically, the method for conducting the first-stage rough search and using the non-dominated sorting genetic algorithm to generate multiple design schemes for the optimization parameters of the relativistic klystron specifically includes:

[0020] Use the non-dominated sorting genetic algorithm to generate multiple design schemes for the optimization parameters of the relativistic klystron, use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to calculate the objective function values of the optimization objectives obtained from the multiple design schemes for the optimization parameters of the relativistic klystron in sequence, take the multiple design schemes for the optimization parameters of the relativistic klystron as the initial population, then conduct non-dominated sorting on the initial population, divide it into different ranks, and calculate the crowding distance to measure the diversity between individuals, and then enter the optimization iteration loop. The specific process of the optimization iteration loop is as follows:

[0021] Generate an offspring population through selection, crossover, and mutation operations. After merging the parent population and the offspring population, perform non-dominated sorting again, and select the first preset number of individuals according to the rank and crowding degree to form a new generation population. Determine whether the number of iterations reaches the preset number of iterations or converges. If so, output the Pareto solution set; otherwise, repeat the optimization iteration loop process again until the preset number of iterations or convergence is reached.

[0022] Specifically, the S230 includes:

[0023] Perform the second-stage fine optimization. Use the design schemes of the optimization parameters of multiple relativistic klystrons generated by the first-stage rough search as part of the initial particles of the multi-objective particle swarm optimization algorithm for the optimal solution, then generate design schemes of the optimization parameters of multiple relativistic klystrons, and select the optimal scheme from multiple design schemes.

[0024] Specifically, the preset physical constraint mechanism includes:

[0025] Apply a physical constraint range to the optimization parameters.

[0026] Specifically, the optimization objectives include at least two of the following parameters:

[0027] The parameters include the length of the electron beam drift tube section, output efficiency, and bandwidth of the relativistic klystron.

[0028] According to the second aspect of the present invention, there is provided a relativistic klystron design system based on multi-objective optimization, including:

[0029] An acquisition module: used to acquire the simulation design parameters and optimization objectives of the relativistic klystron input by the user;

[0030] An optimization module: used to simulate and simulate the relativistic klystron using the klystron large-signal simulation software KlyH based on the one-dimensional disk model according to the simulation design parameters, and optimize the optimization objectives using a hybrid multi-objective optimization algorithm to obtain an optimized relativistic klystron.

[0031] Beneficial effects:

[0032] The present invention provides a design method for a relativistic klystron based on multi-objective optimization. The large-signal simulation software KlyH of the klystron based on a one-dimensional disk model is used to simulate and imitate the relativistic klystron, and a hybrid multi-objective optimization algorithm is used to optimize the optimization objectives. Through the integration of the hybrid multi-objective optimization algorithm and the large-signal simulation software KlyH of the klystron based on a one-dimensional disk model, it is possible to efficiently explore different design schemes of relativistic klystrons in a multi-dimensional parameter space. It can not only optimize multiple key indicators such as efficiency, interaction length, and power simultaneously, but also significantly improve the global search ability, quickly generate high-quality designs that meet complex requirements. The multi-objective optimization method of this application is particularly suitable for the analysis of a large number of coupled variables and nonlinear effects involved in the design of high-efficiency and compact RKA, thereby breaking through the bottlenecks in terms of efficiency, bandwidth, and engineering practicability of traditional methods, and greatly improving the reliability and engineering feasibility of the design of the relativistic klystron in the present invention. Description of the Drawings

[0033] Figure 1 is a flowchart of the design method for a relativistic klystron based on multi-objective optimization provided in the specific embodiment of the present invention;

[0034] Figure 2 is a schematic diagram of the composition of the design system for a relativistic klystron based on multi-objective optimization provided in the specific embodiment of the present invention;

[0035] Figure 3 is a flowchart of the operation of the large-signal simulation software KlyH of the klystron based on a one-dimensional disk model provided in the specific embodiment of the present invention;

[0036] Figure 4 is a flowchart of the implementation of the hybrid multi-objective optimization algorithm provided in the specific embodiment of the present invention;

[0037] Figure 5 is a diagram of the solution set after hybrid optimization provided in the specific embodiment of the present invention;

[0038] Figure 6 are the results of optimizing the efficiency and interaction length respectively using the NSGA-II algorithm and the OMOPSO algorithm provided in the specific embodiment of the present invention;

[0039] Figure 7 Schematic diagram of the comparison of the calculation efficiencies of the NSGA-II algorithm and the OMOPSO algorithm provided in the specific embodiment of the present invention. Detailed Embodiment

[0040] To enable those skilled in the art to better understand the technical solution of the present invention, the technical solution of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the present invention. Based on the embodiments in this application, other similar embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the scope of protection of this application. In addition, the directional terms mentioned in the following embodiments, such as "up", "down", "left", "right", etc., are only references to the directions of the accompanying drawings. Therefore, the directional terms used are for illustration rather than to limit the present invention.

[0041] The present invention will be further described below in conjunction with the accompanying drawings and preferred embodiments.

[0042] Please refer to Figure 1 , this embodiment provides a design method for a relativistic klystron based on multi-objective optimization, including:

[0043] S100. Obtain the simulation design parameters and optimization objectives of the relativistic klystron input by the user;

[0044] S200. According to the simulation design parameters, use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to simulate the relativistic klystron, and use the hybrid multi-objective optimization algorithm to optimize the optimization objectives to obtain an optimized relativistic klystron.

[0045] It should be noted here that the klystron large-signal simulation software KlyH based on the one-dimensional disk model is an independently controllable relativistic klystron simulation software designed by the present invention, which is used to realize the one-dimensional large-signal simulation of the relativistic klystron, can efficiently simulate the working process of the relativistic klystron, and can efficiently perform calculation tasks such as large-signal simulation, bandwidth analysis and efficiency evaluation. Through highly optimized numerical calculation methods, the KlyH software greatly improves the simulation speed, supports fast scanning of large-scale parameter spaces and multi-objective optimization calculations, can accurately simulate the behavior of the relativistic klystron under different working conditions, and generates accurate simulation results, ensuring the scientificity and accuracy of design decisions, and greatly meeting the requirements of high-efficiency and high-precision simulations.

[0046] Furthermore, the klystron large-signal simulation software KlyH based on the one-dimensional disk model has a numerical program developed in Fortran language based on the one-dimensional disk model as its core. It can simulate the modulation process of the electron beam in the cavity and calculate key performance parameters such as electron beam efficiency, interaction length, and gain. Compared with traditional analytical methods, this numerical simulation program can more accurately describe the electron-microwave interaction, especially suitable for the nonlinear analysis of relativistic klystrons. This method can perform simulation calculations based on input parameters such as electron beam parameters, cavity geometry parameters, input microwave power, and resonator modulation parameters, and provide high-precision performance output values, providing a reliable physical basis for optimization.

[0047] During the optimization process, the KlyH software serves as the core for calculating the objective function, providing simulation feedback for the optimization parameter design schemes of different relativistic klystrons. Compared with traditional empirical design methods, this method quantifies the relationships between parameters such as efficiency, gain, and interaction length through numerical calculations, avoiding the high-cost problem of over-reliance on three-dimensional simulations and experiments.

[0048] The core algorithms of KlyH one-dimensional simulation include the following parts:

[0049] (1) One-dimensional electron disk model: For a solid electron beam, the following assumptions are made: It is considered to be composed of many layers of electrons, each layer is an infinitely thin electron disk, and the behavior of each electron on the disk is the same, without transverse motion, and the disk does not deform during the motion, just like a rigid body. It is assumed that the disk is transparent, that is, other disks are allowed to pass through. When two disks are infinitely close, the force between them is zero, the focusing magnetic field is infinitely large, and the electron beam operates in the space-charge-limited current state;

[0050] (2) Position and velocity model of electron disks: Calculated according to kinematic theory through relativistic equations of motion;

[0051] (3) Gap impedance model of resonators: Calculate the gap impedances of the input cavity, output cavity, and intermediate cavities according to the equivalent circuit;

[0052] (4) Induced current and gap voltage model: By performing Fourier expansion on the induced current, the induced current components of each harmonic can be obtained for calculating the induced current, and then the gap voltage can be calculated based on parameters such as gap impedance and induced current;

[0053] (5) Output power and efficiency calculation model: Calculate parameters such as the output power and efficiency of the relativistic klystron based on data such as the high-frequency voltage, current, and external quality factor of the output cavity.

[0054] Specifically, the simulation design parameters include:

[0055] The microwave parameters, calculation parameters, electron beam parameters, and resonator cavity modulation parameters of a relativistic klystron.

[0056] It should be further noted that in the multi-objective optimization design method of the relativistic klystron in this embodiment, the microwave parameters of the relativistic klystron include the input microwave power and the input microwave frequency; the calculation parameters include the number of disks, the number of times of pushing the disks in each high-frequency period, the maximum number of iterations, and the initial value of the total number of times of pushing the disks; the electron beam parameters include the beam voltage, beam current, beam radius, and drift tube radius; the resonator cavity modulation parameters include the resonant frequency of the resonator cavity, the axial position of the resonator cavity, the characteristic impedance of the resonator cavity, the gap width, and the harmonic number. When using the klystron large-signal simulation software KlyH based on the one-dimensional disk model for simulation, the gap voltage and induced current of the resonator cavity are calculated one by one through a loop, and it is judged whether the gap voltage converges. The output cavity is the resonator cavity that finally extracts energy in the relativistic klystron, and its convergence state directly determines the performance index. The convergence result can ensure that the simulation result reflects the steady-state physical process and avoids transient errors. If the gap voltage converges and it is the output cavity, then the output efficiency and output data are calculated.

[0057] In some specific embodiments, the optimization objectives include at least two of the following parameters:

[0058] The parameters include the electron beam drift tube section length, output efficiency, bandwidth, etc. of the relativistic klystron.

[0059] It can be understood that in the simulation of the relativistic klystron, improving the efficiency usually requires a longer interaction length, that is, the electron beam drift tube section length, while shortening the interaction length often reduces the output efficiency. Therefore, the optimized design of the relativistic klystron needs to balance between high efficiency and a compact structure. Therefore, the performance optimization objectives are high efficiency and low tube length. The optimization of the present invention is carried out based on specific relativistic klystron examples to ensure that the optimization results meet physical constraints and engineering feasibility.

[0060] Furthermore, the process of converting the simulation design parameters input by the user into a klystron simulation example is as follows:

[0061] (1) Read the microwave parameters, calculation parameters, electron beam parameters, and resonator cavity modulation parameters input by the user;

[0062] (2) Construct the initialization model data of the relativistic klystron according to these parameters;

[0063] (3) Set the initial parameters of the electron beam, including the beam voltage, beam current, beam radius, drift tube radius, etc.;

[0064] (4) Set the initial parameters of the input microwave, including the input microwave power 、 the input microwave frequency, etc.;

[0065] (5) Set calculation parameters, including the number of disks, the number of times to push the disks in each high-frequency cycle, the maximum number of iterations, the initial value of the total number of times to push the disks, etc.;

[0066] (6) Set the modulation parameters of each resonant cavity, including the characteristic impedance R / Q of the resonant cavity, the resonant frequency, the axial position of the resonant cavity, the external quality factor of the resonant cavity, etc.;

[0067] (7) Transfer the constructed model to the core algorithm program of KlyH one-dimensional simulation for calculation.

[0068] After the core algorithm program of KlyH one-dimensional simulation is calculated, it will return the simulation results of the relativistic klystron, including parameters such as output power, efficiency, gain, bandwidth, etc.

[0069] Specifically, the hybrid multi-objective optimization algorithm includes a first-stage rough search and a second-stage fine optimization. The method for optimizing the optimization objective using the hybrid multi-objective optimization algorithm includes:

[0070] S210. Determine the optimization parameters, where the optimization parameters include but are not limited to the resonant frequency of the resonant cavity, the drift tube length between cavities, the external quality factor of the input and output cavities, etc.;

[0071] S220. Conduct the first-stage rough search, and use the non-dominated sorting genetic algorithm to quickly generate multiple design schemes of the optimization parameters of the relativistic klystron;

[0072] Specifically, the S220 includes:

[0073] Use the non-dominated sorting genetic algorithm to generate multiple design schemes of the optimization parameters of the relativistic klystron, and screen the Pareto optimal solution set of the non-dominated sorting genetic algorithm according to the physical characteristics of the relativistic klystron;

[0074] Specifically, the method for conducting the first-stage rough search and using the non-dominated sorting genetic algorithm to generate multiple design schemes of the optimization parameters of the relativistic klystron specifically includes:

[0075] Use the non-dominated sorting genetic algorithm to generate multiple design schemes of the optimization parameters of the relativistic klystron. Use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to calculate the objective function values of the optimization objectives obtained from multiple design schemes of the optimization parameters of the relativistic klystron in turn. Take multiple design schemes of the optimization parameters of the relativistic klystron as the initial population, and then conduct non-dominated sorting on the initial population, divide it into different levels, and calculate the crowding distance to measure the diversity between individuals, and then enter the optimization iteration loop. The specific process of the optimization iteration loop is as follows:

[0076] Generate the offspring population through selection, crossover, and mutation operations. After merging the parent population and the offspring population, perform non-dominated sorting again, and select the first preset number of individuals according to rank and crowding degree to form a new generation population. Determine whether the number of iterations reaches the preset number of iterations or converges. If so, output the Pareto solution set; otherwise, repeat the optimization iteration loop process again until the preset number of iterations or convergence is reached.

[0077] Furthermore, the first preset number in this embodiment can be selected by the user according to actual needs. For example, the first N individuals can be selected to form a new generation population. The preset number of iterations can also be set according to the user's settings. For example, the preset number of iterations can be set to 100 times.

[0078] The optimization parameters meet the requirements of the preset physical constraint mechanism.

[0079] Specifically, the preset physical constraint mechanism includes: imposing a physical constraint range on the optimization parameters.

[0080] It should be further noted that for the physical constraints of the relativistic klystron in this embodiment, the optimization examples and the selection of optimization parameters are both through physical simulation and effectiveness analysis. Since the design of the relativistic klystron is restricted by many factors, variables cannot be adjusted arbitrarily during the optimization process, but must be optimized within the range that can be realized in engineering. To ensure the feasibility of the optimization results, a boundary constraint mechanism, that is, a preset physical constraint mechanism, is introduced in this method to ensure that the solutions of the optimization parameters do not violate physical laws. Since the physical characteristics such as the electron beam characteristics, resonator parameters, and input microwave power of the relativistic klystron are coupled with each other, the range of the parameters to be optimized cannot be given arbitrarily, but the key design variables need to be adjusted within a reasonable range.

[0081] See Figure 4 , in some specific embodiments, the specific steps of the first-stage coarse search are as follows:

[0082] Step 1: Input stage;

[0083] Determine the basic parameters (simulation design parameters) of the relativistic klystron instance to be optimized, determine the basic parameters of the non-dominated sorting genetic algorithm, and determine that the range of the optimization parameters is within the physically allowed range;

[0084] Step 2: Initialize the population;

[0085] Generate an initial population using the NSGA-II algorithm, and these individuals represent different schemes for the design of the relativistic klystron;

[0086] Step 3: Calculate the objective function value;

[0087] Evaluate each individual, and use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to calculate the objective function values of key indicators such as its output efficiency and interaction length.

[0088] Step 4: Fast non-dominated sorting and crowding degree calculation;

[0089] a. Divide the population into different ranks through non-dominated sorting to determine which individuals are the current optimal solutions;

[0090] b. Calculate the crowding degree of each individual to maintain the diversity of the population.

[0091] Step 5: Selection, crossover, and mutation;

[0092] a. Adopt the elite selection strategy to retain excellent individuals.

[0093] b. Perform crossover and mutation operations to generate a new offspring population.

[0094] Step 6: Output the Pareto solution set;

[0095] After multiple iterations, a preliminary Pareto optimal solution set is obtained. These solutions cover the design space of the relativistic klystron and provide multiple output efficiency - tube length trade-off schemes.

[0096] It can be understood that in this embodiment, the klystron large-signal simulation software KlyH based on the one-dimensional disk model is used to calculate the objective function, avoiding physically infeasible solutions that may occur in pure mathematical optimization. By combining the actual working principle of the relativistic klystron, the optimization results can meet the effectiveness of the electron beam - microwave interaction and the rationality of the cavity resonance frequency and cavity position. This method not only improves the applicability of the optimization algorithm but also enhances the reliability and engineering feasibility of the relativistic klystron design.

[0097] S230. Perform the second-stage fine optimization. Fine-optimize the design schemes of the optimization parameters of multiple relativistic klystrons generated in the first-stage rough search, then generate multiple design schemes of the optimization parameters of the relativistic klystron, and select the optimal scheme from the multiple design schemes.

[0098] Specifically, the S230 includes:

[0099] Perform the second-stage fine optimization. Use the design schemes of the optimization parameters of multiple relativistic klystrons generated in the first-stage rough search as part of the initial particles of the multi-objective particle swarm optimization algorithm for the optimal solution, then generate multiple design schemes of the optimization parameters of the relativistic klystron, and select the optimal scheme from the multiple design schemes.

[0100] It can be understood that the above design schemes of the optimization parameters all meet the requirements of the preset physical constraint mechanism.

[0101] Specifically, the method specifically includes:

[0102] Use the optimal solution multi-objective particle swarm algorithm to generate design schemes for the optimization parameters of multiple relativistic klystrons, initialize their positions and velocities, and use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to calculate the objective function values of the optimization objectives obtained from the design schemes of the optimization parameters of multiple relativistic klystrons in sequence. At the same time, the individual optimal solution of each particle is initialized to its current position, and the global optimal solution is initialized to the optimal particle in the population. The velocity of each particle is adjusted according to the individual and global optimal solutions, and the position is updated accordingly. If the new position of the particle is better than its individual optimal solution, the individual optimal solution is updated. Then, the particles are mutated with a certain probability to increase diversity and avoid falling into local optima. The updated particles need to recalculate the objective function values. Then, through non-dominated sorting, all non-dominated particles in the population are identified and added to the non-dominated set. If the new position of the particle is better than its current individual optimal solution, the individual optimal solution is updated, and the non-dominated particles are stored in the external set to maintain the diversity of the solution set. Then, according to the particles in the external set, the global optimal solution is updated until the preset number of iterations or convergence is reached, and the Pareto solution set is output.

[0103] It can be understood that for the optimization problem of high-efficiency compact klystrons, those skilled in the art can find that the non-dominated sorting genetic algorithm (NSGA-II) has a low time cost and the optimal solution set is relatively dispersed, which is suitable for global fast search. At the same time, the optimal solution multi-objective particle swarm algorithm (OMOPSO) has a high time cost and the optimal solution set is relatively concentrated, which is suitable for local optimal search. Therefore, a hybrid multi-objective optimization algorithm (including NSGA-II and OMOPSO) is used to optimize the optimization objectives, as Figure 6 and Figure 7 shown Figure 6 are the results of optimizing the output efficiency and interaction length using the NSGA-II algorithm and the OMOPSO algorithm respectively. It can also be seen from this that the optimal solution set obtained by the NSGA-II algorithm is relatively dispersed, but the general efficiency is lower than that of the OMOPSO algorithm. The optimal solution set obtained by the OMOPSO algorithm is relatively concentrated, but the total interaction length obtained is shorter; Figure 7The horizontal axis (population size) represents the size of the population during the algorithm's operation, ranging from 0 to 1000 individuals. The vertical axis (computational efficiency) represents the computational time cost of the algorithm under different population sizes, with the unit being minutes. It can be seen from the figure that in the case of 5 generations and 10 generations of the OMOPSO algorithm, as the population size increases, the computational time increases significantly. For example, when there are 1000 individuals, the computational time of OMOPSO_10 generations is close to 200 minutes. Therefore, the time cost of the OMOPSO algorithm is relatively high, especially when the population size is large. While for the NSGA-II algorithm in the case of 5 generations and 10 generations, the computational time is significantly lower than that of the OMOPSO algorithm. Even at the maximum population size (1000 individuals), the computational time of NSGA-II_10 generations is only about 30 minutes. Therefore, the time cost of the NSGA-II algorithm is low and does not change much with the population size.

[0104] Furthermore, according to Figure 6 and Figure 7 , the multi-objective optimization strategy that combines the NSGA-II algorithm and the OMOPSO algorithm can give full play to their respective advantages and achieve efficient global and local optimization. Figure 5 is the optimization result obtained by optimizing the total interaction length and output efficiency using a hybrid multi-objective optimization algorithm that combines the NSGA-II algorithm and the OMOPSO algorithm. Its optimization effect is significantly better than Figure 6 the effects when only the NSGA-II algorithm or only the OMOPSO algorithm is used alone in . The output efficiency has increased by nearly 10%, and at the same time, the interaction length has been shortened. Therefore, the combination of the above two algorithms can give full play to their respective advantages and enable the present application to achieve a balance in the optimization between high efficiency and a compact structure.

[0105] It should be further noted that the multi-objective particle swarm optimization algorithm (MOPSO) was proposed by Carlos Artemio Coello Coello in 2002. This algorithm is an extension based on the particle swarm optimization algorithm and aims to solve multi-objective optimization problems. The optimal multi-objective particle swarm algorithm (OMOPSO algorithm) in the present invention is proposed based on MOPSO. The OMOPSO algorithm constructs a non-dominated solution set based on the Pareto dominance relationship, uses an external set to save the currently found non-dominated solution set, proposes the concept of strong ε-dominance based on ε-dominance, and uses it to update the external set so that the algorithm can maintain good distribution. By improving the selection methods of the global best and individual best in the MOPSO algorithm, at the same time adopting a new particle update strategy to accelerate the convergence of the solution set, adding an adaptive mutation operator to avoid falling into a local non-dominated optimal solution, and proposing a non-dominated set construction method based on quicksort to accelerate the running efficiency of the algorithm. The main framework of the OMOPSO algorithm is as follows:

[0106] (1) Initialization

[0107] Randomly generate a population containing N particles, where each particle represents a solution. Initialize its position and velocity, and calculate the objective function value of each particle. At the same time, initialize the individual best solution of each particle to its current position, and initialize the global best solution to the best particle in the population;

[0108] (2) Iteration process (iteration reaches the preset maximum number of generations or the solution converges)

[0109] The velocity of each particle is adjusted according to the individual and global best solutions, and the position is updated accordingly. If the new position of the particle is better than its individual best solution, update the individual best solution, and then mutate the particle with a certain probability to increase diversity and avoid falling into a local optimum. The updated particle needs to recalculate the objective function value; then, through non-dominated sorting, identify all non-dominated particles in the population and add them to the non-dominated set. If the new position of the particle is better than its current individual best solution, update the individual best solution; store the non-dominated particles in the external set to maintain the diversity of the solution set, and then update the global best solution according to the particles in the external set.

[0110] See Figure 4 , in some specific embodiments, the specific steps of the second-stage fine optimization are as follows:

[0111] Step 1: Pareto solution screening and conversion;

[0112] Convert the Pareto optimal solution set generated from the first-stage rough search into particles in the particle swarm optimization algorithm.

[0113] Step 2: Initialize the particle swarm:

[0114] Use the screened Pareto solutions as part of the initial particles to initialize the particle swarm.

[0115] Step 3: Particle update and local search;

[0116] a. The OMOPSO algorithm uses the information exchange and learning mechanism between particles to update the position and velocity of particles;

[0117] b. Conduct local refinement search in the search space to further improve the quality and accuracy of the solution.

[0118] Step 4: Output Pareto solutions;

[0119] After local optimization, obtain the final Pareto optimal solution set, and these solutions have higher convergence and accuracy.

[0120] Step 5: Output stage

[0121] Import the final Pareto solution set into the KlyH simulation software for verification and visualization analysis to ensure the feasibility and performance of the design scheme.

[0122] It can be understood that in the solution proposed in this embodiment, which combines the KlyH and multi-objective optimization hybrid algorithm, since the NSGA-II algorithm performs excellently in terms of computational output efficiency and global search, while the OMOPSO algorithm performs excellently in local search, a multi-objective optimization method combining the NSGA-II algorithm and the OMOPSO algorithm is adopted. First, the NSGA-II algorithm is used in the first-stage coarse search to construct a preliminary Pareto solution set within the range of optimization parameters of the given klystron instance. The goal of the first stage is to cover the design space of the relativistic klystron as much as possible, ensure the diversity of the optimization scheme, and provide multiple efficiency-length trade-off schemes for subsequent optimization. Then, based on the Pareto optimal solution set generated by the NSGA-II algorithm, the OMOPSO algorithm is used in the second-stage fine optimization to further improve the quality of the solution. As an improved version of the particle swarm optimization algorithm, the OMOPSO algorithm can use the Pareto solutions of the NSGA-II algorithm as guiding particles to further refine the design scheme in the search space, improve the convergence and accuracy of the optimization solution, and at the same time strictly meet the physical constraints of the klystron. This method can not only achieve a balance between high efficiency and a compact structure, but also ensure that the optimization results meet the actual engineering applications. Compared with a single optimization method, it can find a relativistic klystron design scheme that takes into account both high efficiency and a compact structure more efficiently, further improving the feasibility and stability of the relativistic klystron design of the present invention.

[0123] See Figures 3 - 5 , the working principle of the present invention will be described below through specific examples:

[0124] 1. According to the idea of the present invention, design and implement the principle of a relativistic klystron simulation framework based on KlyH numerical simulation. Figure 3 is the working principle diagram of the klystron large-signal simulation software KlyH based on a one-dimensional disk model, and the specific process is as follows;

[0125] First, input simulation design parameters such as DC electron beam parameters, resonant cavity high-frequency parameters, total number of resonant cavities, injected microwave power and frequency, number of disks, maximum number of iterations, etc., and initialize the resonant cavity gap voltage and induced current;

[0126] Then, use the electron load conductance and resonant cavity gap impedance calculation models to provide a basis for the gap voltage calculation. During the calculation process, the gap voltage and induced current of each resonant cavity are solved one by one through a loop, and it is judged whether the gap voltage converges. If it does not converge, continue to iterate until the convergence condition is met. Specifically, it includes:

[0127] (1) Start and input data, where the input data includes: DC electron beam parameters, high-frequency parameters of the resonant cavity, total number of resonant cavities Ncav, injected microwave power and frequency, number of disks, maximum number of iterations itermax, and initial value of the total number of disk pushes;

[0128] (2) Initialize by assigning initial values to the gap voltages and induced currents of each resonant cavity, and initialize the loop variables i = 1 and j = 0, where i is used to traverse the resonant cavities and j is used to control the number of iterations;

[0129] (3) Calculate the electron load conductance and the gap impedances of each resonant cavity to prepare for the subsequent calculation of the gap voltages;

[0130] (4) Calculate the gap voltages and induced currents of the resonant cavities one by one in a loop, and determine whether the gap voltages converge. If they do not converge, push the disks to recalculate the induced currents and gap voltages until the gap voltages converge:

[0131] (5) Determine whether the current cavity is the output cavity. If so, calculate the efficiency and output data; otherwise, return to step (4);

[0132] Finally, if the calculation reaches the output cavity, solve for the electron beam output efficiency and other output parameters; otherwise, return to continue the iterative calculation.

[0133] 2. Hybrid multi-objective optimization method based on NSGA-II algorithm and OMOPSO algorithm, and the specific process is as Figure 4 shown;

[0134] In this example, taking the optimization of a 6-cavity single-injection relativistic klystron as an example, the steps to determine the instance parameters of the relativistic klystron are as follows. The simulation design parameters include calculation parameters: the number of disks is 64, each disk is pushed 32 times per high-frequency cycle, the maximum number of iterations is 200, and the initial value of the total number of disk pushes is 1000; electron beam parameters: beam voltage V0 = 450 kV, beam current I0 = 105 A, beam radius Rb = 2.1 mm, drift tube radius Rt = 3.5 mm; microwave parameters: input microwave power P0 = 100 W, input microwave frequency f0 = 10 GHz; the modulation parameters of each resonant cavity are shown in Table 1, where f is the resonant frequency of the resonant cavity, z is the relative position of the resonant cavity along the axis with respect to the input cavity, Q ext is the external quality factor of the resonant cavity. To meet the requirements of engineering design, the value of Q ext of the bunching cavity is taken as 950000, R / Q is the characteristic impedance of the resonant cavity. Using the klystron large-signal simulation software KlyH based on the one-dimensional disk model, the output efficiency of the relativistic klystron in this example is 55.2%, and the total interaction length (the length of the electron beam drift tube section) is 43.57 cm.

[0135] Table 1 Modulation Parameters of Each Resonator before Optimization

[0136]

[0137] It is determined that the optimization objectives of this example are the electron beam drift tube interval length and output efficiency of the relativistic klystron. The optimization parameters are the resonant frequency (f) of the cavity, the drift tube length (l) between the cavities, and the external quality factor (Q ext ) of the input and output cavities;

[0138] It is determined that the seed number of the NSGA-II algorithm is 1000, the maximum number of iterations is 10, the mutation rate is 0.07, and the crossover rate is 0.9; it is determined that the seed number of the OMOPSO algorithm is 1000, the maximum number of iterations is 10, and the mutation rate is 0.07;

[0139] Use the NSGA-II algorithm for multi-objective optimization. Its general process is to randomly generate multiple relativistic klystron schemes according to the input data, use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to calculate the objective function value, and then start non-dominated sorting and crowding degree calculation. Through selection, crossover, and mutation, non-dominated sorting and crowding degree calculation of the new population are carried out until the number of iterations is reached or convergence occurs, and then the Pareto solution set is output;

[0140] The design schemes in the Pareto solution set are used as part of the initial particle swarm of the OMOPSO algorithm to ensure that the initial solution set has good quality for further optimization;

[0141] After the OMOPSO algorithm optimization is completed, the finally obtained optimal design scheme is output, including parameters such as high efficiency and low tube length, and the Pareto front graph is drawn, Figure 5 which is the solution set graph after hybrid optimization, and the selected optimal design scheme. Table 2 is the table of modulation parameters of each resonator after optimization. KLY5 simulation is carried out, and its efficiency is increased by 26.6%, and the total interaction length is 41.6 cm, which is reduced by 19.7 mm compared with that before optimization. It can be seen that the multi-objective optimization method of the present invention provides an intuitive basis for engineering decision-making.

[0142] Table 2 Modulation Parameters of Each Resonator after Optimization

[0143]

[0144] 3. Optimize the optimization strategy combining variable selection and physical constraints, specifically as follows:

[0145] This example targets the physical constraints of a relativistic klystron. When selecting examples and optimizing parameters for physical analysis, key parameters such as the resonant frequency of the resonator, the length of the drift tube between cavities, and the external quality factor of the input and output cavities are finally selected as optimization parameters. Due to limitations such as the velocity modulation and bunching of the electron beam, the application of microwave resonators, and structural design and manufacturing processes, variables cannot be adjusted arbitrarily during the optimization process, but must be optimized within the range achievable in engineering. To ensure the feasibility of the optimization results, a boundary constraint mechanism is introduced in this application during the optimization process, that is, a definite floating range is set for each optimization parameter to ensure that the optimization solution does not violate physical laws. This method not only improves the applicability of the optimization algorithm, avoids physically infeasible solutions that may occur in pure mathematical optimization, ensures that the final design scheme can effectively guide 3D simulation and practical applications, and improves the reliability and engineering achievability of relativistic klystron design.

[0146] It should be noted here that this embodiment provides a design method for a relativistic klystron based on multi-objective optimization. The klystron large-signal simulation software KlyH based on a one-dimensional disk model is used to simulate and analyze the relativistic klystron, and a hybrid multi-objective optimization algorithm is used to optimize the optimization objectives. Through the integration of the hybrid multi-objective optimization algorithm and the klystron large-signal simulation software KlyH based on a one-dimensional disk model, different design schemes of relativistic klystrons can be efficiently explored in a multi-dimensional parameter space. It can not only optimize multiple key indicators such as efficiency, interaction length, and power simultaneously, but also significantly improve the global search ability, quickly generate high-quality designs that meet complex requirements. The multi-objective optimization method of this application is particularly suitable for the analysis of a large number of coupled variables and non-linear effects involved in the design of high-efficiency and compact RKA, thus breaking through the bottlenecks in terms of efficiency, bandwidth, and engineering practicality of traditional methods, and greatly improving the reliability and engineering achievability of the relativistic klystron design in the present invention.

[0147] Please refer to Figure 2 , this embodiment provides a design system for a relativistic klystron based on multi-objective optimization. The system includes:

[0148] Acquisition module 100: used to acquire the simulation design parameters and optimization objectives of the relativistic klystron input by the user;

[0149] Optimization module 200: used to simulate and analyze the relativistic klystron according to the simulation design parameters by using the klystron large-signal simulation software KlyH based on a one-dimensional disk model, and optimize the optimization objectives by using a hybrid multi-objective optimization algorithm to obtain an optimized relativistic klystron.

[0150] It should be noted that this embodiment provides a relativistic klystron design system based on multi-objective optimization. The large-signal simulation software KlyH of the klystron based on the one-dimensional disk model is used to simulate and analyze the relativistic klystron, and the hybrid multi-objective optimization algorithm is used to optimize the optimization objectives. Through the integration of the hybrid multi-objective optimization algorithm and the large-signal simulation software KlyH of the klystron based on the one-dimensional disk model, different design schemes of the relativistic klystron can be efficiently explored in the multi-dimensional parameter space. It can not only optimize multiple key indicators such as efficiency, interaction length, and power simultaneously, but also significantly improve the global search ability, quickly generate high-quality designs that meet complex requirements. The multi-objective optimization method of this application is particularly suitable for the analysis of a large number of coupled variables and nonlinear effects involved in the design of high-efficiency and compact RKA, thus breaking through the bottlenecks in terms of efficiency, bandwidth, and engineering practicability of traditional methods, and greatly improving the reliability and engineering feasibility of the relativistic klystron design in the present invention.

[0151] It should be noted that the terms "first", "second", etc. in the specification, claims and above-mentioned drawings of this application are used to distinguish similar objects, and do not have to be used to describe a specific order or sequence. It should be understood that such data can be interchanged under appropriate circumstances, so that the embodiments of the present application described here can be implemented in an order other than those illustrated or described here. In addition, the terms "comprising" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or are inherent to these processes, methods, products or devices.

[0152] The above-described technical features can be combined arbitrarily. Although not all possible combinations of these technical features are described, any combination of these technical features should be considered to be covered by this specification as long as such combination does not exist in contradiction.

[0153] The specific embodiments of the present invention described above do not constitute a limitation to the protection scope of the present invention. Any other corresponding changes and deformations made according to the technical concept of the present invention should be included in the protection scope of the claims of the present invention.

Claims

1. A design method of a relativistic klystron based on multi-objective optimization, characterized in that, Including: S100. Obtain the simulation design parameters and optimization objectives of the relativistic klystron input by the user; S200. According to the simulation design parameters, use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to simulate and simulate the relativistic klystron, and use the hybrid multi-objective optimization algorithm to optimize the optimization objectives to obtain an optimized relativistic klystron; The hybrid multi-objective optimization algorithm includes a first-stage coarse search and a second-stage fine optimization; The method of using the hybrid multi-objective optimization algorithm to optimize the optimization objectives includes: S210. Determine the optimization parameters; S220. Conduct the first-stage coarse search, use the non-dominated sorting genetic algorithm to generate multiple design schemes of the optimization parameters of the relativistic klystron, use the klystron large-signal simulation software KlyH based on the one-dimensional disk model to calculate the objective function values of the optimization objectives obtained by multiple design schemes of the optimization parameters of the relativistic klystron in turn, use multiple design schemes of the optimization parameters of the relativistic klystron as the initial population, then perform non-dominated sorting on the initial population, divide it into different ranks, and calculate the crowding distance to measure the diversity between individuals, and then enter the optimization iteration loop; The optimization parameters meet the requirements of the preset physical constraint mechanism; S230. Conduct the second-stage fine optimization, perform fine optimization on multiple design schemes of the optimization parameters of the relativistic klystron generated in the first-stage coarse search, then generate multiple design schemes of the optimization parameters of the relativistic klystron, and select the optimal scheme from multiple design schemes.

2. The design method of a relativistic klystron based on multi-objective optimization according to claim 1, characterized in that The simulation design parameters include: The microwave parameters, calculation parameters, electron beam parameters, and resonator modulation parameters of the relativistic klystron.

3. The design method of a relativistic klystron based on multi-objective optimization according to claim 1, wherein The S220 includes: Use the non-dominated sorting genetic algorithm to generate multiple design schemes of the optimization parameters of the relativistic klystron, and screen the Pareto optimal solution set of the non-dominated sorting genetic algorithm according to the physical characteristics of the relativistic klystron.

4. The design method of a relativistic klystron based on multi-objective optimization according to claim 3, characterized in that, The process of the optimization iteration loop is as follows: Generate the offspring population through selection, crossover, and mutation operations, merge the parent population and the offspring population and then perform non-dominated sorting again, and select the first preset number of individuals according to the rank and crowding degree to form a new generation population, judge whether the number of iterations reaches the preset number of iterations or converges, if so, output the Pareto solution set, otherwise, repeat the optimization iteration loop process again until the preset number of iterations or convergence is reached.

5. The design method of a relativistic klystron based on multi-objective optimization according to claim 4, characterized in that The S230 includes: Conduct the second-stage fine optimization, use multiple design schemes of the optimization parameters of the relativistic klystron generated in the first-stage coarse search as part of the initial particles of the multi-objective particle swarm optimization algorithm for the optimal solution, then generate multiple design schemes of the optimization parameters of the relativistic klystron, and select the optimal scheme from multiple design schemes.

6. The design method of a relativistic klystron based on multi-objective optimization according to claim 5, characterized in that The preset physical constraint mechanism includes: Apply a physical constraint range to the optimization parameters.

7. The design method of a relativistic klystron based on multi-objective optimization according to claim 1, characterized in that, The optimization objectives include at least two of the following parameters: The parameters include the electron beam drift tube interval length, output efficiency, and bandwidth of the relativistic klystron.

8. A relativistic klystron design system based on multi-objective optimization, characterized in that, Using the relativistic klystron design method based on multi-objective optimization according to any one of claims 1 to 7, the system includes: Acquisition module: used to acquire the simulation design parameters and optimization objectives of the relativistic klystron input by the user; Optimization module: used to simulate the relativistic klystron according to the simulation design parameters by using the klystron large-signal simulation software KlyH based on the one-dimensional disk model, and optimize the optimization objectives by using the hybrid multi-objective optimization algorithm to obtain an optimized relativistic klystron.

Images