An automated numerical simulation method, device and equipment for inductive heating type crystal growth

Abstract

This invention provides an automated numerical simulation method, apparatus, and device for induction heating crystal growth, belonging to the field of crystal growth technology. It aims to automate the numerical simulation of crystal growth. The method includes: acquiring a first macro code file, a second macro code file, and a third macro code file of a target model; receiving multiple combinations of operating parameters; and for each combination of operating parameters, executing steps to simulate the growth environment in crystal growth based on the target model, including: modifying the model parameters in the first, second, and third macro code files respectively based on each combination of operating parameters; and sequentially executing the modified first macro code file, the modified second macro code file, and the modified third macro code file in a preset order to obtain simulation results corresponding to each combination of operating parameters. The simulation results include the temperature of the induction furnace and the fluid flow rate.

Description

Technical Field

[0001] This invention relates to the field of crystal growth technology, and in particular to an automated numerical simulation method, apparatus and equipment for induction heating crystal growth. Background Technology

[0002] Silicon carbide (SiC) is a third-generation wide-bandgap semiconductor material with excellent properties such as high breakdown field strength, high thermal conductivity, high saturated electron mobility, high chemical stability, and radiation resistance, making it a promising candidate for high-power electronic devices. The most commonly used growth technique for SiC single crystals is the physical vapor transport (PVT) method. The PVT method utilizes a temperature gradient within the crucible (high temperature at the bottom, low temperature at the top) to transport the SiC vapor obtained from the sublimation of SiC powder at the bottom of the crucible to a lower-temperature SiC seed crystal at the top for recrystallization. Over time, a SiC ingot of a certain thickness crystallizes on the seed crystal. Heating methods for SiC powder are divided into induction and resistance methods. Induction heating uses a high-frequency alternating current passed through an induction coil to inductively heat the crucible and the surrounding graphite cylinder, then heats the SiC powder inside the graphite crucible through heat conduction and radiation. Temperature field changes have a significant impact on SiC vapor transport and ingot defects; however, because SiC crystal growth occurs within a high-temperature, sealed container, it is impossible to directly observe changes in the internal temperature field and gas flow.

[0003] Currently, numerical simulation technology is used to calculate all temperature and flow field information during the growth of SiC crystals. However, the growth process of SiC crystals by induction heating is complex and requires coupling multiple software programs to achieve accurate numerical simulation. When faced with the need to study parameter sensitivity or optimize design using a large amount of simulation data, the existing manual numerical simulation process is often time-consuming, labor-intensive, and inefficient. Summary of the Invention

[0004] Based on the background technology, this invention proposes an automated numerical simulation method, apparatus, and equipment.

[0005] In a first aspect, the present invention provides an automated numerical simulation method for induction heating crystal growth, the method comprising:

[0006] Obtain a first macro code file, a second macro code file, and a third macro code file of the target model; wherein, the first macro code file is used to obtain the mesh file of the target model, the second macro code file is used to obtain the power density of the induction furnace in the target model, and the third macro code file is used to obtain the temperature and fluid flow velocity of the induction furnace in the target model;

[0007] It can receive multiple combinations of operating parameters, and at least one of the operating parameters has a different value range among the different combinations of operating parameters;

[0008] For each combination of operating parameters, a step of simulating the growth environment in crystal growth is performed based on the target model. The step includes: modifying the model parameters in the first macro code file, the second macro code file, and the third macro code file respectively based on each combination of operating parameters.

[0009] According to a preset order, the modified first macro code file, the modified second macro code file, and the modified third macro code file are executed sequentially to obtain simulation results corresponding to each combination of operating parameters. The simulation results include the temperature of the induction furnace and the fluid flow rate.

[0010] Optionally, the step of executing the modified first macro code file, the modified second macro code file, and the modified third macro code file sequentially according to a preset order includes:

[0011] The modified first macro code file is executed to obtain the updated mesh file, and based on the updated mesh file, the modified third macro code file is executed to obtain the node coordinate information of the target model;

[0012] The node coordinate information is loaded into the modified second macro code file, and the modified second macro code file is executed to obtain the node power density corresponding to each node of the induction furnace in the target model.

[0013] The power density of multiple nodes is loaded into the modified third macro code file, and the modified third macro code file is executed to obtain the temperature and fluid flow velocity corresponding to each node of the induction furnace in the target model.

[0014] Optionally, when the power parameters included in the target model are modified, the step of executing the modified second macro code file based on the node coordinate information to obtain the node power density corresponding to each node of the induction furnace in the target model includes:

[0015] Determine the induced current based on the power parameters;

[0016] Based on the induced current and the node coordinate information, the modified second macro code file is executed to obtain the actual power density corresponding to each node of the induction furnace in the target model, and the actual power value is determined according to the multiple actual power densities.

[0017] Match the actual power value with the value of the power parameter;

[0018] When the actual power value does not match the value of the power parameter, the induced current is corrected and the actual power value is re-determined based on the difference between the actual power value and the power parameter.

[0019] When the actual power value matches the value of the power parameter, the actual power density corresponding to the induced current is determined as the node power density corresponding to each node of the induction furnace in the target model.

[0020] Optionally, loading the power densities of the multiple nodes into the modified third macro code file and executing the modified third macro code file to obtain the temperature and fluid flow velocity corresponding to each node of the induction furnace in the target model includes:

[0021] Obtain the first preset node coordinates, which are the node coordinates of the target model after updating the mesh file;

[0022] When executing the modified third macro code file, the temperature and fluid velocity corresponding to the first preset node coordinates are obtained at preset time intervals, and it is determined whether the temperature field and fluid flow field of the target model have both reached a steady state; the determination of whether the temperature field and fluid flow field of the target model have both reached a steady state includes: determining the first difference between the temperature obtained at the current time and the temperature obtained at the previous time, and the second difference between the fluid velocity obtained at the current time and the fluid velocity obtained at the previous time.

[0023] Based on the first difference and the second difference, determine whether the temperature field and fluid flow field of the target model have both reached a steady state;

[0024] If neither the temperature field nor the fluid flow field of the target model has reached a steady state, the steps of obtaining the temperature and fluid velocity corresponding to the first preset node coordinates at preset time intervals and determining whether the temperature field and the fluid flow field of the target model have both reached a steady state continue to be executed.

