root/trunk/matml/src/sphlow/sphlow.html

Revision 288, 8.8 kB (checked in by powell, 4 years ago)

Significant updates to Sphlow HTML page.

  • Property svn:keywords set to Author Date Id Revision
Line 
1<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
2<html>
3<head>
4   <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
5   <meta name="Author" content="Adam C. Powell, IV">
6   <meta name="GENERATOR" content="Mozilla/4.75 [en] (X11; U; Linux 2.2.17 i686) [Netscape]">
7   <title>Sphlow Flotation Calculator 0.2 by Adam Powell</title>
8</head>
9<body text="#000066" bgcolor="#CCFFFF" link="#003300" vlink="#330033" alink="#990000">
10
11<center>
12<h1>
13Sphlow Flotation Calculator Version 0.2</h1></center>
14
15<center><applet code="Sphlow.class" width=800 height=500 codebase=.>You
16are using a Java-challenged browser.  Sorry!</applet><p>
17
18<a href="#targ0">About Sphlow</a> -
19<a href="#targ1">How to use it</a> -
20<a href="#targ2">How it works</a> -
21<a href="#targ3">Discussion</a> -
22<a href="#targ4">Known bugs</a> -
23<a href="#targ5">Credits and source</a></center>
24
25<h3><a NAME="targ0">About Sphlow:</h3>
26Sphlow was written to model metallurgical processes involving flotation of
27inclusions, such as tundishes or electron beam melting hearths. Inclusions
28are second-phase particles, such as oxide or nitride impurities, tungsten
29carbide tool bit chips left over from machining operations, etc. which
30get into metal and must be removed while the metal is molten. A given
31tundish or hearth design and operating parameter set will lead to a rough
32critical rising/sinking velocity for removal of all inclusions, such that
33(nearly) all inclusions rising or sinking faster than that critical velocity
34will be removed in that vessel. So, if inclusions are distributed
35in size and density (<i>D</i> and <i>rho<sub>p</sub></i>), then since velocity
36increases with diameter, those particles in the distribution above the
37critical velocity contour in the graph on the right will be removed from
38the molten metal by floating to the top or sinking to the bottom of the
39vessel.
40<p>Of course, you're free to use it for whatever purpose you like.</div>
41
42<h3><a NAME="targ1">How to use it:</h3>
43To draw the first plot, hit return in the velocity field. You may
44then use the menus and change the value in any field, just make certain
45you press &lt;return&gt; after each entry so that it will register.
46To zoom in on a graph, drag the mouse from the top left corner to the bottom
47right corner of what you'd like to see. To zoom out, drag the mouse
48upward. You can also change the limits in the <b>"D-rho curve controls"</b>
49panel.
50<p>So, type in your fluid viscosity and density, or choose from the liquids
51offered. When you give it a critical velocity, Sphlow will tell
52you the Reynolds number, friction factor and particle diameter and show
53you the relevant point for the given particle density in the friction factor-Reynolds
54number curve on the left. It will also plot diameters which rise/fall
55at the same velocity for a range of particle densities in the curve on
56the right. You may plot up to nine such velocity contours by selecting
57a new curve number and entering a new velocity. If you change liquid
58properties or <i>D-rho</i> graph limits, Sphlow will redraw all of
59the velocity contours.
60<p>The red vertical line on the right indicates the liquid density.
61The starting point corresponds to nickel spheres in water.</div>
62
63<h3><a NAME="targ2">How it works:</h3>
64This applet calculates and displays the relationship between size, density
65and terminal rising/sinking velocity of a spherical particle in a fluid.
66It is based on the relationship between friction factor and Reynolds number
67in laminar flow, from Stokes flow where <i>Re</i>&lt;0.2 to Newton's law
68flow for Reynolds numbers from about 10<sup>3</sup> to 10<sup>5</sup>.
69At each point, it first calculates the ratio <i>f/Re</i> (friction factor
70divided by Reynolds number) given by:
71<br>
72<center><table CELLSPACING=0 CELLPADDING=0 NOSAVE >
73<tr ALIGN=CENTER NOSAVE>
74<td ALIGN=CENTER NOSAVE><i><u>f</u></i></td>
75
76<td ALIGN=CENTER VALIGN=CENTER ROWSPAN="2" NOSAVE>=</td>
77
78<td NOSAVE><u>4 <i>mu g</i> |<i>rho<sub>p</sub></i>-<i>rho<sub>f</sub></i>|</u></td>
79</tr>
80
81<tr ALIGN=CENTER NOSAVE>
82<td ALIGN=CENTER NOSAVE><i>Re</i></td>
83
84<td><i>rho<sub>f</sub></i><sup>2</sup><i>u</i><sup>3</sup></td>
85</tr>
86</table></center>
87where <i>mu</i> and <i>g</i> are the viscosity and gravitational acceleration
88respectively, <i>rho<sub>f</sub></i> and <i>rho<sub>p</sub></i> the fluid
89and particle densities, and <i>u</i> the rising or sinking velocity.
90It then draws a diagonal line of slope 1 in the graph on the left corresponding
91to that value of <i>f/Re</i>, and the intersection with the red curve gives
92the Reynolds number and friction factor. Since
93<br>
94<center><table CELLSPACING=0 CELLPADDING=0 NOSAVE >
95<tr ALIGN=CENTER NOSAVE>
96<td ALIGN=CENTER VALIGN=CENTER ROWSPAN="2" NOSAVE><i>Re</i></td>
97
98<td ROWSPAN="2" NOSAVE>=</td>
99
100<td><i><u>rho<sub>f</sub> u D</u></i> ,</td>
101</tr>
102
103<tr ALIGN=CENTER VALIGN=CENTER NOSAVE>
104<td NOSAVE><i>mu</i></td>
105</tr>
106</table></center>
107the particle diameter <i>D</i> follows straightforwardly.</div>
108
109<h3><a NAME="targ3">Discussion:</h3>
110This model assumes a quiescent liquid. If there is a lot of vertical
111flow, recirculation or turbulence in your problem, particles will be swept
112up and down, so it is wise to be conservative in your estimate of critical
113velocity (<i>i.e.</i> use a <b>much</b> larger velocity than the vessel
114height divided by residence time).
115<p>The model also assumes constant liquid and particle densities.
116If these densities vary strongly with temperature, things will be somewhat
117more complicated: if the liquid coefficient of thermal expansion (CTE)
118is higher than that of the particle, there will be a range of densities
119at which particles will be neutrally buoyant; if the particle CTE is higher
120then it will always float or sink. Density also changes in porous
121inclusions, <i>e.g.</i> porous titanium nitride in molten titanium, where
122dissolving the dense nitride decreases a particle's average density but
123filling a pore increases it. (See the <a href="http://www.tms.org/pubs/journals/MT/B/9712/abstracts-9712.B.html#5">paper
124of Jean-Pierre Bellot and Alec Mitchell</a> in the <a href="http://www.tms.org/pubs/journals/MT/B/9712/contents-9712.B.html">12/1997
125issue</a> of <i><a href="http://www.tms.org/pubs/journals/MT/MT.html">Metallurgical
126and Materials Transactions</a></i> for details. Hey, while you're
127there, check out <a href="http://www.tms.org/pubs/journals/MT/B/9712/abstracts-9712.B.html#29">my
128paper</a> too!)
129<p>Finally, the model assumes spherical inclusion particles. For
130non-spherical particles, the friction factor will usually be smaller than
131that of the smallest sphere containing the particle (always so in the case
132of Stokes flow where <i>Re</i>&lt;0.2), so velocity will be higher, and
133this model will give a conservative velocity estimate for such particles.</div>
134
135<h3><a NAME="targ4">Known bugs:</h3>
136Due to a bug in the ptplot package, the "fill" button will redraw the plot
137with limits set to encompass all points that have ever been drawn, not
138merely those which are currently active. It should not be necessary
139to use this button. Another bug in ptplot limits the minimum size
140of the graphs, so a version smaller than 800x500 will not be feasible anytime
141soon. Also, some of the velocity contours are not drawn to the top
142of the <i>D-rho</i> graph space; I will fix that at some point.
143<p>The javascript layers business seems to be broken on newer versions
144of Netscape. Oh well, I may try to fix it someday.</div>
145
146<h3><a NAME="targ5"><a href="http://www.jars.com/"><img SRC="jarsregistered.gif" ALT="JARS-registered" BORDER=0 height=93 width=141 align=RIGHT></a>Credits
147and source:</h3>
148Sphlow is based in large part on the <a href="http://ptolemy.eecs.berkeley.edu/java/ptplot/index.html">ptplot
149package</a> at UCBerkeley, copyright Regents of the University of California.
150I also learned a lot from the <a href="http://www.phrantic.com/scoop/onjava.html">Scoop
151on Java</a> (where my <a href="http://lyre.mit.edu/~powell/Software/Java/hello.html">first applet</a> came from),
152and the <a href="http://lyre.mit.edu/~powell/Software/Java/lemniscate.html">Lemniscate</a>,
153among the <a href="http://www-history.mcs.st-and.ac.uk/~history/Java/">curves</a>
154from the <a href="http://www.st-andrews.ac.uk/">University of St. Andrews</a>.
155Byte code for most of the classes is due to the SGI JDK 1.1 (DiamCalculator.class
156was recompiled in May, 2000 using guavac 1.2 to fix a bug), but it should
157all be JDK 1.0 compliant.
158The Magic Red Curve in the friction factor graph on the left is based on
159principles of fluid flow around a sphere, and taken from <i><a href="http://www.amazon.com/exec/obidos/ASIN/047107392X">Transport
160Phenomena</a></i> by Bird, Stewart and Lightfoot.
161
162<p>Get the <a href="./">source
163code here</a> (distributed under <a href="http://www.gnu.org/copyleft/gpl.html">GPL</a>).
164<p>Email compliments/complaints/comments, bug reports and suggestions for
165new liquids to <a href="mailto:Adam Powell <apowell@opennovation.com>">Adam
166Powell</a>.</div>
167
168</body>
169</html>
Note: See TracBrowser for help on using the browser.