Search Documents

By H. Singh, Anupam Agrawal, H. S. Grewal

"Erosion in industrial appliances causes huge loss of valuable resources especially when solid particles are present in the working medium. A number of erosion models that are available in the literature claimed to predict the erosion response. In this work, five erosion models have been evaluated using computational fluid dynamics (CFD) methodology. Three most commonly used experimental test rigs (solid particle and unsubmerged slurry jet type and pot type) have been simulated and comparisons with the experimental results have been performed. The study helps in understanding the effect of fluid medium and flow pattern on the erosion process. The knowledge acquired shows that phenomenological based erosion models give better results considering the wide range of operating parameters. The in-depth study of the flow field shows the possibility that erosion mechanism, being parameter dependent, could change with location on the same specimen.1.0 IntroductionErosion is a surface degradation phenomenon which causes the removal of material due to impacts either of solid particles and/or bubbles or gas entrained in a fluid medium. Erosion has various categories depending upon the erosion medium such as, solid particles erosion (solid particles entrained in gaseous medium), slurry erosion (solid particles entrained in liquid medium) or cavitation erosion (erosion caused by exploding gas or vapor bubbles). Humphrey [1] has listed 21 different factors that decide the erosion response of the components under a given working condition. The number of factors involved in the erosion process gives us an idea about its complex nature."

Jan 1, 2012

Numerical simulation of the fluid flow and heat transfer has been carried out for a waste-heat boiler in an Outokumpu type copper flash smelter. In this boiler, SO2-containing off-gas flows into the waste-heat boiler radiation section at about 1350°C. A fraction of cooler off-gas from the electrostatic precipitator at about 200°C was recirculated into the radiation section, and a small amount of ambient air was leaked into the boiler from the bottom. The computational results show that the flow is strongly recirculating and distinctly three dimensional in most of the boiler, particularly in the radiation section. The hot off-gas cooled down rapidly after entering the radiation section due to radiation. The predicted flow pattern and temperature distribution were in a good agreement with laboratory models and industrial measurements.

Jan 1, 1995

By Menglong Bian, Xiaotong Dong

Thermal damage is an urgent problem faced by deep mining. High-temperature rock mass, hot water and electromechanical equipment are the sources of the thermal damage. The increase in air temperature in mine tunnel is the result of fluid–solid coupling heat transfer from various heat sources. This paper takes a roadway with hot water and a mine car as the research object. Based on the Fluent, Static Structural and Steady-State Thermal modules in ANSYS, the thermal-fluid–structure coupling analysis in the roadway is carried out. It is found that the hot water has a significant effect on the thermal environment of the mine. Under the experimental parameters, the maximum heat flux density and equivalent elastic strain at the interface between the hot water and the rock are 9789.3 W/m2 and 0.00054532 m/m, respectively. The heat flux distribution in the rock area near the mine car does not change much, and the equivalent elastic strain value in this area is high. After the air is heated in the hot water area, the temperature curve changes gradually and then rises steadily. By changing the inlet air temperature parameter, it is found that the lower the inlet air temperature, the better the cooling effect, and the change of the inlet air temperature can change the air temperature field in the roadway as a whole. By changing the parameter of inlet air speed, it is found that higher wind speed brings better cooling effect, but with the increase in wind speed, the cooling effect of roadway thermal environment is very limited. When the wind speed increases from 5 m/s to 10 m/s, the air temperature of most measuring points drops by more than 2℃; when the wind speed increases from 15 m/s to 20 m/s, the air temperature of most measuring points drops below 0.7℃. After the wind speed reaches a certain value, the cooling effect of ventilation is very limited.

May 5, 2022

By M. K. Ameen, R. El-Hacha

This paper presents the analytical results obtained from numerical simulation of reinforced concrete slabs strengthened with various FRP systems and subjected to blast loading. A three-dimensional finite element model using AUTODYN software was used to predict the overall behavior of the slabs. The numerical results were validated by comparing with experimental results from a program carried out to investigate the efficiency of different strengthening systems on the blast resistance improvement of small-scale reinforced concrete (RC) slabs. The experimental program consisted of testing nine RC slabs strengthened with three different externally bonded (EB) FRP systems attached at the bottom side of the slabs and one RC slab tested without strengthening for comparison purposes. The EB FRP strengthening systems investigated included two layers of Carbon Fiber Reinforced Polymer (CFRP) sheets, two layers of Glass Fiber Reinforced Polymer (GFRP) sheets, and one layer of Glass-Fiber Resin Polymer (GRP) sheet. Different epoxy resin, patching mortars and impregnating epoxy adhesives were used for the slab surface preparation and FRP system application of the strengthening layers. The aim was to evaluate the residual strength of the test specimens after exposure to blast loads. The blast loading parameters were evaluated for each test using AUTODYN simulation and compared with the computer program CONWEP. All specimens were subjected to same spherical TNT charge weight and stand-off distance with the respect to the slab top surface. Extensive comparison is made between the experimental results and the AUTODYN simulations for the control unstrengthened slab and the strengthened slabs including time history records of the bottom steel strain and mid-span acceleration and the peak dynamic deflection of the slabs.

Jan 1, 2010

By Y. Jin, H. You, H. Pang, B. Lin

Successful application of steam assisted gravity drainage (SAGD) on the recovery of land facies ultra-heavy oil in Xinjiang oilfield depends largely on shortening the preheating period with the aid of the micro-fracturing technology, which is conducted through controlled water injection to a pair of horizontal wells to create homogenously distributed microcracks in the zone between the two wells. The numerical simulation of this mechanical process is yet to be investigated. Up to date, there has been research work on rock mechanics during the SAGD production period of marine deposited ultra-heavy oilsands, but the counterpart research on land facies ultra-heavy oilsands of Xinjiang oilfield is left undone. In this regard, this study is attempted to perform numerical simulation of the fracturing process on land facies unconsolidated oilsands. With this objective, the study first investigates geological and physical identification of the field collected cores, then carries out mechanical analysis on laboratory triaixal test results, based on which numerical simulation of the fractured zone is performed. The simulation is applied through a finite element platform adopting extended Drucker-Prager constitutive model, with the parameters determined by experimental work. It is verified that hydraulic fracturing of unconsolidated oilsands generates a microcrack feature in terms of a dilative zone instead of linear apertures. The extent to which the fractured zone propagates is visually demonstrated based on the magnitude distribution of equivalent plastic strain or volumetric dilation. Evaluation of the connection of injection and production wells can possibly be determined by the expansion of the saturated region around each well. Finally, the implications on the efficiency of hydraulic fracturing in stimulation of the pay zone are investigated.

Jan 1, 2015

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rockmass strength. Past studies have largely used arbitrary scaling approaches to produce rockmass strength from laboratory rock tests. This paper proposes a stochastic approach to produce rockmass strength from laboratory rock tests. Based on the laboratory data, MATLAB with an extreme value stochastic model generated a database for each physico-mechanical property of the rock. Then finite-difference software FLAC2D simulated laboratory-sized rock specimens. Random material properties in MATLAB were generated using the grids developed in the numerical model. The random properties were then applied to the final FLAC2D model. Model runs simulated the approach performed in the laboratory. The results from the model indicate that a stochastic approach produced compressive strengths that were lower than those from the deterministic approach. Failure modes for each specimen were different from those observed in the laboratory. In addition, random density also influenced the failure mode, highlighting the importance of stochastic analysis in rocks.

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rock strength. For accurate stability prediction, intact rock behavior is imperatively included in rock mass behavior. However, past research largely used arbitrary scaling approaches to produce rock strength. This paper considers a stochastic approach in order to determine the strength of a rock. Based on the laboratory data, a MATLAB® with extreme value stochastic model generates a database for each physico-mechanical property. Then FLAC® simulates laboratory-sized rock specimens. The grids developed in the numerical model can in turn develop random material properties in MATLAB®, which researchers then applied to the final FLAC® model. Model runs simulate the approach performed in the laboratory. The results from the model indicate that a stochastic approach produces strengths that are lower than a deterministic approach. Failure modes for each specimen are different, similar to observations in the laboratory. In addition, random density also influences the failure mode, highlighting the importance of stochastic analysis in rocks.

By D. F. Malan, J. A. L. Napier

A number of shallow coal or hard rock mines employ pillar mining systems as a strategy for roof failure control. In certain platinum mine layouts, pillars are designed to ‘crush’ in a stable manner as they become loaded in the panel back area. The correct sizing of pillars demands some knowledge of the pillar strength and the overall layout stress distribution. It is particularly important to understand the impact of the layout geometry on the effective regional ‘stiffness’ of the rock mass around each pillar. An important design strategy is to model relatively detailed layout configurations which include a precise representation of the local pillar layout geometry and to analyse multiple mining scenarios and extraction sequences to select optimal pillar sizes and barrier pillar spacing. Although computational solution techniques are now impressive in terms of run time efficiency, a major difficulty is often encountered in assigning suitable material properties to the pillars and in devising an effective material description of the layered rock strata overlying the mine excavations. This paper outlines an efficient numerical strategy that can be used to assess large-scale pillar layout performance while retaining the ability to modify individual pillar constitutive behaviour. The proposed method is applied to selected layouts to compare estimated average pillar stress values against values determined by detailed modelling and against observed behaviour.

May 3, 2023

By R. L. Lewis

The White Pine Mine is the second largest active underground copper mine in the United States, covering an area of approximately 34 km2 in the Upper Peninsula of Michigan. A numerical flow model of the White Pine Mine was constructed to evaluate engineering controls that would reduce the rate of mine water discharge after closure. Of the passive controls evaluated, the placement of low-flow bulkheads within the mine workings was predicted to offer the greatest benefit, reducing the rate of discharge by as much as 35 percent.

Jan 1, 1996

By Cheng Li, Yanlei Wang, Lingrui Kong, Jianxin Tang, Fang Yuan

Coal mined-out areas have a significant impact on mining safety and affect the pressure distributions of neighboring coal faces. The PFC3D software is used to study the goaf mechanical characteristics, and the core loading test verifies the rationality of the numerical simulation method. The results from the numerical modeling showed that the void ratio, particle size ratio, and particle size distribution of the rock particle sample (RPS) had a significantly remarkable influence on the compression deformation characteristics of the rock block. The parallel-bond tensile strength and parallel-bond cohesion had a notable impact on the compression deformation traits of the rock block only when the axial stress level was lower than 12 MPa. By contrast, when the axial stress exceeded 5 MPa, the effective modulus significantly influenced the compression deformation features. However, the stiffness ratio, loading rate, and parallel-bond friction angle had almost no effect on the compression deformation characteristics. Furthermore, when the stress level was less than 10 MPa, the stress–strain curves from the rock core samples’ numerical simulation and loading tests were approximately linear. In addition, the Pearson correlation coefficients of the stress–strain curves suggest that the goaf zone’s core numerical simulation results had a robust correlation with the loading test results. The findings of this study can serve as a reference for the field application of PFC3D software and safety assessments in goaves.

Jan 18, 2022

By Toshio Suzuki, Machiko Ode, Won Tae Kim, Seong Gyoon Kim

"The Ostwald ripening and growth and coarsening of the secondary dendrite arm are studied using a phase-field model. The phase-field parameters are obtained at a thin interface limit proposed by Kim et al. In Ostwald ripening analysis, the difference of the concentration due ·to· the curvature effect is shown. The third power of the average particle radius Rav is proportional to time and the ripening constant shows a good agreement with the value by Deppers et al. The particle radii distribution shows the transition behavior to attain a self-similar stationary shape. The growth and coarsening behavior of the secondary dendrite arm is also analyzed using a phase-field model. The estimated arm spacing and its dependence on time are in good agreements with experimental data.IntroductionThe ripening or coarsening process controls the microstructure during the late stage of solidification. The ripening phenomena are recognized firstly by the Ostwald [1]. The first quantitative analysis for the Ostwald ripening in the limit of vanishing solid fraction was provided by Lifshitz, Slyozov and Wagner (LSW) [2,3]. In their study, that the mean radius grows proportional to the cubic root of time and the size distribution of precipitations shows a ""self-similar"" shape are shown. Most of the studies for the Ostwald ripening confirm or predict these two general rules above. In case of finite solid fraction, the general rules are observed but the ripening constant is larger and the self-similar shape is wider than the LSW prediction respectively [4-6]. To analyze· the high solid fraction system, the theories that consider the solute field correlation are proposed [4,7-8]. In the theory, however, since the particle is treated as the point, the effect of the concentration field deviation from the monopole field should be taken into account. Recently the numerical simulation, in which not only diffusive transportation of solute but also convection effect is considered, has been done using a phase-field model [9]."

Jan 1, 2000

By Guo Wenbing

"INTRODUCTIONThere is a great amount of coal (about 140 billions tons) left unmined under the surface structures, water bodies, railways (referred to as the “three-body”). Coal mining under “threebody” especially that under surface structures has become a major problem in the mining area. The key problem of mining under surface structures is to control overburden strata and surface movements. At present, the major surface subsidence control methods are strip pillar mining, backfilling mining method, room and pillar mining, grouting of bed separations, harmonic mining, and so on. Wongawilli mining method is a high-efficient shortwall mining method developed based on the room and pillar mining, and named after the first successful trial in Australia “Wongawilli” coal seam. The room and pillar mining equipment such as continuous miner is employed in this method. Its major advantage is flexible face layout. It can be used to mine the small wedge coal blocks and those that are difficult to be mined with retreat longwall mining method. It requires less capital investment, short lead time for production, flexible equipment operation and face move and high productivity. It has been applied in some China’s coal mines.The “Wongawilli strip coal mining technology” is a new coal mining method for mining under surface structures. It combines the strip mining longwall layout and Wongawilli high efficient mining technology. This technology can fully utilize the advantages of both the strip pillar and Wongawilli mining methods, and overcomes the disadvantages such as frequent face move, and low mining efficiency in strip pillar mining, and poor ventilation condition and poor long-term stability problems of coal pillars. In this paper, the Wongawilli strip mining method was designed for Wangtaipu mine and its feasibility was studied using numerical modeling and surface subsidence prediction.Geological and mining conditionsAt present, most of the No.15 coal seam reserves totaling 377.07 million tons in Wangtaipu mine are located under surface structures. In order to mine the coal under the surface structures, the Wongawilli strip mining method was employed. An appropriate panel was selected for trial and the surface movement monitoring lines were set up. The trial area is located in panel XV2214 (East) ( re 1). In this area, the floor elevation of coal seam is +634~646m, the surface elevation is +810~817.5m, and the average depth of coal seam is about 174m.The average panel length is 282m and average panel width is 73.5m. The seam dips 1~3° and is 1.8~2.5m thick with an average of 2.15m. The roof and floor rock characteristics are shown in Table 1."

Jan 1, 2015

By T. G. Sitharam, Resmi Sebastian

"Wave propagation studies across jointed rocks have significance in the fields of rock mechanics and geophysics. The passage of a seismic wave induces strains in rocks and often, these strains are of low magnitudes except at the vicinity of source. These low magnitude strains are generally in the elastic range and hence elastic constants necessary for the design and analysis of strength of rock mass, can be derived from wave velocities. Measurement of travel time velocities and measurement of resonance frequencies are the two methods usually adopted for obtaining wave propagation velocities.Resonant column apparatus measures wave velocities from the resonant frequencies obtained during testing. Depending on the type vibration applied, torsional or flexural, shear and compression wave velocities are secured. By varying the amplitude of vibrations applied, various strain levels can be introduced and corresponding resonant frequencies and wave velocities can be calculated. The resonant column experiments were performed on a model material, hardened Plaster of Paris (POP) samples, in laboratory. The paper explains the preparation of samples and the working mechanism of resonant column apparatus along with the estimation of wave velocity from resonant frequency. The numerical simulation of resonant column apparatus using distinct element method and its validation with the experimental results are presented in the paper.A distinct element code is usually used to simulate the movement of materials with intersecting discontinuities. Finite rotations and displacements of individual blocks representing the material are permitted and new contacts between the blocks are identified as the calculation progresses. Three dimensional Distinct Element Code (3DEC), the software intended for the numerical simulation of rock engineering projects, can well simulate the behavior of jointed rocks under static and dynamic loading conditions. The numerical simulation of resonant column apparatus has been performed with 3DEC considering its efficiency in modeling jointed rocks.The numerical modeling of the resonant column apparatus with application of torsional and flexural vibrations to the sample is elaborated in the paper. The response of samples to a range of frequencies is obtained and the frequency that induces maximum displacement in the sample is measured as the resonant frequency, for any particular strain level. The displacements measured at various points of the sample during application of vibrations are compared and presented. The numerical modeling of resonant column apparatus with an intact sample (without joint) and a sample with a joint incorporated are explained at a particular strain level. The resonant frequencies obtained from laboratory experiments and numerical simulations are compared and illustrated. It has been found that resonant frequencies obtained using numerical and laboratory experiments are in good comparison and can be extended for various strain levels and sample conditions."

Jan 1, 2015

"A numerical procedure based on Smoothed Particle Hydrodynamics (SPH) is developed for rock blasting phenomena involving multiple intersecting joints. In this procedure, an elasto-plastic damage constitutive model has been applied in the framework of SPH to simulate yielding response of the surrounding rock medium. In this paper, a square shape rock medium with centrally located emulsion explosive in a borehole is numerically blasted with and without the presence of two diagonal rock joints. The model is simulated by assuming all four sides of the rock medium as free surfaces. The paper elaborates the propagation of stress wave, expansion dynamics of gaseous products, and development of fracture, fragmentation of the medium and dissipation of energy due to damage.INTRODUCTIONDamage prediction of rock mass and its mean fragmentation size are of important to design blast parameters for a given explosive. Current blasting practice solely depends on empirical relationships between fragmentation sizes to those blast design parameters while mathematical descriptions of the physical laws governing the process of rock blasting are limited. A constitutive relationship based on properties of explosives and rock mass is of special interest because of its capabilities in describing the key features of rock fracture and fragmentation, resulting in understanding of the physical process involved in rock blasting. It is well known that blast-induced damage and fragmentation of rock mass occur due to detonation induced stress wave and product gas driven fracture propagation. Nevertheless, the fundamental mechanism of the failure process and the behavior of rock characteristics under these two different loading are not well understood in many cases, which are closely related to growth of fractures and interaction of cracks under different loading conditions.Continuum damage mechanics is considered to be appropriate to describe the failure process under such high-amplitude stress wave. Based on this approach, several theoretical models have been developed to analyze the dynamic damage of brittle material such as rock and concrete with microstructure (Chen, 1999; Taylor et al., 1986; Wu et al., 2004; Yazdchi et al., 1996). These models are mostly based on the theory of damage mechanics or fracture mechanics with some statistical treatment to incorporate random distribution of microcracks. In the application of fracture mechanics theory, sometimes it is difficult to estimate the micromechanical parameters such as crack density and fracture toughness in different modes. In this study, damage in rock material is evaluated according to the earlier model developed by Grady and Kipp (1980), which is a fracture model coupled with a material description for stress wave propagation to predict quantitatively fracture and fragmentation under explosive loading conditions."

Jan 1, 2015

By B. Panton, J. Carvalho, E. T. Brown

"A major consideration in the design of high capacity tiedown anchors for dam, bridge and tower foundations is the tensile resistance of the rock mass to pullout, typically as a result of overturning moments or hydrostatic uplift. Rock mass pullout capacity for installed anchors is developed from the tensile strength and fracture propagation properties of intact rock, the orientation and physical properties of the discontinuities and anchor confinement at depth. The typically assumed, but generally conservative, design approach is to calculate anchor pullout capacity using the dead weight of a uniformly shaped inverted “cone” with an assumed initiation point and breakout angle. In some cases, some tensile or shear strength is considered across the area of the assumed pullout cone, however, these strength estimates are based primarily on based on limited historic testing data, empirical relationships for rock mass strength developed for other applications, and engineering judgement.As an alternative to the current foundation anchor design method, this paper presents a new approach to the rock cone pullout problem using Discrete Fracture Networks (DFN) combined with numerical simulations. The simulations presented in the paper investigate the influence of fractures in a synthetic rock mass on ultimate anchor strength. In addition, the global influence on fracture network statistics on rock mass pullout failure is also discussed and it is contended that some of the aleatory and epistemic uncertainty inherent in the design of rock anchors can be reduced by considering fracture statistics in anchor design."

Jan 1, 2015

By A. J. Basu

Rock mass properties testing for determining strength parameters is an integral part of geomechanics. These tests are performed either in a laboratory environment or in a field (in-situ). The authors share their experience in using numerical modelling in simulating the laboratory test procedures. This approach helps studying the rock breakage and failure mechanism more efficiently and provides a high degree of visualization. Moreover, imperfections and discontinuities can be easily incorporated in these models. This capability of providing arbitrary fractures and discontinuities facilitates in building a large knowledgebase of rock strength and failure mechanisms. This knowledgebase may then be used to generate heuristics for determining rock strengths for rock mass with fractures and discontinuities.

Jan 1, 1997

By D. P. Roberts

This investigation aimed to simulate the experiments performed by Nic Gay (Gay, N.C.1976. Fracture growth around openings in large blocks of rock subjected to uniaxial and biaxial compression. International journal of Rock Mechanics and Minining Sciences and Geomechical Abstracts, vol. 13. pp. 231-243) on fracture growth around openings in blocks of rock subjected to uniaxial and biaxial compression, using the Elfen discrete element code. The results of the physical experiments indicated a number of trends that were successfully replicated in the numerical simulations, including crushing ahead of the face, the formation of a sub-vertical tensile fracture at centre-span, the formation of groups of fractures originating from the crushed face zone, the formation of large corner-to-corner fractures which precede sample failure, and linking of slot corner and sample corner fractures leading to sample failure. In addition, other features were apparent that were not reported by Gay. It is believed that these took the form of conjugate shear fractures (as reported by Gay in a previous paper) and other features that may have been altered or inhibited by the presence of pre-existing structure in the physical samples. It was found that the samples, though nominally made up of the same rock type, must have had significantly differing properties. The presence of pre-existing structures may have contributed to this discrepancy. Strain measurements from the physical experiments were reasonably well matched, except where the gauge lay directly in the path of advancing fractures.

Jan 1, 2012

By Bernhard Hiebl

Global trends that are constantly challenging the iron and steel producers are related to worsening of the iron ore grades, higher substitution rates of main iron ores to cheaper raw materials due to economic pressure and increased environmental standards. Changing the raw materials in an operating plant can significantly influence the sintering or pelletizing process. Adjusting all process parameters to regain stable and high productive operation takes considerable time and causes production losses. Furthermore changes in raw materials lead to variations in the process gas flows and emission concentrations.Therefore SVAI developed two models for numerical simulation of sintering process and pelletizing process. Amongst others the models are based on Ergun equation and include chemical reactions of the related process. For already operating plants the simulations shall predict changes of the process parameters when raw materials are exchanged or when operational parameters are varied. Furthermore the simulations give advantageous information for the design of green field plants and allow the simulation of different plant configurations.The objective of this paper is to demonstrate what the latest developments are in these numerical simulations and what Siemens VAI is doing to overcome the remaining challenges in this field.

Aug 2, 2018

By Li. Jiarong, Zuojian Liang, Qingyan Xu, Baicheng Liu, Shizhong Liu, Hailong Yuan

Single crystal solidification of turbine blades is a key technology for the production of reliable and high efficiency gas turbines. The solidification process of single crystal investment castings causes a number of problems that have to be solved such as the increase of casting defects and complexity of production process. Computer modeling and simulation is helpful to predict the defects and improve the process design. In the paper, a mathematical model for three dimensional simulation of solidification and single crystal growth of investment castings was developed based on CA, mechanism of crystal growth and basic transfer equations such as heat, mass transfer equations. Many factors including constitutional undercooling, curvature undercooling and anisotropy, which had vital influences on microstructure evolution, were considered in the model. Temperature, crystal growth interface and structure defects of investment castings of single crystal superalloy were investigated at withdrawal rates of 4.5mm/min and 7mm/min. Simulated results showed that microstructure evolution and structure defects were predicted reasonably and calculated results were coincident with experimental ones. The study indicated that three dimensional solidification simulation technique is a powerful tool for understanding the fundamental of solidification and the formation of structure defects in single crystal investment castings and finally for optimizing process design.

Jan 1, 2004

By Miaoyong Zhu, Sen Luo

"A two-dimensional CAFD model was developed for simulating the solute diffusion-controlled dentritic growth in the low Peclet number regime. Mass conservation and local equilibrium at the solid/liquid interface were solved to determine the growth velocity of interface with consideration of constitutional and curvature undercooling. The modified decentered square growth proposed by Wang et al was adopted to capture the first eight neighboring cells. The comparison between the classic Lipton-Glicksman-Kurz (LGK) analytical model for steady state tip variables and the simulated values showed that the present model was capable of quantitatively determining the dendritic growth, and the simulation tests demonstrated that the effects of the mesh anisotropy in growth kinetics and crystallographic orientation were eliminated. The numerical simulations were performed to the dendritic growth ofFe-0.6%C binary alloy in undercooled melt, and the results showed that the present model could predict the dendrite morphology with artificial growth orientation in acceptable level.IntroductionDuring the solidification of metals and alloys, the dendrite grows from the substance or the undercooling melt and solute redistributes between the dendrite and the interdendritic liquid. The solidification microstructure is close related to the phase and solute evenness and has a close linkage with the properties of final production.In the past few decades, various kinds of deterministic and stochastic models have been developed to predict the evolution of microstructures in solidification, with the fast development of computer power. Among these methods, Phase Field (PF) model [1, 2] has been known as one of the most adequate deterministic models for directly simulating the dendritic growth. However, there is a vital drawback that the PF method requires considerable computation time and is limited to simulate very small domains with a few dendrites [3]. On the other hand, stochastic models, such as Monte Carlo (MC) method [4, 5] and Cellular Automaton (CA) method [6-8], have been usually applied to predict the solidification microstructure evolution. Comparing with the MC method, the dendritic growth kinetics is explicitly integrated in the CA model, while the MC model doesn't. Therefore, the CA method is more reasonable to predict the dendritic growth during the solidification of melt."

Jan 1, 2013