[0025] When the temperature field and fluid flow field of the target model reach a steady state, the temperature and fluid flow velocity corresponding to the current time are determined as the temperature and fluid flow velocity corresponding to each node of the induction furnace after the target model updates the mesh file.

[0026] Optionally, determining whether the temperature field and fluid flow field of the target model have both reached a steady state based on the first difference and the second difference includes:

[0027] Compare the first difference with a preset temperature threshold, and the second difference with a preset speed threshold;

[0028] If the first difference is less than the preset temperature threshold and the second difference is less than the preset velocity threshold, it is determined that the temperature field and fluid flow field of the target model have both reached a steady state.

[0029] Optionally, the method further includes:

[0030] Obtain at least one second preset node coordinate, wherein the second preset node coordinate is the node coordinate that affects crystal growth in the target model;

[0031] The temperature and fluid flow velocity corresponding to the second preset node coordinates are determined as the simulation results corresponding to the combination of operating parameters, and a target list is obtained. The target list includes multiple combinations of operating parameters and the simulation results corresponding to each combination of operating parameters.

[0032] Optionally, the method further includes:

[0033] Based on the number of the operating condition parameter combinations, the multiple operating condition parameter combinations are divided into multiple batches;

[0034] For each combination of operating parameters in the same batch, the step of simulating the growth environment in crystal growth is performed simultaneously.

[0035] Optionally, modifying the model parameters in the first macro code file, the second macro code file, and the third macro code file based on each combination of operating condition parameters includes:

[0036] Regular expressions are used to determine the target model parameters from the first macro code file, the second macro code file, and the third macro code file, respectively; wherein different macro code files correspond to slightly different target model parameters.

[0037] Determine the target parameter value of the working condition parameter that matches the target model parameter in the working condition parameter combination;

[0038] Modify the numerical values ​​of the target model parameters to the target parameter values.

[0039] A second aspect of the present invention provides an automated numerical simulation apparatus for induction heating crystal growth, the apparatus comprising:

[0040] The acquisition module is used to acquire a first macro code file, a second macro code file, and a third macro code file of the target model; wherein, the first macro code file is used to acquire the mesh file of the target model, the second macro code file is used to acquire the power density of the induction furnace in the target model, and the third macro code file is used to acquire the temperature and fluid flow velocity of the induction furnace in the target model;

[0041] The receiving module is used to receive multiple combinations of operating parameters, where at least one operating parameter has a different value range among the different combinations of operating parameters.

[0042] The simulation module is used to perform a step of simulating the growth environment in crystal growth based on the target model for each combination of operating parameters. The step includes: modifying the model parameters in the first macro code file, the second macro code file, and the third macro code file respectively based on each combination of operating parameters; and executing the modified first macro code file, the modified second macro code file, and the modified third macro code file in a preset order to obtain simulation results corresponding to each combination of operating parameters. The simulation results include the temperature of the induction furnace and the fluid flow rate.

[0043] A third aspect of the present invention provides an electronic device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executed, performs the steps in the automated numerical simulation of induction heating crystal growth as described in the first aspect.

[0044] The present invention provides a method comprising: acquiring a first macro code file, a second macro code file, and a third macro code file of a target model; wherein the first macro code file is used to acquire a mesh file of the target model, the second macro code file is used to acquire the power density of an induction furnace in the target model, and the third macro code file is used to acquire the temperature and fluid flow velocity of the induction furnace in the target model; receiving multiple combinations of operating parameters, wherein at least one operating parameter has a different numerical range among different combinations of operating parameters; and, for each combination of operating parameters, performing a step of simulating the growth environment in crystal growth based on the target model, the step comprising: modifying the model parameters in the first macro code file, the second macro code file, and the third macro code file respectively based on each combination of operating parameters; and sequentially executing the modified first macro code file, the modified second macro code file, and the modified third macro code file in a preset order to obtain simulation results corresponding to each combination of operating parameters, the simulation results including the temperature and fluid flow velocity of the induction furnace;

[0045] Therefore, this invention automates the numerical simulation process of the target model by using a first macro code file to obtain the mesh file of the target model, a second macro code file to obtain the power density of the induction furnace in the target model, and a third macro code file to obtain the temperature and fluid flow velocity of the induction furnace in the target model. Thus, when modifying the numerical values ​​of the operating parameters, the simulation results for the target model can be obtained by modifying the model parameters in the first, second, and third macro code files and then re-executing the first, second, and third macro code files, without requiring manual parameter modification, achieving a highly efficient numerical simulation process.

[0046] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description

[0047] To more clearly illustrate the technical solutions in the embodiments of the present invention or related technologies, the accompanying drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. It should be noted that the scale in the drawings is for illustration only and does not represent the actual scale.

[0048] Figure 1 A flowchart illustrating the steps of the automated numerical simulation method for induction heating crystal growth provided in an embodiment of the present invention is shown.

[0049] Figure 2 A schematic diagram of the crystal growth environment in an embodiment of the present invention is shown;

[0050] Figure 3 A flowchart illustrating the steps of simulating the crystal growth environment based on a target model in an embodiment of the present invention is shown.

[0051] Figure 4 This diagram illustrates the flowchart of an automated numerical simulation method for induction heating crystal growth provided in an embodiment of the present invention.

[0052] Figure 5 This invention illustrates a flowchart of the execution of the modified second macro code file in an embodiment of the invention.

[0053] Figure 6 This invention illustrates a flowchart of the execution of the modified third macro code file in an embodiment of the present invention.

[0054] Figure 7 A schematic diagram of the structure of the automated numerical simulation device for induction heating crystal growth provided in an embodiment of the present invention is shown. Detailed Implementation

[0055] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0056] The growth process of SiC crystals takes place in a high-temperature, sealed container, making it impossible to directly observe the changes in the internal temperature field and gas flow. Therefore, numerical simulation techniques are currently commonly used to calculate all temperature and flow field information during the growth process of SiC crystals.

