A calculation method, system and computer device for fluid-solid coupling heat transfer

Abstract

The application provides a kind of calculation method, system and computer equipment for fluid-structure coupling heat transfer.For the problem that existing weak coupling algorithm long tail iteration converges slowly in later stage, a two-stage decoupling acceleration strategy is proposed.In the first stage, fluid-solid weak coupling algorithm is iterated, and an early truncation mechanism is introduced: when the global temperature rise has not finally converged, but the convective heat transfer coefficient of the interface reaches a stable convergence state first, stop the fluid-structure two-way iteration and extract the current convective heat transfer coefficient;In the second stage, the convective heat transfer coefficient is applied as a thermal boundary condition to the solid calculation domain, and the flow equation operation is stopped, and the heat conduction equation of the solid is solved in a one-way decoupling state until complete steady state.The application skillfully utilizes the space-time difference of the convergence speed of the characteristic physical quantity, eliminates the redundant iteration without losing the engineering precision, and greatly improves the solving efficiency of heat transfer simulation.

Description

Technical Field

[0001] This invention belongs to the field of computer-aided engineering and computational fluid dynamics numerical simulation, and more specifically, relates to a calculation method, system and computer equipment for fluid-structure interaction heat transfer. Background Technology

[0002] Currently, fluid-solid coupling heat transfer problems are widespread in engineering practice, such as thermal management in the automotive industry, chip heat dissipation in electronic devices, thermal-hydraulic analysis of nuclear reactors, and petrochemical fields. During product research and development, quickly and accurately predicting the heat exchange state in these scenarios is crucial for achieving precise adjustment and control of equipment operating temperature, improving operational efficiency, and extending equipment lifespan.

[0003] In the design evaluation phase, commonly used methods for obtaining temperature field distribution mainly include experimental testing and numerical calculation. Among these, numerical calculation is more widely used in industrial production because it has the advantages of higher iterative efficiency, lower cost, and the ability to obtain full-field data compared to experimental testing.

[0004] Existing numerical methods for fluid-structure interaction heat transfer include two types of algorithms: strongly coupled (implicit update) and weakly coupled (explicit update). Strongly coupled algorithms offer high accuracy and can provide a more precise picture of transient temperature changes, requiring fewer iterations to reach convergence. However, this algorithm couples the equations of the fluid and solid domains into a single large matrix, resulting in extremely complex matrix construction, massive computational cost per step, and numerous limitations in highly compressible flows. Weakly coupled algorithms, on the other hand, have a simpler program framework. The equations of the fluid and solid domains are solved separately, with the boundary temperature updated only once at the fluid-solid interface at the end of each time step. Due to explicit decoupling, their single-step computation speed is faster. However, the disadvantage of weakly coupled algorithms lies in their extremely slow overall convergence speed, often limiting their application to steady-state results. Furthermore, in the later stages of computation, the system requires a massive number of time steps to bridge small residuals to achieve complete convergence.

[0005] When faced with complex problems in modern industrial design, such as high Reynolds number flows and cross-scale heat transfer characteristics, existing strongly or weakly coupled algorithms throughout the entire process consume extremely large amounts of computing resources. Current computing power is often insufficient to meet the needs of rapid iterative product design. Therefore, developing a fluid-solid heat transfer algorithm that balances high efficiency and high accuracy is particularly urgent. Summary of the Invention

[0006] In view of this, the present invention provides a calculation method, system and computer device for fluid-structure interaction heat transfer.

[0007] According to a first aspect of the present invention, a calculation method for fluid-structure interaction heat transfer is provided, the method comprising the following steps: Step S1: Obtain the fluid-structure interaction calculation model that includes the fluid calculation domain and the solid calculation domain, and set the initial velocity field, pressure field and temperature field; Step S2, First stage solution: Within consecutive iteration time steps, the fluid-solid weak coupling heat transfer algorithm is used to solve the fluid control equations of the fluid computational domain and the solid control equations of the solid computational domain respectively, and the temperature field information is updated interactively at the interface between the fluid computational domain and the solid computational domain. Step S3: During the iterative process of solving the first stage, the preset convergence index is monitored in real time; when the convergence index meets the first preset convergence condition, the first stage solution is stopped, and the convective heat transfer coefficient at the current interface is extracted. The first preset convergence condition indicates that the global temperature field of the fluid-structure interaction calculation model has not reached the final fully steady-state convergence state, and the convective heat transfer coefficient has reached the stable convergence state. Step S4, Second stage solution: The extracted convective heat transfer coefficient is applied as a thermodynamic boundary condition to the target computational domain, and the target computational domain is solved by one-way decoupling until the final fully converged state is reached to obtain the steady-state temperature field distribution. The unidirectional decoupling solution refers to the fact that the bidirectional interaction of the temperature field between the fluid computational domain and the solid computational domain is no longer performed during the solution process.

[0008] Alternatively, the target computational domain may be the solid-state computational domain; In step S4, the convective heat transfer coefficient is used as a third type of convective heat transfer boundary condition and applied only to the solid computational domain. During the unidirectional decoupling solution process, the update of the fluid computational domain is frozen and the solution of the fluid control equations of the fluid computational domain is stopped. The solid control equations of the solid computational domain are solved only to obtain the steady-state temperature field distribution of the solid computational domain.

[0009] Alternatively, the fluid governing equation is the incompressible Navier-Stokes equation, and the solid governing equation is the solid heat conduction equation.

[0010] Optionally, in step S3, the first preset convergence condition includes at least one of the following: Condition A: The normalized gradient of the temperature field at the interface or within the computational domain with respect to the iteration time step falls back to the range of 0.01 to 0.001; Condition B: Assume the current iteration time step is The historical regional average value of the global temperature field ave ( T )satisfy And the convective heat transfer coefficient h Gradient of residual variation with iteration time step satisfy .

[0011] Alternatively, the convective heat transfer coefficient h The calculation formula is:

[0012] In the above formula, Q The total convective heat transfer through the interface, A Let be the area of ​​the interface. T w The wall temperature on the solid side of the interface is [temperature value missing]. T ref The reference temperature is the fluid side of the interface.

[0013] Alternatively, the specific expression of the solid heat conduction equation is as follows:

[0014] In the above formula, For solid density, For isobaric specific heat capacity, T For temperature, t For time, Thermal conductivity, For temperature gradient, The divergence of heat flux density, The intensity of the heat source within the volume.

[0015] Alternatively, the specific expression of the incompressible Navier-Stokes equation is as follows:

[0016]

[0017]

[0018] In the above formula, velocity vector For velocity divergence, the volume expansion rate of fluid particles is zero under incompressible conditions; For fluid density, For local acceleration, For convective acceleration, For pressure gradient, For dynamic viscosity, For the velocity Laplace operator, It is a volume force vector.

[0019] Alternatively, both the incompressible Navier-Stokes equation and the solid heat conduction equation can be solved spatially using the finite volume method.

[0020] According to a second aspect of the present invention, a computational system for fluid-structure interaction heat transfer is provided, the system comprising the following functional modules: The model initialization module is used to obtain a fluid-structure interaction calculation model that includes both fluid and solid computational domains, and to set the initial velocity, pressure, and temperature fields. The first-stage solution module is used to solve the fluid control equations of the fluid computational domain and the solid control equations of the solid computational domain respectively within consecutive iteration time steps using a fluid-solid weak coupling heat transfer algorithm, and interactively update the temperature field information at the interface between the fluid computational domain and the solid computational domain. The monitoring and extraction module is used to monitor the preset convergence index in real time during the iterative process of the first stage solution; when the convergence index meets the first preset convergence condition, the first stage solution is stopped and the convective heat transfer coefficient at the current interface is extracted; wherein, the first preset convergence condition indicates that the global temperature field of the fluid-structure interaction calculation model has not reached the final fully steady-state convergence state, and the convective heat transfer coefficient has reached the stable convergence state. The second-stage solution module is used to apply the extracted convective heat transfer coefficient as a thermodynamic boundary condition to the target computational domain, and to perform a one-way decoupled solution on the target computational domain until the final fully converged state is reached to obtain the steady-state temperature field distribution; wherein, the one-way decoupled solution means that the bidirectional interaction of the temperature field between the fluid computational domain and the solid computational domain is no longer performed during the solution process.

