Search Documents

By Leslie C. Taylor, Anthony L. Kee

Computational modeling of explosively loaded saturated sand on a suspended plate is an inherently complex problem. In this study a computational method is used to predict the dynamic pressure load/impulse from a buried charge. A fully coupled code modeling the explosive charge, sand and air with a Eulerian mesh and the plate structure with Lagrange elements, such as Gemini-Dyna, is an appropriate computational tool. In water-saturated sand, the loadings are soil and water-transmitted shock and bubble or blast effects depending on the depth and the charge size. The computations provide the capability to determine the total impulse load placed on a horizontal plate on or above the soil due to an explosion in the saturated sand for different burial depths in saturated sand and different plate stand off distances. The following test case will be presented: (a) 4.5 Kg (10 lbs) TNT, 10.16 cm (4 in) depth of burial to top of charge, and 20.32 cm (8 inch) stand off distance (Test 4).

Jan 1, 2005

By C. Tan, C. Detournay

Wellbore stability is a major concern for drilling in naturally fractured formations. The infiltration of drilling mud into the fractures is a factor, which may compromise wellbore integrity due to the loss of strength experienced by the fractures. The reduction in the fracture frictional and cohesive properties is attributed to weakening of the infill materials and lubrication of fracture asperities by the mud filtrate. One way to investigate the influence of mud infiltration on wellbore stability in fractured formations is to use a numerical modeling approach. A requirement for the numerical model is the ability to predict the region of mud infiltration. This paper describes a numerical technique to model the extent of mud infiltration for application in coupled fluid-mechanical simulations in densely fractured rock masses using the computer code, UDEC. Several factors are considered including fracture pattern, in-situ stresses, reduction in joint frictional angle and cohesion, viscosity of the mud filtrate and fracture plugging. Results of numerical simulations are presented to illustrate the technique and show the influence of these factors on the extent of mud infiltration. Un des problèmes importants posés par le forage en milieu fissuré saturé, est celui du maintient de la stabilité du puit. L’infiltration de boue de forage dans le massif peut avoir pour effet de réduire les propriétés frictionnelles et cohésives des fractures et de compromettre ainsi la stabilité du puit de forage. Le radoucissement des propriétés des fractures est attribué à la dégradation du materiau de remplissage et au pouvoir lubrifiant des boues. L’analyse numérique offre un moyen d’investiguer l’influence jouée par l’infiltration de boue sur la stabilité du puit. Un des critères de sélection du modèle numérique sera son pouvoir de prédiction du domaine d’infiltration de boue. Cet article décrit une technique numérique permettant de modéliser l’infiltration de boue de forage dans le cadre de simulations hydro-mécaniques, en milieu à haute densité de fractures, effectuées avec le logiciel de calcul UDEC. Parmis les facteurs tenus en compte, on trouvera la nature du réseau de fractures, les contraintes in-situ, le radoucissement en friction et cohésion des joints, la viscosité de la boue et le bourrage des fractures. Des résultats de simulations sont présentés pour illustrer la technique et montrer l’influence de ces facteurs sur la taille du domaine d’infiltration de boue. Bohrlochstabilität ist ein wichtiger Aspekt für das Bohren in geklüfteten Gebirgsformationen. Das Eindringen der Bohrlochstützflüssigkeit in die Klüfte kann die Bohrlochintegrität wegen des Verlustes der Festigkeit des Kluftsystems gefährden. Die Reduktion der Reibung und kohäsiven Eigenschaften ist der Schwächung der Kluftfüllung und der Schmierwirkung der eingerungenen Stützflüssigkeit in Bezug auf die Kluft-Asperities geschuldet. Numerische Berechnungen sind ein Weg zur Untersuchung des Einflusses der Infiltration der Bohrlochstützflüssigkeit auf die Bohrlochstabilität in geklüfteten Gebirgsformationen. Eine Voraussetzung für numerische Modelle ist die Fähigkeit den Bereich der Infiltration der Stützflüssigkeit zu prognostizieren. Diese Publikation beschreibt ein numerisches Verfahren, um die Reichweite der Infiltration von Stützflüssigkeiten mittels fluid-mechanisch gekoppelter Simulationen im engständig geklüfteten Gebirge mit dem Programm UDEC zu modellieren. Verschiedene Faktoren, wie Kluftmuster, Gebirgsspannungszustand, Reduktion von Kluftreibungswinkel- und kohäsion, Viskosität der Stützflüssigkeit und Rissabdichtung wurden untersucht. Ergebnisse numerischer Simulationen werden präsentiert, um die Vorgehensweise zu illu

