Search Documents

By Robert B. Johnson

The refraction seismic method has been used as an indirect means of exploring for fluorspar in southern Illinois by the Illinois State Geological Survey. Use of the method has been restricted to those areas where mineralization is associated with faults, a restriction necessitated by the operating principles of the refraction method. The seismic data enable the interpreter to locate fault zones in the subsurface and, under ideal circumstances, to make estimates of dip direction and magnitude and vertical displacement provided the local stratigraphic units are known and are geophysically recognizable. In this instance, refraction seismic exploration has been successful in the surface mapping of well-developed faults, both simple and complex with which fluorspar may or may not be associated. However, less clear results have been obtained where faulting has been less severe and where there are numerous accompanying mineralized joints.

Jan 1, 1964

By Olav Mellgren

Part I. The interaction between potassium ethyl xanthate and lead salts has been studied thermo-chemically. It is shown that ethyl xanthate reacts with lead carbonate, basic carbonate, thiosulfate and sulfate according to principles of simple chemistry, whereas reactions with lead sulfite, metasulfite, basic thiosulfate and hydroxide at high pH are associated with secondary reactions. A new compound "xanthogen thiosulfate" is postulated. Part 11. Adsorption of ethyl xanthate on galena has been studied quantitatively and thennochemi-cally. It is shown that an ion-exchange mechanism governs the adsorption process and that the heat of adsorption of xanthate on galena is equivalent to the heat of reaction of ethyl xanthate with bulk lead salts. It is concluded that the entropy content of the surface "compound" on galena is different from that of a corresponding bulk lead salt. Alkali depression is in part shown to be due to the instability of the xanthate-galena system at high pH and to the formation of basic lead salts. Sulfate and carbonate ions are shown to have a stabilizing effect on the xanthate-galena system. PARTI. THERMOCHEMISTRY OF ETHYL XANTHATE AND OF ITS REACTION WITH LEAD SALTS In the past little attention has been paid to the chemistry of lead xanthate. It is known that xan-thates form slightly soluble compounds with lead and that suspensions of lead xanthates in alkaline solution decompose rather readily.' The solubility of lead ethyl xanthate is of the order of 2.6 x 10-6 moles per 1. At increasing pH the xanthate concentration of a saturated solution contains more xanthate. Thus at pH 11 this xanthate content is about 100-fold larger than at pH 7.1 Xanthates are known to be consumed when added to a suspension of slightly soluble lead salts such as lead sulfate and lead carbonate. It is unknown to what extent this xanthate consumption is due to a simple chemical reaction because, frequently, observations are made which reveal secondary reaction products in such reactions. Thermochemical study of the reaction between xanthate and lead salts offers a means to clarify the significance of these reactions as the enthalpy changes can be accurately predicted from heat of formation of lead salts in pure solution reactions. In the present work the heats of formation of lead ethyl xanthate from potassium ethyl xanthate and lead salts are determined and some thermochemical aspects of xanthate in solution are studied. The results of these investigations are presented in Part I. In Part II the heat of adsorption of ethyl xanthate on galena surfaces is presented. All enthalpy changes were measured directly in a microcalorimeter. Calorimeter: The calorimeter used for measuring enthalpy changes was similar to that used by Redinha and Kitchener2 which was based on a description by Derbyshire and Marshall.3 Other useful references4,5 are listed. A sketch of the calorimeter

Jan 1, 1967

By P. K. Leung, S. R. Allada, P. M. Dranchuk, D. Quon

