Fluid-structure interaction calculation method and apparatus based on wall temperature heat flux correction

Abstract

This invention relates to a fluid-structure interaction (FSI) calculation method and apparatus based on wall temperature heat flux correction, belonging to the field of FSI calculation technology. The method includes: calculating initial wall heat flux and initial surface pressure; using the initial wall heat flux as the thermal boundary condition, performing transient heat conduction calculations on the structure within an initial time interval to obtain a first structural temperature distribution; calculating structural deformation to obtain the first wall heat flux and the first surface pressure; repeating the steps to obtain a second structural temperature distribution; calculating the corrected heat flux using a heat flux interpolation model within the first time interval; and determining whether the corrected heat flux exceeds an accuracy threshold; if not, skipping the structural deformation and fluid domain calculations at the second time step and proceeding to the second coupling step; if so, calculating the structural deformation and fluid domain at the second time step and proceeding to the second coupling step; until the preset total simulation time is reached. This method can ensure calculation accuracy while improving simulation efficiency.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of fluid-structure coupling calculation, in particular to a fluid-structure thermal coupling calculation method and device based on wall temperature heat flow correction. BACKGROUND

[0002] During the flight of an aircraft, aerodynamic heating and aerodynamic force act on the surface of the aircraft together, resulting in significant temperature rise and thermal deformation of the structure. In order to accurately predict the thermal response and geometric change of the structure, multi-field coupling calculation needs to be carried out between the fluid solver, the structure thermal solver and the structure mechanics solver. Since the characteristic time of the fluid domain (millisecond level) is much smaller than the characteristic time of the structure heat transfer and deformation (second level), the existing method generally adopts a loose coupling strategy, that is, in each coupling step, the heat flow and pressure are provided by the fluid solver, and the temperature and deformation are calculated by the structure solver according to the heat flow and pressure.

[0003] However, in the existing improved scheme, some methods adjust the coupling step size by monitoring the structure temperature rise threshold, which can avoid the invalid calculation of fixed small step size, but the wall heat flow is not dynamically corrected throughout the simulation, and only the initial heat flow calculated by the fluid domain is used as a constant load to the solid domain, ignoring the physical correlation between the heat flow and the wall temperature, resulting in the accumulation of heat flow precision deviation in long-time simulation. At the same time, such methods force to calculate the structure deformation and update the fluid domain mesh at the end of each coupling step, without considering the actual influence of deformation on the coupling accuracy, resulting in unnecessary computational overhead. Another scheme adopts a multi-level adaptive step size design, which first uses a fixed small step size to quickly obtain the macroscopic characteristics of the cumulative thermal deformation, and then carries out high-precision calculation according to the regional temperature variation rate. Although this method realizes the regional adaptation of the step size, the heat flow correction is only used as an auxiliary means for the first level calculation, which cannot essentially reduce the number of fluid domain recalculation, and the structure deformation trigger only depends on the temperature rise threshold, without considering the synergistic effect of mechanical load on the deformation, which deviates from the real physical process. SUMMARY

[0004] Therefore, it is necessary to provide a fluid-structure thermal coupling calculation method and device based on wall temperature heat flow correction, which can ensure the calculation accuracy and improve the calculation efficiency of the simulation.

[0005] A fluid-structure thermal coupling calculation method based on wall temperature heat flow correction, the method comprising:

[0006] Step 1, setting a coupling step size expansion factor and a precision threshold; calculating an initial wall heat flow and an initial surface pressure at an initial time;

[0007] Step 2, using the initial wall heat flow as a thermal boundary condition, performing structure transient heat conduction calculation in an initial time interval to obtain a first structure temperature distribution at a first time;

[0008] Step 3, in the initial coupling step, based on the first structure temperature distribution and the initial surface pressure, a structure deformation is calculated, the structure deformation is updated to the fluid domain grid for calculation, and the first wall heat flow and the first surface pressure at the first time are obtained;

[0009] Step 4, in the first coupling step, based on the first wall heat flow, step 2 is repeated to obtain the second structure temperature distribution at the second time; and according to the initial wall heat flow, the initial structure temperature distribution, the first wall heat flow and the first structure temperature distribution, a correction heat flow is calculated through a heat flow interpolation model;

