Mesh deformation and particle deposition layer and wall surface transient thermal coupling temperature prediction method

Abstract

The invention provides a grid deformation and particle deposition layer and wall surface transient thermal coupling temperature prediction method. In the modern industrial field, high-efficiency operation and precise maintenance of high-temperature heat flow equipment are crucial. The traditional simulation method neglects the existence of a sedimentary layer or simplifies the sedimentary layer, so that the prediction of the heat transfer characteristics of the thermal boundary has deviation. Through the user-defined function UDF, the deformation of the deposition layer is dynamically predicted, the heat conductivity coefficient of the deposition layer is synchronously updated and assigned, the deposition layer and the solid wall surface are coupled, and then the influence of the particulate matter deposition layer on the heat exchange characteristic of the high-temperature heat flow equipment is predicted. The prediction accuracy of thermal boundary change caused by particle deposition in high-temperature heat flow equipment is improved, and the problem of prediction deviation in traditional simulation is solved. According to the method, the influence of the particulate matter deposition layer on the heat exchange characteristics of the high-temperature heat flow equipment can be accurately predicted, and the method has certain guiding significance on design, operation and maintenance of the high-temperature heat flow equipment through which the particulate matter flows.

Description

Technical Field

[0001] The present invention relates to the field of particle deposition thermal coupling technology, and in particular to a method for predicting grid deformation and transient thermal coupling temperature of a particle deposition layer and a wall surface. Background Art

[0002] Many high-temperature heat flow devices are traversed by soot generated by combustion or fine particles drawn in from the outside. These particles are transported along with the hot flue gas, collide with walls, bounce off them, and deposit. Particle deposition narrows the flow channel, altering heat transfer, particularly near walls, where the deposited particle layer alters the heat transfer characteristics of the wall. Alternatively, the particle layer can accumulate on component surfaces, causing localized high temperatures and, in severe cases, equipment damage. Due to the high cost and time required to conduct particle deposition tests under actual engine conditions, several algorithms have been developed to predict particle deposition.

[0003] To accurately simulate heat conduction between the deposit and the solid wall, it is necessary to identify the mesh representing the deposit structure within the deformed solid domain mesh and distinguish the deposit mesh from the solid structure mesh. During the unsteady numerical calculation process, the DPM injects a fixed number of particles at each time point and tracks the flow of the particles in the fluid. The critical viscosity criterion programmed in the UDF is used to determine particle deposition. After particle tracking is completed, the particle deposition of each unit grid is converted into grid displacement, causing the blade wall to deform. The coupled grid deformation is then used to drive the synchronous displacement of the solid wall and its boundary layer. Therefore, it is necessary to identify the deformed solid domain mesh in order to distinguish the deposit layer mesh from the blade metal mesh. However, there is no specific function in the UDF that can represent the internal mesh that has undergone displacement. Therefore, it is necessary to propose a deposit mesh identification method based on mesh movement to dynamically track the growth process of the deposit layer. Summary of the Invention

[0004] In response to the aforementioned problem of thermal boundary changes caused by particle transport and particle deposition in high-temperature heat flow equipment, which affects its design, operation, and maintenance, a method for predicting grid deformation and transient thermal coupling temperature between particle deposition layers and wall surfaces is provided. The present invention primarily employs a method for identifying a deposited grid based on grid movement. By dynamically tracking the growth of the deposition layer, a numerical method can accurately identify the deposition area grid in a relatively complex geometric model. Based on extraordinary flow and particle deposition, particle deposition and deposition layer growth are predicted. The effects of the presence of the deposition layer on the flow field and wall convection heat transfer, as well as the heat conduction of the deposition layer on the solid wall, are transiently calculated, enabling dynamic simulation of the deposition effect of particles in high-temperature heat flow equipment.

[0005] The technical means adopted in the present invention are as follows:

[0006] A method for predicting the temperature of a grid deformation and transient thermal coupling between a particle deposition layer and a wall surface, comprising:

[0007] Establish the fluid domain and solid domain of the geometric model, establish the connection between the coincident surfaces of the fluid domain and the solid domain, and mesh the fluid domain and the solid domain respectively;

