水辅助注射成型充填过程中工艺参数影响的数值模拟

基于Hele-Shaw流动模型,对一直管的水辅助注射成型的充填过程的工艺参数的影响进行数值模拟。通过单因素实验法研究了工艺参数对水穿透长度与其厚度分数的影响。结果表明,熔体注射量对水穿透距离影响最大,而熔体温度、模具温度、注水温度、注水压力对水穿透距离与其厚度分数的影响不是很明显。

Numerical simulation of LWD resistivity response of carbonate formation using self-adaptive hp-FEM

Most of the carbonate formation are highly heterogeneous with cavities of different sizes, which makes the prediction of cavity-filled reservoir in carbonate rocks difficult. Large cavities in carbonate formations pose serious threat to drilling operations. Logging-whiledrilling （LWD） is currently used to accurately identify and evaluate cavities in reservoirs during drilling. In this study, we use the self-adaptive hp-FEM algorithm simulate and calculate the LWD resistivity responses of fracture-cavity reservoir cavities. Compared with the traditional h-FEM method, the self-adaptive hp-FEM algorithm has the characteristics of the self-adaptive mesh refinement and the calculations exponentially converge to highly accurate solutions. Using numerical simulations, we investigated the effect of the cavity size, distance between cavity and borehole, and transmitted frequency on the LWD resistivity response. Based on the results, a method for recognizing cavities is proposed. This research can provide the theoretical basis for the accurate identification and quantitative evaluation of various carbonate reservoirs with cavities encountered in practice.

PREDICTION OF RIME ICE ACCRETION AND RESULTING EFFECT ON AIRFOIL PERFORMANCE

The roughness effect based on the wall function method is introduced into the numerical simulation of the rime ice accretion and the resulting effect on the aerodynamic performance of the airfoil. Incorporating the two-phase model of air/super-cooled droplets in the Eulerian coordinate system, this paper presents the simulation of the rime ice accretion on the NACA 0012 airfoil. The predicted rime ice shape is compared with those results of measurements and simulations by other icing codes. Also the resulting effects of rime ice on airfoil aerodynamic performance are discussed. Results indicate that the rime ice accretion leads to the loss of the maximum lift coefficient by 26%, the decrease of the stall angle by about 3° and the considerable increase of the drag coefficient.

NUMERICAL SIMULATIONS OF TEMPERATURE FIELD IN DIRECT METAL LASER SINTERING PROCESS

A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radiation, and convection. A formula for the calculation of the heat conductivity of a sintering system containing solid phase, liquid phase, and gas phase is given. Due to the continuous movement of the laser beam, a local coordinate system centered on the laser beam is used to simplify the analytical calculation. Assuming that it is approximately a Gaussian laser beam, the heat conduction equation is resolved based on the assumption of the thermal insulating boundary conditions and the fixed thermal physical parameters. The FORTRAN language is employed to compile the program to simulate the temperature field in the direct copper powder sintering process. It shows a good agreement with the preliminary experimental results.[KH3/4D]

Finite-difference numerical modeling with even-order accuracy in two-phase anisotropic media

To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.

Construction of reasonable pillar group for undersea mining in metal mine

In order to achieve the safe mining in Sanshandao Gold Mine,five schemes of secure pillar group are designed.Using the method of the renormalization group,the failure mechanism of the pillar group is explored,and the safety factor of the pillar system is also obtained.The displacement characteristics,stress-strain laws,distribution of plastic zone and damage range of different pillar group are analyzed using numerical calculation software FLAC3D.To determine a reasonable pillar group scheme,the pillar stability and roof deformation are utilized to evaluate the safety of the pillar group.In addition,the theory of fuzzy comprehensive evaluation is adopted to verify the optimal scheme.The pillar group with the lowest roof deformation value is chosen as the optimal plan,which renders a factor of safety of 2.06 for the pillar group.According to this scheme,pillars with the width of 10 m are set along the strike of undersea deposit with the interval of 50 m.Rib pillars of 15 m in width are set at the location of the exploration line of 127,151 and 167.The analysis can be used to provide guidance for optimal design of pillar structures in undersea mining.

Numerical Study of Solid-Liquid Two-Phase Flow in StirredTanks with Rushton Impeller(Ⅱ) Prediction of Critical Impeller Speed

Three-dimensional solid-liquid flow is mathematically formulated by means of the 'two-fluid' approach and the two-phase k-ε-AP turbulence model. The turbulent fluctuation correlations appearing in the Reynolds time averaged governing equations are fully incorporated. The solid-liquid flow field and solid concentration distribution in baffled stirred tanks with a standard Rushton impeller are numerically simulated using an improved 'inner-outer'iterative procedure. The flow pattern is identified via the velocity vector plots and a recirculation loop with higher solid concentration is observed in the central vicinity beneath the impeller. Comparison of the simulation with experimental data on the mean velocities and the turbulence quantities of the solid phase is made and quite reasonable agreement is obtained except for the impeller swept volume. The counterpart of liquid phase is presented as well.The predicted solid concentration distribution for three experimental cases with the average solid concentration up to 20% is also found to agree reasonably with the experimental results published in the literature.