[0010] Step 5, it is judged whether the correction heat flow exceeds the precision threshold; if not, the structure deformation and the fluid domain calculation at the second time are skipped, the second wall heat flow and the second surface pressure are obtained, and then the second time interval corresponding to the second coupling step is updated according to the coupling step length expansion factor, and the second coupling step is entered; if yes, step 3 is performed to calculate the structure deformation and the fluid domain at the second time, the second wall heat flow and the second surface pressure are obtained, and then the second time interval corresponding to the second coupling step is updated according to the coupling step length expansion factor, and the second coupling step is entered;

[0011] Step 6, steps 4 to 5 are repeated until a preset total simulation time is reached.

[0012] In another aspect, a fluid-structure thermal coupling calculation device based on wall temperature heat flow correction is also provided, comprising:

[0013] An initial calculation module is configured to set a coupling step length expansion factor and a precision threshold; and calculate an initial wall heat flow and an initial surface pressure at an initial time;

[0014] A structure temperature distribution calculation module is configured to perform structure transient heat conduction calculation in an initial time interval with the initial wall heat flow as a thermal boundary condition to obtain a first structure temperature distribution at a first time;

[0015] A first time calculation module is configured to, in the initial coupling step, calculate a structure deformation based on the first structure temperature distribution and the initial surface pressure, update the structure deformation to the fluid domain grid for calculation, and obtain a first wall heat flow and a first surface pressure at the first time;

[0016] A second time calculation module is configured to, in the first coupling step, repeat the process of the structure temperature distribution calculation module based on the first wall heat flow to obtain a second structure temperature distribution at a second time; and calculate a correction heat flow through a heat flow interpolation model according to the initial wall heat flow, the initial structure temperature distribution, the first wall heat flow and the first structure temperature distribution;

[0017] a judgment module, configured to judge whether the corrected heat flow exceeds the precision threshold value; if not, skip the structure deformation and fluid domain calculation at the second time, obtain the second wall heat flow and the second surface pressure, and then update the second time interval corresponding to the second coupling step according to the coupling step expansion factor, and enter the second coupling step; if yes, perform the process of the first time calculation module, calculate the structure deformation and the fluid domain at the second time, obtain the second wall heat flow and the second surface pressure, and then update the second time interval corresponding to the second coupling step according to the coupling step expansion factor, and enter the second coupling step;

[0018] a simulation calculation module, configured to repeat the process from the second time calculation module to the judgment module until the preset total simulation time is reached.

[0019] Compared with the prior art, the fluid-solid thermal coupling calculation method and device based on wall temperature heat flow correction provided by the application have the following beneficial effects:

[0020] 1. Compared with the traditional technology of adopting "fixed heat flow load" or "heat flow correction as an auxiliary means", the application fully utilizes the strong linear relationship between the wall heat flow and the temperature, calculates the corrected heat flow in the time interval corresponding to the coupling step through the heat flow interpolation model, realizes the dynamic adaptation of the heat flow under the non-equidistant step, effectively avoids the accumulation of heat flow precision deviation in long-time simulation, significantly improves the heat flow calculation precision, and makes the heat flow calculation more consistent with the real physical law.

[0021] 2. When the corrected heat flow does not exceed the precision threshold value, the active expansion strategy of updating the next coupling step length by the coupling step length expansion factor is adopted, without maintaining a fixed small step length, but dynamically adjusting the step length according to the heat flow correction precision, which not only reduces the number of invalid coupling iterations, but also avoids the additional overhead caused by frequent fluid domain recalculation, compared with the traditional fixed step length or passive adjustment step length method, the calculation efficiency is greatly optimized, and the engineering application cost is reduced.

[0022] 3. The structure deformation is triggered only at the initial coupling step and when the heat flow correction exceeds the precision threshold value, and the structure deformation calculation combines the structure temperature distribution and the surface pressure, which not only ensures the accuracy of the deformation calculation, but also maximizes the unnecessary calculation consumption, and realizes the balance between precision and efficiency. BRIEF DESCRIPTION OF DRAWINGS

[0023] In order to more clearly illustrate the technical solutions of the embodiments of the application, the following will briefly introduce the drawings needed in the embodiments. It should be understood that the following drawings only show some embodiments of the application, and those skilled in the art can obtain other related drawings according to these drawings without creating labor.

[0024] Figure 1 It is a flowchart of the fluid-solid thermal coupling calculation method based on wall temperature heat flow correction in an embodiment.

[0025] Figure 2 a structural block diagram of a fluid-structure thermal coupling calculation device based on wall temperature heat flow correction in an embodiment;

[0026] Figure 3 an internal structure diagram of a computer device in an embodiment.

[0027] The purposes, functional features and advantages of the present application will be further described with reference to the embodiments and the accompanying drawings. DETAILED DESCRIPTION

[0028] The technical solutions in the embodiments of the present application will be clearly and completely described with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the protection scope of the present application.

[0029] It should be noted that the description of “first”, “second” and the like in the present application is only for the purpose of description, and cannot be understood as indicating or implying the relative importance of the indicated technical features or implicitly indicating the number of the indicated technical features. Therefore, the features limited by “first”, “second” can explicitly or implicitly include at least one of the features. In the description of the present application, the meaning of “plurality” is at least two, for example, two, three, etc., unless otherwise specifically limited.

[0030] It can be understood that the technical solutions of various embodiments of the present application can be combined with each other, but it must be based on the fact that a person skilled in the art can realize it, and when the combination of technical solutions appears contradictory or unachievable, it should be considered that the combination of technical solutions does not exist, nor within the protection scope required by the present application.

[0031] The embodiments of the present application will be described in detail with reference to the drawings in the embodiments of the present application.

[0032] Embodiment 1

[0033] The present application aims to overcome the problems of low coupling efficiency, limited time step and ignoring structural thermal deformation feedback in the prior art, and provides a fluid-structure thermal coupling calculation method based on wall temperature heat flow correction. Under the premise of considering the influence of structural thermal deformation on the geometry of the fluid domain, the method fully utilizes the strong linear relationship between the wall heat flow and the wall temperature, dynamically expands the coupling time step, combines the temperature-driven heat flow interpolation strategy, and cooperates with the grid deformation, to significantly improve the calculation efficiency of long-time simulation while ensuring the calculation accuracy.

[0034] As Figure 1As shown, the fluid-structure-thermal coupling calculation method based on wall temperature heat flux correction provided in this embodiment includes the following steps:

[0035] Step 1: Set the coupling step size expansion factor and accuracy threshold; calculate the initial wall heat flux and initial surface pressure at the initial moment.

[0036] Step 2: Using the initial wall heat flow as the thermal boundary condition, perform transient heat conduction calculations on the structure within the initial time interval to obtain the first structural temperature distribution at the first moment.

[0037] Step 3: In the initial coupling step, the structural deformation is calculated based on the first structural temperature distribution and the initial surface pressure. The structural deformation is then updated to the fluid domain mesh for calculation to obtain the first wall heat flow and the first surface pressure at the first moment.

[0038] Step 4: In the first coupling step, based on the first wall heat flow, repeat step 2 to obtain the second structure temperature distribution at the second time; calculate the corrected heat flow through the heat flow interpolation model based on the initial wall heat flow, the initial structure temperature distribution, the first wall heat flow and the first structure temperature distribution.

[0039] Step 5: Determine whether the corrected heat flux exceeds the accuracy threshold. If not, skip the calculation of structural deformation and fluid domain at the second time step, obtain the second wall heat flux and the second surface pressure, and then update the second time interval corresponding to the second coupling step according to the coupling step size expansion factor, and enter the second coupling step. If yes, proceed to step 3, calculate the structural deformation and fluid domain at the second time step, obtain the second wall heat flux and the second surface pressure, and then update the second time interval corresponding to the second coupling step according to the coupling step size expansion factor, and enter the second coupling step.

[0040] Step 6: Repeat steps 4 to 5 until the preset total simulation time is reached.

