Curator Welcome
Time Dependent Wave Equation 1D One Dimensional Optical Resonator Modes
Gas Physics Introduction
Gas Physics Introduction Chapter 1:Energy Distribution Evolution Energy Diffusion from Mono-energetic of a Three Dimensional Gas Rotating Gas in a 3D Toroid Gas Velocity Equilibration 2D Gas Velocity Equilibration 3D Circulating Gas in a Wind Tunnel Density Modulated Rotating Gas in a 3D Toroid Position Diffusion (Mixing) of a Three Dimensional Gas Energy Distribution 3D with MathJax Equations and Embedded Canvases Gas Energy Distributions 2D Rigid Rotors Energy Distribution 2D 1D Gas Energy and Momentum Distribution Rotor Energies 2D Rotor Energies 3D H2O Energies 3D CO2 Energies 3D 2D Rigid Arrays with Rotation and Translation 3D Rigid Arrays with Rotation and Translation Gas Energy Distributions 1D, 2D, 3D Chapter 2:Particle Scattering Effects
Inelastic Collision Conversion to Internal Energy Brownian Motion Animation 2D Using Hard Sphere Scattering Large Disc Kinetics 2D Condensation of a Gas on a 2D Solid Maxwell's Demon Energy Seperator 2D Minimize Final Energy 2D and 3D Gas Diffusion 2D Gas Expansion 2D Linear Gas Expansion 2D Constant Pressure Gas Expansion 2D Circular Gas Explosion 2D Circular MagnetoCaloric Effect 2D Convection Loops in a Numerical Gas Chapter 3: Gas Energy Equipartition Energy Equipartition of a Two Dimensional Gas with Two Different Masses Energy Equipartition of a Three Dimensional Gas with Two Different Masses Gas Energy Equipartition 3D
Chapter 2.1
Chapter 2.2
Chapter 2.3
Chapter 3.1
Chapter 3.2
Chapter 3.3
Chapter 4.1
Chapter 4.2
Chapter 4.3
Photoelectric Effect Alpha Particle Emission from a Large Nucleus 2D Black Body Emission from Charged Harmonic Oscillators Quantum Animation of Energy Band Gap 1D Animation of Quantum Wave Packets in a pn Junction 1D Quantum Mirror Focusing by Iterating the Schrodinger Equation Quantum Lens Focusing by Iterating the Schrodinger Equation Quantum Two Mirror Resonator by Iterating the Schrodinger Equation Entangled Quantum Object Propagation in a Parabolic Potential Classic and Quantum Object Propagation in a Parabolic Potential Accurate Wave Packet Propagation in a Parabolic Potential Comparison of Quantum to Classical Physics 1D-Time Dependent Accurate Wave Packet Propagation in a Swaged Potential Series Solution of Wave Packet Propagation in a Infinite Square Potential Propagation of a Wave Packet in a Parabolic Potential Propagation of a 2D Wave Packet in a 2D Parabolic Potential Propagation of an Asymmetric Shaped Wave Packet in a 2D Parabolic Potential Two Slit Diffraction by Iterating the Schrodinger Equation Two Potential Slit Diffraction by Iterating the Schrodinger Equation Propagation of an Electron Wave Packet in an Infinite Square Well (Series Solution) Propagation of a Free Wave Packet (Series Solution) Series and Algebraic Solutions for Propagation of a Wave Packet in a Parabolic Potential Comparison of Quantum to Classical Physics 1D-Time Independent Mexican Hat Stationary Quantum States of a Finite Square Potential Well Various Wave Packet Motion in a Finite Square Potential Well Various Wave Packet Motion in a Swaged Potential Well Comparison of Quantum to Classical Physics 1D-Time Independent Cosine/Parabola Electromagnetic Modes of a 2D Cavity Black Body Radiation Experiment and Theory Emission of the Inner Walls of a Black Body Cavity Multiple Slit Diffraction Statistics of Indistinguishable Dice or Particles Quantum Field Collapse Visualization Matrix Operations and Eigenmodes
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Tumbling Block on Inclined Plane New Roller on Inclined Plane Sliding Block on Inclined Plane Classic Objects Motion 1D Generation of Curve for Fastest Time between Two Points Curve for Fastest Time between Two Points Bead Sliding On a Stiff Wire Sound Wave in Two Solid Media Mechanical Properties of a 2D Numerically Modeled Lattice Vibrating Reed Numerical Model 2D Laminar Flow of a Fluid in a 2D Bearing Animation of a Fountain 2D Gravity Leveling of a 2D Liquid Longitudinal and Shear Wave Sound Speed in a 2D Lattice Mechanical Properties of a 2D Numerically Modeled Lattice Vibration Frequencies of a 2D Numerically Modeled Lattice Damped Harmonic Oscillator Motion and Modes of a Particle in a Mexican Hat Potential Forced Oscillator Roller Bearing Cam and Roller Follower Simple Hybrid Gear Train Ratchet Animation Child's Swing Height Increase Mechanics Roller Chain Animation Elasticity of 3D Gas Driving a Piston Free to Forced Harmonic Oscillator Transition Physics of Mounting a Tire on a Rim Physics of a Pendulum Clock Escapement Physics of a Gas Pressure Gauge
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Physics of a Gas Pressure Gauge Ion Drift in Neutral Particle Sea Diffusion Equation Solved by Finite Difference Method Linear Particle Diffusion Heat Flow 1D with Time Dependent Boundary Solid Cube in Liquid Column Drag on a Large Sphere Due to a Digital Gas Animation of a Gas Centrifuge Laminar Flow of a Gas in an Annular Space Laminar Flow of a Fluid in a 2D Bearing Incompressible Fluid Flow through a Constriction Fluid Flow Due To Pressure Gradients
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Parallel Resistors Gas Discharge Animation Oscilloscope Animation Permanent Magnet in Solenoid Cross Product of Two Vectors Magnetic Braking with Cylindrical Geometry Stopping a Falling Magnet by Eddy Currents Magnetic Brake Magnetic Field of a Long Solenoid Coil Light Propagation in a Transparent Dielectric Material Dipole Array Animations Angular Distribution of Emission of Arrays of Excited Dipoles Geometric Derivation of Snell's Refraction Angle Law Electromagnetic (EM) Wave Response to a Row of Dipoles Electromagnetic (EM) Wave Response to a Single Dipole Optics of a Transparent Plate by Iterating the Maxwell Wave Equation Diffraction by Iterating the Maxwell or Schrodinger Wave Equations Lens Focusing by Iterating the Maxwell or Schrodinger Wave Equations Optical Wave Packet Propagation by Finite Difference Matrix 1D Potential(x,y) Using Laplace's Equation Div(gradV)=0 with Variable Boundary Geometry Potential(x,y) Using Laplace's Equation Div(epsilon*gradV)=0 Laplace's Equation Div(gradV)=0 with Mixed BCs Potential(x,y) for Surface Electrodes with Dielectric Slab Dipolar Media Field Plots LCR Oscillator Animation Beam Trajectory of a Moving Laser Beam Trajectory in a Square Ring Laser Doppler Effect Converting Magnetic Forces to Electric Forces Converting Permanent Magnet Orbital Moments to Solenoid Currents Electron Beam Collimation by Axial Magnetic Field Electron Buncher for Klystron Charged Particle Spatial and Speed Bunching Electrical Conduction-Particle Picture Electromagnetic Bell Triode Vacuum Tube
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
E=mc^2 and Mexican Jumping Bean Radar Pulse Return Time Interval Experiment to Determine Gamma Current Loop with a Rotating Test Particle Moving Fabry-Perot Interferometer Space Time Invariant Intervals of a Moving Clock Relativity Transformations by Requiring Simulataneity Accelerated Space Trip Clock Times and Rates Minguzzi Invariant Inertial Frame Time Accelerated Space Trip with Light Pulse Signaling
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Solar Power Generation Planet Formation Illumination of Planets Change of Orbit due to Asteroid Impact Earth Time Zones Day Lengths Astronomical Aberration Kepler Planetary and Sun Orbits Model of Ocean Tides Gravity Course Iteration Spacecraft Speed Assistance by Gravity Hohmann Transfer from One Circcular Orbit to Another Star Distances Using the Parallax Method
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Semiconductor Electron Energies and Band Gap 1D Semiconductor Electron Energies and Band Gap 2D Simple Fermi Surface p-nJunctionDiode Depletion Width Temperature Dependence of 2D Fermi-Dirac Particles
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Assembly of 3D Parts Using Javascript Web Graphics Library 3D Compared to Javascript 3D Graphics 3D Graphics Tutorial Beginnings Proof of Pythagoras a^2+b^2=c^2
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Multiple Slider Animation Template for Plotting with HTML Enhanced Sliders or Ranges Collapsible Dropdown Menu Example Arrow Functions Fast Fourier Transform (FFT) Demonstration
Animation of a Progress Bar
Template for Scientific Plots
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Time Dependent Wave Equation 1D One Dimensional Optical Resonator Modes
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Probability of Birthday Match Random Creation of Probability Function Profiles Two Dies Probability Distribution Three Dies Probability Distribution
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Light Rays in a Transparent Isosceles Prism Light Ray in a Series of Wedges of Various Index of Refraction Snells Law of Refraction Refraction Due to Dipoles in Two Media Refraction and Reflection Two Media Animation of Concave Mirror Reflection Animation of Convex Lens Refraction Prism Refraction-Wave Picture Thick Lens Ray Trace Monochromator (Czerny-Turner Type) Wavelets from Concave Mirror> Thick Lens Focusing
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Sagnac Effect Fiber Optical Gyro (FOG) Ring Laser Gyro (RLG)
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Paticle Motion in a Venturi: Wall Reflections Paticle Motion in a Venturi Potential Elliptical Ion Trap Using Boundary Potential Distribution Prediction of Hard Disc Collision Time and Position Rigid Rotor Collisions
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Acoustic Wave with Built-in Pressure Acoustic Waves in 2D Numerically Modeled Solid Wave 1D on a Linear Array Waves on Discrete Model of Vibrating String Linear Acoustic Wave in a 3D Gas Spherical Acoustic Wave in a 3D Gas Laterally Perturbed Two Dimensional Gas with Lennard Jones Repelling Forces Longitudinally Perturbed Two Dimensional Gas with Lennard Jones Repelling Forces
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Bicep and Forearm Animation Blood Pressure Measurement Multiple Dose Residual Medicine Vs Time The Eye and its Resolution Breathing Animation Scoliosis Animation Muscle Animation
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Scattering from a Crystalline Target 2D Scattering from a Crystalline Target 3D
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12
Energy Conduction in a 2D Numerically Modeled Solid Lattice Unbalanced Gas Dynamics:Energy Flow in a Gas
Chapter 4
Chapter 5
Chapter 6
Chapter 7
Chapter 8
Chapter 9
Chapter 10
Chapter 11
Chapter 12