[0021] According to a third aspect of the present invention, a computer device is provided, comprising a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of any of the above-described calculation methods for fluid-structure interaction heat transfer.

[0022] The beneficial effects of this invention are as follows: The computational method for fluid-structure interaction (FSI) heat transfer provided by this invention creatively employs a two-stage decoupling strategy, based on the physical characteristic that the convective heat transfer coefficient (HTC) converges earlier than the absolute temperature field in FSI systems. By introducing a first preset convergence condition to prematurely truncate the time-consuming weak coupling process, and then transitioning to pure solid-state heat conduction calculations based on unidirectional HTC boundary conditions in the second stage, the method significantly reduces the number of invalid iterations in the Navier-Stokes equations while maintaining the accuracy of the final temperature distribution prediction, thereby reducing the overall simulation time by several times. Using the computational method for FSI heat transfer provided by this invention can effectively accelerate the iteration cycle of related engineering product development and has significant engineering application value.

[0023] Other features and advantages of the present invention will be described in detail in the following detailed description section. Attached Figure Description

[0024] The present invention can be better understood by referring to the following description taken in conjunction with the accompanying drawings, in which the same or similar reference numerals are used throughout the drawings to denote the same or similar parts.

[0025] Figure 1 A flowchart illustrating the implementation of a computational method for fluid-structure interaction heat transfer according to an embodiment of the present invention is shown. Figure 2 A schematic diagram of the overall three-dimensional configuration of a thermal management test example for an electronic device according to an embodiment of the present invention is shown. Figure 3 The steady-state temperature distribution contour map obtained by using the traditional weakly coupled algorithm is shown. Figure 4 The steady-state temperature distribution cloud map obtained by using the two-stage fast algorithm of this invention is shown. Figure 5 A schematic diagram comparing the temperature convergence curve using a traditional weakly coupled algorithm with the stage cutoff point of an embodiment of the present invention is shown.

[0026] Figure 6 A structural block diagram of a computational system for fluid-structure interaction heat transfer according to an embodiment of the present invention is shown. Detailed Implementation

[0027] To enable those skilled in the art to more fully understand the technical solutions of the present invention, exemplary embodiments of the present invention will be described more comprehensively and in detail below with reference to the accompanying drawings. Obviously, the one or more embodiments of the present invention described below are merely one or more specific ways to implement the technical solutions of the present invention, and are not exhaustive. It should be understood that other ways belonging to a general inventive concept can be used to implement the technical solutions of the present invention, and should not be limited to the embodiments described exemplary. Based on one or more embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0028] Example: Addressing the issues of extremely slow convergence speed in the later stages of calculations in existing weakly coupled fluid-structure interaction (FSI) methods, and the complexity and computational limitations of strongly coupled methods, this invention provides an efficient two-stage calculation method for fluid-solid coupled heat transfer while maintaining computational accuracy. This invention aims to significantly accelerate the calculation speed of fluid-structure interaction heat transfer-related designs by eliminating long-tail iterations through a decoupling algorithm.

[0029] It should be noted that the coupling in the embodiments of the present invention refers to the information transfer and solution of a single physical field of temperature field between multiple computational domains in numerical calculation, and does not involve the coupling of other multi-physical fields such as electromagnetic field and stress field.

[0030] Figure 1 A flowchart illustrating the implementation of the calculation method for fluid-structure interaction heat transfer according to an embodiment of the present invention is shown. (Refer to...) Figure 1 The calculation method for fluid-structure interaction heat transfer according to embodiments of the present invention includes the following steps: Step S1: Obtain the fluid-structure interaction calculation model that includes the fluid calculation domain and the solid calculation domain, and set the initial velocity field, pressure field and temperature field; Step S2, First stage solution: Within consecutive iteration time steps, the fluid-solid weak coupling heat transfer algorithm is used to solve the fluid control equations in the fluid computational domain and the solid control equations in the solid computational domain, and the temperature field information is updated interactively at the interface between the fluid computational domain and the solid computational domain. Step S3: During the iterative process of solving the first stage, monitor the preset convergence index in real time; when the convergence index meets the first preset convergence condition, stop solving the first stage and extract the convective heat transfer coefficient at the current interface. Among them, the first preset convergence condition indicates that the global temperature field of the fluid-structure interaction calculation model has not reached the final fully steady-state convergence state, and the convective heat transfer coefficient has reached the stable convergence state. Step S4, Second stage solution: The extracted convective heat transfer coefficient is applied as a thermodynamic boundary condition to the target computational domain. The target computational domain is solved by unidirectional decoupling until the final fully converged state is reached, and the steady-state temperature field distribution is obtained. Among them, unidirectional decoupling solution refers to the elimination of bidirectional temperature field interaction between the fluid computational domain and the solid computational domain during the solution process.