[0041] In the specific implementation of step 1, the calculation parameters of the fluid domain and the solid domain are first initialized. These parameters include the computational grid, physical property parameters, initial conditions, and boundary conditions for both domains. The range of values ​​for the coupling step size scaling factor is as follows: The preferred value is 1.5, which balances the stability and acceleration of the coupled computation; precision threshold. To correct the relative deviation threshold between the heat flux and the actual wall heat flux in the previous coupling step, a threshold of 5% is preset and used to determine the effectiveness of the heat flux correction.

[0042] At the initial moment At that time, the high-precision Navier-Stokes fluid solver is invoked, and based on the initial and boundary conditions set in step 1, the initial structural temperature distribution is given. As boundary conditions of the fluid domain, steady-state fluid domain calculations are performed on the fluid domain at the initial moment to obtain the initial wall heat flux. Initial surface pressure .

[0043] This step lays the foundation for the orderly advancement of subsequent coupled calculations by clarifying the initial calculation parameters and core control thresholds, avoiding uncontrolled calculation accuracy or low efficiency due to parameter ambiguity. Simultaneously, it acquires the initial load data for the coupled calculations. The initial wall heat flux provides the basic thermal boundary conditions for subsequent structural heat conduction calculations; the initial surface pressure provides the initial mechanical loads for subsequent structural deformation calculations, ensuring that the coupled calculations start from a true initial state; and the initial structural temperature distribution is a natural product of fluid domain calculations, supporting the integrity of the overall data chain.

[0044] In the specific implementation of step 2, the initial coupling time step is first set. The initial time interval corresponding to the initial coupling step is determined based on the initial coupling time step size. Among them, the first moment It is worth noting that the first moment is both the end and the beginning of the initial coupling step. The value of the initial coupling time step must satisfy the following: The change in wall heat flux within the interval is less than 5% of the initial wall heat flux, ensuring that no heat flux correction is needed to guarantee calculation accuracy within this interval.

[0045] Then, several intermediate times are set within the initial time interval. ( ), The value is a positive integer and is relatively small. The transient heat conduction structure solver is invoked, with the initial wall heat flux as the reference value. For the thermal boundary conditions of the solid domain, from the initial time Beginning, based on several intermediate moments By progressively advancing the transient heat conduction calculations of the structure, the structural temperature distribution at each intermediate moment is obtained. Finally, based on the structural temperature distribution at several intermediate moments... The temperature distribution of the first structure at the first moment is obtained by calculating the transient heat conduction integral. .

[0046] This step, as the first structural heat conduction calculation, does not perform heat flow correction due to the lack of historical heat flow data, which is consistent with the actual data conditions. At the same time, it achieves accurate integration by subdividing intermediate time points, providing reliable thermal load data for subsequent structural deformation calculations.

[0047] In the specific implementation of step 3, that is, in the initial coupling step, the structural mechanics solver is first invoked, based on the temperature distribution of the first structure. With initial surface pressure Solve the thermoelastic equation to calculate the structural deformation displacement field. The calculation process also considers thermal load (the structural temperature distribution at the first moment). ) and mechanical load (initial surface pressure) This aligns with the true physical causes of structural deformation during flight; subsequently, the structural deformation displacement field was analyzed. The data is transferred to the fluid domain mesh, and a hybrid method combining the spring method and the background mesh is used to deform and update the fluid domain mesh, ensuring that the fluid domain geometry is consistent with the actual geometry of the structure after deformation. Finally, the Navier-Stokes fluid solver is rerun on the updated fluid domain mesh, and the first structural temperature distribution calculated in the initial coupled step is applied. Set the boundary conditions of the fluid domain to obtain the first wall heat flux at the first moment. and first surface pressure .

[0048] This step achieved the first closed-loop feedback from structural deformation to fluid domain geometry update to fluid domain data recalculation. The updated data from the first coupling step provided crucial historical data for subsequent thermal flow correction, and at the same time verified the interaction between the structure and the fluid domain under the initial load.

[0049] In the specific implementation of steps 4 and 5, that is, in the first coupling step, the coupling step size expansion factor is first determined. Calculate the coupling time step of the first coupling step Then, based on the coupling time step of the first coupling step... The second moment at the end of the first coupling step is determined as: Then, based on the first and second moments, the time interval corresponding to the first coupling step is determined as follows: It is worth noting that in the second moment... It is both the end time of the first coupling step and the start time of the second coupling step.