[0057] However, existing numerical simulations of SiC induction heating crystal growth mostly rely on various software programs for manual model building, mesh creation, boundary condition setting, and parameter input. This is not only inefficient but also prone to human error, making it difficult to meet the demands of large-scale simulation tasks with multiple parameters and operating conditions. In existing technologies, due to the complexity of SiC induction heating furnace models, when furnace or process parameters need to be modified, manually drawing and modifying the model mesh using traditional methods requires a significant amount of time to modify each point, line, and surface element, redraw the mesh, reload the mesh, assign material parameters, boundary conditions, and set example parameters. Similarly, when modifying SiC induction heating power parameters, a considerable amount of time is required to modify the current on each coil and recalculate the power density to achieve the desired power.

[0058] Furthermore, traditional simulation methods typically focus on single computational tasks, lacking comprehensive automated control mechanisms and making it difficult to achieve batch simulations and automatic result extraction under multiple sets of process parameters. Since numerical simulations of SiC induction heating crystal growth usually involve using Gambit to create the mesh, Maxwell to calculate the induction power density, and Fluent to calculate heat transfer and flow within the SiC furnace, multiple software programs need to be coupled together for computation. Traditional manual methods require modifying the model in Maxwell to calculate the induction power density, modifying the model and mesh in Gambit, importing the data into Fluent to modify boundary conditions and material parameters, calculating again, and then exporting the results. Due to the complexity of file interaction and processing between multiple software programs, traditional manual methods can only achieve software interaction by manually importing, modifying, and exporting these files one by one. Moreover, when Fluent is calculating, the criteria for convergence often rely on human experience, making accurate judgment difficult. This hinders automated integration and interaction among these software programs, resulting in a lack of high-throughput computing capabilities in manual methods.

[0059] Furthermore, post-processing after numerical simulation requires analysis of furnace temperature variations and gas flow patterns. Temperature field variations primarily include the top and bottom temperatures of the furnace and the radial temperature gradient of the seed crystal. Gas flow analysis necessitates obtaining the magnitude and direction of the argon gas flow velocity. Traditional manual methods require adding different anchor points for each case to derive the corresponding temperature points, and then deriving the flow velocity from the gas region. Subsequently, before studying parameter sensitivity, significant time is spent integrating process parameters and results from all operating conditions. This manual approach makes post-processing output and data organization time-consuming and labor-intensive, hindering data management and subsequent analysis, and severely restricting downstream tasks such as parameter sensitivity analysis, optimization design, and machine learning-driven modeling.

[0060] In view of this, the present invention provides an automated numerical simulation method, apparatus and equipment for induction heating crystal growth. By recording multiple macro code files and executing the macro code files in a preset order, the numerical simulation process is automated. When it is necessary to modify the operating parameters of the target model, the parameter values ​​to be modified can be directly input, and the adjustment and numerical simulation can be automatically realized, thereby improving the efficiency of numerical simulation.

[0061] Reference Figure 1 , Figure 1 This invention illustrates an automated numerical simulation method for induction heating crystal growth, as provided in an embodiment of the present invention. Figure 1 As shown, the method specifically includes:

[0062] Step S101: Obtain the first macro code file, the second macro code file, and the third macro code file of the target model.

[0063] The first macro code file is used to obtain the mesh file of the target model, the second macro code file is used to obtain the power density of the induction furnace in the target model, and the third macro code file is used to obtain the temperature and fluid flow velocity of the induction furnace in the target model.

[0064] In this embodiment, a first macro code file, a second macro code file, and a third macro code file can be pre-written or recorded. Different macro codes can be recorded using different software. For example, the first macro code file is used to obtain the mesh file of the target model. It can be recorded using macros provided by Gambit software and the target model to modify model parameters and redraw the mesh. The second macro code file is used to obtain the power density of the induction furnace in the target model. It can be recorded using macros provided by Maxwell software to modify model parameters and physical property parameters and calculate the power density of the induction furnace in the target model. The third macro code file is used to obtain the temperature and fluid flow velocity of the induction furnace in the target model. It can be recorded using macros provided by Fluent to modify model parameters and physical property parameters and calculate the temperature and fluid flow velocity of the induction furnace in the target model.

[0065] Step S102: Receive multiple combinations of operating parameters, where at least one operating parameter has a different value range among the different combinations of operating parameters.

[0066] Specifically, multiple combinations of operating parameters can be directly input by the user. Of course, the user can also input the value range of each operating parameter. In this case, a list of operating parameter combinations can be generated based on the value range of each operating parameter, resulting in multiple combinations of operating parameters.

[0067] The target model is used to simulate the induction heating crystal growth process, with reference to... Figure 2 , Figure 2 A cross-sectional view of an induction furnace for growing SiC crystals in related technologies is shown, such as... Figure 2 As shown, the induction heating furnace includes a crucible lid 1, a seed crystal 2, a soft felt 3, a graphite heating cylinder 4, a crucible 5, an induction coil 6, a porous graphite plate 7, and SiC powder 8. The induction heating furnace heats the crucible 5 and the graphite heating cylinder 4 outside the crucible by passing a high-frequency alternating current through the induction coil 6. Then, the SiC powder 8 inside the graphite crucible is heated through heat conduction and radiation. As the temperature rises and the soft felt 3 provides insulation, the SiC powder 8 sublimates into SiC gas, which passes through the porous graphite plate 7 and is transported to the lower-temperature SiC seed crystal 2 on the crucible lid 1 for recrystallization. Over time, a SiC ingot of a certain thickness crystallizes on the seed crystal.

[0068] Therefore, the parameters for simulating crystal growth can include power, bottom felt thickness, coil axial movement distance, graphite cap axial distance, powder height, top hole size, graphite conductivity, powder thermal conductivity, felt thermal conductivity, and graphite thermal conductivity. In other words, the operating parameters can include power, bottom felt thickness, coil axial movement distance, graphite cap axial distance, powder height, top hole size, graphite conductivity, powder thermal conductivity, felt thermal conductivity, and graphite thermal conductivity.

[0069] Step S103: For each combination of operating parameters, perform the steps of simulating the growth environment in crystal growth based on the target model.

[0070] Specifically, refer to Figure 3 , Figure 3 The following steps are illustrated in an embodiment of the present invention for simulating the growth environment during crystal growth based on a target model: Figure 3 As shown, step S103 may include:

[0071] Step S1031: Based on each combination of operating condition parameters, modify the model parameters in the first macro code file, the second macro code file, and the third macro code file respectively.