Effect of electromagnetic force on turbulent flow of molten metal in aluminum electrolysis cells

The standard k-ε model was adopted to simulate the flow field of molten metal in three aluminum electrolysis cells with different anode risers. The Hartman number, Reynolds number and the turbulent Reynolds number of molten metal were calculated quantitatively. The turbulent Reynolds number is in the order of 103 , and Reynolds number is in the order of 104 if taking the depth of molten metal as the characteristic length. The results show that the molten metal flow is the turbulence of high Reynolds number, the turbulent Reynolds number is more appropriate than Reynolds number to be used to describe the turbulent characteristic of molten metal, and Hartman number displays very well that electromagnetic force inhibits turbulent motion of molten metal.

Numerical Simulation of the Whole Three-Dimensional Flow in a Stirred Tank with Anisotropic Algebraic Stress Model

In accordance to the anisotropic feature of turbulent flow, ananisotropic algebraic stress model is adopted to predict theturbulent flow field and turbulent characteristics generated by aRushton disc turbine with the improved inner-outer iterativeprocedure. The predicted turbulent flow is compared with experimentaldata and the simulation by the standard k-ε turbulence model. Theanisotropic algebraic stress model is found to give better predictionthan the standard k-ε turbulence model. The predicted turbulent flowfield is in accordance to experimental data and the trend of theturbulence intensity can be effectively reflected in the simulation.

Numerical simulation of influence of Langmuir adsorption constant on gas drainage radius of drilling in coal seam

To determine reasonable distance of gas pre-drainage drillings in coal seams, a solid–gas coupling model that takes gas adsorption effect into account was constructed. In view of different adsorption constants,the paper conducted the numerical simulation of pre-drainage gas in drillings along coal seam, studied the relationship of adsorption constants and permeability, gas pressure, and effective drainage radius of coal seams, and applied the approach to the layout of pre-drainage gas drillings in coal seams. The results show that the permeability of coal seams is on the gradual increase with time, which is divided into three sections according to the increase rate: the drainage time 0–30 d is the sharp increase section;30–220 d is the gradual increase section; and the time above 200 d is the stable section. The permeability of coal seams is in negative linear and positive exponent relation with volume adsorption constant VLand pressure adsorption constant PL, respectively. The effective drainage radius is in negative linear relation with VLand in positive exponent relation with PL. Compared with the former design scheme, the engineering quantity of drilling could be reduced by 25%.

Effect of heterogeneity on mechanical and acoustic emission characteristics of rock specimen

The influence of heterogeneity on mechanical and acoustic emission characteristics of rock specimen under uniaxial compress was studied with numerical simulation methods.Weibull distribution function was adopted to describe the mesoscopic heterogeneity of rocks.The failure process of heterogeneous rock specimen under uniaxial loading was simulated using FLAC 3D software.Five schemes were adopted to investigate the influence of heterogeneity.The results demonstrate that as the homogeneity increases,the peak strength and brittleness of rocks increase,and the macro elastic modulus improves as well.Heterogeneity has great influence on macro elastic modulus and strength when the homogeneity coefficient is less than 20.0.The volume expansion is not so obvious when the homogeneity increases.As the homogeneity coefficient increases the acoustic emissions modes change from swarm shock to main shock.When the homogeneity coefficient is high,the cumulative acoustic emission events-axial strain curve is gentle before the rock failure.The numerical results agree with the previously numerical results and earlier experimental measurements.

Mathematical modelling and numerical optimization of particle heating process in copper flash furnace

A mathematical model of the particle heating process in the reaction shaft of flash smelting furnace was established and the calculation was performed.The results indicate that radiation plays a significant role in the heat transfer process within the first 0.6 m in the upper part of the reaction shaft,whilst the convection is dominant in the area below 0.6 m for the particle heating.In order to accelerate the particle ignition,it is necessary to enhance the convection,thus to speed up the particle heating.A high-speed preheated oxygen jet technology was then suggested to replace the nature gas combustion in the flash furnace,aiming to create a lateral disturbance in the gaseous phase around the particles,so as to achieve a slip velocity between the two phases and a high convective heat transfer coefficient.Numerical simulation was carried out for the cases with the high-speed oxygen jet and the normal nature gas burners.The results show that with the high-speed jet technology,particles are heated up more rapidly and ignited much earlier,especially within the area of the radial range of R=0.3−0.6 m.As a result,a more efficient smelting process can be achieved under the same operational condition.

Kinetic Study on Extraction of Red Pepper Seed Oil with Supercritical CO_2

