CN108804744A - The method for numerical simulation of suspension material atomization film forming - Google Patents
The method for numerical simulation of suspension material atomization film forming Download PDFInfo
- Publication number
- CN108804744A CN108804744A CN201810318556.3A CN201810318556A CN108804744A CN 108804744 A CN108804744 A CN 108804744A CN 201810318556 A CN201810318556 A CN 201810318556A CN 108804744 A CN108804744 A CN 108804744A
- Authority
- CN
- China
- Prior art keywords
- liquid
- atomizing disk
- numerical simulation
- phases
- phase
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Application Of Or Painting With Fluid Materials (AREA)
Abstract
The present invention provides a kind of method for numerical simulation of suspension material atomization film forming, and the relationship of thickness of liquid film and each parameter is obtained by dimensional analysis and linear least-squares regression analysis.Liquid jet continuously impinges upon on atomizing disk, and to form film, which leaves atomizing disk to liquid jet along tangential direction, and is crushed and forms ligament or drop for extension in atomization panel surface.Numerical computations are carried out to the process using computational fluid dynamics software FLUENT, high cost and wastage of material caused by avoiding experiment or blindly designing have certain directive significance to the atomization mechanism of feed liquid.
Description
Technical field
The present invention relates to the method for numerical simulation to form a film on suspension material atomizing disk.
Background technology
Atomizing disk is widely used in producing in spraying, drop, particle and powder.In atomization process, liquid is formed on disk
Film, liquid film are broken into ligament or drop after leaving plate edge.Atomizing disk is widely used in spray drying device, according to this method
Application, these generated drops can obtain various rulers with hot setting to generate particle or powder using this method
Very little drop or powder.The size of drop depends on being formed in the fluid dynamic of liquid film on disk, and with the direct phase of thickness of liquid film
It closes.
It is extremely difficult to the measurement of thickness of liquid film and liquid-drop diameter since feed liquid atomization process is sufficiently complex, and large-scale work
Industry makeup, which sets and can not directly make a table apparatus, to be for experiment, and is summarized mostly according to the experience of engineering staff and forefathers at present
Semiempirical formula is predicted so that the spray drying device that actual design goes out always has such or such problem.Use calculating
Machine software emulates atomization process, and the guidance of the atomization mechanism of feed liquid is had a very important significance.
Invention content
In view of the above-mentioned deficiencies in the prior art, it is an object of the present invention to which the present invention provides the numbers of suspension material atomization film forming
It is worth analogy method.The present invention caused by also avoiding experiment to a certain extent using the method for numerical simulation or blindly designing it is high at
Sheet and technical risk, and have important directive significance for the atomization mechanism of feed liquid.
In order to achieve the above object, the present invention is achieved by the following technical programs：
A kind of method for numerical simulation of suspension material atomization film forming procedure, concrete operation step are as follows：
(1), using the 3-D geometric model figure of three-dimensional drawing Software on Drawing atomizing disk：
Atomizing disk threedimensional model is built using the modeling module (Design Model) in ANSYS Workbench
Mould；
(2), it establishes the finite element model of atomizing disk and is solved, concrete principle is as follows：
A, the assumed condition of numerical simulation is set；
A) do not consider to conduct heat between gas-liquid；
B) liquid is toppled over continuously along the central shaft of disk；
C) in order to reduce calculation amount, nearby the limited part of air is included atomizing disk；
B, it is modeled using Euler-Euler's method, gas-liquid two-phase is used as continuous phase processor, in order to describe multiphase
Between mutually run through and continuous Multiphase Flow, it is proposed that phase volume fraction concept uses αqIndicate the volume fraction shared by q phases,
Each mutually to meet quality and momentum conservation equation alone, the equation of satisfaction is as follows：
The mass-conservation equation of q phases is：
The momentum conservation equation of q phases is：
In formula, μqAnd λqExpression is shearing and the bulk viscosity (Pas) of q phases, F respectivelyqFor external volume power, Flift,qFor
Lift, FVm,qFor virtual mass power, RpqInteraction force between phase, p are the pressure that all phases are shared；
C, using k- ε turbulence models, tubulence energy k and tubulence energy dissipative shock wave ε formula are as follows：
In formula, μ is molecular viscosity (Pas), μtFor turbulent viscosity (Pas), PkIt is turbulent shear output item [kg/ (m
s3)], Cε1、Cε2、σk、σεFor constant, respectively 1.44,1.92,1,1.3；
D, in the FLUENT of ANSYS Workbench import atomizing disk 3-D geometric model, step A, B, C,
Assuming that and the rational computational domain of foundation in theoretical foundation and physical model, parameters are set, and simulation calculates the thickness of liquid film
Angle value obtains the relationship of thickness of liquid film and each parameter by dimensional analysis and linear least-squares regression analysis；
E, atomization test model is designed, and the numerical simulation result of step D is compared and analyzed with experimental result, is verified
The applicability of method for numerical simulation.
The present invention compared with prior art, have following obvious prominent substantive distinguishing features and notable technology into
Step：
Present invention employs the method for numerical simulation that a kind of suspension material is atomized film forming procedure, using calculating fluid dynamic
It learns software FLUENT and numerical simulation is carried out to suspension material atomization film forming procedure, the present invention is also existed using the method for numerical simulation
High cost and technical risk caused by avoiding experiment to a certain extent or blindly designing, and have weight for the atomization mechanism of feed liquid
The directive significance wanted.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, lower mask body is introduced the present invention and is implemented
Technical solution in example, it is clear that described embodiments are some of the embodiments of the present invention, instead of all the embodiments.Base
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts it is all its
His embodiment, shall fall within the protection scope of the present invention.
Embodiment one：
The method for numerical simulation of this suspension material atomization film forming, concrete operation step are as follows：
(1), using the 3-D geometric model figure of three-dimensional drawing Software on Drawing atomizing disk：
Atomizing disk threedimensional model is modeled using the modeling module (Design Model) in ANSYS Workbench
(2), it establishes the finite element model of atomizing disk and is solved, the specific method is as follows：
A, the assumed condition of numerical simulation is set；
A) do not consider to conduct heat between gas-liquid；
B) liquid is toppled over continuously along the central shaft of disk；
C) in order to reduce calculation amount, nearby the limited part of air is included atomizing disk；
B, it is modeled using Euler-Euler's method, gas-liquid two-phase is used as continuous phase processor, in order to describe multiphase
Between mutually run through and continuous Multiphase Flow, it is proposed that phase volume fraction concept uses αqIndicate the volume fraction shared by q phases,
Each mutually to meet quality and momentum conservation equation alone, the equation of satisfaction is as follows：
The mass-conservation equation of q phases is：
The momentum conservation equation of q phases is：
In formula, μqAnd λqExpression is shearing and the bulk viscosity (Pas) of q phases, F respectivelyqFor external volume power, Flift,qFor
Lift, FVm,qFor virtual mass power, RpqInteraction force between phase, p are the pressure that all phases are shared；
C, using k- ε turbulence models, tubulence energy k and tubulence energy dissipative shock wave ε formula are as follows：
In formula, μ is molecular viscosity (Pas), μtFor turbulent viscosity (Pas), PkIt is turbulent shear output item [kg/ (m
s3)], Cε1、Cε2、σk、σεFor constant, respectively 1.44,1.92,1,1.3；
D, in the FLUENT of ANSYS Workbench import atomizing disk 3-D geometric model, step A, B, C,
Assuming that and the rational computational domain of foundation in theoretical foundation and physical model, parameters are set, and simulation calculates the thickness of liquid film
Angle value obtains the relationship of thickness of liquid film and each parameter by dimensional analysis and linear least-squares regression analysis；
E, atomization test model is designed, and the numerical simulation result of step D is compared and analyzed with experimental result, is verified
The applicability of method for numerical simulation
Embodiment two：
The present embodiment specifically introduces the method for numerical simulation that molten metal prepares the atomization process of shot, and this method is at least wrapped
Include following steps：
First, atomizing disk threedimensional model is modeled using the modeling module (DM) in ANSYS Workbench：
Diameter 75mm, rotating speed 1500rpm, feed liquid flow 2.5kg/min, density 2590kg/m3, viscosity 0.7Pas.Choosing
With FLUENT software fluid volume function models, to capture the interface between liquid phase and gas phase.New material is created in FLUENT,
Relevant parameter is inputted according to the physical property of material, i.e., molten metal is modeled.
Secondly, establish the finite element model of atomizing disk and solved, it should be noted that the present invention in for based on
The method for numerical simulation of FLUENT in ANSYS Workbench, basic skills and step and traditional method are similar to step,
It will not be described here, and the hypothesis and theoretical foundation of the present invention is described below.
First, molten metal is atomized film forming procedure complexity, and method for numerical simulation reproduces melt atomization film forming procedure ratio completely
It is more difficult, it is not easy to obtain desired result.Therefore, the numerical simulation of molten metal of the invention atomization film forming procedure is based primarily upon
It is assumed hereinafter that：
(1) do not consider to conduct heat between gas-liquid.
(2) liquid is toppled over continuously along the central shaft of disk.
Second, gas-liquid two-phase is modeled using Euler-Euler's method, gas-liquid two-phase is used as continuous phase processor, is
It can describe mutually to run through between multiphase and continuous Multiphase Flow, it is proposed that phase volume fraction concept uses αqIndicate q phases institute
The volume fraction accounted for, each mutually to meet quality and momentum conservation equation alone, the equation of satisfaction is as follows：
The mass-conservation equation of q phases is：
The momentum conservation equation of q phases is：
In formula, μqAnd λqExpression is shearing and the bulk viscosity (Pas) of q phases, F respectivelyqFor external volume power, Flift,qFor
Lift, FVm,qFor virtual mass power, RpqInteraction force between phase, p are the pressure that all phases are shared；
Third, using k- ε turbulence models, tubulence energy k and tubulence energy dissipative shock wave ε formula are as follows：
In formula, μ is molecular viscosity (Pas), μtFor turbulent viscosity (Pas), PkIt is turbulent shear output item [kg/ (m
s3)], Cε1、Cε2、σk、σεFor constant, respectively 1.44,1.92,1,1.3.
The 3-D geometric model that atomizing disk is imported in the FLUENT of ANSYS Workbench, in above hypothesis and reason
The rational computational domain of foundation on the basis of and physical model, are arranged parameters, and simulation calculates the thickness value of liquid film, passes through
Dimensional analysis and linear least-squares regression analysis obtain the relationship of thickness of liquid film and each parameter；Devise atomization film forming experiment
Model, and above-mentioned numerical simulation result is compared and analyzed with experimental result, demonstrate the applicability of method for numerical simulation.
It should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although ginseng
According to previous embodiment, invention is explained in detail, it will be understood by those of ordinary skill in the art that：It still can be with
Technical scheme described in the above embodiments is modified or equivalent replacement of some of the technical features；And
These modifications or replacements, the spirit and model of various embodiments of the present invention technical solution that it does not separate the essence of the corresponding technical solution
It encloses.
Claims (1)
1. a kind of method for numerical simulation of suspension material atomization film forming, which is characterized in that concrete operation step is as follows：
(1), using the 3-D geometric model figure of three-dimensional drawing Software on Drawing atomizing disk：
Atomizing disk threedimensional model is modeled using the modeling module in ANSYS Workbench
(2), it establishes the finite element model of atomizing disk and is solved, the specific method is as follows：
A, the assumed condition of numerical simulation is set；
A) do not consider to conduct heat between gas-liquid；
B) liquid is toppled over continuously along the central shaft of disk；
C) in order to reduce calculation amount, nearby the limited part of air is included atomizing disk；
B, it is modeled using Euler-Euler's method, gas-liquid two-phase is used as continuous phase processor, in order to describe between multiphase
Mutually run through and continuous Multiphase Flow, it is proposed that phase volume fraction concept uses αqIndicate the volume fraction shared by q phases, each phase
Quality and momentum conservation equation can be met alone, the equation of satisfaction is as follows：
The mass-conservation equation of q phases is：
The momentum conservation equation of q phases is：
In formula, μqAnd λqExpression is shearing and the bulk viscosity (Pas) of q phases, F respectivelyqFor external volume power, Flift,qTo rise
Power, FVm,qFor virtual mass power, RpqInteraction force between phase, p are the pressure that all phases are shared；
C, using k- ε turbulence models, tubulence energy k and tubulence energy dissipative shock wave ε formula are as follows：
In formula, μ is molecular viscosity (Pas), μtFor turbulent viscosity (Pas), PkIt is turbulent shear output item [kg/ (m
s3)], Cε1、Cε2、σk、σεFor constant, respectively 1.44,1.92,1,1.3；
D, the 3-D geometric model that atomizing disk is imported in the FLUENT of ANSYS Workbench, in step A, B, C, hypothesis
And the rational computational domain of foundation in theoretical foundation and physical model, setting parameters, simulation calculate the thickness value of liquid film,
The relationship of thickness of liquid film and each parameter is obtained by dimensional analysis and linear least-squares regression analysis；
E, atomization test model is designed, and the numerical simulation result of step D is compared and analyzed with experimental result, demonstrates number
It is worth the applicability of analogy method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810318556.3A CN108804744B (en) | 2018-04-11 | 2018-04-11 | Numerical simulation method for atomization film formation of suspension material |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810318556.3A CN108804744B (en) | 2018-04-11 | 2018-04-11 | Numerical simulation method for atomization film formation of suspension material |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108804744A true CN108804744A (en) | 2018-11-13 |
| CN108804744B CN108804744B (en) | 2022-07-12 |
Family
ID=64094883
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810318556.3A Active CN108804744B (en) | 2018-04-11 | 2018-04-11 | Numerical simulation method for atomization film formation of suspension material |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108804744B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110414141A (en) * | 2019-07-30 | 2019-11-05 | 辽宁工程技术大学 | Drop during compressible fluids Transonic Flowss is atomized Three-dimensional Numerical Simulation Method |
| CN112329169A (en) * | 2020-11-03 | 2021-02-05 | 华南农业大学 | Numerical simulation analysis method for flow and heat transfer process of hot air drum type phoenix Dancong tea green removing machine |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106528971A (en) * | 2016-10-31 | 2017-03-22 | 烟台科力博睿地震防护科技有限公司 | A method and a system for determining the relationship between the size of a throat and the thickness of a liquid film of a crude oil refining nozzle |
| CN107038284A (en) * | 2017-03-20 | 2017-08-11 | 上海大学 | Multi-cavity rotary furnace and the method for numerical simulation for carrying out catalyst granules heating |
-
2018
- 2018-04-11 CN CN201810318556.3A patent/CN108804744B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106528971A (en) * | 2016-10-31 | 2017-03-22 | 烟台科力博睿地震防护科技有限公司 | A method and a system for determining the relationship between the size of a throat and the thickness of a liquid film of a crude oil refining nozzle |
| CN107038284A (en) * | 2017-03-20 | 2017-08-11 | 上海大学 | Multi-cavity rotary furnace and the method for numerical simulation for carrying out catalyst granules heating |
Non-Patent Citations (7)
| Title |
|---|
| NARAYANA P. RAYAPATI等: "Eulerian multiphase population balance model of atomizing, swirling flows", 《INTERNATIONAL JOURNAL OF SPRAY AND COMBUSTION DYNAMICS》 * |
| YUAN ZHU等: "Mechanism study and numerical simulation of Uranium nitriding induced by high energy laser", 《E3S WEB OF CONFERENCES》 * |
| 徐让书等: "离心喷嘴内部流动与液膜初级破碎的耦合模拟", 《沈阳工业大学学报》 * |
| 杨国华等: "基于树形自适应网格的旋流液膜雾化过程仿真", 《推进技术》 * |
| 毛传林等: "气泡雾化喷嘴雾化射流场性能仿真", 《计算机仿真》 * |
| 赵辉等: "基于铀的激光氮化机理研究和数值模拟", 《材料导报》 * |
| 陈雁等: "空气喷涂平面成膜的双流体模型模拟", 《后勤工程学院学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110414141A (en) * | 2019-07-30 | 2019-11-05 | 辽宁工程技术大学 | Drop during compressible fluids Transonic Flowss is atomized Three-dimensional Numerical Simulation Method |
| CN112329169A (en) * | 2020-11-03 | 2021-02-05 | 华南农业大学 | Numerical simulation analysis method for flow and heat transfer process of hot air drum type phoenix Dancong tea green removing machine |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108804744B (en) | 2022-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108804744A (en) | The method for numerical simulation of suspension material atomization film forming | |
| Morales-Higa et al. | Ladle shroud as a flow control device for tundish operations | |
| Vu et al. | Numerical investigation of dynamic behavior of a compound drop in shear flow | |
| Abbas et al. | Numerical simulation of electrohydrodynamic jet and printing micro-structures on flexible substrate | |
| Li et al. | Effects of different nozzle materials on atomization results via CFD simulation | |
| Liu et al. | Mathematical modeling of multi-sized argon gas bubbles motion and its impact on melt flow in continuous casting mold of steel | |
| Zhang | Mathematical simulation of fluid flow in gas-stirred liquid systems | |
| Wang et al. | Motion of single bubble and interactions between two bubbles in liquid steel | |
| Ge et al. | Gas-liquid-solid swirling flow polishing and bubble collapse impact characteristics | |
| CN108897902A (en) | The method for numerical simulation of material evaporation in spray drying tower | |
| Yu et al. | CFD simulation and experimental studies of suspension flow field in ultrasonic polishing | |
| Kwiatoń et al. | Multiphase modeling of droplet fall using Euler-Euler approach | |
| Kowitwarangkul et al. | CFD Simulation of Molten Steel Flow with Isothermal Condition in the Continuous Casting Tundish | |
| Minghao et al. | Numerical simulation of three-dimensional transient flow characteristics for a dual-fluid atomizer | |
| Menshutina et al. | CFD analysis of the dispersed phase behavior for micropowders production via spray drying and ultrasonic atomization | |
| Zhou et al. | Study on the Cavity Forming Induced by a Gas Jet Impinging on a Liquid Surface Based on a Deformed Mesh Method | |
| Duan et al. | Numerical simulation and experimental study of the characteristics of packing feature size on liquid flow in a rotating packed bed | |
| Kotapati et al. | Time-accurate three-dimensional simulations of synthetic jets in quiescent air | |
| Zhou et al. | Numerical investigation of a gas–solid turbulent jet flow with Reynolds number of 4500 using lattice Boltzmann method | |
| Yu et al. | Effect of nozzle inlet parameters on the dry granulation atomization process of Si3N4 ceramic bearing balls | |
| Cui et al. | 3-D CFD validation of an axisymmetric jet in cross-flow (NASA Langley Workshop Validation: Case 2) | |
| Guo et al. | Numerical study on the formation and solidification of LMPA microdroplet in a microfluidic device | |
| Liu et al. | A new general equation of mean particle size for different atomization processes | |
| Morenko | Laminar incompressible viscous flow past a cylinder performing rotary oscillations | |
| Aniszewski et al. | Numerical simulations of planar jets stripping of liquid coatings |
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 |