[0050] Then, the heat flow from the first wall surface As the thermal boundary condition for the solid domain, step 2 is repeated to perform transient heat conduction calculations on the structure within the first time interval, yielding the result at the second time interval. Second structure temperature distribution Then, the corrected heat flux is calculated, and it is determined whether the corrected heat flux exceeds the accuracy threshold, so as to determine whether structural deformation and fluid domain calculations need to be performed at the second time step.

[0051] Specifically, based on the initial wall heat flow Initial structural temperature distribution Heat flow at the first wall Temperature distribution of the first structure The corrected heat flux is calculated using a heat flux interpolation model.

[0052] The expression for the heat flux interpolation model is:

[0053] (1)

[0054] In the formula, Indicates the first In the coupling step The wall heat flow at an intermediate moment; Indicates time The heat flow on the wall; Indicates time The heat flow on the wall; Indicates time The structural temperature distribution; Indicates time The structural temperature distribution; Indicates the first In the coupling step The structural temperature distribution at each intermediate moment.

[0055] Then, the following formula is used to determine whether the corrected heat flux exceeds the accuracy threshold:

[0056] (2)

[0057] In the formula, Indicates the first In the coupling step The wall heat flow at an intermediate moment; Indicates time The heat flow on the wall; Indicates the precision threshold.

[0058] If formula (2) holds true, it means that the corrected heat flux has not exceeded the accuracy threshold, indicating that the corrected heat flux meets the accuracy requirements. As the current solid-domain thermal boundary condition, we continue to advance the transient heat conduction calculation of the structure, and directly skip the structural deformation and fluid domain calculations at the second time interval until the end of the first time interval, thus obtaining the second time interval. Second wall heat flow and second surface pressure Then, the second time interval corresponding to the second coupling step is updated according to the coupling step size expansion factor, and the second coupling step is entered. It should be noted that the magnitude of the second surface pressure is still equal to... ,Right now .

[0059] If formula (2) does not hold, it means that the corrected heat flux exceeds the accuracy threshold, indicating that the accuracy of the corrected heat flux is insufficient. Step 4 needs to be repeated to calculate the structural deformation and fluid domain at the second moment and obtain the heat flux at the second wall. With the second surface pressure Then, the second time interval corresponding to the second coupling step is updated according to the coupling step size expansion factor, and the second coupling step is entered. It should be noted that the magnitude of the second surface pressure at this time is not equal to... ,Right now .

[0060] In the second coupling step, according to the coupling step size expansion factor Update the second time interval corresponding to the second coupling step, specifically by: adjusting the coupling step size expansion factor. Calculate the coupling time step of the second coupling step. Then, based on the coupling time step size of the second coupling step, the third moment at which the second coupling step ends is determined as follows: The time interval corresponding to the second coupling step is determined based on the second and third time points. .

[0061] The coupling time step is calculated based on the coupling step size expansion factor, and the expression is as follows:

[0062] (3)

[0063] In the formula, Indicates the first The coupling time step of the coupling step; Indicates the coupling step size expansion factor; Indicates the first The coupling time step of the coupling step.

[0064] As can be understood, the core of step 4 lies in utilizing the strong linear relationship between wall heat flux and temperature to construct an interpolation model based on historical data to dynamically correct the heat flux. This replaces the frequent fluid domain recalculations in traditional methods, significantly reducing computational overhead. Simultaneously, it supports accurate integration under non-isochronous steps, balancing computational accuracy and flexibility. Step 5 then achieves dynamic adjustment of the coupling step size and on-demand triggering of fluid recalculation through accuracy criteria. This ensures timely correction when computational accuracy is insufficient and maximizes computational efficiency by increasing the step size, avoiding redundant calculations.

[0065] In the specific implementation of step 6, the heat flow correction and structural heat conduction of step 4 and the accuracy judgment, step size adjustment and fluid recalculation of step 5 are executed cyclically. Each cycle is based on the historical data of the first two coupled steps to perform heat flow correction. The coupled step size is dynamically expanded according to the accuracy. Structural deformation and fluid domain recalculation are triggered only when the corrected heat flow exceeds the threshold, forming a complete closed loop of data accumulation, heat flow correction, accuracy judgment and step size adjustment / recalculation.

