Search Documents

By Sang-Joon Kim, Hae-Geon Lee, Hyun-Jin Cho

"Decarburization phenomena on the surface of liquid iron have been simulated by numerical method and the results were compared to the experimental data. During the DeC process in a converter, a large amount of splashed droplets are generated from the top surface of the bath and consequently the reaction surface area is increasing tremendously. Since the decarburization mainly occurs at the interface, it is important to understand the reactions on the liquid iron droplets. The flow and species concentration fields around liquid droplets taking into account the reaction of carbon with oxygen at the interface were calculated by a commercial CFD code FLUENT which solving multiphase Navier-Stokes equation with Volume Of Fluid method, species transport equations and energy transport equation. The initial carbon concentration was fixed to 4 wt% and the concentrations of oxygen in the atmosphere were varied from 5 to 50 wt%. As comparing the numerical calculation with experimental data written in the reference, the results have a good agreement for decarburization rate.IntroductionIn a primary steelmaking process, one of the main functions is to decarburize from hot metal and the 4% of carbon is decreased to about 0.5% of carbon. Even though there are some possible reactions leading to the removal of carbon from hot metal, the main reaction occurs on the free surface by direct oxidation of carbon at the gas-liquid interface with producing the carbon monoxide gas. In many metallurgical processes involving a jet impinging onto a steel bath surface, especially in a Basic Oxygen Furnace (BOF) converter, it is important to understand the behavior of molten iron droplets splashing out of the steel. The iron droplets are generated by an oxygen or air jet blowing and it increases the total surface area of liquid iron tremendously. Since the decarburization reaction is mainly affected by the surface area, it is necessary to verify the decarburization phenomena around the splashed droplets from liquid iron in order to allow an analysis ofBOF process possible."

Jan 1, 2011

By K. Krishnapisharody

A mathematical model has been developed to predict the interface dynamics and droplet generation of the buoyant phase in a two-fluid system. The model employs the well-known SIMPLE algorithm for obtaining flow fields and the Volume of Fluid (VOF) method, for interface tracking. The implementation of the model in water-oil systems, where the oil layer is placed above a recirculating water bath, is discussed. The evolution of the water-oil interface leading to the formation of oil droplets is computed realistically.

Jan 1, 2008

By Guo Chen, Yifeng Shi, Xiaobiao Shang, Nanshang Shen, Bangqi Zhang, Jinhui Peng, Junruo Chen

"The energy is introduced directly into the volume of the honeycomb ceramic regenerator due to the complex interaction of microwaves with minerals in microwave gas heating system. A three-dimensional numerical model is developed to simulate the electromagnetic field by solving Maxwell's equations using finite-difference time-domain method. Computational result is presented for the microwave power absorption of the regenerator. Effects of the dielectric constant and the effective loss factor on the microwave power absorption in the regenerator are studied as well. Numerical results show that the heating pattern strongly depends on the dielectric property of the regenerator. The demonstration of the model of electromagnetic field distribution can be applied effectively and efficiently to optimize the fundamental parameters in the microwave heating processing.IntroductionIn recent years, microwave energy has been utilized in many fields, such as food processing, chemical engineering and mineral processing [l]. In contrast to conventional heating, where heat is transferred from a surface to the interior, microwave energy is introduced directly into the volume of minerals due to the complex interaction of the microwave with the mineral. Based on volumetric heating effects, microwave heating has several advantages such as pollution free, rapid heating technique, internal heating, high efficiency over conventional heating [ 1-4]."

Jan 1, 2013

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

Based on the commercial computational fluid dynamics software, Fluent, the fluid flow and temperature field in an iron-bath reactor were simulated numerically, and the effects of top-blown and side-blown operating parameters on the mixing time of iron bath, distribution range of oxygen in top space, and post-combustion ratio were investigated. The simulated results showed that the mixing time of iron bath increased initially and then decreased with the increase of inclined degree and immersion depth of side-blown lance. The distribution range of oxygen became wider as pulling the top lance away from the iron bath. Furthermore, post-combustion ratio was reduced given the top lance was too close the bath. The simulation results provide theoretical base for the process optimization of iron-bath reactor and further development of new type iron-bath smelting reduction process.

Jan 1, 2009

By Bryan J. Travis, Thomas L. Cook

