Multi-phase flow simulation method and device for cooling equipment, equipment and medium

Abstract

The embodiment of the invention provides a multiphase flow simulation method and device for cooling equipment, equipment and a medium, and the method comprises the steps: obtaining a geometric structure model and input working condition parameters of to-be-simulated cooling equipment, and carrying out the data preprocessing; performing simulation pretreatment according to preset structural characteristics of the cooling equipment, the pretreated geometric structure model and the input working condition parameters to construct a multiphase flow simulation model of the to-be-simulated cooling equipment; and performing simulation calculation and post-processing by adopting the multiphase flow simulation model to obtain a simulation result. According to the embodiment of the invention, through the finite element mesh model constructed based on the geometric structure parameters and the input working condition parameters, the simulation model can be quickly and automatically established, and the simulation efficiency and accuracy are improved.

Description

TECHNICAL FIELD

[0001] Embodiments of the present application relate to the technical field of simulation, in particular to a multi-phase flow simulation method and device for a cooling device, equipment and medium. BACKGROUND

[0002] The cooling device is a core component in the battery thermal management system of an electric vehicle. The cooling device is constructed by stacking multiple plate heat exchanger plates, and is divided into a refrigerant circuit and a coolant circuit. By introducing refrigerant from an air conditioning system and evaporating when flowing through an expansion valve, the refrigerant and the coolant flow in the cooling device in a convection form. During the convection process, the refrigerant absorbs the heat of the coolant in the battery cooling circuit through heat exchange, thereby achieving effective heat exchange and cooling the battery.

[0003] The cooling efficiency of the power battery of an electric vehicle is affected by the power of the cooling device, the power of the water pump, the flow rate of the refrigerant, and the flow rate of the coolant, etc. The multi-phase flow simulation of the cooling device can provide a large amount of result data in a short time. When properly set and verified, the simulation can provide relatively accurate results. Compared with the test method, the simulation can effectively predict parameters that are difficult to measure in the test, such as flow rate distribution, temperature distribution, and pressure distribution. In addition, the simulation can evaluate multiple working conditions and designs in a short time, thereby reducing the cost required for the test.

[0004] However, the existing simulation method for the cooling device requires a simulation engineer to manually set the simulation process, including steps such as geometry modeling, boundary condition setting, and mesh division. In particular, for high-efficiency cooling devices (chiller devices) with a multi-layer flow channel plate structure, the pre-processing process before simulation is complex, and parameters need to be set separately for different flow channel layers, which results in tedious operation and is prone to errors.

[0005] This manual simulation method has problems such as low efficiency, large amount of repetitive work, and high risk of human error. At the same time, due to the lack of standardized simulation process, the simulation results of different engineers may have large differences, which affects the accuracy and reliability of the design scheme.

[0006] In the field of multi-phase flow simulation of cooling devices, the existing technology mainly faces the following technical problems: 1) low simulation efficiency, the traditional method requires a large amount of manual operation and repetitive work; 2) insufficient simulation accuracy, it is difficult to accurately predict the multi-phase flow characteristics in the flow channel; 3) low standardization of simulation process, the simulation results of different engineers have large differences. These problems seriously restrict the research and development efficiency and product quality of the cooling device. SUMMARY

[0007] Embodiments of the present application provide a multiphase flow simulation method, device, equipment and medium for a cooling device, to solve the problem that the existing simulation method for the cooling device needs a simulation engineer to manually set a simulation process operation, increases errors caused by human input, and the simulation process is inefficient.

[0008] According to an aspect of the present application, a multiphase flow simulation method for a cooling device is provided, comprising:

[0009] obtaining a geometric structure model and input working condition parameters of a cooling device to be simulated, and performing data preprocessing;

[0010] performing pre-simulation processing according to preset structural features of the cooling device, the preprocessed geometric structure model and the input working condition parameters, to construct a multiphase flow simulation model of the cooling device to be simulated;

[0011] performing simulation calculation and post-processing using the multiphase flow simulation model, to obtain a simulation result;

[0012] The preset structural features of the cooling device include a plurality of layers of flow channel plates stacked between a top plate and a bottom plate, a refrigerant inlet and a refrigerant outlet opened in the top plate or the bottom plate, and a cooling liquid inlet and a cooling liquid outlet opened in the top plate or the bottom plate; the refrigerant inlet forms a refrigerant flow channel inner cavity through refrigerant flow channels in the layers of flow channel plates to the refrigerant outlet; and the cooling liquid inlet forms a cooling liquid flow channel inner cavity through cooling liquid flow channels in the layers of flow channel plates to the cooling liquid outlet.

[0013] According to another aspect of the present application, a multiphase flow simulation device for a cooling device is provided, comprising:

[0014] a preprocessing module configured to obtain a geometric structure model and input working condition parameters of a cooling device to be simulated, and perform data preprocessing;

[0015] a model construction module configured to perform pre-simulation processing according to preset structural features of the cooling device, the preprocessed geometric structure model and the input working condition parameters, to construct a multiphase flow simulation model of the cooling device to be simulated;

[0016] a result obtaining module configured to perform simulation calculation and post-processing using the multiphase flow simulation model, to obtain a simulation result;

[0017] The preset structural features of the cooling device include: a plurality of layers of flow channel plates stacked between a top plate and a bottom plate; a refrigerant inlet and a refrigerant outlet formed in the top plate or the bottom plate, and a cooling liquid inlet and a cooling liquid outlet formed in the top plate or the bottom plate; the refrigerant inlet forms a refrigerant flow channel inner cavity by passing through the refrigerant flow channel in each layer of flow channel plates to the refrigerant outlet; and the cooling liquid inlet forms a cooling liquid flow channel inner cavity by passing through the cooling liquid flow channel in each layer of flow channel plates to the cooling liquid outlet.

[0018] According to another aspect of the present application, there is provided a computer program product comprising a computer program which, when executed by a processor, implements the multiphase flow simulation method of the cooling device according to any one of the embodiments of the present application.

[0019] According to another aspect of the present application, there is provided an electronic device comprising:

[0020] at least one processor; and a memory connected to the at least one processor in communication; wherein the memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor to enable the at least one processor to perform the multiphase flow simulation method of the cooling device according to any one of the embodiments of the present application.

[0021] According to another aspect of the present application, there is provided a computer readable storage medium storing computer instructions for causing a processor to implement the multiphase flow simulation method of the cooling device according to any one of the embodiments of the present application when executed by the processor.

[0022] The embodiments of the present application can quickly and conveniently construct a finite element grid model based on geometric structure parameters and input working condition parameters through the preset structural features of the cooling device, and can truly reflect the physical characteristics of the device; by configuring physical model parameters and post-processing acquisition parameters, the simulation process is ensured to be consistent with the actual working principle of the device, thereby obtaining a high-credibility multiphase flow simulation result, providing reliable data support for performance optimization, fault diagnosis and energy efficiency improvement of the cooling device, adapting to various changes and adjustments of the cooling device, significantly shortening the research and development cycle and reducing the experimental cost.

[0023] It should be understood that the content described in this part is not intended to identify key or important features of the embodiments of the present application, nor to limit the scope of the present application. Other features of the present application will become apparent from the following description. BRIEF DESCRIPTION OF DRAWINGS

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the description of the embodiments. Obviously, the drawings described below only show some of the embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor based on these drawings.

[0025] Figures 1A-1C and Figure 1E Three views of the cooling device suitable for the embodiments of the present application;

[0026] Figure 1D Flowchart of the multiphase flow simulation method of the cooling device provided by the embodiments of the present application;

[0027] Figure 2 Flowchart of the geometric structure model preprocessing in the multiphase flow simulation method of the cooling device provided by the embodiments of the present application;

[0028] Figure 3 Flowchart of the input working condition parameter preprocessing in the multiphase flow simulation method of the cooling device provided by the embodiments of the present application;

[0029] Figure 4 Flowchart of the simulation preprocessing in the multiphase flow simulation method of the cooling device provided by the embodiments of the present application;

[0030] Figure 5 Flowchart of the simulation calculation and post-processing in the multiphase flow simulation method of the cooling device provided by the embodiments of the present application;

[0031] Figure 6 Flowchart of the function of the plug-in corresponding to the multiphase flow simulation method of the cooling device provided by the embodiments of the present application;

[0032] Figure 7 Structural schematic diagram of the multiphase flow simulation device of the cooling device provided by the embodiments of the present application;

[0033] Figure 8 Structural schematic diagram of the electronic device provided by the embodiments of the present application. DETAILED DESCRIPTION

[0034] In order to make the person skilled in the art better understand the present application, the following will combine the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments only show some of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor should be within the scope of protection of the present application.

[0035] It is to be understood that the terms "first", "second", and the like used in the description and the claims of the present application and the above drawings are used to distinguish similar objects, and are not necessarily used to describe a particular sequential or chronological order. It should be understood that the data thus used can be interchanged, where appropriate, so that the embodiments of the application described herein can be carried out in other than the order shown or described herein. Furthermore, the terms "comprise" and "have", and any variations thereof, are intended to cover non-exclusive inclusion, for example, a process, method, system, product, or apparatus that comprises a list of steps or units is not necessarily limited to those steps or units that are clearly listed, but can include other steps or units that are not clearly listed or inherent to such processes, methods, products, or apparatuses.