[0066] Specifically, let the current coupling step be... ,in, In the first In the coupling step, the corresponding time interval is , will the Wall heat flow at the end of the coupling step As the thermal boundary condition for the solid domain, repeat step 2 to obtain the... Structural temperature distribution at the end of the coupling step Then, the corrected heat flux is calculated, and it is determined whether the corrected heat flux exceeds the accuracy threshold to determine whether structural deformation and fluid domain calculations are needed at the second time step.

[0067] Specifically, it needs to be based on the first Structural temperature distribution at the end of the coupling step Wall heat flow Passing the exam Structural temperature distribution at the end of the coupling step Wall heat flow The corrected heat flux is calculated using the heat flux interpolation model of formula (1). During the calculation, the time interval... Several intermediate times are set inside. Then, the corrected heat flow is determined by formula (2) to see if it exceeds the accuracy threshold.

[0068] If formula (2) holds true, it means that the corrected heat flux has not exceeded the accuracy threshold, and the first step is skipped directly. Coupling step end time The structural deformation and fluid domain calculations will As the current thermal boundary condition for the solid domain, we will continue to advance the calculation of transient heat conduction of structures until the [missing information]. The coupling step time interval ends, and the first step is obtained. Wall heat flow at the end of the coupling step and surface pressure Then update the first step according to the coupling step size expansion factor. The time interval corresponding to the coupling step, and then enter the first... Coupling step.

[0069] If formula (2) does not hold, it means that the corrected heat flux exceeds the accuracy threshold, and step 3 needs to be repeated to calculate the second step. Coupling step end time The structural deformation and fluid domain are used to obtain the first... Wall heat flow at the end of the coupling step and surface pressure Then update the first step according to the coupling step size expansion factor. The time interval corresponding to the coupling step, and then enter the first... Coupling step.

[0070] This step ensures the continuous advancement of long-term fluid-structure-thermal coupling simulation of aircraft. Through the cyclical application of core innovative strategies, high fidelity and high computational efficiency are maintained throughout the entire simulation cycle, effectively supporting the need for thermal response prediction of aircraft with long flight times of hours.

[0071] It should be understood that, although this embodiment Figure 1 The steps are shown sequentially as indicated by the arrows, but they are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order in which these steps are performed; they can be executed in other orders. Figure 1 At least some of the steps in the process may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least some of the sub-steps or stages of other steps.

[0072] Example 2

[0073] Based on the fluid-structure-thermal coupling calculation method based on wall temperature heat flow correction in Example 1, this example discloses a fluid-structure-thermal coupling calculation device based on wall temperature heat flow correction, such as... Figure 2 As shown, the fluid-structure-thermal coupling calculation device based on wall temperature heat flux correction includes: an initial calculation module 401, a structural temperature distribution calculation module 402, a first-moment calculation module 403, a second-moment calculation module 404, a judgment module 405, and a simulation calculation module 406, wherein:

[0074] The initial calculation module 401 is used to set the coupling step size expansion factor and accuracy threshold; and to calculate the initial wall heat flux and initial surface pressure at the initial moment.

[0075] The structural temperature distribution calculation module 402 is used to perform transient heat conduction calculation of the structure within the initial time interval using the initial wall heat flow as the thermal boundary condition, so as to obtain the first structural temperature distribution at the first moment.

[0076] The first-moment calculation module 403 is used to calculate the structural deformation based on the first structural temperature distribution and the initial surface pressure in the initial coupling step, update the structural deformation to the fluid domain mesh for calculation, and obtain the first wall heat flux and the first surface pressure at the first moment.

[0077] The second-moment calculation module 404 is used in the first coupling step to repeat the process of the structure temperature distribution calculation module based on the first wall heat flow to obtain the second structure temperature distribution at the second moment; and to calculate the corrected heat flow through the heat flow interpolation model based on the initial wall heat flow, the initial structure temperature distribution, the first wall heat flow and the first structure temperature distribution.

