Fusion longitudinal field coil design method based on MATLAB and COMSOL

Abstract

The invention discloses a fusion longitudinal field coil design method based on MATLAB and COMSOL, belongs to the technical field of coil design, and not only can make full use of strong data processing and optimization design capabilities of MATLAB, but also can make full use of the advantages of strong electromagnetic calculation capability and accurate electromagnetic calculation result of COMSOL to form a joint simulation model with complementary advantages. Compared with a conventional method for designing a fusion longitudinal field coil by utilizing a magnetic field formula with a simplified structure, the method has the advantage that the defect that the actual operation and theoretical design of the method have large errors is overcome by utilizing the high accuracy of MATLAB and finite element simulation calculation results. Compared with a conventional method for optimally designing the fusion longitudinal field coil by independently utilizing the MATLAB, the method has the advantage that the defects of the MATLAB in the aspect of calculation can be overcome by utilizing the advantages of the finite element in the aspects of electromagnetic and solid mechanics calculation.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of coil design, and more particularly to a fusion toroidal field coil design method based on MATLAB and COMSOL. BACKGROUND

[0002] The use of clean energy, fusion energy, is one of the dreams of mankind, but the challenges faced in its generation are very complex. Magnetic confinement fusion is considered one of the most likely ways to achieve the application of fusion energy, and the tokamak has become the mainstream research form of magnetic confinement fusion devices. The toroidal field (TF) coil, as the core component of the magnetic system of the magnetic confinement fusion device, plays a decisive role in the magnetic field strength of the fusion power and the high-parameter plasma operation. The toroidal field coil works in extreme conditions of high magnetic field and large current, and the coil itself and the associated support structure need to withstand strong electromagnetic force, which brings great challenges to the design of the coil structure.

[0003] The toroidal field coil, also known as the toroidal field coil, is an important part of the magnetic system of the fusion reactor, and its role is to generate a toroidal magnetic field for confining the plasma, and its stable operation is related to the safe operation of the entire fusion reactor. The shape design of the fusion toroidal field coil currently mainly adopts Princeton-D shape, that is, a set of 12-20 toroidal arranged toroidal field coils are simplified as a ring-shaped curved surface.

[0004] In the field of engineering practice, in order to facilitate the processing and manufacturing of the toroidal field coil, the coil skeleton for winding the cable is usually divided into multiple circular arcs and straight lines, and the size of these circular arcs and straight line segments is fitted from the actual engineering curve and the theoretical analytical curve. Due to the existence of the simplification error of the theoretical analytical curve in calculating the magnetic field and the fitting error of the actual engineering curve, this method cannot meet the design requirements of the toroidal field coil of the future compact and strong magnetic field fusion device. SUMMARY

[0005] In view of the above defects or improvement needs of the prior art, the application provides a fusion toroidal field coil design method based on MATLAB and COMSOL, which aims to solve the technical problem that the prior art cannot meet the design requirements of the toroidal field coil of the future compact and strong magnetic field fusion device.

[0006] To achieve the above-mentioned purpose, according to one aspect of the application, a fusion toroidal field coil design method based on MATLAB and COMSOL is provided, comprising: S1: initializing the geometric parameters and geometric variables of the toroidal field coil center line in MATLAB; updating and iterating the initialized geometric variables based on the initialized geometric parameters using the set genetic algorithm, obtaining the current geometric variables of the toroidal field coil center line, and passing the current geometric variables to COMSOL; S2: establishing a finite element model in COMSOL by using the received current geometric variable to calculate the tension parameter on the central line of the longitudinal field coil and transmitting the tension parameter on the central line of the longitudinal field coil to MATLAB; S3: in MATLAB, judging whether the tension parameter on the central line of the longitudinal field coil and the iteration number meet the iteration stop condition; if not, updating the current geometric variable by using the genetic algorithm and transmitting it to COMSOL, and returning to S2; if yes, outputting the current geometric variable to design the fusion longitudinal field coil.

[0007] Further, the tension parameter on the central line of the longitudinal field coil is a tension distribution on the central line of the longitudinal field coil.

[0008] Further, the tension distribution on the central line of the longitudinal field coil is expressed as: ; wherein, is the local tension of the longitudinal field coil, is the local magnetic induction intensity, is the operating current of the longitudinal field coil, is the local radius of curvature.

[0009] Further, the tension parameter on the central line of the longitudinal field coil is a tension fluctuation value on the central line of the longitudinal field coil.