The alternating direction explicit procedure (ADEP) makes use of the boundary conditions to reduce multi-dimensional problems to a series of one-dimensional problems. The method, previously applied to reservoirs containing only an under-saturated oil, has now been extended to cover the case of natural gas reservoirs. l his involves solving a non-linear partial differential equation, application of the procedure is straightforward and no calcuIational problems were encountered. INTRODUCTION There has been a growing interest in formulating mathematical models of petroleum and natural gas reservoirs — models which permit the engineer to examine and evaluate the physical and economic consequences of various alternative production policies. The tremendous reduction in the cost of solving such models in recent years has made possible their use as an almost routine management tool. This reduction is the result not only of improved computer hardware but also of the development of more efficient mathematical techniques. The present paper is concerned with the application of a recently proposed numerical method (ADEP) to two-dimensional gas reservoirs. STATEMENT OF THE PROBLEM Given a two-dimensional Region R, bounded by a closed Curve C (Fig. 1) such that the behavior on Curve C is known, the differential mass balance for each fluid phase in R, neglecting gravitational effects and assuming Darcy's law for fluid transport, can be written as:

Jan 1, 1967

By Charles Pack

Extensive progress has been made in the metallurgy of alloys for die castings. Enthusiastic proponents of some alloys are inclined to make extravagant claims for their materials, which may be justified by physical tests on test bars made from the respective alloys but are seldom justified by the results obtained with commercial die castings made from those alloys. Chemical analysis of the test bars and castings fails to determine the cause of these variations, since they are due to other physical factors entering into the commercial production of die castings. In this paper, an attempt will be made to outline some of the factors in commercial die-casting practice to which these variable results may be attributed, in the hope that the manufacturer of die castings, recognizing the importance of these factors, will devote more attention to them in the design and construction of dies and equipment. In 1914,l the author defined die castings as "metal castings made by forcing molten metal, under pressure, into metallic molds or dies." This definition was given in order to differentiate die castings from other molded metal products such as "permanent-mold" castings, L'slush" castings and similar products that are sometimes erroneously termed die castings. The term "die casting" could be interpreted to include all castings made in metal molds or dies; but since this term, like most of our technical nomenclature, has been developed by commercial practice we must look to this practice for its proper definition. The Germans, with their usual devotion to desecriptive names, have adopted the name of sprilzguss (squirt cast) for die castings. Undoubtedly this name is more descriptive of the die-casting process, where the metal is literally squirted into the die and serves to differentiate this process from other metal mold processes, where the metal is poured into the die. In the earlier paper mentioned, the author also emphasized the fact that the definition given for die castings indicated clearly that the satisfactory operation of the process must comprise three factors, all of which must be considered of vital importance to the successful operation of the die-casting process:

Jan 1, 1931

By Richard Rice

TURBO blowers for blast-furnace blowing have now been in use for some years, and a review of the experience gained and the present state of progress may be interesting. FIG. 1.-STEAM TURBINE DIRECT CONNECTED TO AIR COMPRESSOR. The first blast-furnace blowing engine of this type to be installed in America was built by the General Electric Co. and is installed at the oxford Furnace of the Empire Iron & Steel Co. This machine was put into operation March 8, 191.0. A duplicate unit is situated at the furnace of the Northern Iron Co. at Port Henry, N. Y., which unit was put into opera ion June 30, 1910 (Fig. 1).

Jan 5, 1914

By A. B. Kinzel

Since the exposition of the behavior of certain iron alloys by Sykesl involving the existence of an austenite loop and the discovery of such a loop in the chrome-iron system by Bain,2 there has been much conjecture both on his part and on the part of metallurgists in general as to the exact shape and location of this loop. In previous work commercial steel has been used and the microscope has been employed to locate the critical points so that the position of the loop has been quite approximate. Numerous attempts have been made to measure the critical points of the carbon-free alloys but all of these for one reason or another have failed. A satisfactory method for studying these phenomena is presented here. MateRial The first difficulty in this study was to obtain pure chromium-iron alloys. This has been overcome by melting in a vacuum electrolytic iron and chromium analyzing 99.84 Cr, Fe and C. A pure magnesite crucible was used and the melt allowed to solidify in the crucible. A typical analysis is chromium, 10.10 per cent., carbon, 0.006 per cent., silicon and manganese nil. To overcome the problem presented by segregation, ingots were made into wire 0.04 in. in diameter. This was accomplished by hand-forging the ingots into rods, using an oxyhydrogen flame for heating, in order to avoid carbon. Subsequently the rods were swaged to wire, the heating being carried out in hydrogen. Thus local segregation was eliminated although differences in analysis between portions 50 ft. apart were noted. An electric-resistance method was first tried to determine the critical points, but due to the physical defects in the wire this failed. Bain had previously shown the presence or absence of critical points in a qualitative way by observing the sag in a catenary. The author has devised the following method which gives an accurate determination of the critical points from expansion-temperature curves of suspended wire, a special telescopic dilatometer being used for the purpose.