[0008] Set different time steps for the fluid domain and solid domain to simulate and solve the flow equations, energy equations and turbulence equations in the unsteady flow field;

[0009] Initialize using the user-defined function UDF and record the initial boundary surface mesh center coordinate position P (X, Y, Z) and normal n (x, y, z) direction;

[0010] Release particles in the fluid domain and determine their motion trajectories by tracking the particle trajectories using the Lagrangian method;

[0011] Set UDF custom boundary conditions and determine the particle deposition state through the particle deposition model;

[0012] When particles are deposited on the solid surface, the physical properties of the particles are recorded in the corresponding grid cells and the height of the deposited layer is calculated;

[0013] The sediment is evenly distributed to each node of the grid unit. Based on the dynamic mesh technology, the grid nodes at the interface between the fluid domain and the solid domain are moved simultaneously to achieve simultaneous boundary deformation of the fluid domain and the solid domain.

[0014] Mark the sediment layer mesh after boundary deformation in the solid domain, solve the fluid flow and heat transfer in the fluid domain and solid domain, and obtain the temperature of the sediment and solid domain;

[0015] The UDF is used to traverse the grid, extract and calculate the temperature of the grid where the sediment layer contacts the original solid, and solve the original boundary surface temperature under the sediment layer.

[0016] Furthermore, the overlapping surfaces between the fluid domain and the solid domain are connected through shared topology, the grid division adopts unstructured grid, the grid nodes at the junction of the fluid domain and the solid domain are completely corresponding, and boundary layer encryption is set on both sides of the junction.

[0017] Furthermore, the releasing of particles in the fluid domain and determining particle motion trajectories by tracking particle trajectories using the Lagrangian method specifically include:

[0018] Particles are released into the fluid domain through the discrete phase model (DPM) and their migration trajectories are calculated. The Lagrangian discrete phase model is used to track the motion trajectories of the particles, and the Navier-Stokes equations are solved to treat the gas as a continuous phase. When calculating the migration trajectories, the drag, gravity, Saffman force, pressure gradient force, thermophoretic force, added mass force and Brownian force of the particles are taken into account, and the mutual influence between turbulence and particles is considered through a two-way coupling method.

[0019] Furthermore, determining the deposition state of the particles using the particle deposition model specifically includes:

[0020] The particle deposition model calculates the softening temperature T according to the particle composition. s , and the critical viscosity μ corresponding to the softening temperature critical ; When the calculated viscosity reaches or falls below the critical viscosity, or the particle T p When the temperature reaches or exceeds the softening temperature, the adhesion probability of the particles is assumed to be 1, otherwise, the adhesion probability P s (T p ) is calculated as follows:

[0021]

[0022] Among them, μ Tp is the viscosity of the particle, μ soft is the particle viscosity corresponding to the particle softening temperature.

[0023] Furthermore, the calculation method of the sedimentation layer height is as follows:

[0024]

[0025] Where n(x) represents the particle size in area A. cell The amount of sediment in the grid cell, m particle represents the mass of a single particle and ρ particle represents the density of a single particle; C represents the magnification factor.

[0026] Furthermore, the dynamic mesh technology realizes the movement of mesh nodes on the interface between the fluid domain and the solid domain after particle deposition through mesh smoothing and reconstruction. The distance and direction of the mesh node movement are controlled by a user-defined function (UDF).

[0027] Furthermore, the original boundary surface temperature T face The calculation method is as follows:

[0028]

[0029] in, is the thermal conductivity of the sedimentary layer grid; is the distance from the center point of the sedimentary layer grid to the surface; is the temperature of the center point of the sedimentary layer grid; is the thermal conductivity of the raw material grid; is the distance from the center point of the raw material grid to the surface; T2 is the temperature of the center point of the raw material grid.

[0030] Compared with the prior art, the present invention has the following advantages:

[0031] The method for predicting the temperature of the grid deformation and transient thermal coupling between the particle deposition layer and the wall provided by the present invention not only simulates the deposition law of particulate matter, but also takes into account the influence of the sediment on the flow field, and further introduces the coupling effect between the particulate matter and the unsteady flow field and the thermal resistance effect of the sediment. By introducing an amplification factor in the particle concentration and the sediment height, the dynamic deposition process of particulate matter in the high-temperature heat flow equipment can be simulated more realistically. In addition, the present invention innovatively proposes a multi-grid identification technology based on boundary layer movement, which can accurately identify the deformed solid domain grid and accurately assign the thermal conductivity of the sediment layer, thereby more accurately calculating the particle deposition and the influence of the sediment on the thermal boundary in the high-temperature heat flow equipment.

[0032] Based on the above reasons, the present invention can be widely promoted in the field of particle deposition thermal coupling technology. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative labor.

[0034] Figure 1 This is a flow chart of the method for predicting the temperature of the grid deformation and the transient thermal coupling between the particle deposition layer and the wall surface of the present invention.

[0035] Figure 2 Schematic diagram of the particle-plate impact model and grid in an embodiment of the present invention.

[0036] Figure 3 This is a cloud diagram comparing the experimental and calculated results of the particle deposition characteristics on the target surface and the temperature on the back of the target after the flat plate impact in an embodiment of the present invention.

[0037] Figure 4 This is a comparison chart of experimental and numerical results of the average temperature within circles of different diameters on the back of the target just below the center line in an embodiment of the present invention. DETAILED DESCRIPTION

[0038] It should be noted that, in the absence of conflict, the embodiments and features of the embodiments of the present invention can be combined with each other. The present invention will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0039] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. The following description of at least one exemplary embodiment is actually only illustrative and is in no way intended to limit the present invention and its application or use. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0040] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or combinations thereof.

[0041] Unless otherwise specified, the relative arrangement of the parts and steps, numerical expressions and numerical values ​​set forth in these embodiments do not limit the scope of the present invention. At the same time, it should be clear that, for ease of description, the sizes of the various parts shown in the accompanying drawings are not drawn according to actual proportional relationships. The technology, methods and equipment known to those of ordinary skill in the relevant art may not be discussed in detail, but in appropriate cases, the technology, methods and equipment should be considered as a part of the specification. In all examples shown and discussed here, any specific value should be interpreted as being merely exemplary, rather than as a limitation. Therefore, other examples of the exemplary embodiments can have different values. It should be noted that similar numbers and letters represent similar items in the following drawings, and therefore, once an item is defined in an accompanying drawing, it does not need to be further discussed in subsequent drawings.

[0042] In the description of the present invention, it should be understood that the directions or positional relationships indicated by directional words such as "front, back, up, down, left, right", "horizontal, vertical, vertical, horizontal" and "top, bottom" are usually based on the directions or positional relationships shown in the accompanying drawings. They are only for the convenience of describing the present invention and simplifying the description. Unless otherwise specified, these directional words do not indicate or imply that the device or element referred to must have a specific direction or be constructed and operated in a specific direction. Therefore, they cannot be understood as limiting the scope of protection of the present invention: the directional words "inside and outside" refer to the inside and outside relative to the outline of each component itself.

[0043] For ease of description, spatially relative terms such as "above", "above", "on the upper surface of", "above", etc. may be used herein to describe the spatial positional relationship of a device or feature to other devices or features as shown in the figures. It should be understood that spatially relative terms are intended to include different orientations of the device in use or operation in addition to the orientation described in the figures. For example, if the device in the drawings is inverted, the device described as "above other devices or structures" or "above other devices or structures" will be positioned as "below other devices or structures" or "below their position devices or structures". Thus, the exemplary term "above" can include both "above" and "below". The device can also be positioned in other different ways (rotated 90 degrees or in other orientations), and the spatially relative descriptions used here are interpreted accordingly.

[0044] In addition, it should be noted that the use of terms such as "first" and "second" to limit components is only for the convenience of distinguishing the corresponding components. Unless otherwise stated, the above terms have no special meaning and therefore cannot be understood as limiting the scope of protection of the present invention.

[0045] like Figure 1 As shown, the present invention provides a method for predicting the temperature of a grid deformation and a particle deposition layer and a wall transient thermal coupling, comprising:

[0046] Establish the fluid domain and solid domain of the geometric model, establish the connection between the coincident surfaces of the fluid domain and the solid domain, and mesh the fluid domain and the solid domain respectively;

[0047] In specific implementation, as a preferred embodiment of the present invention, the overlapping surfaces between the fluid domain and the solid domain are connected through shared topology, the grid division adopts unstructured grid, the grid nodes at the junction of the fluid domain and the solid domain are completely corresponding, and boundary layer encryption is set on both sides of the junction.

[0048] Different time steps are set for the fluid and solid domains to simulate and solve the flow, energy, and turbulence equations in the unsteady flow field. During implementation, an appropriate turbulence model is selected based on the operating conditions of the high-temperature equipment, along with the DPM, material properties, and boundary conditions. Different time steps are set for the fluid and solid domains to perform unsteady numerical simulations. After the solution converges, a certain number of UDMs (user-defined memory) are allocated through a UDF to record the initial boundary surface mesh center coordinates and normal directions, with the recording time being t = 0.

[0049] Initialize using the user-defined function UDF and record the initial boundary surface mesh center coordinate position P (X, Y, Z) and normal n (x, y, z) direction;

[0050] Release particles in the fluid domain and determine their motion trajectories by tracking the particle trajectories using the Lagrangian method;

[0051] In a specific implementation, as a preferred embodiment of the present invention, the releasing of particles in the fluid domain and determining the particle motion trajectory by tracking the particle trajectory using the Lagrangian method specifically include:

[0052] Particles are released into the fluid domain through the discrete phase model (DPM) and their migration trajectories are calculated. The Lagrangian discrete phase model is used to track the motion trajectories of the particles, and the Navier-Stokes equations are solved to treat the gas as a continuous phase. When calculating the migration trajectories, the drag, gravity, Saffman force, pressure gradient force, thermophoretic force, added mass force and Brownian force of the particles are taken into account, and the mutual influence between turbulence and particles is considered through a two-way coupling method.

[0053] During implementation, an appropriate turbulence model is selected according to the calculation requirements, and the initial temperature of the particles is initialized through a user-defined function (UDF).

[0054] Set UDF custom boundary conditions and determine the particle deposition state through the particle deposition model;

[0055] In specific implementation, as a preferred embodiment of the present invention, the determination of the deposition state of particles by using a particle deposition model specifically includes:

[0056] The particle deposition model calculates the softening temperature T according to the particle composition. s , and the critical viscosity μ corresponding to the softening temperature critical ; When the calculated viscosity reaches or falls below the critical viscosity, or the particle T p When the temperature reaches or exceeds the softening temperature, the adhesion probability of the particles is assumed to be 1, otherwise, the adhesion probability P s (T p ) is calculated as follows:

[0057]

[0058] Among them, μ Tp is the viscosity of the particle, μ soft is the particle viscosity corresponding to the particle softening temperature.

[0059] The deposition model is embedded into the ANSYS Fluent solver via UDF.

[0060] When particles are deposited on a solid surface, the physical properties of the particles are recorded in the corresponding grid cells, and the height of the deposited layer is calculated. The storage of data on grid cells, surfaces, and nodes, as well as the transfer between them, are achieved through user-defined memory macros.

[0061] In specific implementation, as a preferred embodiment of the present invention, the method for calculating the height of the sedimentary layer is as follows:

[0062]

[0063] Where n(x) represents the particle size in area A. cell The amount of sediment in the grid cell, m particle represents the mass of a single particle and ρ particle represents the density of a single particle; C represents the amplification factor. In order to simulate the actual particle deposition phenomenon on a large time scale, the amplification factor C needs to be multiplied when calculating the grid deformation height.

[0064] The sediment is evenly distributed to each node of the grid unit. Based on the dynamic mesh technology, the grid nodes at the interface between the fluid domain and the solid domain are moved simultaneously to achieve simultaneous boundary deformation of the fluid domain and the solid domain.

[0065] In specific implementation, as a preferred embodiment of the present invention, the dynamic mesh technology realizes the movement of mesh nodes on the interface between the fluid domain and the solid domain after particle deposition through mesh smoothing and reconstruction, and the distance and direction of the mesh node movement are controlled by a user-defined function (UDF).

