CN113935257A - Gas-liquid-solid three-phase flow simulation method - Google Patents

Gas-liquid-solid three-phase flow simulation method Download PDF

Info

Publication number
CN113935257A
CN113935257A CN202111164848.4A CN202111164848A CN113935257A CN 113935257 A CN113935257 A CN 113935257A CN 202111164848 A CN202111164848 A CN 202111164848A CN 113935257 A CN113935257 A CN 113935257A
Authority
CN
China
Prior art keywords
simulation
particles
phase
liquid
gas
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
Application number
CN202111164848.4A
Other languages
Chinese (zh)
Other versions
CN113935257B (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.)
Yanshan University
Original Assignee
Yanshan University
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 Yanshan University filed Critical Yanshan University
Priority to CN202111164848.4A priority Critical patent/CN113935257B/en
Publication of CN113935257A publication Critical patent/CN113935257A/en
Application granted granted Critical
Publication of CN113935257B publication Critical patent/CN113935257B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a gas-liquid-solid three-phase flow simulation method, which is used for truly simulating the motion track and the deposition position of solid particles in a hydraulic oil tank. The gas-liquid-solid three-phase flow simulation method comprises the following steps: setting simulation parameters of a liquid phase and a gas phase in the FLEUNT to perform gas-liquid two-phase steady state simulation; setting solid phase simulation parameters in the EDEM; and the gas-liquid two-phase simulation data are kept unchanged, the simulation mode is switched to be transient, and the CFD-DEM coupling interface is accessed. And then connecting FLUENT and EDEM through a coupling interface to perform gas-liquid-solid three-phase flow simulation until the simulation is completely stable, and ending the simulation. The invention has the advantages of simple parameter setting, high accuracy, intuitive and understandable simulation result, short simulation period, high efficiency, wide application range, strong universality and the like.

Description

