CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION

Similar documents
= z 20 z n. (k 20) + 4 z k = 4

Uncertainty in measurements of power and energy on power networks

Lecture 21: Numerical methods for pricing American type derivatives

Difference Equations

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Lecture 12: Discrete Laplacian

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

EEE 241: Linear Systems

Solution Thermodynamics

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

The Geometry of Logit and Probit

Advanced Quantum Mechanics

Convexity preserving interpolation by splines of arbitrary degree

Errors for Linear Systems

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

Structure and Drive Paul A. Jensen Copyright July 20, 2003

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

THE SUMMATION NOTATION Ʃ

NUMERICAL DIFFERENTIATION

6.3.7 Example with Runga Kutta 4 th order method

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

Finite Element Modelling of truss/cable structures

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

One-sided finite-difference approximations suitable for use with Richardson extrapolation

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

Mathematical Preparations

Module 9. Lecture 6. Duality in Assignment Problems

5 The Rational Canonical Form

Pulse Coded Modulation

x = , so that calculated

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

Fundamental loop-current method using virtual voltage sources technique for special cases

Foundations of Arithmetic

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Physics 181. Particle Systems

Assignment 4. Adsorption Isotherms

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Formal solvers of the RT equation

Report on Image warping

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Some modelling aspects for the Matlab implementation of MMA

Differentiating Gaussian Processes

Inductance Calculation for Conductors of Arbitrary Shape

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

18.1 Introduction and Recap

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

2 Finite difference basics

Numerical Heat and Mass Transfer

Note 10. Modeling and Simulation of Dynamic Systems

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)

Notes on Frequency Estimation in Data Streams

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation

Lecture Notes on Linear Regression

PHYS 705: Classical Mechanics. Calculus of Variations II

Credit Card Pricing and Impact of Adverse Selection

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

8.6 The Complex Number System

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Indeterminate pin-jointed frames (trusses)

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS)( ) 8/25/2010

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

Important Instructions to the Examiners:

A new Approach for Solving Linear Ordinary Differential Equations

1 Convex Optimization

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

DUE: WEDS FEB 21ST 2018

Linear Regression Analysis: Terminology and Notation

Relevance Vector Machines Explained

Grid Generation around a Cylinder by Complex Potential Functions

Grover s Algorithm + Quantum Zeno Effect + Vaidman

Army Ants Tunneling for Classical Simulations

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Curve Fitting with the Least Square Method

Poisson brackets and canonical transformations

Chapter 13: Multiple Regression

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering /

Min Cut, Fast Cut, Polynomial Identities

10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization

Supplementary Notes for Chapter 9 Mixture Thermodynamics

CALCULUS CLASSROOM CAPSULES

The Feynman path integral

Color Rendering Uncertainty

Transcription:

Górnctwo Geonżynera Rok 0 Zeszyt / 006 Igor Konstantnovch Mladetskj * Petr Ivanovch Plov * Ekaterna Nkolaevna Kobets * Tasya Igorevna Markova * CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION The separaton characterstc of the scheme s determned by separaton characterstcs of each operaton ther number and sequence of connecton. Let s assume that before mxture the quantty of partcles wth separatng sgn X was P and after mxture t was P (Fg. ) and was equal to the sum of quanttes of partcles before mxture and entered nto the gven operaton wth crculatng loadng.e. m + j j Р Р P where j s a seral number of crculatng loadng. The crculaton s formed by consecutve connecton of devces.e. n P P E j + Then m n P P + P E j + Whence n m P P E j + * Natonal Mnng Unversty of Ukrane Dnepropetrowsk Ukrana 9

Thus the mxture operaton separaton characterstc s P CM P m n /( P j + E Let s apply the receved formulas on examples. Fg.. Consecutve connecton of separaton and mxture operatons 0

We assume there s a technologcal block for separaton (Fg. ). Fg.. The scheme of a technologcal separaton block Let s make the separaton characterstc of ths connecton. Accordng to laws of consecutve connecton P β PP OCP For the feedback technologcal connecton the characterstc s Р ОС РР Makng dentcal transformatons we shall receve the requred formula: РР РР РР Р β РР Р( Р) Р+ РР Let s complcate the separaton devces connecton (Fg. ). The sequence of separaton characterstcs for such connecton s as follows Р Р РР Р Р β FB FB and the feedback characterstc s Р FB РР Р РР Р FB FB

Makng dentcal transformatons we shall receve: PPP P β PP 5 PP 5 We shall not carry out the further transformatons as they are trval. The connecton separaton characterstc turns out rather complex but there are dentcal transformatons resultng n an obvous knd expresson. Fg.. The scheme of a complcated technologcal separaton block However for the crcut shown n Fgure the dentcal transformatons become so bulky that the probablty of a mstake comes near to one. Therefore for such scheme t s necessary to make the equatons of fractons balance for each separaton recepton: ( ) ( ) р Р р + р + р ; р Р р + р + р ; 0 0 ( ) ( ) р Р р + р ; р Р р + р ; р Р р ; р Р р ; () р Р р ; р Р р ; p p + p OC To determne outputs of all products t s necessary to solve a system of equatons (). As t has a hgh order the soluton s carred out numercally for example by a smple teratons method. Numercal solutons demand specal formal transformaton to obtan an algorthmc form.

Fg.. The scheme of a technologcal separaton block wth deep feedback For a system of lnear equatons t looks lke: uur() uur( ) ur Х АХ + В () where B ur s a vector of free terms; A s a matrx of factors; Х uur s a vector of requred roots approxmatons. The matrx of factors A wll be as follows Р Р Р Р Р Р Р Р Р Р Р Р uur and the vector of free terms s Х ( РХ0 Р Х0000000). For the analyss of the separaton scheme separaton characterstcs of all devces are set. As the separaton s bnary each devce has two separaton characterstcs. As the result the scheme shown n Fgure s characterzed by eght separaton characterstcs. They are

represented as dscrete set values for example by sx values. Thus the ntal data on the separaton scheme wll be as a matrx of 6 8 sze. The calculaton of the separaton crcut can be done n the followng way. The frst column of the separaton characterstcs matrx s chosen and the system of the equatons () s solved as a result we shall receve the vector of the frst fracton of partcles separaton. Than the second column of the separaton characterstcs matrx s chosen and agan the system of the equatons s solved. The vector of the second fracton of partcles separaton s receved. The soluton s made as many tmes as many columns are n the separaton characterstcs matrx. As a result we shall receve a new matrx of 6 8 sze. To determne the -th product output t s necessary to sum the elements of the correspondng -th soluton matrx row. To determne qualty ndces t s necessary to multply row elements and the fracton qualty parameters and to sum the receved products. So for the crcut n Fgure 5 the soluton s obtaned as the matrx 0.00 0. 0.00 0. 0.00 0. 0.00 0. 0.09 0.8 0.05 0.97 0.0 0.95 0.06 0.09 0.5 0. 0. 0. 0.07 0.6 0.0 0.07 0. 0.6 0. 0. 0. 0.09 0.066 0.0 0. 0.08 0.8 0.07 0.8 0.0 0.088 0.0 0.5 0.0 0. 0.0 0.6 0. 0. 0.006 0.6 0.6 The ffth product s the enrched product and the eghth one s depleted. The sum of outputs of these products s equal to one as t s receved as a result of calculatons. The technologcal balance of concentraton processes s connected wth the determnaton of output parameters of each technologcal pont. When the technology has deep recycles such calculatons become nconvenent and are connected wth teratve procedures to select the requred parameters resultng n materal balance. The purpose of the gven argumentatons s to fnd formal algorthm whch would allow makng calculatons by the computer. Let s assume that we have some separaton crcut (Fg. 5). Fg. 5. The separaton crcut wth the recycle

Let s wrte down a system of balance equatons for the gven crcut: Q ( Q + Q ); Q ( Q + Q ); 0 0 Q Q ; Q Q The number of the equatons s equal to the number of unknown values therefore the soluton s sngle-valued and takng nto account the low order of the system t can be solved by a substtuton. We shall determne an output of the enrched product: Q Q 0 ; Q Q0 Thus on the bass of known outputs on separate operatons and known laws of products transformaton n separaton crcuts t s possble to determne the output of any product n relaton to an ntal feedng. Indvdual outputs are calculated when requred by the Q IJ rato of charges IJ. Q0 Nowadays the mathematcal approach s developed allowng calculaton of estmatng accuracy of a product output for a rather smple technology of mneral wealth enrchment []. Such estmaton s carred out wth the help of partal dfferentaton of the product output functon. The more complcated s the crcut the more complcated s the operaton of dfferentaton. Therefore t s necessary to have a method allowng to execute the mentoned estmaton of accuracy for the crcut of any complexty. So f to calculate outputs at known parameters of the product preparaton and then to gve ncrements to these parameters and agan make calculatons t s possble to receve an estmaton of output dervatves on correspondng parameters that serves as a measure of accuracy of nterested parameters estmaton. Separaton characterstcs Р(α) and functons of partcles dstrbuton on fractons F(α) can be receved expermentally. For ths reason each fracton defnton error one can consder as known: σ σ...... σ σ σ...... σ р α р α р αп F α F α F αn The complex crcut product output wll be calculated on the bass of the repeated soluton of system of the lnear equatons of products balance. At the begnnng we assume that n the result the certan value s obtaned. Knowng t we gve an ncrement to the functon Р(α). For example we shft t upwards or downwards and agan calculate an output. Now t wll be. 5

For the error defnton t s necessary for functons F(α) and Р(α) to have dervatves on all the range α. In ths case the dfferental of the functon wll be: ΔP( α ) P( α ) P ( α ) ΔΔF( α ) ΔF( α ) ΔF ( α ) Then the dfferental of the argument s: Δ ( α ) P( α) ΔF( α) P( α) ΔF( α ) Δ P F( α ) P( α ) ΔF ( α ) P( α ) ΔF ( α ) The functon dervatves can be wrtten down as follows: P Δ P( α) ; P ΔP( α ) F Δ F ( α ) F ΔΔF( α ) Then the error from each varable change wll be formed from errors of each product fracton defnton as functons F(α) and Р(α) are ntegral: n ( ) Δ p α σ P σpα ; ΔP( α ) n ( ) Δ p α σ F σfα ΔF( α ) As a result the general error of the numercal output estmaton can be wrtten by the followng expresson: σ σ +σ Р F As an example we shall gve the calculatons of the output error for the separaton crcut shown n Fgure 6. ) We determne the product output of the complex crcut on the bass of the repeated soluton of the system of lnear equatons of products balance. As a result we have receved the certan value. 6 Knowng t we shft the functon Р(α) downwards and agan calculate the output.