[0036] The embodiments of the present application are proposed to solve the above technical problems, and a multiphase flow simulation method of a cooling device is provided, which realizes a technical breakthrough through the following core improvement points:

[0037] 1) Geometric model preprocessing technology is adopted to automatically extract flow channel structure characteristics;

[0038] 2) An automatic selection mechanism for multiphase flow models is developed;

[0039] 3) Automatic processing of the simulation process is realized.

[0040] These improvements significantly improve the simulation efficiency and accuracy, while reducing the errors caused by manual intervention.

[0041] The key terms involved in the embodiments of the present application are explained as follows:

[0042] 1) Multiphase flow simulation: refers to a technology for numerically simulating the flow of fluid containing two or more different phases (such as gas-liquid two-phase). In a cooling device, it is mainly used to simulate the gas-liquid phase change process of refrigerant.

[0043] 2) Geometric structure model: a digital model describing the three-dimensional geometry of a cooling device, including flow channel plates, inlet and outlet pipes, and other structural features, which can fully represent the internal and external three-dimensional contour structure of the cooling device.

[0044] 3) Input working condition parameters: the working condition parameters of the cooling device to be simulated during operation, which are the boundary condition parameters required for simulation calculation, including flow rate, temperature, pressure, and other physical quantities.

[0045] 4) Multiphase flow model: a mathematical model used to describe multiphase flow, mainly involving mixed multiphase flow model (MMP), two-phase thermal equilibrium model (TPTE), and Euler multiphase separation flow model (EMP) in the present application.

[0046] A typical cooling device includes a Chiller device for refrigerating a battery system. The three views of the cooling device are as followsFigures 1A-1C As shown, the cooling device includes a plurality of flow channel plates 10 arranged in layers, such as Figures 1A-1C As shown in the middle, the cooling device includes a first flow channel plate 1, a second flow channel plate 2, and a plurality of flow channel plates 10 arranged below, and further includes a top plate 8 and a bottom plate 3, wherein the top plate 8 is provided with a cooling liquid inlet 4 and a refrigerant inlet 6, and further provided with a refrigerant outlet 5 and a cooling liquid outlet 7. In fact, the cooling liquid inlet and outlet, and the refrigerant inlet and outlet can be provided on one of the top plate or the bottom plate or on both of them according to the needs. Various patterns of flow channels are provided in each flow channel plate 10, and the flow channels of adjacent layers of flow channel plates are buckled together to form flow channels. The flow channels respectively communicating with the inlet and outlet of the cooling liquid form cooling liquid flow channels, and the flow channels respectively communicating with the inlet and outlet of the refrigerant form refrigerant flow channels. Optionally, the flow channels between the first flow channel plate and the second flow channel plate form cooling liquid flow channels if they communicate with the cooling liquid inlet, the flow channels between the second flow channel plate and the third flow channel plate form refrigerant flow channels if they communicate with the refrigerant inlet, the flow channels between the third flow channel plate and the fourth flow channel plate form the next cooling liquid flow channels, and the flow channels between the fourth flow channel plate and the fifth flow channel plate form the next refrigerant flow channels. In this way, the refrigerant flow channels communicate with each other, the cooling liquid flow channels communicate with each other, and the two types of flow channels form two side flow channels and are spaced apart from each other to facilitate heat conduction.

[0047] The technical scheme of the embodiment of the present application realizes efficient and accurate calculation of the multiphase flow simulation of the cooling device through the innovative means of automatically processing the geometric structure model, intelligently selecting the multiphase flow model, optimizing the grid division, and automatically processing the configuration.

[0048] Figure 1D is a flow diagram of the multiphase flow simulation method of the cooling device provided by the embodiment of the present application. The embodiment of the present application can be applied to the multiphase flow simulation of the cooling device, and can be executed by a cooling device multiphase flow simulation device. The device can be configured in any electronic device with computing function. Typically, it can be executed by the STAR-CCM+ simulation software platform. Based on the multiphase flow simulation software, the plug-in for automatic simulation modeling and execution is developed to complete the multiphase flow simulation processing of the cooling device of the embodiment of the present application. The steps shown below will be described. Figure 1D As shown in Figure 1D The method includes the following steps.

[0049] In step 110, the geometric structure model of the cooling device to be simulated and the input working condition parameters are obtained, and data preprocessing is performed.

[0050] STAR-CCM+ simulation software platform is a powerful CFD simulation tool that supports multi-physical field coupling calculation, has perfect meshing, solver and post-processing functions. In practical application, STAR-CCM+ can run on a server or a high-performance computing workstation.

[0051] The geometric structure model refers to a digitized model describing the three-dimensional geometric characteristics of the cooling equipment, including the top plate, the bottom plate, the multi-layer flow channel plate, and the inlet and outlet pipelines of the refrigerant and the cooling liquid and the like structures. Specifically, the geometric structure model can be established by a CAD software (such as SolidWorks, CATIA, etc.) and imported into the simulation platform in a general format such as STEP or IGES. The geometric structure model can be marked with different functional areas such as the top plate, the flow channel plate, the bottom plate, and the cooling liquid inlet and outlet, the refrigerant inlet and outlet.

[0052] The input working condition parameters include the physical property parameters of the refrigerant and the cooling liquid, the boundary condition parameters and the like. The physical property parameters of the refrigerant and the cooling liquid mainly include density, viscosity, specific heat capacity, thermal conductivity and the like, which can affect the flow, heat conduction and the like performance changes of the fluid; the boundary condition parameters include the flow rate, temperature, pressure and the like of the refrigerant and the cooling liquid inlet. These parameters are usually determined by experimental measurement or engineering experience, and can be set according to the simulation requirements.

[0053] The data preprocessing includes geometric model simplification and working condition parameter adjustment processing. The geometric model simplification mainly removes small features such as chamfers, small holes and the like that do not affect the simulation results; the working condition parameter adjustment can include parameter standardization, converting parameters in different units into standard units recognizable by the simulation software, and can also be calculating the required parameters for model simulation according to the provided physical property parameters and boundary condition parameters and the like. Through preprocessing, the efficiency and stability of subsequent simulation calculation can be improved, and preparation is made for the configuration of model construction, pre-processing and post-processing.

[0054] In actual implementation, the complete geometric model of the cooling equipment can be first exported from the CAD system, and then the geometric repair tool of STAR-CCM+ is used to check and repair the model to ensure the geometric closure and continuity. At the same time, the working condition parameters are arranged into a structured data format for subsequent calling.

[0055] In step 120, according to the preset structure characteristics of the cooling equipment, the pre-processed geometric structure model and the input working condition parameters, the simulation pre-processing is performed to construct a multiphase flow simulation model of the cooling equipment to be simulated.

[0056] The preset structural features refer to common structural features of the cooling device, including the plurality of layers of flow channel plates 10 stacked between the top plate 8 and the bottom plate 3, the refrigerant inlet 6 and the refrigerant outlet 5 formed in the top plate 8 or the bottom plate 3, and the cooling liquid inlet 4 and the cooling liquid outlet 7 formed in the top plate 8 or the bottom plate 3, and the refrigerant flow channel formed from the refrigerant inlet 6 to the refrigerant outlet 5 through the refrigerant flow channels in the plurality of layers of flow channel plates 10 to form a refrigerant flow channel inner cavity, the cooling liquid flow channel formed from the cooling liquid inlet 4 to the cooling liquid outlet 7 through the cooling liquid flow channels in the plurality of layers of flow channel plates 10 to form a cooling liquid flow channel inner cavity, and the like, as shown in FIG. 8. These features are common structural features of the cooling device, and the number of layers of flow channel plates, the size and position of the inlet and outlet, and the pattern of the flow channel inner cavity in the flow channel plate of different cooling devices are different, but the common structural features are consistent, and the finite element grid of the multiphase flow simulation model of the cooling device can be automatically constructed based on the common structural features. Figure 1E The common structural features are shown in FIG. 8. These features are common structural features of the cooling device, and the number of layers of flow channel plates, the size and position of the inlet and outlet, and the pattern of the flow channel inner cavity in the flow channel plate of different cooling devices are different, but the common structural features are consistent, and the finite element grid of the multiphase flow simulation model of the cooling device can be automatically constructed based on the common structural features.

[0057] Based on the preset structural features and the geometric structure model, the grid division is performed, the physical model is selected, and then the model grid is configured based on the input working condition parameters and the like for simulation pretreatment configuration, so as to construct the multiphase flow simulation model. For the grid division, the strategy of polyhedral grid cooperating with boundary layer grid can be used in the embodiments of the present application. The polyhedral grid is suitable for complex geometry and has high calculation efficiency; the boundary layer grid can accurately capture the near-wall flow characteristics.

[0058] For the selection of the physical model, the embodiments of the present application provide three kinds of multiphase flow models: mixed multiphase flow (MMP) model, two-phase thermal equilibrium (TPTE) model, and Euler multiphase separated flow (EMP) model. The MMP model has high calculation efficiency and is suitable for rapid evaluation; the TPTE model can accurately predict the phase change process; and the EMP model can describe the interaction between phases in detail. Through multiple simulation calculations or comparison of experimental data, the optimal model can be selected.