[0078] The judgment module 405 is used to determine whether the corrected heat flux exceeds the accuracy threshold. If not, the calculation of structural deformation and fluid domain at the second time moment is skipped, the second wall heat flux and the second surface pressure are obtained, and then the second time interval corresponding to the second coupling step is updated according to the coupling step size expansion factor, and the second coupling step is entered. If yes, the process of the first time moment calculation module is performed to calculate the structural deformation and fluid domain at the second time moment, obtain the second wall heat flux and the second surface pressure, and then the second time interval corresponding to the second coupling step is updated according to the coupling step size expansion factor, and the second coupling step is entered.

[0079] The simulation calculation module 406 is used to repeat the process from the calculation module to the judgment module at the second time step until the preset total simulation time is reached.

[0080] In this embodiment, the specific working process and working principle of the initial calculation module 401, the structural temperature distribution calculation module 402, the first moment calculation module 403, the second moment calculation module 404, the judgment module 405, and the simulation calculation module 406 are the same as those in Embodiment 1, and therefore will not be described again in this embodiment. Each unit module can be implemented entirely or partially through software, hardware, or a combination thereof. Each unit module can be embedded in the processor of the computer device in hardware form or independent of it, or it can be stored in the memory of the computer device in software form, so that the processor can call and execute the operations corresponding to the above unit modules.

[0081] Example 3

[0082] like Figure 3 The diagram illustrates a terminal device disclosed in this embodiment, comprising a transmitter, a receiver, a memory, and a processor. The transmitter transmits instructions and data, the receiver receives instructions and data, the memory stores computer-executed instructions, and the processor executes the computer-executed instructions stored in the memory to implement the method described in Embodiment 1 above.

[0083] It is important to note that the aforementioned memory can be either standalone or integrated with the processor. When the memory is set up independently, the terminal device also includes a bus for connecting the memory and the processor.

[0084] Example 4

[0085] This embodiment discloses a computer-readable storage medium storing computer-executable instructions. When a processor executes the computer-executable instructions, it implements the method in Embodiment 1 above.

[0086] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0087] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0088] The embodiments described above are merely examples of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention.

Claims

1. A fluid-structure thermal coupling calculation method based on wall surface temperature heat flow correction, characterized in that, The method comprises: Step 1, setting a coupling step length expansion factor and an accuracy threshold; calculating an initial wall heat flow and an initial surface pressure at an initial time; Step 2, performing structure transient heat conduction calculation in an initial time interval with the initial wall heat flow as a thermal boundary condition to obtain a first structure temperature distribution at a first time; Step 3, in the initial coupling step, calculating structure deformation based on the first structure temperature distribution and the initial surface pressure, updating the structure deformation to the fluid domain grid for calculation to obtain a first wall heat flow and a first surface pressure at the first time; Step 4, in the first coupling step, based on the first wall heat flow, repeating step 2 to obtain a second structure temperature distribution at a second time; and calculating a correction heat flow through a heat flow interpolation model according to the initial wall heat flow, the initial structure temperature distribution, the first wall heat flow and the first structure temperature distribution; Step 5, judging whether the correction heat flow exceeds the accuracy threshold; if not, skipping the structure deformation and fluid domain calculation at the second time to obtain a second wall heat flow and a second surface pressure, and then updating a second time interval corresponding to the second coupling step according to the coupling step length expansion factor to enter the second coupling step; if yes, performing step 3 to calculate the structure deformation and the fluid domain at the second time to obtain the second wall heat flow and the second surface pressure, and then updating the second time interval corresponding to the second coupling step according to the coupling step length expansion factor to enter the second coupling step; Step 6, repeating steps 4 to 5 until a preset total simulation time is reached; The expression of the heat flow interpolation model is: ； wherein represents the wall heat flux at the first intermediate time in the coupling step second intermediate time in the coupling step represents the wall heat flux at time represents the wall heat flux at time represents the wall heat flux at time represents the wall heat flux at time represents the structure temperature distribution at time represents the structure temperature distribution at time represents the structure temperature distribution at time represents the structure temperature distribution at time represents the structure temperature distribution at the first intermediate time in the coupling step second intermediate time in the coupling step 2. The method for fluid-structure thermal coupling calculation based on wall surface temperature heat flow correction according to claim 1, characterized in that, In step 2, the structure transient heat conduction calculation in the initial time interval with the initial wall heat flow as the thermal boundary condition to obtain the first structure temperature distribution at the first time comprises: Setting an initial coupling time step, and determining an initial time interval according to the initial coupling time step; Setting a plurality of intermediate times in the initial time interval, and starting from the initial time, performing structure transient heat conduction calculation based on the plurality of intermediate times to obtain structure temperature distributions at the plurality of intermediate times with the initial wall heat flow as the thermal boundary condition; Obtaining the first structure temperature distribution at the first time according to the structure temperature distributions at the plurality of intermediate times.