[0072] In this embodiment, if the value of at least one of the multiple operating parameters included in the operating parameter combination differs from that in the target model, it is necessary to first modify the model parameters in the first macro code file, the second macro code file, and the third macro code file according to the value of the operating parameter in the operating parameter combination. This is to modify the mesh file of the target model, the power density corresponding to the target model, and the temperature and fluid flow velocity corresponding to the target model according to the operating parameter combination, so as to obtain the heat transfer and flow field information of the target model under each operating parameter combination.

[0073] Specifically, regular expressions can be used to match the operating parameters that need to be modified in the first, second, and third macro code files. Then, based on the values ​​corresponding to the operating parameters that need to be modified, the model parameters in the first, second, and third macro code files are modified to determine the temperature field and fluid velocity field results of crystal growth under that combination of operating parameters. For example, taking the modification of model parameters in the first macro code file as an example, the model parameters that need to be modified can be matched using regular expressions. If the value of the bottom soft felt thickness needs to be modified, then the value of the moving edge in the recorded moving soft felt macro code is modified to achieve the modification of the model parameters.

[0074] Step S1032: Execute the modified first macro code file, the modified second macro code file, and the modified third macro code file in a preset order to obtain the simulation results corresponding to each combination of operating parameters. The simulation results include the temperature of the induction furnace and the fluid flow rate.

[0075] Specifically, after modifying the first, second, and third macro code files, the modified first, second, and third macro code files are executed in a preset order. It is important to note that these modified macro code files are not executed sequentially; their execution order can be determined based on the solution order of numerical values ​​such as temperature and fluid velocity in a typical crystal growth simulation.

[0076] When executing the modified first, second, and third macro code files, it is necessary to determine whether the temperature and velocity fields of the induction furnace in the target model have reached a steady state to ensure the reliability of the simulation results, such as the temperature of the induction furnace and the fluid flow velocity. Therefore, during the calculation of the temperature of the induction furnace and the fluid flow velocity, it is necessary to continuously monitor the changes in the calculated temperature and fluid flow velocity to facilitate the determination that the temperature and velocity fields of the induction furnace have reached a steady state and obtain accurate simulation results. The simulation results include the temperature of the induction furnace and the fluid flow velocity. The temperature of the induction furnace represents the crystal growth temperature, and the fluid flow velocity in the induction furnace represents the crystal growth rate. These two results are helpful for tasks such as parameter sensitivity analysis, optimization design, and machine learning-driven modeling of crystal growth.

[0077] The automated numerical simulation method for induction heating crystal growth provided in this invention uses macro code files to modify the model parameters of the target model, as well as the mesh file, power density, induction furnace temperature, and fluid flow velocity. This enables an automated numerical simulation process. Furthermore, when modifying the numerical parameters of the operating conditions, the model parameters of the first, second, and third macro code files can be directly modified, and the modified macro code files can be executed to obtain the numerical simulation results corresponding to each combination of operating parameters, thereby improving the efficiency of numerical simulation.

[0078] In one embodiment, the specific execution process of the modified first macro code file, the modified second macro code file, and the modified third macro code file includes: First, executing the modified first macro code file to obtain an updated mesh file, and based on the updated mesh file, executing the modified third macro code file to obtain the node coordinate information of the target model; then, loading the node coordinate information into the modified second macro code file and executing the modified second macro code file to obtain the node power density corresponding to each node of the induction furnace in the target model; subsequently, loading multiple node power densities into the modified third macro code file and executing the modified third macro code file to obtain the temperature and fluid flow velocity corresponding to each node of the induction furnace in the target model.

[0079] In this embodiment, the second macro code file in Maxwell models the target model and calculates the power density. Since Maxwell performs power density calculations after simplifying the target model, it's necessary to import the node coordinate information of the target model into the second macro code file. This allows the calculated node power density to be used in subsequent heat transfer and fluid flow velocity calculations. Therefore, before executing the modified second macro code file, the updated mesh file needs to be loaded into the third macro code file to obtain the node coordinate information of the target model. This node coordinate information can consist of multiple node coordinates, each representing a minimum unit in the mesh file. After loading the node coordinate information, executing the modified second macro code file determines the power density corresponding to each node coordinate of the induction furnace in the target model.

[0080] It is important to note that when importing the updated mesh file into the modified third macro code file and executing the modified third macro code file, the modified third macro code file does not complete execution due to the lack of power density parameters. It is necessary to export the obtained node coordinate information, load the node coordinate information into the modified second macro code file, obtain the node power density corresponding to each node of the induction furnace in the target model, and then continue to execute the modified third macro code file based on multiple node power densities to obtain the temperature and fluid flow velocity corresponding to each node of the induction furnace in the target model.

[0081] This can be achieved by writing a Python program to couple the acquisition of mesh files, the calculation of power density, and the calculation of temperature and fluid flow fields. The necessary conditions for these processes can be set, such as setting the task name, software location, project file location, and defining a file list. Setting the task name facilitates task management during execution, and setting the software and project file locations facilitates file retrieval during calculation. For example, setting the installation location of the software that acquires mesh files allows for automatic retrieval of the software. The file list defines the necessary filenames and locations for easy retrieval, such as defining the name and location of the mesh files, allowing for direct retrieval from the list when needed.

[0082] It is important to note that while power parameters are used to limit the intensity of the AC electric field generated by the coil during actual crystal growth, the intensity of the AC electric field generated by the coil needs to be calculated based on the induced current during numerical simulation. Therefore, when the power parameters included in the target model are modified, it is necessary to first calculate the induced current corresponding to the power parameters and determine whether the actual power value calculated based on the induced current is consistent with the power parameter value, so as to ensure the accuracy of the power density value obtained from the numerical simulation. Specifically, this process may include: first, determining the induced current based on the power parameters; then, based on the induced current and node coordinate information, executing the modified second macro code file to obtain the actual power density corresponding to each node of the induction furnace in the target model, and determining the actual power value based on multiple actual power densities; then, matching the actual power value with the power parameter value; if the actual power value and the power parameter value do not match, correcting the induced current and re-determining the actual power value based on the difference between the actual power value and the power parameter value; when the actual power value and the power parameter value match, determining the actual power density corresponding to the induced current as the node power density corresponding to each node of the induction furnace in the target model.

