[go: up one dir, main page]

CN119940051B - Prediction method of particle distribution characteristics inside valve under wet solid conditions - Google Patents

Prediction method of particle distribution characteristics inside valve under wet solid conditions

Info

Publication number
CN119940051B
CN119940051B CN202411932905.2A CN202411932905A CN119940051B CN 119940051 B CN119940051 B CN 119940051B CN 202411932905 A CN202411932905 A CN 202411932905A CN 119940051 B CN119940051 B CN 119940051B
Authority
CN
China
Prior art keywords
particle
valve
particles
tab
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202411932905.2A
Other languages
Chinese (zh)
Other versions
CN119940051A (en
Inventor
林哲
乐琦
明友
陶俊宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Machinery Key Core Basic Component Innovation Center Anhui Co ltd
Hefei General Machinery Research Institute Co Ltd
Zhejiang Sci Tech University ZSTU
Original Assignee
General Machinery Key Core Basic Component Innovation Center Anhui Co ltd
Hefei General Machinery Research Institute Co Ltd
Zhejiang Sci Tech University ZSTU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by General Machinery Key Core Basic Component Innovation Center Anhui Co ltd, Hefei General Machinery Research Institute Co Ltd, Zhejiang Sci Tech University ZSTU filed Critical General Machinery Key Core Basic Component Innovation Center Anhui Co ltd
Priority to CN202411932905.2A priority Critical patent/CN119940051B/en
Publication of CN119940051A publication Critical patent/CN119940051A/en
Application granted granted Critical
Publication of CN119940051B publication Critical patent/CN119940051B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及数值模拟技术领域,公开了一种湿气固条件下阀门内部颗粒分布特征预测方法。方法包括建立阀门三维模型,对模型进行网格划分,结合湿气固条件下颗粒碰撞反弹特性,修正切向恢复系数和法向恢复系数,利用修正后的碰撞反弹模型开展数值模拟,获取模拟数据。通过数据提取和图像处理,生成阀门内部流动分布特征图和颗粒分布特征图,依据特征图判断颗粒积聚区域,并提取颗粒数密度,量化分析阀门内部高浓度颗粒分布。基于量化分析结果,提取颗粒捕集率参数,获得湿气固条件下由壁面液膜引起的颗粒捕集特征。本方法能够有效预测阀门内部颗粒分布及积聚特性,为阀门设计与优化提供依据。

The present invention relates to the field of numerical simulation technology, and discloses a method for predicting the particle distribution characteristics inside a valve under wet solid conditions. The method comprises establishing a three-dimensional model of the valve, meshing the model, combining the particle collision and rebound characteristics under wet solid conditions, correcting the tangential restitution coefficient and the normal restitution coefficient, and using the corrected collision and rebound model to carry out numerical simulation and obtain simulation data. Through data extraction and image processing, a flow distribution characteristic map and a particle distribution characteristic map inside the valve are generated, the particle accumulation area is judged based on the characteristic map, and the particle number density is extracted to quantitatively analyze the distribution of high-concentration particles inside the valve. Based on the quantitative analysis results, the particle capture rate parameter is extracted to obtain the particle capture characteristics caused by the wall liquid film under wet solid conditions. This method can effectively predict the particle distribution and accumulation characteristics inside the valve, and provide a basis for valve design and optimization.

Description

Valve internal particle distribution characteristic prediction method under wet gas-solid condition
Technical Field
The application relates to the technical field of numerical simulation, in particular to a valve internal particle distribution characteristic prediction method under a wet gas-solid condition.
Background
The valve is used for regulating and controlling medium flow in industrial production, and is widely applied to the fields of chemical industry, petroleum, electric power and the like. With the development of industrial technology, the application of gas-solid two-phase flow is increasing, for example, when industrial gas containing solid particles is conveyed under high temperature and high pressure, the particles may be accumulated or aggregated due to the wetting action of the wall surface, and the stability of fluid flow and valve performance are affected.
Currently, two-phase flow under wet solids conditions is under relatively limited investigation. The existing research is mainly aimed at dry environments, and the influence of a liquid film on particle collision and aggregation behaviors is difficult to accurately describe by a related model, so that the prediction accuracy of particle distribution characteristics is insufficient. In addition, the existing numerical simulation method has low applicability under the moisture-solid condition, and is difficult to meet the requirements of valve design and maintenance on particle distribution characteristic prediction.
Therefore, research is conducted on particle distribution characteristics under the moisture-solid condition, effective data support can be provided for optimizing valve performance, and the method has important significance in improving the efficiency and safety of industrial processes.
Disclosure of Invention
The embodiment of the application aims to provide a method for predicting the internal particle distribution characteristics of a valve under a wet gas-solid condition, so as to solve the problem that the collision rebound and accumulation behaviors of particles are difficult to accurately predict under the wet gas-solid condition in the prior art, and particularly, the prediction accuracy of the particle accumulation occurrence area and the capture rate is limited under the influence of a wall liquid film. The method effectively makes up the defects and provides a reliable basis for the optimal design of the valve.
According to the embodiment of the application, a method for predicting the distribution characteristics of particles in a valve under a wet gas-solid condition is provided, which comprises the following steps:
s1, establishing a valve three-dimensional model;
s2, carrying out grid division on the valve three-dimensional model;
S3, correcting tangential recovery coefficients and normal recovery coefficients in the particle collision rebound model by combining the three-dimensional valve model subjected to grid division with particle collision rebound characteristics under the moisture-solid condition, and carrying out numerical simulation on the corrected particle collision rebound model to obtain simulation data;
s4, carrying out data extraction on the simulation data, and then carrying out image processing to obtain a flow distribution characteristic diagram and a particle distribution characteristic diagram;
S5, judging a particle accumulation occurrence area according to the flow distribution characteristic diagram and the particle distribution characteristic diagram;
s6, extracting particle number density according to the particle accumulation occurrence area, and quantitatively analyzing the distribution of high-concentration particles in the valve under the moisture-solid condition by adopting a particle number density analysis method;
s7, extracting particle trapping rate parameters according to the quantitative analysis result, and further obtaining particle trapping characteristics caused by the wall liquid film under the moisture-solid condition;
The tangential recovery coefficient The expression of (2) is as follows:
;
The normal recovery coefficient The expression of (2) is as follows:
;
Wherein, the In order to fit the coefficients of the coefficients,At the critical angle of the beam,The impact angle is inclined for the particles.
The technical scheme provided by the embodiment of the application can comprise the following beneficial effects:
According to the embodiment, the three-dimensional model construction, the grid division and the corrected particle collision rebound model are adopted, the tangential recovery coefficient and the normal recovery coefficient are corrected by combining the particle motion characteristics under the moisture-solid condition, the technical problem that the traditional model under the dry wall surface can not accurately reflect the particle collision rebound and accumulation behaviors under the moisture-solid condition is solved, and the particle distribution prediction accuracy under the moisture-solid condition is improved. By introducing the quantitative analysis method of particle number density, the problem of insufficient analysis precision of particle accumulation areas and trapping rate in the prior art is solved, and accurate identification and quantitative analysis of the particle accumulation areas are further realized. Finally, the influence characteristics of the wetting wall on particle distribution are further obtained by extracting the particle trapping rate parameters, so that a more accurate technical basis is provided for valve design optimization.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application as claimed.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the application and together with the description, serve to explain the principles of the application.
FIG. 1 is a flow chart illustrating a method for predicting internal particle distribution characteristics of a valve under wet gas-solid conditions, according to an exemplary embodiment.
FIG. 2 is a three-dimensional model diagram of a ball valve, according to an example embodiment.
FIG. 3 is a flow chart illustrating meshing of a ball valve according to an exemplary embodiment.
FIG. 4 is a schematic diagram illustrating a ball valve internal flow path calculation region, according to an example embodiment.
FIG. 5 is a partial mesh encryption schematic of a three-dimensional model of a ball valve, according to an example embodiment.
Fig. 6 is a schematic diagram illustrating particle-wall collision bounce according to an exemplary embodiment.
Fig. 7 is a flow chart illustrating a particle impact rebound model boundary condition UDF file under wet conditions, according to an exemplary embodiment.
FIG. 8 is a grid cell arrangement illustrating a particle density solution according to an exemplary embodiment.
Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples do not represent all implementations consistent with the application.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used in this specification and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any or all possible combinations of one or more of the associated listed items.
FIG. 1 is a flow chart illustrating a method for predicting internal particle distribution characteristics of a valve under moisture solids conditions, according to an exemplary embodiment, as shown in FIG. 1, the method comprising the steps of:
in the implementation of step S1, a valve three-dimensional model is established;
Specifically, a three-dimensional model is built according to the actual geometric dimensions of the valve through operations such as stretching, array, mirror image, cutting, rotation and the like of a three-dimensional modeling software SOLIDOORKS, and as shown in fig. 2, the valve rod 1, the bracket 2, the valve cover 3, the ball 4, the valve seat 5 and the valve body 6 are respectively modeled and assembled into a complete ball valve model. And exporting the assembled ball valve three-dimensional model into STEP format, so that the subsequent model is conveniently imported into ANSYS SPACECLAIM software for model pretreatment.
In the implementation of the step S2, grid division is carried out on the valve three-dimensional model;
Specifically, the three-dimensional model of the valve is led into ANSYS SPACECLAIM software for model pretreatment, then grid division is carried out in FLUENT MESHING software, the number of grids is set, local grid encryption treatment is carried out on the valve cavity area inside the valve, and high-precision grids of the interaction area of fluid and particles are ensured so as to reflect the flow characteristics more accurately.
The flow chart of the grid division is shown in figure 3 and comprises the steps of importing a model and preprocessing, defining a grid area, generating a grid, checking the quality of the grid and the flow coefficientRefinement of grid and flow coefficientComparing flow coefficients, comparing minimum mesh sizeAnd particle diameterAnd solving the numerical value in the next step.
In the process of importing a model and preprocessing, in order to reduce the calculation amount of subsequent grid division and improve the grid precision, the valve three-dimensional model is simplified in ANSYS SPACECLAIM software, only parts related to a flow channel are reserved, irrelevant parts such as a bracket are deleted, holes and gaps irrelevant to the flow channel are filled, chamfers and fillets are deleted, concave-convex areas are filled and the like. In order to ensure that the gaseous medium in the ball valve pipeline and the solid particles in the medium are fully mixed, the length of the upstream pipeline of the ball valve is prolonged to 10 times of the diameter of the pipeline by a pulling function, the length of the downstream pipeline is prolonged to 10 times of the diameter of the pipeline,
In the defined grid area, the wall surface, the inlet and the outlet of the valve and the valve cavity area inside the valve to be grid encrypted are named, then the named inlet surface and outlet surface are selected through volume extraction in the preparation tab to extract the fluid area inside the valve and named fluid, and finally the conversion to FLUENT MESHING software is selected in the Workbench tab.
In the process of generating grids, the workflow of 'WATERTIGHT GEOMETRY' is selected in 'FLUENT MESHING' software, after a geometric model is imported, local dimensions are added to the inner surface of a named valve to be encrypted by the grids, so that the local grids are encrypted, then, a surface grid is generated, when the geometric structure is described, the 'geometric model is checked to be composed of only fluid areas without gaps', boundary conditions are updated to be a pressure outlet and a pressure inlet, a boundary layer is added, finally, a polyhedral grid 'poly-hexcore' generator grid is selected, a domain diagram is calculated for the internal flow channel of the ball valve in a dotted line frame of fig. 4, and fig. 5 is a schematic diagram for encrypting the local grids in the valve.
In checking the grid quality, clicking "execute grid check" in the "grid" tab, the grid quality parameters are visible on the underlying console, with emphasis on whether there is a negative volume and an orthogonal quality, and if there is a negative volume or an orthogonal quality below 0.5, it is necessary to return to the "model import and pre-process" step for re-modification.
After the grid quality inspection is completed, the flow coefficient is usedAs a target variable, grid independence verification starts. Specifically, an initial grid is used for carrying out numerical solution to obtain a flow coefficient. Then, by adjusting the minimum size of the gridTo refine the number of grids, and re-performing numerical solution to obtain flow coefficient. Comparing the errors of the two flow coefficients, if the errors are more than 5%, continuously refining grids and solving again, gradually optimizing grid division until the errors of the flow coefficients are less than 5%, and if the errors are less than 5%, indicating that the numerical solution is basically irrelevant to the grid division. Next, the minimum dimensions of the grid are comparedAnd particle diameter. If it isThe current grid division passes the independence verification and can be used for subsequent numerical solution, ifReturning to the "grid generation" step, the grid density is readjusted to ensure grid independence and accuracy of the particle distribution characteristics.
The grid quality of the exemplary embodiment in the disclosure finally reaches more than 0.7, meets the requirement of numerical simulation calculation, and can be solved in the next step.
In the implementation of step S3, the three-dimensional valve model after meshing is combined with the particle collision rebound characteristic under the moisture-solid condition, and the tangential recovery coefficient in the particle collision rebound model is calculatedAnd normal coefficient of restitutionCarrying out correction, carrying out numerical simulation on the corrected particle collision rebound model, and obtaining simulation data;
in particular, when the gas-solid two-phase flow calculation is directly performed, stable convergence of the calculation or ideal calculation accuracy is often difficult to achieve due to the strong coupling action of particles and gas in the flow field and the complexity of particle movement. Therefore, in order to improve the calculation efficiency and the accuracy of the result, the method is required to be carried out in two steps, namely, firstly, a stable gas flow field is established through numerical simulation, the basic characteristics of the gas flow are ensured to reach a stable state, then, on the basis, particles are introduced, and the movement track, the distribution characteristic and the interaction of the particles with the wall surface of the gas flow field are further calculated.
More specifically, firstly, a three-dimensional valve model after grid division is imported into FLUENT software, and unified grid units are set in a grid scaling tab, a single-phase gas flow field is established, steady state and gravity acceleration are set in a universal tab, energy and viscosity are opened in a model tab, and the flow field is setSST turbulence model and wall function, setting single-phase gas as 'air' in a 'materials' tab, setting fluid domain material as 'air' in a 'unit area condition' tab, defining boundary conditions of 'outlet' and 'inlet' as pressure in a 'boundary condition' tab, setting calculation method as Coupled and control method as defaults in a 'solving' tab, setting convergence residual value in a 'calculation monitoring' tab, setting an initialization method as 'standard initialization' in an 'initialization' tab, clicking an 'initialization' button to finish initialization setting, setting iteration times and time steps in a 'calculation setting' tab, and clicking 'start calculation'.
Next, by the experiment of collision and rebound between particles and wall with liquid film and analyzing the normal component and tangential component of the velocity of particles during collision, as shown in FIG. 6, the tangential recovery coefficient is obtainedAnd normal coefficient of restitutionAnd the tangential recovery coefficient and the normal recovery coefficient of the particle impact rebound model under the moisture-solid condition are used for correcting;
Finally, after the stable flow field is obtained after the calculation is completed, setting ' transient ' and ' gravity acceleration ' in a ' universal ' tab, opening ' discrete phase ' in a ' model ' tab, checking ' interaction with a continuous phase ', creating ' injection source ', selecting ' particle type ' as ' inertia ', selecting ' injection source type ' as ' file ', wherein ' file ' is an externally imported user-defined particle package file, selecting ' function ' in the ' user-defined ' tab, loading a reflection boundary condition UDF function of a particle impact rebound model under the condition of moisture solid containing correction, setting ' discrete phase boundary type ' in DPM of all wall surfaces in the ' boundary condition ' tab as ' user-defined ', setting ' discrete phase BC function ' as a user-defined file ' bc_reflection: libduf ', setting the iteration times and time step in a ' running calculation ' tab, clicking ' to start calculation, and obtaining data of CASE files and DATE files after the calculation.
By constructing the gas flow field in advance before particle calculation, the problems that two-phase flow is difficult to converge and the error is large in the traditional method when two-phase flow is directly calculated are effectively solved. The numerical calculation method of the embodiment of the application obviously improves the accuracy and the reliability of the result, and simultaneously provides a more reliable and innovative solution for the numerical simulation of the two-phase flow.
In numerical computation, the residual is one of the key indicators for measuring whether the computation converges, and the judgment condition for the termination of the computation generally includes the residual and the number of iteration steps. When the residual reaches a preset threshold or the iteration number reaches an upper limit, the calculation is automatically ended. In order to improve the accuracy of flow field simulation as much as possible, the residual error convergence standard is set to 10 -20 times and the iteration number is 2000 at this time, so as to ensure the reliability and accuracy of the calculation result.
A flow chart for compiling a reflection boundary condition UDF function of a particle impact rebound model under modified moisture-solid conditions is shown in fig. 7, and the specific steps are as follows:
(1) Initializing variables including the angle between the particle velocity and the normal vector of the wall Angle of reflectionNormal direction of particlesCritical speed of particleCoefficient of normal recoveryCoefficient of tangential recovery;
(2) Calculating a vector, judging whether the rotation condition is axisymmetric according to the value of rp_axi_ swirl, wherein rp_axi_ swirl is a global variable in Fluent and is used for identifying whether an axisymmetric rotation flow model is started, if so, calculating a three-dimensional vector to avoid numerical instability, and if not, directly using a two-dimensional vector;
(3) Checking the particle type, judging whether the particle is an inert particle, if yes, continuing the next processing, if not, skipping the particle processing, wherein the particle type is set to be inert in the exemplary embodiment of the application, so that the next processing can be directly performed;
(4) Comparing the normal velocity of the particles And critical speedThe particle velocity and the normal vector algorithm are multiplied by the velocityCalculating critical velocity of particles using empirical formula;
(5) If it isThe particle velocity is determined to be 0, and the particle velocity of 0 is set, meaning that the particles are trapped by the wall liquid film;
If it is Then entering wall reflection logic, firstly calculating a reflection angle, and calculating a tangential recovery coefficient according to the size of the reflection angle and a fitting coefficient and a formula obtained by combining experimentsAnd normal coefficient of restitutionFinally, the particle speed is readjusted by using two recovery coefficients, and the initial particle speed is updated to continue moving.
The impact rebound speed and angle of the particles and the wall surface with the liquid film are changed under the influence of the liquid film, and meanwhile, the existence of the liquid film is equivalent to the addition of a protective layer on the wall surface, so that the impact speed of the particles and the wall surface is reduced, and the movement track of the particles is changed.
Therefore, the fitting coefficient and the formula obtained by the experiment are obtained by collision and rebound experiments of particles and a wall surface with a liquid film, and the normal component and the tangential component of the velocity of the particles during collision are analyzed to obtain the tangential recovery coefficientAnd normal coefficient of restitutionAnd the method is used for writing a reflection boundary condition UDF function of the particle impact rebound model under the modified moisture-solid condition in the subsequent numerical simulation so as to simulate the particle impact rebound behavior more truly.
Through the obtained experimental data, tangential recovery coefficients can be obtained by fitting in Origin softwareAnd normal coefficient of restitutionAngle of impact with particle inclinationThe changing functional relation:
(1) Tangential recovery coefficient
(2) Normal coefficient of restitution
Wherein, the In order to fit the coefficients of the coefficients,Is a critical angle.
In an exemplary embodiment of the present application, the fitting coefficients after fitting by experimental data are:
In the implementation of step S4, the flow distribution feature map and the particle distribution feature map are obtained after the data extraction is performed on the analog data and then the image processing is performed.
Specifically, the calculated analog data is subjected to graphic processing, and a characteristic section is required to be created for displaying an internal flow field distribution characteristic diagram better. Selecting "cloud" in the "results" tab in FLUENT and clicking on "new face" to create an appropriate feature section, sequentially selecting "visibility" and "Velocity Magnitude" in the "coloring variables" tab, and then clicking on "save/display"; likewise, a cloud image is created again on the created feature section with respect to the internal Pressure distribution, and "Pressure" and "Static Pressure" are selected in sequence in the "coloring variable" tab, and then "save/display" is clicked;
the "coloring variables" in the "particle track" tab select "Particle Variables" and "Particle Velocity Magnitude" in sequence and then click "track" and "save/display" so that a map of the movement track of the interior particle can be obtained.
In order to better display the movement of particles in the valve or the distribution condition of particles trapped by the wetted wall surface, the valve grid and the particle track are displayed in a graph, the transparency of the valve grid is properly set, namely, the grids are checked in a scene option tab, the grid is Particle Velocity, the grid transparency is 60%, and finally the particle distribution characteristic graph is obtained.
In the implementation of step S5, the particle accumulation occurrence area is determined according to the flow distribution feature map and the particle distribution feature map.
Specifically, the flow distribution characteristic diagram and the particle distribution characteristic diagram are utilized to visually display the flow field condition in the valve, the interaction condition of particles and a wall liquid film under the moisture-solid condition is observed, and the locking speed is reduced to the area where the particles are located. By analyzing the particle mass concentration distribution and the volume concentration distribution of these regions, a specific region where particle accumulation occurs is labeled.
In the implementation of step S6, the particle number density is extracted according to the particle accumulation occurrence area, and the quantitative analysis of the high-concentration particle distribution inside the valve under the moisture-solid condition is performed by adopting a particle number density analysis method.
The particle number density is used as a characterization parameter, the particle number density of a particle accumulation occurrence area is obtained by calculating the equivalent particle number in each grid unit and summing up the probability fractions of each particle reaching the grid unit according to the particle number density result, the spatial distribution characteristics of the particles under the moisture-curing condition are further quantified according to the particle number density result, the quantification standard comprises the sum of the particle numbers in each grid unit and the distribution uniformity among the grid units, the concentration gradient and the accumulation characteristics of the particles on the spatial distribution are evaluated according to the spatial distribution characteristics of the particles under the moisture-curing condition by counting the equivalent number of the particles in the grid units and combining the probability distribution of the particles reaching the grid units, and the distribution range and the position of the particles with high concentration are further confirmed, wherein the concentration of the particles in the accumulation area is larger than that of the particles at an inlet. The particle number density is defined as the number of equivalent particles in a given grid (s×s), the grid cell arrangement to be used is shown in fig. 8, and the probability score calculation formula is as follows:
Wherein, the For each probability score that a particle falls on a respective grid,AndIs the coordinates of the particles and,AndFor the coordinates of the center points of the respective grids,AndAt the center point of two adjacent gridsShaft and method for producing the sameDistance in the axial direction.
By the method, the distribution density degree of the particles in different areas can be intuitively quantified. Areas of high particle count density are often closely related to particle accumulation and these areas are often critical sites for valve operation that are susceptible to wear or clogging. Based on analysis of particle number density, a reference basis can be provided for optimizing a flow field under a moisture-solid condition, and a new thought is provided for researching particle distribution characteristics.
In the implementation of step S7, the particle capturing rate parameter is extracted according to the result of the quantitative analysis, so as to further obtain the particle capturing characteristic caused by the wall liquid film under the moisture-solid condition.
Specifically, the particle trapping rate is defined as the trapping probability fraction of particles in the wall liquid film area, and the trapping characteristic of the particles in the valve is obtained by calculating the trapping rate difference of the particles in different wetting wall areas and constructing a particle trapping distribution model. The calculation formula of the particle trapping rate is as follows:
Wherein, the In order to achieve the particle trapping rate, the particle size distribution,In order for the number of particles to be trapped by the wall,The total particle number in the flow field is calculated.
The capture distribution characteristics of the particles in different wetted wall areas were analyzed by numerical simulation and data extraction. In certain specific areas inside the valve, the wall liquid film may lead to higher particle trapping rates due to flow characteristics or geometry effects. For these regions, potential particle accumulation regions can be identified by spatially distributed features of the trapping rate, thereby optimizing the valve cavity shape design and reducing the risk of particle deposition during operation.
As can be seen from the above embodiments, the embodiment of the present application adopts three-dimensional model construction, meshing and modified particle impact rebound model, and combines particle motion characteristics under moisture-solid condition to tangential recovery coefficientAnd normal coefficient of restitutionThe method solves the technical problem that the model under the traditional drying wall surface can not accurately reflect the collision rebound and accumulation behaviors of the moisture-solid condition on the particles, and further improves the accuracy of particle distribution prediction under the moisture-solid condition. By introducing the quantitative analysis method of particle number density, the problem of insufficient analysis precision of particle accumulation areas and trapping rate in the prior art is solved, and accurate identification and quantitative analysis of the particle accumulation areas are further realized. Finally, the influence characteristics of the wetting wall on particle distribution are further obtained by extracting the particle trapping rate parameters, so that a more accurate technical basis is provided for valve design optimization.
Other embodiments of the application will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure herein. This application is intended to cover any variations, uses, or adaptations of the application following, in general, the principles of the application and including such departures from the present disclosure as come within known or customary practice within the art to which the application pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the application being indicated by the following claims.
It is to be understood that the application is not limited to the precise arrangements and instrumentalities shown in the drawings, which have been described above, and that various modifications and changes may be effected without departing from the scope thereof. The scope of the application is limited only by the appended claims.

Claims (7)

1.一种湿气固条件下阀门内部颗粒分布特征预测方法,其特征在于,包括:1. A method for predicting particle distribution characteristics inside a valve under wet solid conditions, comprising: S1:建立阀门三维模型;S1: Establish a three-dimensional model of the valve; S2:对所述阀门三维模型进行网格划分;S2: Meshing the three-dimensional model of the valve; S3:将网格划分后的阀门三维模型结合湿气固条件下颗粒碰撞反弹特性,对颗粒碰撞反弹模型中的切向恢复系数和法向恢复系数进行修正,对修正后的颗粒碰撞反弹模型开展数值模拟,获得模拟数据;S3: Combine the meshed 3D valve model with the particle collision and rebound characteristics under wet solid conditions, modify the tangential restitution coefficient and normal restitution coefficient in the particle collision and rebound model, and perform numerical simulation on the modified particle collision and rebound model to obtain simulation data; S4:对所述模拟数据进行数据提取再经过图像处理,得到流动分布特征图和颗粒分布特征图;S4: extracting data from the simulation data and then performing image processing to obtain a flow distribution characteristic map and a particle distribution characteristic map; S5:根据所述流动分布特征图和颗粒分布特征图,判断颗粒积聚发生区域;S5: determining a particle accumulation region based on the flow distribution characteristic map and the particle distribution characteristic map; S6:根据所述颗粒积聚发生区域,提取颗粒数密度,采用颗粒数密度分析方法对湿气固条件下阀门内部高浓度颗粒分布进行量化分析;S6: extracting the particle number density according to the particle accumulation occurrence area, and using a particle number density analysis method to quantitatively analyze the distribution of high-concentration particles inside the valve under wet solid conditions; S7:根据所述量化分析的结果,提取颗粒捕集率参数,进一步获得湿气固条件下由壁面液膜造成的颗粒捕集特征;S7: extracting a particle capture rate parameter based on the results of the quantitative analysis, and further obtaining the particle capture characteristics caused by the wall liquid film under wet solid conditions; 所述切向恢复系数的表达式如下:The tangential restitution coefficient The expression is as follows: ; 所述法向恢复系数的表达式如下:The normal restitution coefficient The expression is as follows: ; 其中,为拟合系数,为临界角度,为颗粒倾斜冲击角度。in, is the fitting coefficient, is the critical angle, is the particle impact angle. 2.根据权利要求1所述的方法,其特征在于,对所述阀门三维模型进行网格划分,包括:2. The method according to claim 1, wherein meshing the valve three-dimensional model comprises: 将阀门的三维模型导入Ansys SpaceClaim软件中进行模型预处理,预处理后再导入至FLUENT MESHING软件中进行网格划分,设置网格数量并针对阀门内部的阀腔区域进行局部网格加密处理。The 3D model of the valve was imported into Ansys SpaceClaim software for model preprocessing. After preprocessing, it was imported into FLUENT MESHING software for meshing. The number of meshes was set and local mesh encryption was performed on the valve cavity area inside the valve. 3.根据权利要求1所述的方法,其特征在于,将网格划分后的阀门三维模型结合湿气固条件下颗粒碰撞反弹特性,对颗粒碰撞反弹模型中的切向恢复系数和法向恢复系数进行修正,对修正后的颗粒碰撞反弹模型开展数值模拟,获得模拟数据,包括:3. The method according to claim 1 is characterized in that the tangential restitution coefficient in the particle collision rebound model is calculated by combining the three-dimensional valve model after meshing with the particle collision rebound characteristics under wet solid conditions. and normal restitution coefficient Correction is made, and numerical simulation is carried out on the corrected particle collision and rebound model to obtain simulation data, including: (1)将网格划分后的阀门三维模型导入至FLUENT软件中,并在“网格缩放”选项卡中设置统一网格单位;建立气体流场,在“通用”选项卡中设置“稳态”及“重力加速度”;在“模型”选型卡中打开“能量”和“粘性”,并设置湍流模型和壁面函数;在“材料”选项卡中设置单相气体为“空气”;在“单元区域条件”选项卡中设置流体域材料为“空气”;在“边界条件”选项卡中定义“出口”、“进口”的边界条件;在“求解”选项卡中设置计算方法及控制方法;在“计算监控”选项卡中设置收敛残差数值;在“初始化”选项卡中设置成初始化方法为“标准初始化”,并点击“初始化”按钮完成初始化设置;在“计算设置”选项卡中设置迭代次数与时间步长,并点击“开始计算”;(1) Import the meshed valve 3D model into FLUENT software and set a uniform mesh unit in the "Mesh Scaling" tab; establish the gas flow field and set "Steady State" and "Gravity Acceleration" in the "General" tab; open "Energy" and "Viscosity" in the "Model" tab and set the turbulence model and wall function; set the single-phase gas to "Air" in the "Material" tab; set the fluid domain material to "Air" in the "Unit Region Condition" tab; define the boundary conditions of "Outlet" and "Inlet" in the "Boundary Condition" tab; set the calculation method and control method in the "Solution" tab; set the convergence residual value in the "Calculation Monitoring" tab; set the initialization method to "Standard Initialization" in the "Initialization" tab and click the "Initialize" button to complete the initialization settings; set the number of iterations and time step in the "Calculation Settings" tab and click "Start Calculation"; (2)通过颗粒与带液膜的壁面碰撞反弹实验,并分析颗粒碰撞时速度的法向分量和切向分量,得到切向恢复系数和法向恢复系数,并用于修正湿气固条件下的颗粒碰撞反弹模型的切向恢复系数和法向恢复系数;(2) Through the particle collision and rebound experiment with the wall with liquid film, the normal component and tangential component of the particle velocity during collision are analyzed to obtain the tangential restitution coefficient and normal restitution coefficient , and is used to correct the tangential restitution coefficient and normal restitution coefficient of the particle collision rebound model under wet solid conditions; (3)在上述计算完成后,继续再建立颗粒流场,在“通用”选项卡中设置“瞬态”及“重力加速度”;在“模型”选项卡中打开“离散相”,勾选“与连续相的交互”,并创建“喷射源”,选择“喷射源类型”为“file”,其中“file”为外部导入的用户自定义的颗粒包文件;在“用户自定义”选项卡中加载修正后的湿气固条件下的颗粒碰撞反弹模型的反射边界条件UDF函数,并设置“边界条件”选项卡中壁面的DPM相关参数为所加载的UDF函数;在“运行计算”选项卡中设置迭代次数与时间步长,并点击“开始计算”;最后,得到CASE文件和DATE文件的模拟数据。(3) After the above calculations are completed, continue to establish the particle flow field, set "Transient" and "Gravity Acceleration" in the "General" tab; open "Discrete Phase" in the "Model" tab, check "Interaction with Continuous Phase", and create a "Jet Source", select "Jet Source Type" as "file", where "file" is the user-defined particle package file imported externally; load the modified reflection boundary condition UDF function of the particle collision rebound model under wet gas-solid conditions in the "User Defined" tab, and set the DPM related parameters of the wall in the "Boundary Conditions" tab to the loaded UDF function; set the number of iterations and time step in the "Run Calculation" tab, and click "Start Calculation"; finally, obtain the simulation data of the CASE file and DATE file. 4.根据权利要求3所述的方法,其特征在于,对所述模拟数据进行数据提取后再经过图像处理,得到流动分布特征图和颗粒分布特征图,包括:4. The method according to claim 3, wherein the step of extracting the simulation data and then performing image processing to obtain a flow distribution characteristic map and a particle distribution characteristic map comprises: 在FLUENT软件中“结果”选项卡中选择“云图”,绘制阀门内部特征截面上的流动分布特征图,其包括了气体速度云图、气体压力云图,在“颗粒轨迹”选项卡中绘制湿润壁面上的颗粒分布特征图。In the FLUENT software, select "Cloud Map" in the "Results" tab to draw the flow distribution characteristic map on the characteristic section inside the valve, which includes the gas velocity cloud map and the gas pressure cloud map. In the "Particle Trajectory" tab, draw the particle distribution characteristic map on the wetted wall surface. 5.根据权利要求1所述方法,其特征在于,根据所述流动分布特征图和颗粒分布特征图,判断颗粒积聚发生区域,包括:5. The method according to claim 1, wherein determining the particle accumulation region based on the flow distribution characteristic map and the particle distribution characteristic map comprises: 通过所述流动分布特征图和颗粒分布特征图显示,观察在阀门内部被壁面液膜捕捉后速度为0的颗粒的位置区域,分析并标记颗粒位置,通过流场数值对比分析颗粒数量和密集程度,确定湿气固条件下颗粒积聚的区域。The flow distribution characteristic diagram and the particle distribution characteristic diagram are used to observe the position area of particles with a velocity of 0 after being captured by the wall liquid film inside the valve, analyze and mark the particle positions, analyze the number and density of particles through flow field numerical comparison, and determine the area where particles accumulate under wet-solid conditions. 6.根据权利要求5所述方法,其特征在于,根据所述颗粒积聚发生区域,提取颗粒数密度,采用颗粒数密度分析方法对湿气固条件下阀门内部高浓度颗粒分布进行量化分析,包括:6. The method according to claim 5, characterized in that the particle number density is extracted based on the particle accumulation occurrence area, and a particle number density analysis method is used to quantitatively analyze the distribution of high-concentration particles inside the valve under wet solid conditions, comprising: S61:引入颗粒数密度作为表征参数,通过计算每个网格单元内的当量颗粒数量,根据流场中每个颗粒到达网格单元的概率分数求和,得到颗粒积聚发生区域的颗粒数密度;S61: Introducing particle number density as a characterization parameter, by calculating the equivalent number of particles in each grid cell and summing the probability fractions of each particle reaching the grid cell in the flow field, the particle number density in the area where particle accumulation occurs is obtained; S62:根据颗粒数密度结果,进一步量化颗粒在湿气固条件下的空间分布特征,量化标准包括每个网格单元内颗粒数量的总和及其在网格单元间的分布均匀性;S62: Based on the particle number density results, further quantify the spatial distribution characteristics of particles under wet solid conditions. The quantification criteria include the total number of particles in each grid cell and the uniformity of their distribution among grid cells. S63:根据颗粒在湿气固条件下的空间分布特征,通过统计颗粒在网格单元内的当量数量,结合各网格单元中颗粒到达的概率分布,评估颗粒在空间分布上的浓度梯度和积聚特征,进一步确认高浓度颗粒的分布范围和位置,所述高浓度颗粒指积聚区域内颗粒的浓度大于进口处颗粒的浓度。S63: Based on the spatial distribution characteristics of particles under wet solid conditions, by counting the equivalent number of particles in the grid unit and combining the probability distribution of particle arrival in each grid unit, the concentration gradient and accumulation characteristics of the particles in the spatial distribution are evaluated to further confirm the distribution range and position of high-concentration particles, which refers to particles with a concentration in the accumulation area greater than the concentration at the inlet. 7.根据权利要求6所述方法,其特征在于,根据所述量化分析的结果,提取颗粒捕集率参数,进一步获得湿气固条件下由壁面液膜造成的颗粒捕集特征,包括:7. The method according to claim 6, wherein the particle capture rate parameter is extracted based on the results of the quantitative analysis, and the particle capture characteristics caused by the wall liquid film under wet solid conditions are further obtained, comprising: 根据所述量化分析的高浓度颗粒分布结果,进一步分析颗粒在湿润壁面上的分布情况,定义颗粒捕集率为颗粒在壁面液膜区域的捕集概率分数,通过计算颗粒在不同湿润壁面区域的捕集率差异,得到阀门内部颗粒的捕集特征,从而实现湿气固条件下阀门内部颗粒分布特征的预测。Based on the high-concentration particle distribution results of the quantitative analysis, the distribution of particles on the wetted wall surface is further analyzed, and the particle capture rate is defined as the capture probability fraction of particles in the liquid film area of the wall surface. By calculating the difference in capture rate of particles in different wetted wall areas, the capture characteristics of particles inside the valve are obtained, thereby realizing the prediction of the particle distribution characteristics inside the valve under wet-solid conditions.
CN202411932905.2A 2024-12-26 2024-12-26 Prediction method of particle distribution characteristics inside valve under wet solid conditions Active CN119940051B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411932905.2A CN119940051B (en) 2024-12-26 2024-12-26 Prediction method of particle distribution characteristics inside valve under wet solid conditions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411932905.2A CN119940051B (en) 2024-12-26 2024-12-26 Prediction method of particle distribution characteristics inside valve under wet solid conditions

Publications (2)

Publication Number Publication Date
CN119940051A CN119940051A (en) 2025-05-06
CN119940051B true CN119940051B (en) 2025-10-24

Family

ID=95534214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411932905.2A Active CN119940051B (en) 2024-12-26 2024-12-26 Prediction method of particle distribution characteristics inside valve under wet solid conditions

Country Status (1)

Country Link
CN (1) CN119940051B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102460527A (en) * 2009-05-01 2012-05-16 爱克斯崔里斯科技有限公司 Improvements to Particle Detectors
CN116884883A (en) * 2023-09-01 2023-10-13 山西创芯光电科技有限公司 A method for reducing air bubbles in the bottom filling of infrared detectors

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5964985A (en) * 1994-02-02 1999-10-12 Wootten; William A. Method and apparatus for converting coal to liquid hydrocarbons
DE102018010197A1 (en) * 2018-12-18 2020-06-18 GRID INVENT gGmbH Electronic element and electrically controlled display element

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102460527A (en) * 2009-05-01 2012-05-16 爱克斯崔里斯科技有限公司 Improvements to Particle Detectors
CN116884883A (en) * 2023-09-01 2023-10-13 山西创芯光电科技有限公司 A method for reducing air bubbles in the bottom filling of infrared detectors

Also Published As

Publication number Publication date
CN119940051A (en) 2025-05-06

Similar Documents

Publication Publication Date Title
Lyu et al. Fast and versatile fluid-solid coupling for turbulent flow simulation
CN113010977B (en) Gear oil injection lubrication optimization simulation method
Nirschl et al. Three-dimensional calculations of the simple shear flow around a single particle between two moving walls
WO2012092502A1 (en) Generating inviscid and viscous fluid-flow simulations over an aircraft surface
CN112069689A (en) Simulation method and system for fuel atomization characteristic of aircraft engine
CN114638131B (en) A crystal plasticity finite element modeling method based on electron backscatter diffraction analysis and Gmesh software
CN112434475B (en) Numerical simulation calculation result post-processing method for pressurized water nuclear reactor pressure vessel
Xia et al. Particle swarm optimization of aerodynamic shapes with nonuniform shape parameter–based radial basis function
CN118296990A (en) Modeling method for submerged boundary and grid treatment of aero-engine combustion chamber
CN1388444A (en) Simulation method of casting cavity filling process
CN119940051B (en) Prediction method of particle distribution characteristics inside valve under wet solid conditions
CN112685976A (en) Numerical simulation method of natural circulation system under motion condition
CN116824015A (en) PINNs-based flow field data generation method, device, equipment and medium
CN112257312B (en) A Parametric Modeling Method for the Erosion Model of Fine Flake Particles in Turbine Blade Materials
Cummings et al. Supersonic, turbulent flow computation and drag optimization for axisymmetric afterbodies
CN115563900A (en) Method for suppressing cavitation of cryogenic butterfly valve, method for verifying effectiveness, and electronic equipment
CN118313314B (en) Particle analog free incoming flow boundary processing method, device, equipment and medium
CN114330080A (en) Prediction method of aircraft surface-symmetric cross-basin flow field
CN120012649A (en) A numerical prediction method for dynamic characteristics of coarse particle solid-liquid two-phase flow in vertical lifting pipeline
CN113627099B (en) A method and device for predicting cavitation flow field
CN113505544A (en) Bicycle motion virtual numerical value wind tunnel system based on finite volume method
CN114048696B (en) Pneumatic simulation analysis method for flapping-rowing-torsion flapping wing aircraft
CN112182762B (en) A CFD Calculation Method and System of Aircraft Wake Vortex Based on Adaptive Grid
CN115964891A (en) Simulation analysis method and device for automobile fan and storage medium
CN119692254B (en) Molten pool width prediction method for 3D printing selective laser melting technology

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant