Assignment title: Information
1
CS408 Assignment 5
Due: March 26, 2017
CS 808 Students: For assignment 5, you can do either CS 408 Assignment 5 or CS 408
Assignment 6A. Assignment 6 will be to do a small project related to your research paper.
Create a program to display a 2D view of grid-based simulation of gas moving in a flat space.
See on-line notes on Gaseous Phenomena. The method for performing simulation of gaseous
phenomena is explained in the “Explanation of Gas Simulation Program”.
http://www2.cs.uregina.ca/~anima/408/Notes/Gas/index.htm
For your program, you will need to create a graphical display of a grid (large 2D array), by
displaying each location as a rectangle of some colour. The colours are used to represent density.
Although you may choose your own colouring scheme, the following way of setting the R, G, B
colour components gives a nice display where high densities correspond to bright colours.
R = density
G = density * density * 0.05
B = density * density * density * 0.0001
The size of the grid should be controlled by symbolic constants. Initially, try 200 x 200. Each cell
in the grid is drawn as a rectangle (GL_QUAD). Use a double ‘for’ loop to display 200 x 200
rectangles. It is straightforward to calculate the x and y locations of each rectangle. Use an
orthogonal projection, rather than the usual perspective projection, by replacing the
gluPerspective command with a
gluOrtho2D(0, WIDTH, 0, HEIGHT);
In your program you should choose some way of representing the current simulation state, which
has a current density for each cell and a current (x, y) velocity for each cell. You can use three
2D arrays or one 2D array of structs or classes, or any other way of structuring your data that you
want. You will also require an identical set of values for the previous state of the simulation.
To initialize the simulation, set all densities and velocities to zero (0). Every cell should have a
1/20 chance of having some nonzero amount of gas in it. If there is gas, choose a random mass
for it between 0 and 100 and random initial velocity, with x and y each between -3.0 and +3.0.
During each simulation step, you will have an old grid and a new grid. Here is a simple approach.
Initially set every cell in a new grid to zero. Go through the old grid cell by cell. For each cell,
calculate where the given amount of gas (as expressed by its density) should move according to
the current velocity recorded for the cell. Update the corresponding four cells in the new grid by
changing their densities and velocities. After updating is complete, display the new grid. Then set
the old grid to the new grid and repeat. This approach is described in the “Explanation of Gas
Simulation Program”.notes.2
For your program, increase efficiency by avoiding copying the new grid to the old grid at the end
of each frame. Instead of doing this copying, you can simply switch which grid you are using for
each purpose. One way of doing so is to have two grids in a 2-element array instead of having an
old grid and a new grid. For frame i, use grid [i % 2] as the old grid and grid [1 – i % 2] as the
new grid.
Creative Feature: Add a creative feature of your choice to your simulation at least 1/5 as large
in scope as the above features. Several ideas may result from thinking about which aspects of the
users might want to change at runtime.
.