[0083] In this embodiment, the power density is actually calculated based on the induced current. Therefore, when modifying the power parameters, it is necessary to first calculate the induced current according to the preset empirical formula, and then use the induced current to obtain the power density of each node in the induction furnace. In order to ensure the accuracy of the obtained power density, after calculating the induced current, the node power density of each node in the induction furnace can be obtained based on the induced current. Then, the power densities of components such as the graphite heating cylinder, crucible, soft felt, induction coil, and porous graphite plate are integrated and summed to obtain the actual power value. The actual power value is then matched with the value of the power parameters to determine whether the induced current is accurate.

[0084] Specifically, the error between the actual power value and the power parameter value can be determined. If the error value is less than or equal to the set threshold, it indicates that the induced current is relatively accurate and the calculated power density can reflect the actual power density. In this case, multiple actual power densities can be exported as the node power density corresponding to each node of the induction furnace in the target model. However, if the error value is greater than the set threshold, it indicates that the induced current error is large. The previously obtained induced current can be corrected according to the magnitude of the error value, and the process of determining the actual power value and matching can be repeated until the actual power value obtained by the induced current can match the power parameter value.

[0085] In one embodiment, when simulating the crystal growth environment using a target model, the induction furnace in the target model includes a process where temperature and fluid flow velocity continuously change until a steady state is reached. However, when obtaining simulation results, the temperature and fluid flow velocity values ​​of the induction furnace obtained during the temperature and fluid flow velocity change phases do not reflect the crystal growth environment. Therefore, the temperature and fluid flow velocity in the induction furnace can be continuously monitored until the growth environment reaches a steady state. At this point, the temperature and fluid flow velocity in the induction furnace can be used as simulation results. This process specifically includes: first, obtaining the first preset node coordinates, which are the node coordinates after the target model's mesh file is updated; then, when executing the modified third macro code file, obtaining the temperature and fluid velocity corresponding to the first preset node coordinates at preset time intervals, and determining whether the temperature field and fluid flow field of the target model have both reached a steady state. Specifically, determining whether the temperature field and fluid flow field of the target model have reached a steady state is crucial. The steps for determining whether the temperature field and fluid flow field of the target model have both reached a steady state include: determining a first difference between the temperature obtained at the current time and the temperature obtained at the previous time, and a second difference between the fluid velocity obtained at the current time and the fluid velocity obtained at the previous time; then, determining whether the temperature field and fluid flow field of the target model have both reached a steady state based on the first and second differences; if the temperature field and fluid flow field of the target model have not reached a steady state, continuing to execute the steps of obtaining the temperature and fluid velocity corresponding to the first preset node coordinates at preset time intervals, and determining whether the temperature field and fluid flow field of the target model have both reached a steady state; if the temperature field and fluid flow field of the target model have both reached a steady state, determining the temperature and fluid flow corresponding to the current time as the temperature and fluid flow velocity corresponding to each node of the induction furnace after the target model mesh file is updated.

[0086] In this embodiment, any node in the induction furnace can be used as a monitoring point to continuously acquire the temperature and fluid flow velocity of the monitoring point in order to determine whether the temperature field and fluid flow field in the induction furnace remain stable. Once it is determined that the temperature field and fluid flow field of the induction furnace have reached a steady state, the temperature and fluid flow velocity under the steady state are used as the simulation results.

[0087] Specifically, the temperature and fluid flow velocity at the monitoring point can be continuously acquired at preset time intervals. These time intervals can be determined based on actual conditions, such as 1 second, 2 seconds, or 5 seconds; this embodiment is not limited to these intervals. Each time the temperature and fluid flow velocity at the monitoring point are acquired, the acquired temperature and fluid flow velocity can be compared with those acquired previously to determine whether the temperature field and fluid flow field of the target model have reached a steady state.

[0088] Understandably, when the temperature field and fluid flow field of the target model reach a steady state, the temperature and fluid flow velocity of the target model no longer change or change only slightly. Then, the magnitude of the first difference between the obtained fluid flow velocity and the fluid flow velocity obtained at the previous time, and the magnitude of the second difference between the obtained fluid flow velocity and the fluid flow velocity obtained at the previous time can be determined. Based on the magnitude of the first difference and the second difference, it can be determined whether the temperature field and fluid flow field of the target model have reached a steady state.

[0089] Specifically, the magnitude of the first difference can be compared with the magnitude of the preset temperature threshold, and the magnitude of the second difference can be compared with the magnitude of the preset speed threshold. The preset temperature threshold and the preset speed threshold can be determined according to the actual situation, and this embodiment does not limit them.

[0090] Specifically, if the first difference is less than a preset temperature threshold and the second difference is less than a preset velocity threshold, it is determined that the temperature field and fluid flow field of the target model have reached a steady state, and the temperature and fluid flow velocity of each node in the induction furnace at this time can be used as the simulation result; however, if the first difference is greater than or equal to the preset temperature threshold, or the second difference is greater than or equal to the preset velocity threshold, it indicates that the temperature field and / or fluid flow field of the target model have not reached a steady state and cannot be used as the final simulation result. In this case, the temperature and fluid flow velocity of the first preset node are monitored until it is determined that the temperature field and fluid flow field of the target model have reached a steady state.

[0091] In one embodiment, considering that different locations within the induction furnace have varying degrees of influence on crystal growth, the temperature and fluid flow velocity of nodes that significantly affect crystal growth can be obtained as simulation results to facilitate subsequent analysis and research. In this case, at least one second preset node coordinate can be obtained first, which are the node coordinates that affect crystal growth in the target model. Then, the temperature and fluid flow velocity corresponding to the second preset node coordinates are determined as the simulation results corresponding to the combination of operating parameters, and a target list is obtained. The target list includes multiple combinations of operating parameters and the simulation results corresponding to each combination of operating parameters.

[0092] In this embodiment, multiple second preset node coordinates can be pre-determined based on experience. These second preset node coordinates can be coordinates of nodes at locations that significantly influence crystal growth, such as the top hole in the furnace, the center of the seed crystal, and the edge of the seed crystal. The temperature and fluid flow velocity corresponding to these second preset node coordinates are then determined. These temperatures and fluid flow velocities are used as the simulation results for the corresponding combinations of operating parameters. This eliminates the need for manual screening, improving the efficiency of subsequent analysis of the simulation results. After obtaining the simulation results corresponding to each of the multiple combinations of operating parameters, a target list can be generated to facilitate the display of the simulation results and subsequent post-processing of the simulation results corresponding to each combination of operating parameters. The first preset node coordinates can be the same as or different from the second preset node coordinates.

