Difference between revisions of "Cs382"

From Earlham CS Department
Jump to navigation Jump to search
(Modified Images)
 
(14 intermediate revisions by 2 users not shown)
Line 2: Line 2:
 
  | __TOC__
 
  | __TOC__
 
  |}
 
  |}
This page documents the work of CS382 - Scientific Computing, Fall 2007
+
 
----
 
 
==Using the enVision Tabletop Groundwater Simulator==
 
==Using the enVision Tabletop Groundwater Simulator==
  
Line 21: Line 20:
 
[[Image:parabolic.jpg|thumb|right|Parabolic contaminant flow model]]
 
[[Image:parabolic.jpg|thumb|right|Parabolic contaminant flow model]]
  
*Experiments
+
Experiments
**Demonstrating porosity
+
*Demonstrating porosity
***model water flow unconfined aquifier
+
**model water flow unconfined aquifier
**Illustrating groundwater flow in a confined aquifer
+
*Illustrating groundwater flow in a confined aquifer
***We will use a cellular automata model where at the lowest level, a cell is either fresh water or contaminated. We see this problem split into two concepts - speed and direction.  
+
**We will use a cellular automata model where at the lowest level, a cell is either fresh water or contaminated. We see this problem split into two concepts - speed and direction.  
****Direction: The illustration to the right demonstrates our assumptions about how the water will move through the material. The simulation will calculate a new direction at each generation based on it's position relative to the known locations of water input and output.
+
***Direction: The illustration to the right demonstrates our assumptions about how the water will move through the material. The simulation will calculate a new direction at each generation based on it's position relative to the known locations of water input and output.
****Speed: Remains constant throughout generations for a given run. The "speed" value represents a combination of speed of water flow and material porosity, and in terms of the simulation is the possibility that a a neighboring cell in the flow direction becomes contaminated.
+
***Speed: Remains constant throughout generations for a given run. The "speed" value represents a combination of speed of water flow and material porosity, and in terms of the simulation is the possibility that a a neighboring cell in the flow direction becomes contaminated.
 
+
*Describing recharge, transition and discharge areas
**Describing recharge, transition and discharge areas
+
**modeling behavior of water recharge, discharge in wells, lake, etc
***modeling behavior of water recharge, discharge in wells, lake, etc
 
  
*Computational Tools
+
Computational Tools
**C
+
*C
***+Very fast
+
**+Very fast
***+Libraries are available
+
**+Libraries are available
***+Good distributed Libraries
+
**+Good distributed Libraries
***-Potentially difficult to use
+
**-Potentially difficult to use
***-no graphics libraries
+
**-no graphics libraries
**Netlogo
+
*Netlogo
***+Fancy Graphics
+
**+Fancy Graphics
***+Fun to use
+
**+Fun to use
***+Available examples/code
+
**+Available examples/code
***-Slow
+
**-Slow
***-Small problem size
+
**-Small problem size
***-No Distributed processing
+
**-No Distributed processing
  
 
=== Peter and Mikio ===
 
=== Peter and Mikio ===
* Experiment
+
Experiment
** Describing the model
+
* Describing the model
*** Describing the various parts of the Groundwater Simulator by attaching tags: Key words -- wells, artesian wells, lake, underground storage tank, septic tank, springs, vegetative layer, river/ocean, recharge area, discharge area, aquifers, confining layer, clay layers
+
** Describing the various parts of the Groundwater Simulator by attaching tags: Key words -- wells, artesian wells, lake, underground storage tank, septic tank, springs, vegetative layer, river/ocean, recharge area, discharge area, aquifers, confining layer, clay layers
** Illustrating and Calculating Porasity of different types of earth materials
+
* Illustrating and Calculating Porasity of different types of earth materials
** Determining how it is easy for ground water to move in different earth materials.
+
* Determining how it is easy for ground water to move in different earth materials.
* Computetional Tool
+
Computetional Tool
** NetLogo for computatinal experiment
+
* NetLogo for computatinal experiment
  
 
=== Brad and Nate ===
 
=== Brad and Nate ===
 
 
Our goal is an incremental approach towards illustrating groundwater contamination in a confined aquifer.  The confined aquifer, viewed between wells 1 and 8, offers an environment within the groundwater simulator with the fewest variables.  The first 4 experiments are an effort to illustrate the behavior and underlying science that must be understood and demonstrated in the final experiment.
 
Our goal is an incremental approach towards illustrating groundwater contamination in a confined aquifer.  The confined aquifer, viewed between wells 1 and 8, offers an environment within the groundwater simulator with the fewest variables.  The first 4 experiments are an effort to illustrate the behavior and underlying science that must be understood and demonstrated in the final experiment.
  
