CN108804744B - Numerical simulation method for atomization film formation of suspension material - Google Patents

Numerical simulation method for atomization film formation of suspension material Download PDF

Info

Publication number
CN108804744B
CN108804744B CN201810318556.3A CN201810318556A CN108804744B CN 108804744 B CN108804744 B CN 108804744B CN 201810318556 A CN201810318556 A CN 201810318556A CN 108804744 B CN108804744 B CN 108804744B
Authority
CN
China
Prior art keywords
liquid
numerical simulation
phase
model
atomization
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
CN201810318556.3A
Other languages
Chinese (zh)
Other versions
CN108804744A (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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201810318556.3A priority Critical patent/CN108804744B/en
Publication of CN108804744A publication Critical patent/CN108804744A/en
Application granted granted Critical
Publication of CN108804744B publication Critical patent/CN108804744B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 invention provides a numerical simulation method for atomization film forming of suspension liquid materials, and the relation between the thickness of a liquid film and each parameter is obtained through dimensional analysis and a linear least square regression analysis method. The liquid jet continuously impinges on the atomizing disk, where it spreads to form a film, which exits the disk tangentially and breaks up to form ligaments or droplets. The numerical calculation is carried out on the process by adopting computational fluid dynamics software FLUENT, so that high cost and raw material waste caused by experiments or blind design are avoided, and certain guiding significance is provided for the atomization mechanism of the feed liquid.

Description

Numerical simulation method for atomization film formation of suspension material
Technical Field
The invention relates to a numerical simulation method for film forming on a suspension material atomizing disc.
Background
Atomizing disks are widely used in the production of sprays, droplets, granules and powders. During the atomization process, a liquid film forms on the disc, which breaks up into ligaments or droplets after leaving the edge of the disc. Atomizing disks are widely used in spray drying equipment and, depending on the application of the method, the droplets produced can be solidified at high temperatures to produce particles or powders, which can be used to obtain droplets or powders of various sizes. The size of the droplets depends on the hydrodynamic force of the liquid film formed on the disc and is directly related to the thickness of the liquid film.
Because the feed liquid atomization process is very complicated, the measurement of the liquid film thickness and the liquid drop diameter is extremely difficult, and a large-scale industrialized device cannot be directly manufactured into a device for testing, at present, the measurement is mostly predicted according to the experience of engineering personnel and a semi-empirical formula summarized by predecessors, so that the actually designed spray drying equipment always has the problems. The computer software is used for simulating the atomization process, and the method has very important significance for guiding the atomization mechanism of the feed liquid.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a numerical simulation method for atomization film formation of suspension materials. The invention also avoids high cost and technical risk caused by experiment or blind design to a certain extent by adopting a numerical simulation method, and has important guiding significance for the atomization mechanism of the feed liquid.
In order to achieve the purpose, the invention is realized by the following technical scheme:
a numerical simulation method for a suspension material atomization film forming process comprises the following specific operation steps:
(1) drawing a three-dimensional geometric model diagram of the atomizing disc by using three-dimensional drawing software:
modeling a three-dimensional Model of the atomizing disc by using a modeling module (Design Model) in ANSYS Workbench;
(2) establishing a finite element model of the atomizing disc and solving, wherein the specific principle is as follows:
A. setting a hypothetical condition of numerical simulation;
a) heat transfer between gas and liquid is not considered;
b) the liquid is poured continuously along the central axis of the tray;
c) to reduce the number of calculations, a limited portion of the air near the atomizing disk is included;
B. modeling is carried out by using an Euler-Euler method, gas-liquid two phases are treated as continuous phases, and in order to describe the multiphase flow which is mutually penetrated and continuous, a phase volume fraction concept is provided and alpha is usedqAnd (3) representing the volume fraction occupied by the q phase, wherein each phase can independently satisfy mass and momentum conservation equations as follows:
the conservation of mass equation for q-phase is:
Figure BDA0001624558380000021
the q-phase conservation of momentum equation is:
Figure BDA0001624558380000022
in the formula, muqAnd λqRespectively, represents the shear and bulk viscosity (Pa · s) of the q-phase, FqExternal volumetric force, Flift,qIs a lifting force, FVm,qTo virtual mass force, RpqIs the interaction force between the phases, p is the pressure shared by all phases;
C. using a k-epsilon turbulence model, the formula of the turbulence energy k and the dissipation ratio epsilon of the turbulence energy is as follows:
Figure BDA0001624558380000023
Figure BDA0001624558380000024
wherein μ is a molecular viscosity (Pa · s) and μtIs turbulent viscosity (Pa · s), PkIs the turbulent shear yield term [ kg/(m.s)3)]，Cε1、Cε2、σk、σεIs constant and is 1.44, 1.92, 1, 1.3 respectively;
D. introducing a three-dimensional geometric model of the atomizing disc into FLUENT of ANSYS Workbench, establishing a reasonable calculation domain and a physical model on the assumption and theoretical basis of step A, B, C, setting various parameters, simulating and calculating the thickness value of the liquid film, and obtaining the relation between the thickness of the liquid film and the parameters by a dimensional analysis and linear least square regression analysis method;
E. and (4) designing an atomization experiment model, and comparing and analyzing the numerical simulation result of the step D and the experiment result, so that the applicability of the numerical simulation method is verified.
Compared with the prior art, the invention has the following obvious and prominent substantive characteristics and remarkable technical progress:
the numerical simulation method for the atomization film-forming process of the suspension material is adopted, numerical simulation is carried out on the atomization film-forming process of the suspension material by adopting computational fluid dynamics software FLUENT, high cost and technical risk caused by experiments or blind design are avoided to a certain extent by adopting the numerical simulation method, and the numerical simulation method has important guiding significance for the atomization mechanism of the feed liquid.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention are specifically described below, and it is obvious that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The first embodiment is as follows:
the numerical simulation method for the atomization film forming of the suspension liquid material comprises the following specific operation steps:
(1) drawing a three-dimensional geometric model diagram of the atomizing disc by using three-dimensional drawing software:
modeling is carried out on the three-dimensional Model of the atomizing disc by utilizing a modeling module (Design Model) in ANSYS Workbench
(2) Establishing a finite element model of the atomizing disc and solving the finite element model, wherein the specific method comprises the following steps:
A. setting a hypothetical condition of numerical simulation;
a) heat transfer between gas and liquid is not considered;
b) the liquid is poured continuously along the central axis of the tray;
c) to reduce the number of calculations, a limited portion of the air near the atomizing disk is included;
B. modeling was performed using the Euler-Euler method, with both gas and liquid phases treated as continuous phases, in order to be able to describe the interpenetrated and continuous multiphase flow between the phasesThe phase volume fraction concept is presented, using alphaqAnd (3) representing the volume fraction occupied by the q phase, wherein each phase can independently satisfy mass and momentum conservation equations as follows:
the conservation of mass equation for q-phase is:
Figure BDA0001624558380000031
the q-phase conservation of momentum equation is:
Figure BDA0001624558380000032
in the formula, muqAnd λqRespectively, represents the shear and bulk viscosity (Pa · s) of the q-phase, FqExternal volumetric force, Flift,qIs a lifting force, FVm,qTo virtual mass force, RpqIs the interaction force between the phases, p is the pressure shared by all phases;
C. using a k-epsilon turbulence model, the formula of the turbulence energy k and the dissipation ratio epsilon of the turbulence energy is as follows:
Figure BDA0001624558380000033
Figure BDA0001624558380000034
wherein μ is a molecular viscosity (Pa · s) and μtIs turbulent viscosity (Pa · s), PkIs the turbulent shear yield term [ kg/(m.s)3)]，Cε1、Cε2、σk、σεIs constant and is 1.44, 1.92, 1, 1.3 respectively;
D. introducing a three-dimensional geometric model of the atomizing disc into FLUENT of ANSYS Workbench, establishing a reasonable calculation domain and a physical model on the assumption and theoretical basis of step A, B, C, setting various parameters, simulating and calculating the thickness value of the liquid film, and obtaining the relation between the thickness of the liquid film and the parameters by a dimensional analysis and linear least square regression analysis method;
E. designing an atomization experiment model, and comparing and analyzing the numerical simulation result of the step D and the experiment result, thereby verifying the applicability of the numerical simulation method
Example two:
this example specifically describes a numerical simulation method for an atomization process of preparing a steel shot from a molten metal, which at least includes the following steps:
firstly, a modeling module (DM) in ANSYS Workbench is utilized to model a three-dimensional model of the atomizing disc:
75mm in diameter, 1500rpm in rotation speed, 2.5kg/min in feed liquid flow, 2590kg/m in density3Viscosity 0.7 pas. The FLUENT software fluid volume function model is selected to capture the interface between the liquid phase and the gas phase. And (3) creating a new material in FLUENT, and inputting related parameters according to the physical properties of the material, namely modeling the molten metal.
Secondly, establishing a finite element model of the atomizing disc and solving, wherein it should be noted that the method is a numerical simulation method based on FLUENT in ANSYS Workbench, the basic method and steps are similar to those of the traditional method and steps, which are not repeated herein, and the hypothesis and theoretical basis of the invention are introduced below.
Firstly, the atomization film-forming process of molten metal is complex, and it is difficult to completely reproduce the atomization film-forming process of molten metal by a numerical simulation method, and an ideal result is not easy to obtain. Therefore, the numerical simulation of the atomization film-forming process of molten metal according to the present invention is based on the following assumptions:
(1) without regard to heat transfer between the gas and liquid.
(2) The liquid continuously pours along the central axis of the tray.
Secondly, modeling gas-liquid two phases by using an Euler-Euler method, wherein the gas-liquid two phases are treated as continuous phases, and in order to describe the multiphase flow which is mutually penetrated and continuous, a phase volume fraction concept is provided and alpha is usedqAnd (3) representing the volume fraction occupied by the q phase, wherein each phase can independently satisfy mass and momentum conservation equations as follows:
the conservation of mass equation for q-phase is:
Figure BDA0001624558380000041
the q-phase conservation of momentum equation is:
Figure BDA0001624558380000042
in the formula, muqAnd λqRespectively, represents the shear and bulk viscosity (Pa · s) of the q-phase, FqExternal volumetric force, Flift,qIs a lifting force, FVm,qTo virtual mass force, RpqIs the interaction force between the phases, p is the pressure shared by all phases;
thirdly, using a k-epsilon turbulence model, the formula of the turbulence energy k and the dissipation rate epsilon of the turbulence energy is as follows:
Figure BDA0001624558380000043
Figure BDA0001624558380000044
wherein μ is a molecular viscosity (Pa · s) and μtIs turbulent viscosity (Pa · s), PkIs the turbulent shear yield term [ kg/(m.s)3)]，Cε1、Cε2、σk、σεIs a constant, 1.44, 1.92, 1, 1.3 respectively.
Introducing a three-dimensional geometric model of an atomizing disc into FLUENT of ANSYS Workbench, establishing a reasonable calculation domain and a physical model on the basis of the assumptions and theories, setting various parameters, simulating and calculating the thickness value of the liquid film, and obtaining the relation between the thickness of the liquid film and the parameters by a dimensional analysis and linear least square regression analysis method; an atomization film forming experiment model is designed, the numerical simulation result and the experiment result are compared and analyzed, and the applicability of the numerical simulation method is verified.
It should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting the same; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (1)