[0066] The sedimentary layer mesh after boundary deformation is marked in the solid domain, and the fluid flow and heat transfer in the fluid and solid domains are solved to obtain the temperature of the sediment and solid domains. At this time, the original solid surface has been covered by the sediment layer, and its temperature cannot be extracted. Therefore, a UDF is used to traverse the mesh and extract and calculate the temperature of the mesh where the sediment layer contacts the original solid to solve the original boundary surface temperature under the sediment layer.

[0067] After the solid wall mesh is deformed, sediment is identified by using a user-defined function (UDF) to determine whether the mesh cells containing particle deposits have moved. If so, the solid boundary layer mesh that follows the movement is traversed to determine whether it is a sedimentary layer mesh. The thermal conductivity of the sediment is then assigned using a user-defined function (UDF) of physical property parameters.

[0068] Use UDF to traverse the boundary surface mesh and extract the center point position P of the boundary surface mesh * , and compare it with the coordinate position P of the recorded initial wall grid center point. If the positions are inconsistent, it is determined that the grid node has displaced. Along the normal direction n of the deformed boundary surface grid, find the internal grid of the solid domain connected to the deformed grid, mark the internal grid as cell1, and extract the coordinates of the grid center point Use the vector direction determination method to determine whether the grid belongs to the sedimentary layer grid. The determination method is as follows:

[0069] n2=(X1-X,Y1-Y,Z1-Z)

[0070]

[0071] If n*n2>0, it means that the center point of cell1 has moved to the outer layer of the original blade solid. At this time, the network is marked as a sedimentary layer grid. If n*n2<0, it means that the center point of cell1 is still in the inner layer of the original blade solid.

[0072] If cell1 is determined to be a sedimentary layer grid, the search for the next layer of grid will continue with cell1 as the starting point until the line vector connecting the coordinates of the grid center point inside a solid domain and P(X, Y, Z) is less than 0.

[0073] The particle thermal conductivity is assigned to the grid marked as the deposition layer, while the other grids remain unchanged. Then, in the time step t = t0 + Δt, the updated grid is calculated to solve the effect of the deposition layer on the flow field and the heat conduction in the solid domain.

[0074] In specific implementation, as a preferred embodiment of the present invention, the original boundary surface temperature T face The calculation method is as follows:

[0075]

[0076] in, is the thermal conductivity of the sedimentary layer grid; is the distance from the center point of the sedimentary layer grid to the surface; is the temperature of the center point of the sedimentary layer grid; is the thermal conductivity of the raw material grid; is the distance from the center point of the raw material grid to the surface; T2 is the temperature of the center point of the raw material grid.

[0077] The particle release process is repeated until t reaches the total time T of the numerical prediction simulation.

[0078] Example

[0079] like Figure 2 As shown, in this embodiment, a particle plate impact model and a grid schematic diagram are designed. Figure 2 (a) is a cross-section of the grid, and (b) is a schematic diagram of the structure of the fluid-solid coupling boundary layer.

[0080] The transient thermal coupling temperature of the wall is predicted by the method of the present invention. The particle deposition characteristics on the target surface after the flat plate impact and the target back surface temperature are compared with the experimental and calculated results. Figure 3 shown. Figure 3(a) is the particle deposition height, (b) is the experimental result after the experiment using the method of the present invention, (c) is a schematic diagram of the thermal conductivity of the deposition layer, and (d) is the actual calculation result. From the comparison of the experimental results and the actual calculation results in Figure (3), it can be seen that the method of the present invention can predict the wall temperature relatively accurately.

[0081] like Figure 4 As shown, the method of the present invention is used to predict the average temperature within the range of circles of different diameters on the back of the target just below the center line, and the experimental results are compared with the actual temperature values. The prediction results of the method of the present invention have a small error.