Gas-liquid-solid three-phase flow simulation method
Technical Field
The invention belongs to the technical field of hydraulic simulation, and particularly relates to a gas-liquid-solid three-phase flow simulation method.
Background
With the development of computer technology, Computational Fluid Dynamics (CFD) is widely used for Fluid motion research. Parameters such as speed, temperature and pressure of the flow field can be effectively analyzed through CFD simulation, and the method has important guiding significance on design and optimization of the flow field structure.
When oil in a hydraulic system enters a hydraulic oil tank, solid impurities and bubbles are often mixed to form multiphase flow. To study multiphase flow in a hydraulic tank, it is often necessary to perform simulation calculations.
The current simulation methods for analyzing multiphase flow mainly include an Euler-Euler method and an Euler-Lagrange method. The euler-euler method looks at a certain spatial point in the flow field, considers the research object as a continuous phase, and researches the motion condition of fluid particles passing through the spatial point. The euler-lagrange method looks at the motion of fluid particles, treats the object of study as a discrete phase, tracks each particle, and observes and analyzes the motion history of each discrete particle. Among them, the euler-lagrange Method includes a Discrete Particle Method (DPM) and a Discrete Element Method (DEM).
However, most of the existing multiphase flow simulation methods only perform gas-liquid two-phase flow or solid-liquid two-phase flow calculation, and gas-liquid-solid three-phase flow is not considered at the same time, so that the actual conditions in the hydraulic oil tank are different, and the simulation result only has a certain reference value.
In addition, solid particles mixed with oil in the oil tank often contain metal particles and nonmetal particles, the particles are different in size and irregular in shape, and the stress condition of the solid particles in the oil is very complex, and the solid particles can be subjected to the buoyancy, drag force, pressure gradient force and other forces of the oil to the particles, the contact force between the particles and the wall surface and the like. The existing three-phase flow simulation method is mainly used for simulating by utilizing a two-phase flow model and a Discrete Particle Model (DPM) in CFD simulation software FLUENT, the simulation method cannot completely consider the shape and the size and the stress condition of particles in practice, the particles are used as regular spheres for simulation, and although the shape coefficient can be set according to a mathematical method, the shape and the size of the particles cannot be close to the shape and the size of the actual particles. And the DPM model only considers the force of oil liquid to solid particles, neglects the contact force between particles and wall surfaces, and is not in line with the practice. In addition, from simulation results, the method can only roughly simulate the movement track of the particles in the hydraulic oil tank, and cannot accurately display and predict the final deposition positions of the particles.
Therefore, the existing multiphase flow simulation method cannot truly simulate the real movement condition of particles in the hydraulic oil tank. How to consider the influence between gas phase and particles, particle to particle and particle to constraint in multiphase flow simulation and to accurately simulate the motion trail and deposition position of particles as much as possible is a problem to be solved urgently in multiphase flow simulation.
Disclosure of Invention
Aiming at the technical problems in the existing simulation method, the invention provides a gas-liquid-solid three-phase flow simulation method to truly simulate the movement condition of particles in a hydraulic oil tank. The simulation method is based on a CFD-DEM coupling method, fully considers the interaction between gas phase and liquid phase and solid phase, the actual shape and size of particles and the stress condition of the particles in a flow field, and can accurately and effectively predict the motion track and the deposition position of solid particles in a hydraulic oil tank.
Specifically, the invention provides a gas-liquid-solid three-phase flow simulation method, which comprises the following steps:
step S1, setting simulation parameters of liquid phase and gas phase in simulation software FLEUNT, and performing gas-liquid two-phase steady state simulation;
the gas-liquid two-phase steady-state simulation adopts an Euler-Euler method, a liquid phase and a gas phase are used as continuous phases for simulation, and simulation parameters of the liquid phase and the gas phase are set in simulation software FLUENT;
then setting the liquid phase simulation and the gas phase simulation as a steady state simulation, and performing step S2 after the liquid phase simulation result and the gas phase simulation result enter the steady state;
s2, setting solid phase simulation parameters in simulation software EDEM, keeping gas-liquid two-phase simulation data unchanged, switching a simulation mode into transient simulation, and then accessing a CFD-DEM coupling interface;
the CFD-DEM coupling interface in the step S2 includes two coupling interfaces, the two coupling interfaces are respectively a coupling interface based on multiphase flow and a coupling interface based on DPM, the two coupling interfaces respectively include two calculation and setting methods according to the density of particles in the liquid, and the two calculation and setting methods are respectively a calculation and setting method of volume fraction of particles to be considered and a calculation and setting method of volume fraction of particles not to be considered;
the particle volume fraction is defined as follows:
Figure BDA0003291354170000031
in the formula, alphapAnd alphalVolume fractions of particles and liquid are respectively, and when the volume fraction eta of the particles is more than 10%, the volume fraction of the particles needs to be considered in a corresponding simulation method; when the volume fraction η of the particles is less than or equal to 10%, the corresponding simulation method does not need to consider the volume fraction of the particles;
step S3, connecting FLUENT and EDEM through the selected coupling interface, performing gas-liquid-solid three-phase flow simulation, and ending the simulation until the simulation is completely stable;
in the step S3, the euler-lagrange method is adopted for gas-liquid-solid three-phase simulation, the liquid phase and the gas phase are regarded as continuous phases for simulation, the solid particles are regarded as discrete phases for simulation, solid-phase simulation parameters are set in the EDEM, the solid particles include forces applied to the solid particles during simulation, and the stress model of the solid particles is as follows:
Figure BDA0003291354170000032
Figure BDA0003291354170000033
in the formula, mpAnd IpMass and inertia tensor for the particles, respectively; u. ofpAnd ωpRespectively the linear and angular velocity of the particles;Ffis the force of the fluid on the particles; fcIs the contact force to which the particles are subjected; t iscIs the contact torque to which the particles are subjected.
Preferably, if the volume fraction of the particles is not considered in the coupling interface based on the multiphase flow, an Euler model needs to be started, after the model is started, gas-liquid two-phase flow simulation is performed in FLUENT, and the particles and the fluid interact through a self-defined source term;
if the multiphase flow-based coupling interface considers the volume fraction of particles, after an Euler model is started, the Euler phase number is set to be 3 phases, gas-liquid-solid three-phase flow simulation is performed in FLUENT, and the particles and the fluid interact through a custom source term.
Preferably, if the volume fraction of the particles is not considered, the coupling interface based on the DPM starts an euler model and a DPM model in FLUENT, sets gas-liquid two-phase simulation in the Eulerian model, sets solid-phase simulation in the DPM model, and initializes DPM information in the current step by parameters such as the position, the volume, the speed and the like of the particles in the EDEM;
if the volume fraction of the particles is considered, the Eulerian model and the DPM model are started in FLUENT, the DDPM model is activated, and the rest parameter settings are consistent with the condition that the volume fraction of the particles is not considered.
The forces applied to the particles by the movement of the particles in the oil tank can be divided into two types, namely the acting force applied to the particles by the flow field, the acting force between the particles and the wall surface. The existing two-phase flow simulation method only considers the force applied to particles by a partial flow field, and does not consider the acting force between the particles and the wall surface. The invention fully considers the actual motion situation of the particles in the oil tank, improves the existing simulation method, and adds the action of gas phase on the particles to form three-phase flow simulation; the force applied to the particles is compensated according to the actual movement, which is more practical. In the existing CFD-DEM coupling interface, only the floating force, drag force and gravity force of the particles in the flow field are considered, but in addition to the forces, the particles moving in the actual flow field are also subjected to Saffman force, baseset force, virtual mass force, pressure gradient force, Magnus force and the like. The invention rewrites the existing coupling interface program and adds the residual force of the flow field on the particles, so that the particles can better meet the actual situation of particle motion.
Preferably, wherein FcThe contact force to which the particles are subjected is expressed as:
Fc＝Fc,n+Fc,t
in the formula, Fc,nAnd Fc,tRespectively representing normal force and tangential force applied in the process of contacting the particles and the wall surface;
Tcthe contact torque to which the particles are subjected is expressed as:
Tc＝Tt+Tr
in the formula, TtAnd TrRepresenting contact torques generated by tangential contact force and rolling friction, respectively;
Ffthe acting force of the fluid on the particles is expressed as
Ff＝FG+FB+FP+FDrag+FVR+FSaff
Figure BDA0003291354170000051
Figure BDA0003291354170000052
Figure BDA0003291354170000053
In the formula ufIs the velocity of the fluid; fGThe particles are subjected to gravity; fBIs buoyancy; fpIs a pressure gradient force; fDragIs the fluid drag force; fVRIs a virtual mass force; fSaffIs Saffman lifting force; cDIs the drag coefficient; dpIs the particle diameter; rhofAnd ρpAre respectively fluid andthe density of the particles; μ is the hydrodynamic viscosity; reIs the particle Reynolds number.
The invention has the following beneficial effects:
(1) the gas-liquid-solid three-phase flow simulation method has the advantages of simple parameter setting, high accuracy, intuitive and understandable simulation result, short simulation period, high efficiency, wide application range, strong universality and the like.
(2) All parameters of the gas-liquid-solid three-phase flow simulation method are set in simulation software FLUENT and EDEM software, and the operation is simple; the gas-liquid-solid three-phase flow simulation method can be used for modeling according to the actual shape and size of particles, setting parameters such as density, Poisson's ratio and elastic modulus of the particles and really restoring the intrinsic properties of the particles; contact parameters such as friction coefficient, coefficient of restitution and the like between the particles and the contact surface can be set according to actual working conditions; the contact action between particles and between the particles and the wall surface is considered, so that the method is more practical and has high simulation accuracy;
(3) the gas-liquid-solid three-phase flow simulation method rewrites the existing coupling interface program, supplements the force applied by the particles in the flow field, better accords with the actual situation and has high simulation accuracy; the result obtained by the gas-liquid-solid three-phase flow simulation method directly shows the specific deposition position of the particles in the simulation model, and intuitively shows the motion condition of the particles;
(4) the multiphase flow simulation method has the advantages of short simulation period, less consumed computing resources and high computing efficiency; the gas-liquid-solid three-phase flow simulation method is suitable for not only an oil tank, but also numerical simulation of any flowing liquid doped with solid phase and gas phase in theory, and has wide application range and strong universality.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate exemplary embodiments of the invention and, together with the description, serve to explain the invention and are not intended to limit the invention. In the drawings:
FIG. 1 is a flow chart of simulation disclosed in the embodiments of the present invention;
FIG. 2 is a diagram of a simulation application model and process disclosed in an embodiment of the present invention;
FIG. 3 is a flow chart of CFD-DEM coupled simulation disclosed in the embodiments of the present invention;
FIG. 4a and FIG. 4b are both simulation coupling model improvement algorithms disclosed in the embodiments of the present invention;
fig. 5a to 5c are graphs of simulation results of three-phase flow simulation of a hydraulic oil tank according to an embodiment of the present invention, respectively, where fig. 5a is a graph of oil velocity simulation, fig. 5b is a graph of gas volume fraction simulation, and fig. 5c is a graph of particle deposition position simulation;
fig. 6a and fig. 6b are comparison diagrams of particle deposition positions of three-phase flow and two-phase flow of a hydraulic oil tank according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail below with reference to the accompanying drawings.
DETAILED DESCRIPTION OF EMBODIMENT (S) OF INVENTION
In order to make the technical solutions and advantages of the present invention more apparent, exemplary embodiments of the present invention are described in further detail below with reference to the accompanying drawings. It is clear that the described embodiments are only a part of the embodiments of the invention, and not an exhaustive list of all embodiments. And the embodiments and features of the embodiments may be combined with each other without conflict.
FIGS. 1 and 2 schematically show a gas-liquid-solid three-phase flow simulation method according to an embodiment of the invention, comprising the steps of:
step S1: the liquid and gas phases were subjected to steady state simulation in the simulation software FLUENT. And the gas-liquid simulation adopts an Euler-Euler model, the gas phase and the liquid phase are both regarded as continuous phases, and relevant simulation parameters of the gas phase and the liquid phase are set in FLUENT to perform steady-state simulation.
Step S2: and setting solid phase simulation parameters in the EDEM, keeping the gas-liquid two-phase simulation data unchanged, switching the simulation mode to be transient, and accessing the CFD-DEM coupling interface.
The CFD-DEM coupling interface in the step S2 includes two coupling interfaces, the two coupling interfaces are respectively a coupling interface based on multiphase flow and a coupling interface based on DPM, the two coupling interfaces respectively include two calculation and setting methods according to the density of particles in the liquid, and the two calculation and setting methods are respectively a calculation and setting method of volume fraction of particles to be considered and a calculation and setting method of volume fraction of particles not to be considered;
the particle volume fraction is defined as follows:
Figure BDA0003291354170000071
in the formula, alphapAnd alphalVolume fractions of particles and liquid are respectively, and when the volume fraction eta of the particles is more than 10%, the volume fraction of the particles needs to be considered in a corresponding simulation method; when the particle volume fraction η is less than or equal to 10%, the corresponding simulation method does not need to consider the volume fraction of the particles.
Step S3: and the CFD-DEM coupling interface is connected with FLUENT and EDEM to perform gas-liquid-solid three-phase flow simulation. The gas-liquid-solid simulation part is transient simulation, calculation is carried out by adopting an Euler-Lagrange method, a gas phase and a liquid phase are taken as continuous phases, a solid phase is taken as a discrete phase, EDEM is connected in FLUENT, and simulation is started. And after the simulation is finished, storing the simulation data and checking the simulation result.
In a preferred embodiment, as shown in FIG. 3. FLUENT and EDEM coupling process is as follows: firstly, performing bubble flow simulation in FLUENT, performing steady state simulation on gas-liquid phases by adopting an Euler-Euler simulation method, performing transient flow field calculation, and accessing an EDEM coupling interface after the steady state is waited; and setting particle attributes and defining contact models and parameters in the EDEM, and setting the particle number and time for throwing the particles. After FLUENT and EDEM coupling calculation is started, the coupling interface can transfer flow field force borne by the particles to the EDEM, the EDEM calculates resultant force borne by the particles according to a Newton second law, the position and the speed of the particles are updated, the information is converted into momentum action of the particles on the flow field, and the momentum action is transferred to the FLUENT through the coupling interface. And repeating the calculation until FLUENT reaches the set time step, finishing the coupling calculation and finishing the simulation.
Fig. 4a and 4b schematically show a simulated coupling algorithm modification according to an embodiment of the invention. Since the programming architecture and content of the different forces to which the particles are subjected in the flow field are similar, fig. 4a and 4b only schematically show the programming content of the Saffman force to which the particles are subjected in the coupling interface. The Saffman force is related to the rotation angular velocity of the particles, the algorithm firstly defines the components of the rotation angular velocity, defines other variables used in the Saffman force calculation formula, and finally writes the force calculation formula.
In a preferred embodiment, a simulation result of gas-liquid-solid three-phase flow of a hydraulic oil tank is shown in FIGS. 5a-5 c. The oil liquid mixed with gas and solid particles uniformly enters from an oil return opening of the oil tank, passes through a partition plate in the middle of the oil tank and flows out from an oil suction opening of the oil tank. The oil speed distribution cloud chart shows that the maximum oil speed is distributed on one side of an oil return pipe of an oil tank, namely the left side of the partition plate; the distribution cloud chart of the gas volume fraction shows that the position with the maximum gas content in the oil is positioned on one side of an oil return pipe of an oil tank, and the oil on one side of an oil suction pipe of the oil tank hardly contains gas; the deposition position diagram of the particles in the hydraulic oil tank shows that the particles are almost deposited on one side of an oil return opening of the oil tank, the two corners on the left side of an oil return pipe are provided with gathered particles, and the particles on the right side of the oil return pipe are distributed in an arc shape in simulation.
In a preferred embodiment, a comparison of particle deposition positions of a gas-liquid-solid three-phase flow and a gas-liquid two-phase flow of a hydraulic oil tank is shown in fig. 6a and 6 b. The simulation results of the gas-liquid-solid three-phase flow and the gas-liquid two-phase flow are similar in overall trend, namely, the two corners on the left side of the oil return pipe are provided with gathered particles, and the particles on the right side of the oil return pipe are distributed in an arc shape in the simulation.
The difference between the two is that the overall position of the particles in the gas-liquid-solid three-phase flow simulation result is farther than that of the gas-liquid two-phase flow, namely closer to the wall surface of the oil tank, so that the difference degree between the gas-liquid-solid three-phase flow simulation result and the gas-liquid two-phase flow simulation result is larger, and the importance of considering the action of gas on solid particles when performing fluid simulation on the hydraulic oil tank is reflected.
The above-mentioned embodiments are merely illustrative of the preferred embodiments of the present invention, and do not limit the scope of the present invention, and various modifications and improvements made to the technical solution of the present invention by those skilled in the art without departing from the spirit of the present invention shall fall within the protection scope defined by the claims of the present invention.

