Search Documents

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 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 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

By Min Chen

A numerical simulation of the stress field of the solidified shell in the mold during continuous casting of peritectic steel slab was calculated with commercial software. The results showed that the maximal stress on wide face of the solidified shell was located at the corner and near middle of the slab, and it increased with increasing the width of slab. In addition, the stress was also influenced by operating parameters such as superheat degree, drawing speed and cooling intensity of the mold, and the maximal stress value and its position changed with these parameters. For the slab with section size of 3200mm×150mm, the proper drawing speed was around 1.2m/min, with cooling intensity of 5500L/min while the superheat degree was 15°C to 25°C in order to prevent longitudinal crack occurring during peritectic steel continuous casting.

Jan 1, 2009

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 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 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 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 B. Pesic

Northwest Alloys, a wholly owned subsidiary of Alcoa, produces magnesium from dolomite by reduction with ferro silicon and aluminum. This process produces alkaline by-products in quantities that allow their consideration for use as commercial by-products. Two such products are dicalcium silicate slag, Di-Cal, and magnesium oxide condensate, MgO. The alkaline properties of these materials offer acid neutralizing capability comparable to the traditionally used lime, therefore, warranting a study for their use in applications where the neutralization of acidic waters or process solutions is required, such as treatment of mine waters. In this ongoing study, three mine waters (Bunker Hill Mine, Berkeley Pit, and Blackbird Mine) were selected for treatment with Di-Cal and MgO. Lime was used as the background control. The solid precipitates were subsequently contacted with pure water at four different pH conditions for stability testing of prepared metal precipitates (expanded TCLP). It was found that neutralization with Di-Cal results in precipitates which are of equivalent to superior stability than those from lime treatment.

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 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 Kwang-Yeom Kim, Dong-Gyou Kim, Hyo-Kwan Yoon

This study was intended to simulate the blasting model using notched-hole with numerical approach. Numerical analysis approach was applied to compare and evaluate the blasting model using notched-hole and circular hole. To implement non-linear dynamic blasting model, blasting analysis was carried out with AUTODYN code. Behavior of two models was compared and evaluated using AUTODYN and visual analysis and quantitative evaluation based on pressure and stress to verify the effect to induce the destruction territory to the specific direction when blasting with notched hole was attempted.

Jan 1, 2010

By Jin Zhao, Weijie Wei, Jinwang Zhang, Xiaohang Wan, Shengli Yang

Reasonable mining-retention ratio (ratio of room width to pillar width) is vital to recovery rate and stope stability in underground mining of solid potassium salt deposits. Based on the engineering background of a solid underground potassium salt mine in Khammouan Province, characteristics of the plastic zone of the pillar, and the variation of recovery rate and stope stability under different mining-retention ratio conditions are discussed. The results show that compared with current room width (8 m) and pillar width (10 m), the optimization method of improvement of stope stability. The research results are helpful to recover as much solid potash resources as possible under the premise of ensuring the stability of stope.

Jun 25, 2024