From GeoMod
'''Initial:'''
from visual import *
scene.forward = (0, -0.5,1) scene.up = (0,0,1) scene.forward = (0, 1,0) scene.lights = [vector(0, -10, 0)]
xaxis = curve(pos=[ (10, 0, 0), (0, 0,0)], color=(1,0,0)) yaxis = curve(pos=[ (0, 10, 0), (0, 0,0)], color=(0,1,0)) zaxis = curve(pos=[ (0, 0, 10), (0, 0,0)], color=(0,0,1))
ball = sphere(pos=(1,0,10), color=(1,0,0),radius=1, vel=0) ball2 = sphere(pos=(1,5,10), color=(4,1,0),radius=1.5, vel=0)
g = -9.8 dt = 0.01 t = 0.0 s = .5 *(ball.vel * (t*t)) + 10 ct = 0
while 1:
#rate(50)
t = t + dt
ball.vel = ball.vel + g * dt
ball2.vel = ball2.vel + g * dt
ball.pos.z = ball.pos.z + ball.vel * dt
ball2.pos.z = ball2.pos.z + ball2.vel * dt
if ball.pos.z < 1:
ball.pos.z = 1
if ball2.pos.z < 1.5:
ball2.pos.z = 1.5
ball.vel = -ball.vel - 2.25
ball2.vel = -ball2.vel - .5
s = .5 *(ball.vel * (t*t)) + 10
ball.pos.x = ball.pos.x + .09 + ball.pos.z
ball2.pos.x = ball2.pos.x + .0009 + ball2.pos.z
if ball.pos.x >= 9:
break
if ball2.pos.x >= 9:
break
ct = ct + 1
print t,ball.vel,ball.pos.z,s,ball.pos.x,ball2.vel,ball2.pos.z,s,ball2.pos.x,ct
'''The Final:'''
from visual import * from random import uniform, seed
scene.forward = (0, -0.5,1) scene.up = (0,0,1) scene.forward = (0, 1,0) scene.lights = [vector(0, -10, 0)]
xaxis = curve(pos=[ (10, 0, 0), (0, 0,0)], color=(1,0,0)) yaxis = curve(pos=[ (0, 10, 0), (0, 0,0)], color=(0,1,0)) zaxis = curve(pos=[ (0, 0, 10), (0, 0,0)], color=(0,0,1))
seed (2)
balls = [] for i in range(5):
rad = uniform(.5, 2)
vx = uniform(.005, .01)
balls.append (sphere(radius = rad, velx = vx, pos = (1,0,10), velz = 0))
g = -9.8 dt = 0.01 vel = 0.0 t = 0
while 1:
t +=dt
for i in balls:
i.velz = i.velz + g * dt
i.pos.z = i.pos.z + i.velz * dt
if i.pos.z < i.radius:
i.pos.z = i.radius
i.velz = -i.velz * 0.75
i.velx = i.velx * 0.75
i.pos.x = i.pos.x + i.velx
=='''08 February 08'''==
Steady State Flow Equation:
<math>frac{\nabla^2h=0</math>}
where:
frac{\delta^2h/delta_x^2+0}
area= frac{\delta_v/delta_t}=frac{\delta(delta_s/delta_t)/delta_t}
gravity=delta^2"s"/delta"t"^2
Substituting our original formula:
h_i = sq rt(frac\(h_i/-1)^2 + frac\(h_i+1)^2/2)
This is the Final formula:
h<sub>i</sub> = ((delta"h"/delta"x")sub"i"+1/2 - (delta"h"/delta"x")sub"i"-1/2)/delta"x"
[[Image:HPIM1689.JPG|thumb|400px]]
This is the start of the equation
[[Image:HPIM1690.JPG|thumb|400px]]
This is the next step
[[Image:HPIM1691.JPG|thumb|400px]]
This is another step in the substituion of the equation
[[Image:HPIM1686.JPG|thumb|400px]]
{{break}}
=='''15 February 08'''==
The Excel Spreadsheet & Charts of Darcy's Law for the Dam and for Temperature Gradients
[[Image:Darcys_Law_Temp.xls]]
[[Image:Darcys_Law_Water.xls]]
=='''18 February 08'''==
Class discussion included Dr. Cox of the Geology Departement who discussed the formation of the Mississippi Embayment. This was caused by the passing over of the lower crust of the earth over a moving volcanic Plume. This Plume contained enough heat and upwelling force that it caused a hot spot swelling at the embayment. The embayment was susceptable to this typw of heat trasnfer because of the makeup of the existing material. A makeup of of lesser material set between the two mountain ranges located in and around Denver, Colorado and the Appalachian. As this swell formed and the outcome was a formation of hills, mountains, etc., as quickly as these landforms emerged, the inverse occurred because of the makeup within these features. The features include not only the weakened or lesser material, but deposits of the material from the Plume. As this passed over the Mississippi area, what was eventually left was a depression of an elipsiodal shape. As the depression took on a shape, other depositions gathered in the new formed area. This also became a colection place for water.
The class discussed how this type of movement could be characterized and studied in a model in the Python program. The parameters of the model will include a grid of particular size (2Km X 2Km) that will include the embayment (size 300m X 400m) with each cell within the grid at 50 Km. The 2D model will eventually include the erosion, heat transfer calculation, the expansion and contraction of the depression, the density, rate of expansion and the effects at the crust. There will have to be an agreement of the the different layer sizes and compostion
Once the grid is constructed, another model will include the new grid and show the depression of the embayment. After this has been formualted, the class will begin to incorporate the other aspects of this project.
=='''25 February 08'''==
The class discussed the simplification of the formulas that will be used in the next modeling exercise. There was a homework assignment regarding simplifying the solutions for the next portion of the input formula; the explicit and implicit calculations.
=='''27 February 08'''==
The class was given the correct set up for the formulas and there was a demonstration in excel to determine the limits with which the formula would continue to compute before becoming erratic. The final solution was to use the more simpler formula; the explicit because of the steps needed to perform calculations and the limits were larger for calculations.
=='''29 February 08'''==
There were two presentations given to explain the formation of the which yielded the results for the embayment. The presentations were made to explain how the formulas were applied in order to determine how the crust was effected and what conclusions can be made from the information. The results will be used for the second part of the model which will show the formation of the depression as the crust passes over the volcanic plume.
==''' 05 March 08'''==
The parameters of the model were discussed again in terms of looking at the processes in three deminsion. However, after a peliminary model was reviewed it was determined that the limits have to be revised. In addition, all supplemental information associated with finding the formulas and processes to be used in the model will be posted as links for future studies and reviews.
==''' 07 March 08'''==
The thermal conductivity formula as well as the previous steady state flow formula will be used in determining the amount of uplift at each cell absorbed by the plume below it with a particular temperature. Once the heat energy is absorbed by the cell, the conductivity of the cell will move the heat energy through the cell and out through the top of the cell wiht some loss to the pyhsical makeup of the cell and the amount of heat energy retained by the cell. Then the cell must dissipate heat energy laterally into other adjoining cells where those affected cells will then repeat the same initial process. However, the previous cell will be subjected to some additional absorption given by the new affected cells and the process is again repeated until a certain amount of cells have been affected by this plume and all heat intially given into the cells have been absorbed and dissipated.
