Search Documents

By Dong-dong Liu, Jian-ping Niu, Qi Liu, Shu-sen Cheng

A numerical model was developed to simulate the three-dimensional heat transfer in the belly and lower shaft region of the blast furnace wall consisting of copper steel composite stave, packing and shell. The temperature distribution of copper steel composite stave without the slag layer or refractory lining was investigated and then compared with the copper stave. The results showed that copper steel composite stave can satisfy the requirement of overheatingfree and had an analogous heat transfer behavior as copper stave but more stable cooling intensity. This model had been validated by the good agreement between the computed and measured temperature in the thermal test.

Jan 1, 2015

By Dengfu Chen, Huabiao Chen, Yunwei Huang, Mujun Long, Youguang Ma, Lintao Gui

In the current work, a three-dimensional model coupling turbulent flow, heat and solute transport was developed to investigate the solute transport and redistribution in the mold of slab continuous casting where the Reynolds number is very high and the molten steel flows strongly. The conservation equations of momentum, energy and species for a multicomponent system which includes solute element C, Si, Mn, P, S were solved with the commercial software ANSYS Fluent. The parameters used in the model are based on the actual continuous casting. The fluid flow, temperature distribution and solute elements distribution were analyzed. Results showed that negative segregation occurred near the strand surface, and the species concentration reached a peak value at the solidification front during solidification. The segregation degree of solute element S is higher than the other solute elements because of its lower equilibrium partition coefficient. The species concentration in the liquid pool is homogeneous because of the turbulent flow effect.

Jan 1, 2015

By Y. Ren, B. Yang

The magnetic field characteristics due to electromagnetic stirring in the secondary cooling zone of an extra-thick slab during continuous casting are numerically determined using Maxwell's equations, with particular emphasis on the distribution of magnetic induction intensity and electromagnetic force, as well as the influence of process parameters on the magnetic field. The orientation of the electromagnetic force is the same when the electromagnetic stirrers at the 30° and 60° posistions of the arc-shaped slab are powered in the same direction; otherwise, the orientation is opposite. The magnetic induction at the centre of the section of the stirrer's central axis at the 30° positio, as well as the current, frequency, and fitting formulae, are as follows: B = {0.075I.(1/A)-0.3} mT and B = {-1.75f.(1/Hz)+22} mT. The electromagnetic force and current fitting formula at the central point of the B-B section is F = {0.03168I2.(1/A2)-03.762I.(1/A)+326.7} N/m3, and the fitting formula of electromagnetic force and frequency at the centre point is divided into two segments. The first takes the form of a parabola: = {-21((f-4.3).(1/Hz)2)+970} N/m3, and the second segment is negative linear: F = {-1.65f.(1/Hz)+51.5} N/m3.

Apr 16, 2024

By Tingyao Liu, Lihui Han, Qing Liu, Bao Wang, Teng Yang, Yong Sheng

This paper mainly investigates the mass transfer phenomena of rare-earth concentrate from leaching process in Rushton stirred tank reactor with rotate speed of 100-600 rpm. In the research, a certain amount of sphere particles with 50-450 μm in diameter is treated as the mineral particles. Comparison of laboratory experiment data with numerical simulation is performed. The results show that the critical impeller speed for complete dispersion is 300 rpm. And the maximum value of turbulence kinetic energy appeared 15 mm away from paddle, where the mixing time is decreased 55% and 80%, respectively, when mineral particles are injected into stirred tank through one or four filling ports. Moreover, with mixing process going on, the numbers of particles near the paddle region accounts for as much as 58% of the total amount. Meanwhile, following the particle diameter increasing, the particle concentration of the upper part of stirred tank reactor is increased.

Jan 1, 2015

By Yuting Xue, Brijes Mishra

The relaxation test on sandstone specimens showed typical behavior in pre-failure region and stepwise behavior in the postfailure region. Overall, more significant stress relaxation occurred within the failed specimens than the intact ones. Numerical simulations were conducted with pre-defined failure plane and with Voronoi tessellation to visualize the relaxation behavior. The model with pre-defined failure plane showed the key role of failure plane and asperity in simulating the step-wise post-failure relaxation behavior. The inhomogeneous stress distribution within failed specimen and the viscous deformation of intact rocks induced high stress concentration at the asperities. The observed step-wise relaxation initiated from the failure of asperity. Furthermore, the simulations with Voronoi tessellation showed time-dependent fracture development during relaxation in post-failure region. The presence of fractures completely changed the stress distribution. Stress concentration occurred at the front area of fractures and at the interacting areas between fractures. During relaxation, fractures still developed with time and the sudden significant increase in the “damage” coincided with the acceleration of stress relaxation leading to step-wise relaxation. Finally, the results showed the possibility of using residual strength as the long-term strength of failed rock.