[0031] Furthermore, in this embodiment of the invention, the target computational domain is a solid-state computational domain; In step S4, the convective heat transfer coefficient is used as the third type of convective heat transfer boundary condition and is applied only to the solid computational domain. During the unidirectional decoupling solution process, the update of the fluid computational domain is frozen and the solution of the fluid control equations of the fluid computational domain is stopped. The solid control equations of the solid computational domain are solved only to obtain the steady-state temperature field distribution of the solid computational domain.

[0032] Furthermore, in this embodiment of the invention, the fluid governing equation is the incompressible Navier-Stokes equation, and the solid governing equation is the solid heat conduction equation.

[0033] Furthermore, in step S3 of this embodiment of the invention, the first preset convergence condition includes at least one of the following: Condition A: The normalized gradient of the temperature field at the interface or within the computational domain with respect to the iteration time step falls back to the range of 0.01 to 0.001; Condition B: Assume the current iteration time step is Historical regional average of the global temperature field ave ( T )satisfy And the convective heat transfer coefficient h Gradient of residual variation with iteration time step satisfy .

[0034] Furthermore, in this embodiment of the invention, the convective heat transfer coefficient h The calculation formula is:

[0035] In the above formula, Q The total convective heat transfer through the interface, A Let be the area of ​​the interface. T w The wall temperature on the solid side of the interface is [temperature value missing]. T ref The reference temperature is the fluid side of the interface.

[0036] Furthermore, in this embodiment of the invention, the specific expression of the solid heat conduction equation is as follows:

[0037] In the above formula, For solid density, For isobaric specific heat capacity, T For temperature, t For time, Thermal conductivity, For temperature gradient, The divergence of heat flux density, The intensity of the heat source within the volume.

[0038] Furthermore, in this embodiment of the invention, the specific expression of the incompressible Navier-Stokes equation is as follows:

[0039]

[0040]

[0041] In the above formula, velocity vector For velocity divergence, the volume expansion rate of fluid particles is zero under incompressible conditions; For fluid density, For local acceleration, For convective acceleration, For pressure gradient, For dynamic viscosity, For the velocity Laplace operator, It is a volume force vector.

[0042] Furthermore, in this embodiment of the invention, both the incompressible Navier-Stokes equation and the solid heat conduction equation are solved spatially using the finite volume method.

[0043] Specifically, the underlying physical and mathematical mechanism by which the calculation method for fluid-structure interaction heat transfer in this invention can significantly improve computational efficiency lies in: In the early stages of fluid-solid weakly coupled heat transfer calculations, the flow field is established quickly, the system temperature gradient is large, and the temperature rise is extremely rapid. The temperature rise in this stage can account for 80% or more of the total temperature rise, and the time steps consumed typically account for only 40% or less of the total time steps. However, in the later stages of the calculation, the system needs to consume several times more time steps to slowly compensate for small temperature differences; this is known as inefficient long-tail iteration. The embodiments of this invention directly truncate this inefficient iteration.

[0044] In fluid-solid convective heat transfer, the formula for the convective heat transfer coefficient (HTC) is: Since the heat flow characteristics in the numerator and the temperature difference characteristics in the denominator are both dynamically determined by the local gradient of the fluid boundary layer, they develop synchronously. This results in the convergence rate of the physical parameter $h$ being much faster than the convergence rate of the absolute temperature field itself. In the first stage before complete convergence (temperature rise of about 80%), the heat transfer capacity of the interface has essentially approached the true final steady-state constant.