[0082] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the above embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for predicting the temperature of grid deformation and transient thermal coupling between particle deposition layer and wall surface, characterized in that: include: Establish the fluid domain and solid domain of the geometric model, establish the connection between the coincident surfaces of the fluid domain and the solid domain, and mesh the fluid domain and the solid domain respectively; Set different time steps for the fluid domain and solid domain to simulate and solve the flow equations, energy equations and turbulence equations in the unsteady flow field; Initialize using the user-defined function UDF and record the initial boundary surface mesh center coordinate position P (X, Y, Z) and normal n (x, y, z) direction; Release particles in the fluid domain and determine their motion trajectories by tracking the particle trajectories using the Lagrangian method; Set UDF custom boundary conditions and determine the particle deposition state through the particle deposition model; When particles are deposited on the solid surface, the physical properties of the particles are recorded in the corresponding grid cells and the height of the deposited layer is calculated; The sediment is evenly distributed to each node of the grid unit. Based on the dynamic mesh technology, the grid nodes at the interface between the fluid domain and the solid domain are moved simultaneously to achieve simultaneous boundary deformation of the fluid domain and the solid domain. Mark the sediment layer mesh after boundary deformation in the solid domain, solve the fluid flow and heat transfer in the fluid domain and solid domain, and obtain the temperature of the sediment and solid domain; The UDF is used to traverse the grid, extract and calculate the temperature of the grid where the sediment layer contacts the original solid, and solve the original boundary surface temperature under the sediment layer.

2. The method for predicting temperature of grid deformation and transient thermal coupling of particle deposition layer and wall according to claim 1, characterized in that: The overlapping surfaces between the fluid domain and the solid domain are connected through shared topology. The grid division adopts unstructured grid. The grid nodes at the junction of the fluid domain and the solid domain are completely corresponding, and boundary layer encryption is set on both sides of the junction.

3. The method for predicting temperature of grid deformation and transient thermal coupling of particle deposition layer and wall according to claim 1, characterized in that: The releasing of particles in the fluid domain and determining particle motion trajectories by tracking particle trajectories using the Lagrangian method specifically include: Particles are released into the fluid domain through the discrete phase model (DPM) and their migration trajectories are calculated. The Lagrangian discrete phase model is used to track the motion trajectories of the particles, and the Navier-Stokes equations are solved to treat the gas as a continuous phase. When calculating the migration trajectories, the drag, gravity, Saffman force, pressure gradient force, thermophoretic force, added mass force and Brownian force of the particles are taken into account, and the mutual influence between turbulence and particles is considered through a two-way coupling method.

4. The method for predicting temperature of grid deformation and transient thermal coupling of particle deposition layer and wall according to claim 1, characterized in that: Determining the deposition state of particles using a particle deposition model specifically includes: The particle deposition model calculates the softening temperature T according to the particle composition. s , and the critical viscosity μ corresponding to the softening temperature critical ; When the calculated viscosity reaches or falls below the critical viscosity, or the particle T p When the temperature reaches or exceeds the softening temperature, the adhesion probability of the particles is assumed to be 1, otherwise, the adhesion probability P s (T p ) is calculated as follows: Among them, μ Tp is the viscosity of the particle, μ soft is the particle viscosity corresponding to the particle softening temperature.

5. The method for predicting grid deformation and transient thermal coupling temperature of particle deposition layer and wall according to claim 1, characterized in that: The calculation method of the sediment layer height is as follows: Where n(x) represents the particle size in area A. cell The amount of sediment in the grid cell, m particle represents the mass of a single particle and ρ particle represents the density of a single particle; C represents the magnification factor.

6. The method for predicting grid deformation and transient thermal coupling temperature of particle deposition layer and wall according to claim 1, characterized in that: The dynamic mesh technology realizes the movement of mesh nodes on the interface between the fluid domain and the solid domain after particle deposition through mesh smoothing and reconstruction. The distance and direction of the mesh node movement are controlled by a user-defined function (UDF).

7. The method for predicting grid deformation and transient thermal coupling temperature of particle deposition layer and wall according to claim 1, characterized in that: The original boundary surface temperature T under the sediment layer face The calculation method is as follows: in, is the thermal conductivity of the sedimentary layer grid; is the distance from the center point of the sedimentary layer grid to the surface; is the temperature of the center point of the sedimentary layer grid; is the thermal conductivity of the raw material grid; is the distance from the center point of the raw material grid to the surface; T2 is the temperature of the center point of the raw material grid.

Images