Jun 25, 2020

By Yang Shuxin, Guo Qiliang, He Ping, Lu Lizhi, Yao Rui

"We combined stepwise regression and constrained regression statistical method to obtain the best boundary conditions for the calculation of in-situ stress field in Wushaoling tunnel using finite element method. Our numerical results reveal that: the main geologic factors contributing to the initial background stress are gravity and tectonics; the principal stress values would drop in fault zones; the angle between the maximal horizontal stress and the axis of Wushaoling tunnel is acute; Wushaoling tunnel tends to become stable when this tunnel has 2:1 ratio of the height to its width; the possibility of rock-burst exists under high stress.INTRODUCTIONStress state of engineering region is an important factor affecting the stability of project. Geological disasters in the project area are closely related to their stress environment. Wushaoling tunnel which is 20,050 meters long, lying along Lanzhou-Wuwei segment of Lan-Xin Railway in China, is the key project of Lan-Wu double-track railway. The region surrounding this tunnel is under a high stress environment. Accurate identification of the stress environment of the engineering area and the interaction between geomaterials and geotechnical engineering construction holds significant scientific and practical value for engineering design, construction, and operation.Many scholars explored and discussed the present-day tectonic stress through investigations of tectonic evolution history (Ding, & Lu, 1986), crustal deformation (Zoback, 1992), focal mechanisms (Yang et al., 2003, Xie, Cui, & Zhao, 2004), and so forth. However, the above methods would only provide the approximate direction of the principal stress in a broad region, and could not give the absolute stress values. In-situ stress measurement is one of the primary means to directly gauge the current crustal stress. However, in-situ stress measurement can only be applied in small number limited by measurement conditions."

Jan 1, 2015

By B. Q. Li

A numerical study of the developing fluid flow and thermal aspects during the staning phase of aluminum continuous casting processes is presented. The deforming finite element method, coupled with a quasi-Lagrangian formulation, is used to represent the growth of an ingot and the dynamic fluid flow and thermal behavior that evolves in the growing ingot. The mathematical formulation of the problem and the numerical treatment are discussed. The algorithm solves the free moving nodes as part of total unknowns with an automatic time step control for time integration. This permits an optimal computing efficiency while still keeping the computational errors under control. A selection of the results obtained using the algorithm is presented. The results are concerned with the evolution of fluid flow, temperature and solidification phenomena during a DC round ingot casting process, as it evolves from the beginning of casting operation to a steady state. These results show that the time-varying, complex fluid flow and thermal phenomena during continuous casting processes can be very well represented by the computational algorithm. The results also illustrate that, during the casting process, ~e fluid flow and temperature fields change drastically at the initial stage but evolve slowly afterwards

Jan 1, 1993

By Ruirui Wei, Mingrui Yang, Xuewei Lv

Calcium ferrite, as the main bonding phase in the sintering process, greatly influences the property of sintered ore. Ultrasonic vibration has been widely applied to the treatment of metals to improve the properties of metals. In this work, ultrasonic vibration was introduced into the calcium ferrite and numerical simulation was carried out by solving the Rayleigh-Plesset equation to investigate the influence of ultrasonic conditions, including frequency, power and initial radii of cavitation bubbles on the cavitation bubbling behavior in calcium ferrite melt. The results indicated that lower ultrasonic frequency, higher power, and an appropriate initial radius of bubble promotes ultrasonic cavitation.

Mar 1, 2018

By M. Patel, G. Djambazov, J. Shrimpton