A mathematical model for extraction of red pepper seed oil with supercritical CO2was proposed.Some factors influencing the process were investigated,including operation pressure,temperature and extraction yield Xe(%).The model was solved by the method of weighted residuals and used to simulate the process numerically.The kinetic equation is expressed as Xe16.8606exp(t/9.98177)16.95457 and the simulation results are in excellent agreement with the experimental data.The optimal operating parameters are 30 MPa,318 K and 60 min.

Approach to Generalized Synchronization with Application to Chaos-Based Secure Communication

A constructive theorem is established for generalized synchronization (GS) related to C1 diffeomorphic transformations of unidirectionally coupled dynamical arrays. The theorem provides some interpretations about the underlying mechanism of various GS phenomena in nature. As a direct application of the theorem, a chaos-based secure Internet communication scheme is proposed. Moreover, a cellular neural network (CNN) of Chen's chaotic circuits with GS property is designed and studied. Numerical simulation shows that this Chen's CNN has high security and is fast and reliable for secure Internet communications.

Comparison of numerical simulations and experiments in conical gas–solid spouted bed

Flow behavior of gas and particles in conical spouted beds is experimentally studied and simulated using the twofluid gas-solid model with the kinetic theory of granular flow.The bed pressure drop and fountain height are measured in a conical spouted bed of 100 mm I.D.at different gas velocities.The simulation results are compared with measurements of bed pressure drop and fountain height.The comparison shows that the drag coefficient model used in cylindrical beds under-predicted bed pressure drop and fountain height in conical spouted beds due to the partial weight of particles supported by the inclined side walls.It is found that the numerical results using the drag coefficient model proposed based on the conical spouted bed in this study are in good agreement with experimental data.The present study provides a useful basis for further works on the CFD simulation of conical spouted bed.

NORMALIZED DIFFERENCE SNOW INDEX SIMULATION FOR SNOW-COVER MAPPING IN FOREST BY GEOSAIL MODEL

The snow-cover mapping in forest area is always one of the difficult points for optical satellite remote sensing. To investigate reflectance variability and to improve the mapping of snow in forest area, GeoSail model was used to simulate the reflectance of a snow-covered forest. Using this model, the effects of varying canopy density, solar illumination and view geometry on the performance of the MODIS (Moderate-resolution Imaging Spectroradiometer) snow-cover mapping algorithm were investigated. The relationship between NDSI (Normalized Difference Snow Index), NDVI (Normalized Difference Vegetation Index) and snow fraction was discussed in detail. Results indicated that the weak performance would be achieved if fixed criteria were used for different regions especially in the complicated land cover components. Finally, some suggestions to MODIS SNOWMAP algorithm were put forward to improve snow mapping precision in forest area based on the simulation, for example, new criteria should be used in coniferous forest, that is, NDSI greater than 0.3 and NDVI greater than zero. Otherwise, a threshold on view zenith angle may be used in the criteria such as 45°.

Finite Element Numerical Simulation of Ground Subsidence in Liangjia Colliery

Being aimed at the ground subsidence due to underground coal mining,a numerical model of rock was established and an appropriate method of numerical simulation was put forward.Using the measured subsidence data on the ground,the equivalent mechanical parameters of the rock stratums can be back-calculated by the properly treatment of coal excavation area,then the ground subsidence of other coal mining area can be predicted by FFM.It provided reference for the treatment of the buildings on the ground of this colliery.

Turbulence Model Investigations on the Boundary Layer Flow with Adverse Pressure Gradients

In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients （APGs） is conducted by using Reynolds-averaged Navier-Stokes （RANS） equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox （2006） k-co are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.

THE INFLUENCE OF SCREEN PARAMETERS ON PERFORMANCE CHARACTERISTICS

A series of laboratory screening tests have been performed with the objective of evaluating the effect of feed rate, screen length and screen motion on the quality of separation. The separation performance has been evaluated by the use of proposed ISO Standard measures which are based on the partition curve.The laboratory rig, the test programme and results derived from it are discussed. The effect of screen length, feed rate and screen motion on the levels of screen performance as indicated by these measures and a normalised separation size parameter are discussed in relation to the assumed screening conditions for the tests reported and conclusions are summarised.

Numerical and experimental investigations of weld pool geometry in GTA welding of pure aluminum

A 2-D numerical model was developed to predict the shape of weld pool in stationary GTA welding of commercial pure aluminium, without considering fluid flow in the weld pool. A Gaussian current density and heat input distribution on the surface of the workpiece were considered. The parameters of Gaussian distribution were modified by comparing calculated results with experimental ones. It was found that these distribution parameters are fimctions of applied current and arc length. Effects of arc length, applied current and welding time on the geometry of the weld pool were investigated. To check the validity of the model, a series of experiments were also conducted. In general, the agreement between calculated overall shape of the weld pool and the experimental one was acceptable, especially in low applied currents. Therefore, it can be concluded that in pure aluminium, the heat conduction is dominant mechanism of heat transfer in the weld pool.