[0010] Further, the tension fluctuation value on the central line of the longitudinal field coil is expressed as: ; wherein, is the tension fluctuation value on the central line of the longitudinal field coil, is the maximum tension on the center, is the minimum tension on the center.

[0011] Further, the geometric parameters include: the inner and outer radii of the central line of the longitudinal field coil. Further, the geometric variables include: the height of the straight line segment, the number of circular arc segments, the coordinates of the centers of the circular arc segments, and the number of central angles.

[0012] According to another aspect of the present application, a fusion longitudinal field coil design device based on MATLAB and COMSOL is provided, comprising: An initialization module is configured to initialize the geometric parameters and geometric variables of the central line of the longitudinal field coil in MATLAB; update the initialized geometric variables based on the initialized geometric parameters by using a set genetic algorithm to obtain the current geometric variables of the central line of the longitudinal field coil, and transmit the current geometric variables to COMSOL; a calculation module for establishing a finite element model in COMSOL by using the received current geometric variable to calculate the tension parameter on the central line of the longitudinal field coil and transmitting the tension parameter on the central line of the longitudinal field coil to MATLAB; a circulation module for judging in MATLAB whether the tension parameter on the central line of the longitudinal field coil and the iteration number satisfy an iteration stopping condition; if not, updating the current geometric variable by using a genetic algorithm and transmitting the current geometric variable to COMSOL, and returning to the calculation module to re-execute; if yes, outputting the current geometric variable to design the fusion longitudinal field coil.

[0013] According to another aspect of the present application, an electronic device is provided, comprising a memory and a processor, the memory storing a computer program, and the processor implementing the steps of the fusion longitudinal field coil design method based on MATLAB and COMSOL when executing the computer program.

[0014] According to another aspect of the present application, a computer readable storage medium is provided, storing a computer program, and the computer program implementing the steps of the fusion longitudinal field coil design method based on MATLAB and COMSOL when executed by a processor.

[0015] In general, the above technical solutions conceived by the present application can achieve the following beneficial effects compared with the prior art: The present application provides a fusion longitudinal field coil design method based on MATLAB and COMSOL. This method can make full use of the powerful data processing and optimization design capabilities of MATLAB, and make full use of the strong electromagnetic calculation capabilities and accurate electromagnetic calculation results of COMSOL, forming a complementary joint simulation model. Compared with the conventional method of designing fusion longitudinal field coils by using simplified magnetic field formulas, the high accuracy of MATLAB and finite element simulation calculation results overcomes the above-mentioned shortcomings of large errors between actual operation and theoretical design. Compared with the conventional method of optimizing and designing fusion longitudinal field coils by using MATLAB alone, the present application can overcome the shortcomings of MATLAB in calculation by taking advantage of the advantages of finite element in electromagnetic and solid mechanics calculation. BRIEF DESCRIPTION OF DRAWINGS

[0016] Figure 1 A flowchart of a fusion longitudinal field coil design method based on MATLAB and COMSOL provided for embodiment 1 of the present application; Figure 2 A flowchart of another fusion longitudinal field coil design method based on MATLAB and COMSOL provided for embodiment 1 of the present application; Figure 3A schematic diagram of the fusion field coil center line and geometric parameters and variables provided for the embodiment 1 of the present application is shown in the figure; Figure 4 A schematic diagram of the geometric model established in COMSOL provided for the embodiment 1 of the present application is shown in the figure; Figure 5 A simulation magnetic field distribution diagram for the fusion field coil center line is shown in the figure; Figure 4 Figure 6 A simulation magnetic field distribution diagram for the fusion field coil center line is shown in the figure; Figure 5 A tension distribution diagram of the fusion field coil center line after data processing is shown in the figure. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solutions and advantages of the present application clearer and more apparent, the present application is further described in detail below in combination with the figures and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and do not limit the present application. In addition, the technical features involved in each embodiment of the present application described below can be combined with each other as long as they do not conflict with each other.

[0018] Embodiment 1 The embodiment provides a fusion field coil design method based on MATLAB and COMSOL, as shown in the figures Figure 1 and Figure 2 , comprising S1-S3.

[0019] S1: initializing the geometric parameters and geometric variables of the fusion field coil center line in MATLAB; updating and iterating the initialized geometric variables based on the initialized geometric parameters by using the set genetic algorithm, obtaining the current geometric variables of the fusion field coil center line, and transmitting the current geometric variables to COMSOL. Figure 3 A schematic diagram of the fusion field coil center line and geometric parameters and variables is shown in the figure; the specific process in MATLAB is explained in detail as follows.