"Removing zinc by distillation can leave the lead bullion virtually free of zinc and also produces pure zinc crystals. Batch distillation is considered in a hemispherical kettle with water-cooled lid, under high vacuum (50 Pa or less). Sufficient zinc concentration at the evaporating surface is achieved by means of a mechanical stirrer.The numerical model is based on the multi-physics simulation package PHYSICA. The fluid flow module of the code is used to simulate the action of the stirring impeller and to determine the temperature and concentration fields throughout the liquid volume including the evaporating surface. The rate of zinc evaporation and condensation is then modelled using Langmuir's equations. Diffusion of the zinc vapour through the residual air in the vacuum gap is also taken into account.Computed results show that the mixing is sufficient and the rate-limiting step of the process is the surface evaporation driven by the difference of the equilibrium vapour pressure and the actual partial pressure of zinc vapour. However, at higher zinc concentrations, the heat transfer through the growing zinc crystal crust towards the cold steel lid may become the limiting factor because the crystallization front may reach the melting point.The computational model can be very useful in optimising the process within its safe limits.IntroductionLead bullion is commonly refined by removing various impurities (e.g. copper, silver, bismuth, antimony, etc) sequentially. This in normally achieved by the addition of reagents which selectively react with the impurities and the resulting compounds removed by various means. To aid efficient mixing (ensure uniform reactions), batch mixing of materials via the aid of impellers is common practice in a number of industrial sectors. Excess zinc, added during the desilvering process can be removed by one of three methods: (1) Vacuum dezincing - During vacuum dezincing, a vacuum is drawn on the agitated molten lead within a hemispherical vessel (Figure 1). Vaporized zinc condenses on the inner top surface of the vessel. The temperature at which zinc commences to vaporize lowers owing to heating under a vacuum. This enables vaporized zinc to be recovered in the form of high purity crystals without oxidation. The zinc crystals are recovered once the vessel lid is removed and the vacuum is broken. (2) Chlorine dezincing - In chlorine dezincing, molten desilvered lead is reacted with chlorine gas, forming a surface layer of zinc chloride contaminated with small amounts of lead chloride."

Jan 1, 2004

By S. Viswanathan, A. S. Sabau

"Dimensional changes between a pattern die and its corresponding investment casting part occur as a result of complex phenomena such as thermal expansion/contraction and hot deformation (elastic, plastic, and creep) during the processing of the pattern material (wax), mold material (shell) and solidifying alloy. Determining the pattern tooling dimensions is crucial to the dimensional control of the investment casting process. To date, computational methodologies have not been applied to the prediction of dimensional changes in investment castings.This paper deals with the evaluation of wax pattern dimensions, which is one important factor in determining pattern tooling dimensions in investment casting. Cerita 29-51, an industrial wax, is considered in this study. The wax pattern dimensions are determined by considering its thermo physical and thermomechanical properties, restraint of geometrical features by the metal die, and process parameters such as dwell time, platen temperature, injection pressure, and injection temperature.IntroductionThe investment casting process involves making a disposable wax pattern by injecting wax into a metal mold, building a ceramic shell mold around the wax pattern by the application of a series of ceramic coatings to the wax pattern, dewaxing the ceramic shell mold, and casting the alloy into the shell mold. Dimensional changes between the pattern tooling and its corresponding cast part occur as a result of thermal expansion, shrinkage, hot deformation, and creep of the pattern .material (wax), mold material (shell), and solidifying alloy during the processing [1]."

Jan 1, 2001

In this study, numerical simulations of yield pillars employing the finite element method with time-dependent Drucker-Prager material model are conducted. The stress control mechanism of yield pillars on weak floor strata is studied. The development of pillar yield zone, the distribution of pillar stress and the development of failure zone in the weak floor strata are investigated. The study indicates that a properly selected yield pillar size may significantly reduce the peak stress induced in the weak floor strata and effectively control the development of failure zone in the floor.

Jan 1, 1999

By Congling Shi, Marek Borowski, Yue Wang, Jianwei Cheng, Xixi Zhang, Xuan Xu

As the first barrier to prevent accidents caused by explosions in underground coal mines, the integrity and stability of mine seals have great effects on mine workers’ lives and underground production safety. In this paper, the mechanical responding performance of typical mine seals under the impact load are studied by using numerical simulation method. By such simulation studies, the safety of different seals subjected to the explosion load is evaluated. Seals’stability are also investigated under effects caused by various influence factors including thickness of mine seal, cutting depth into surrounding rocks, rock type, and the roof-tofloor convergence under given the explosion load. The research results can provide theoretical basis and reasonable improvements for mine seal’s construction practices for improving its anti-explosion performance.