Fg. 6. The technologcal separaton block ) Usng the gven equatons we have receved ther solutons whch are regstered n Table. ΔP( α ) P( α ) P ( α ) ΔΔF( α ) ΔF( α ) ΔF ( α ) TABLE α 0 0.5 0.75 0.65 0.875 Р(α) 0.0 0.0 0.06 0.0 0.0 0.05 F(α) 0.08 0.0 0.0 0.0 0 0.0 We calculate an argument dfferental. The receved values are regstered n Table. Δ ( α ) P( α ) ΔF( α ) P ( α ) ΔF( α ) P Δ ( α ) P( α ) ΔF( α ) P( α ) ΔF ( α ) F TABLE α 0 0.5 0.75 0.65 0.875 p (α) 0.006 0.0065 0.0 0.009 0.006 0.058 F (α) 0.0096 0.005 0.008 0.007 0 0.08 7

Functon dervatves are wrtten down as follows: P Δ P( α) ; P ΔP( α ) F Δ F ( α ) F ΔΔF( α ) The values receved at ther calculaton are gven n Table. TABLE α 0 0.5 0.75 0.65 0.875 p (α)/ Р 0. 0.5 0.8 0.0 0. 0.76 F (α)/ F 0. 0.5 0. 0.7 0 0.8 The error from each varable change s calculated as follows: n ( ) Δ p α σ P σpα ; ΔP( α ) n ( ) Δ p α σ F σfα ΔF( α ) We have receved the values of each varable and have wrtten them n Table. TABLE σ P 0.000005 0.0000 0.00005 0.000009 0.00000 0.00005 σ F 0.00000 0.00000 0.00006 0.00009 0.000000 0.000067 Now t s possble to calculate the general error of numercal output estmaton: σ σ +σ Р F σ 0 00009 + 0 0006 0 0008 σ 5% 8

Ths crcut s not too complex and so we stll can calculate an output error by a classcal method wth the purpose to check the offered numercal approach. Devces separaton characterstcs are dentcal and are gven n Table 5. There also the dstrbuton functon of jonts of raw materal gong to separaton s gven. TABLE 5 α 0 0.5 0.75 0.65 0.875 Р 0. 0.5 0. 0.7 0.78 0.8 ΔF 0.7 0.6 0.7 0.08 0. 0. Frst t s necessary to determne outputs of all products. For ths purpose the separaton characterstcs for each technology pont are set. The feedback separaton characterstc s Р FB РРР The feedback separaton characterstc s Р FB РРРР For each crcut pont the separaton characterstc s formed by the scheme of multplcaton of separaton characterstcs from the nput up to a pont. Pont : Р Ι РРР РРРР + РРРРРРР The denomnator at all separaton characterstcs wll be the same therefore further we shall wrte down: pont : Р РР Ι pont : Р РР Р Ι pont : Р V РР Р Р Ι pont 5: РV РРРРР Ι pont 6: Р V Р Р Ι pont 7: Р pont 8: Р pont 9: Р V РР Р Ι V РР Р Р Ι РРРРР Ι Х 9

