Fluent
Fluent
Fluent
Consider fluid flowing through a circular pipe of contant cross-section. The pipe diameter
D=0.2 m and length L=8 m. The inlet velocity Vin=1 m/ s. Consider the velocity to be
constant over the inlet cross-section. The fluid exhausts into the ambient atmosphere
which is at a pressure of 1 atm. Take density ρ=1 kg/ m3 and coefficient of viscosity µ= 2
x 10-3 kg/(ms). The Reynolds number Re based on the pipe diameter is
where Vavg is the average velocity at the inlet, which is 1m/s in this case.
Solve this problem using FLUENT. Plot the centerline velocity, wall skin-friction
coefficient, and velocity profile at the outlet. Validate your results.
Note: The values used for the inlet velocity and flow properties are chosen for
convenience rather than to reflect reality. The key parameter value to focus on is the
Reynolds no.
Preliminary Analysis
We expect the viscous boundary layer to grow along the pipe starting at the inlet. It will
eventually grow to fill the pipe completely (provided that the pipe is long enough). When
this happens, the flow becomes fully-developed and there is no variation of the velocity
profile in the axial direction, x (see figure below). One can obtain a closed-form solution
to the governing equations in the fully-developed region. You should have seen this in
the Introduction to Fluid Mechanics course. We will compare the numerical results in the
fully-developed region with the corresponding analytical results. So it's a good idea for
you to go back to your textbook in the Intro course and review the fully-developed flow
analysis. What are the values of centerline velocity and friction factor you expect in the
fully-developed region based on the analytical solution? What is the solution for the
velocity profile?
We'll create the geometry and mesh in GAMBIT which is the preprocessor for FLUENT,
and then read the mesh into FLUENT and solve for the flow solution.
If you would prefer to skip the mesh generation steps, you can create a working directory
(see below), download the mesh from here (right click and save as pipe.msh) into the
working directory and go straight to step 4.
In order to create the rectangle, we will first create the vertices at the four corners. We'll
then join adjacent vertices by straight lines to form the "edges" of the rectangle. Lastly,
we'll create a "face" corresponding to the area enclosed by the edges. In Step 2, we'll
mesh the face i.e. the rectangle. Note that in 3D problems, you'll have to form a "volume"
from faces. So the hierarchy of geometric objects in GAMBIT is vertices -> edges ->
faces -> volumes.
Create a folder called pipe in a convenient location. We'll use this as the working folder
in which files created during the session will be stored.
Note for ACCEL computer lab users: Each user gets his/her own 100 MB of disk
space under S: at ACCEL. You can put your files in S: and it'll be accessible from any
computer. This is where you should put files that you want to keep and access later on.
Start GAMBIT
Start > Programs > Lab Apps > Fluent Inc Products > Gambit 2.3.16 > Gambit 2.3.16
This brings up the GAMBIT startup window. Click Browse and select the folder that you
just created. Enter -id pipe in the options box to tell GAMBIT to use pipe as the default
file prefix, then click Run.
In Windows, the Exceed X-server starts up before the GAMBIT interface comes up.
Exceed is a third-party application needed to render the interface in Windows (GAMBIT
was originally developed under Unix). To make best use of screen real estate, move the
windows and resize them so that you approximate this screen arrangement. This way you
can read instructions in the browser window and implement them in GAMBIT.
You can resize the text in the browser window to your taste and comfort:
In Internet Explorer: Menubar > View > Text Size, then choose the appropriate font size.
In Netscape: Menubar > View > Increase Font or Menubar > View > Decrease Font.
Note that the job name pipe appears after ID: in the title bar of the Utility Menu.
• Operation Toolpad:
We'll more or less work our way across the Operation Toolpad as we go through
the solution steps. Notice that as each of the top buttons is selected, a different
"sub-pad" appears. The Geometry sub-pad is shown in the above snaphot.
• Global Control Toolpad:
The Global Control Toolpad has options such as Fit to Screen and Undo
that are very handy during the course of geometry and mesh creation.
• GAMBIT Graphics:
This is the window where the graphical results of operations are displayed.
• GAMBIT Description Panel:
The Description Panel contains descriptions of buttons or objects that the mouse
is pointing to. Move your mouse over some buttons and notice the corresponding
text in the Description Panel.
• GAMBIT Transcript Window:
This is the window to which output from GAMBIT commands is written and
which provides feedback on the actions taken by GAMBIT as you perform
operations. If, at some point, you are not sure you clicked the right button or
entered a value correctly, this is where to look to figure out what you just did.
You can click on the arrow button in the upper right hand corner to make the
Transcript window full-sized. You can click on the arrow again to return the
window to its original size. Go ahead, give this a try.
Select Solver
If the window titlebar does not say the solver is FLUENT 5/6, then you need to specify:
Verify this has been done by looking in the Transcript Window where you should see:
The boundary types that you'll be able to select in the third step depends on the solver
selected.
We can assume that the flow is axisymmetric. The problem domain is:
We will put the origin of the coordinate system at the lower left corner of the rectangle.
The coordinates of the corners are shown in the figure below:
We will first create four vertices at the four corners and join adjacent vertices to get the
edges of the rectangle. We will then form a face that covers the area of the rectangle.
Create Vertices
Find the buttons described below by pointing the mouse at each of the buttons and
reading the Description Window.
Operation Toolpad > Geometry Command Button > Vertex Command Button >
Create Vertex
Notice that the Create Vertex button has already been selected by default. After you
select a button under a sub-pad, it becomes the default when you go to a different sub-pad
and then come back to the sub-pad.
In the Transcript window, GAMBIT reports that it "Created vertex: vertex.1". The
vertices are numbered vertex.1, vertex.2 etc. in the order in which they are created.
Vertex 2: (0,0.1,0)
Vertex 3: (8,0.1,0)
Vertex 4: (8,0,0)
Note that for a 2D problem, the z-coordinate can always be left to the default value of 0.
Operation Toolpad > Global Control > Fit to Window Button
This fits the four vertices of the rectangle we have created to the size of the Graphics
Window.
Create Edges
We'll now connect appropriate pairs of vertices to form edges. To select any entity in
GAMBIT, hold down the Shift key and click on the entity.
Operation Toolpad > Geometry Command Button > Edge Command Button >
Create Edge
Select two vertices that make up an edge of this rectangle by holding down the Shift
button and clicking on the corresponding vertices. As each vertex is picked, it will appear
red in the Graphics Window. Then let go of the Shift button. We can check the selected
vertices by clicking on the up-arrow next to Vertices:.
This will bring up a window containing the vertices that have been selected. Vertices can
be moved from the Available and Picked lists by selecting them and then pressing the left
or right arrow buttons.
After the correct vertices have been selected, click Close, then click Apply in the Create
Straight Edge window.
Create Face
Operation Toolpad > Geometry Command Button > Face Command Button >
Form Face
To form a face out of the area enclosed by the four lines, we need to select the four ledges
that enclose this area. This can be done by holding down the Shift key, clicking on each
line (notice that the currently selected line appears red), and then releasing the Shift key
after all four lines have been selected.
Alternatively, an easier way to do this would be to click on the up arrow next to edges:
This will bring up the Edge List window. Click on All-> to select all of the edges at once.
Click Close.
We'll now create a mesh on the rectangular face with 100 divisions in the axial direction
and 5 divisions in the radial direction. We'll first mesh the four edges and then the face.
The desired grid spacing is specified through the edge mesh.
Mesh Edges
Operation Toolpad > Mesh Command Button > Edge Command Button > Mesh
Edges
Shift-click or bring up the Edge List window and select both the vertical lines. If this is
difficult, one can zoom in on an edge by holding down the Ctrl button, clicking and
dragging the mouse to specify an area to zoom in on, and releasing the Ctrl button. To
return to the main view, click on the Global Control Toolpad > Fit to Window Button
again. You can also hold down Ctrl and double-click in the window to zoom out to a
fitting window. To pan the view, hold down the middle mouse button and drag the
mouse.
Once a vertical edge has been selected, select Interval Count from the drop down box that
says Interval Size in the Mesh Edges Window. Then, in the box to the left of this combo
box, enter 5 for the interval count.
Click Apply. Nodes appear on the edges showing that they are divided into 5.
Repeat the same process for the horizontal edges, but with an interval count of 100.
Now that the edges are meshed, we are ready to create a 2-D mesh for the face.
Mesh Face
Operation Toolpad > Mesh Command Button > Face Command Button > Mesh
Faces
Shift left-click on the face or use the up arrow next to Faces to select the face. Click
Apply.
We'll next set the boundary types in GAMBIT. The left edge is the inlet of the pipe, the
right edge the outlet, the top edge the wall, and the bottom edge the axis.
Operation Toolpad > Zones Command Button > Specify Boundary Types Command
Button
This will bring up the Specify Boundary Types window on the Operation Panel. We will
first specify that the left edge is the inlet. Under Entity:, pick Edges so that GAMBIT
knows we want to pick an edge (face is default).
Now select the left edge by Shift-clicking on it. The selected edge should appear in the
yellow box next to the Edges box you just worked with as well as the Label/Type list
right under the Edges box.
For Type:, select VELOCITY_INLET. (Note: Sometimes all the items in a dropdown
menu will not be visible. If you cannot find the VELOCITY_INLET option in the Type
menu, try maximizing the window. If it is still not visible, try auto-hiding your taskbar.
Right-click on the taskbar and go to properties.)
Click Apply. You should see the new entry appear under Name/Type box near the top of
the window.
Repeat this process for the other three edges according to the following table:
Edge
Name Type
Position
Left inlet VELOCITY_INLET
Right outlet PRESSURE_OUTLET
Top wall WALL
Bottom centerline AXIS
You should have the following edges in the Name/Type list when finished:
Save and Export
Type in pipe.msh for the File Name:. Select Export 2d Mesh since this is a 2 dimensional
mesh. Click Accept.
Check pipe.msh has been created in your working directory (the box will be filled in red).
The "2ddp" option is used to select the 2-dimensional, double-precision solver. In the
double-precision solver, each floating point number is represented using 64 bits in
contrast to the single-precision solver which uses 32 bits. The extra bits increase not only
the precision but also the range of magnitudes that can be represented. The downside of
using double precision is that it requires more memory.
Import Grid
Navigate to the working directory and select the pipe.msh file. This is the mesh file that
was created using the preprocessor GAMBIT in the previous step. FLUENT reports the
mesh statistics as it reads in the mesh:
Check the number of nodes, faces (of different types) and cells. There are 500
quadrilateral cells in this case. This is what we expect since we used 5 divisions in the
radial direction and 100 divisions in the axial direction while generating the grid. So the
total number of cells is 5*100 = 500.
Also, take a look under zones. We can see the four zones inlet, outlet, wall, and centerline
that we defined in GAMBIT.
First, we check the grid to make sure that there are no errors.
Any errors in the grid would be reported at this time. Check the output and make sure that
there are no errors reported. Check the grid size:
Translation: The grid can be translated in any direction by holding down the Left Mouse
Button and then moving the mouse in the desired direction.
Zoom In: Hold down the Middle Mouse Button and drag a box from the Upper Left Hand
Corner to the Lower Right Hand Corner over the area you want to zoom in on.
Zoom Out: Hold down the Middle Mouse Button and drag a box anywhere from the
Lower Right Hand Corner to the Upper Left Hand Corner.
Use these operations to zoom into the grid to obtain the view shown below.
Note: The zooming operations cannot be performed without a middle mouse button.
You can also look at specific parts of the grid by choosing the boundaries you wish to
view under Surfaces (click to select and click again to deselect a specific boundary).
Click Display again when you have selected your boundaries. For example, the wall,
outlet, and centerline boundaries have been selected in the following view:
These options will display the graph:
For convenience, the button next to Surfaces selects all of the boundaries while the
deselects all of the boundaries at once.
Choose Axisymmetric under Space. We'll use the defaults of pressure based
("segregated", in older versions) solver, implicit formulation, steady flow and absolute
velocity formulation. Click OK.
Laminar flow is the default. So we don't need to change anything in this menu. Click
Cancel.
For incompressible flow, the energy equation is decoupled from the continuity and
momentum equations. We need to solve the energy equation only if we are interested in
determining the temperature distribution. We will not deal with temperature in this
example. So leave the Energy Equation unselected and click Cancel to exit the menu.
Change Density to 1.0 and Viscosity to 2e-3. These are the values that we specified under
Problem Specification. We'll take both as constant.
Click Change/Create. Close the window.
For all flows, FLUENT uses gauge pressure internally. Any time an absolute pressure is
needed, it is generated by adding the operating pressure to the gauge pressure. We'll use
the default value of 1 atm (101,325 Pa) as the Operating Pressure.
We'll now set the value of the velocity at the inlet and pressure at the outlet.
We note here that the four types of boundaries we defined are specified as zones on the
left side of the Boundary Conditions Window. The centerline zone should be selected by
default. Make sure it is, then make sure the Type of this boundary is selected as axis and
click Set.... Notice that there is nothing to set for the axis. Click OK.
Move down the list and select inlet under Zone. Note that FLUENT indicates that the
Type of this boundary is velocity-inlet. Recall that the boundary type for the "inlet" was
set in GAMBIT. If necessary, we can change the boundary type set previously in
GAMBIT in this menu by selecting a different type from the list on the right.
Click on Set.... Enter 1 for Velocity Magnitude. Click OK. This sets the velocity of the
fluid entering at the left boundary.
The (absolute) pressure at the outlet is 1 atm. Since the operating pressure is set to 1 atm,
the outlet gauge pressure = outlet absolute pressure - operating pressure = 0. Choose
outlet under Zone. The Type of this boundary is pressure-outlet. Click on Set.... The
default value of the Gauge Pressure is 0. Click Cancel to leave the default in place.
Lastly, click on wall under Zones and make sure Type is set as wall. Click on each of the
tabs and note that only momentum can be changed under the current conditions. This will
not be so under later exercises so make a note of the location of these options. Click OK.
Go to Step 5: Solve!
Step 5: Solve!
In the Solution Initialization menu that comes up, choose inlet under Compute From. The
Axial Velocity for all cells will be set to 1 m/s, the Radial Velocity to 0 m/s and the
Gauge Pressure to 0 Pa. These values have been taken from the inlet boundary condition.
Change the residual under Convergence Criterion for continuity, x-velocity, and y-
velocity, all to 1e-6.
Also, under Options, select Plot. This will plot the residuals in the graphics window as
they are calculated.
Click OK.
Type in pipe.cas for Case File. Click OK. Check that the file has been created in your
working directory. If you exit FLUENT now, you can retrieve all your work at any time
by reading in this case file.
The residuals for each iteration is printed out as well as plotted in the graphics window as
they are calculated.
The residuals fall below the specified convergence criterion of 1e-6 in about 46 iterations.
Actual number of convergence steps may vary slightly.
Enter pipe.dat for Data File and click OK. Check that the file has been created in your
working directory. You can retrieve the current solution from this data file at any time.
Centerline Velocity
We'll plot the variation of the axial velocity along the centerline.
Under Y Axis Function, pick Velocity... and then in the box under that, pick Axial
Velocity.
Please note that X Axis Function and Y Axis Function describe the x and y axes of the
graph, which should not be confused with the x and y directions of the pipe.
Finally, select centerline under Surfaces since we are plotting the axial velocity along the
centerline. This finishes setting up the plotting parameters.
Click Plot.
This brings up a plot of the axial velocity as a function of the distance along the
centerline of the pipe.
(Click picture for larger image)
In the graph that comes up, we can see that the velocity reaches a constant value beyond a
certain distance from the inlet. This is the fully-developed flow region.
Change the axes extents: In the Solution XY Plot window, click on Axes.... Under
Options, deselect Auto Range. The boxes under Range should now be activated. Select X
under Axis. Enter 1 for Minimum and 3 for Maximum under Range.
We'll turn on the grid lines to help estimate where the flow becomes fully developed.
Check the boxes next to Major Rules and Minor Rules under Options. Click Apply.
Now, pick Y under Axis and once again deselect Auto Range under Options, then enter
1.8 for Minimum and 2.0 for Maximum under Range. Also select Major Rules and Minor
Rules to turn on the grid lines in the Y direction. We have now finished specifying the
range for each axes, so click Apply and then Close.
Go back to the Solution XY Plot menu and click Plot to replot the graph with the new
axes extents. We can see that the fully-developed region starts at around x=3m and the
centerline velocity in this region is 1.93 m/s.
In the Solution XY Plot Window, check the Write to File box under Options. The Plot
button should have changed to Write.... Click on Write.... Enter vel.xy as the XY File
Name and click OK. Check that this file has been created in your FLUENT working
directory.
Leave the Solution XY Plot Window and the Graphics Window open and in the main
FLUENT window click on:
EPS - if you have a postscript viewer, this is the best choice. EPS allows you to save the
file in vector mode, which will offer the best viewable image quality. After selecting
EPS, choose Vector from under File Type.
TIFF - this will offer a high resolution image of your graph. However, the image file
generated will be rather large, so this is not recommended if you do not have a lot of
room on your storage device.
JPG - this is small in size and viewable from all browsers. However, the quality of the
image is not particularly good.
After selecting your desired image format and associated options, click on Save...
Enter vel.eps, vel.tif, or vel.jpg depending on your format choice and click OK.
Verify that the image file has been created in your working directory. You can now copy
this file onto a disk or print it out for your records.
FLUENT provides a large amount of useful information in the online help that comes
with the software. Let's probe the online help for information on calculating the
coefficient of skin friction.
Click on S in the links on top and scroll down to skin friction coefficient. Click on the
second 965 link (normally, you would have to go through each of the links until you find
what you are looking for). We can see an excerpt on the skin coefficient as well as the
equation for calculating it.
Click on the link for Reference Values panel, which tells us how to set the reference
values used in calculating the skin coefficient.
Select inlet under Compute From to tell FLUENT to calculate the reference values from
the values at inlet. Check that density is 1 kg/m3 and velocity is 1 m/s. (Alternately, you
could have just typed in the appropriate values). Click OK.
Go back to the Solution XY Plot menu. Uncheck Write to File under Options since we
want to plot to the window right now. We can leave the other Options and Plot Direction
as is since we are still plotting against the x distance along the pipe.
Under the Y Axis Function, pick Wall Fluxes..., and then Skin Friction Coefficient in the
box under that.
Reset axes ranges: Go to Axes... and re-select Auto-Range for the Y axis. Click Apply.
Set the range of the X axis from 1 to 8 by selecting X under Axis, entering 1 under
Minimum, and 8 under Maximum in the Range box (remember to de-select Auto-Range
first if it is checked).
Click Apply, Close, and then Plot in the Solution XY Plot Window.
(Click picture for larger image)
We can see that the fully developed region is reached at around x=3.0m and the skin
friction coefficient in this region is around 1.54. Compare the numerical value of 1.54
with the theoretical, fully-developed value of 0.16.
Save the data from this plot: Pick Write to File under Options and click Write.... Enter
cf.xy for XY File and click OK.
Velocity Profile
We'll next plot the velocity at the outlet as a function of the distance from the center of
the pipe. To do this, we have to set the y axis of the graph to be the y axis of the pipe (the
radial direction).
To plot the position variable on the y axis of the graph, uncheck Position on X Axis under
Options and choose Position on Y Axis instead. To make the position variable the radial
distance from the centerline, under Plot Direction, change X to 0 and Y to 1. To plot the
axial velocity on the x axis of the graph, for X Axis Function, pick Velocity... and Axial
Velocity under that.
Since we want to plot this at the outlet boundary, pick outlet under Surfaces.
Change both the x and y axes to Auto-Range. (Don't forget to click apply before selecting
a different axis)
Uncheck Write to File under Options so that we can see the graph. Click Plot.
(Click picture for larger image)
Save the data from this plot: Pick Write to File under Options and click Write.... Enter
profile.xy for XY File and click OK.
To see how the velocity profile changes in the developing region, let us add the profiles
at x=0.6m (x/D=3) and x=0.12m (x/D=6) to the above plot. First, create a line at x=0.6m
using the Line/Rake tool:
We'll create a straight line from (x0,y0)=(0.6,0) to (x1,y1)=(0.6,0.1). Select Line Tool
under Options. Enter x0=0.6, y0=0, x1=0.6, y1=0.1. Enter line1 under New Surface
Name. Click Create.
To see the line just created, select
Note that line1 appears in the list of surfaces. Select all surfaces except default-interior.
Click Display. This displays all surfaces but not the mesh cells. Zoom into the region
near the inlet to see the line created at x=0.6m. (Click here to review the zoom
functionality discussion in step 4.) line1 is the white vertical line to the right in the figure
below.
Now we can plot the velocity profiles at x=0.6m (x/D=3) and x=0.12m (x/D=6) along
with the outlet profile. In the Solution XY plot menu, use the same settings as above.
Under Surfaces, in addition to outlet, select line1 and line2. Make sure Node Values is
selected under Options. Click Plot. Your symbols might be different from the ones below.
You can change the symbols and line styles under the Curves... button. Click on Help in
the Curves menu if you have problems figuring out how to change these settings.
The profile three diameters downstream is fairly close to the fully-developed profile at
the outlet. If you redo this plot using the fine grid results in the next step, you'll see that
this is not actually the case. The coarse grid used here doesn't capture the boundary layer
development properly and underpredicts the development length.
In FLUENT, you can choose to display the computed cell-center values or values that
have been interpolated to the nodes. By default, the Node Values option is turned on, and
the interpolated values are displayed. Node-averaged data curves may be somewhat
smoother than curves for cell values.
Velocity Vectors
One can plot vectors in the entire domain, or on selected surfaces. Let us plot the velocity
vectors for the entire domain to see how the flow develops downstream of the inlet.
Zoom into the region near the inlet. (Click here to review the zoom functionality
discussion in step 4.) The length and color of the arrows represent the velocity
magnitude. The vector display is more intelligible if one makes the arrows shorter as
follows: Change Scale to 0.4 in the Vectors menu and click Display.
You can reflect the plot about the axis to get an expanded sectional view:
Under Mirror Planes, only the axis surface is listed since that is the only symmetry
boundary in the present case. Select axis and click Apply. Close the Views window.
The velocity vectors provide a picture of how the flow develops downstream of the inlet.
As the boundary layer grows, the flow near the wall is retarded by viscous friction. Note
the sloping arrows in the near wall region close to the inlet. This indicates that the
slowing of the flow in the near-wall region results in an injection of fluid into the region
away from the wall to satisfy mass conservation. Thus, the velocity outside the boundary
layer increases.
By default, one vector is drawn at the center of each cell. This can be seen by turning on
the grid in the vector plot: Select Draw Grid in the Vectors menu and then click Display
in the Grid Display as well as the Vectors menus. Velocity vectors are the default, but
you can also plot other vector quantities. See section 27.1.3 of the user manual for more
details about the vector plot functionality.
It is very important to assess the dependence of your results on the mesh used by
repeating the same calculation on different meshes and comparing the results. We will re-
do the previous calculation on a 100x10 mesh and compare the results with the 100x5
mesh used previously. If you prefer to skip the GAMBIT steps for modifying the mesh,
download the 100x10 mesh (by right-clicking on the link) and go directly to the FLUENT
analysis discussed below.
The 100x5 mesh is saved as pipe.dbs in your working folder. Copy and paste the file in
the same folder. Rename Copy of pipe.dbs to pipe2.dbs. We will work with pipe2.dbs in
order to retain pipe.dbs as is. Launch GAMBIT and browse to where pipe2.dbs is saved.
Notice that under Session ID, pipe2 is now listed. Select this and click Run. Note in the
main menu bar that pipe2 is the ID of this job. Files created during this session will have
that prefix.
We will delete the face mesh, modify the edge meshes for the vertical edges and remesh
the face. To delete the original face mesh, choose
Operation Toolpad > Mesh Command Button > Face Command Button > Delete Face
Meshes
In the Delete Face Meshes Window that comes up, uncheck the Remove unused lower
mesh box. This tells GAMBIT to remove the face mesh only and keep the edge meshes
associated with the face mesh. Since we will be changing the mesh on only two edges of
the rectangle, there is no need to redo the meshes for all four edges.
Select the only face of the rectangle by shift-clicking on it and then click Apply.
To change the number of divisions on the vertical edges from 5 to 10, choose:
Operation Toolpad > Mesh Command Button > Edge Command Button > Mesh Edges
Select the two vertical edges by holding down the Shift button, clicking on each in turn,
and then releasing the Shift button. Select Interval count from the box under Spacing that
says Interval size. Change the number in the box next to the Interval count box from 5 to
10.
Make sure that the Remove old mesh box is checked under Options. This will make sure
that the old edge meshes are erased before the new edge meshes are created.
Click Apply.
Remember that you can zoom in by holding down Ctrl, dragging a box across the area
you want to zoom in on, and then releasing Ctrl. Do this now and make sure that the
vertical edges have 10 divisions.
Operation Toolpad > Mesh Command Button > Face Command Button > Mesh Faces
Shift-click on the face in the Graphics Window to select it. Click Apply.
Type in pipe2.msh for the File Name:. Select Export 2d Mesh option. Click Accept.
Repeat steps 4 and 5 of this tutorial with the 100x10 mesh (a tad on the repetitious side
but consider it good practice).
One you obtain the solution, plot the variation of the centerline velocity along the x-
direction as described in step 6. Compare this result with that obtained on the previous
mesh which is stored in the vel.xy file created earlier. To do this, after centerline velocity
has been plotted, click on Load File... in the Solution XY Plot window. Navigate to your
working folder if necessary and click on vel.xy and OK. Click Plot.
In the graphics window, we can see both of the lines plotted in the same window. Adjust
the axes so that you can zoom in on the beginning of the fully developed region.
In the centerline velocity plot above, the white and red symbols represent the results on
the 100x10 mesh and 100x5 meshes, respectively. The centerline velocity in the fully-
developed region for the finer mesh is 1.98 m/s. This value agrees better with the
analytical value of 2 m/s that the value of 1.93 m/s obtained on the coarser mesh. Save
the data for this plot as vel2.xy. The velocity result gets more accurate on refining the
mesh as expected.
Plot the skin friction coefficient as described in step 6. Compare the result with that
obtained on the 100x5 mesh by loading it from cf.xy.
(Click here for larger image)
The finer mesh provides a skin friction coefficient of 0.159 in the fully-developed region,
which is much closer to the theoretical value of 0.16 than the corresponding coarser mesh
value of 0.154. Save the data for this plot as cf2.xy.
Similarly, plot the velcoity profile at the outlet and compare with the coarser grid result in
out.xy. The two results compare well with the greatest deviation occurring near the
centerline. Save the data for this plot as out2.xy.
If you repeat the calculation on a 100x20 mesh, you'll see that the results on the two
finest meshes are grid-independent to a high level of accuracy. In the plots below, the
white, red and green symbols correspond to the 100x20, 100x10 and 100x5 meshes,
respectively.
Velocity along centerline:
Skin Coefficient:
Go to Problem 1
3 Fuel 22.81
4 Organics 45.42
5 Inert 30.22
6 Recyclable 1.32
7 Other 0.19
Mineral matter 15 % - 25 %
Carbon 35 % - 40 %
Hydrogen 5%-8%
Nitrogen 1 % - 1.5 %
Sulphur 0.2 % - 0.5 %
Oxygen 25 % - 30 %
Gross Calorific Value of 2,600 kcal / kg
RDF (Avg)
Pyrolysis
Conventional gasification
The char that would have been left behind in pyrolysis is converted to
additional syngas. Inorganic materials end up as either bottom ash (low-
temperature gasification) or as an inert, glassy slag (high-temperature
gasification). Gasification has been in use for more than 200 years. It has
been used worldwide to gasify coal to produce `town gas¿ for heating,
cooking and lighting. MSW gasification has been implemented worldwide in
the past few years, with some facilities producing steam and others
electricity.
Plasma gasification
ENVIRONMENTAL BENEFITS
The fluff is expected to have a gross calorific value (GCV) of 2600 kcal/kg to
3,000 Kcal/kg.
THE CONVERSION OF MSW INTO RDF
The conversion of MSW into RDF involves the following operations
M anual Segregation.
Shredding.
Screening to separate both, fine inerts and some percentage of bio-degradable
matter.
Rotary conveying and as per requirements of drying system.
Fines Screening.
Density Separator (Ballistic separation)
The schematic process flow diagram for conversion of MSW to RDF is given in
Figure
Bhide A D and Sunderesan B B. (1983) Processing method for the future solid waste management
in developing countries. pp. 124–134, Indian National Science Documentation, New Delhi
T Rand J. Haukohl U AMarxen WORLD BANK TECHNICAL PAPER NO. 462 June 2000 Municipal Solid Waste Incineration
Requirements for a Successful Project
Stephen D. Jenkins Conversion technologies NOV-2005new alternative for MSW management Waste management
world
Buekens A., Huang, H. (1998), “Comparative evaluation of techniques for controlling the
formation and emission of chlorinated dioxins/furans in municipal solid waste incineration”. Journal of Hazardous
Materials.Columbia University (2001).
Life After Fresh Kills: Moving Beyond New York City’s Current
Waste Management Plan. Policy, Technical and Environmental Considerations.
Denison, R.A. (1996) “Environmental life-cycle comparisons of recycling, landfilling, and Incineration” Annual Reviews of Energy
and Environment, Vol. 21 191-237.
Elliott, D.C. (1988). “Relation of Reaction Time and Temperature to Chemical Composition of Pyrolysis Oils,” ACS Symposium
Series 376, Pyrolysis Oils from Biomass.
ETIS Web site (1999). “Characterization of Alkali Metal Emission from Fuels and Samples Collected from Fluidized Bed
Gasification.” Paper No. FBC99-0125
Evans, R.J., Milne, T.A (1997). “Chemistry of Tar Formation and Maturation in the
Thermochemical Conversion of Biomass.” Developments in Thermochemical Biomass
Conversion, Vol. 2.
Feldmann, H.F. (1973) “A process for converting solid wastes to pipeline gas.” US Patent
3,3773,187
Korsgren, J.G. et al. (1999). “Characterization of Alkali Metal emission from fuels and samples collected from Fluidized bed
gasification.” Proceedings of the15th annual conference on fluidized bed combustion. Savannah, GA.
Krigmont, H (1999). “IBGCC power generation concept: A gateway for a cleaner future. Allied
Milne, T.A., Evans, R.J., (1998). “Biomass Gasifier Tars: Their Nature, Formation, and
Conversions.” National Renewable Energy Laboratory, Golden, CO. NREL/TP-570-25357.
Miranda, M., Hale, B. (1999). “Re-Covering All the Bases: A Comparison of Landfills andResource Recovery Facilities in Puerto
Rico.” Nickolas School of the Environment, DukeUniversity.
Morris, M. (1998). “Electricity Production from Solid Waste Fuels Using advanced Gasification Technologies.” TPS Termiska
Processer AB. www.tps.se
Niessen, W., Markes., C et al. (1996). “Evaluation of Gasification and Novel Thermal
Processes for the treatment of Municipal Solid Waste.” National Renewable Energy Laboratory,
Schultz, H. (2000). “A Bright Future for Renewable Energy Based on Waste.” International Directory of Waste Management.
http://www.jxj.com/yearbook/iswa/2000/bright_future_schultz.html
Scottish Agricultural Website (2002)
http://www.sac.ac.uk/envsci/External/WillowPower/Conversn.htm
SCS Engineers. (1992) New York City Waste Composition Study,1989-1990
Spliethoff, H (2001). “Status of Biomass Gasification for Power Production.” IFRF Combustion Journal. Article Number 200109.
Tchobanoglous, G., Theisen, H., and Vigil, S. (1993) “Integrated Solid Waste Management” Chapter 4, McGraw-Hill, New York.
Themelis, N.J., Kim., Y.H., Brady M.H. (2002). “Energy Recovery from New York City Solid Wastes.” IWSA
Themelis, N.J., Gregory, A.F. (2001). “Sources and Material Balance of Mercury in the New York-New Jersey Harbor.” Report to the
NY Academy of Sciences.
Uil, H., Mozafarrian, M., et. al. (2000) “New and Advanced Processes for Biomass
Gasification.” Netherlands Energy Research Foundation (ECN).
Investigation Into Municipal Solid Waste Gasification – PUB Draft May 27, 2004
Prepared for: Alameda Power & Telecom Draft to the Public Utilities Board Prepared by: Advanced Energy Strategies, Inc.
May 27, 2004
Dear Sir,
Kindly get two sets of question paper prepared from two different
experts. You can share a copy of letter and annexure provided to
you.Finally
01 set need to be moderated by you. After completion, please call me
so that I will depute my person to pick up all the three papers. Please
pack all papers separately. Kindly get the keys signed by the experts.
Please ensure to send me back the Annexure-III.
Regards.