1. A numerical simulation method for atomization film formation of suspension liquid materials is characterized by comprising the following specific operation steps:
(1) drawing a three-dimensional geometric model diagram of the atomizing disc by using three-dimensional drawing software:
modeling is carried out on three-dimensional model of atomizing disc by utilizing modeling module in ANSYSTEMS workbench
(2) Establishing a finite element model of the atomizing disc and solving, wherein the specific method comprises the following steps:
A. setting a hypothetical condition of numerical simulation;
a) heat transfer between gas and liquid is not considered;
b) the liquid is poured continuously along the central axis of the tray;
c) to reduce the number of calculations, a limited portion of the air near the atomizing disk is included;
B. modeling is carried out by using an Euler-Euler method, gas-liquid two phases are treated as continuous phases, and in order to describe the multiphase flow which is mutually penetrated and continuous, a phase volume fraction concept is provided and alpha is usedqAnd (3) representing the volume fraction occupied by the q phase, wherein each phase can independently satisfy mass and momentum conservation equations as follows:
the conservation of mass equation for q-phase is:
Figure FDA0003552623910000011
the q-phase conservation of momentum equation is:
Figure FDA0003552623910000012
in the formula, muqAnd λqRespectively, represents the shear and bulk viscosity (Pa · s) of the q-phase, FqExternal volumetric force, Flift,qIs a lifting force, FVm,qTo virtual mass force, RpqIs the interaction force between the phases, p is the pressure shared by all phases;
C. using a k-epsilon turbulence model, the formula of the turbulence energy k and the dissipation ratio epsilon of the turbulence energy is as follows:
Figure FDA0003552623910000013
Figure FDA0003552623910000014
wherein μ is a molecular viscosity (Pa · s) and μtIs the turbulent viscosity (Pa s), PkIs the turbulent shear yield term [ kg/(m.s)3)]，Cε1、Cε2、σk、σεIs constant and is 1.44, 1.92, 1, 1.3 respectively;
D. introducing a three-dimensional geometric model of the atomizing disc into FLUENT of ANSYS Workbench, establishing a reasonable calculation domain and a physical model on the assumption and theoretical basis of the step A, B, C, setting various parameters, simulating and calculating the thickness value of the liquid film, and obtaining the relation between the thickness of the liquid film and the parameters through dimensional analysis and a linear least square regression analysis method;
E. and (4) designing an atomization experiment model, and comparing and analyzing the numerical simulation result of the step D and the experiment result, so that the applicability of the numerical simulation method is verified.
CN201810318556.3A 2018-04-11 2018-04-11 Numerical simulation method for atomization film formation of suspension material Active CN108804744B (en)

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 CN108804744A (en) 2018-11-13
CN108804744B true 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)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110414141B (en) * 2019-07-30 2022-11-04 辽宁工程技术大学 Three-dimensional numerical simulation method for liquid drop atomization in process of transonic flow of compressible fluid
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
Eulerian multiphase population balance model of atomizing, swirling flows;Narayana P. Rayapati等;《International journal of spray and combustion dynamics》;20111231;第14-44页 *
Mechanism study and numerical simulation of Uranium nitriding induced by high energy laser;Yuan Zhu等;《E3S Web of Conferences》;20180131;第1-4页 *
基于树形自适应网格的旋流液膜雾化过程仿真;杨国华等;《推进技术》;20180117(第03期);第81-89页 *
基于铀的激光氮化机理研究和数值模拟;赵辉等;《材料导报》;20160725(第14期);第157-162页 *
气泡雾化喷嘴雾化射流场性能仿真;毛传林等;《计算机仿真》;20130915(第09期);第219-223页 *
离心喷嘴内部流动与液膜初级破碎的耦合模拟;徐让书等;《沈阳工业大学学报》;20111215(第06期);第63-68、103页 *
空气喷涂平面成膜的双流体模型模拟;陈雁等;《后勤工程学院学报》;20151130(第06期);第86-90页 *