Let s calculate a denomnator of separaton characterstcs (Table 6): TABLE 6 α 0 0.5 0.75 0.65 0.875 Р 0. 0.5 0. 0.7 0.78 0.8 Denomnator 0.9858 0.978 0.869 0.765 0.7757 0.797 Let s wrte down all separaton characterstcs calculated accordng to the above-stated expressons (Table 7). TABLE 7 α Р Р Р Р Р 5 Р 6 Р 7 Р 8 Р 9 5 6 7 8 9 0 0.0 0.7 0.06 0.007 0.000 0.896 0.07 0.08 0.005 0.5.0 0.5 0.00 0.00 0.0005 0.8690 0.0 0.095 0.009 0.75.50 0.60 0.80 0.076 0.09 0.690 0.760 0.0 0.0 0.65.068 0.97 0.60 0.8 0.7 0.90 0.7 0.9 0. 0.875.89.0055 0.78 0.67 0.77 0.86 0. 0.75 0.5.60.057 0.89 0.696 0.570 0.7 0.86 0.58 0.5 We calculate the output and qualty n each pont of technology. Pont (Table 8). TABLE 8 α Р ΔF 5 6 0.0 0.7 0.7 0 0 0.5.0 0.6 0.658 0.09 0.05 0.75.50 0.7 0.06 0.675 0.05 0.65.068 0.08 0.05 0.065 0.05 0.875.89 0. 0.569 0.557 0.05.6 0. 0.605 0.605 0.. 0.8 0.89 0

α 0.89 0.9 n β 0. Smlarly we calculate for other ponts of technology. Therefore we have: 0.5 β 0.6 0.9 β 0.76 0. β 0.8 5 0.5 β 0.85 6 0.65 β 6 0.6 7 0.9 β 7 0.5 8 0.09 β 8 0.6 9 0.05 β 9 0.7 Let s calculate partal outputs for each separaton recepton. The frst recepton: 0. 0.6 0. 0.6 0.6 0.58 The second recepton: 0.6 0.5 0.6 0.76 0.5 0.9 The thrd recepton: 0.76 0.6 0.7 0.8 0.6 0.

The fourth recepton: 0.8 0.7 0.77 0.85 0.7 0. The functon of the concentrate output s made smlarly to the separaton characterstcs and looks lke: ( ) ( ) + ( )( ) Let s generate expressons for partal dervatves. Accordng to the calculated partal outputs we shall wrte down the expressons: 0.9 β 0.8 + 0.08 0.79 β 0.9 0.9( ) 0.6 β 0.9 + 0.0085 0.97 β 0.759 + 0.8 Let s calculate dervatves on correspondng outputs. 0.9( 0.8 + 0.08 ) 0.9 ( 0.8+ 0.06 ) 0.8 + ( 0.8 0.08 ) 0 6( 0.9 + 0 0085 ) 0.6 ( 0.9 + 0 07 ) 0.9 + ( 0.9 0.0085 ) 097(0 759 + 08 ) 0.8 097 0.9 + (0759 0.8 ) 079(0.9 0.9( )) 0.79 0.9 0.69 (0.9 0.9( ))

α ν Let s further calculate partal dervatves from partal outputs on correspon- β ν dng varables. 0.8. α α 0.8 0.9 0.9 α α 0. 0.9 0.95 α α 0. 0.69. α α 0. 0 0 6 β 0 6 0.6 0.8 0.6. β β ( β 0.6) Let s assume σ α σ β σν. To determne ther numercal values whch wll be the same as for numercal determnaton of errors; let s do the next actons (Fg. 7). Fg. 7. The crcut of recepton of output functon arguments ncrements

Accordng to functons ncrements we can determne an argument ncrement Δα for each accepted dscrete value of the functon usng the plot. Then we can fnd the average value of such ncrement on all a range of α change. It wll be the requred Δα. In ths case t we have 0.6 0.5 β 0.5 Δα 0.0 0.9 0.9 0.57 β β ( β 0.5) 0.76 0.6 β 0.6 0. 0.9 0.665 β β ( β 0.6) 08 07 β 07 0 ν 06 ν 06 ν 076 ν (0.76 ν)( ) + (0.6 ν) 0.9 0.6 ν ν (0.76 ν) 076 ν 08 ν (0.8 ν)( ) + (0.76 ν) 0.9 0.85 ν ν (0.8 ν) 08 ν 085 ν (0.85 ν)( ) + (0.8 ν) 0.69 0. ν ν (0.85 ν)

α β ν σ σ + σ + σ α β ν σ (. 0.0) + ( 0.6 0.0) + ( 0.67 0.0) 0 00 σ 0.0 α β ν σ σ + σ + σ α β ν σ 0.0 α β ν σ σ + σ + σ α β ν σ 0.0 α β ν σ σ + σ + σ α β ν σ 0.0 The error functon s: σ σ + σ 0.08 + σ + σ Though the dvergence s about 5% we beleve that the approach of the error s numercal analyss s comprehensble as opportuntes of such analyss are practcally unlmted. Concluson The offered method of the output determnaton of the accuracy estmaton for enrchment products can be appled for enterprse technologcal balance formng. REFERENCES [] Kozn V.Z.: Approbaton and control of technologcal processes of enrchment. М. Nedra 985 9 5