INTRODUCTION Our growing concern for adequate and secure sources of energy and minerals has stimulated vigorous exploration for new sources, research toward a better understanding of geological processes, and development of new extraction technologies. The need for control, or at least prediction, of subsurface fluid flow is important for many of these technologies for example: primary, secondary and tertiary recovery of oil; ground water and waste management; in-situ fossil energy extraction (oil shale, coal, tar sands); and solution mining of uranium, copper and other minerals. These technologies, especially the last two, are characterized by highly complex systems. A partial list of the physical processes occurring would include flow in porous/fractured media, multi-phase and multi-component flow with heat transport, chemically active fluids and soils, tracers, diffusion and dispersion, fracturing and dissolution. A great deal of understanding of how these processes behave and interact can be obtained from models. MODELS Theoretical models are valuable because they: 1) provide a frame of reference for interpreting results of laboratory and field experiments; 2) once validated by experiments, allow a variety of geometries, injection/production strategies, etc., to be examined (at relatively low cost) for efficiency and stability; 3) can provide guidance to the design of experiments and field operations. Most models are based on the fundamental principles of mass, momentum and energy balance. But from this starting point, many paths can be taken. For example, there is the theory of Payatakes (1973) which concentrates on the microscale dynamics. In this approach, the rock or soil matrix is represented by a complex of characteristic channels (such as periodically constricted smooth tubes). Detail of flow within the characteristic channel is calculated very accurately. A difficulty with this model is that description of a representative channel can require several parameters which may not be easily measured. Also, it is not clear how the channels can be combined practically to model a large scale flow. Another type of model is the "global" one described by Bear (1972). Here, the continuum equations for conservation of mass, momentum and energy are averaged over a distribution of pore channels, resulting in a set of conservation equations in which the small scale structure of the medium is replaced by quantities such as porosity, permeability, dispersion coefficients and tortuosity. This approach allows computation of large scale flows. However, additional constitutive equations are needed which relate averaged quantities such as permeability to observables such as local saturation, particle size distributions, and others. An important difference between models is the way they handle the momentum equation. Payatakes' model solves the full equation. In others, it is replaced by a simpler relation such as Darcy's law (valid for slow flow rates) or by Forchheimer's equation which extends Darcy's law to higher flow rates. These simple relations can nevertheless match a great deal of experimental data (Dullien 1975). The permeability term which appears in Darcy's law and in Forchheimer's expression has been related to other quantities such as porosity and particle surface area. The Kozeny-Carmen equation is a well-known example, valid for some ranges of porosity and particle sizes and for some materials. Several other semi-theoretical, semi-empirical formulations have been devised, none of which are entirely satisfactory. One goal of researchers has been the ability to predict permeability of a porous material from basic measurable quantities such as grain size distributions without recourse to adjustable parameters. This effort has proceeded from consideration of distributions of idealized, non-intersecting channels, to intersecting channels, to consideration of both distributed pore radii and pore neck radii. This last approach has been used successfully to predict permeabilities (ranging over several orders of magnitude) in sandstones (Dullien 1975). Additionally, studies on explicit networks of channels (e.g., Fatt 1956) and percolation theory (e.g., Larson et. al. 1981) have been used to examine interconnectivity effects in porous media with the hope of eventually being able to predict permeability accurately. To be of more than theoretical interest, a model must be compared s controlled laboratory conditions where boundary and initial conditions and relevant material properties can be accurately determined. Also, extraneous forces can be eliminated so that the interactions between processes of interest can be clearly seen and flow models can be vigorously tested. In contrast to the well-defined environment and

Jan 1, 1982

By T. W. Doe, J. Noorishad

The problem of fluid flow from a well into a horizontal fracture has been studied using a nonlinear coupled hydromechanical finite element model. Such an approach is required due to the strong nonlinearity and anisotropy of deformation and fluid flow moduli introduced by the presence of the fracture. The resulting nonuniform total stress field invalidates the application of the conventional fluid flow techniques to the problem. Stable convergent solution from the above analysis are verified by a coupled steady state algorithm.

Jan 1, 1982

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 T. F. Adams, L. G. Margolin

The Bedded Crack Model (BCM) is a constitutive model for brittle materials. It is based on effective modulus theory and makes use of a generalized Griffith criterion for crack growth. It is used in a solid dynamic computer code to simulate stress wave propagation and fracture in rock. A general description of the model is given and then the theoretical basis for it is presented. Some effects of finite cell size in numerical simulations are discussed. The use of the BCM is illustrated in simulations of explosive fracture of oil shale. There is generally good agreement between the calculations and data from field experiments.

Jan 1, 1982

By C H. Ryu

This paper presents a simplified two-dimensional numerical model of block fragmentation in a Jointed rock mass during the late or throw stage of a blast after wave effects have subsided. The model is based on the distinct element method. Examples illustrate the approach and allow for a comparison of effects of Joint patterns on fragmentation and resulting block motion.

Jan 1, 1986

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 K. S. Rao, R. P. Tiwari, V. K. Shrivastava