[0093] In one embodiment, simulation results corresponding to multiple operating condition combinations can be obtained simultaneously based on the computing power of the device to improve the efficiency of numerical simulation. Specifically, multiple operating condition parameter combinations can be divided into multiple batches according to the number of operating condition parameter combinations. Then, for each operating condition parameter combination in the same batch, the steps of simulating the growth environment in crystal growth can be executed simultaneously.

[0094] Considering that simulating crystal growth for only one combination of operating parameters each time is inefficient, the operating parameter combinations can be divided into multiple batches according to the total number. The number of operating parameter combinations in each batch can be determined according to the computing power of the equipment used. That is, when the computing power of the equipment is high, the number of operating parameter combinations in each batch is large, while when the computing power of the equipment is low, the number of operating parameter combinations in each batch is small.

[0095] After dividing multiple combinations of operating parameters into multiple batches, the crystal growth simulation process can be performed simultaneously for each combination of operating parameters in each batch. Specifically, a separate case folder can be generated for each batch of multiple operating parameters. Then, the target model, the first macro code file, the second macro code file, the third macro code file, and other file lists can be copied to each case folder. Then, steps such as modifying the macro code file and executing the modified macro code file can be performed according to different combinations of operating parameters.

[0096] For any of the above embodiments, when modifying the model parameters of the macro code file, regular expressions can be used for matching to identify the model parameters that need to be modified. Specifically, regular expressions can be used to determine the target model parameters from the first macro code file, the second macro code file, and the third macro code file, respectively; wherein different macro code files correspond to slightly different target model parameters; then, the target parameter values ​​of the working condition parameters that match the target model parameters in the working condition parameter combinations are determined; and finally, the values ​​of the target model parameters are modified to the target parameter values.

[0097] For example, the first macro code file can be retrieved, and the target model parameters that need to be modified can be matched according to the combination of working condition parameters using regular expressions. Then, the values ​​of the target model parameters are replaced with the values ​​of the corresponding working condition parameters in the working condition parameter combination. Since the first macro code file includes macro code for modifying the working condition parameters, after replacing the values ​​of the target model parameters, executing the first macro code file will automatically modify the values ​​of the working condition parameters and draw the mesh file. The modification methods for the second macro code file are the same as those for the first macro code file, and will not be elaborated here.

[0098] The automated numerical simulation method for crystal growth during induction heating provided in this invention uses macro code files to modify the model parameters of the target model, as well as the mesh file, power density, induction furnace temperature, and fluid flow velocity. This enables an automated numerical simulation process. Furthermore, when modifying the numerical parameters of the operating conditions, the model parameters of the first, second, and third macro code files can be directly modified, and the modified macro code files can be executed to obtain the numerical simulation results corresponding to each combination of operating parameters, thereby improving the efficiency of numerical simulation.

[0099] The automated numerical simulation method for induction heating crystal growth provided by the embodiments of the present invention will be described in detail below with reference to examples:

[0100] Reference Figure 4 , Figure 4 This diagram illustrates a flowchart of an automated numerical simulation method for induction heating crystal growth provided in an embodiment of the present invention. The first macro code file is obtained based on macro recording provided by Gambit software, the second macro code file is obtained based on Maxwell software, and the third macro code file is obtained based on Fluent software. Figure 4As shown, firstly, the first macro code file, the second macro code file, and the third macro code file are obtained. The first macro code file uses Gambit software to draw the mesh file, the second macro code file uses Maxwell software to calculate the power density of the target model, and the third macro code file uses Fluent software to calculate the temperature and fluid flow velocity of the target model.

[0101] Next, multiple combinations of operating parameters are received. Before this, operating parameters can be initialized, defining various combinations. Specifically, based on important operating parameters in the induction heating crystal growth process, several operating parameters are defined as follows: power parameter, bottom felt thickness parameter, coil axial movement distance parameter, graphite cap axial distance parameter, powder height parameter, top hole size parameter, graphite conductivity, powder thermal conductivity, felt thermal conductivity, and graphite thermal conductivity. Each received combination of operating parameters then includes these 10 parameters. Specifically, the value range corresponding to each operating parameter can be received, and multiple different combinations of operating parameters can be obtained by iteratively traversing the value range of each operating parameter.

[0102] Next, multiple combinations of operating parameters are divided into multiple batches according to computing power. The crystal growth simulation is performed simultaneously on multiple combinations of operating parameters in each batch. During this process, a separate case folder is generated for each combination of operating parameters in the same batch. The basic file list is copied to multiple case folders. Then, regular expressions are used to determine the model parameters that need to be modified from the first macro code file, the second macro code file, and the third macro code file, respectively, based on the combination of operating parameters.

[0103] Next, the modified first macro code file, the modified second macro code file, and the modified third macro code file are executed in a specific order to obtain the simulation results. Specifically, refer to... Figures 4-6 , Figure 5 A schematic diagram illustrating the process of obtaining power density when power parameters are modified is shown. Figure 6The diagram illustrates the process of obtaining the temperature and fluid velocity parameters of the induction furnace in the target model. Gambit is invoked to read the baseline project file and load the modified first macro code. Gambit automatically executes the first macro code to modify the target model and redraw the mesh, obtaining the mesh file. Then, Fluent is launched and reads the baseline project file. The updated mesh file is loaded into the modified third macro code file, and the modified third macro code file is executed to obtain the node coordinate information. Finally, the node coordinate information is loaded into the second macro code file, and the modified second macro code file is executed to obtain the nodal power density of each node of the induction furnace in the target model. When the operating parameter to be modified in the second macro code file is the power parameter, since the power parameter is calculated by Maxwell and then integrated and summed, in order to calculate the required power value, the corresponding power induction current is first preset by experience. Maxwell is started and the baseline project file is read. After executing the modified second macro code file to calculate the node power density of multiple nodes, the actual power value is obtained by integrating and summing the power densities of graphite parts, soft felt and stainless steel parts. The actual power value is then compared with the value of the power parameter to determine whether the error meets the set threshold. If the condition is met, the node power density corresponding to each node in the induction furnace is exported. If the condition is not met, the induction current used in the previous calculation is corrected according to the error value.