Jan 1, 1962

By M. L. Fuller

MOST commercial alloys undergo changes in phase composition after casting. This is a natural result of the fact that the alloys are not in a state of phase equilibrium as cast and phase changes will therefore take place, with more or less rapidity, in an attempt to attain equilibrium. The practical aspect of aging in zinc-alloy die castings concerns the changes in physical properties and in dimensions that may accompany aging phase changes. Harmful alterations have been reduced to a minimum and zinc-alloy die castings are being produced today that are so stable with respect to dimensions and physical properties that they are entirely suitable for most practical purposes1. In the first paper2 of this series, and in this paper, aging phase changes of zinc alloy die castings have been investigated and the possible influence of these changes on dimensional changes of the castings has been evalu-ated and compared with the dimensional changes actually experienced. The change of dimensions has been selected for practical comparison in these papers because the effect of phase changes on the dimensions of a casting is directly calculable; from X-ray crystal-analysis data of the phase structure. The die-casting alloys considered contain, in addition to zinc, 4 per cent of aluminum, about 0:03 per cent magnesium, and from 0 to 3 per cent of copper. It was shown in the former paper that the eutectoidal decomposition of the beta aluminum-zinc phase, while intrinsically capable of causing a considerable shrinkage, can have very little effect on the observed shrinkage of die castings because it takes place within a few days of aging, whereas the observed shrinkage continues for five weeks. It was therefore concluded that the other existing phase change,

Jan 1, 1935

By R. Schuhmann

Jan 1, 1965

By K. Leutwyler

The key to realistic casing stress analysis in thermal recovery installations is accurate knowledge of the temperatures involved. Much information leading to prediction of heat losses from tubing string to wellbore has been published. This .same information permits calculation of casing temperature under certain conditions. This paper reviews the basis for most of the published work, the line source solution to the diflusivity equation and examines some of its application criteria in the light of unsteady-state conditions. An alternate finite differences method is presented for correlation, and experimental data pertaining to the heat transfer mechanism in the annulus is reported. In the light of the information presented in the paper, it appears that the assumption of a line source for the tubing-wellbore system is valid even for reasonably short periods of heat flow, .such as "steam soak" operations. The empirically modified equation presented earlier matches computer derived results of the finite differences approach lor periods as short as 50 hours. Previously predicted cas ing temperatures and the associated stress levels are thus verified based on these correlative computations and the reported experimental information. Fall out from the computer output resulting in some studies of the influence of insulating cement properties on casing temperature, and experimentally studied effects of nitrogen placed in the annulus at various pressures., are described. INTRODUCTION Generally speaking, predictions of average casing temperatures have been based on an idealized model involving a centralized tubing string at uniform constant temperature radiating energy towards the casing under steady-state conditions.'. ' Heat is transferred away from the casing to the formation by unsteady-state conduction across thermal barriers consisting of cement and the mud sheath. Heat balances written on the casing permit prediction of instantaneous casing temperatures. During early periods of injection, casing temperature varies appreciably as a function of time; however, for long injection times this variation decays as heat flow from the wellbore to the formation stabilizes, and it becomes virtually negligible as conditions surrounding the wellbore become quasi-steady-state (Fig. 1). Carslaw and Jaeger -escribe a solution to the diffusivity equation for a line source (Fig. 2). A constant supply of heat is introduced into the medium from the line placed into the origin of the x-y coordinate system and parallel to the z-axis. Heat spreads outward in the infinite cylindrical solid and the temperature increase on a line going through point (x, y) parallel to the z-axis at a distance r can be predicted by: The following criterion must be considered: the Carslaw and Jaeger equation is obtained on the basis of a large