[0020] Firstly, the geometric parameters and variables of the fusion field coil center line are initialized in MATLAB, wherein the geometric parameters are the inner and outer radii of the fusion field coil center line, which need to be given in the initial conditions; the geometric variables are the parameters changed in the iteration process during the optimization process, mainly including the height of the straight line segment, the number of circular arc segments, the coordinates of the centers of the circular arc segments, and the number of central angles; the target optimization variable is the tension fluctuation value on the fusion field coil center line.

[0021] Then, the main parameters of the genetic algorithm are set in the MATLAB genetic algorithm function, and the main parameters of the genetic algorithm are set: the initial population size is 100, the genetic termination generation MAXGEN is 100; the selection operator adopts the tournament mode, and the selection probability is 4; the crossover operator adopts the middle crossover operator, and the crossover probability is 1; the mutation operator adopts the uniform mutation, and the mutation probability is 1.​

[0022] Finally, the MATLAB genetic algorithm generates the geometric variables of the central line of the longitudinal field coil, and passes the variable parameters to COMSOL: in MATLAB, the genetic algorithm is run, the population is initialized, the height of the straight segment, the number of circular arc segments, the coordinates of the center of each circular arc segment, and the number of degrees of the center are used as the design geometric variables, and MATLAB uses COMSOL Multiphysics with MATLAB data interface to pass the initial conditions and design geometric variables to COMSOL for modeling and simulation calculation.

[0023] S2: In COMSOL, a finite element model is established using the received current geometric variables to calculate the tension parameter on the central line of the longitudinal field coil, and the tension parameter on the central line of the longitudinal field coil is transmitted to MATLAB.

[0024] As an optional implementation, the tension parameter on the central line of the longitudinal field coil is the tension distribution on the central line of the longitudinal field coil. Wherein, the tension distribution on the central line of the longitudinal field coil is represented as: ; wherein, is the local tension of the longitudinal field coil, is the local magnetic induction, is the operating current of the longitudinal field coil, is the local radius of curvature.

[0025] As an optional implementation, the tension parameter on the central line of the longitudinal field coil is the tension fluctuation value on the central line of the longitudinal field coil. Further, the tension fluctuation value on the central line of the longitudinal field coil is represented as: ; wherein, is the tension fluctuation value on the central line of the longitudinal field coil, is the maximum tension on the center, is the minimum tension on the center.

[0026] It should be noted that, Figure 4 is a schematic diagram of the geometric model established in COMSOL by the present example. Figure 5 is a simulation magnetic field distribution diagram of Figure 4 ; Figure 6 is a simulation magnetic field distribution diagram of Figure 5 ; is a fusion longitudinal field coil central line tension distribution diagram after data processing.

[0027] S3: In MATLAB, it is judged whether the tension parameter on the central line of the longitudinal field coil and the iteration number satisfy the iteration stop condition; if not, the current geometric variables are updated using the genetic algorithm and transmitted to COMSOL, and then returned to S2; if yes, the current geometric variables are output to design the fusion longitudinal field coil.

[0028] Embodiment 2 The embodiment provides a fusion poloidal field coil design device based on MATLAB and COMSOL, and the device comprises the following modules: An initialization module is configured to initialize geometric parameters and geometric variables of a poloidal field coil center line in MATLAB, and to update and iterate the initialized geometric variables based on the initialized geometric parameters by using a set genetic algorithm, so as to obtain current geometric variables of the poloidal field coil center line and transmit the current geometric variables to COMSOL. A calculation module is configured to establish a finite element model in COMSOL by using the received current geometric variables, to calculate tension parameters on the poloidal field coil center line, and to transmit the tension parameters on the poloidal field coil center line to MATLAB. A loop module is configured to judge, in MATLAB, whether the tension parameters on the poloidal field coil center line and the number of iterations meet an iteration stop condition; if not, the current geometric variables are updated by using the genetic algorithm and transmitted to COMSOL, and the loop module returns to the calculation module to be executed again; if yes, the current geometric variables are output, so as to design a fusion poloidal field coil.

[0029] Embodiment 3 The application further relates to an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above method when executing the computer program.

[0030] The electronic device can be a desktop computer, a notebook computer, a palm computer, a cloud server and the like. The processor can be a central processing unit (CPU), and can also be another general-purpose processor, a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or another programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, and the like. The memory can be used to store a computer program and / or a module, and the processor realizes various functions of the electronic device by running or executing the computer program and / or the module stored in the memory, and calling data stored in the memory.

[0031] Embodiment 4 The application further relates to a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the steps of the above method.