Hover over the menu bar to pick a physics animation.

Circulating Gas in a Wind Tunnel with a Venturi

A venturi is a tapered down region (see Canvas 1) of a hollow tube where a fluid moving in the tube speeds up. I will examine the dynamics of a GAS in a venturi where the learner has control over the geometry and gas parameters. A gas is more interesting than a liquid since a gas is much more compressible than a liquid. A venturi is what was used in carburetors of internal combustion engines. The increased speed of the air in the venturi reduced the air pressure and the reduced pressure "sucked" the gasoline vapor (and sometimes liquid) out of the carburetor float bowl. A venturi is often used in wind tunnels where high wind speeds are required to measure aircraft body lift and drag. Many interesting variables like velocity, pressure, and density Vs position around the wind tunnel will be plotted.

For an incompressible fluid the speed, `v`, in the venturi is just

`v^2=A_0/A v_0^2 \ \ \(1)`

where `v_0` is the speed before entering the tapered section and venturi,`A` is the cross-sectional area of the venturi and `A_0` is the cross-sectional area of the section prior to the tapered section. Equation 1 stems simply from the fact that the rate of flow must be constant along the incompressible fluid path (no sinks or sources exist inside). The Bernoulli Principle provides for a pressure, P, reduction of the fluid

`P/P_0=v_0^2/v^2=A/A_0\ \ \(2)`

