Electromagnetic, heat and stress bidirectional coupling analysis method in complex electromagnetic environment

Abstract

The invention discloses an electromagnetic, heat and stress bidirectional coupling analysis method in a complex electromagnetic environment, and the method comprises the steps: building a fusion reactor cladding coupling model based on three-dimensional modeling design software; establishing an electromagnetic field distribution analysis model, and calculating to obtain an electromagnetic field, current distribution, Joule heat and electromagnetic force; establishing a temperature field distribution analysis model, and calculating temperature field distribution in the cladding by solving a fluid-solid conjugate heat transfer equation; establishing a stress-strain analysis model, solving a solid mechanical stress-strain equation, obtaining stress and strain mechanical performance distribution of the cladding, performing cladding performance analysis, feeding back to the electromagnetic field distribution analysis model and the temperature field distribution analysis model, and performing iterative analysis calculation to obtain a final analysis result. According to the method, the defect that mutual influence among an electromagnetic field, a temperature field and a stress field is neglected in an existing analysis method is overcome, and the electromagnetic performance, the thermal hydraulic performance and the mechanical performance of the cladding in a complex electromagnetic environment are accurately analyzed.

Description

Technical Field

[0001] The present invention relates to the technical field of magnetic confinement fusion reactors, and in particular to an electromagnetic, thermal and stress bidirectional coupling analysis method under a complex electromagnetic environment. Background Art

[0002] Magnetic confinement fusion reactor is one of the main methods of controllable confinement fusion that can be used at present, among which the tokamak magnetic confinement device is the focus of current research. The main principle of the tokamak magnetic confinement fusion device is to confine and control plasma through a strong magnetic field. The entire device is a complex electromagnetic system. During the operation of the device, the device inside the vacuum chamber will always be under the influence of a strong magnetic field. The fusion reactor blanket is an important component in the vacuum chamber and is the key to the application of fusion energy. The blanket system mainly provides the following functions: breeding tritium, absorbing plasma radiation and particle heat flux, and providing shielding to reduce the heat and neutron radiation of the vacuum container and the components inside the container.

[0003] In a magnetic confinement fusion reactor, the blanket is inside the vacuum chamber system and is in a complex multi-physical field environment, involving physical fields, electromagnetic fields, temperature fields, flow fields, stress fields, etc., which interact and are deeply coupled. The blanket is subjected to the surface heat flux density from plasma radiation, as well as the nuclear thermal power density generated by the 14MeV fast neutrons produced by the fusion reaction interacting with the structural materials in the blanket. The blanket structure has a temperature field distribution, and the thermal stress generated by the temperature difference causes the stress field to change. Considering the mechanical loads such as the internal coolant pressure and the material's own gravity on the blanket, thermomechanical stress is generated.

[0004] However, when operating in a strong magnetic field environment, these thermal and mechanical performance parameters will have large deviations, affecting the correct judgment of structural safety. The candidate structural materials for the cladding are mostly low-activated ferrite / martensite (RAFM) steel. The electromagnetic effect between the ferromagnetic structural material and the strong magnetic field leads to a huge electromagnetic load, including static magnetization effect and transient electromagnetic induction effect, which brings additional electromagnetic force and torque to the cladding structure and is loaded during the entire operation of the cladding. The effect of the electromagnetic field changes the performance of the cladding, and the temperature field and stress field change greatly. The thermal-stress field in turn affects the electromagnetic field distribution and changes the size and distribution of the electromagnetic load. With the change of the temperature field distribution, the saturation magnetic induction intensity of the cladding structural material changes, which in turn affects the distribution of the electromagnetic field. At the same time, its conductivity also changes with temperature, which means that the temperature will affect the size of the induced current in the cladding module, thereby affecting the size and distribution of the electromagnetic field and the electromagnetic load of the cladding; the distribution of the stress field and strain field changes the eddy current distribution inside the structure, which has a feedback effect on the electromagnetic field. Summary of the invention

[0005] The purpose of the present invention is to overcome the shortcomings of the prior art. To achieve the above purpose, a bidirectional coupling analysis method of electromagnetic, thermal and stress in a complex electromagnetic environment is adopted to solve the problems raised in the above background technology.

[0006] A method for bidirectional coupling analysis of electromagnetic, thermal and stress in a complex electromagnetic environment includes the following steps:

[0007] Step S1, based on the three-dimensional modeling design software and according to the design scheme of the helium-cooled solid breeder blanket, a fusion reactor blanket coupling model is established, wherein the model includes a magnet system, a plasma, a blanket system, and a vacuum chamber;

[0008] Step S2, establishing an electromagnetic field distribution analysis model, and calculating the electromagnetic field, current distribution, Joule heat, and electromagnetic force;

[0009] Step S3, establishing a temperature field distribution analysis model, and calculating the temperature field distribution in the cladding by solving the fluid-solid conjugate heat transfer equation;

[0010] Step S4, establish a stress-strain analysis model, solve the solid mechanics stress-strain equation, obtain the stress and strain mechanical property distribution of the cladding, perform performance analysis of the cladding, and feed it back to the electromagnetic field distribution analysis model and the temperature field distribution analysis model for iterative analysis and calculation to obtain the final analysis result.

[0011] As a further solution of the present invention: the specific steps in step S1 include:

[0012] The design scheme of the helium-cooled solid breeder blanket uses the 3D modeling design software CATIA to establish a fusion reactor blanket coupling model, which includes a magnet system, plasma, blanket system, and vacuum chamber;

[0013] The cladding system includes a first wall, a protective layer, a reinforcing plate, a cover plate, a lithium silicate ball bed, a beryllium ball bed, a cooling plate, and a back plate;

[0014] Three complete sectors containing seven inner and eight outer envelopes around a typical envelope are modeled to capture the effect of ferromagnetism on the envelope.

[0015] As a further solution of the present invention: the specific steps in step S2 include:

[0016] Step 21, establishing an electromagnetic field distribution analysis model, obtaining the physical property parameters of the saturation magnetization intensity, relative magnetic permeability and resistivity of the loaded ferromagnetic material, and obtaining the influence of temperature on the physical properties of the material;

[0017] Step 22: Use the Mesh function of the Maxwell software to select an appropriate mesh size for division according to the spatial scale of each component;

[0018] Step 23, according to the design parameters, the current density on the cross section of the magnet coil and the plasma is uniformly loaded;

[0019] Step 24, setting the solver to magnetic transient, and evaluating the influence of induced current on electromagnetic force based on eddy current effect;

[0020] Step 25: Set the calculation time according to the analysis conditions and adjust the time step to ensure the accuracy and convergence of the calculation results;

[0021] Step 26: Finally, the electromagnetic field and current distribution are calculated, and Joule heat and electromagnetic force are obtained.

[0022] As a further solution of the present invention: the specific steps in step S3 include:

[0023] Step 31, establish and simplify the temperature field distribution analysis model, load the physical property parameters of the materials including the structure, coolant, and proliferation agent, and determine the density, thermal conductivity, and specific heat capacity for research based on the temperature influence;

[0024] Step 32, using the Fluid Line function of the simulation software to simulate the coolant fluid, and to perform heat exchange between the coolant and the structural material through coupling between the fluid and the wall unit;

[0025] Step 33, loading neutron reaction nuclear heat deposition, plasma heat flux, and electromagnetic Joule heat heat sources as body power density and heat flux boundary conditions;

[0026] Step 34, setting the contact thermal resistance between the ball bed and the structural material, using the Yagi-Kunii model, and setting the coolant inlet and outlet boundaries according to the design parameters;

[0027] Step 35: Use Ansys Mesh module to build a thermal analysis grid and encrypt the key areas of the flow channel;

[0028] Step 36: Calculate the temperature distribution and save the results.

[0029] As a further solution of the present invention: the specific steps in step S4 include:

[0030] Step 41: Establish a stress-strain analysis model, and load material property parameters including thermal expansion coefficient, Young's modulus, and Poisson's ratio;

[0031] Step 42: Set constraints on the upper and lower surfaces of the analysis model and behind the back plate to prevent the model from unexpected displacement or rotation;

[0032] Step 43: The temperature field obtained by thermal analysis, the electromagnetic force field obtained by electromagnetic analysis, and the mechanical load caused by the high-pressure coolant helium and the purge gas are loaded as loads into the stress-strain analysis model;

[0033] Step 44: Solve the solid mechanics stress-strain equation, analyze the cladding performance, and feed back the deformation to the electromagnetic field distribution analysis model and the temperature field distribution analysis model, perform iterative analysis and calculation, and obtain the final analysis result.

[0034] Compared with the prior art, the present invention has the following technical effects:

[0035] Adopting the above technical scheme, through the complex environment of strong magnetic field, high heat flux density and internal coolant pressure faced by the fusion reactor cladding, a processing method suitable for the bidirectional coupling analysis of electromagnetic field, temperature field and stress field is proposed, a basic coupling scheme is formulated, and the coupling interface of each physical field solving program is developed. Based on the unified data mapping format and the unified ICoCo packaging concept, the process control, data interaction and mapping between programs, and the coupling iterative algorithm are studied;

[0036] Based on the coupling and cross-fertilization of multiple physical fields, the electromagnetic field, temperature field and stress field are coupled, and the mutual influence between the electromagnetic field, temperature field and stress field is quantified to obtain more accurate numerical simulation results;

[0037] The mutual influence between electromagnetics, heat and stress is taken into account at the same time, and the results obtained are more in line with the operating status of the blanket in the actual environment, providing a reference for the design and manufacture of future fusion reactor blankets. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] The specific embodiments of the present invention are described in detail below in conjunction with the accompanying drawings:

[0039] Figure 1 A schematic diagram of the steps of the bidirectional coupling analysis method disclosed in the present application;

[0040] Figure 2 A flow chart of a bidirectional coupling analysis method according to an embodiment disclosed in this application;

[0041] Figure 3 A schematic diagram of a fusion reactor model and a blanket model for an embodiment disclosed in the present application;

[0042] Figure 4 A schematic diagram of an electromagnetic analysis grid according to an embodiment disclosed in the present application;

[0043] Figure 5 A schematic diagram of a temperature field distribution analysis model of an embodiment disclosed in this application;

[0044] Figure 6This is a schematic diagram of a stress analysis model of an embodiment disclosed in this application. DETAILED DESCRIPTION

[0045] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0046] Please refer to Figure 1 and Figure 2 In an embodiment of the present invention, a method for analyzing electromagnetic, thermal and stress bidirectional coupling in a complex electromagnetic environment includes the following steps:

[0047] Step S1, based on the three-dimensional modeling design software and according to the design scheme of the helium-cooled solid breeder blanket, a fusion reactor blanket coupling model is established, wherein the model includes a magnet system, a plasma, a blanket system, and a vacuum chamber.

[0048] According to the above steps, one of the specific implementation methods is that the electromagnetic-thermal-mechanical bidirectional coupling scheme adopts the Supervisor unified platform solution based on the Salome platform. Each physical program needs to develop a corresponding coupling interface based on the ICoCo (Interface for Code Coupling) concept, and use the MED library (Model Exchange Data) to support the program to store grid / field data in memory or files. The Salome platform serves as the Supervisor main program to control the operation of each physical program and the data transfer between programs. The ICoCo interface concept requires that each program can implement a series of pre-defined interface functions to complete the corresponding functions respectively. These functions include: initializing the program, calculating the time step, initializing the calculation with the time step, solving / iterative solving the time step, extracting the output field, setting the input field and terminating the program. Data exchange between programs is an important part of program coupling. The input and output data may be large-capacity, complex and diverse in data structure. The MED library solution is used for data transfer; this solution supports the program to store grid / field data in memory or files, and supports field data mapping on non-matching grids.

[0049] Step S1 specifically includes:

[0050] The design scheme of the helium-cooled solid breeder blanket uses the 3D modeling design software CATIA to establish a fusion reactor blanket coupling model, which includes a magnet system, plasma, blanket system, and vacuum chamber;

[0051] The cladding system includes a first wall, a protective layer, a reinforcing plate, a cover plate, a lithium silicate ball bed, a beryllium ball bed, a cooling plate, and a back plate;

[0052] Three complete sectors containing seven inner and eight outer envelopes around a typical envelope are modeled to capture the effect of ferromagnetism on the envelope.

[0053] like Figure 3 As shown, the figure is a schematic diagram of fusion reactor modeling and blanket model;

[0054] One of the specific implementation methods is to establish a fusion reactor blanket coupling model based on the design plan of the helium-cooled solid breeder blanket of the fusion engineering test reactor using the three-dimensional modeling design software CATIA. The model includes a magnet system (including a longitudinal field magnet system, a poloidal field magnet system and a correction field magnet system), plasma, a blanket system, and a vacuum chamber.

[0055] The cladding system includes a first wall, a protective layer, a reinforcement plate, a cover plate, a lithium silicate ball bed, a beryllium ball bed, a cooling plate, and a back plate.