Claims (4)

1. A gas-liquid-solid three-phase flow simulation method is characterized by comprising the following steps: which comprises the following steps:
step S1, setting simulation parameters of liquid phase and gas phase in simulation software FLEUNT, and carrying out gas-liquid two-phase steady state simulation;
the gas-liquid two-phase steady-state simulation adopts an Euler-Euler method, a liquid phase and a gas phase are used as continuous phases for simulation, and simulation parameters of the liquid phase and the gas phase are set in simulation software FLUENT;
then setting the liquid phase simulation and the gas phase simulation as a steady state simulation, and performing step S2 after the liquid phase simulation result and the gas phase simulation result enter the steady state;
s2, setting solid phase simulation parameters in simulation software EDEM, keeping gas-liquid two-phase simulation data unchanged, switching a simulation mode into transient simulation, and then accessing a CFD-DEM coupling interface;
the CFD-DEM coupling interface in the step S2 includes two coupling interfaces, the two coupling interfaces are respectively a coupling interface based on multiphase flow and a coupling interface based on a discrete particle method, the two coupling interfaces respectively include two calculation and setting methods according to the density of particles in the liquid, and the two calculation and setting methods are respectively a calculation and setting method of volume fraction of particles to be considered and a calculation and setting method of volume fraction of particles not to be considered;
wherein the volume fraction of the particles is defined by the formula:
Figure FDA0003291354160000011
in the formula, alphapAnd alphalVolume fractions of particles and liquid, respectively;
when the volume fraction eta of the particles is more than 10 percent, the volume fraction of the particles needs to be considered in a corresponding simulation method; when the volume fraction η of the particles is less than or equal to 10%, the corresponding simulation method does not need to consider the volume fraction of the particles;
step S3, connecting FLUENT and EDEM through the selected coupling interface respectively, and performing gas-liquid-solid three-phase flow simulation until the simulation is completely stable and the simulation is finished;
in the step S3, an euler-lagrange method is adopted for gas-liquid-solid three-phase simulation, a liquid phase and a gas phase are taken as continuous phases for simulation, solid particles are taken as discrete phases for simulation, solid phase simulation parameters are set in the EDEM, the solid particles include forces applied to the solid particles during simulation, and the stress model of the solid particles is as follows:
Figure FDA0003291354160000021
Figure FDA0003291354160000022
in the formula, mpAnd IpMass and inertia tensor for the particles, respectively; u. ofpAnd ωpLinear and angular velocities of the particles, respectively; ffIs the force of the fluid on the particles; fcIs the contact force to which the particles are subjected; t iscIs the contact torque to which the particles are subjected.
2. The gas-liquid-solid three-phase flow simulation method according to claim 1, characterized in that: if the multiphase flow-based coupling interface does not consider the volume fraction of the particles, the calculation and setting method specifically comprises the following steps: an Euler model is required to be started, after the Euler model is started, gas-liquid two-phase flow simulation is carried out in FLUENT, and solid particles and fluid interact through a user-defined source term;
if the multiphase flow-based coupling interface considers the volume fraction of the particles, the calculation and setting method specifically comprises the following steps: after the Euler model is started, the Euler phase number is set to be 3 phases, gas-liquid-solid three-phase flow simulation is performed in FLUENT, and solid particles and fluid interact through a custom source item.
3. The gas-liquid-solid three-phase flow simulation method according to claim 1, characterized in that: if the volume fraction of the particles is not considered by the coupling interface based on the discrete particle method, the calculating and setting method specifically comprises the following steps: starting an Eulerian model and a DPM model in FLUENT, setting gas-liquid two-phase simulation in the Eulerian model, setting solid-phase simulation in the DPM model, and initializing DPM information in the current step by the position, volume and speed parameters of particles in EDEM;
if the volume fraction of the particles is considered in the coupling interface based on the discrete particle method, the calculating and setting method specifically comprises the following steps: after the Eulerian model and DPM model were started in FLUENT, the DDPM model was activated and the remaining parameter settings were consistent without regard to the volume fraction of the particles.
4. The gas-liquid-solid three-phase flow simulation method according to claim 1, characterized in that:
wherein, FcThe contact force to which the particles are subjected is expressed as:
Fc＝Fc,n+Fc,t
in the formula, Fc,nAnd Fc,tRespectively representing the normal force and the tangential force applied in the process of contacting the particles and the wall surface;
Tcthe contact torque to which the particles are subjected is expressed as:
Tc＝Tt+Tr
in the formula, TtAnd TrRepresenting contact torques generated by tangential contact force and rolling friction, respectively;
Ffthe acting force of the fluid on the particles is expressed as
Ff＝FG+FB+FP+FDrag+FVR+FSaff
Figure FDA0003291354160000031
Figure FDA0003291354160000032
Figure FDA0003291354160000033
In the formula ufIs the velocity of the fluid; fGThe particles are subjected to gravity; fBIs buoyancy; fpIs a pressure gradient force; fDragIs the fluid drag force; fVRIs a virtual mass force; fSaffIs Saffman lifting force; cDIs the drag coefficient; dpIs the particle diameter; rhofAnd ρpDensity of the fluid and particles, respectively; μ is the hydrodynamic viscosity; reIs the particle reynolds number.
CN202111164848.4A 2021-09-30 2021-09-30 Gas-liquid-solid three-phase flow simulation method Active CN113935257B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111164848.4A CN113935257B (en) 2021-09-30 2021-09-30 Gas-liquid-solid three-phase flow simulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111164848.4A CN113935257B (en) 2021-09-30 2021-09-30 Gas-liquid-solid three-phase flow simulation method