Where `P_0` is the initial pressure.  The reason for this is that fluid pressure represents a potential energy that decreases when the `v^2` increases. Equation 2 is what gives rise to lift in aircraft wings. This pressure reduction can be observed as the reduction of the force of the particles against the walls of the venturi.

It is essential to define how we compute a pressure in a numerical gas that has a variable drift (wind) speed, `v_("drift")`, along the wind tunnel axis.

In a gas that has NO drift velocity, the pressure is the average force per unit area on an imaginary surface due to the rebound of gas particles.

Hereafter, a bold character denotes a vector which has [x,y,z] components. Also, a vector with a "hat" is a unit vector.  

The change in velocity, `bb v`, on rebound is always along the normal, `bb hat "n"`, to this surface and the force, `bb "F"_p`, per particle is

`F_p= 2 m_p (bb hat "n" cdot bb "v") bb hat "n" \ \ \ (2a)`

In a gas that HAS a drift velocity, the pressure is STILL the average force on an imaginary surface due to the rebound of gas particles. However, with a drift velocity, the force causing the pressure would be huge so we subtract the LOCAL drift velocity off of the total velocity before computing the pressure. The pressure due to the flow of the gas is called dynamic pressure and should be kept separate from the static pressure. The static pressure is then due to the average velocity change along the normal to the imaginary surface and the force per particle is