[0059] In actual implementation, first, the pretreated geometric model is imported into STAR-CCM+, then the grid parameters are set and the calculation grid is generated. Then the continuum is created, the multiphase flow model is selected and the physical property parameters are set. Finally, the calculation domain is created and the boundary conditions are set. The whole process can be realized automatically through secondary development scripts.

[0060] In step 130, simulation calculation and post-processing are performed by using the multiphase flow simulation model to obtain simulation results.

[0061] The simulation calculation refers to the process of solving the multiphase flow control equation by using a numerical method. For example, the finite volume method can be used for discrete solution in the embodiments of the present application.

[0062] The post-processing can be configured in advance, mainly including flow field visualization and quantitative analysis. The distribution of parameters such as velocity, temperature and pressure can be observed by creating cross sections. In particular, for the refrigerant side, the gas-liquid volume fraction distribution needs to be analyzed, which directly reflects the phase change process. Quantitative analysis includes calculating key performance indicators such as heat transfer and pressure drop.

[0063] For example, in actual implementation, first set the solver parameters and start the calculation. Monitor the residual curve and key parameter changes during the calculation to ensure convergence. After the calculation is completed, create the required cross sections and reports, and extract the key performance parameters. Finally, compare the results with the experimental data for verification.

[0064] The cooling equipment multiphase flow simulation method provided by the embodiments of the present application realizes efficient and accurate simulation analysis through innovative means such as optimization of geometric structure model processing, intelligent selection of multiphase flow model, and automated simulation process. Compared with traditional methods, the present scheme shortens the simulation time while improving the simulation accuracy, providing reliable technical support for the optimization design of cooling equipment. Specifically, the standardization processing of pre-set structure characteristics ensures the consistency of the simulation model; the intelligent selection of the multiphase flow model improves the calculation accuracy; the automated process reduces human error and improves work efficiency.

[0065] Figure 2 The flowchart for geometric structure model preprocessing in the cooling equipment multiphase flow simulation method provided by the embodiments of the present application. In the above embodiments, the function of obtaining the geometric structure model of the cooling equipment to be simulated and preprocessing the data can be realized by the following operations:

[0066] Step 111: According to the functional area annotation information in the geometric structure model, determine the number of flow channel plates, and extract the local structure model of the following parts: the first two layers of flow channel plates, the top plate, the bottom plate, the refrigerant inlet, the refrigerant outlet, the cooling liquid inlet and the cooling liquid outlet.

[0067] The geometric model refers to a three-dimensional digital model of the cooling device established by CAD software, which contains all the geometric features and structural information of the device. The model is usually constructed using parametric modeling methods, which facilitates subsequent modification and optimization. Each component in the geometric model has functional area annotation information, which identifies the functional attributes of different components, such as flow channel plates, inlets, and outlets. The multi-layer flow channel plate is the core heat exchange component in the cooling device, composed of multiple layers of metal plates, with complex internal flow channel structures for the flow and heat exchange of refrigerant and cooling liquid. The local structure model refers to a simplified model of a specific component extracted from the complete geometric model, which retains the main functional features of the component and removes unnecessary details. The first two layers of flow channel plates are extracted because in practice, the cooling device is composed of two types of flow channel plates stacked alternately. When the number of flow channel plates increases, due to the complex surface patterns and numerous interference points, the loading time of the digital model is long, and other structural software is difficult to process the internal cavity. Therefore, the first two layers of flow channel plates are used as representatives, which can be directly copied and stacked in the simulation software, reducing the time spent in other software and saving computer resources, allowing the simulation of multi-layer plates to continue.

[0068] By extracting the local structure model of these key parts, the calculation amount can be greatly reduced while ensuring the simulation accuracy, and the simulation efficiency can be improved. At the same time, the extraction of the local structure model also provides convenience for subsequent optimization processing.

[0069] Step 112: Optimize the flow-independent features of each local structure model.

[0070] Flow-independent features refer to geometric features that have little effect on fluid flow characteristics, such as small fillets, chamfers, and minor protrusions. These features may exist in actual manufacturing, but they increase the number of grids in simulation calculations without significantly affecting the calculation results. Optimization includes simplifying or deleting these flow-independent features, as well as necessary repair and fairing processing of the model. Specific optimization methods may include, for example: deleting holes with a diameter less than 1mm, simplifying fillets with a radius less than 0.5mm, merging adjacent surfaces with a spacing less than 0.3mm, etc. The optimized model retains key geometric features that affect flow, such as flow channel cross-sectional shape, flow direction changes, etc., while removing unnecessary small features.

[0071] This optimization can significantly reduce the number of grids, improve calculation efficiency, and ensure the accuracy of simulation results.

[0072] In actual implementation, the key components are identified and the local model is extracted by marking the information of the functional area, which provides a basis for subsequent optimization processing; and the optimization processing of small features further improves the simulation efficiency of the extracted local model. The two steps together constitute the core link of data preprocessing, laying a good foundation for subsequent simulation calculation.

[0073] Specifically, the basic size parameters of mesh division can be obtained from the geometric size of the actual Chiller product, the setting parameters of the boundary layer are calculated from the inlet flow rate and the inlet radius, and the relative position parameters of the cross section required for post-processing are obtained from the number and size of the flow channels, the distance between the two adjacent flow channels of the same side fluid, and the distance from the inlet of each side fluid to the first layer flow channel.

[0074] In the embodiments of the application, the local structure model of the key part is extracted and the flow-independent feature optimization processing is performed. In this way, the calculation amount can be greatly reduced while ensuring the simulation accuracy, the simulation efficiency is improved, so that accurate simulation results can be quickly obtained, and then reliable basis can be provided for the optimization design of the cooling equipment.

[0075] In some embodiments, the input working condition parameters can include: refrigerant, coolant and medium of the cooling equipment; refrigerant side inlet flow rate and temperature, refrigerant outlet pressure, pressure before the expansion valve, and gas-liquid volume fraction of the refrigerant inlet and outlet; inlet flow rate, temperature and pressure of the coolant side.

[0076] The refrigerant refers to the working medium used for heat exchange in the Chiller system, and common ones are R134a, R404A, etc. The refrigerant realizes heat transfer through the phase change process (evaporation and condensation) in the system, and its physical property parameters directly affect the heat exchange efficiency of the system. The refrigerant side inlet flow rate refers to the mass flow rate of the refrigerant entering the system per unit time, and the temperature reflects the thermal state of the refrigerant entering the system. The refrigerant outlet pressure is the pressure parameter of the refrigerant flowing out of the system, and the pressure before the expansion valve is the pressure state of the refrigerant before the expansion valve. These parameters jointly determine the flow and phase change behavior of the refrigerant in the system. The gas-liquid volume fraction of the refrigerant inlet and outlet reflects the phase distribution of the refrigerant at the inlet and outlet, and is an important index for evaluating the performance of the system.

[0077] The coolant refers to the liquid medium used for heat transfer in the battery thermal management system, and common ones are water, ethylene glycol solution, etc. The coolant side inlet flow rate refers to the volume flow rate of the coolant entering the system per unit time, the temperature reflects the thermal state of the coolant entering the system, and the pressure is the pressure parameter of the coolant entering the system. These parameters jointly determine the flow characteristics and heat exchange effect of the coolant in the system.

[0078] The medium of the cooling device refers to the material constituting the structure of the Chiller, which usually includes metal materials such as aluminum, copper, etc. The heat conductivity of these materials directly affects the overall heat exchange efficiency of the system. The selection of the medium needs to consider factors such as corrosion resistance, thermal conductivity, and mechanical strength.

[0079] In actual implementation, the parameters of the refrigerant side and the cooling liquid side need to be collected and input synchronously to ensure the completeness and accuracy of the system simulation. The refrigerant side parameters are mainly used to analyze the performance of the refrigeration cycle, while the cooling liquid side parameters are used to evaluate the effect of the thermal management system. Both of them together constitute the basic input conditions for the simulation of the Chiller system.

[0080] In the embodiments of the present application, the specific content of the input working condition parameters is defined in detail. In this way, comprehensive and accurate input data can be obtained during the simulation process, thereby improving the accuracy of the simulation results, and further providing reliable data support for the design and optimization of Chiller products.

[0081] Figure 3 The flow chart for pre-processing the input working condition parameters in the multiphase flow simulation method of the cooling device provided in the embodiments of the present application is shown in the figure. In some embodiments, the function of obtaining the input working condition parameters of the cooling device to be simulated and pre-processing the data can be realized by the following operations:

[0082] Step 113: According to the refrigerant, the cooling liquid, and the medium of the cooling device, the physical property parameters of the refrigerant, the cooling liquid, and the cooling device are obtained respectively;