[0104] Then, the node power density corresponding to each node can be loaded into the third macro code file. When the modified third macro code file is executed, the temperature and fluid flow velocity corresponding to each node can be calculated based on the power density of multiple nodes, thereby obtaining the temperature and fluid flow velocity of each node of the induction furnace in the target model.

[0105] At this point, the coordinates of the first preset node can be obtained from the basic engineering file to continuously monitor the temperature and fluid flow velocity of the first preset node coordinates. Specifically, the temperature and fluid flow velocity of the first preset node coordinates can be acquired every 1 second, and a first difference between the temperature of the currently acquired first preset node coordinates and the temperature of the first preset node coordinates acquired at the previous time can be determined, as well as a second difference between the fluid flow velocity of the currently acquired first preset node coordinates and the fluid flow velocity of the first preset node coordinates acquired at the previous time can be determined. Then, the magnitude of the first difference and a preset temperature threshold, and the magnitude of the second difference and a preset velocity threshold can be compared. If the first difference is less than the preset temperature threshold and the second difference is less than the preset velocity threshold, it indicates that the temperature and velocity fields of the induction furnace in the target model have reached a steady state, and the temperature and velocity field data obtained at this time can be used as the simulation result. However, if the first difference is greater than or equal to the preset temperature threshold or the second difference is greater than or equal to the preset velocity threshold, it indicates that the temperature and velocity fields of the induction furnace in the target model have not yet reached a steady state, and it is necessary to continue monitoring the temperature and velocity of the first preset node coordinates and comparing the magnitude of the first difference and the preset temperature threshold, and the magnitude of the second difference and the preset velocity threshold.

[0106] To facilitate post-processing of the data, the basic engineering file also includes the coordinate information of the locations to be monitored, namely the second preset node coordinates. These can be locations that have a significant impact on crystal growth, such as the top hole in the furnace, the center of the seed crystal, and the edge of the seed crystal. The temperature and fluid flow velocity corresponding to the second preset node coordinates can be used as the simulation results corresponding to the combination of operating parameters. Then, a target list can be generated based on the simulation results corresponding to each combination of operating parameters. Subsequently, tasks such as parameter sensitivity analysis, optimization design, and machine learning-driven modeling can be performed directly based on the target list, thereby improving the efficiency of data post-processing.

[0107] Based on the same inventive concept, embodiments of the present invention also provide an automated numerical simulation device for induction heating crystal growth, referring to... Figure 7 , Figure 7 A schematic diagram of the structure of the automated numerical simulation device for induction heating crystal growth provided in an embodiment of the present invention is shown, as follows: Figure 7 As shown, the device specifically includes:

[0108] The acquisition module 201 is used to acquire the first macro code file, the second macro code file, and the third macro code file of the target model; wherein, the first macro code file is used to acquire the mesh file of the target model, the second macro code file is used to acquire the power density of the induction furnace in the target model, and the third macro code file is used to acquire the temperature and fluid flow velocity of the induction furnace in the target model.

[0109] The receiving module 202 is used to receive multiple combinations of operating parameters, where at least one operating parameter has a different value range among the different combinations of operating parameters.

[0110] The simulation module 203 is used to perform the steps of simulating the growth environment in crystal growth based on the target model for each combination of operating parameters. The steps include: modifying the model parameters in the modified first macro code file, the modified second macro code file, and the modified third macro code file respectively based on each combination of operating parameters; and executing the modified first macro code file, the modified second macro code file, and the modified third macro code file in a preset order to obtain the simulation results corresponding to each combination of operating parameters. The simulation results include the temperature of the induction furnace and the fluid flow rate.

[0111] Based on the same inventive concept, embodiments of the present invention also provide an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the steps of the automated numerical simulation method for induction heating crystal growth as described in any of the above embodiments.

[0112] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.

[0113] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0114] The above provides a detailed description of an automated numerical simulation method, apparatus, and equipment for induction heating crystal growth provided by the present invention. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

[0115] Other embodiments of the invention will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This invention is intended to cover any variations, uses, or adaptations of the invention that follow the general principles of the invention and include common knowledge or customary techniques in the art not disclosed herein. The specification and examples are to be considered exemplary only, and the true scope and spirit of the invention are indicated by the following claims.

[0116] It should be understood that the present invention is not limited to the precise structure described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of the invention is limited only by the appended claims.

[0117] The terms "an embodiment," "embodiment," or "one or more embodiments" as used herein mean that a particular feature, structure, or characteristic described in connection with an embodiment is included in at least one embodiment of the invention. Furthermore, please note that the examples of the phrase "in one embodiment" do not necessarily all refer to the same embodiment.

[0118] Numerous specific details are set forth in the specification provided herein. However, it will be understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures, and techniques have not been shown in detail so as not to obscure the understanding of this specification.

[0119] In the claims, any reference signs placed between parentheses should not be construed as limiting the claims. The word "comprising" does not exclude the presence of elements or steps not listed in the claims. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention can be implemented by means of hardware comprising a plurality of different elements and by means of a suitably programmed computer. In a unit claim enumerating a plurality of means, several of these means may be embodied by the same item of hardware. The use of the words first, second, and third, etc., does not indicate any order. These words may be interpreted as names.

[0120] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. An automated numerical simulation method for induction heating crystal growth, characterized in that, The method includes: Obtain a first macro code file, a second macro code file, and a third macro code file of the target model; wherein, the first macro code file is used to obtain the mesh file of the target model, the second macro code file is used to obtain the power density of the induction furnace in the target model, and the third macro code file is used to obtain the temperature and fluid flow velocity of the induction furnace in the target model; It can receive multiple combinations of operating parameters, and at least one of the operating parameters has a different value range among the different combinations of operating parameters; For each combination of operating parameters, a step of simulating the growth environment in crystal growth is performed based on the target model. The step includes: modifying the model parameters in the first macro code file, the second macro code file, and the third macro code file respectively based on each combination of operating parameters. The modified first macro code file, the modified second macro code file, and the modified third macro code file are executed sequentially according to a preset order to obtain simulation results corresponding to each combination of operating parameters. The simulation results include the temperature of the induction furnace and the fluid flow velocity. The step of executing the modified first macro code file, the modified second macro code file, and the modified third macro code file sequentially according to a preset order includes: The modified first macro code file is executed to obtain the updated mesh file, and based on the updated mesh file, the modified third macro code file is executed to obtain the node coordinate information of the target model; The node coordinate information is loaded into the modified second macro code file, and the modified second macro code file is executed to obtain the node power density corresponding to each node of the induction furnace in the target model. The power density of multiple nodes is loaded into the modified third macro code file, and the modified third macro code file is executed to obtain the temperature and fluid flow velocity corresponding to each node of the induction furnace in the target model.