`F_p= 2 m_p bb hat "n" cdot (bb "v"-bb "v"_("drift")) bb hat "n" \ \ \ \ (2b)`

where `m_p` is the particle mass, here taken to be unity to make the numerical calculations simpler. Since, after subtracting off `bb"v"_("drift")` the velocity distribution becomes isotropic, the orientation of the imaginary surface makes no difference. For this animation the pressure is computed by summing up the momenta changes of the particles hitting the inside wall of a slice of the the wind tunnel during one iteration time period. This sum, which is the Force, `F_s`, is then divided by the wall area, `A_s`, of the slice which yields the pressure on that slice.

`P_s=F_s/A_s \ \ \ (2c)`

Recirculating Wind Tunnel Geometry

I have chosen to animate a Recirculating Wind Tunnel because, with recirculation, it is possible to show conservation of energy in a vey convincing manner. The Recirculating Wind Tunnel here is composed of several objects: a 180 degree bend on the right (+y) side, a pipe at the bottom (+x) side a 180 degree bend at the left (-y) side, a tapered reducing section to the left (-y) of the center, the venturi centered (at y=0) at the top (-x) side and a tapered increasing radius section to the right (+y) of the venturi and connecting to the right hand 180 degree bend. A 180 degree bend is an object that is characterized by two parameters: 1. A median arc of radius `R_s` (similar to the equator of a sphere). 2. A set of tubular circles of radius `r_s` centered at this median .
To locate any point on the tubular circle one also needs to define the azimuth, `theta (0-180\ degrees)` pointing to a point on the median arc as well as the altitude `phi (0-360\ degrees)` of the point on the tube. Thus, an equation in Cartesian coordinates that defines the boundary point can be written

`bb "r"_v=[r_x,r_y,r_z]=[ R_s cos theta + r_s cos theta cos phi,\ \ \ R_s sin theta + r_s sin theta cos phi,\ \ \ r_s sin phi] \ \ \ \ \ (3)`

The pipes are hollow cylinders that have longitudinal parameters instead of the azimuth angle `theta \ \ ` The tapered sections are similar to pipes except the radius at the two ends differs. For viewing on the canvas, these objects can be rotated by multiplying each bend vertex position vector by a 3`x`3 rotation matrix, M.   Particles are added to the bend in groups at intervals along the `theta` angles, and to the pipes and tapered sections at intervals along the longitudinal axes.

Initial Particle Velocities

The particles will be given a combination of random velocities as well as drift speeds.   For this animation, in the bends, the drift speeds will be along the median or `theta` direction (i.e. they will be locally parallel to the bend median line).  For display only, the particles also need to be rotated with M in order to remain inside the wind tunnel.

The initial velocity vectors of the particles are defined by a speed `v_drift` and the random speed `v_r`. The speed in the bend must be matched to that in the pipes that connect to the bend. Also the drift speed in the tapered sections and the venturi should be increased to conform to equation 1. For the bends the speeds are defined by

`bb "v"=v_("drift") hat bb theta+(v_(rx)hat bb "x"+v_(ry)hat bb "y"+v_(rz)hat bb "z")`

where

`hat bb theta=-sin theta hat bb "x"+cos theta hat bb "y"`