[0083] The physical property parameters refer to the parameters describing the physical properties of the substance, including density, viscosity, specific heat capacity, thermal conductivity, etc. For the refrigerant, the physical property parameters can also include saturation pressure, latent heat of vaporization, and other phase change characteristic parameters; for the cooling liquid, the physical property parameters can include thermal physical property parameters; for the cooling device material, it can include thermal conductivity, density, etc. These parameters are the basic data for multiphase flow simulation, and directly affect the accuracy of the simulation results. In actual implementation, the physical property parameters can be obtained by consulting the physical property database or experimental measurement, such as the density of refrigerant R134a at 25℃ is 1206 kg / m³, and the thermal conductivity is 0.081 W / (m·K).

[0084] Step 114: According to the pre-expansion valve pressure, the refrigerant inlet gas-liquid volume fraction is obtained;

[0085] The gas-liquid volume fraction refers to the volume proportion of the gas phase and the liquid phase of the refrigerant in a specific state. The pressure before the expansion valve is a key parameter for determining the state of the refrigerant. By consulting the thermodynamic property table or the equation of state of the refrigerant, the saturation temperature and the proportion of the gas-liquid two-phase under the pressure can be determined. For example, when the pressure before the expansion valve of R134a is 1.5 MPa, the corresponding saturation temperature is 55°C, at this time the refrigerant is in a subcooled liquid state, and the gas-liquid volume fraction is 0 (full liquid state). When the pressure is reduced to 0.3 MPa, the gas-liquid volume fraction may be 0.3 (gas phase) and 0.7 (liquid phase). This parameter is crucial for the calculation of the phase change process in the subsequent multi-phase flow simulation.

[0086] Step 115: Calculate the inlet velocity of the cooling liquid side according to the inlet flow rate of the cooling liquid side.

[0087] The inlet velocity refers to the average flow rate of the cooling liquid when it enters the flow channel, that is, the flow rate at the inlet of the cooling liquid. It is calculated by dividing the inlet flow rate by the cross-sectional area of the flow channel. The inlet flow rate is usually given by the working condition parameters, with a unit of kg / s or m³ / s; the cross-sectional area of the flow channel is obtained according to the geometric model measurement. For example, when the cooling liquid flow rate is 0.5 kg / s, the cross-sectional area of the flow channel is 0.001 m², and the density of the cooling liquid is 1000 kg / m³, the inlet velocity is 0.5 m / s. The inlet velocity is an important parameter for setting boundary conditions, and directly affects the calculation accuracy of the flow field distribution and heat transfer effect.

[0088] In actual implementation, the first obtained physical property parameters provide basic physical characteristic data for subsequent calculation; the gas-liquid volume fraction determined by the pressure before the expansion valve reflects the state characteristics of the refrigerant; and the inlet velocity calculated based on the inlet flow rate provides a boundary condition for the flow field simulation. The three steps together constitute the pre-processing process of the multi-phase flow simulation, ensuring the completeness and accuracy of the simulation input.

[0089] In the embodiments of the present application, the physical property parameters of the refrigerant, the cooling liquid and the equipment are obtained respectively, the gas-liquid volume fraction of the refrigerant is determined, and the inlet velocity of the cooling liquid is calculated. In this way, accurate input parameters can be provided for the multi-phase flow simulation, so as to ensure the reliability of the simulation calculation, and thus the accuracy of the performance evaluation of the Chiller product can be improved.

[0090] Figure 4 The flow chart for the simulation pre-processing in the multi-phase flow simulation method of the cooling equipment provided in the embodiments of the present application is shown. In some embodiments, the simulation pre-processing can include the following operations according to the pre-set structural characteristics of the cooling equipment, the pre-processed geometric structure model and the input working condition parameters:

[0091] Step 121: Import the extracted local structure models of the refrigerant and coolant inlet and outlet, the first two layers of flow channel plates, the top plate, and the bottom plate from the pretreated geometric model.

[0092] The local structure model (Parts) refers to the key component model extracted from the complete geometric model, including the flow channel inlet and outlet, the flow channel plate, and other core structures. These local models retain the geometric characteristics and topological relationships of the original structure, but remove the small features that do not affect the simulation results. The extraction of local structure models can significantly reduce the computational load while ensuring the simulation accuracy.

[0093] Step 122: Surface repair of the imported local structure models.

[0094] Surface repair refers to the process of topological inspection and defect repair of the imported geometric model. Common repair operations include filling holes, repairing self-intersecting surfaces, and eliminating small gaps. Surface repair can ensure the smooth progress of subsequent meshing and simulation calculations.

[0095] Step 123: Establish a Cartesian coordinate system at the inlet and outlet of the refrigerant and coolant, respectively, where the X and Y axes of the Cartesian coordinate system are the inlet and outlet planes, and the Z axis is the flow direction of the flow channel.

[0096] The Cartesian coordinate system is a reference system used to define the spatial position and direction of the flow channel. The X-Y plane is parallel to the flow channel inlet and outlet cross-section, and the Z axis extends along the flow direction. This coordinate system setting facilitates the subsequent definition of flow channel arrays and boundary conditions. In practical applications, independent coordinate systems need to be established for the refrigerant side and the coolant side.

[0097] For example, at the inlet and outlet of the refrigerant and coolant flow channels on both sides, a Cartesian coordinate system is established, with the z-axis direction being the flow direction and perpendicular to the flow channel plate. The refrigerant inlet coordinate is named "csys_5", the refrigerant outlet coordinate is named "csys_10", the coolant inlet coordinate is named "csys_11", and the coolant inlet coordinate is named "csys_15". According to the actual number of flow channel plates, the first two layers of flow channel plates are arrayed along the z-direction of the refrigerant inlet coordinate system csys_5 to the actual structure state of the product, and the Part is then surface repaired and named "Chiller". The internal cavities of the refrigerant and coolant flow regions on both sides are then extracted and named "Refrigerant" and "Coolant", respectively.

[0098] Step 124: Along the Z-axis direction, according to the number of flow channel plates, the local structure models of the first two layers of flow channel plates are replicated and spliced to construct the simulation internal cavity structure model of the cooling equipment to be simulated.

[0099] The simulation inner cavity structure model is a complete flow channel system constructed by replicating and splicing the basic flow channel plate. The model retains the laminated structure characteristics of the actual product, while reducing the modeling workload. The replicating and splicing operation is performed along the Z-axis direction to ensure the consistency of the flow channel direction. In actual implementation, multiple layers of flow channel structures can be quickly generated through array operation. Specifically, the flow channel inner cavity structure between the first two layers of flow channel plates is used to replace the flow channel inner cavity structure between other flow channel plates.

[0100] Step 125: According to each local structure model and the simulation inner cavity structure model, extract the inner cavity structure model of the fluid region on both sides of the refrigerant and the coolant, respectively.

[0101] The inner cavity structure model refers to the spatial region model through which the fluid actually flows. Through the extraction operation, a pure fluid calculation domain can be obtained, removing the solid structure part. The refrigerant side and the coolant side need to be extracted separately to form independent calculation domains. These models will serve as the basis for subsequent mesh division and physical model definition. Specifically, the refrigerant pipeline between the refrigerant inlet and outlet extracted from the geometry model structure can be connected with the refrigerant side flow channel inner cavity structure in the flow channel plate to form the refrigerant inner cavity structure model; similarly, the coolant pipeline between the coolant inlet and outlet extracted from the geometry model structure can be connected with the coolant side flow channel inner cavity structure in the flow channel plate to form the coolant inner cavity structure model.

[0102] Step 126: Select a mesh division model and set mesh parameters.

[0103] The mesh division model determines the generation method and quality of the calculation mesh. Common models include polyhedral mesh, hexahedral mesh, etc. Mesh parameters include base size, boundary layer setting, etc., which directly affect the calculation accuracy and efficiency. In actual implementation, the optimal parameter combination can be determined through mesh independence research. Optionally, the mesh division selects a polyhedral mesh model, uses thin wall layer mesh for the solid side flow channel plate, and selects boundary layer mesh for the fluid region; when setting the mesh parameters, it can specifically include: keeping other parameter settings unchanged, gradually refining the mesh, comparing the calculation results of different numbers of meshes, and obtaining the mesh parameters when the mesh independent solution is obtained.

[0104] Step 127: Create a physical continuum for the inner cavity structure models of the refrigerant and the coolant, and select a multiphase flow model as the multiphase flow simulation model of the cooling device to be simulated.

[0105] The physical continuum defines the physical model and material properties used in the simulation. The multiphase flow model is used to describe the gas-liquid two-phase flow of the refrigerant, and common models include mixed multiphase flow model, Eulerian multiphase flow model, etc. Model selection needs to consider the balance between calculation accuracy and efficiency.

[0106] Step 128: Set physical property parameters for the multiphase flow simulation model according to the preprocessed input operating condition parameters.

[0107] The physical property parameters include physical property parameters such as density, viscosity, and specific heat of the fluid. These parameters need to be set according to the actual operating conditions and directly affect the accuracy of the simulation results. The parameter sources can be experimental measurement data or standard physical property databases. Specifically, the corresponding physical property parameters can be set for each regional grid in the multiphase flow simulation model so as to participate in subsequent simulation calculation.

[0108] Step 129: Create a calculation domain for the multiphase flow simulation model and set boundary condition parameters of the calculation domain according to operating condition requirements.