Jan 1, 2003

By R. Hogg, A. C. Maffei, D. T. Ray

Based on standard population balance models for agglomeration and size reduction processes, a simplified model has been developed to describe simultaneous growth and breakage in the flocculation of agitated suspensions. The results of semi-batch flocculation tests have been used for estimation of the model parameters and empirical relationships have been obtained to describe the effects of process variables and to establish scale-up criteria. Application of the model to the design, operation and control of continuous flocculation processes is discussed.

Jan 1, 1990

By S. P. Mehrotra, P. C. Kapur

A phase-species classification is suggested for the mathematical models of flotation kinetics. Important models are discussed and the structure of their rate constants and parameters is analyzed. Split or enhancement factors are formulated for different kinetic models. The steady state dynamics of the flotation circuits can always be described in terms of these factors. For simulation of circuits, feed is characterized in a formal manner starting from a multivariate distribution of particles in flotation rate, size and mineral content. Synthesis of optimum circuits which simultaneously searches for the best process variables and structural parameters from a generalized master network is described. Trends in the areas of search algorithms, role of middlings and feed identification etc. are indicated.

Jan 1, 1989

By Akul Joshi, Ming C. Leu, Robert Guttenberg, Susan L. Murray

The aerospace and automotive industries are facing a significant risk for cumulative trauma disorders from high-repetition, long-duration tasks. Additional risk factors such as shocks, vibrations and sustained uncomfortable postures oftentimes contribute to musculoskeletal, neurological injuries associated with the hand-arm system. The power tools used for the operation reduce the comfort and working efficiency of the operators, thus lowering their health and safety and the quality of operation. This paper investigates the fastening operation on the assembly line of a major aerospace company for quantifying hand-arm vibrations with the objective of developing a dynamic model of the hand-arm system. The model will be used to ascertain the effect of the various risk factors on the quality of the operation. Experiment Setup A system of three tri-axial accelerometers is used to collect vibrations entering the hand-arm system and their positioning follows the ISO-5349 standard. Reaction force during the operation is measured using a force sensor. The posture of the hand-arm system is measured using two goniometers and one torsiometer. The pilot study conducted has two subjects and one pistol-grip power tool. The experiment setup and positions of various sensors is depicted in Figure 1. [ ] Analysis Method The analysis consists of transmissibility of vibrations and frequency weighted acceleration based on the recommendations of the ISO-5349 standard [1,2]. The transmissibility [3] was deduced in frequency domain after measuring data in time domain from the

Jan 6, 2006

By Jan K. Wolski

This paper presents a simple method for determining air temperature gradients along a tunnel with time-dependent intake air properties. The variable intake air properties may be the flow direction and rate, air temperature and humidity, and the carbon dioxide content. The transient intake air properties are simulated by a series of steady intake air conditions, using a superposition method. The mathematical model was verified by reduced-scale experiments in the Waldo mine. The method may be usefull in conditions where the computation speed is more important than the precision of the results.

Jan 1, 1995

By Kenneth P. Michalek, James E. Kelly, Jonathan A. Dantzig

"Many investigations have been reported regarding the simulation of continuous casting with emphasis on in-mold heat transfer. These analyses have been limited, however, because they require either estimation or measurement of the characteristic heat transfer between the ingot and the mold. In this paper, a self-consistent model 'is presented for calculating the heat transfer between a solidifying continuous casting and the mold, in the context of the continuous casting of steel. Using both thermal and stress analyses, the air gap formed during continuous casting is calculated and used to compute the heat transfer coefficients. Iteration between the thermal and stress analyses is used to determine a self-consistent solution. IntroductionThe continuous casting process for steel is illustrated schematically in Fig. 1. Molten steel is tapped into a tundish from which it flows at a controlled rate into the continuous casting mold. The initial shell of the casting forms in the mold by transfer of heat from the metal through the mold wall to the cooling water. The remainder of the casting is solidified by application of water sprays, and/or radiation cooling below the mold.Heat transfer in the mold portion of the process is important because it controls the formation of various defects which are deleterious to subsequent processing. Examples of defects formed in this region are surface irregularities, such as cold folding and oscillation marks, and more serious defects, such as longitudinal or transverse cracks. A more catastrophic problem which may occur is a breakout, in which the shell ruptures and molten steel is discharged in an uncontrolled way."

Jan 1, 1986

By Mostafa O. El-Bealy

"To estimate the interdendritic cracking tendency in casting processes, a new model for interdendritic cracking criteria has been proposed. The model combines a concept of the alloying element effect on the interdendritic strain generation and therefore on the susceptibility of interdendritic cracking phenomena. A susceptibility of interdendritic cracking is modeled by using El-Bealy's interdendritic cracking susceptibility function ""BJCSF' and the opening displacement between the dendrites ""EIA "". These concepts are modeled based on the previous experimental observations and ElBealy' approach. The model is implemented numerically for a one-dimensional model problem with some similarities to the casting of steel. The variation of some key parameters that are known to influence the frequency and severity of the interdendritic cracking is then studied.I. INTRODUCTIONInterdendritic cracking ""hot tearing"" commonly encountered in both ferrous and nonferrous castings (1,2]. It is generally believed that interdendritic cracks start to nucleate and develop in the interdendritic strain region within the mushy zone [l-5]. The interdendritic strain region is then definitely coherent, but continuous films of liquid still exit (4,5]. This distinguishes interdendritic cracks from cold cracks forming during the casting process after complete solidification (6,7].Much effort has been put into the understanding of the interdendritic cracking phenomena, and several theories have been proposed. Fellini [8] stated that interdendritic crack ""hot tear"" will result if the material is subjected to too high of accumulated strain within so called vulnerable part of the solidification interval. Kinosita et. al. [9] and Campbell (10] also emphasized the role of tensile stresses in the forming of interdendritic cracks. Most interdendritic cracks criteria simply consider the size of solidification interval, stating that solidification range is associated with a lager interdendritic cracking susceptibility than a short range. Clyne and Davies (11] formulated a more refined such criteria based upon the time spent in different regimes of the solidification interval. They defined a vulnerable region in which thin continuous films of interdendritic liquid exist, and the permeability is low. When thermal strains are induced in this region, the film is not able to sustain the stresses, and an interdendritic crack will form."

Jan 1, 2001

By Gordon Irons, Kevin Graham, Krishnakumar Krishnapisharody

Ladle or secondary metallurgy in steelmaking is an essential process step for the production of high-quality steel. Professor David Robertson was a pioneer in the modeling of these processes, in particular gas-metal fluid dynamic interactions and the rates of slag-metal reactions. This is now a highly active area of research, and some of the recent developments in this area will be reviewed.

Jan 1, 2014

By Vaughan R. Voller

In many high temperature processes involving liquid melts the common practice is to solidify a fraction of the melt adjacent to the containment wall; a practice that can provide protection of the structural components of the process. The position, shape and transient behavior of these so called "skulls" and/or "ledges" are directly dependent on the thermal process conditions. This paper presents a model of a generic ledge process that can predict the steady shape of a ledge. Two numerical implementations are presented and compared: (1) based on a space grid of continuously deforming finite elements (2) based on a conjugate grid of fixed control volumes and deforming boundary elements.

Jan 1, 1996

By R. Hamar, G. Lamanthe