ABSRACT: In this paper numerical simulation of jointed block mass behaviour obtained in an experimental study has been presented. Singh et al. (2002) conducted extensive tests on the jointed model materials with several joint configurations. The results of this experimental study were used to validate the strength behaviour of numerical block mass model developed using professionally available software. This software basically contains different kinds of material models. From the results of numerical simulation the existence of strength anisotropy was observed. The strength anisotropy so obtained has been compared with the experimental values reported by Singh et al. (2002). The anisotropic curves in both cases are U-shaped with flat base. They show good agreement in sliding mode of failure zone in which inclination ß (= 90° - ?) of joint set-I lie between 20° to 40° for all cases of interlocking. ZUSAMMENFASSUNG: In dieser numerischen Papiersimulation des verbundenen Blockmasse Verhaltens, das in einer experimentellen Studie erhalten wird, ist dargestellt worden. Singh et al. (2002) leiteten umfangreiche Tests auf den verbundenen vorbildlichen Materialien mit einigen gemeinsamen Konfigurationen. Die Resultate dieser experimentellen Studie wurden verwendet, um das Stärke Verhalten des numerischen Blockmasse Modells zu validieren, das mit professionell vorhandener Software entwickelt wurde. Diese Software enthält im Allgemeinen unterschiedliche Arten der materiellen Modelle. Von den Resultaten der numerischen Simulation wurde das Bestehen von Stärke Anisotrophie beobachtet. Die Stärke Anisotrophie also erreicht ist mit den experimentellen Werten verglichen worden, die von Singh et al. (2002) berichtet werden. Die anisotropen Kurven in beiden Fällen sind mit flacher Unterseite u-förmig. Sie zeigen gute Vereinbarung in gleitender Erscheinungsform des Ausfalls Zone, in der Neigung ß (= 90° - ?) von der Verbindung stellenSie Lüge zwischen 20° auf 40° für alle Fälle vom Ineinandergreifen ein. RESUME: Dans cette simulation numérique de papier du comportement joint de la masse de bloc obtenu en étude expérimentale a été présenté. Singh et autres (2002) ont effectué les essais étendus sur les matériaux modèles joints avec plusieurs configurations communes. Les résultats de cette étude expérimentale ont été employés pour valider le comportement de résistance du modèle numérique de la masse de bloc développé en utilisant le logiciel professionnellement disponible. Ce logiciel contient fondamentalement différents genres de modèles matériels. Des résultats de la simulation numérique on a observé l'existence de l'anisotropie de force. L'anisotropie de force ainsi obtenu a été comparée aux valeurs expérimentales rapportées par Singh et autres (2002). Les courbes anisotropes dans les deux cas sont en U avec la base plate. Elles montrent la bonne concordance dans le mode coulissant de la zone d'échec en lequel la inclination ß (= 90° - ?) du joint placez-Je le mensonge entre 20° à 40° pour tous les cas d'enclencher.

Jan 1, 2003

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

"An advanced casting simulation-based optimization approach was applied in this study to assist in the improvement of the mold design of aerospace components made of A3 56 alloy. By using this approach, mold filling and solidification related defects (including macro-shrinkage and shrinkage porosity) were significantly minimized and hence it helped in cost reduction, performance enhancement and quality assurance of complex A3 56 cast parts.An experimental validation and detailed calibration procedures of the models for prediction of macro-shrinkage and shrinkage porosities were performed using A3 5 6 plates cast in fiiran-silica sand molds using the Prometal's Rapid Casting Technology (RCT) - mold printing technology. Correlations between Niyama values and the pore percentage were also developed. Thus, the severity level of shrinkage porosity can be determined via the Niyama criterion. Predictions were then compared with the macro-shrinkage and porosity measurements in plates of various plate thicknesses and in other commercial A3 56 casting components.IntroductionThis paper describes the process modeling approach used in the development of the mold design technology used to produce assist in the improvement of the mold design of aerospace components made of A356 alloy. Nova Flow&Solid™ software [3] was applied in this study. This casting simulation software was previously validated for several alloys including Ti, steel and superalloys [1-4]. The objective of this paper is to calibrate and validate the simulation software for A356 alloy sand mold castings. For this purpose, an experimental validation and detailed calibration procedures of the models for prediction of macro-shrinkage and microshrinkage (e.g., shrinkage porosities) were performed using A356 plates cast in fiiran-silica sand molds using the Prometal Rapid Casting Technology (RCT) mold printing technology. The details of the RCT technology are shown in Fig. 1 [5]. Prometal's RCT technology is used for rapid production of sand molds and cores directly from CAD files. The RCT technology relies upon ink jet printer technology to selectively bind the grains of sand to create three-dimensional sand molds and cores. Mold filling and solidification modeling as well as stress analysis are typically performed to properly design the casting and the rig system of each prototype component."

Jan 1, 2012

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