![]() |
|
![]() |
||||
|
|
||||||
Be the first pioneers to continue the Astronomy Discussions at our new Astronomy meeting place... The Space and Astronomy Agora |
| Computer Modelling.
Forum List | Follow Ups | Post Message | Back to Thread Topics | In Response To Posted by Alexander on April 6, 2001 21:18:52 UTC |
|
Computer modelling is not difficult, but it requires to know some calculus and some programming language. Most laws of physics become differential equations when it comes to calculate motion or flow of something. Say, if your task is to model comet's orbit around Sun - then you take Newton law of gravity: F=-GMm/r2 and write a differential equation: m(d2t/dt2)=-GMm/r2). But unlike in this simple case, many differential equations do not have analytical solution, thus have to be solved numerically. (Also, if there is one more gravitating body around your asteroid or on its way, then there is no analytical solution for such gravitating system too.) Numerical differentiation (and integration) is to split your process into small time intervals - say, from ti to ti+1 and replace derivatives within this interval by ratios of finite differencies like: dr/dt = (ri+1-ri)/(ti+1-ti) and d^2r/dt^2=(ri+1 -2ri+ri-1)/((ti+1-ti)(ti-ti-1), thus instead of differential equation you have a simple algebraic equation. Computer can solve it in a millisecond or so (using C, C++ or fortran), so you should write a simple code which starts from given initial or boundary values of r and dr/dt and determines ri, and then repeats it step by step using previous values to calculate next ones, thus showing, for example,dot-to-dot orbit of asteroid around mass M or perturbation of its orbit when passing by Earth. As r is a vector in 3-D, it has 3 components, so you actually should have a set of 3 algebraic equations for x,y, and z components of r.
|
|
|
| Additional Information |
|---|
| About Astronomy Net | Advertise on Astronomy Net | Contact & Comments | Privacy Policy |
|
Unless otherwise specified, web site content Copyright 1994-2026 John Huggins All Rights Reserved Forum posts are Copyright their authors as specified in the heading above the post. "dbHTML," "AstroGuide," "ASTRONOMY.NET" & "VA.NET" are trademarks of John Huggins |