[0056] The finely modeled equatorial envelope is taken as a typical object for analysis. In order to consider the influence of ferromagnetism on the envelope, three complete sectors containing seven inner envelopes and eight outer envelopes are modeled around the typical envelope.

[0057] In order to simplify the calculation, the cladding other than the typical cladding on the equatorial plane adopts a simplified structure, that is, no structural division is made within the cladding, and the material properties are loaded as the properties of the cladding components converted in proportion.

[0058] Step S2: Establishing an electromagnetic field distribution analysis model to calculate the electromagnetic field, current distribution, Joule heat, and electromagnetic force. The specific steps include:

[0059] Calculate the electromagnetic field, induced current and electromagnetic force distribution of the cladding, and use the interactive interface of Maxwell software to set them up;

[0060] Step 21, establishing an electromagnetic field distribution analysis model, obtaining the physical property parameters of the saturation magnetization intensity, relative magnetic permeability and resistivity of the loaded ferromagnetic material, and obtaining the influence of temperature on the physical properties of the material;

[0061] like Figure 4 As shown, the figure is a schematic diagram of the electromagnetic analysis grid;

[0062] The specific implementation method is to load material property parameters. The structural material is a ferromagnetic material, and the saturation magnetization intensity, relative magnetic permeability and resistivity need to be given. The material properties are involved in considering the influence of temperature.

[0063] Step 22: Use the Mesh function of the Maxwell software to select an appropriate mesh size for division according to the spatial scale of each component, specifically:

[0064] Use Maxwell's Mesh function to divide the mesh and select the appropriate mesh size according to the spatial scale of different components;

[0065] Step 23: According to the design parameters, the current density on the cross section of the magnet coil and the plasma is uniformly loaded, specifically:

[0066] According to the design parameters, the magnet coil and plasma are loaded with current, and the current density is evenly distributed on the cross section. The coil current parameters under the typical plasma configuration are shown in the following table:

[0067] Coil Number of turns Single turn current (kA) Total current (kA) CS1U 374 -56 -20944 CS2U 374 -13 -4862 CS3U 374 13 4862 CS1L 374 -56 -20944 CS2L 374 4 1496 CS3L 374 13 4862 PF1U 616 15 9240 PF2U 324 -6 -1944 PF3U 324 -14 -4536 PF1L 616 20 12320 D1 324 0 0 D2 324 0 0 PF2L 324 -1 -324 PF3L 324 -20 -6480 TF 132 67.5 8910

[0068] Step 24: Set the solver to magnetic transient and evaluate the influence of induced current on electromagnetic force based on eddy current effect. Specifically:

[0069] The force and torque parameters of the cladding are set as eddy current effector to consider the influence of induced current on electromagnetic force;

[0070] Step 25: Set the calculation time according to the analysis conditions and adjust the time step to ensure the accuracy and convergence of the calculation results. Specifically:

[0071] Set the calculation time and time step. The transient time refers to the working condition to be analyzed, and the time step is adjusted until it does not affect the calculation results and convergence.

[0072] Step 26: Finally, the electromagnetic field and current distribution are calculated, and Joule heat and electromagnetic force are obtained, specifically:

[0073] The electromagnetic field and current distribution are calculated, and the Joule heat and electromagnetic force are calculated, extracted and saved to provide data for subsequent temperature field and stress field calculations.

[0074] Step S3: Establish a temperature field distribution analysis model, and calculate the temperature field distribution in the cladding by solving the fluid-solid conjugate heat transfer equation. The specific steps include:

[0075] like Figure 5 As shown, the figure is a schematic diagram of the thermal analysis model;

[0076] Step 31: Establish and simplify the temperature field distribution analysis model, load the physical parameters of the materials including the structure, coolant, and proliferation agent, and determine the density, thermal conductivity, and specific heat capacity based on the temperature influence for research, specifically:

[0077] Load material physical parameters, including structural materials, coolant materials, proliferation materials, etc. Given material density, thermal conductivity, specific heat capacity and other physical parameters, consider the influence of temperature;

[0078] Step 32: Use the Fluid Line function of the simulation software to simulate the coolant fluid, and conduct heat exchange between the coolant and the structural material through the coupling between the fluid and the wall unit. Specifically:

[0079] Use the Fluid Line function to simulate the coolant fluid, and simulate the convection heat transfer between the fluid and the wall through the coupling between the fluid unit and the surface unit, so as to realize the heat exchange between the coolant fluid and the structural material;