2. The automated numerical simulation method for induction heating crystal growth according to claim 1, characterized in that, When the power parameters included in the target model are modified, the step of executing the modified second macro code file based on the node coordinate information to obtain the node power density corresponding to each node of the induction furnace in the target model includes: Determine the induced current based on the power parameters; Based on the induced current and the node coordinate information, the modified second macro code file is executed to obtain the actual power density corresponding to each node of the induction furnace in the target model, and the actual power value is determined according to the multiple actual power densities. Match the actual power value with the value of the power parameter; When the actual power value does not match the value of the power parameter, the induced current is corrected and the actual power value is re-determined based on the difference between the actual power value and the power parameter. When the actual power value matches the value of the power parameter, the actual power density corresponding to the induced current is determined as the node power density corresponding to each node of the induction furnace in the target model.

3. The automated numerical simulation method for induction heating crystal growth according to claim 1, characterized in that, The step of loading the power density of multiple nodes into the modified third macro code file and executing the modified third macro code file to obtain the temperature and fluid flow velocity corresponding to each node of the induction furnace in the target model includes: Obtain the first preset node coordinates, which are the node coordinates of the target model after updating the mesh file; When executing the modified third macro code file, the temperature and fluid velocity corresponding to the first preset node coordinates are obtained at preset time intervals, and it is determined whether the temperature field and fluid flow field of the target model have both reached a steady state; the determination of whether the temperature field and fluid flow field of the target model have both reached a steady state includes: determining the first difference between the temperature obtained at the current time and the temperature obtained at the previous time, and the second difference between the fluid velocity obtained at the current time and the fluid velocity obtained at the previous time. Based on the first difference and the second difference, determine whether the temperature field and fluid flow field of the target model have both reached a steady state; If the temperature field or fluid flow field of the target model has not reached a steady state, continue to execute the steps of obtaining the temperature and fluid velocity corresponding to the first preset node coordinates at preset time intervals, and determining whether the temperature field and fluid flow field of the target model have both reached a steady state. When the temperature field and fluid flow field of the target model reach a steady state, the temperature and fluid flow velocity corresponding to the current time are determined as the temperature and fluid flow velocity corresponding to each node of the induction furnace after the target model updates the mesh file.

4. The automated numerical simulation method for induction heating crystal growth according to claim 3, characterized in that, The step of determining whether the temperature field and fluid flow field of the target model have both reached a steady state based on the first difference and the second difference includes: Compare the first difference with a preset temperature threshold, and the second difference with a preset speed threshold; If the first difference is less than the preset temperature threshold and the second difference is less than the preset velocity threshold, it is determined that the temperature field and fluid flow field of the target model have both reached a steady state.

5. The automated numerical simulation method for induction heating crystal growth according to claim 4, characterized in that, The method further includes: Obtain at least one second preset node coordinate, wherein the second preset node coordinate is the node coordinate that affects crystal growth in the target model; The temperature and fluid flow velocity corresponding to the second preset node coordinates are determined as the simulation results corresponding to the combination of operating parameters, and a target list is obtained. The target list includes multiple combinations of operating parameters and the simulation results corresponding to each combination of operating parameters.

6. The automated numerical simulation method according to claim 1, characterized in that, The method further includes: Based on the number of the operating condition parameter combinations, the multiple operating condition parameter combinations are divided into multiple batches; For each combination of operating parameters in the same batch, the step of simulating the growth environment in crystal growth is performed simultaneously.

7. The automated numerical simulation method for induction heating crystal growth according to claim 1, characterized in that, The step of modifying the model parameters in the first macro code file, the second macro code file, and the third macro code file based on each combination of operating condition parameters includes: Regular expressions are used to determine the target model parameters from the first macro code file, the second macro code file, and the third macro code file, respectively; wherein different macro code files correspond to slightly different target model parameters. Determine the target parameter value of the working condition parameter that matches the target model parameter in the working condition parameter combination; Modify the numerical values ​​of the target model parameters to the target parameter values.

8. An automated numerical simulation device for induction heating crystal growth, characterized in that, The device includes: The acquisition module is used to acquire a first macro code file, a second macro code file, and a third macro code file of the target model; wherein, the first macro code file is used to acquire the mesh file of the target model, the second macro code file is used to acquire the power density of the induction furnace in the target model, and the third macro code file is used to acquire the temperature and fluid flow velocity of the induction furnace in the target model; The receiving module is used to receive multiple combinations of operating parameters, where at least one operating parameter has a different value range among the different combinations of operating parameters. The simulation module is used to perform a step of simulating the growth environment in crystal growth based on the target model for each combination of operating parameters. The step includes: modifying the model parameters in the first macro code file, the second macro code file, and the third macro code file respectively for each combination of operating parameters; and executing the modified first macro code file, the modified second macro code file, and the modified third macro code file sequentially in a preset order to obtain simulation results corresponding to each combination of operating parameters. The simulation results include the temperature and fluid flow rate of the induction furnace. The step of executing the modified first macro code file, the modified second macro code file, and the modified third macro code file sequentially in a preset order includes: The modified first macro code file is executed to obtain the updated mesh file, and based on the updated mesh file, the modified third macro code file is executed to obtain the node coordinate information of the target model; The node coordinate information is loaded into the modified second macro code file, and the modified second macro code file is executed to obtain the node power density corresponding to each node of the induction furnace in the target model. The power density of multiple nodes is loaded into the modified third macro code file, and the modified third macro code file is executed to obtain the temperature and fluid flow velocity corresponding to each node of the induction furnace in the target model.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes, it implements the steps in the automated numerical simulation method for induction heating crystal growth as described in any one of claims 1-7.

Images