*Experiments
+
Experiments
** [[Cs382/Diffusion_Experiment|Diffusion]]
+
* [[Cs382/Diffusion_Experiment|Diffusion]]
*** Show diffusion without groundwater movement.
+
** Show diffusion without groundwater movement.
** [[Cs382/Flow_Rate_Experiment|Flow Rate]]
+
* [[Cs382/Flow_Rate_Experiment|Flow Rate]]
*** Show the leading edge of groundwater contamination as a indicator of flow rate (related to section 5 and 13 in manual)  
+
** Show the leading edge of groundwater contamination as a indicator of flow rate (related to section 5 and 13 in manual)  
** [[Cs382/Plume_Length_Experiment|Contaminant Plume Length]]
+
* [[Cs382/Plume_Length_Experiment|Contaminant Plume Length]]
*** Determine whether contaminant plume length is affected by flow rate for a given amount of dye
+
** Determine whether contaminant plume length is affected by flow rate for a given amount of dye
** [[Cs382/Soil_Density_Experiment|Soil Density]]
+
* [[Cs382/Soil_Density_Experiment|Soil Density]]
*** Use displacement method and measurements of aquifer component to determine the density of the soil. We can use this value in silico.
+
** Use displacement method and measurements of aquifer component to determine the density of the soil. We can use this value in silico.
** Illustrate laminar flow in a confined aquifer (Activity 7-1)
+
* Illustrate laminar flow in a confined aquifer (Activity 7-1)
*** Show laminar flow between wells 1 and 8.
+
** Show laminar flow between wells 1 and 8.
  
* Computational Tools
+
Computational Tools
** NetLogo
+
* NetLogo
*** + Visualization built in
+
** + Visualization built in
*** + Agent and cell based simulation structure built in
+
** + Agent and cell based simulation structure built in
*** - Possible limitation on world size / agent count in RAM
+
** - Possible limitation on world size / agent count in RAM
*** - Possible run time slower than groundwater simulator at higher flow rates
+
** - Possible run time slower than groundwater simulator at higher flow rates
*** - Not parallel
+
** - Not parallel
** Python and MYMPI
+
* Python and MYMPI
*** + Parallelizable
+
** + Parallelizable
*** + Faster than NetLogo in serial code ?
+
** + Faster than NetLogo in serial code ?
*** + Visualization software exists
+
** + Visualization software exists
**** TKInter - easy to install; seemingly easy to use
+
* TKInter - easy to install; seemingly easy to use
*** - Visualization software must be integrated
+
** - Visualization software must be integrated
*** - MYMPI is untested
+
** - MYMPI is untested
**** Need to compile stuff.
+
** Need to compile stuff.
  
 
=== Plume Tracking - Bryan and Brad ===
 
=== Plume Tracking - Bryan and Brad ===
* Setup
+
==== Setup ====
** physical simulator setup approximately 16 inches away and perpendicular to the line of sight of a web enabled camera.
+
* physical simulator setup approximately 16 inches away and perpendicular to the line of sight of a web enabled camera.
** A script was used to capture output of the output of the camera from the server at a rate of one every two seconds. A faster rate may be possible, but the current script did not have time to get the image and rename it within a 1 second interval.
+
* A script was used to capture output of the output of the camera from the server at a rate of one every two seconds. A faster rate may be possible, but the current script did not have time to get the image and rename it within a 1 second interval.
* Procedure
+
==== Procedure ====
** set pump flow rate at maximum and allow water table to equalize
+
* set pump flow rate at maximum and allow water table to equalize
** start image capture script
+
* start image capture script
** inject a full pipette bulb into well number 1
+
* inject a full pipette bulb into well number 1
** remove pipette before allowing bulb to reinflate
+
* remove pipette before allowing bulb to reinflate
** allow simulator to run for approximately 5 minutes or until the majority of the dye in the system has been discharged
+
* allow simulator to run for approximately 5 minutes or until the majority of the dye in the system has been discharged
** stop image capture script
+
* stop image capture script
  
We did three complete runs, each with a different dye colors.  We used blue, purple and green because we thought they would give the most contrast for edge detection.
+
==== Raw Images ====
 +