[0109] The boundary condition parameters define the boundary conditions of the calculation domain, including the settings of flow rate, pressure, temperature, and other parameters. Reasonable boundary conditions are the key to obtaining accurate simulation results. In actual implementation, accurate settings need to be made according to the actual operating conditions. In this step, the boundary conditions of the inlet and outlet can be specifically set for the boundary layer grid of each calculation domain.

[0110] Optionally, setting the boundary condition parameters of the calculation domain according to operating condition requirements includes: keeping other parameter settings unchanged, using a fluid calculation domain and a porous medium calculation domain for the refrigerant calculation domain respectively, comparing the calculation results of different refrigerant calculation domains, and setting the boundary conditions of the refrigerant.

[0111] In actual implementation, the above steps constitute a complete simulation pre-processing flow. First, a basic geometric model is established through local structure extraction and repair, then a complete flow channel system is constructed through coordinate system definition and model replication, then a calculation domain is established through fluid region extraction and grid division, and finally simulation preparation is completed through physical model definition and parameter setting. Each step is closely linked to ensure the accuracy and calculation efficiency of the simulation model.

[0112] In the embodiments of the present application, through the systematic simulation pre-processing flow, a high-quality Chiller multiphase flow simulation model can be quickly constructed. In this way, the simulation efficiency and accuracy can be significantly improved, so as to provide reliable data support for product design and optimization, and thus shorten the product development cycle.

[0113] In some embodiments, the multiphase flow model can specifically include a mixed multiphase flow model, a two-phase thermal equilibrium model, and an Eulerian multiphase separated flow model.

[0114] The mixed multi-phase flow model (MMP) is a computational model that treats different phase fluids as a single mixture, suitable for flow conditions with strong inter-phase interactions and uniform phase distribution. By solving the momentum and energy equations of the mixture, the model simplifies the calculation process and is suitable for scenarios with rapid phase change of refrigerants. In practical applications, the model can accurately predict the overall flow characteristics of refrigerants during the evaporation process.

[0115] The two-phase thermal equilibrium model (TPTE) assumes that the gas-liquid two-phase is in a state of thermodynamic equilibrium, and simulates the phase change process by solving the phase change equation. This model is suitable for situations where the phase change rate is fast and the inter-phase heat exchange is sufficient, and can accurately predict the rapid evaporation process of refrigerants after the expansion valve. In Chiller simulation, the model can well simulate the gas-liquid two-phase flow and heat exchange of refrigerants.

[0116] The Eulerian multi-phase separated flow model (EMP) treats each phase as an independent continuous medium and solves the conservation equations of each phase separately, suitable for situations where the phase distribution is uneven and there is obvious slip between phases. This model can more accurately capture the dynamic changes of the phase interface and is suitable for simulating the complex interaction between the cooling liquid and the refrigerant. In Chiller simulation, the model can describe in detail the distribution and movement of refrigerant bubbles in the cooling liquid.

[0117] In actual implementation, the three models can be used in combination. For example, the MMP model can be used first for overall flow characteristic analysis, then the TPTE model can be used to simulate the phase change process in specific areas, and finally the EMP model can be used to study the phase separation phenomenon in detail. This combined application method can fully utilize the advantages of each model.

[0118] The three models provide solutions for different flow characteristics: the MMP model is suitable for overall flow analysis, the TPTE model focuses on the phase change process, and the EMP model is good at handling phase separation phenomena. By comparing the calculation results of the three models, the most suitable simulation method for specific Chiller structure and working conditions can be selected.

[0119] In the embodiments of the present application, three different types of multi-phase flow models are provided. In this way, the most suitable model can be selected according to the actual simulation requirements, so that the complex multi-phase flow phenomena in the Chiller can be more accurately simulated, and the reliability and engineering application value of the simulation results can be improved.

[0120] In some embodiments, the grid division model can be a polyhedral grid model, the grid division model of the solid side flow channel plate can be a thin wall layer grid model, and the grid division model of the fluid region can be a boundary layer grid.

[0121] The polyhedral mesh model is a type of mesh composed of polyhedral elements, suitable for the discretization of complex geometries. In Chiller simulation, polyhedral mesh can better adapt to the irregular geometric features of the flow channel plate, improving the calculation accuracy. Compared with traditional tetrahedral mesh, polyhedral mesh has fewer elements and better convergence, especially suitable for areas containing complex flow features. In actual implementation, the system will automatically generate polyhedral mesh according to the geometry, ensuring that the mesh quality meets the calculation requirements.

[0122] The thin-wall layer mesh model is a mesh division method designed specifically for thin-walled structures, which accurately describes the geometric characteristics of thin-walled structures by arranging a small number of layers in the thickness direction. For the solid side flow channel plate in Chiller, thin-wall layer mesh can accurately capture the temperature gradient in the thickness direction while reducing the consumption of computing resources. Thin-wall layer mesh usually uses prismatic elements or hexahedral elements, and setting 3-5 layers of mesh in the thickness direction can meet the engineering precision requirements. This mesh division method is particularly suitable for the simulation of thin-walled structures such as plate heat exchangers.

[0123] The boundary layer mesh model is a high-resolution mesh for the near-wall region of fluid, which accurately captures the velocity gradient and temperature gradient in the boundary layer by arranging dense mesh layers in the wall normal direction. In Chiller simulation, boundary layer mesh is crucial for accurately predicting heat transfer and flow resistance between fluid and wall. This mesh division method can significantly improve the calculation accuracy of near-wall flow and heat transfer.

[0124] In actual implementation, there is a synergistic relationship between the three mesh division models. Polyhedral mesh is used to handle overall complex geometry, thin-wall layer mesh is optimized for the thin-wall characteristics of the solid side flow channel plate, and boundary layer mesh focuses on the near-wall flow details in the fluid region. This combined division strategy can fully leverage the advantages of each mesh type, achieving the best balance between calculation accuracy and efficiency.

[0125] In the embodiments of the present application, a combined division strategy of polyhedral mesh, thin-wall layer mesh and boundary layer mesh is adopted. In this way, the optimal mesh division method can be adopted for the structural characteristics of different regions of Chiller, thereby significantly improving the calculation accuracy and efficiency, and further obtaining more accurate simulation results.

[0126] In some embodiments, setting the mesh parameters can specifically include: controlling other parameter settings unchanged, gradually refining the mesh, comparing the calculation results of different numbers of meshes, and obtaining the mesh parameters when the mesh-independent solution is obtained.

[0127] Grid parameters refer to the various settings used to control the quality and density of mesh generation during CFD simulation. These parameters include base grid size, boundary layer mesh layers, and grid growth rate. The setting of grid parameters directly affects the accuracy and efficiency of the calculation results. For example, the base grid size determines the initial density of the mesh, the boundary layer mesh layers affect the accuracy of capturing near-wall flow, and the grid growth rate controls the smoothness of the transition from dense to sparse mesh. In practical applications, the base grid size can be set to different values such as 0.5mm, 1mm, or 2mm, and the boundary layer mesh layers can be 5, 10, or 15, depending on the simulation accuracy requirements and the limitations of computing resources.

[0128] Grid-independent solution refers to the state where the impact of increasing mesh density on the calculation results can be ignored after the mesh is refined to a certain extent. Achieving grid-independent solution means that the calculation results are accurate enough and no longer affected by mesh density. The method of judging grid-independent solution is usually by comparing the key calculation results (such as pressure drop, temperature distribution, etc.) under different mesh densities. When the difference between the results of two consecutive mesh refinements is less than a pre-set threshold (such as 1%), it is considered that the grid-independent solution has been reached. For example, in Chiller simulation, the changes in refrigerant outlet temperature or pressure drop with mesh refinement can be monitored. When the change rate of these parameters is less than 1%, it can be determined that the current grid parameters are the optimal settings.

[0129] In actual implementation, there is a close relationship between grid parameters and grid-independent solution. By systematically adjusting grid parameters and observing the trend of calculation results, the grid parameter settings required to achieve grid-independent solution can be scientifically determined. This process needs to balance the calculation accuracy and efficiency, avoiding unnecessary consumption of computing resources caused by excessive mesh refinement.

[0130] In the embodiments of the present application, by controlling other parameters unchanged, gradually refining the mesh and comparing the calculation results, the optimal grid parameters can be scientifically determined. In this way, it can ensure that the simulation results meet the accuracy requirements and do not waste computing resources due to excessive mesh refinement, thereby significantly improving the simulation efficiency and providing reliable data support for subsequent optimization design.

[0131] In actual implementation, the grid parameter setting process is an iterative optimization process. First, initial grid parameters need to be set for preliminary calculation, and then whether the grid-independent solution is reached is judged according to the calculation results. If not, adjust the grid parameters (such as reducing the base grid size or increasing the boundary layer mesh layers) for the next round of calculation until the difference between the results of two consecutive calculations is less than the pre-set threshold. This process ensures the reliability of the simulation results while avoiding unnecessary waste of computing resources.