The random (thermal) speed is defined as

`sqrt(v_(rx)^2+v_(ry)^2+v_(rz)^2)=v_r`

and the relative values of the three components are chosen using the random number generator: To start, define:

`bb "v"_t=[v_(xt),v_(yt),v_(zt)]=2[rnd()-1/2,rnd()-1/2,rnd()-1/2]`

where each `rnd()` yields a different random number between 0 and 1. Now divide `bb "v"_t` by its magnitude to make it a unit vector:

`hat bb "v"_t=bb "v"_t/|bb "v"_t|`

finally we have the vector version of `v_r`:

`bb "v"_r=v_r hat bb "v"_t`

The initial speeds in the pipes and tapered sections are defined similarly to those in the bends with both a drift velocity along the axis and a random velocity except the drift speed is increased in the venturi and tapered sections to agree using equation 1.

Ratio of Drift Speed to the Speed of Sound

The ratio of the drift speed `v_(drift)` to the speed of sound in still (non-moving) air is called the Mach number, M. The equation for the speed of sound in an ideal gas is usually given as:

`v_(sound)=sqrt(gamma (k T_K)/m)`

Here `gamma=C_P/C_V` is the ratio of the heat capacity at constant pressure to that at constant volume, `k` is Boltzmann's constant, `T_K` is the Kelvin (absolute) temperature, and m is the mass of the particle. We have seen for all thermal problems that `k T_K` is exactly the same as the average "random" kinetic energy `k T_K=E_(avera"ge")`. So the equation becomes:

`v_(sound)=sqrt(gamma (E_(avera"ge"))/m)`

where m is the particle mass, taken to be 1 in this program. Now the "random" kinetic energy is the same as the thermal kinetic energy in still air i.e. it does not include the circulating energy. The running circulating energy is a bit smaller than the starting circulating energy because the average particle distance from the bend center increases due to the need to make circles in the bend. Our initial number, M, outside the venturi can be written as

`M=v_("drift")/sqrt(gamma (v_r^2)/2)`

For our simple spherical particles `gamma=5/3`. The reason that `gamma` depends on the form of the gas molecules is that for polyatomic molecules, their rotational energy and vibrational energies do not contribute to the gas pressure.

Explanation of the Parameter Sliders and Radio Buttons

Angle Sliders

For display only, the (rx,ry,rz) angle sliders adjust the rotation of the wind tunnel about its (x,y,z) axes. With (rx,ry,rx)=(0,0,0), the z axis is along the vertical of the screen and the y axis is along the horizontal of the screen. That leaves the x axis that points outward from the screen surface. The initial setup is with ry=90 degrees so the y axis is still horizontal and the x axis points downward with the z axis pointing outward from the screen surface.

Initial Gas Velocities

The gas drift velocity is initially along the median line of the wind tunnel. The gas thermal velocity is superimposed on the drift velocity and is entirely random in direction but initially all the particles have the same velocity magnitude (speed).

Venturi Parameters

The Throat Length is the sum of the lengths of the two tapered sections (funnels) and the venturi length. The Venturi Radius Ratio is the ratio of the venturi tube radius to that of the un-tapered sections of the wind tunnel tube.

Particle Radius

Particle Radius affects the time between collisions and therefore the time before the gas velocities become random in direction and also the time before the energies reach their expected Maxwell-Boltzmann distribution. In a gas that has an initial drift speed, I call this process "thermalization". Larger radius means sooner thermalization.

SLice Thickness

In order to display the magnitudes of the Density, Speed, and the Pressure, I chose to color sections (slices) of the wind tunnel according to the color bar that is displayed at the center of the wind tunnel. A slice thickness is chosen so that there is an integer number of equally thick slices that comprise the wind tunnel.

Plot Radio Buttons and Pause Button

As might be expected, only one of the the Density, Speed, or Pressure can be displayed at once so the Plot Buttons enable selection between these. An important note is, while the sliders always required a restart from the selected initial conditions, changing the button selection does not cause a restart since all three of the slice results are continuously updated at no significant speed loss.
The Start and Pause buttons also do not require a restart.