Also Published As

Publication number Publication date
CN108804744A (en) 2018-11-13

Similar Documents

Publication Publication Date Title
Li et al. Investigation on the multiphase sink vortex Ekman pumping effects by CFD-DEM coupling method
Sheikholeslami et al. Nanofluid flow and heat transfer in a rotating system in the presence of a magnetic field
CN108804744B (en) Numerical simulation method for atomization film formation of suspension material
Liu et al. Modeling of quasi-four-phase flow in continuous casting mold using hybrid Eulerian and Lagrangian approach
Zhang et al. Large eddy simulation on flow structure in a dissipative ladle shroud and a tundish
Morales-Higa et al. Ladle shroud as a flow control device for tundish operations
Lopez et al. Recent developments of a numerical model for continuous casting of steel: Model theory, setup and comparison to physical modelling with liquid metal
Li et al. Multiscale mathematical model with discrete–continuum transition for gas–liquid–slag three-phase flow in gas-stirred ladles
Liu et al. Mathematical modeling of multi-sized argon gas bubbles motion and its impact on melt flow in continuous casting mold of steel
Yang et al. Hydrodynamic Modeling and Mathematical Simulation of Flow Field and Temperature Profile for Molten Stainless Steel in an Asymmetrical T‐Type Single‐Strand Continuous Casting Tundish with Arch or Round Hole (s) at Dam Bottom
Liu et al. Numerical Investigation on Motion and Removal of Inclusions in Continuous Casting Tundish with Multiorifice Filter
Hsu et al. Large eddy simulations of turbulent Couette–Poiseuille and Couette flows inside a square duct
Yan et al. Numerical study on the effect of a novel swirling flow generator for submerged entry nozzle in tundish
Mishra et al. Numerical modelling of Sen and mould for continuous slab casting
CN108345737B (en) Design method of bloom continuous casting rotational flow water gap
Tran et al. CFD pre-study of Nozzle reactor for fast hydrothermal liquefaction
Vakhrushev et al. Experimental verification of a 3-phase continuous casting simulation using a water model
Kowitwarangkul et al. CFD Simulation of Molten Steel Flow with Isothermal Condition in the Continuous Casting Tundish
CN109376452A (en) A kind of method for numerical simulation of gathering line hot-washing wax remover phase-change heat transfer
Zhou et al. Gas-liquid two phase flow modelling of incompressible fluid and experimental validation studies in vertical centrifugal casting
Guo et al. Numerical study on the formation and solidification of LMPA microdroplet in a microfluidic device
Mao et al. Uniform flow field design in porous media filter tower and experimental verification
Zhou et al. Study on the Cavity Forming Induced by a Gas Jet Impinging on a Liquid Surface Based on a Deformed Mesh Method
Wang et al. Numerical simulation of effect of various parameters on atomization in an annular slit atomizer
Xia et al. Flow and heat transfer performance of slag and matte in the settler of a copper flash smelting furnace

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