[0132] In the technical solutions of the above embodiments, grid parameter setting, multiphase flow model setting, and refrigerant calculation domain setting are performed respectively, the grid parameters, the multiphase flow model, and the refrigerant calculation domain type can be controlled separately, simulation calculation is performed, and the optimal simulation simulation method is selected by comparing the calculation results with experimental data.

[0133] Figure 5 The flowchart for simulation calculation and post-processing in the multiphase flow simulation method of the cooling equipment provided in the embodiments of the present application, in some embodiments, the function of obtaining simulation results by using the multiphase flow simulation model for simulation calculation and post-processing can be realized by the following operations:

[0134] Step 131: Simulation calculation is performed by using the multiphase flow simulation model.

[0135] The multiphase flow simulation model refers to a calculation model for simulating the flow and heat exchange of the refrigerant and the cooling liquid in the Chiller flow channel. It can be constructed through the foregoing steps.

[0136] In the embodiments of the present application, by using the multiphase flow simulation model for simulation calculation, the flow and heat exchange process of the refrigerant and the cooling liquid in the Chiller can be accurately simulated, thereby providing reliable data support for subsequent optimization design.

[0137] Step 132: According to the post-processing configuration, a cross section of each layer of the flow channel on the refrigerant and cooling liquid side is respectively created along the direction perpendicular to the flow direction in the flow channel, the mass flow on the cross section is obtained, and the flow consistency of the two fluids in the flow channel plate is judged.

[0138] The post-processing configuration refers to parameter setting for processing and analyzing the result data after simulation, including cross section creation position, data extraction type, and secondary analysis processing of result data. The flow channel refers to the channel for the flow of the refrigerant and the cooling liquid in the Chiller, which is stacked by multiple flow channel plates. The mass flow refers to the mass of the fluid passing through a certain cross section per unit time, which is an important parameter for measuring the flow characteristics of the fluid. The flow consistency refers to whether the mass flow of the refrigerant and the cooling liquid in the corresponding flow channel matches, which reflects the balance of heat exchange.

[0139] In the embodiments of the present application, by creating the flow channel cross section and analyzing the mass flow, the flow matching of the two fluids can be intuitively evaluated, thereby judging the heat exchange efficiency and providing a basis for flow channel structure optimization.

[0140] Step 133: Perpendicular cross sections and horizontal cross sections of at least three sampling points are respectively and uniformly distributed along the flow direction in the flow channel plate, the velocity distribution on the perpendicular cross sections and the horizontal cross sections of each sampling point, and the gas-liquid volume fraction distribution of the refrigerant are obtained.

[0141] The sampling points refer to the position points selected according to certain rules in the flow channel for detailed analysis. The vertical section refers to the section perpendicular to the flow channel plate, which is used to analyze the distribution characteristics of the fluid in the thickness direction. The horizontal section refers to the section parallel to the flow channel plate, which is used to analyze the distribution characteristics of the fluid in the plane direction. The velocity distribution refers to the change of the velocity size and direction of the fluid at different positions. The gas-liquid volume fraction distribution refers to the spatial change of the volume proportion of the gas phase and the liquid phase in the refrigerant. The specific positions of the sampling points can be artificially set or selected according to the equal division points of the flow channel plate.

[0142] In the embodiments of the present application, by creating sampling sections in multiple directions and obtaining detailed flow and phase distribution data, the flow characteristics in the flow channel can be comprehensively evaluated, potential flow unevenness or phase separation problems can be found, and then the optimization design of the flow channel structure is guided.

[0143] In actual implementation, there is a close cooperative relationship between the above steps. First, the basic simulation data is obtained through the multiphase flow simulation model, then the preliminary flow analysis is performed through the post-processing configuration, and finally the more detailed flow characteristic data is obtained through the sampling point analysis. The three steps together constitute a complete simulation analysis process, which provides comprehensive data support for the optimization design of the Chiller flow channel structure.

[0144] In some embodiments, before the simulation calculation is performed, the following operations can also be included:

[0145] Step 141: According to the geometric structure model, the number of flow channel plates of the cooling equipment to be simulated, the thickness of each flow channel plate, the distance between adjacent flow channels of the same side fluid, and the distance from the inlet to the first flow channel plate of each side fluid are obtained as post-processing configuration data.

[0146] The number of flow channel plates refers to the number of plate-shaped structures in the cooling equipment for separating the flow channels of the refrigerant and the cooling liquid. The distance between adjacent flow channels of the same side fluid refers to the spacing distance between two adjacent flow channels of the same fluid (refrigerant or cooling liquid), for example, from the top plate to the bottom plate, assuming that the inner cavity between the first layer flow channel plate and the second layer flow channel plate is the cooling liquid flow channel, named the first cooling liquid flow channel, then the inner cavity between the second layer flow channel plate and the third layer flow channel plate is the refrigerant flow channel, the inner cavity between the third layer flow channel plate and the fourth layer flow channel plate is the second cooling liquid flow channel, and the distance between the adjacent flow channels of the same side fluid is the distance between the first cooling liquid flow channel and the second cooling liquid flow channel. The distance from the inlet to the first flow channel plate of each side fluid refers to the path length passed by the fluid from the inlet pipe to the first flow channel plate, for example, the distance between the cooling liquid inlet and the upper surface of the first layer flow channel plate, or the distance between the refrigerant inlet and the upper surface of the first flow channel plate belonging to the refrigerant flow channel. These parameters together constitute the post-processing configuration data, which is used to guide the subsequent simulation result analysis and visualization processing.

[0147] In actual implementation, the system automatically extracts these parameters from the geometric model. For example, the thickness and spacing of the flow channel plate can be obtained by measuring tools, and the distance from the inlet to the first flow channel plate can be determined by path tracking. These data are stored in a specific data structure for subsequent steps.

[0148] Step 142: Post-processing configuration according to the post-processing configuration data.

[0149] Post-processing configuration refers to setting the analysis and display mode of the simulation results according to the extracted structure parameters. This includes determining the cross-section positions to be created, setting visualization parameters, defining analysis indicators, etc. Post-processing configuration directly affects the visualization effect and analysis depth of the simulation results.

[0150] In actual implementation, the system automatically creates a corresponding number of cross-sections according to the number of flow channel plates, sets appropriate display scales according to the thickness of the flow channel plates, and adjusts the visualization density according to the spacing between the flow channels. For example, for a device with 10 layers of flow channel plates, the system will create cross-sections at each layer of flow channel plate position; for flow channel plates with smaller thickness, magnification display will be used to ensure clarity. Typically, the number of flow channel plates can be tens of layers.

[0151] In the embodiments of the present application, the structure parameters required for post-processing are automatically extracted and configured. In this way, the accuracy and consistency of post-processing analysis can be ensured, thereby avoiding errors caused by manual measurement and setting, and thus improving the credibility and analysis efficiency of the simulation results.

[0152] The technical scheme of the embodiments of the present application can be based on the STAR-CCM+ platform for secondary development to generate a plug-in capable of simulating Chiller multi-phase flow, as shown in the following table: Figure 6 The functions of the plug-in can include the following operations, which are implemented by corresponding functional modules:

[0153] S601, check the geometric model of the cooling equipment to be simulated and write in the data of the input working condition parameters, which are used to check whether the imported structure entity (Parts) representing the calculation domain and the naming are accurate, and write in the required simulation data such as input working condition parameters and other configuration parameters needed to control the simulation process and post-processing;

[0154] S602, meshing based on the geometric model and setting meshing parameters, which are used to set the meshing model and set the meshing parameters;

[0155] S603, setting continuum and selecting multi-phase flow model, and setting property parameters, which are used to create a continuum model according to the material properties and assign the calculation domain; continuum refers to medium, refrigerant and coolant, for example, setting the material properties of the medium;

[0156] S604, creating calculation regions and setting boundary conditions according to working conditions, for setting boundary conditions of each calculation region;

[0157] The calculation region refers to a logical unit containing a calculation grid and defining a physical field solving range, and is a core area for CFD simulation calculation. For example, the fluid calculation region is the inner cavity of the cooling liquid or refrigerant, and the solid calculation region is the flow channel plate or pipeline structure.

[0158] S605, creating a report, for creating a report about pressure, temperature, heat exchange, etc.

[0159] S606, setting solving parameters, for setting parameters for changing the solver;

[0160] S607, creating a section, for creating a section required for post-processing, including the section of each layer of the flow channel of the refrigerant and the cooling liquid, and three vertical and horizontal sections created uniformly in the flow direction respectively;

[0161] S608, saving and running, for saving the simulation file and performing running calculation.

[0162] After the simulation calculation is performed, the simulation analysis result can be obtained based on the post-processing configuration and saved.

[0163] The technical scheme of each embodiment of the present application can be implemented through a plug-in, and each functional module in the plug-in can be implemented through script code writing. Specifically, the secondary development script writing can include:

[0164] Writing code for checking the geometric model and parameter input;

[0165] Writing code for selecting a grid model and setting grid parameters;

[0166] Writing code for establishing a continuum and selecting a multiphase flow model and setting physical property parameters;

[0167] Writing code for creating regions and setting boundary condition parameters according to working conditions;

[0168] Writing code for creating a report;

[0169] Writing code for setting solver parameters;