[0045] After HTC converges, it is used to replace the original nonlinear fluid-structure interaction boundary, and is applied unidirectionally as a constant boundary condition to the solid region for solution. This is equivalent to reducing the complex non-fluid coupled system to a simple pure heat conduction solution system. Therefore, the calculation can overcome obstacles and achieve complete convergence instantly.

[0046] The following is a more detailed explanation of the calculation method for fluid-structure interaction heat transfer in this embodiment of the invention, based on a specific example.

[0047] like Figure 2 As shown, the embodiments of the present invention use a typical electronic device heat dissipation model for verification.

[0048] The model's 3D geometry comprises three main components: the external chassis 1, the internal PCB 2, and the core heat-generating component, CPU 3. After meshing, the air channels within chassis 1 constitute the fluid computation domain, while PCB 2 and CPU 3 form the solid computation domain. Boundary conditions are set as follows: the left Xmin surface of the chassis is designated as the velocity inlet with a given inlet velocity of 1 m / s; the right Xmax surface of the chassis is designated as the pressure outlet; and the remaining external surfaces of chassis 1 are all designated as non-slip solid walls. A heat dissipation source is added inside CPU 3, with a total heat output set to 2.8 W.

[0049] The specific implementation flow of the calculation method for fluid-structure interaction heat transfer in this embodiment of the invention is as follows.

[0050] Phase 1: Truncation-based weakly coupled computation.

[0051] Set the initial flow field and temperature field of the system, and enter the time step iteration loop.

[0052] 1) First, perform incompressible Navier-Stokes equations calculations within the fluid computational domain. Solve the momentum prediction equation to obtain the temporary velocity field; construct and solve the pressure correction equation to obtain the pressure correction value; update the temporary velocity field using the pressure correction value, and update the pressure field simultaneously.

[0053] 2) Using the updated velocity and pressure fields, solve the fluid energy equation in the fluid domain; at the same time, solve the solid heat conduction equation in the solid computational domain.

[0054] 3) After each single-step solution is completed, the temperature information at the fluid-solid interface is updated interactively. Throughout the above loop, the program background continuously monitors the residual change in the convective heat transfer coefficient (HTC) at the interface and the overall system temperature gradient. When the program reaches the 105th iteration time step (e.g., ...), ... Figure 5 When the position is marked by the vertical dashed line in the middle, the system determines that the physical quantity at this time satisfies the first preset convergence condition: At this point, the normalized temperature gradient has dropped to around 0.01, the temperature rise has reached approximately 80% of the steady-state total temperature rise, and it meets the criteria for the historical average temperature difference. More importantly, the gradient residual of HTC extracted from the interface at this time satisfies... This indicates that the physical quantity of the heat transfer coefficient has been completely stabilized.

[0055] The system then issues a stop command, forcibly interrupting the time-consuming weakly coupled fluid-structure bidirectional iteration, and saves the current high-precision HTC data field.

[0056] Second stage: One-way decoupling and rapid convergence.

[0057] The system extracts the HTC field and reference temperature data saved at the end of the first stage.

[0058] The HTC is applied as a third-order convective thermal boundary condition corresponding to Newton's law of cooling, on the outer surface of the solid-state computing domain (PCB and CPU). At this point, the system completely freezes the fluid computing domain (stops solving all fluid flow equations and ceases data transfer), solving only the pure solid-state steady-state heat conduction equations within the isolated solid-state computing domain. Since there is no need to wait for data matching with the flow field, the temperature in the solid region rises instantaneously within a few tiny computational steps and locks into the final steady-state convergence state. Figure 5 The part to the right of the dotted line is completed instantly in this stage, and the time taken is negligible.

[0059] Accuracy comparison of verification results: Figure 5 The convergence process of the traditional weakly coupled method (solid line curve) is shown, which requires about 600 iterations to achieve absolute convergence through slow, long-tail approximation. In contrast, the algorithm of this invention completes the core flow field calculation in only about 105 steps and then converges instantly, reducing the overall computation time by about 80%.

[0060] Figure 3 The baseline steady-state temperature distribution cloud map is obtained after 600 steps by the traditional weakly coupled algorithm, with the highest temperature point being approximately 391.52 K (118.37 ℃). Figure 4 This is a steady-state temperature distribution cloud map obtained using the algorithm of this embodiment of the invention. (Comparison) Figure 3 and Figure 4It is evident that the overall thermal diffusion pattern and temperature gradient distribution of the two cloud maps maintain a high degree of consistency. Precise numerical comparison reveals that the maximum temperature difference between the core maximum temperature obtained by this algorithm and the traditional benchmark algorithm is only 0.81 K. This result fully demonstrates that the calculation method for fluid-structure interaction heat transfer proposed in this embodiment of the invention significantly reduces simulation calculation time while perfectly ensuring industrial-grade accuracy, making it a highly efficient and practical engineering thermal simulation calculation technique.

[0061] Accordingly, based on the calculation method for fluid-structure interaction heat transfer in the embodiments of the present invention, the embodiments of the present invention also propose a calculation system for fluid-structure interaction heat transfer.

[0062] Figure 6 A structural block diagram of a computational system for fluid-structure interaction heat transfer according to an embodiment of the present invention is shown. (Refer to...) Figure 6 The computational system for fluid-structure interaction heat transfer according to embodiments of the present invention includes the following functional modules: The model initialization module is used to obtain a fluid-structure interaction calculation model that includes both fluid and solid computational domains, and to set the initial velocity, pressure, and temperature fields. The first-stage solution module is used to solve the fluid control equations in the fluid computational domain and the solid control equations in the solid computational domain using a fluid-solid weak coupling heat transfer algorithm within consecutive iteration time steps, and to interactively update the temperature field information at the interface between the fluid computational domain and the solid computational domain. The monitoring and extraction module is used to monitor the preset convergence index in real time during the iteration process of the first stage solution. When the convergence index meets the first preset convergence condition, the first stage solution is stopped and the convective heat transfer coefficient at the current interface is extracted. The first preset convergence condition indicates that the global temperature field of the fluid-structure interaction calculation model has not reached the final fully steady-state convergence state, and the convective heat transfer coefficient has reached the stable convergence state. The second-stage solution module is used to apply the extracted convective heat transfer coefficient as a thermodynamic boundary condition to the target computational domain, and to perform unidirectional decoupling solution on the target computational domain until the final fully converged state is reached to obtain the steady-state temperature field distribution; wherein, unidirectional decoupling solution means that there is no longer bidirectional interaction of the temperature field between the fluid computational domain and the solid computational domain during the solution process.

[0063] Accordingly, based on the calculation method for fluid-structure interaction heat transfer in the embodiments of the present invention, the embodiments of the present invention also propose a computer device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the calculation method for fluid-structure interaction heat transfer as proposed in the embodiments of the present invention.

[0064] While one or more embodiments of the present invention have been described above, those skilled in the art will recognize that the present invention can be implemented in any other form without departing from its spirit and scope. Therefore, the embodiments described above are illustrative and not restrictive, and many modifications and substitutions will be apparent to those skilled in the art without departing from the spirit and scope of the invention as defined in the appended claims.

Claims

1. A calculation method for heat transfer in fluid-structure interaction, characterized in that, Includes the following steps: Step S1: Obtain the fluid-structure interaction calculation model that includes the fluid calculation domain and the solid calculation domain, and set the initial velocity field, pressure field and temperature field; Step S2, First stage solution: Within consecutive iteration time steps, the fluid-solid weak coupling heat transfer algorithm is used to solve the fluid control equations of the fluid computational domain and the solid control equations of the solid computational domain respectively, and the temperature field information is updated interactively at the interface between the fluid computational domain and the solid computational domain. Step S3: During the iterative process of solving the first stage, the preset convergence index is monitored in real time; when the convergence index meets the first preset convergence condition, the first stage solution is stopped, and the convective heat transfer coefficient at the current interface is extracted. The first preset convergence condition indicates that the global temperature field of the fluid-structure interaction calculation model has not reached the final fully steady-state convergence state, and the convective heat transfer coefficient has reached the stable convergence state. Step S4, Second stage solution: The extracted convective heat transfer coefficient is applied as a thermodynamic boundary condition to the target computational domain, and the target computational domain is solved by one-way decoupling until the final fully converged state is reached to obtain the steady-state temperature field distribution. The unidirectional decoupling solution refers to the fact that the bidirectional interaction of the temperature field between the fluid computational domain and the solid computational domain is no longer performed during the solution process.