By Zengwu Zhao, Peng Xu, Jieliang Wang, Qingwei Bai, Xinyang Liu, Xiao Lei

The paper investigates the capture and recovery of rare earth minerals from Bayan Obo ore using superconducting magnetic separation. A finite element model (FEM) was established to study the influence of various factors on the capture of rare earth minerals under a high-gradient superconducting magnetic field. The study aims to improve the efficiency of rare earth mineral separation and recovery.

Oct 1, 2023

By Mark S. Rayson, Jessy Lee, David O. Johnson, Matthew R. Rawls

The ability to manipulate the structure of liquid chemicals in an explosive is important in mining applications, where precision blasting is crucial for the safe, reliable, economical and eco-friendly extraction of ores. Static mixers allow explosives, emulsion-based in particular, to be custom-made as they offer precise and reproducible mixing without the need for moving mechanical parts. In this paper, we perform numerical simulations on the laminar mixing of two nitrate-based fluids in a Kenics static mixer unit for various mixer parameters. Two numerical methods were investigated in the study of the mixing pattern: solving the mass transfer equation and using particle tracking. The former introduces artificial numerical diffusion which can lead to spurious smearing of the component boundaries, a problem that can be avoided with the latter technique. We have also demonstrated that the particle tracing method is capable of qualitatively predicting the blending patterns under the conditions studied.

Jan 1, 2014

By M. Kamyab

The Deep Exploration Technology project, as part of the Australian Cooperative Research Centre, is developing technologies to drill deep boreholes for mineral exploration purposes. The aim of the initiative is to transfer drilling technologies from oil and gas into mining applications. Several modifications have to be pondered due to this technology exchange: smaller cuttings sizes because of using impregnated diamond bits; narrower annulus space; higher drilling fluid velocity and likelihood of turbulent flow regime in the annulus; and higher pressure loss. In particular the results of our lab experiments have shown that the increase in concentration of small size cuttings can significantly change the rheology of the mud when drilling in hard rocks. This, in return, could have a noticeable increase in the pressure loss encountered along the annulus space. As a continuation of the above study, numerical simulations were conducted to investigate the cuttings transport behaviour in small size annulus space. Fluid velocity, annulus size, cuttings concentration, cuttings density, and cuttings size were some of the parameters that were varied to see their effects on cuttings transportation and the pressure loss along the annulus space. The results indicate that, due to the differences between oil and gas and mining drilling, conventional oil well planning is not applicable for mineral exploration drilling and numerical simulation would help in this matter.

Aug 1, 2013

By Valerio Viti

In High Pressure Die Casting a molten metal is injected through a thin gate into the cast cavity. High injection pressures and high gate velocities create atomization phenomena which can negatively affect the final quality of the cast. In order to control atomization process it is essential to understand the atomization patterns and the two-phase flowfield that exist at the gate exit. In the present work, 2D numerical simulations of the flow of molten Magnesium through a high aspect-ratio rectangular gate are performed using a water analogue for an open and closed cavity. The numerical simulations made use of VOF-type (Eulerian) physical models and Lagrangian models. Further, a sub-grid scale model was implemented that, in conjunction with VOF-type equations can efficiently predict the general atomization pattern without the need for high-resolution grids. The numerical results were compared to experimental data for validation.

Jan 1, 2006

By S. O. Choi, H. -S. Shin, W. -K. Song

Nuclear waste should be disposed with the best warranty for the national health and welfare. Among the several repository systems, the deep underground disposal is considered to be a most suitable system in Korean peninsula. However it demands basically the comprehensive stability analysis such as the behaviour of discontinuities in rock mass, the behaviour of ground water, and the reaction of rock stress according to the heat from the nuclear waste. For these purposes, the underground nuclear waste disposal repository was simulated in two- and three- dimensions. It shows that emphasis will be given to the thermal conductivity of discontinuities and the thermal convection of ground water as well as the thermal stress generated by the nuclear waste for design and construction of repository system, and they should be obtained from the field measurements. Déchets nucléaires devront tre déposés avec la meilleure garantie pour la santé et le bien public. Parmi les plusieurs systèmes de dépôt, la disposition en souterrin profond est considérée au système le plus approprié en Corée. Pourtant ce système demande essentiellement l’analyse éten-due de stabilité comme le comportement des discontinues en masse rocheuse, le comportement de l’eau souterrin et la réaction des contraintes de roche selon la chaleur depuis les déchets nucléaires. Pour ce but, le dépôt des déchets nucléaires en sous-sol a étè simulé en deux et trois dimensions. Il montre que l’accent sera donné à la conductivité thermale des discontinues et la convection thermale de l’eau souterrin aussi bien que la contrainte thermale générée par les déchets nucléaires pour le plan et la construction du systèmes de dépôt, et ils devront tre obtenus depuis les mesures sur place. Die Kernwüste muß mit der besten Gewähr für die nationale Gesundheit und die Wohlfahrt angeordnenet werden. Unter einigem Fundgrubesystem, das tief unterirdisch Wegwerfen werde in Betracht das passendest System in der koreanisch Halbinsel gezogen. Jedoch, es verlangt grundzätzlich die umfassend Stabilitätsberechnung so wie das Behnehmen mangelnder Konti-nuität in der Felsmasse, das Behnehmen des Grundwasser, und die Reaktion der Felsenbelastung nach der Hitze von der Kernwüste. Für diesen Zwecke, die unterirdische Wegferfenfundgrube der Kernwüste wurde in der zwei und drei Dimensionen gesimuliert. Es beweist daß die Betonung wird zur thermisch Leitfähigkeit mangelnder Kontinuität und thermisch Konvektion des Grund-wasser sowie die thermisch Belastung ergezeugt durch die Kernwüste für den Entwurf und die Konstruktion des Fundgrubesystem gegeben werden, und sie müssen von den Feldmaße ergelangt werden .

Jan 1, 2003

By M Zoorabadi

Overcoring is a common technique to measure in situ rock stress using the concept of stress relief. The strains measured during overcoring are used to calculate the stress tensor assuming continuous, homogeneous, isotropic and linear elastic behaviour. This approach does not recognise any effect of the thickness of the rock annulus produced during overcoring instead assuming that the stresses are entirely a function of the difference in strain between the start of the overcoring and the end with no regard to the transient strains produced. Three-dimensional numerical modelling was used to simulate the overcoring process and investigate the impact of core diameter on the transient stresses and strains during overcoring. This study showed that in thicker core, maximum transient axial stresses and strains occur when the drill bit approaches the measuring point and greater tensile axial stresses were produced when the annulus was thinner. CITATION:Zoorabadi, M, 2018. Numerical study of annulus thickness on overcoring stress measurements, in Proceedings The Fourth Australasian Ground Control in Mining Conference (AusRock), pp 274–278 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 30, 2018

By Ramana G. Reddy, Yudong Li

Numerical modeling of the thermal plasma process was carried out based on the thermal plasma reactor in our lab and confirmed using experimental data. The inlet boundary conditions of a non-transferred DC arc thermal plasma reactor were used in modeling the temperature and fluid flow distribution in the reactor. Different mesh grid sizes were used to confirm the model is independent of grid size. Temperature profile and gas flow distribution in the thermal plasma reactor were developed by the computational fluid dynamics (CFD) with ANSYS Fluent. The predicted temperatures are in good agreement with the experimentally measured temperatures in the reactor. The influence of plasma torch input power as well as the plasma gas flux on the temperature distribution was investigated using this model. The influence of power input and gas flow rate on temperature and velocity distributions are not independent. Generally, higher power input and lower gas flow rate will give rise to the temperature increase in the reactor.

Mar 1, 2017

By Liming Yuan

A computational fluid dynamics (CFD) study was conducted to model effects of coal properties on the potential for spontaneous heating in longwall gob (mined out) areas. A two longwall panel district using a bleeder ventilation system was simulated. The permeability and porosity profiles for the longwall gob were generated from a geotechnical model and were used as inputs for the three dimensional CFD modeling. The spontaneous heating is modeled as the low temperature oxidation of coal in the gob using kinetic data obtained from previous laboratory scale spontaneous combustion studies. Heat generated from coal oxidation is dissipated by convection and conduction, while oxygen and oxidation products are transported by convection and diffusion. Unsteady state simulations were conducted for three different US coals and simulation results were compared with some available test results. The effects of coal surface area and heat of reaction on the spontaneous heating process were also examined.

Jan 1, 2008