[0032] In particular, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, at least one disk storage device, a flash memory device, or other volatile solid-state memory device.

[0033] Embodiment 5 The embodiment of the present application provides a computer program product or computer program, which comprises computer instructions stored in a computer readable storage medium. A processor of a computer device reads the computer instructions from the computer readable storage medium, and the processor executes the computer instructions, so that the computer device executes steps of the method of the above-mentioned embodiment of the present application.

[0034] The technical features of the above-mentioned embodiments can be combined in any manner. In order to make the description simple, all possible combinations of the technical features in the above-mentioned embodiments are not described, however, as long as the combinations of the technical features do not exist contradictory, it should be considered that the combinations are within the scope of the present application. It should be noted that the "in an embodiment of the present application", "for example", "for instance" and the like in the present application are intended to illustrate the present application, and are not used to limit the present application.

[0035] The above-mentioned embodiments only express several implementation manners of the present application, the description is more specific and detailed, but it should not be understood as a limitation to the patent application scope. It should be pointed out that, for those skilled in the art, without departing from the concept of the present application, a number of modifications and improvements can be made, which are all within the protection scope of the present application.

Claims

1. A method for designing fusion longitudinal field coils based on MATLAB and COMSOL, characterized in that, include: S1: Initialize the geometric parameters and geometric variables of the centerline of the longitudinal field coil in MATLAB; use the set genetic algorithm to update and iterate the initialized geometric variables based on the initialized geometric parameters to obtain the current geometric variables of the centerline of the longitudinal field coil, and pass the current geometric variables to COMSOL; S2: In COMSOL, a finite element model is established using the received current geometric variables to calculate the tension parameters on the center line of the longitudinal field coil, and the tension parameters on the center line of the longitudinal field coil are transmitted to MATLAB. S3: In MATLAB, determine whether the tension parameters and iteration number on the center line of the longitudinal field coil meet the iteration stopping condition; if not, use a genetic algorithm to update the current geometric variables and transmit them to COMSOL, then return to S2; if they meet the condition, output the current geometric variables to design the fusion longitudinal field coil.

2. The fusion longitudinal field coil design method based on MATLAB and COMSOL as described in claim 1, characterized in that, The tension parameter on the center line of the longitudinal field coil is the tension distribution on the center line of the longitudinal field coil.

3. The fusion longitudinal field coil design method based on MATLAB and COMSOL as described in claim 2, characterized in that, The tension distribution along the center line of the longitudinal field coil is expressed as: ； in, The local tension of the longitudinal field coil. The local magnetic flux density. This is the operating current of the longitudinal field coil. Let be the local radius of curvature.

4. The fusion longitudinal field coil design method based on MATLAB and COMSOL as described in claim 1, characterized in that, The tension parameter on the center line of the longitudinal field coil is the tension fluctuation value on the center line of the longitudinal field coil.

5. The fusion longitudinal field coil design method based on MATLAB and COMSOL as described in claim 4, characterized in that, The tension fluctuation value along the center line of the longitudinal field coil is expressed as: ； in, This represents the tension fluctuation value along the center line of the longitudinal field coil. The maximum tension at the center, The minimum tension is at the center.

6. The fusion longitudinal field coil design method based on MATLAB and COMSOL as described in claim 1, characterized in that, The geometric parameters include the inner and outer radii of the center line of the longitudinal field coil.

7. The fusion longitudinal field coil design method based on MATLAB and COMSOL as described in claim 6, characterized in that, The geometric variables include: the height of the straight line segment, the number of arc segments, the coordinates of the center of each arc segment, and the angle of the center.

8. A fusion longitudinal field coil design device based on MATLAB and COMSOL, characterized in that, include: The initialization module is used to initialize the geometric parameters and geometric variables of the centerline of the longitudinal field coil in MATLAB; it uses a set genetic algorithm to update and iterate the initialized geometric variables based on the initialized geometric parameters to obtain the current geometric variables of the centerline of the longitudinal field coil, and then passes the current geometric variables to COMSOL; The calculation module is used to establish a finite element model in COMSOL using the received current geometric variables, to calculate the tension parameters on the center line of the longitudinal field coil, and to transmit the tension parameters on the center line of the longitudinal field coil to MATLAB. The loop module is used in MATLAB to determine whether the tension parameters and iteration count on the center line of the longitudinal field coil meet the iteration stopping condition. If not, the current geometric variables are updated using a genetic algorithm and transmitted to COMSOL, and then the calculation module is returned to re-execute. If the condition is met, the current geometric variables are output to design the fusion longitudinal field coil.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

Images