Thursday, 11 October 2012

Simplified Motion Part Deux

I have factored in a drag force of the form:  F_D\, =\, \tfrac12\, \rho\, v^2\, C_d\, A,

However, according to my current code, this generates a drag force which exceeds the net acceleration of the bubble as it floats upwards. This doesn't make sense, so I'll have to work on that.

There was a small victory on the plot front. I am now capable of saving values from a loop into an array. Huzzah!

Here is an example of the newly refined (yet still simple) "uvast.py":



In [146]: % run -i uvast.py
Mass (Mkg): 50
Radial Fraction: 0.7
Temperature (MK): 2.5
That took 1672.400 seconds


Coming soon: drag force which obeys laws of physics.




No comments:

Post a Comment