[0080] Step 33: Load the heat sources of neutron reaction nuclear heat deposition, plasma heat flux, and electromagnetic Joule heat as the body power density and heat flux boundary conditions, specifically:

[0081] Loading heat sources, including nuclear heat deposition from neutron-material reactions, loaded into the structure as body power density based on the results of neutron calculations, and heat flux from the plasma loaded on a first wall facing the plasma as a heat flux boundary, and Joule heat obtained from electromagnetic calculations as body power density;

[0082] Step 34: Set the contact thermal resistance between the ball bed and the structural material, adopt the Yagi-Kunii model, and set the coolant inlet and outlet boundaries according to the design parameters, specifically:

[0083] Set the contact thermal resistance between the ball bed and the structural material. The relationship adopts the internationally commonly used Yagi-Kunii model. Set the coolant inlet and outlet boundaries. The inlet flow rate and temperature are based on the design parameters.

[0084] Step 35: Use Ansys Mesh module to build thermal analysis mesh and encrypt the key areas of the flow channel. Specifically:

[0085] Use Ansys Mesh module to build thermal analysis mesh, set body size and encrypt areas such as flow channels;

[0086] Step 36: Calculate the temperature distribution and save the results, specifically:

[0087] The temperature distribution is calculated, extracted and saved to provide material property calculation input for electromagnetic and stress calculations, and to provide input for stress analysis to calculate thermal stress.

[0088] Step S4, establish a stress-strain analysis model, solve the solid mechanics stress-strain equation, obtain the stress and strain mechanical property distribution of the cladding, perform cladding performance analysis, and feed back to the electromagnetic field distribution analysis model and the temperature field distribution analysis model, perform iterative analysis and calculation, and obtain the final analysis results. The specific steps include:

[0089] like Figure 6 As shown, the figure is a schematic diagram of the stress analysis model;

[0090] In this embodiment, since the mechanical interaction between the ball bed and the wall is extremely small and the protective layer does not bear any mechanical function, the ball bed and the protective layer do not need to be considered in the stress analysis model, and the rest is consistent with the thermal analysis model. The specific stress analysis steps are:

[0091] Step 41: Establish a stress-strain analysis model, and load material property parameters including thermal expansion coefficient, Young's modulus, and Poisson's ratio;

[0092] Step 42: Set constraints on the upper and lower surfaces of the analysis model and behind the back plate to prevent the model from unexpected displacement or rotation. Specifically:

[0093] Load constraint conditions. Load coupling constraint boundary conditions on the upper surface of the computational model to make all nodes on the upper surface move in parallel in the polar direction at the same time. Use symmetric constraint boundary conditions on the lower surface. The four faces behind the back plate use constraint conditions that the radial displacement is 0 (that is, the back plate cannot be displaced in the radial direction). Select two nodes on the first wall and make the circumferential displacement of these two nodes 0 to prevent the computational model from rotating.

[0094] Step 43: The temperature field obtained by thermal analysis, the electromagnetic force field obtained by electromagnetic analysis, and the mechanical load caused by the high-pressure coolant helium and the purge gas are loaded as loads into the stress-strain analysis model, specifically:

[0095] Loading loads include thermal loads, mechanical loads and electromagnetic loads. The temperature field obtained by thermal analysis is imported to load thermal loads, and the electromagnetic force field obtained by electromagnetic analysis is imported to load electromagnetic loads. Mechanical loads are mainly caused by the high-pressure coolant helium in the flow channels of the cooling components of the blanket and the purge gas in the pebble bed, and are distributed and loaded on the inner wall of the flow channel and the wall surface of the components in contact with the purge gas.

[0096] Step 44: Solve the solid mechanics stress-strain equation, analyze the cladding performance, and feed back the deformation to the electromagnetic field distribution analysis model and the temperature field distribution analysis model, perform iterative analysis and calculation, and obtain the final analysis results. One of the specific implementation steps is:

[0097] Solve the solid mechanics stress-strain equations to obtain the stress and strain mechanical property distribution of the cladding, perform cladding performance analysis, and feed the deformation back to the electromagnetic analysis and thermal analysis models for iterative calculations.

[0098] Although the embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention. The scope of the present invention is defined by the appended claims and their equivalents, and all should be included within the scope of protection of the present invention.

Claims