Jan 1, 1967

By Scott, Turner

MY whole life has been spent in the mining business, PO I naturally tend to address my remarks particularly to the newly-graduated mining and metallurgical engineers among you. To a certain extent, all forms of engineering principles are involved in the training of the mining engineer, and most sorts of engineering are certainly utilized by him, so every one of this graduatins class is concerned with the field I propose to discuss.

Jan 1, 1932

By E. T. Turkdogan

Jan 1, 1962

By Jasper L. Stuckey

PYROPHYLLITE, first identified as soapstone,1 later as agalmatolite,2 and finally as pyrophyllite,3 has been known to occur in North Carolina for more than 130 years and has been produced intermittently in commercial amounts for the past 100 years. Pyrophyllite mining, however, has become a major industry in North Carolina only during the past 25 years. North Carolina and California are the chief producing states. Current annual production averages 125,000 tons, of which slightly more than 90 pct comes from North Carolina. Distribution and Geologic Setting: All known pyrophyllite deposits in North Carolina are associated with rocks of the Carolina Slate Belt,4 which is composed of a vast series of metamorphosed volcanic and sedimentary materials, occurring as two zones extending in a northeast-southwest direction across the state. The larger of these zones, lying in the central Piedmont region, varies from 8 to 50 miles wide. The smaller zone lies along the eastern edge of the Piedmont region and is partly covered by coastal plain sediments.

Jan 1, 1958

By R. L. Wilcox, M. L. Fuller

MOST commercial alloys undergo changes in phase composition after casting. This is a natural result of the fact that the alloys are not in a state of phase equilibrium as cast and phase changes will therefore take place, with more or less rapidity, in an attempt to attain equilibrium. The practical aspect of aging in zinc-alloy die castings concerns the changes in physical properties and in dimensions that may accompany aging phase changes. Harmful alterations have been reduced to a minimum and zinc-alloy die castings are being produced today that are so stable with respect to dimensions and physical properties that they are entirely suitable for most practical purposes1. In the first paper2 of this series, and in this paper, aging phase changes of zinc alloy die castings have been investigated and the possible influence of these changes on dimensional changes of the castings has been evaluated and compared with the dimensional changes actually experienced. The change of dimensions has been selected for practical comparison in these papers because the effect of phase changes on the dimensions of a casting is directly calculable from X-ray crystal-analysis data of the phase structure. The die-casting alloys considered contain, in addition to zinc, 4 per cent of aluminum, about 0.03 per cent magnesium, and from 0 to 3 per cent of copper. It was shown in the former paper that the eutectoidal decomposition of the beta aluminum-zinc phase, while intrinsically capable of causing a considerable shrinkage, can have very little effect on the observed shrinkage of die castings because it takes place within a few days of aging, whereas the observed shrinkage continues for five weeks. It was therefore concluded that the other existing phase change,

Jan 1, 1935

By M. L. Fuller

