Molten pool characteristic analysis method, device, computer equipment and storage medium
By obtaining the natural and forced convection intensity values of the molten pool and determining the convection influencing factors and characteristics, the problem of the lack of molten pool characteristic analysis under the swinging state is solved, ensuring that the molten material is retained in the pressure vessel and reducing the risk of radioactive material release.
Patent Information
- Application Number
- CN202310143458.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-02-14
- Publication Date
- 2025-09-19
- Estimated Expiration
- 2043-02-14
AI Technical Summary
The lack of effective melt pool characteristic analysis methods, especially those for melt pool characteristics in a swinging state, affects the integrity of pressure vessels and the risk of radioactive material release.
By obtaining the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state, the convection influencing factor is determined. The presence of forced convection in a swinging state is judged based on the preset threshold, and the convection characteristics of the molten pool are determined in combination with the free liquid surface height and acceleration.
It provides an accurate method for analyzing the characteristics of the molten pool, especially the characteristics analysis under the swing state, to help determine whether the molten material is retained in the pressure vessel and reduce the risk of radioactive material release.
Smart Images

Figure CN116202730B_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the field of computer technology, and in particular to a method, apparatus, computer equipment, and storage medium for analyzing molten pool characteristics. Background Art
[0002] In the event of a severe accident at a conventional pressurized water reactor nuclear power plant, if the reactor core cannot be adequately cooled, the fuel assemblies and structural materials will melt under the continued heating of decay heat, gradually migrating to the bottom of the pressure vessel's lower head, forming a high-temperature, highly radioactive molten pool. Because the molten pool continuously transfers decay heat to the lower head wall, if the decay heat cannot be fully dissipated, the high-temperature molten pool may melt through the lower head wall, seriously threatening the integrity of the pressure vessel.
[0003] Therefore, it is necessary to find a way to retain the molten material in the lower head of the pressure vessel to avoid subsequent accidents and prevent the release of radioactive materials into the external environment. The success of the melt retention technology is closely related to the characteristics of the molten pool (such as convection characteristics). However, there is currently a lack of effective methods for analyzing the characteristics of the molten pool, especially the lack of methods for analyzing the characteristics of the molten pool in a swinging state (for example, the molten pool formed when an accident occurs in a small offshore nuclear reactor). Summary of the Invention
[0004] Based on this, it is necessary to provide a method, device, computer equipment and storage medium for analyzing the characteristics of the molten pool under a swinging state to address the above technical problems.
[0005] In a first aspect, the present application provides a method for analyzing melt pool characteristics. The method comprises:
[0006] Obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0007] Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0008] According to the relationship between the convection influencing factor and the preset threshold, it is determined whether there is forced convection in a swaying state inside the molten pool;
[0009] If present, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state.
[0010] In one embodiment, obtaining the natural convection intensity value of the molten pool in a static state includes:
[0011] The natural convection intensity value of the molten pool in a static state is determined based on the first fluid parameters of the molten pool; wherein the first fluid parameters include the thermal expansion coefficient, the difference between the maximum temperature in the molten pool and the phase interface temperature, the height of the molten pool and the kinematic viscosity of the melt.
[0012] In one embodiment, obtaining a forced convection intensity value of the molten pool in a swaying state includes:
[0013] Determining a fluid momentum equation based on second fluid parameters of the molten pool; wherein the second fluid parameters include a velocity of the molten pool, a fluid density, a pressure, a kinematic viscosity of the melt, and a rocking angular velocity;
[0014] The fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation;
[0015] According to the dimensionless momentum equation, the forced convection intensity value of the molten pool in the swing state is determined.
[0016] In one embodiment, the fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation, including:
[0017] Determine dimensionless parameters based on characteristic parameters and third fluid parameters of the molten pool; wherein the characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameters include the radius vector, motion velocity, time, fluid density, sway angular velocity, pressure, and gravitational acceleration of the molten pool; and the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration;
[0018] According to the characteristic parameters and dimensionless parameters, the fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation.
[0019] In one embodiment, determining the convection influencing factor of the molten pool according to the natural convection intensity value and the forced convection intensity value includes:
[0020] The forced convection intensity value is squared and multiplied by the natural convection intensity value to obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt.
[0021] In one embodiment, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state, including:
[0022] Constructing a partial differential equation of the relative velocity potential function of the free liquid surface and boundary conditions of the partial differential equation based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state;
[0023] Based on the boundary conditions, the partial differential equations are solved to obtain the convection characteristics of the molten pool.
[0024] In one embodiment, the method further includes:
[0025] Obtain the Nusselt number of the heat transfer in the molten pool in a static state and the Nusselt number of the heat transfer in a swinging state;
[0026] The heat transfer characteristics of the molten pool are determined based on the heat transfer Nusselt number of the molten pool in a static state, the heat transfer Nusselt number of the molten pool in a swinging state, and the convection influence factor of the molten pool.
[0027] In a second aspect, the present application also provides a molten pool characteristic analysis device. The device comprises:
[0028] An acquisition module is used to obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0029] A first determination module is used to determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0030] A second determination module is used to determine whether there is forced convection in a swaying state inside the molten pool based on the relationship between the convection influencing factor and a preset threshold;
[0031] The third determining module, if present, determines the convection characteristics of the molten pool based on the height and lateral span of the free liquid surface in the molten pool.
[0032] In a third aspect, the present application further provides a computer device. The computer device includes a memory and a processor. The memory stores a computer program. When the processor executes the computer program, the following steps are performed:
[0033] Obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0034] Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0035] According to the relationship between the convection influencing factor and the preset threshold, it is determined whether there is forced convection in a swaying state inside the molten pool;
[0036] If present, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool.
[0037] In a fourth aspect, the present application further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the following steps:
[0038] Obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0039] Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0040] According to the relationship between the convection influencing factor and the preset threshold, it is determined whether there is forced convection in a swaying state inside the molten pool;
[0041] If present, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool.
[0042] In a fifth aspect, the present application further provides a computer program product. The computer program product includes a computer program that, when executed by a processor, implements the following steps:
[0043] Obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0044] Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0045] According to the relationship between the convection influencing factor and the preset threshold, it is determined whether there is forced convection in a swaying state inside the molten pool;
[0046] If present, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool.
[0047] The above-mentioned molten pool characteristic analysis method, device, computer equipment and storage medium, the method of the present application obtains the natural convection intensity value of the molten pool in a static state, and the forced convection intensity value of the molten pool in a swinging state, and then determines the convection influence factor of the molten pool based on the natural convection intensity value and the forced convection intensity value, and finally determines whether there is forced convection in a swinging state inside the molten pool based on the relationship between the convection influence factor and the preset threshold; if so, then the convection characteristics of the molten pool in the swinging state can be determined based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in the swinging state, which solves the current problem of lack of effective methods for molten pool characteristic analysis, especially the lack of methods for determining the molten pool characteristics in a swinging state. BRIEF DESCRIPTION OF THE DRAWINGS
[0048] Figure 1 A diagram illustrating the application environment of the molten pool characteristic analysis method provided in this embodiment;
[0049] Figure 2 A schematic flow chart of the first molten pool characteristic analysis method provided in this embodiment;
[0050] Figure 3 A schematic diagram of a process for determining the forced convection intensity value of a molten pool in a swinging state provided in this embodiment;
[0051] Figure 4 A schematic diagram of a process for determining the convection characteristics of a molten pool provided in this embodiment;
[0052] Figure 5 A schematic flow chart of the second molten pool characteristic analysis method provided in this embodiment;
[0053] Figure 6 A schematic flow chart of a third method for analyzing molten pool characteristics provided in this embodiment;
[0054] Figure 7 A structural block diagram of the first molten pool characteristic analysis device provided in this embodiment;
[0055] Figure 8 A structural block diagram of the second molten pool characteristic analysis device provided in this embodiment;
[0056] Figure 9 A structural block diagram of the third molten pool characteristic analysis device provided in this embodiment;
[0057] Figure 10 This is a diagram of the internal structure of the computer device provided in this embodiment. DETAILED DESCRIPTION
[0058] In order to make the purpose, technical solutions and advantages of this application more clear, the following further describes this application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not intended to limit this application.
[0059] In one embodiment, a computer device is provided. The computer device may be a server, and its internal structure diagram may be as follows: Figure 1As shown. The computer device includes a processor, a memory, and a network interface connected via a system bus. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store acquisition data of abnormal data. The network interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, a method for analyzing molten pool characteristics is implemented.
[0060] In one embodiment, a method for analyzing melt pool characteristics is provided, such as Figure 2 As shown, the following steps are included:
[0061] S201, obtaining a natural convection intensity value of the molten pool in a static state and a forced convection intensity value in a swinging state.
[0062] The natural convection intensity value characterizes the strength of natural convection in a stationary melt pool and can be expressed as the Grashof number. Natural convection refers to convection in a stationary melt pool driven by decaying internal heat sources. The forced convection intensity value characterizes the strength of forced convection caused by the swaying of the melt pool and can be expressed as the Ekman number. Forced convection refers to convection caused by the swaying of the melt pool.
[0063] In this embodiment, an optional implementation method for obtaining the natural convection intensity value of the molten pool in a static state is:
[0064] The natural convection intensity value of the molten pool in a static state is determined according to the first fluid parameter of the molten pool.
[0065] Among them, the first fluid parameters are used to determine the parameters required for the natural convection intensity value of the molten pool in a static state, mainly including the thermal expansion coefficient of the molten pool, the difference between the maximum temperature in the molten pool and the phase interface temperature, the height of the molten pool and the kinematic viscosity of the melt.
[0066] Specifically, in this embodiment, the Grashof number can be used to characterize the natural convection intensity value of the molten pool in a static state. The Grashof number can be determined by the following formula (1):
[0067]
[0068] In formula (1), g represents the acceleration due to gravity; β represents the thermal expansion coefficient; ΔT represents the difference between the maximum temperature in the molten pool and the phase interface temperature; h represents the height of the molten pool; and v represents the kinematic viscosity of the melt.
[0069] In this embodiment, an optional implementation method for obtaining the forced convection intensity value of the molten pool in the swinging state is: based on the fluid parameters of the molten pool and the swing parameters in the swinging state, combined with the fluid momentum equation, the forced convection intensity value of the molten pool in the swinging state can be obtained.
[0070] S202: Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value.
[0071] The convection influence factor refers to the correlation factor used to correlate the natural convection intensity value and the forced convection intensity value, and is used to determine whether the molten pool is dominated by natural convection or forced convection.
[0072] An optional implementation of this embodiment is: the natural convection intensity value and the forced convection intensity value can be input into a trained neural network model, and the neural network model determines the convection influence factor of the molten pool based on the natural convection intensity value and the forced convection intensity value.
[0073] Another optional implementation of this embodiment is: corresponding weights can be assigned to the natural convection intensity value and the forced convection intensity value respectively, and the product of the natural convection intensity value and the corresponding weight is calculated to obtain the natural convection intensity influence value, and the product of the forced convection intensity value and the corresponding weight is calculated to obtain the forced convection intensity influence value. Finally, the convection influence factor of the molten pool is determined by taking the sum of the natural convection intensity influence value and the forced convection intensity influence value.
[0074] Another optional implementation of this embodiment is: the convection influence factor of the molten pool can be determined by calculating the product of the natural convection intensity value and the forced convection intensity value.
[0075] S203 , determining whether there is forced convection in a swaying state inside the molten pool based on the relationship between the convection influencing factor and a preset threshold.
[0076] The preset threshold refers to a pre-set threshold used to determine whether there is forced convection in a swaying state inside the molten pool.
[0077] Optionally, in this embodiment, the convection influence factor can be compared with a preset threshold to determine a comparison result, and based on the comparison result, it can be determined whether there is a threshold for forced convection in a swinging state inside the molten pool.
[0078] For example, the preset threshold is 10, and the convection influence factor can be compared with the preset threshold 10. If the convection influence factor is greater than 10, it is determined that natural convection is dominant in the molten pool, and the existence of forced convection is not considered; if the convection influence factor is less than or equal to 10 and greater than 1, it is determined that mixed convection is dominant in the molten pool, that is, natural convection and forced convection compete with each other; if the convection influence factor is less than or equal to 1, it is determined that forced convection is dominant in the molten pool. In this embodiment, both molten pools dominated by mixed convection and molten pools dominated by forced convection are considered to have forced convection in the molten pool. That is, if the convection influence factor is less than or equal to 10, it is determined that forced convection in a swinging state exists in the molten pool.
[0079] S204: If so, determine the convection characteristics of the molten pool according to the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in the swaying state.
[0080] Among them, this embodiment approximately simplifies the swinging behavior of the molten pool into the forced swinging behavior of a single layer of fluid in a two-dimensional rectangular container. The height of the free liquid surface refers to the height between the free liquid surface of the fluid in the two-dimensional rectangular container and the bottom of the two-dimensional rectangular container in a static state, and the horizontal span refers to the width between the two sides of the two-dimensional rectangular container; the convection characteristics refer to the characteristics exhibited by the convection inside the molten pool, mainly including the sloshing frequency and sloshing waveform.
[0081] An optional implementation of this embodiment is: the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in a swinging state can be input into a trained neural network model, and the convection characteristics of the molten pool can be determined by the neural network model.
[0082] This embodiment obtains the natural convection intensity value of the molten pool in a static state and the forced convection intensity value of the molten pool in a swinging state, and then determines the convection influence factor of the molten pool based on the natural convection intensity value and the forced convection intensity value. Finally, based on the relationship between the convection influence factor and the preset threshold value, it is determined whether there is forced convection in the swinging state inside the molten pool; if so, the convection characteristics of the molten pool in the swinging state can be determined based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in the swinging state, which solves the problem of the current lack of effective methods for analyzing molten pool characteristics, especially the lack of methods for determining the characteristics of the molten pool in the swinging state.
[0083] In one embodiment, in order to accurately obtain the forced convection intensity value of the molten pool in the swing state, such as Figure 3As shown, in S201, an optional implementation method for obtaining the forced convection intensity value of the molten pool in the swinging state includes:
[0084] S301, determining a fluid momentum equation according to a second fluid parameter of the molten pool.
[0085] Among them, the second fluid parameters refer to the relevant parameters of the molten pool fluid used to determine the fluid momentum equation, including at least the movement speed of the molten pool, fluid density, pressure, kinematic viscosity of the melt, rocking angular velocity, gravitational acceleration and the radius vector of the fluid particle relative to the rotation axis; the fluid momentum equation refers to the momentum equation determined based on the second fluid parameters of the molten pool.
[0086] In this embodiment, the swinging motion of the molten pool is regarded as a rotational motion with a variable angular velocity. From the perspective of relative motion, the fluid momentum equation in the rotational motion reference system is established. The fluid momentum equation can be expressed by the following formula (2):
[0087]
[0088] In formula (2), represents the speed of movement; ρ represents the density of the fluid; p represents the pressure; ν represents the viscosity of the fluid; represents the sway angular velocity; represents the vector gravitational acceleration; Represents the radius vector of the fluid particle relative to the rotation axis.
[0089] S302, performing dimensionless processing on the fluid momentum equation to obtain a dimensionless momentum equation.
[0090] The dimensionless momentum equation refers to the momentum equation obtained by substituting dimensionless parameters into the fluid momentum equation.
[0091] Optionally, in this embodiment, the dimensionless parameter is determined based on the characteristic parameter and the third fluid parameter of the molten pool.
[0092] Among them, characteristic parameters refer to parameters used to determine dimensionless parameters in combination with the third fluid parameters of the molten pool, and include at least characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration. The third fluid parameters refer to parameters of the molten pool selected for determining dimensionless parameters, and include at least the molten pool's radius vector, motion velocity, time, fluid density, sway angular velocity, pressure, and vector gravitational acceleration. The molten pool's radius vector refers to the radius vector of the fluid particle relative to the rotation axis. Dimensionless parameters refer to parameters determined based on the characteristic parameters and the third fluid parameters of the molten pool, and include at least dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration.
[0093] Optionally, in this embodiment, the dimensionless parameter may be determined based on the ratio of the characteristic parameter to the corresponding third fluid parameter of the molten pool.
[0094] Specifically, the expression of dimensionless parameters is as follows: dimensionless velocity dimensionless time dimensionless gravitational acceleration dimensionless density Dimensionless pressure dimensionless radius Dimensionless swing angular velocity
[0095] Among them, L0 represents the characteristic length, U0 represents the characteristic velocity, T0 represents the characteristic time, ρ0 represents the characteristic density, Ω0 represents the characteristic angular velocity, and g0 represents the characteristic gravitational acceleration; represents the radius vector of the molten pool, represents the speed of movement, t represents time, ρ represents the fluid density, represents the angular velocity of the swing, p represents the pressure and represents the vector acceleration due to gravity.
[0096] According to the characteristic parameters and dimensionless parameters, the fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation.
[0097] Optionally, in this embodiment, the characteristic parameters and dimensionless parameters are substituted into the fluid momentum equation to obtain a dimensionless momentum equation. The dimensionless momentum equation can be expressed by the following formula (3):
[0098]
[0099] S303: Determine the forced convection intensity value of the molten pool in the swinging state according to the dimensionless momentum equation.
[0100] In this embodiment, the Ekman number is determined based on the dimensionless momentum equation, and the forced convection intensity value of the molten pool in the swinging state is determined based on the Ekman number.
[0101] In this embodiment, the Ekman number can be expressed using the following formula (4):
[0102]
[0103] Where EK is the Ekman number; ν is the kinematic viscosity of the fluid; Ω is the maximum angular velocity of the swing; and L0 is the characteristic length, which is the vertical height of the free surface of the molten pool relative to the swing axis in the static state.
[0104] In this embodiment, the fluid momentum equation is determined based on the second fluid parameter of the molten pool; the fluid momentum equation is then dimensionlessly processed to obtain a dimensionless momentum equation; and finally, based on the dimensionless momentum equation, the forced convection intensity value of the molten pool in a swinging state can be accurately determined.
[0105] On the basis of the above embodiment, in order to more accurately determine the convection impact factor, an optional implementation of S202 is:
[0106] The forced convection intensity value is squared and multiplied by the natural convection intensity value to obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt. The convection influence factor can be calculated according to the following formula (5):
[0107]
[0108] Wherein, EK in formula (5) represents the Ekman number used to characterize the forced convection intensity value of the molten pool in the swinging state, Gr represents the Grashof number used to characterize the natural convection intensity value of the molten pool in the static state, g represents the acceleration of gravity; β represents the thermal expansion coefficient; ΔT represents the difference between the maximum temperature in the molten pool and the phase interface temperature; h represents the height of the molten pool; Ω represents the maximum swing angular velocity; L0 represents the characteristic length, which is taken as the vertical height of the free liquid surface of the molten pool relative to the swing axis in the static state.
[0109] In this embodiment, the forced convection intensity value is squared and the squared result is multiplied by the natural convection intensity value, so as to accurately obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt, thereby improving the accuracy of judging whether forced convection exists in the molten pool.
[0110] In one embodiment, in order to accurately determine the convection characteristics of the molten pool, such as Figure 4 As shown, an optional implementation of S204 includes:
[0111] S401: Construct a partial differential equation of the relative velocity potential function of the free surface and boundary conditions of the partial differential equation based on the height and lateral span of the free surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in a swaying state.
[0112] Among them, the partial differential equation of the relative velocity potential function refers to the partial differential equation related to the free liquid surface constructed based on the height and lateral span of the free liquid surface in the molten pool; the boundary condition refers to the constraint condition used to solve the partial differential equation.
[0113] Optionally, in this embodiment, based on the height and lateral span of the free liquid surface in the molten pool, combined with the non-inertial coordinate system (with the lateral span as the x-axis and the height of the free liquid surface as the y-axis), the partial differential equation of the relative velocity potential function of the free liquid surface and the boundary conditions of the partial differential equation are constructed.
[0114] The partial differential equation of the relative velocity potential function of the free liquid surface in this embodiment can be expressed by the following formula (6):
[0115]
[0116] In this embodiment, the boundary conditions of the partial differential equation can be expressed by the following formula (7):
[0117]
[0118]
[0119]
[0120]
[0121] Among them, Φ in formula (6) and formula (7) r (x, y, t) represents the relative velocity potential function of liquid sloshing; h(x, t) represents the wave height of the free liquid surface relative to the initial liquid surface; It represents the horizontal acceleration caused by the swaying motion; It represents the vertical acceleration caused by the rocking motion; g represents the acceleration due to gravity.
[0122] S402, based on the boundary conditions, solving the partial differential equation to obtain the convection characteristics of the molten pool.
[0123] Optionally, in this embodiment, the boundary conditions and partial differential equations may be input into a solution model, and the solution model solves the partial differential equations based on the boundary conditions to obtain the convection characteristics of the molten pool. The convection characteristics of the molten pool include sloshing frequency and sloshing waveform.
[0124] This embodiment constructs a partial differential equation of the relative velocity potential function of the free liquid surface and the boundary conditions of the partial differential equation based on the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal and vertical accelerations of the molten pool in a swinging state. Based on the boundary conditions, the partial differential equation is solved to accurately obtain the convection characteristics of the molten pool. The determined convection characteristics of the molten pool provide a basis for the subsequent structural design of the lower head.
[0125] In one embodiment, to determine the heat transfer characteristics of the molten pool, such as Figure 5As shown, an optional embodiment of a melt pool characteristic analysis method includes:
[0126] S501, obtaining a natural convection intensity value of the molten pool in a static state and a forced convection intensity value in a swinging state.
[0127] S502: Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value.
[0128] S503 , obtaining the heat transfer Nusselt number of the molten pool in a static state and the heat transfer Nusselt number of the molten pool in a swinging state.
[0129] Among them, the heat transfer Nusselt number refers to the criterion number of the intensity of convective heat transfer, and also represents the ratio of the thermal conductivity resistance of the bottom layer of the fluid laminar flow to the convective heat transfer resistance.
[0130] Optionally, in this embodiment, the structure and operating conditions of the molten pool may be determined first, and based on the structure and operating conditions of the molten pool, the heat exchange Nusselt number of the molten pool in a stationary state and the heat exchange Nusselt number in a swinging state may be determined.
[0131] S504 , determining the heat transfer characteristics of the molten pool according to the heat transfer Nusselt number of the molten pool in a stationary state and the heat transfer Nusselt number of the molten pool in a swinging state, as well as the convection influence factor of the molten pool.
[0132] Optionally, in this embodiment, the general heat transfer expression of the molten pool can be solved based on the heat transfer Nusselt number of the molten pool in a static state and the heat transfer Nusselt number of the molten pool in a swinging state, as well as the convection influence factor of the molten pool, and then the heat transfer characteristics of the molten pool can be determined.
[0133] Optionally, an alternative implementation of constructing a general heat transfer expression for the molten pool in this embodiment includes:
[0134] In this embodiment, the Nusselt number amplification factor is used to correlate the heat transfer Nusselt number of the molten pool in the static state with the heat transfer Nusselt number of the molten pool in the swing state, thereby determining the first heat transfer characteristic expression. The first heat transfer characteristic expression can be expressed by the following formula (8):
[0135]
[0136] In formula (8), f roll represents the Nusselt number amplification factor, Nu st The Nusselt number represents the heat transfer of the molten pool at a stationary state, Nu roll It represents the heat transfer Nusselt number of the molten pool in the rocking state.
[0137] In this embodiment, the second heat transfer characteristic expression is determined based on the Nusselt number amplification factor and the convection influence factor. The second heat transfer characteristic expression can be expressed by the following formula (9):
[0138] f roll =f(Ri r ) (9)
[0139] Among them, f in formula (9) roll represents the Nusselt number amplification factor, Ri r represents the Richardson number, which is the convection influence factor.
[0140] In this embodiment, the third heat transfer characteristic expression is determined based on the heat transfer Nusselt number of the molten pool in a static state. The third heat transfer characteristic expression can be expressed by the following formula (10):
[0141]
[0142] In formula (10), C and n are constants; g represents the acceleration due to gravity; Ra′ represents the modified Rayleigh number, a dimensionless quantity; β represents the thermal expansion coefficient; q v represents the power density of the heat source in the molten pool; h represents the height of the molten pool; α represents the thermal diffusivity of the molten pool; ν represents the kinematic viscosity of the melt; and λ represents the thermal conductivity of the melt.
[0143] According to the first heat transfer characteristic formula, the second heat transfer characteristic formula and the third heat transfer characteristic formula, a general heat transfer expression of the molten pool is determined. The general heat transfer expression of the molten pool can be expressed by the following formula (11):
[0144]
[0145] This example accurately determines the heat transfer characteristics of the melt pool by obtaining the Nusselt numbers for both the stationary and swinging states of the melt pool, along with the convection factor of the melt pool. This determined heat transfer characteristics provide a foundation for the subsequent structural design of the lower end cap.
[0146] In one embodiment, Figure 6 As shown, an optional implementation method of the melt pool characteristic analysis method is:
[0147] S601: Determine a natural convection intensity value of the molten pool in a static state based on first fluid parameters of the molten pool, wherein the first fluid parameters include a thermal expansion coefficient, a difference between a maximum temperature in the molten pool and a phase interface temperature, a height of the molten pool, and a kinematic viscosity of the melt.
[0148] S602: Determine a fluid momentum equation based on a second fluid parameter of the molten pool.
[0149] S603: Determine dimensionless parameters based on characteristic parameters and third fluid parameters of the molten pool. The characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameters include the molten pool's radius vector, velocity, time, fluid density, sway angular velocity, pressure, and gravitational acceleration; and the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration.
[0150] S604 , performing dimensionless processing on the fluid momentum equation according to the characteristic parameters and the dimensionless parameters to obtain a dimensionless momentum equation.
[0151] S605: Determine the forced convection intensity value of the molten pool in the swinging state according to the dimensionless momentum equation.
[0152] S606 , square the forced convection intensity value, and multiply the squared result by the natural convection intensity value to obtain a convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt.
[0153] S607: Determine whether there is forced convection in a swaying state inside the molten pool based on the relationship between the convection influencing factor and the preset threshold.
[0154] S608, if it exists, construct a partial differential equation of the relative velocity potential function of the free liquid surface and the boundary conditions of the partial differential equation based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in the swinging state.
[0155] S609, based on the boundary conditions, solving the partial differential equation to obtain the convection characteristics of the molten pool.
[0156] S6010, obtaining the heat transfer Nusselt number of the molten pool in a stationary state and the heat transfer Nusselt number in a swinging state.
[0157] S6011: Determine the heat transfer characteristics of the molten pool based on the heat transfer Nusselt number of the molten pool in a static state, the heat transfer Nusselt number of the molten pool in a swinging state, and the convection influence factor of the molten pool.
[0158] This embodiment obtains the natural convection intensity value of the molten pool in a static state and the forced convection intensity value of the molten pool in a swinging state, and then determines the convection influence factor of the molten pool based on the natural convection intensity value and the forced convection intensity value. Finally, based on the relationship between the convection influence factor and the preset threshold value, it is determined whether there is forced convection in the swinging state inside the molten pool; if so, the convection characteristics of the molten pool in the swinging state can be determined based on the height and lateral span of the free liquid surface in the molten pool, which solves the problem of the current lack of effective methods for analyzing molten pool characteristics, especially the lack of methods for determining the characteristics of the molten pool in the swinging state.
[0159] It should be understood that, although the various steps in the flowcharts involved in the various embodiments described above are displayed in sequence according to the instructions of the arrows, these steps are not necessarily executed in sequence in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order restriction on the execution of these steps, and these steps can be executed in other orders. Moreover, at least a portion of the steps in the flowcharts involved in the various embodiments described above can include multiple steps or multiple stages, and these steps or stages are not necessarily executed and completed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a portion of steps or stages in other steps.
[0160] Based on the same inventive concept, embodiments of the present application also provide a melt pool characteristic analysis device for implementing the melt pool characteristic analysis method described above. The solution provided by this device is similar to the solution described in the method described above. Therefore, the specific limitations of one or more melt pool characteristic analysis device embodiments provided below can be found in the limitations of the melt pool characteristic analysis method described above and will not be repeated here.
[0161] In one embodiment, Figure 7 As shown, a molten pool characteristic analysis device 1 is provided, comprising: a first acquisition module 10, a first determination module 20, a second determination module 30 and a third determination module 40, wherein:
[0162] The first acquisition module 10 is used to obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0163] A first determining module 20 is used to determine a convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0164] A second determination module 30 is configured to determine whether there is forced convection in a swaying state inside the molten pool based on a relationship between the convection influencing factor and a preset threshold;
[0165] The third determining module 40, if present, determines the convection characteristics of the molten pool based on the height and transverse span of the free liquid surface in the molten pool.
[0166] In one embodiment, the Figure 7 The first acquisition module 10 is specifically used to determine the natural convection intensity value of the molten pool in a static state based on the first fluid parameters of the molten pool; wherein the first fluid parameters include the thermal expansion coefficient, the difference between the maximum temperature in the molten pool and the phase interface temperature, the height of the molten pool and the kinematic viscosity of the melt.
[0167] In one embodiment, the Figure 7 On the basis of Figure 8 As shown, Figure 7 The first acquisition module 10 further includes:
[0168] A first determining unit 101 is configured to determine a fluid momentum equation based on second fluid parameters of the molten pool, wherein the second fluid parameters include a velocity of the molten pool, a fluid density, a pressure, a kinematic viscosity of the melt, and a rocking angular velocity;
[0169] The second determining unit 102 is configured to perform dimensionless processing on the fluid momentum equation to obtain a dimensionless momentum equation;
[0170] The third determining unit 103 is configured to determine the forced convection intensity value of the molten pool in the swaying state according to the dimensionless momentum equation.
[0171] In one embodiment, the Figure 8 The second determination unit 102 is specifically used to: determine the dimensionless parameters according to the characteristic parameters and the third fluid parameters of the molten pool; wherein the characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameters include the radius vector, motion velocity, time, fluid density, swing angular velocity, pressure, and gravitational acceleration of the molten pool; the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration; and perform dimensionless processing on the fluid momentum equation according to the characteristic parameters and the dimensionless parameters to obtain the dimensionless momentum equation.
[0172] In one embodiment, the Figure 7 The first determination module 20 is specifically used to: square the forced convection intensity value, and multiply the squared result by the natural convection intensity value to obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt.
[0173] In one embodiment, the Figure 7 On the basis of Figure 9 As shown, Figure 7 The third determining module 40 further includes:
[0174] A construction unit 401 is used to construct a partial differential equation of a relative velocity potential function of the free liquid surface and boundary conditions of the partial differential equation according to the height and lateral span of the free liquid surface in the molten pool;
[0175] The fourth determining unit 402 is configured to solve the partial differential equation based on the boundary conditions to obtain the convection characteristics of the molten pool.
[0176] In one embodiment, the Figure 7 A molten pool characteristic analysis device 1, further comprising:
[0177] The second acquisition module 50 is used to obtain the heat exchange Nusselt number of the molten pool in a static state and the heat exchange Nusselt number in a swinging state;
[0178] The fourth determining module 60 is configured to determine the heat transfer characteristics of the molten pool according to the heat transfer Nusselt number of the molten pool in a stationary state, the heat transfer Nusselt number of the molten pool in a swinging state, and the convection influence factor of the molten pool.
[0179] Each module in the aforementioned melt pool characteristic analysis device may be implemented in whole or in part through software, hardware, or a combination thereof. Each module may be embedded in or independent of a processor in a computer device in the form of hardware, or may be stored in a computer device memory in the form of software, so that the processor can call and execute the corresponding operations of each module.
[0180] In one embodiment, a computer device is provided. The computer device may be a server, and its internal structure diagram may be as follows: Figure 10 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O) and a communication interface. The processor, the memory and the input / output interface are connected via a system bus, and the communication interface is connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store relevant data of the molten pool. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, a molten pool characteristic analysis method is implemented.
[0181] Those skilled in the art will understand that Figure 10 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.
[0182] In one embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, the following steps are implemented:
[0183] Obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0184] Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0185] According to the relationship between the convection influencing factor and the preset threshold, it is determined whether there is forced convection in a swaying state inside the molten pool;
[0186] If it exists, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in the swaying state.
[0187] In one embodiment, when the processor executes the computer program, the processor further implements the following steps: obtaining the natural convection intensity value of the molten pool in a static state, including:
[0188] The natural convection intensity value of the molten pool in a static state is determined based on the first fluid parameters of the molten pool; wherein the first fluid parameters include the thermal expansion coefficient, the difference between the maximum temperature in the molten pool and the phase interface temperature, the height of the molten pool and the kinematic viscosity of the melt.
[0189] In one embodiment, when the processor executes the computer program, the processor further implements the following steps: obtaining a forced convection intensity value of the molten pool in a swinging state, including:
[0190] Determining a fluid momentum equation based on second fluid parameters of the molten pool; wherein the second fluid parameters include a velocity of the molten pool, a fluid density, a pressure, a kinematic viscosity of the melt, and a rocking angular velocity;
[0191] The fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation;
[0192] According to the dimensionless momentum equation, the forced convection intensity value of the molten pool in the swing state is determined.
[0193] In one embodiment, when the processor executes the computer program, the processor further implements the following steps: performing dimensionless processing on the fluid momentum equation to obtain the dimensionless momentum equation, including:
[0194] Determine dimensionless parameters based on characteristic parameters and third fluid parameters of the molten pool; wherein the characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameters include the radius vector, motion velocity, time, fluid density, sway angular velocity, pressure, and gravitational acceleration of the molten pool; and the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration;
[0195] According to the characteristic parameters and dimensionless parameters, the fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation.
[0196] In one embodiment, when the processor executes the computer program, the processor further implements the following steps: determining the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value, including:
[0197] The forced convection intensity value is squared and multiplied by the natural convection intensity value to obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt.
[0198] In one embodiment, when the processor executes the computer program, the processor further implements the following steps: determining the convection characteristics of the molten pool based on the height and lateral span of the free liquid surface in the molten pool, including:
[0199] According to the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in the swaying state, the partial differential equation of the relative velocity potential function of the free liquid surface and the boundary conditions of the partial differential equation are constructed;
[0200] Based on the boundary conditions, the partial differential equations are solved to obtain the convection characteristics of the molten pool.
[0201] In one embodiment, when the processor executes the computer program, the method further implements the following steps:
[0202] Obtain the Nusselt number of the heat transfer in the molten pool in a static state and the Nusselt number of the heat transfer in a swinging state;
[0203] The heat transfer characteristics of the molten pool are determined based on the heat transfer Nusselt number of the molten pool in a static state, the heat transfer Nusselt number of the molten pool in a swinging state, and the convection influence factor of the molten pool.
[0204] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented:
[0205] Obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0206] Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0207] According to the relationship between the convection influencing factor and the preset threshold, it is determined whether there is forced convection in a swaying state inside the molten pool;
[0208] If it exists, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in the swaying state.
[0209] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: obtaining the natural convection intensity value of the molten pool in a static state, including:
[0210] The natural convection intensity value of the molten pool in a static state is determined based on the first fluid parameters of the molten pool; wherein the first fluid parameters include the thermal expansion coefficient, the difference between the maximum temperature in the molten pool and the phase interface temperature, the height of the molten pool and the kinematic viscosity of the melt.
[0211] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: obtaining a forced convection intensity value of the molten pool in a swinging state, including:
[0212] Determining a fluid momentum equation based on second fluid parameters of the molten pool; wherein the second fluid parameters include a velocity of the molten pool, a fluid density, a pressure, a kinematic viscosity of the melt, and a rocking angular velocity;
[0213] The fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation;
[0214] According to the dimensionless momentum equation, the forced convection intensity value of the molten pool in the swing state is determined.
[0215] In one embodiment, when the computer program is executed by a processor, the computer program further implements the following steps: performing dimensionless processing on the fluid momentum equation to obtain a dimensionless momentum equation, including:
[0216] Determine dimensionless parameters based on characteristic parameters and third fluid parameters of the molten pool; wherein the characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameters include the radius vector, motion velocity, time, fluid density, sway angular velocity, pressure, and gravitational acceleration of the molten pool; and the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration;
[0217] According to the characteristic parameters and dimensionless parameters, the fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation.
[0218] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: determining the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value, including:
[0219] The forced convection intensity value is squared and multiplied by the natural convection intensity value to obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt.
[0220] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: determining the convection characteristics of the molten pool based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state, including:
[0221] According to the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in the swaying state, the partial differential equation of the relative velocity potential function of the free liquid surface and the boundary conditions of the partial differential equation are constructed;
[0222] Based on the boundary conditions, the partial differential equations are solved to obtain the convection characteristics of the molten pool.
[0223] In one embodiment, when the computer program is executed by a processor, the method further implements the following steps:
[0224] Obtain the Nusselt number of the heat transfer in the molten pool in a static state and the Nusselt number of the heat transfer in a swinging state;
[0225] The heat transfer characteristics of the molten pool are determined based on the heat transfer Nusselt number of the molten pool in a static state, the heat transfer Nusselt number of the molten pool in a swinging state, and the convection influence factor of the molten pool.
[0226] In one embodiment, a computer program product is provided, comprising a computer program, which, when executed by a processor, implements the following steps:
[0227] Obtain the natural convection intensity value of the molten pool in a static state and the forced convection intensity value in a swinging state;
[0228] Determine the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value;
[0229] According to the relationship between the convection influencing factor and the preset threshold, it is determined whether there is forced convection in a swaying state inside the molten pool;
[0230] If it exists, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in the swaying state.
[0231] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: obtaining the natural convection intensity value of the molten pool in a static state, including:
[0232] The natural convection intensity value of the molten pool in a static state is determined based on the first fluid parameters of the molten pool; wherein the first fluid parameters include the thermal expansion coefficient, the difference between the maximum temperature in the molten pool and the phase interface temperature, the height of the molten pool and the kinematic viscosity of the melt.
[0233] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: obtaining a forced convection intensity value of the molten pool in a swinging state, including:
[0234] Determining a fluid momentum equation based on second fluid parameters of the molten pool; wherein the second fluid parameters include a velocity of the molten pool, a fluid density, a pressure, a kinematic viscosity of the melt, and a rocking angular velocity;
[0235] The fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation;
[0236] According to the dimensionless momentum equation, the forced convection intensity value of the molten pool in the swing state is determined.
[0237] In one embodiment, when the computer program is executed by a processor, the computer program further implements the following steps: performing dimensionless processing on the fluid momentum equation to obtain a dimensionless momentum equation, including:
[0238] Determine dimensionless parameters based on characteristic parameters and third fluid parameters of the molten pool; wherein the characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameters include the radius vector, motion velocity, time, fluid density, sway angular velocity, pressure, and gravitational acceleration of the molten pool; and the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration;
[0239] According to the characteristic parameters and dimensionless parameters, the fluid momentum equation is dimensionlessly processed to obtain the dimensionless momentum equation.
[0240] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: determining the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value, including:
[0241] The forced convection intensity value is squared and multiplied by the natural convection intensity value to obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt.
[0242] In one embodiment, when the computer program is executed by the processor, the following steps are further implemented: determining the convection characteristics of the molten pool based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state, including:
[0243] According to the height and lateral span of the free liquid surface in the molten pool, as well as the horizontal acceleration and vertical acceleration of the molten pool in the swaying state, the partial differential equation of the relative velocity potential function of the free liquid surface and the boundary conditions of the partial differential equation are constructed;
[0244] Based on the boundary conditions, the partial differential equations are solved to obtain the convection characteristics of the molten pool.
[0245] In one embodiment, when the computer program is executed by a processor, the method further implements the following steps:
[0246] Obtain the Nusselt number of the heat transfer in the molten pool in a static state and the Nusselt number of the heat transfer in a swinging state;
[0247] The heat transfer characteristics of the molten pool are determined based on the heat transfer Nusselt number of the molten pool in a static state, the heat transfer Nusselt number of the molten pool in a swinging state, and the convection influence factor of the molten pool.
[0248] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, database or other media used in the embodiments provided in this application may include at least one of non-volatile and volatile memory. Non-volatile memory may include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory may include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The database involved in the various embodiments provided herein may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, distributed databases based on blockchains. The processor involved in the various embodiments provided herein may be, but are not limited to, a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic unit, a data processing logic unit based on quantum computing, and the like.
[0249] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.
[0250] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present application. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present application, and these modifications and improvements fall within the scope of protection of the present application. Therefore, the scope of protection of the present application shall be determined by the appended claims.
Claims
1. A method for analyzing molten pool characteristics, characterized in that: The method comprises: determining a natural convection intensity value of the molten pool in a static state based on first fluid parameters of the molten pool, wherein the first fluid parameters include a thermal expansion coefficient, a difference between a maximum temperature in the molten pool and a phase interface temperature, a height of the molten pool, and a kinematic viscosity of the melt; Determining a fluid momentum equation based on second fluid parameters of the molten pool; performing dimensionless processing on the fluid momentum equation to obtain a dimensionless momentum equation; and determining a forced convection intensity value of the molten pool in a swaying state based on the dimensionless momentum equation; wherein the second fluid parameters include a molten pool velocity, fluid density, pressure, kinetic viscosity of the melt, and swaying angular velocity; Determining a convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value; determining whether there is forced convection in the swaying state within the molten pool according to a relationship between the convection influencing factor and a preset threshold; If present, the convection characteristics of the molten pool are determined based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state.
2. The method according to claim 1, characterized in that The dimensionless processing of the fluid momentum equation to obtain the dimensionless momentum equation includes: Determine the dimensionless parameters based on the characteristic parameters and the third fluid parameters of the molten pool; wherein the characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameters include the radius vector, motion velocity, time, fluid density, sway angular velocity, pressure, and gravitational acceleration of the molten pool; the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration; the radius vector of the molten pool refers to the radius vector of the fluid particle relative to the rotation axis; The fluid momentum equation is dimensionally processed according to the characteristic parameters and the dimensionless parameters to obtain a dimensionless momentum equation.
3. The method according to claim 1, characterized in that Determining the convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value includes: The forced convection intensity value is squared, and the squared result is multiplied by the natural convection intensity value to obtain the convection influence factor of the molten pool after eliminating the kinematic viscosity of the melt.
4. The method according to claim 1, wherein Determining the convection characteristics of the molten pool according to the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state includes: Constructing a partial differential equation of the relative velocity potential function of the free liquid surface and boundary conditions of the partial differential equation based on the height and lateral span of the free liquid surface in the molten pool, and the horizontal acceleration and vertical acceleration of the molten pool in a swaying state; Based on the boundary conditions, the partial differential equation is solved to obtain the convection characteristics of the molten pool.
5. The method according to any one of claims 1 to 4, characterized in that The method further comprises: Obtaining the heat transfer Nusselt number of the molten pool in a static state and the heat transfer Nusselt number in a swinging state; The heat transfer characteristics of the molten pool are determined according to the heat transfer Nusselt number of the molten pool in a stationary state and the heat transfer Nusselt number of the molten pool in a swinging state, as well as the convection influence factor of the molten pool.
6. A molten pool characteristic analysis device, characterized in that: include: an acquisition module for determining a natural convection intensity value of the molten pool in a static state based on first fluid parameters of the molten pool, wherein the first fluid parameters include a thermal expansion coefficient, a difference between a maximum temperature in the molten pool and a phase interface temperature, a height of the molten pool, and a kinematic viscosity of the melt; The acquisition module further includes a first determination unit, a second determination unit, and a third determination unit; wherein the first determination unit is used to determine a fluid momentum equation based on a second fluid parameter of the molten pool; the second fluid parameter includes a velocity of the molten pool, a fluid density, a pressure, a kinematic viscosity of the melt, and a rocking angular velocity; the second determination unit is used to perform dimensionless processing on the fluid momentum equation to obtain a dimensionless momentum equation; and the third determination unit is used to determine a forced convection intensity value of the molten pool in a rocking state based on the dimensionless momentum equation; A first determining module is used to determine a convection influence factor of the molten pool according to the natural convection intensity value and the forced convection intensity value; a second determining module, configured to determine whether there is forced convection in the swaying state inside the molten pool based on a relationship between the convection influencing factor and a preset threshold; The third determining module, if present, determines the convection characteristics of the molten pool according to the height and lateral span of the free liquid surface in the molten pool.
7. The device according to claim 6, characterized in that The second determination unit is further specifically used to: determine dimensionless parameters based on characteristic parameters and a third fluid parameter of the molten pool; wherein the characteristic parameters include characteristic length, characteristic velocity, characteristic time, characteristic density, characteristic angular velocity, and characteristic gravitational acceleration; the third fluid parameter includes the radius vector, motion velocity, time, fluid density, swing angular velocity, pressure, and gravitational acceleration of the molten pool; the dimensionless parameters include dimensionless radius, dimensionless velocity, dimensionless time, dimensionless density, dimensionless angular velocity, dimensionless pressure, and dimensionless gravitational acceleration; according to the characteristic parameters and the dimensionless parameters, perform dimensionless processing on the fluid momentum equation to obtain a dimensionless momentum equation; the radius vector of the molten pool refers to the radius vector of the fluid particle relative to the rotation axis.
8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the melt pool characteristic analysis method according to any one of claims 1 to 5 are implemented.
9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the melt pool characteristic analysis method according to any one of claims 1 to 5 are implemented.
10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the melt pool characteristic analysis method according to any one of claims 1 to 5 are implemented.
Citation Information
Patent Citations
Large nuclear reactor molten pool natural-convection heat transfer test system and method
CN105911087A
Heat exchange calculation method and device for melting pool and electronic equipment
CN112964397A