1. A method for bidirectional coupling analysis of electromagnetic, thermal and stress in a complex electromagnetic environment, characterized in that: The following steps are involved: Step S1, based on the three-dimensional modeling design software and according to the design scheme of the helium-cooled solid breeder blanket, a fusion reactor blanket coupling model is established, wherein the model includes a magnet system, a plasma, a blanket system, and a vacuum chamber; Step S2, establishing an electromagnetic field distribution analysis model, and calculating the electromagnetic field, current distribution, Joule heat, and electromagnetic force; Step S3, establishing a temperature field distribution analysis model, and calculating the temperature field distribution in the cladding by solving the fluid-solid conjugate heat transfer equation; Step S4, establish a stress-strain analysis model, solve the solid mechanics stress-strain equation, obtain the stress and strain mechanical property distribution of the cladding, perform performance analysis of the cladding, and feed it back to the electromagnetic field distribution analysis model and the temperature field distribution analysis model for iterative analysis and calculation to obtain the final analysis result.

2. According to the method of electromagnetic, thermal and stress bidirectional coupling analysis under complex electromagnetic environment in claim 1, it is characterized in that: The specific steps in step S1 include: The design scheme of the helium-cooled solid breeder blanket uses the 3D modeling design software CATIA to establish a fusion reactor blanket coupling model, which includes a magnet system, plasma, blanket system, and vacuum chamber; The cladding system includes a first wall, a protective layer, a reinforcing plate, a cover plate, a lithium silicate ball bed, a beryllium ball bed, a cooling plate, and a back plate; Three complete sectors containing seven inner and eight outer envelopes around a typical envelope are modeled to capture the effect of ferromagnetism on the envelope.

3. According to the method of electromagnetic, thermal and stress bidirectional coupling analysis under complex electromagnetic environment in claim 1, it is characterized in that: The specific steps in step S2 include: Step 21, establishing an electromagnetic field distribution analysis model, obtaining the physical property parameters of the saturation magnetization intensity, relative magnetic permeability and resistivity of the loaded ferromagnetic material, and obtaining the influence of temperature on the physical properties of the material; Step 22: Use the Mesh function of the Maxwell software to select an appropriate mesh size for division according to the spatial scale of each component; Step 23, according to the design parameters, the current density on the cross section of the magnet coil and the plasma is uniformly loaded; Step 24, setting the solver to magnetic transient, and evaluating the influence of induced current on electromagnetic force based on eddy current effect; Step 25: Set the calculation time according to the analysis conditions and adjust the time step to ensure the accuracy and convergence of the calculation results; Step 26: Finally, the electromagnetic field and current distribution are calculated, and Joule heat and electromagnetic force are obtained.

4. According to the method of electromagnetic, thermal and stress bidirectional coupling analysis under complex electromagnetic environment in claim 1, it is characterized in that: The specific steps in step S3 include: Step 31, establish and simplify the temperature field distribution analysis model, load the physical property parameters of the materials including the structure, coolant, and proliferation agent, and determine the density, thermal conductivity, and specific heat capacity for research based on the temperature influence; Step 32, using the Fluid Line function of the simulation software to simulate the coolant fluid, and to perform heat exchange between the coolant and the structural material through coupling between the fluid and the wall unit; Step 33, loading neutron reaction nuclear heat deposition, plasma heat flux, and electromagnetic Joule heat heat sources as body power density and heat flux boundary conditions; Step 34, setting the contact thermal resistance between the ball bed and the structural material, using the Yagi-Kunii model, and setting the coolant inlet and outlet boundaries according to the design parameters; Step 35: Use Ansys Mesh module to build a thermal analysis grid and encrypt the key areas of the flow channel; Step 36: Calculate the temperature distribution and save the results.

5. According to claim 1, a method for analyzing electromagnetic, thermal and stress bidirectional coupling under complex electromagnetic environment, characterized in that: The specific steps in step S4 include: Step 41: Establish a stress-strain analysis model, and load material property parameters including thermal expansion coefficient, Young's modulus, and Poisson's ratio; Step 42: Set constraints on the upper and lower surfaces of the analysis model and behind the back plate to prevent the model from unexpected displacement or rotation; Step 43: The temperature field obtained by thermal analysis, the electromagnetic force field obtained by electromagnetic analysis, and the mechanical load caused by the high-pressure coolant helium and the purge gas are loaded as loads into the stress-strain analysis model; Step 44: Solve the solid mechanics stress-strain equation, analyze the cladding performance, and feed back the deformation to the electromagnetic field distribution analysis model and the temperature field distribution analysis model, perform iterative analysis and calculation, and obtain the final analysis result.

Images