"To realize complex shaped pieces, like turbine blades, superalloys are cast with the lost-wax process. The end properties of these parts depend upon their microstructure as well as upon the occurence of casting defects such as micropores and hot tears. Such features are governed by the local solidification conditions which may be describe by the calculation of thermal transfers between alloy and mould. It needs the evaluation of numerous parameters which are fitted from experimental measurements in real castings. In further simulations, any parameter may be changed, in order to determine its influence. Such a single modification is quite impossible to realize in real castings, which shows that numerical modeling is a very good way to study casting processes.This modeling technique of casting is applied here to the real study of a typical low pressure turbine blade of a civil engine.IntroductionAircraft engine manufacturers use superalloys to produce turbine blades. They have to obtain a given microstructure and to avoid solidification defects to give good mechanical properties to theirs castings. For example, microporosity is induced by a large shrinkage that cannot be completely fed with the remaining liquid alloy .The optimization of casting processes needs a good Knowledge of the solidification process, that must be controlled at every place in the casting, and at every time of its freezing. Material and heat transfers that occur during the solidification, can be estimated from rather complex calculations. To model these physical phenomena, the evaluation of many parameters must be made. In this aim, some casting experiments are fitted with numerous thermal captors. The identification of experimental and computed thermograms let us estimate these parameters.A long-range scientific cooperation between SNECMA and INPG has enabled the realization of such studies on this topic. From all of them, this paper describes an experimental and computed investigation which has been led on a typical low pressure turbine blade of a civil engine."

Jan 1, 1986

By M. K. Alam

Chemical Vapor Deposition (CVD) processes are currently attracting special attention in materials processing. CVD processes are now used to produce special materials in the form of thin films, coatings, composites, and optical fibers. This process produces material of high purity, in complex shapes, and in a variety of microstructures. It is a capital intensive process, involving use of high cost raw materials. Modeling simulations have the potential of reducing the time for CVD . process design, and the need for experimentation. Modeling would, therefore, appear to be an attractive option for producing a product of acceptable quality. At Ohio University, progress has been made in modeling of CVD processes, specially in the areas of optical fiber preform fabrication, and fiber coating.

Jan 1, 1990

By S. Niyomwas, R. G. Reddy

"A mathematical model was developed to describe the plasma gas and particle dynamics and conversion yields. The velocity and temperature fields in the plasma were calculated by solving 1 the Navier-Stokes equations and the thermal energy balance equations. Conservation equations for the system species consider the multi-component diffusion and chemical reactions. Through the application of this model, temperature distribution, flow fields, and production rates were computed. The model was used to study the thermal decomposition of ilmenite in the non-transferred arc plasma reactor. The computed production rates of TiN from ilmenite using nitrogen and methane gas mixtures, are in the range of 5 to 100 x 104 kg/m3/s at 2500 KIntroductionThe widespread interest in thermal plasma processing has stimulated a great deal of research effort into the various aspects of this process, including plasma' spraying, thermal plasma chemical vapor deposition, plasma synthesis of fine powders, thermal plasma decomposition, plasma melting, and plasma extractive metallurgy.In the plasma synthesis of fine powders, the earlier experimental and theoretical studies aimed at predicting the final products1-5. Then it was followed by more recent investigation concerning the profile of plasma gas velocity, plasma temperature and concentration. In this regard, useful plasma processing models have been reported by numerous investigation6 -11 which were increased by both theoretical analysis and measurement.As a result of this work a reasonable plasma-processing model is proposed. It is to be noted that all the thermal analysis that has been made neglected the fluid -flow effects. Also, in the analysis aimed at the quantitative description of the fluid-flow field, the effect of heat transfer was not taken into account."

Jan 1, 2000

By B. Q. Li, S. P. Song

"Numerical models are developed to represent complex electromagnetic, free surface deformation and fluid flow phenomena in floating zone refining processes with radio frequency inductive heating. The computational methodology is based on the coupled finite/boundary element method for the electromagnetic field distribution and the Galerkin finite element method for magnetically driven fluid flow and free surface profiles of the molten zone. With the models, the complex transport and free surface phenomena in a floating zone system can be studied as a function of various operating conditions including applied current, frequency, inductor position and shape, surface tension, floating zone diameter and height. Results show that the equilibrium shape is sensitive to the height of the molten zone and that the free surface shapes can have a significant effect on the melt flow in the molten zone. Our extensive simulations also suggest that to ensure accurate results, the same mesh should be used, whenever possible, for both the electromagnetic forces and fluid flow calculations.1. IntroductionThe floating zone processes have been widely used in the materials industry to grow ultra-pure rare earth metal or semiconductor single crystals for sensor or electronic applications that have a tough demand on purity. Developed based on the principle of magnetic melting/levitation, the processes possess a main advantage that the molten zone is supported in air by the electromagnetic forces generated by the coil (see Figure 1), and therefore is free from contamination from a container. Consequently, the final products made from the floating zone processes have higher inherent purity than those by other techniques [1-4].Complex physical phenomena occur in the floating zone processing of materials. These include free surface deformation, magnetically-driven melt flow, which may be turbulent or mildly turbulent [5] depending on the applied conditions, and solidification. Thermal stresses also develop in the single crystals, which are affected by the flow and temperature field and which can have direct implications to qualities of single crystals [6]. An understanding of the fundamentals governing these physical phenomena can be of critical value in elucidating the physics that controls the properties of the metal or semiconductor crystals being produced, thereby providing a rational basis for process improvement and quality control."