2. The calculation method for fluid-structure interaction heat transfer according to claim 1, characterized in that, The target computational domain is the solid-state computational domain; In step S4, the convective heat transfer coefficient is used as a third type of convective heat transfer boundary condition and applied only to the solid computational domain. During the unidirectional decoupling solution process, the update of the fluid computational domain is frozen and the solution of the fluid control equations of the fluid computational domain is stopped. The solid control equations of the solid computational domain are solved only to obtain the steady-state temperature field distribution of the solid computational domain.

3. The calculation method for fluid-structure interaction heat transfer according to claim 2, characterized in that, The fluid governing equations are the incompressible Navier-Stokes equations, and the solid governing equations are the solid heat conduction equations.

4. The calculation method for fluid-structure interaction heat transfer according to claim 3, characterized in that, In step S3, the first preset convergence condition includes at least one of the following: Condition A: The normalized gradient of the temperature field at the interface or within the computational domain with respect to the iteration time step falls back to the range of 0.01 to 0.001; Condition B: Assume the current iteration time step is The historical regional average value of the global temperature field ave ( T )satisfy Where T represents temperature and the convective heat transfer coefficient is... h Gradient of residual variation with iteration time step satisfy .

5. The calculation method for fluid-structure interaction heat transfer according to claim 4, characterized in that, The convective heat transfer coefficient h The calculation formula is: In the above formula, Q The total convective heat transfer through the interface, A Let be the area of ​​the interface. T w The wall temperature on the solid side of the interface is [temperature value missing]. T ref The reference temperature is the fluid side of the interface.

6. The calculation method for fluid-structure interaction heat transfer according to claim 5, characterized in that, The specific expression for the solid heat conduction equation is as follows: In the above formula, For solid density, For isobaric specific heat capacity, T For temperature, t For time, Thermal conductivity, For temperature gradient, The divergence of heat flux density, The intensity of the heat source within the volume.

7. The calculation method for fluid-structure interaction heat transfer according to claim 6, characterized in that, The specific expression for the incompressible Navier-Stokes equation is as follows: In the above formula, velocity vector For velocity divergence, the volume expansion rate of fluid particles is zero under incompressible conditions; For fluid density, For local acceleration, For convective acceleration, For pressure gradient, For dynamic viscosity, For the velocity Laplace operator, It is a volume force vector.

8. The calculation method for fluid-structure interaction heat transfer according to claim 7, characterized in that, Both the incompressible Navier-Stokes equation and the solid heat conduction equation are solved spatially using the finite volume method.

9. A computational system for fluid-structure interaction heat transfer, characterized in that, include: The model initialization module is used to obtain a fluid-structure interaction calculation model that includes both fluid and solid computational domains, and to set the initial velocity, pressure, and temperature fields. The first-stage solution module is used to solve the fluid control equations of the fluid computational domain and the solid control equations of the solid computational domain respectively within consecutive iteration time steps using a fluid-solid weak coupling heat transfer algorithm, and interactively update the temperature field information at the interface between the fluid computational domain and the solid computational domain. The monitoring and extraction module is used to monitor the preset convergence index in real time during the iterative process of the first stage solution; when the convergence index meets the first preset convergence condition, the first stage solution is stopped and the convective heat transfer coefficient at the current interface is extracted; wherein, the first preset convergence condition indicates that the global temperature field of the fluid-structure interaction calculation model has not reached the final fully steady-state convergence state, and the convective heat transfer coefficient has reached the stable convergence state. The second-stage solution module is used to apply the extracted convective heat transfer coefficient as a thermodynamic boundary condition to the target computational domain, and to perform a one-way decoupled solution on the target computational domain until the final fully converged state is reached to obtain the steady-state temperature field distribution; wherein, the one-way decoupled solution means that the bidirectional interaction of the temperature field between the fluid computational domain and the solid computational domain is no longer performed during the solution process.

10. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the calculation method for fluid-structure interaction heat transfer as described in any one of claims 1 to 8.

Images