[0170] Writing code for creating a section required for post-processing;

[0171] Writing code for saving a file and running calculation.

[0172] STAR-CCM+ is a powerful CFD simulation platform supporting multi-field modeling and simulation. Secondary development based on STAR-CCM+ refers to the process of customizing the software by developers using programming languages such as Java, C++, and Python under the STAR-CCM+ framework. Developers can write control logic and algorithms to extend platform functionality, improve simulation efficiency, and optimize simulation processes. Through secondary development of Chiller multiphase flow simulation, an automated CFD (Computational Fluid Dynamics) analysis process for Chiller products can be realized. By automating repetitive tasks, reducing manual operation time, and avoiding errors caused by manual input, the simulation speed is accelerated, and work efficiency is improved.

[0173] The beneficial effects of the embodiments of the present application are as follows: For any structure of Chiller product, the first two layers of flow channel plates and other structures are extracted, the structure is appropriately optimized, and the required parameters for Chiller multiphase flow simulation are sorted according to the input working condition information, including grid setting parameters, physical property parameters of refrigerant and coolant, boundary condition parameters, solver parameters, and structure parameters required for post-processing. The geometric model is imported and pre-processing operations are performed, and the subsequent operations and input of simulation are completed using the secondary development plug-in one-key operation, reducing a large amount of manual operation, avoiding errors caused by human input, and quickly, accurately, and efficiently realizing the process of Chiller multiphase flow simulation.

[0174] Figure 7 A structure diagram of a multiphase flow simulation device for a cooling device is provided by the embodiments of the present application. The device includes a pre-processing module 710, a model construction module 720, and a result acquisition module 730.

[0175] The pre-processing module 710 is configured to obtain a geometric structure model of a cooling device to be simulated and input working condition parameters, and perform data preprocessing. The model construction module 720 is configured to perform pre-simulation processing according to preset structural features of the cooling device, the pre-processed geometric structure model, and the input working condition parameters, to construct a multiphase flow simulation model of the cooling device to be simulated. The result acquisition module 730 is configured to perform simulation calculation and post-processing using the multiphase flow simulation model, and obtain a simulation result. The preset structural features of the cooling device include a plurality of layers of flow channel plates stacked between a top plate and a bottom plate, a refrigerant inlet and a refrigerant outlet formed in the top plate or the bottom plate, and a coolant inlet and a coolant outlet formed in the top plate or the bottom plate. The refrigerant inlet passes through a refrigerant flow channel in each layer of flow channel plate to the refrigerant outlet to form a refrigerant flow channel inner cavity. The coolant inlet passes through a coolant flow channel in each layer of flow channel plate to the coolant outlet to form a coolant flow channel inner cavity.

[0176] Optionally, the preprocessing module 710 is specifically configured to:

[0177] According to the functional area annotation information in the geometric structure model, the number of flow channel plates is determined, and the local structure models of the following parts are extracted: the first two layers of flow channel plates, the top plate, the bottom plate, the refrigerant inlet, the refrigerant outlet, the cooling liquid inlet, and the cooling liquid outlet;

[0178] Optimization processing of flow-independent features is performed on each of the local structure models.

[0179] Optionally, the input working condition parameters include: refrigerant, cooling liquid, and cooling device medium; refrigerant inlet flow rate and temperature, refrigerant outlet pressure, pre-expansion valve pressure, and refrigerant inlet and outlet gas-liquid volume fraction; cooling liquid inlet flow rate, temperature, and pressure.

[0180] Optionally, the preprocessing module 710 is further specifically configured to:

[0181] According to the refrigerant, cooling liquid, and cooling device medium, the physical property parameters of the refrigerant, cooling liquid, and cooling device are obtained respectively;

[0182] According to the pre-expansion valve pressure, the refrigerant inlet gas-liquid volume fraction is obtained;

[0183] According to the cooling liquid side inlet flow rate, the cooling liquid side inlet velocity is calculated.

[0184] Optionally, the model construction module 720 is specifically configured to:

[0185] Import the extracted local structure models of the refrigerant and cooling liquid inlets and outlets, the first two layers of flow channel plates, the top plate, and the bottom plate from the preprocessed geometric structure model;

[0186] Surface repair is performed on each of the imported local structure models;

[0187] A Cartesian coordinate system is established at the inlet and outlet of the refrigerant and cooling liquid respectively, wherein the X and Y axes of the Cartesian coordinate system are the inlet and outlet planes, and the Z axis is the flow direction of the flow channel on this side;

[0188] According to the number of flow channel plates, the local structure models of the first two layers of flow channel plates are replicated and spliced along the Z axis direction to construct a simulation internal cavity structure model of the to-be-simulated cooling device;

[0189] According to each local structure model and the simulation internal cavity structure model, the internal cavity structure models of the fluid regions on both sides of the refrigerant and cooling liquid are extracted respectively;

[0190] A mesh division model is selected and mesh parameters are set;

[0191] The physical continuum is created for the internal cavity structure model of the refrigerant and the coolant, and a multiphase flow model is selected as a multiphase flow simulation model of the cooling device to be simulated.

[0192] According to the pretreated input working condition parameters, physical property parameters are set for the multiphase flow simulation model.

[0193] The calculation domain is created for the multiphase flow simulation model, and the boundary condition parameter of the calculation domain is set according to the working condition requirements.

[0194] Optionally, the multiphase flow model includes a mixed multiphase flow model, a two-phase thermal equilibrium model, and an Eulerian multiphase separated flow model.

[0195] Optionally, setting the boundary condition parameter of the calculation domain according to the working condition requirements includes:

[0196] While controlling other parameter settings unchanged, the fluid calculation domain and the porous medium calculation domain are used for the refrigerant calculation domain respectively, and the calculation results of different refrigerant calculation domains are compared to set the refrigerant boundary condition.

[0197] Optionally, the grid division model is a polyhedral grid model, the grid division model of the solid side flow channel plate is a thin wall layer grid model, and the grid division model of the fluid region is a boundary layer grid.

[0198] Optionally, setting the grid parameter includes:

[0199] While controlling other parameter settings unchanged, the grid is gradually refined, the calculation results of different numbers of grids are compared, and the grid parameter when the grid independent solution is obtained.

[0200] Optionally, the result acquisition module 730 is specifically configured to:

[0201] The multiphase flow simulation model is used for simulation calculation;

[0202] According to the post-processing configuration, the cross sections of each layer of the flow channel on the refrigerant side and the coolant side are respectively created along the direction perpendicular to the flow direction in the flow channel, the mass flow on the cross section is obtained, and the flow consistency of the two sides of the fluid in the flow channel plate is judged.

[0203] Along the flow direction in the flow channel plate, the vertical cross sections and the horizontal cross sections of at least three sampling points are respectively and uniformly distributed, the velocity distribution on the vertical cross section and the horizontal cross section of each sampling point, and the gas-liquid fraction distribution of the refrigerant are obtained.

[0204] Optionally, the device further includes:

[0205] The post-processing configuration module is configured to, before simulation calculation and post-processing are performed by using the multiphase flow simulation model, acquire, according to the geometric structure model, the number of flow channel plates of the cooling device to be simulated, the thickness of each flow channel plate, the distance between adjacent flow channels of the same side fluid, and the distance from the inlet to the first flow channel plate of each side fluid as post-processing configuration data; and perform post-processing configuration according to the post-processing configuration data.

[0206] The multiphase flow simulation device of the cooling device can also perform the multiphase flow simulation method of the cooling device provided by any embodiment of the present application, and has the corresponding function modules and beneficial effects of performing the method.

[0207] Figure 8 A structural schematic diagram of an electronic device that can be used to implement embodiments of the present application is shown. The electronic device is intended to represent various forms of digital computers, such as laptops, desktops, tablets, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers. The electronic device can also represent various forms of mobile devices, such as personal digital processors, cellular telephones, smart phones, wearable devices (e.g., headsets, glasses, watches, etc.), and other similar computing devices. The components shown here, their connections, and their functions, as described and demonstrated in this document, are meant to be examples only, and are not intended to limit the implementations of the present application described and / or claimed in this document.

[0208] As shown in the figure, the electronic device 80 includes at least one processor 81, and a memory, such as a read-only memory (ROM) 82, a random access memory (RAM) 83, etc., which is communicatively connected to the at least one processor 81, wherein the memory stores a computer program that can be executed by the at least one processor. The processor 81 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 82 or the computer program loaded from the storage unit 88 into the random access memory (RAM) 83. In the RAM 83, various programs and data required for the operation of the electronic device 80 can also be stored. The processor 81, the ROM 82, and the RAM 83 are connected to each other through a bus 84. An input / output (I / O) interface 85 is also connected to the bus 84.

[0209] A plurality of components in the electronic device 80 are connected to the I / O interface 85, including: an input unit 86, such as a keyboard, a mouse, etc.; an output unit 87, such as various types of displays, speakers, etc.; a storage unit 88, such as a magnetic disk, an optical disk, etc.; and a communication unit 89, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 89 allows the electronic device 80 to exchange information / data with other devices through a computer network, such as the Internet, and / or various telecommunications networks.