MOST commercial alloys undergo changes in phase composition after casting. This is a natural result of the fact that the alloys are not in a state of phase equilibrium as cast and phase changes will therefore take place, with more or less rapidity, in an attempt to attain equilibrium. The practical aspect of aging in zinc-alloy die castings concerns the changes in physical properties and in dimensions that may accompany aging phase changes. Harmful alterations have been reduced to a minimum and zinc-alloy die castings are being produced today that are so stable with respect to dimensions and physical properties that they are entirely suitable for most practical purposes1. In the first paper' of this series, and in this paper, aging phase changes of zinc alloy die castings have been investigated and the possible influence of these changes on dimensional changes of the castings has been evalu-ated and compared with the dimensional changes actually experienced. The change of dimensions has been selected for practical comparison in these papers because the effect of phase changes on the dimensions of a casting is directly calculable; from X-ray crystal-analysis data of the phase structure. The die-casting alloys considered contain, in addition to zinc, 4 per cent of aluminum, about 0.03 per cent magnesium, and from 0 to 3 per cent of copper. It was shown in the former paper that the eutectoidal decomposition of the beta aluminum-zinc phase, while intrinsically capable of causing a considerable shrinkage, can have very little effect on the observed shrinkage of die castings because it takes place within a few days of aging, whereas the observed shrinkage continues for five weeks. It was therefore concluded that the other existing phase change,

Jan 1, 1936

By W. J. Reagan

In a study of the life of ingot molds, it is essential to eliminate all of the variables. In the commercial manufacture of steel this is almost an impossibility. In this study many of the variables have been eliminated, and because of the unusually consistent furnace practice, with small variations in the type of steel manufactured, and a similar type of mold design for all ingot sizes, it is possible to determine with some degree of accuracy the cause of variations in ingot-mold life. All of the ingot molds in this study were used in the production of basic open-hearth killed steel with chemical analyses within the following range: C, 0.50-0.85 per cent; Mn, 0.60-0.80; Si, 0.20-0.30; S, 0.04 max.; P, 0.04 max. All ingot molds were of the closed-bottom type and used for bottom casting only. None of the molds in the study failed from cutting of the sides due to top casting; in fact, on only rare occasions are any molds used for top casting. The general type of this mold is shown in Fig. 1. All sizes have rather a heavy band of metal, about 6 in. wide, at the top to give mechanical strength at this point, and all have cast lugs. Cross sections of the various sizes from top and bottom sections are shown in Fig. 2. Ratios of wall thickness to ingot weight, ratio of ingot diameter to ingot length, and other figures of interest are shown in Table 1. In the majority of cases, the cause of failure in this type of mold is due to heavy fire cracking on the "flats" at a point from 3 to 6 in. from the bottom radius of the mold. Fig. 3 pictures this type of failure. The start of failure at these points usually begins with a very small section of the interior surface becoming involved. As a rule, a section of 1 or 2 sq. in. is the first noticeable indication that mold failure is taking place. Gradually the small section increases until a sizable cavity develops, causing a protuberance on the bottom of the ingot, which eventually becomes too large to be stripped. The mold is then removed from service. At the time of failure, this cavity often reaches a size of about 5 sq. in. and about 1/2in. deep. Sometimes cavities are found on only two or three sides and at other times all twelve sides appear to be invoIved.

Jan 1, 1937

By W. J. Reagan

IN a study of the life of ingot molds, it is essential to eliminate all of the variables. In the commercial manufacture of steel this is almost an impossibility. In this study many of the variables have been eliminated, and because of the unusually consistent furnace practice, with small variations in the type of steel manufactured, and a similar type of mold design for all ingot sizes, it is possible to determine with some degree of accuracy the cause of variations in ingot-mold life. All of the ingot molds in this study were used in the production of basic open-hearth killed steel with chemical analyses within the following range: C, 0.50-0.85 per cent; Mn, 0.60-0.80; Si, 0.20-0.30; S, 0.04 max.; P, 0.04 max. All ingot molds were of the closed-bottom type and used for bottom casting only. None of the molds in the study failed from cutting of the sides due to top casting; in fact, on only rare occasions are any molds used for top casting. The general type of this mold is shown in Fig. 1. All sizes have rather a heavy band of metal, about 6 in. wide, at the top to give mechanical strength at this point, and all have cast lugs. Cross sections of the vari-ous sizes from top and bottom sections are shown in Fig. 2. Ratios of wall thickness to ingot weight, ratio of ingot diameter to ingot length, and other figures of interest are shown in Table 1. In. the majority of cases, the cause of failure in this type of mold is due to heavy fire cracking on the "flats" at a point from 3 to 6 in. from the bottom radius of the mold. Fig. 3 pictures this type of failure. The start of failure at these points usually begins with a very small section of the interior surface becoming involved. As a rule, a section of 1 or 2 sq. in. is the first noticeable 'indication that mold failure is taking place. Gradually the small section increases until a sizable cavity develops, causing a protuberance on the bottom of the ingot, which eventually becomes too large to be stripped. The mold is then removed from service. At the time of failure, this cavity often reaches a size of about 5 sq. in. and about 1/2 in. deep. Sometimes cavities are found on only two or three sides and at other times all twelve sides appear to be involved.