Jan 1, 1997

By K. P. Nishad

The origin of many casting defects can be traced at the mold filling stage. Entrapped air/gas pockets and oxide film formation are notable examples. In the present work, the mold-filling process is simulated using a mathematical model, built on a multi-physics based commercial package PHYSICA. The simulation code has the capability to solve multi-phase problems. The numerical model addresses a complicated dynamics of filling operation including free surface flow and subsequent flow behavior and is based on conservation of mass and momentum. The Free surface is modeled using the scalar equation method (SEM). The model is validated with benchmark experimental results from the literature. The role of turbulence and the pouring rate on variables affecting sites prone to oxide film formation are investigated using a simple but novel technique. Key findings of the work are presented.

Jan 1, 2006

By Sushil Bhandari

A computer model is developed for the problem of dispersal of dust (suspended particulate matter) resulting from blasting operations in surface mines. Undulated surface (lower boundary) simulation was carried through sixth order polynomial by least square approximation. Segmentwise terrain simulation was linearly added to build ground feature extending either side of the blast location. Assuming the initial dust concentration in the vertical column to. be uniform and treating the finite vertical line sources as a series of point sources to height, H varying from zero to top of the dust mushroom, integration was done of the concentration effects from H = 0 toL tan[0] (L = horizontal distance of the blast location and [0] is the angle of the top of dust cloud from the point of observation). Using frame of reference moving with the mean wind, the ground level 'concentration (GLC) has been computed over the undulated terrain. Standard deviation (s) of concentration has been parameterized based on mixed layer similarity theory in convective boundary layer.

Jan 1, 2002

By B. Farouk, D. Apelian

"A numerical model has been developed to predict the temperature history of particles injected in a low pressure d.c. plasma jet. The temperature and velocity fields of the plasma jet are predicted as a free jet by solving the parabolized compressible Navier-Stokes equations using a spatial marching scheme. Particle trojectories and heat transfer characteristics are calculated using the predicted plasma jet temperature and velocity fields. Correction factors have been introduced to take into account non-continuum effects encountered in the low pressure environment. The plasma jet profiles as well as the particle/plasma interactions under different jet pressure ratios (from underexpanded to overexpanded) have been investigated.I. IntroductionIn the plasma spraying process a d.c. plasma jet is used as a heat source to melt and accelerate powder particles which subsequently impinge and solidify on a given substrate. During low pressure plasma deposition (LPPD) prooess, the plasma jet is expanded into a low pressure environment (40-80 Corr) which results in high plasma jet velocities with Mach numbers ranging from 2 to 3. Previous efforts to predict heat transfer and fluid flow in plasmas have made use of incompressible flow formulation employing turbulence models [1-3] for atmospheric plasma spraying. However, these theoretical calculations have been limited to the subsonic flow regime; hence, the compressibility effects and viscous heat dissipation were neglected. In the LPPD process, the plasma jet is operated at the supersonic flow regime, therefore, the compressibility effects and viscous heat dissipation are no longer negligible. This is the main motivation of this work.In addition, production of high quality dense deposits requires that a large fraction of the injected particles be in a molten state when they impact the substrate. An additional motivation for this work is that it is important to understand the plasma/particle heat and momentum transfer since these phenomena dictate both the particle thermal history and particle trajectory.Early computational works in supersonic nozzle field flow [4,5,6] consisted of patching methods that divided the flow field into an invisoid free stream and a viscous boundary and. mixing layers. Each was analyzed independently and coupled through appropriate boundary conditions. Later schemes [7,8] involved solving time-dependent, compressible Navier-Stokes equations over the entire computational domain. These schemes show a great promise in predicting flows with complex structure, and where the viscous effects became prevalent."