Publications (2)

Publication Number Publication Date
CN113935257A true CN113935257A (en) 2022-01-14
CN113935257B CN113935257B (en) 2023-05-05

Family

ID=79277866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111164848.4A Active CN113935257B (en) 2021-09-30 2021-09-30 Gas-liquid-solid three-phase flow simulation method

Country Status (1)

Country Link
CN (1) CN113935257B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116108631A (en) * 2022-12-13 2023-05-12 徐州徐工挖掘机械有限公司 Simulation method and calibration device for transient degassing process of hydraulic oil tank

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063248A (en) * 2018-06-26 2018-12-21 广东工业大学 It is a kind of for laser impact intensified gas-liquid-solid coupling calculation
CN109840375A (en) * 2019-01-25 2019-06-04 华北电力大学 A kind of confirmation method of liquid-solid fluid bed CFD drag force model
WO2020042821A1 (en) * 2018-08-29 2020-03-05 华南理工大学 Method for calculating pressure loss of parallel r-type automobile shock absorber
CN111144003A (en) * 2019-12-27 2020-05-12 安徽信息工程学院 Numerical simulation test method for double-droplet impact on flat solid surface based on VOF and dynamic grid self-adaption
CN111680366A (en) * 2020-06-16 2020-09-18 郑州大学 Amesim-based aircraft fuel tank mass and mass center simulation calculation method
CN112380788A (en) * 2020-11-06 2021-02-19 天津大学 Semi-analytic calculation method for bidirectional coupling of super-ellipsoid particles and flow field

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063248A (en) * 2018-06-26 2018-12-21 广东工业大学 It is a kind of for laser impact intensified gas-liquid-solid coupling calculation
WO2020042821A1 (en) * 2018-08-29 2020-03-05 华南理工大学 Method for calculating pressure loss of parallel r-type automobile shock absorber
CN109840375A (en) * 2019-01-25 2019-06-04 华北电力大学 A kind of confirmation method of liquid-solid fluid bed CFD drag force model
CN111144003A (en) * 2019-12-27 2020-05-12 安徽信息工程学院 Numerical simulation test method for double-droplet impact on flat solid surface based on VOF and dynamic grid self-adaption
CN111680366A (en) * 2020-06-16 2020-09-18 郑州大学 Amesim-based aircraft fuel tank mass and mass center simulation calculation method
CN112380788A (en) * 2020-11-06 2021-02-19 天津大学 Semi-analytic calculation method for bidirectional coupling of super-ellipsoid particles and flow field

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
纪辉;聂松林;胡锦国;: "基于超重力净油装置的三相分离数值仿真研究" *
颜世东;刘砚山;冯艳;: "基于FLUENT的固体发动机内流场数值仿真" *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116108631A (en) * 2022-12-13 2023-05-12 徐州徐工挖掘机械有限公司 Simulation method and calibration device for transient degassing process of hydraulic oil tank
CN116108631B (en) * 2022-12-13 2024-03-01 徐州徐工挖掘机械有限公司 Simulation method and calibration device for transient degassing process of hydraulic oil tank