We did three complete runs, each with a different dye colors.  We used blue, purple and green because we thought they would give the most contrast between the dye and sand.
  
 
* [http://cs.earlham.edu/~carrick/run1.tgz run1.tgz] (blue)
 
* [http://cs.earlham.edu/~carrick/run1.tgz run1.tgz] (blue)
 
* [http://cs.earlham.edu/~carrick/run2.tgz run2.tgz] (purple)
 
* [http://cs.earlham.edu/~carrick/run2.tgz run2.tgz] (purple)
 
* [http://cs.earlham.edu/~carrick/run3.tgz run3.tgz] (green)
 
* [http://cs.earlham.edu/~carrick/run3.tgz run3.tgz] (green)
 +
 +
==== Modified Images ====
 +
Each image was adjusted to fix lens distortion, cropped and had its colors inverted for higher contrast. A ruler was added and the webcam's timestamp was preserved.  The image processing was done via batch jobs in Photoshop.
 +
 +
The batch actions are as follows:
 +
# select time portion of image
 +
# cut
 +
# paste into new layer
 +
# select the original background layer
 +
# Apply a lens correction
 +
## distort amount: +6
 +
## rotate: -0.77 degrees
 +
## vertical perspective: -4
 +
# open image file of ruler
 +
# copy all
 +
# copy
 +
# close file
 +
# paste as new layer
 +
# move the current layer (ruler) to the correct final location under the gws
 +
# crop image
 +
# select the time layer
 +
# move it to the proper location
 +
# select background layer
 +
# invert colors
 +
# auto levels
 +
# auto contrast
 +
# auto colors
 +
# save file
 +
 +
Many of these processes could be replicated through the freely available, command-line driven program ImageMagick.  The inversion, cropping, rotation and merging of photo files are well within IM's scope.  The distorting and cropping are very specific to the position of the camera in relation to the groundwater simulator and would need to be adjusted each time the camera was moved.  If this was to be a regular occurrence, it may be beneficial to have a flat piece of cardboard with a grid on it that could be placed directly in front of the simulator.  A photo of this would give a reference for the distortion caused by the camera's lens in its given location.
 +
 +
* [http://cs.earlham.edu/~carrick/run1_mod.tgz run1_mod.tgz] (inverted blue)
 +
* [http://cs.earlham.edu/~carrick/run2_mod.tgz run2_mod.tgz] (inverted purple)
 +
* [http://cs.earlham.edu/~carrick/run3_mod.tgz run3_mod.tgz] (inverted green)
 +
 +
==== Movies ====
 +
The movies were created with Quicktime Pro's "open image sequence".  QT does not appear to have the capability to have custom framerates outside of their standard choices.  This means that the actual simulation and the movie of the simulation run at different speeds.
 +
* [http://cs.earlham.edu/~carrick/run1.mov run1.mov] (inverted blue)
 +
** Actual run time - 5:18 (1.59 seconds/frame average)
 +
** Movie run time - 6:40 (2 seconds/frame )
 +
** Frames - 200
 +
* [http://cs.earlham.edu/~carrick/run2.mov run1.mov] (inverted purple)
 +
** Actual run time - 4:40 (2.12 seconds/frame average)
 +
** Movie run time - 4:24 (2 seconds/frame)
 +
** Frames - 132
 +
* [http://cs.earlham.edu/~carrick/run3.mov run1.mov] (inverted green)
 +
** Actual run time - 5:24 (2.16 seconds/frame average)
 +
** Movie run time - 5:00 (2 seconds/frame)
 +
** Frames - 150

Latest revision as of 21:01, 11 December 2007

Using the enVision Tabletop Groundwater Simulator

General Instructions

  • Setup
  • Teardown and cleaning
  • Packing and travelling

Instructions for Demonstrations

  • First one
  • Second one
  • etc.

Computational Groundwater Simulations

Fitz, Bryan and Mikio

Confined aquifier simulation
Parabolic contaminant flow model

Experiments

  • Demonstrating porosity
    • model water flow unconfined aquifier
  • Illustrating groundwater flow in a confined aquifer
    • We will use a cellular automata model where at the lowest level, a cell is either fresh water or contaminated. We see this problem split into two concepts - speed and direction.
      • Direction: The illustration to the right demonstrates our assumptions about how the water will move through the material. The simulation will calculate a new direction at each generation based on it's position relative to the known locations of water input and output.
      • Speed: Remains constant throughout generations for a given run. The "speed" value represents a combination of speed of water flow and material porosity, and in terms of the simulation is the possibility that a a neighboring cell in the flow direction becomes contaminated.
  • Describing recharge, transition and discharge areas
    • modeling behavior of water recharge, discharge in wells, lake, etc

Computational Tools

  • C
    • +Very fast
    • +Libraries are available
    • +Good distributed Libraries
    • -Potentially difficult to use
    • -no graphics libraries
  • Netlogo
    • +Fancy Graphics
    • +Fun to use
    • +Available examples/code
    • -Slow
    • -Small problem size
    • -No Distributed processing

Peter and Mikio

Experiment

  • Describing the model
    • Describing the various parts of the Groundwater Simulator by attaching tags: Key words -- wells, artesian wells, lake, underground storage tank, septic tank, springs, vegetative layer, river/ocean, recharge area, discharge area, aquifers, confining layer, clay layers
  • Illustrating and Calculating Porasity of different types of earth materials
  • Determining how it is easy for ground water to move in different earth materials.

Computetional Tool

  • NetLogo for computatinal experiment

Brad and Nate

Our goal is an incremental approach towards illustrating groundwater contamination in a confined aquifer. The confined aquifer, viewed between wells 1 and 8, offers an environment within the groundwater simulator with the fewest variables. The first 4 experiments are an effort to illustrate the behavior and underlying science that must be understood and demonstrated in the final experiment.

Experiments

  • Diffusion
    • Show diffusion without groundwater movement.
  • Flow Rate
    • Show the leading edge of groundwater contamination as a indicator of flow rate (related to section 5 and 13 in manual)
  • Contaminant Plume Length
    • Determine whether contaminant plume length is affected by flow rate for a given amount of dye
  • Soil Density
    • Use displacement method and measurements of aquifer component to determine the density of the soil. We can use this value in silico.
  • Illustrate laminar flow in a confined aquifer (Activity 7-1)
    • Show laminar flow between wells 1 and 8.

Computational Tools

  • NetLogo
    • + Visualization built in
    • + Agent and cell based simulation structure built in
    • - Possible limitation on world size / agent count in RAM
    • - Possible run time slower than groundwater simulator at higher flow rates
    • - Not parallel
  • Python and MYMPI
    • + Parallelizable
    • + Faster than NetLogo in serial code ?
    • + Visualization software exists
  • TKInter - easy to install; seemingly easy to use
    • - Visualization software must be integrated
    • - MYMPI is untested
    • Need to compile stuff.

Plume Tracking - Bryan and Brad

Setup

  • physical simulator setup approximately 16 inches away and perpendicular to the line of sight of a web enabled camera.
  • A script was used to capture output of the output of the camera from the server at a rate of one every two seconds. A faster rate may be possible, but the current script did not have time to get the image and rename it within a 1 second interval.

Procedure

  • set pump flow rate at maximum and allow water table to equalize
  • start image capture script
  • inject a full pipette bulb into well number 1
  • remove pipette before allowing bulb to reinflate
  • allow simulator to run for approximately 5 minutes or until the majority of the dye in the system has been discharged
  • stop image capture script

Raw Images

We did three complete runs, each with a different dye colors. We used blue, purple and green because we thought they would give the most contrast between the dye and sand.

Modified Images

Each image was adjusted to fix lens distortion, cropped and had its colors inverted for higher contrast. A ruler was added and the webcam's timestamp was preserved. The image processing was done via batch jobs in Photoshop.

The batch actions are as follows:

  1. select time portion of image
  2. cut
  3. paste into new layer
  4. select the original background layer
  5. Apply a lens correction
    1. distort amount: +6
    2. rotate: -0.77 degrees
    3. vertical perspective: -4
  6. open image file of ruler
  7. copy all
  8. copy
  9. close file
  10. paste as new layer
  11. move the current layer (ruler) to the correct final location under the gws
  12. crop image
  13. select the time layer
  14. move it to the proper location
  15. select background layer
  16. invert colors
  17. auto levels
  18. auto contrast
  19. auto colors
  20. save file

Many of these processes could be replicated through the freely available, command-line driven program ImageMagick. The inversion, cropping, rotation and merging of photo files are well within IM's scope. The distorting and cropping are very specific to the position of the camera in relation to the groundwater simulator and would need to be adjusted each time the camera was moved. If this was to be a regular occurrence, it may be beneficial to have a flat piece of cardboard with a grid on it that could be placed directly in front of the simulator. A photo of this would give a reference for the distortion caused by the camera's lens in its given location.

Movies

The movies were created with Quicktime Pro's "open image sequence". QT does not appear to have the capability to have custom framerates outside of their standard choices. This means that the actual simulation and the movie of the simulation run at different speeds.

  • run1.mov (inverted blue)
    • Actual run time - 5:18 (1.59 seconds/frame average)
    • Movie run time - 6:40 (2 seconds/frame )
    • Frames - 200
  • run1.mov (inverted purple)
    • Actual run time - 4:40 (2.12 seconds/frame average)
    • Movie run time - 4:24 (2 seconds/frame)
    • Frames - 132
  • run1.mov (inverted green)
    • Actual run time - 5:24 (2.16 seconds/frame average)
    • Movie run time - 5:00 (2 seconds/frame)
    • Frames - 150