Jan 1, 1986

By J. F. Shao, G. Duveau, W. Z. Chen, O. Ozanam

In this paper, we present a sensitivity study for the numerical modelling of hydromechanical coupling around a shaft in partially saturated claystone. Two representative constitutive models have been used: an elastoplastic damage model (IDEPLML) and an elastic perfect plastic model based on the Hoek and Brown failure criterion (H&B). A series of comparative studies have been performed with respect to main parameters involved in coupling process (mechanical, hydraulic and hydric properties). The emphasis is put on the evaluation of induced plastic, damaged and failed zone due to excavation and drying process. Dans ce papier, nous présentons une étude numérique paramétrique du comportement hydromécanique couplé d’une argilite partiellement saturé, lors du creusement et de la maintenance d’un puit L’étude paramétrique porte à la fois sur le type de modèle et sur les principaux paramètres (mécaniques hydrauliques et hydriques) intervenant dans le couplage. L’accent est principalement porté sur l’évaluation des zones plastifiées autour de l’ouvrage lors de la désaturation. In unserem Beitrag geht es um eine digitale parametrische Untersuchung über das hydromechanische Verhalten eines teilweise durchtränkten Tons beim Abteufen und dem Betrieb von Schächten. Die parametrische Analyse konzentriert sich auf den Modelltypus und auf die wesentlichen bei der Verbindung wirksamen Parameter, d.h. die mechanischen, hydraulischen und hydrologischen. Der Schwerpunkt liegt auf der Bewährung der plastifizierten Bereiche des Bauwerkes während der Austrocknung.

Jan 1, 2003

By Murali Thirukkonda, Rahul Shah, Srinivas Tangrila, Shankar Rachakonda, Jan Kajuch, Arthur Gurson

"High-temperature superconductors have potential for a variety of applications such as power generators, superconducting magnets for mine sweepers or ship propulsion motors, and magnetic levitation transportation systems. Interest in the use of superconducting wires and tapes increased dramatically after the discovery of high-critical-temperature (high-Tc) oxide superconducting compounds, which offer operating temperatures above that of liquid nitrogen. Several techniques have been proposed to produce commercially viable high-Tc oxide based superconducting wires. One of the promising methods involves multi-pass powder-in-tube (PIT) wire drawing followed by rolling and heat treatment. This work focuses on the development of a numerical model to simulate the PIT drawing process for fabrication of long lengths of silver sheathed Bi-2212 superconducting wires. The finite element based model was used to predict densification of the oxide powder, wire drawing forces, and silver sheath thickness during drawing. A cap-type pressure dependent constitutive equation was implemented in the model to simulate the powder behavior. The model incorporated experimentally obtained material data for the silver and powder. Results indicate that powder densification is increased by decreasing the initial silver thickness and increasing the silver strength. Data from wire drawing experiments were used to verify model predictions. IntroductionSuperconducting wires and tapes can be used to carry large dc currents with no measurable loss due to resistivity. When used in superconducting magnets for power generators, motors, levitation devices, magnetic detectors, or energy storage devices, they could potentially result in significantly reduced energy losses. In the past, most prototypes were tested in liquid helium at 4.2K [1]. However, the necessity to provide liquid-helium cooling systems proved to be a serious limitation in commercialization. In 1986 high-critical-temperature (high-Tc) oxide superconductors were discovered [2], which could allow for higher operating temperatures above liquid nitrogen (77 K). Since then, both academia and industry have concentrated their efforts on the development of high-Tc oxide superconducting wires, tapes and coils. As yet, large scale fabrication of long lengths of these wires with the required current carrying capacity and mechanical properties remains an enormous challenge."

Jan 1, 1994

By E. Kozeschnik

A new model is presented for the simulation of the evolution of precipitates in a multi-component multi-particle multi-phase environment The mathematical formulation as well as the physical background of the model are briefly introduced and its capabilities are discussed on the example of a precipitation sequence in a complex steel. Keywords: Modeling, precipitation kinetics, nucleation, growth, coarsening

Jan 1, 2003