x angle: 0 degrees
y angle: 90 degrees
z angle: 0 degrees
Gas Drift Speed: 1.5
Gas Thermal Speed: 1.5
Venturi Radius Ratio: 0.5
Throat Length: 600
Slice Thickness: 30
Particle Radius: 2
Funnel Length: 200
Canvas 1:3D Drawing of Wind Tunnel and Particles Animation

Dynamics Inside the Wind Tunnel

The particles often collide with the inside of the tunnels. When they do they are reflected about the local normal of the inside wall of the tunnel without losing kinetic energy or anglular momentum. To find the local normal, one must use some vector algebra. For the two bends, the particle's position can be characterized by two vectors: Its azimuth , `bb "R"_theta\ \ `, with respect to bend's arc center and its altitude `bb "r"_phi`.   The actual particle position vector will be called `bb "r"_v` and it will be the sum of the two vectors:

`bb "r"_v=bb "R"_(median)+bb "r"_phi \ \ \ \ \ (4)`

where `bb "R"_(median)` is a vector that points to the median axis of the tunnel and `bb "r"_phi` is a vector that points from the median to the particle. To determine whether the particle is close enough to the inner wall to be reflected, we need the magnitude of `bb "r"_phi`. We therefore solve equation 4 for `bb "r"_phi` and obtain:

`bb "r"_phi=bb "r"_v-bb "R"_(median) \ \ \ \ (5)`

Now we need the value of vector `bb "R"_(median)`. We know that this vector meets the tunnel median line axis at some position [x,y,0]. Note that . Let `bb "r"_(xy)` be a vector that points to `[r.x,r.y,0]`. A unit vector pointing to `bb"r_(xy)` is then `bb "r"_(xy)/|r_(xy)|` Then

`bb "R"_(median) = R_theta bb "r"_(xy)/|r_(xy)| \ \ \ \ (6)`

From equations 5 and 6 one obtains `bb "r"_phi`. The unit normal `hat bb "n"` to the insides of the bends is parallel to `bb "r"_phi`.

In the pipe sections which have x offset `R_(theta)` the vector `bb "r"_(phi) ` from the median axis to the particle can be written

` bb "r"_phi =bb "r"_v-[|R_(theta)|,0,0]\ \ \ \ \ (7)`

In the tapered sections, the median is the same as in the adjacent venturi section but distance, `r_r`, from the median to the taper wall varies linearly with the `y` coordinate. The variation is expressed as

` r_t =r_V+(|y|-L_v/2)/L_T(r_T-r_V)\ \ \ \ \ (8)`

where `L_V` is the venturi length, `L_T` is the taper length, `r_T` is the radius of the remainder of the tunnel and `r_V` is the radius of the venturi. The normals of the tapered sections have combinations of components along the tunnel median (axis) as well as along the directions perpendicular to the axis. The taper normal component parallel to the median axis is related to the slope `m_T` of the taper.



` m_T=(r_T-r_V)/L_T \ \ \ \ \ (9)`

Since `m_T` is the tangent of the slope angle, the sine and cosine of the slope angle`psi` are

` cos psi=1/sqrt(1+m_T^2) \ \ \ \ \ (10)`

` sin psi=m_T/sqrt(1+m_T^2) \ \ \ \ \ (11)`

The unit normal to the taper is composed of the following vector:

` n_T=[-cos phi cos psi,sin psi, -sin phi cos psi] \ \ \ \ \ (12)`

where 'phi' the angle from the median axis to the particle

` phi=arctan (r_(v.z)/|r_v|) \ \ \ \ \ (13)`

The normals in the objects that are not tapered are perpendicular to the tunnel axis.

When the particle collides with the inside, the velocity increment:

`delta bb "v"=2 (hat bb "n" cdot bb "v") hat bb "n" \ \ \ \ \ (8)`

is added to the incident velocity in order to conserve momentum and energy.

Explanation of the Graphics

The average kinetic energy of particles is listed on Canvas 1. The angular momentum cannot be conserved because the particles linear momentum is reduced when they hit the left funnel tapered section. With this model it is difficult to see the Bernoulli principle's reduction of the pressure in the venturi. This is because the gas is becoming rapidly ds and the wind speed energy gets converted to thermal speed so the speeds in the venturi are not nearly as large as expected by equation 1. We do see a large increase in density at the right funnel because the particles are being reflected back by the taper of that funnel.