Jan 1, 1936

By W. G. Davenport, M. Paneni

Experiments in sheets of mercury and water have confirmed that the dynamics of bubble rise are simi-lur in Liquid metal and aqueous systems. RECOGNITION that bubble/metal interactions are important in pyrometallurgical refining processes (e.g., steelmaking, vacuum degassing, deoxidation of copper) has led to a number of investigations into the dynamic behavior of bubbles rising in liquid metals. Terminal velocities of bubbles rising in mercury (20°C)' and silver (1000" to 1050"c)1,2 have been measured. Good agreement between aqueous and liquid metal systems has been obtained. Shapes of air bubbles rising in mercury have been determined experimentally by Davenport, Bradshaw, and Richardsonl who used a system of electric probes to measure maximum horizontal and vertical dimensions, Fig. 5. In the present work the motion and shape of nitrogen bubbles rising in a thin sheet of mercury (i.e., a "two-dimensional" environment) have been studied with the view of confirming dynamic similarity between aqueous and metallic systems. The use of thin sheets of liquid permitted direct visual and photographic comparisons of bubble behavior in clear liquids and in mercury. EXPERIMENTAL METHOD The thin sheets of mercury or water were held in a clear acrylic plastic container, Fig. 1. The sheets of liquid were 80 cm high by 30 cm wide and 0.47 cm thick. Single bubbles were formed by injecting the desired volume of nitrogen into an inverted semicircular cup (stainless steel). The cup was then rotated into the upright position using an external handle thus releasing a single bubble of measured volume. Nitrogen was injected into the cup from a nitrogen cylinder through the external handle and a hollow rotating axle. The volume of a bubble (volume range 0.01 to 10 cu cm) at any instant was obtained by connecting a soap film meter to a small gas space above the liquid. By this technique any change in the volume of a gas bubble in the liquid (i.e., during cup filling or during rise) produces an equivalent displacement from gas space A into the buret B. This displacement is determined by measuring the movement of a horizontal soap film in the buret. Bubble shapes were determined using single-frame photography. Rising velocities were determined by cinephotography with a 1/100-sec stopwatch in the field of view. In both cases dimensions were referenced to a rule attached to the container surface, also in the field of view. Linear dimensions were measured on projections of the single-frame negatives. Rising velocities were calculated from measured heights of rise (10 to 20 cm) and times of rise. The effects of disturbances caused by cup rotation and by bubble release were minimized by photographing bubbles at a point not less than 25 cm above the cup. MATERIALS Mercury: Technical-grade purified: minimum 99 pct Hg; Johnson Matthey and Mallory Limited, Toronto. Nitrogen: Technical-grade purified nitrogen: 99.99 pct Nz, 0.01 pct 02, trace argon; Air Liquide Canada Limitee, Montreal.

Jan 1, 1970