3. The method for fluid-structure thermal coupling calculation based on wall surface temperature heat flow correction according to claim 1, characterized in that, The coupling step enlargement factor has a value range of: .

4. The fluid-structure thermal coupling calculation method based on wall temperature heat flow correction according to any one of claims 1 to 3, characterized in that, In step 5, judging whether the correction heat flow exceeds the accuracy threshold, the expression is: ； wherein represents the first coupling step the wall heat flux at the intermediate time instant; represents the wall heat flux at the time instant represents the wall heat flux at the time instant represents the precision threshold.

5. The method for fluid-structure thermal coupling calculation based on wall surface temperature heat flow correction according to claim 4, characterized in that, In step 5, updating the second time interval corresponding to the second coupling step according to the coupling step length expansion factor comprises: calculating a coupling time step for the second coupling step according to the coupling step size expansion factor ; and then determining a third time instant at which the second coupling step ends as: ; According to the second time and the third time, the time interval corresponding to the second coupling step is determined as .

6. The method for fluid-structure thermal coupling calculation based on wall surface temperature heat flow correction according to claim 5, characterized in that, Calculating the coupling time step of the coupling step according to the coupling step length expansion factor, the expression is: ； wherein denotes the coupling time step of the coupling step denotes the coupling time step of the coupling step denotes the coupling step length expansion factor denotes the coupling time step of the coupling step denotes the coupling time step of the coupling step 7. The fluid-structure thermal coupling calculation method based on wall temperature heat flow correction according to any one of claims 1 to 3, characterized in that, The structure deformation is triggered only in step 3 and when the correction heat flow exceeds the accuracy threshold, and both the thermal load generated by the structure temperature distribution and the mechanical load generated by the surface pressure are considered.

8. The fluid-structure thermal coupling calculation method based on wall surface temperature heat flow correction according to any one of claims 1 to 3, characterized in that, The grid updating adopts a spring method and a background grid hybrid method.

9. A fluid-structure thermal coupling calculation device based on wall temperature heat flow correction, characterized in that, The device comprises: An initial calculation module for setting a coupling step length expansion factor and an accuracy threshold; calculating an initial wall heat flow and an initial surface pressure at an initial time; a structure temperature distribution calculation module, configured to perform a structure transient heat conduction calculation in an initial time interval with the initial wall heat flux as a thermal boundary condition to obtain a first structure temperature distribution at a first time; a first time calculation module, configured to, in an initial coupling step, calculate a structure deformation based on the first structure temperature distribution and the initial surface pressure, update the structure deformation to a fluid domain grid for calculation to obtain a first wall heat flux and a first surface pressure at the first time; a second time calculation module, configured to, in a first coupling step, repeat the process of the structure temperature distribution calculation module based on the first wall heat flux to obtain a second structure temperature distribution at a second time, and calculate a correction heat flux by a heat flux interpolation model according to the initial wall heat flux, the initial structure temperature distribution, the first wall heat flux and the first structure temperature distribution; a judgment module, configured to judge whether the correction heat flux exceeds the precision threshold; if not, skip the structure deformation and the fluid domain calculation at the second time to obtain a second wall heat flux and a second surface pressure, then update a second time interval corresponding to a second coupling step according to a coupling step length expansion factor and enter the second coupling step; if yes, perform the process of the first time calculation module to calculate the structure deformation and the fluid domain at the second time to obtain the second wall heat flux and the second surface pressure, then update the second time interval corresponding to the second coupling step according to the coupling step length expansion factor and enter the second coupling step; a simulation calculation module, configured to repeat the process from the second time calculation module to the judgment module until a preset total simulation time is reached.

Images