The Shell Momentum Balance:
Step-by-Step Method
CBE 320, DJK, Fall 2013
Step 1: Draw a (crude) picture and list assumptions
This will give a differential equation for the momentum flux (the
relevant component of ij or ij , depending how much simplifying you did in the previous step).
If you havent already, make sure you have substituted for all
components of ij , and simplified as much as possible.
What velocity components are nonzero?
What spatial variables does the velocity depend on?
Choose coordinate system.
Start a list of assumptions.
Step 6: Integrate once, if possible
Step 2: Select system
This will give a constant of integration.
If you have a boundary condition for the momentum flux component ij , apply the boundary and solve for the constant.
You will not evaluate the constant here if both of your boundary
conditions are for the velocity.
Draw a second picture, showing the shell in detail.
Shell:
a box with sides k or to the velocity v.
box faces should lie on coordinate surfaces (you choose your
coordinate system to guarantee this).
Make shell thin in the direction the velocity is varying.
You will develop intuition about how to do this as you work problems.
Step 7: Insert Newtons Law of Viscosity
Use Appendix B.1 (pp. 843844) to replace the components of
ij with the appropriate expressions in terms of velocity components.
Simplify as much as possible:
Appendix B.1 is general; some velocity components and
derivatives will be zero.
Make use of your assumptions, and add to the list of assumptions if new ones occur to you.
This will lead to a differential equation for the desired velocity
component (First order if you completed step 6, second order if
you skipped step 6).
Step 3: Apply conservation of momentum (component of interest)
{ { - { {+{ {=0
Rate at which
momentum is
transported
into
the system
Rate at which
momentum is
transported
out of
the system
fij = p dij + tij + rvivj
Force of gravity
acting on
the system
add other body
forces here (e.g.,
electrical, magnetic,
etc.)
if unsteady, then
the RHS is the time
rate of change of
momentum
Step 8: Solve the differential equation
This will produce another unknown constant.
Use a second boundary condition to get this constant.
The result of this step will be the velocity profile (for example,
vz (r)).
Add arrows to second picture indicating where flux components
transport momentum in and out.
Write out the conservation of momentum equation (momentum
balance) for the system.
Step 9: Use the velocity profile to engineer
Step 4: Simplify/eliminate terms if possible
Using the velocity (or stress) profile, you can calculate a variety
of quantities of interest, e.g.,
force on a surface,
torque on an object,
volumetric flow rate,
maximum flow rate,
etc.
Insert appropriate definitions for the ij components (i.e., ij =
pij + ij + vi vj ).
Use assumptions to eliminate terms that are zero, or terms that
cancel.
Simplify here as much as possible; you will have more opportunities to simplify later.
Step 5: Let Shell thickness 0
Let the small dimension of the box 0.
1