[0210] The processor 81 can be various general-purpose and / or special-purpose processing components having processing and computing capabilities. Some examples of the processor 81 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, and the like. The processor 81 performs various methods and processes described above, such as the multiphase flow simulation method of a cooling device.

[0211] In some embodiments, the multiphase flow simulation method of a cooling device can be implemented as a computer program tangibly embodied in a computer readable storage medium, such as the storage unit 88. In some embodiments, part or all of the computer program can be loaded and / or installed onto the electronic device 80 via the ROM 82 and / or the communication unit 89. When the computer program is loaded onto the RAM 83 and executed by the processor 81, one or more steps of the multiphase flow simulation method of a cooling device described above can be performed. Alternatively, in other embodiments, the processor 81 can be configured to perform the multiphase flow simulation method of a cooling device by any other suitable means, such as by means of firmware.

[0212] Various implementations of the systems and techniques described above can be realized in digital electronic circuitry, integrated circuitry, a field programmable gate array (FPGA), an application specific integrated circuit (ASIC), a system on a chip (SOC), a programmable logic device (CPLD), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include implementation in one or more computer programs that are executable and / or interpretable on a programmable system including at least one programmable processor, which can be special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device.

[0213] Computer programs used to implement the methods of the application can be written in any combination of one or more programming languages. These computer programs can be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the computer program, when executed by the processor of the machine, implements the functions / acts specified in the flowcharts and / or block diagrams. The computer program can be executed entirely on a machine, partially on a machine, partially on a machine and partially on a remote machine or entirely on a remote machine or server.

[0214] In the context of this invention, a computer-readable storage medium can be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, apparatus, or device. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0215] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0216] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or middleware components (e.g., application servers), or frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0217] The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system, to solve the defects of large management difficulty and weak business scalability in traditional physical host and VPS service.

[0218] It should be understood that the various forms of flow shown above can be reordered, added to, or have steps deleted. For example, the steps described in the present application can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solutions of the present application can be achieved, which are not limited herein.

[0219] The above detailed description does not constitute a limitation on the protection scope of the present application. Those skilled in the art should understand that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A multiphase flow simulation method for a cooling device, characterized in that, include: Obtain the geometric structure model and input operating parameters of the cooling equipment to be simulated, and perform data preprocessing; Based on the preset structural features of the cooling equipment, the pre-processed geometric model, and the input operating parameters, pre-simulation processing is performed to construct a multiphase flow simulation model of the cooling equipment to be simulated. The multiphase flow simulation model was used for simulation calculations and post-processing to obtain simulation results. The cooling device includes a pre-defined structural feature comprising multiple layers of flow channel plates stacked between a top plate and a bottom plate; a refrigerant inlet and a refrigerant outlet on the top plate or bottom plate, and a coolant inlet and a coolant outlet on the top plate or bottom plate; the refrigerant inlet flows through the refrigerant channels in each layer of flow channel plates to the refrigerant outlet to form a refrigerant flow channel cavity; the coolant inlet flows through the coolant channels in each layer of flow channel plates to the coolant outlet to form a coolant flow channel cavity.

2. The method according to claim 1, characterized in that, Obtaining the geometric model of the cooling equipment to be simulated and performing data preprocessing includes: Based on the functional area labeling information in the geometric structure model, determine the number of flow channel plates and extract the local structural models of the following parts: the flow channel plates of the first two layers, the top plate, the bottom plate, the refrigerant inlet, the refrigerant outlet, the coolant inlet, and the coolant outlet; The flow-independent features of each of the local structural models are optimized.

3. The method according to claim 1, characterized in that, The input operating parameters include: refrigerant, coolant and cooling equipment medium; refrigerant inlet flow rate and temperature, refrigerant outlet pressure, expansion valve pressure and refrigerant inlet and outlet gas-liquid fractions; coolant inlet flow rate, temperature and pressure.

4. The method according to claim 3, characterized in that, Obtaining the input operating parameters of the cooling equipment to be simulated and performing data preprocessing includes: Based on the refrigerant, coolant, and cooling equipment medium, obtain the physical property parameters of the refrigerant, coolant, and cooling equipment respectively; The gas-liquid integral of the refrigerant at the inlet is obtained based on the pressure before the expansion valve. Calculate the inlet velocity on the coolant side based on the inlet flow rate on the coolant side.

5. The method according to any one of claims 1-4, characterized in that, Based on the preset structural features of the cooling equipment, the preprocessed geometric model, and the input operating parameters, the simulation preprocessing includes: Import the extracted refrigerant and coolant inlets and outlets, the first two flow channel plates, and the local structural models of the top and bottom plates from the preprocessed geometric model. Perform surface repair on each imported local structural model; Cartesian coordinate systems are established at the inlet and outlet of the refrigerant and coolant, respectively, where the X-axis and Y-axis of the Cartesian coordinate system are the inlet plane and outlet plane, and the Z-axis is the flow direction of the flow channel on that side. Along the Z-axis, based on the number of flow channel plates, the local structural models of the first two flow channel plates are copied and spliced ​​together to construct the simulated internal cavity structural model of the cooling equipment to be simulated. Based on the local structural models and the simulated internal cavity structural models, the internal cavity structural models of the fluid regions on both sides of the refrigerant and coolant are extracted respectively. Select the mesh generation model and set the mesh parameters; A physical continuum is created for the internal cavity structure model of the refrigerant and coolant, and a multiphase flow model is selected as the multiphase flow simulation model for the cooling device to be simulated. Based on the preprocessed input operating condition parameters, set the physical property parameters for the multiphase flow simulation model; A computational domain is created for the multiphase flow simulation model, and the boundary condition parameters of the computational domain are set according to the operating conditions.

6. The method according to claim 5, characterized in that, The multiphase flow models include the mixed multiphase flow model, the two-phase thermal equilibrium model, and the Eulerian multiphase separation flow model.

7. The method according to claim 5, characterized in that, Setting the boundary condition parameters of the computational domain according to the operating conditions includes: Keeping other parameter settings unchanged, the calculation results of the fluid calculation domain and the porous medium calculation domain are compared to set the refrigerant boundary conditions.

8. The method according to claim 5, characterized in that, The meshing model is a polyhedral mesh model, the meshing model of the solid side channel plate is a thin-walled layer mesh model, and the meshing model of the fluid region is a boundary layer mesh.

9. The method according to claim 5, characterized in that, Setting grid parameters includes: By keeping other parameter settings constant, the mesh is gradually refined, and the calculation results of different numbers of meshes are compared to obtain the mesh parameters when the solution is mesh-independent.

10. The method according to any one of claims 1-4, characterized in that, The simulation results obtained by using the multiphase flow simulation model for simulation calculation and post-processing include: The multiphase flow simulation model was used for simulation calculations. Based on the post-processing configuration, cross sections of each channel on the refrigerant and coolant sides are created along the flow direction perpendicular to the channel, the mass flow rate on the cross section is obtained, and the consistency of the flow rate of the fluid on both sides in the channel plate is determined. Along the flow direction inside the flow channel plate, at least three sampling points are uniformly distributed to create vertical and horizontal cross sections to obtain the velocity distribution on the vertical and horizontal cross sections of each sampling point, as well as the gas-liquid ratio distribution of the refrigerant.

11. The method according to claim 10, characterized in that, Before performing simulation calculations and post-processing using the multiphase flow simulation model, the following steps are also included: Based on the geometric structure model, the number of flow channel plates, the thickness of each flow channel plate, the distance between adjacent flow channels on the same side, and the distance between the inlet and the first flow channel plate of the fluid on each side are obtained as post-processing configuration data for the cooling equipment to be simulated. Perform post-processing configuration based on the post-processing configuration data.

12. A multiphase flow simulation device for a cooling equipment, characterized in that, include: The preprocessing module is used to obtain the geometric structure model and input operating parameters of the cooling equipment to be simulated, and to perform data preprocessing. The model building module is used to perform pre-simulation processing based on the preset structural features of the cooling equipment, the pre-processed geometric structure model, and the input operating parameters, so as to build a multiphase flow simulation model of the cooling equipment to be simulated. The result acquisition module is used to perform simulation calculations and post-processing using the multiphase flow simulation model to obtain simulation results. The cooling device includes a pre-defined structural feature comprising multiple layers of flow channel plates stacked between a top plate and a bottom plate; a refrigerant inlet and a refrigerant outlet on the top plate or bottom plate, and a coolant inlet and a coolant outlet on the top plate or bottom plate; the refrigerant inlet flows through the refrigerant channels in each layer of flow channel plates to the refrigerant outlet to form a refrigerant flow channel cavity; the coolant inlet flows through the coolant channels in each layer of flow channel plates to the coolant outlet to form a coolant flow channel cavity.

13. An electronic device, characterized in that, The electronic device includes: At least one processor; and a memory communicatively connected to the at least one processor; The memory stores a computer program that can be executed by the at least one processor, which is then executed by the at least one processor to enable the at least one processor to perform the multiphase flow simulation method for the cooling device according to any one of claims 1-11.

14. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions that cause a processor to execute the multiphase flow simulation method for the cooling device according to any one of claims 1-11.

Images