Also Published As

Publication number Publication date
CN113935257B (en) 2023-05-05

Similar Documents

Publication Publication Date Title
Wan et al. Direct numerical simulation of particulate flow via multigrid FEM techniques and the fictitious boundary method
Wandel et al. Spline-pinn: Approaching pdes without data using fast, physics-informed hermite-spline cnns
CN112069689B (en) Simulation method and system for fuel atomization characteristic of aircraft engine
EP2008250A1 (en) Method of simulating detailed movements of fluids using derivative particles
CN107341315A (en) The method being predicted using mixed-scale multiphase flow model to multiphase flow
Chawner et al. Progress in geometry modeling and mesh generation toward the CFD vision 2030
CN113935257A (en) Gas-liquid-solid three-phase flow simulation method
Bodnár et al. Numerical simulation of turbulent free-surface flow in curved channel
Baxter et al. Haptic Interaction with Fluid Media.
Onishi et al. Vehicle aerodynamics simulation for the next generation on the K computer: part 2 use of dirty CAD data with modified Cartesian grid approach
Weaver et al. Fluid Simulation by the Smoothed Particle Hydrodynamics Method: A Survey.
Li et al. An efficient Cartesian mesh generation strategy for complex geometries
Watanabe et al. The CFD application for efficient designing in the automotive engineering
JP2010108183A (en) Simulation method and simulation device
Baek et al. Muddy water animation with different details
Zheng The new direction of computational fluid dynamics and its application in industry
Gil et al. Parallel multigrid detached eddy simulation algorithm for three-dimensional unsteady incompressible flows on unstructured meshes
CN103632016A (en) Simulation method based on dynamic mesh technology for hydrostatic bearing pressure field
Ferziger et al. Introduction to numerical methods
Wang et al. Fundamentals of Computational Fluid Dynamics
Yang et al. Motion prediction using wall-resolved and wall-modeled approaches on a Cartesian grid
Jacobs Modelling of multiphase flows on adaptive unstructured meshes with applications to the dynamics of volcanic ash plumes
Yimsiriwatana Airflow simulation of particle suction in hard disk drive manufacturing process
Jackson Truncation Error Based Mesh Adaptation and Its Application to Multi-Mesh CFD
Antepara Zambrano Adaptive mesh refinement method for CFD applications

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