By E. H. Herron, J. H. Perry

Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and to predict reservoir behavior, thereby leading to the most economically desirable form of exploitation. In addition, simulation can be used as a tool for reservoir description to learn more about the physical nature of the reservoir and the mode of primary recovery. This use is essential in most reservoir studies and represents one of the more significant applications of simulation. Prediction of reservoir behavior provides information concerning displacement efficiency, optimum well locations, and the comparison of alternative producing processes. Since oil, water and gas all commonly occur in many reservoirs, a simulation that takes into account the simultaneous flow of three immiscible phases is a prerequisite for obtaining such information. Furthermore, a simulation that describes this flow in two dimensions provides the capability of evaluating combination gas- and water-drive reservoirs with either an areal or a cross-sectional description. Many reservoir problems require such a simulation. The objective of the research described here was to develop a mathematical model and associated computer program for accurate, efficient, and economical prediction of the reservoir flow of three phases in two-dimensional geometry. Several authors have discussed multidimensional reservoir simulation of two phases.l-3 Gottfried et al.' presented an analy- sis of three-phase simulation in one dimension. Fagin and Stewart" and Garrett6 have discussed the simulation of three-phase flow in two dimensions, but their solution techniques were different from the method discussed in this paper. In particular, previous models and computer programs have usually involved either fewer capabilities or inefficient and less rigorous descriptions resulting in limited applicability. In the following sections, we present a description of the model, an evaluation of the method, and the application to reservoir problems. Additional and supporting material is presented in the Appendix. The Three-Phase Model General Description The three-phase model consists of the mathematical description of the flow of oil, water and gas in a reservoir. Conceptually this description considers flow in all three space dimensions, but in the subject simulator it is assumed that no flow occurs in the third dimension. The model includes the effects on reservoir behavior of fluid and rock compressibility, viscosity, gravity, capillary pressure, relative permeability and gas solubility. Within the two-dimensional restriction, the reservoir considered may be completely heterogeneous and anisotropic. The differential equations governing multiphase, multidimensional flow in a reservoir are too complex to be solved analytically, so they are approximated by a system of finite difference equations and solved

Jan 1, 1970

By M. Cross

INTRODUCTION The results of the Japanese dissective investigations have [l-4] transformed our understanding of how the blast furnace operates. These studies revealed an internal structure of the furnace burden with five main zones, as illustrated in Figure 1. In the granular zone the charged material retains its structure whereas in the cohesive zone the iron ore particles are at advanced stages of reduction and form a semi-molten mass with a relatively high flow resistance. In the active coke zone molten iron and slag flow through a matrix of loosely packed coke. The raceway is a relatively small void space generated by displacement of coke due to the hot blast momentum, whilst the hearth is packed with coke, slag and molten iron. [ ] It is now well appreciated that furnace control is largely achieved and maintained by using the burden charge distribution to determine the gas glow distribution within the furnace and, hence, the shape and extent of the cohesive zone. In fact, Wakayama et a1 [5] of Nippon Steel Corporation state that optimum furnace operation is synonomous with the formation of a desirable cohesive zone shape, which in turn is almost solely determined by the burden charge distribution and its mechanical properties [6]. Much of the development work on operating furnaces since the dissective investigations has been oriented towards the identification of external furnace measures which reliably predict the shape of the cohesive zone. Currently the most popular measures are top gas temperature, velocity and CO/CO, ratio profiles [5]. More recently pressure distributions as evaluated from a series of wall measurements down the height of the furnace are also being utilised [7] to identify the position and shape of the cohesive zone. A great deal of emphasis has been placed upon achieving a consistent quality in the burden charge materials, establishing a suitable set of operating conditions and only making changes when forced to, ie., the operation for some reason moves outside a desired operating range or the production requirements change. Operational control schemes based upon this philosophy (eg., Kawasaki's GO-STOP system [8]) have been remarkably successful at achieving both high production rates and relatively low fuel costs. Mathematical models have long been used as tools to assist in the analysis of the blast furnace process. Although, the models have been useful in post-operative analysis and for identifying trends, etc., they have not proved reliable for detailed design criteria or as the basis for effective control schemes. Such mathematical models have moved from simple statistical correlation, through overall and stagewise heat-and-mass balance calcu-

Jan 1, 1984