Summary of the invention:
The object of the invention is to the defect overcoming prior art existence, seek to design the method for numerical simulation that a kind of fully mechanized coal face dust movement track is provided, the movement locus that analog computation goes out different-grain diameter dust can be specified, thus take measurements for reduction targetedly.
To achieve these goals, the method for numerical simulation that the present invention relates to specifically comprises following processing step:
(1) according to fluid mechanics viscous sublayer district's flowing law and Buddhist nun La Zi experiment experience formula, in conjunction with the central point of fully mechanized coal face, derive and show that fully mechanized coal face inlet section viscous sublayer district wind speed computing formula is:
Wherein, v
x---fully mechanized coal face inlet section viscous sublayer district wind speed, v
x>=0, unit m/s;
---fully mechanized coal face inlet section mean wind speed,
unit m/s; R---arbitrfary point, fully mechanized coal face inlet section viscous sublayer district (A) to the distance of fully mechanized coal face central point, unit m; r
0---fully mechanized coal face central point to workplace border vertical range and through arbitrfary point, fully mechanized coal face inlet section viscous sublayer district (A), r
0>=r, unit m.
(2) distinguished and admirable according to fully mechanized coal face inlet section average speed value calculates Reynolds number, to judge distinguished and admirable fluidised form:
If 1. Reynolds number is less than 2000, be Laminar Flow, so according to the wind speed formula (7) in fully mechanized coal face inlet section viscous sublayer district, adopt conventional computer C language to write fully mechanized coal face inlet section wind speed self-defining function;
If 2. Reynolds number is greater than 2000, it is Turbulence Flow, the wind speed formula (7) in so comprehensive fully mechanized coal face inlet section viscous sublayer district and non-circular workplace turbulent region wind speed profile formula (8), adopt conventional computer C language to write fully mechanized coal face inlet section wind speed self-defining function;
Wherein, v
x---fully mechanized coal face turbulent region wind speed, v
x> 0, unit m/s;
---fully mechanized coal face inlet section mean wind speed,
unit m/s; R---arbitrfary point, fully mechanized coal face turbulent region (A) to the distance of fully mechanized coal face central point, unit m; r
0---fully mechanized coal face central point to workplace border vertical range and through arbitrfary point, fully mechanized coal face turbulent region (A), r
0> r, unit m;
(3) fully mechanized coal face inlet section wind speed self-defining function is made an explanation in fluid analysis software FLUENT, and be arranged to the boundary condition of fully mechanized coal face inlet section wind speed; Then grain size analysis software EDEM and fluid analysis software FLUENT is intercoupled, and derive dust granules coordinate data to obtain the fitting function of different-grain diameter dust movement track, thus simulation draws the movement locus of fully mechanized coal face different-grain diameter dust.
The present invention compared with prior art, EDEM software is coupled with FLUENT software can simulate the Multiphase Flow of particle scale, this coupling enhances the ability of FLUENT software simulation polyphasic flow, can comprise the feature as particles such as distribution of particles, particle habit, geometric attribute and surface forces, above-mentioned polyphasic flow attribute all has obvious impact to particle and fluid; The method realizes easy, and principle is reliable, design science, and data are accurate, with low cost.
Embodiment:
By embodiment, further description is made to the present invention below in conjunction with accompanying drawing.
Embodiment:
The method for numerical simulation of the fully mechanized coal face dust movement track that the present embodiment relates to comprises the following steps:
(1) derivation fully mechanized coal face inlet section viscous sublayer district wind speed computing formula:
According to fluid mechanics viscous sublayer district's flowing law and Buddhist nun La Zi experiment experience formula, in conjunction with the central point of fully mechanized coal face, derive and draw fully mechanized coal face inlet section viscous sublayer district wind speed computing formula.
Lineal layout is thought by approximate for the velocity distribution of wind speed, that is: in fully mechanized coal face inlet section viscous sublayer district (being called for short viscous sublayer district)
In formula,
---velocity gradient;
V
x---viscous sublayer district random spot speed, unit m/s, v
x>=0;
Y---arbitrfary point, viscous sublayer district to the vertical range on workplace border, unit m, y >=0;
By stress approximate representation be:
In formula, δ---the thickness in viscous sublayer district, unit m;
τ---viscous shearing stress, unit N/m
2;
μ---viscosity, unit Pas.
Can be obtained by formula (2):
In formula, v
*---resistance speed, unit m/s;
ν---kinematic viscosity, unit m
2/ s.
Order
In formula, the density of ρ---air, unit K g/m
3.
Formula (4) is derived according to Bernoulli equation and Darcy's law and is drawn following formula:
In formula, λ---Darcy coefficient, zero dimension;
α---coefficient of frictional resistance, units/kg/m
3;
Learn according to Ni Gulazi experiment and meet following formula at laminar region λ:
In formula, d---equivalent diameter, unit m;
R
e---Reynolds number, zero dimension.
Formula (5), (6) are substituted into formula (3), uses after arranging
r/r
0the wind speed formula in the fully mechanized coal face inlet section viscous sublayer district represented is as follows:
In formula, v
x---fully mechanized coal face inlet section viscous sublayer district wind speed, v
x>=0, unit m/s;
---fully mechanized coal face inlet section mean wind speed,
unit m/s; R---arbitrfary point, viscous sublayer district (A) to the distance of fully mechanized coal face central point, unit m; r
0---fully mechanized coal face central point to workplace border vertical range and through arbitrfary point, viscous sublayer district (A), r
0>=r, unit m.
(2) write fully mechanized coal face inlet section wind speed self-defining function: for fully mechanized coal face in reality, calculate Reynolds number with the distinguished and admirable fluidised form of interpretation by the average velocity that fully mechanized coal face inlet section is distinguished and admirable:
If 1. Reynolds number is less than 2000, be Laminar Flow, according to the wind speed formula (7) in fully mechanized coal face inlet section viscous sublayer district, adopt C-language Programming Design to write fully mechanized coal face inlet section wind speed self-defining function;
If 2. Reynolds number is greater than 2000, it is Turbulence Flow, the wind speed formula (7) in comprehensive fully mechanized coal face inlet section viscous sublayer district and non-circular workplace turbulent region wind speed profile formula (8), adopt C-language Programming Design to write fully mechanized coal face inlet section wind speed self-defining function;
In formula, v
x---fully mechanized coal face turbulent region wind speed, v
x> 0, unit m/s;
---fully mechanized coal face inlet section mean wind speed,
unit m/s;
R---arbitrfary point, fully mechanized coal face turbulent region (A) to the distance of fully mechanized coal face central point, unit m;
R
0---fully mechanized coal face central point to workplace border vertical range and through arbitrfary point, fully mechanized coal face turbulent region (A), r
0> r, unit m.
Described writing of fully mechanized coal face inlet section wind speed self-defining function with operating procedure is:
1. fully mechanized coal face inlet section wind speed self-defining function write step
User-Defined Functions header file;
Fully mechanized coal face inlet section wind speed self-defining function DEFINE_PROFILE;
Variable (comprising array variable x [ND_ND], coordinate variable a, b, face data variable face_t f) is set;
Loop statement begin_f_loop (f, thread), to be used for scanning in whole workplace region all;
Function F _ CENTROID (x, f, thread), obtains the barycenter in each face, and the coordinate figure assignment of barycenter to x;
Assignment a=x [0] is carried out to coordinate variable, b=x [1];
Use if, else statement, according to the different situations of coordinate, fully mechanized coal face inlet section air current state is different, gives corresponding computing formula to F_PROFILE (f, thread, position);
End loop statement f_loop (f, thread)
2. the operating procedure of fully mechanized coal face inlet section wind speed self-defining function
In FLUENT software, click menu define-user defined – function – interpreted;
In the dialog box interpreted UDFs ejected, click Browse ..., the fully mechanized coal face User-Defined Functions of writing is imported;
Button click interprete.
Certain colliery 3 is chosen as with the geometric model of fully mechanized coal face in the present embodiment
on1110 fully mechanized coal faces are that example is described, according to 3
onthe on-the-spot actual measurement of 1110 fully mechanized coal faces, the pre-processing software GAMBIT adopting FLUENT software to carry sets up physical model, is respectively by the length of model and is set as 20.0m, 10.0m, 3.0m;
First, in the present embodiment 3
onthe distinguished and admirable fluidised form of 1110 fully mechanized coal faces differentiates, when fully mechanized coal face inlet section mean wind speed is 2m/s, calculates Reynolds number:
In formula, the area of S---fully mechanized coal face, unit m
2;
The girth of C---fully mechanized coal face, unit m.
Be greater than 2000 by formula (9) known Reynolds number, then the flow state of fully mechanized coal face is Turbulence Flow.
Secondly, according to wind speed formula (7), non-circular workplace turbulent region wind speed profile formula (8) and 3 in fully mechanized coal face inlet section viscous sublayer district
onthe actual conditions of 1110 fully mechanized coal faces, adopt C language to write fully mechanized coal face inlet section wind speed self-defining function udfvelocity_inlet, specific code is as follows:
(3) numerical simulation movement locus: above-mentioned fully mechanized coal face inlet section wind speed self-defining function udf velocity_inlet is made an explanation in fluid analysis software FLUENT, and is arranged to the boundary condition of fully mechanized coal face inlet section wind speed; Then grain size analysis software EDEM and fluid analysis software FLUENT is intercoupled, and derive dust granules coordinate data to obtain the fitting function of different-grain diameter dust movement track, thus simulation draws the movement locus of fully mechanized coal face different-grain diameter dust, concrete numerical simulation is as follows:
1. grain size analysis software EDEM is coupled after installing with fluid analysis software FLUENT, fully mechanized coal face grid model file is imported in fluid analysis software FLUENT, then computation model, discrete type parameter, boundary condition and sources of particles major parameter are set, and solve parameter;
2. in grain size analysis software EDEM, fully mechanized coal face grid model file is imported, then set environment parameter, solid, the parameter of the dust such as particle diameter, material, quantity, calculating parameter;
3. in fluid analysis software FLUENT, open the EDEM-uncoupled module under Models, the coupling process of EDEM and FLUENT is set wherein, and EDEM file is called in;
4. particle diameter is set in grain size analysis software EDEM and (comprises 1 × 10
3um, 1.6 × 10
2um, 1.0 × 10
2um, 50um, 8um), each particles number is respectively 100, when the computing time that is coupled is 0.3579s, the movement locus of each particle diameter dust as Figure 4-8, less by the known dust size of accompanying drawing content, the identical time is larger to the scope of outlet diffusion, less with z-axis negative direction angle.
The fluid analysis software FLUENT related in the present embodiment is CFD software package general at present, can be used for simulating the Complex Flows be never compressible in highly compressible scope, it includes front processor, solver and preprocessor, all can download in domestic each fatware websites such as software download center, the Pacific Ocean and 121 download websites; Involved grain size analysis software EDEM is the general CAE software being used for Simulation and analysis particle disposal and production operation of a modernization discrete element analysis tech design; it is the flagship product of DEM-Solutions company of Britain, and Beijing sea base technology & development Co. is EDEM China Software commission merchant.