|
|
| Line 1: |
Line 1: |
| '''Orbit modeling''' is the process of creating mathematical models to simulate motion of a massive body as it moves in orbit around another massive body due to [[gravity]]. Other forces such as gravitational attraction from tertiary bodies, [[Drag (physics)|air resistance]], [[Solar wind|solar pressure]], or thrust from a [[propulsion]] system are typically modeled as secondary effects. Directly modeling an orbit can push the limits of [[machine precision]] due to the need to model small perturbations to very large orbits. Because of this, [[Perturbation (astronomy)|perturbation]] methods are often used to model the orbit in order to achieve better accuracy.
| | In decreasing term policy, the level of [https://Www.Flickr.com/search/?q=coverage coverage] will reduce by each year. It is vital that you disclose as much information about yourself as possible or else your policy may be invalidated. It is a bit kooky feeling but that is how term life insurance works. Because of this, [http://www.google.Co.uk/search?hl=en&gl=us&tbm=nws&q=choosing&gs_l=news choosing] your life insurance terms is one of the most important decisions you can make for your family. You will feel more relaxed if you own a 15 year term policy which would provide sufficient income to take care of the family until the last child graduates. " Andy was out of earshot reading a Bugs Bunny book up in his room when this particularly callous remark was uttered. This policy is for one year and has an option to renew it for a specific number of years. <br> <br> It makes it a lot easier for you to decide which policy you want when you can get a cheap online term life insurance quote. Figuring out how much life insurance to buy can be tricky. This category is ideal suited to late-middle-aged and above. Account Value: This is the accumulated gross value of all the investments contributed to the policy which include the income after deducting all the current monthly expenses. However, one should have a clear understanding of payroll deductions in order to distinguish the taxable from the non-taxable. We provide the following types of life insurance policies:. In this type of policy, you are led to believe that you are paying for both full coverage on your life and you are also investing your money into a savings type of fund. <br> <br> Without LTCI you also cannot force a nursing home to take you in and reiterate that apart from your functional limitations, your hypertension has gotten worse and your doctor's recommendation is that you stay in a nursing home. Most people with life insurance have term life insurance. You can perform all your proceedings online and sign up your registration form through the RBC insurance company's highly protected electronic-signature technology, which offers an electronic method for you to sign your forms and testify to the validness of the information you have given. It is important to note that rates do slowly rise as you age. The reason being that APRA has the responsibility of administering: registration of life insurance companies and authorized non-operating holding companies, life insurance company statutory funds, financial management in life companies, judicial management, winding up, among other issues. Once you annuitize, your gains are taxed at your normal tax rate. Everybody understands the solution: comparison-shopping. <br> <br> Let's first look at the reasons your life insurance will not pay a claim. Insurers would provide bowel life insurance for those suffering from IBS in most cases. When you understand these 2 basic types, then you will have an easier time understanding any of the variations. You can surrender an annuity contract which is still in the accumulation phase at any time. We will look at how taxes may affect your insurance pay outs, and how you can go about making sure your loved ones receive all the benefits they are entitled to after your death. These term life insurance companies have spent a lot of money advertising their message during prime time. Life insurance should be kept for the duration of your existence. <br> <br> You could convert your term policy to a permanent policy at that time to guarantee that you always have coverage. Many different insurance companies have different critical illness insurance policies and so it is imperative that you know very well what is going to be covered in your policy and what won't be covered. They will most likely not be able to afford the premiums, and would most likely quality for Medicaid, which does cover long term care. You see if you buy a 25 or 30 year term policy from the outset then all is taken care of. Even if your term has finished, you can renew it every year according to your own wish, so that they can be able to get the coverage later on. So, if you are considering a senior life policy, try to qualify for simplified issue instead of guaranteed issue if you can. The Breatharian Diet' All I Need is the Air That I Breath. <br> <br> Be sure to pick a term length that will last until your dependents are no longer in need of your income. He looks into the variety of [http://www.termlifepolicy.com/insurance-companies/reviews/aig-sunamerica/ different Life Insurance policies] available on the market. Those people who have a personal business can add life insurance to their normal insurance coverage plans. If you are shopping for life insurance products for some time now, you might have already been lured by the appealing and somewhat clichd verbiages that insurance companies use more often than not. To save some money consider switching to annual life insurance premiums instead of monthly life insurance premiums. Before you ask for an extended term insurance offer or perhaps compare life insurance companies, determine that you are doing require one simply because a real insurance policy is just not for everybody since most of them can be instead highly-priced. In return for the term life insurance, the insurer pays the insurance company a premium. <br> <br> In due time, we all have to face the very end leaving all the things we toiled for so many years to our beloved families. It used to be when you were interested in purchasing term life insurance you would get an agent to come and see you, tell him or her what you have in mind and get a few illustrations showing the death benefits and the premium costs. Use an online service such as to do comparison-shopping. Despite simplicity and limitations, term life insurance coverage is still sensible among many customers. That bailout of AIG also violated the Social Security Act and the FEMA Act, disproving the United States Senate understands prohibitions against organized crime. A family with young children may need a 20 year or a 25 year term policy to protect the family through to the end of the college years but because of cost they use the 10 year term policy until they can afford the more costly policies. However, you can find plenty of term life insurance policies and providers who are willing to work with you, despite your age. |
| | |
| == Background == | |
| The study of orbital motion and mathematical modeling of orbits began with the first attempts to predict planetary motions in the sky, although in ancient times the causes remained a mystery. [[Isaac Newton|Newton]], at the time he formulated his laws of [[Newton's laws of motion|motion]] and of [[Newton's law of universal gravitation|gravitation]], applied them to the first analysis of perturbations,<ref name="moulton"/> recognizing the complex difficulties of their calculation.<ref name=moulton>
| |
| {{cite book
| |
| |url=http://books.google.com/books?id=jqM5AAAAMAAJ&printsec=frontcover&source=gbs_ge_summary_r&cad=0#v=onepage&q&f=false
| |
| |title=An Introduction to Celestial Mechanics, Second Revised Edition
| |
| |last=Moulton
| |
| |first=Forest Ray
| |
| |year=1914}} chapter IX. (at [http://books.google.com/books Google books])</ref>
| |
| Many of the great mathematicians since then have given attention to the various problems involved; throughout the 18th and 19th centuries there was demand for accurate tables of the position of the Moon and planets for purposes of navigation at sea.
| |
| | |
| The complex motions of orbits can be broken down. The hypothetical motion that the body follows under the gravitational effect of one other body only is typically a [[conic section]], and can be readily modeled with the methods of [[geometry]]. This is called a [[two-body problem]], or an unperturbed [[Kepler orbit|Keplerian orbit]]. The differences between the Keplerian orbit and the actual motion of the body are caused by [[Perturbation (astronomy)|perturbations]]. These perturbations are caused by forces other than the gravitational effect between the primary and secondary body and must be modeled to create an accurate orbit simulation. Most orbit modeling approaches model the two-body problem and then add models of these perturbing forces and simulate these models over time. Perturbing forces may include gravitational attraction from other bodies besides the primary, solar wind, drag, magnetic fields, and propulsive forces.
| |
| | |
| Analytical solutions (mathematical expressions to predict the positions and motions at any future time) for simple two-body and [[three-body problem]]s exist; none have been found for the [[n-body problem|''n''-body problem]] except for certain special cases. Even the two-body problem becomes insoluble if one of the bodies is irregular in shape.<ref name="roy">{{cite book |last = Roy |first = A.E. | title = Orbital Motion |publisher = Institute of Physics Publishing |edition = third |isbn = 0-85274-229-0 |year=1988}}, chapters 6 and 7.</ref>
| |
| | |
| Due to the difficulty in finding analytic solutions to most problems of interest, computer [[modeling and simulation]] is typically used to analyze orbital motion. Commercial software applications such as [[Satellite Tool Kit]] have been created for the specific purpose of simulating orbits and trajectories of spacecraft.
| |
| | |
| == Keplerian orbit model ==
| |
| {{main|Kepler orbit}}
| |
| In its simplest form, an orbit model can be created by assuming that only two bodies are involved, both behave as spherical point-masses, and that no other forces act on the bodies. For this case the model is simplified to a [[Kepler orbit]].
| |
| | |
| Keplerian orbits follow [[conic sections]]. The mathematical model of the orbit which gives the distance between a central body and an orbiting body can be expressed as:
| |
| | |
| :<math> r(\nu) = \frac{a(1-e^2)}{1+e\cos(\nu)} </math>
| |
| | |
| Where:
| |
| :<math>r</math> is the distance
| |
| :<math>a</math> is the [[semi-major axis]], which defines the size of the orbit
| |
| :<math>e</math> is the [[orbital eccentricity|eccentricity]], which defines the shape of the orbit
| |
| :<math>\nu</math> is the [[true anomaly]], which is the angle between the current position of the orbiting object and the location in the orbit at it is closest to the central body (called the [[periapsis]])
| |
| Alternately, the equation can be expressed as:
| |
| | |
| :<math> r(\nu) = \frac{p}{1+e\cos(\nu)} </math>
| |
| | |
| Where <math>p</math> is called the [[conic section#Features|semi-latus rectum]] of the curve. This form of the equation is particularly useful when dealing with parabolic trajectories, for which the semi-major axis is infinite.
| |
| | |
| An alternate approach uses [[Isaac Newton]]'s [[Newton's law of universal gravitation|law of universal gravitation]] as defined below:
| |
| | |
| :<math>F = G \frac{m_1 m_2}{r^2}</math>
| |
| | |
| where:
| |
| | |
| :<math>F</math> is the magnitude of the gravitational force between the two point masses
| |
| :<math>G</math> is the [[gravitational constant]]
| |
| :<math>m_1</math> is the mass of the first point mass
| |
| :<math>m_2</math> is the mass of the second point mass
| |
| :<math>r</math> is the distance between the two point masses
| |
| | |
| Making an additional assumption that the mass of the primary body is much greater than the mass of the secondary body and substituting in Newton's [[Newton's laws of motion|second law of motion]], results in the following differential equation
| |
| :<math> \ddot{\mathbf{r}} = \frac{G m_1}{r^2} \mathbf{\hat{r}}</math>
| |
| | |
| Solving this differential equation results in Keplerian motion for an orbit.
| |
| In practice, Keplerian orbits are typically only useful for first-order approximations, special cases, or as the base model for a perturbed orbit.
| |
| | |
| == Orbit simulation methods ==
| |
| {{main|Perturbation (astronomy)}}
| |
| Orbit models are typically propagated in time and space using special [[Perturbation (astronomy)|perturbation]] methods. This is performed by first modeling the orbit as a Keplerian orbit. Then perturbations are added to the model to account for the various perturbations that affect the orbit.<ref name="moulton"/>
| |
| Special perturbations can be applied to any problem in [[celestial mechanics]], as it is not limited to cases where the perturbing forces are small.<ref name="roy"/> Special perturbation methods are the basis of the most accurate machine-generated [[Fundamental ephemeris|planetary ephemerides]].<ref name="moulton">
| |
| see, for instance, [[Jet Propulsion Laboratory Development Ephemeris]]</ref>
| |
| | |
| === Cowell's method ===
| |
| [[File:Cowells method.png|thumb|Cowell's method. Forces from all perturbing bodies (black and gray) are summed to form the total force on body ''i'' (red), and this is numerically integrated starting from the initial position (the ''epoch of osculation'').]]
| |
| | |
| Cowell's method is perhaps the simplest of the special perturbation methods;<ref>
| |
| So named for [[Philip Herbert Cowell|Philip H. Cowell]], who, with A.C.D. Cromellin, used a similar method to predict the return of Halley's comet.
| |
| {{cite book
| |
| |last1 = Brouwer
| |
| |first1 = Dirk
| |
| |last2 = Clemence
| |
| |first2 = Gerald M.
| |
| | title = Methods of Celestial Mechanics
| |
| |publisher = Academic Press, New York and London
| |
| |year=1961}}, p. 186.
| |
| </ref>
| |
| mathematically, for <math>n</math> mutually interacting bodies, [[Newton's law of universal gravitation|Newtonian]] forces on body <math>i</math> from the other bodies <math>j</math> are simply summed thus,
| |
| | |
| : <math>\mathbf{\ddot{r}}_i = \sum_{\underset{j \ne i}{j=1}}^n {Gm_j (\mathbf{r}_j-\mathbf{r}_i) \over r_{ij}^3}</math>
| |
| | |
| where
| |
| :<math>\mathbf{\ddot{r}}_i</math> is the [[acceleration]] vector of body <math>i</math>
| |
| :<math>G</math> is the [[gravitational constant]]
| |
| :<math>m_j</math> is the [[mass]] of body <math>j</math>
| |
| :<math>\mathbf{r}_i</math> and <math>\mathbf{r}_j</math> are the [[position vector]]s of objects <math>i</math> and <math>j</math>
| |
| :<math>r_{ij}</math> is the distance from object <math>i</math> to object <math>j</math>
| |
| with all [[Euclidean_vector#Physics|vectors]] being referred to the [[Center_of_mass#Astronomy|barycenter]] of the system. This equation is resolved into components in <math>x</math>, <math>y</math>, <math>z</math> and these are integrated numerically to form the new velocity and position vectors as the simulation moves forward in time. The advantage of Cowell's method is ease of application and programming. A disadvantage is that when perturbations become large in magnitude (as when an object makes a close approach to another) the errors of the method also become large.<ref name="danby">
| |
| {{cite book
| |
| |last = Danby
| |
| |first = J.M.A.
| |
| | title = Fundamentals of Celestial Mechanics
| |
| |publisher = Willmann-Bell, Inc.
| |
| |edition = second |isbn = 0-943396-20-4
| |
| |year=1988}}, chapter 11.</ref>
| |
| Another disadvantage is that in systems with a dominant central body, such as the [[Sun]], it is necessary to carry many [[Significant figures|significant digits]] in the [[arithmetic]] because of the large difference in the forces of the central body and the perturbing bodies.<ref>
| |
| {{cite book
| |
| |last = Herget
| |
| |first = Paul
| |
| | title = The Computation of Orbits
| |
| |publisher = privately published by the author
| |
| |year=1948}}, p. 91 ff.</ref>
| |
| | |
| === Encke's method ===
| |
| | |
| [[File:Enckes method.PNG|thumb|Encke's method. Greatly exaggerated here, the small difference δ'''r''' (blue) between the osculating, unperturbed orbit (black) and the perturbed orbit (red), is numerically integrated starting from the initial position (the ''epoch of osculation'').]]
| |
| | |
| Encke's method begins with the [[osculating orbit]] as a reference and integrates numerically to solve for the variation from the reference as a function of time.<ref>
| |
| So named for [[Johann Franz Encke]];
| |
| {{cite book
| |
| |last = Battin
| |
| |first = Richard H.
| |
| | title = An Introduction to the Mathematics and Methods of Astrodynamics, Revised Edition
| |
| |publisher = American Institute of Aeronautics and Astronautics, Inc.
| |
| |isbn = 1-56347-342-9
| |
| |year=1999}}, p. 448</ref>
| |
| Its advantages are that perturbations are generally small in magnitude, so the integration can proceed in larger steps (with resulting lesser errors), and the method is much less affected by extreme perturbations than Cowell's method. Its disadvantage is complexity; it cannot be used indefinitely without occasionally updating the osculating orbit and continuing from there, a process known as ''rectification''.<ref name="danby"/>
| |
| <ref>Battin (1999), sec. 10.2.</ref>
| |
| | |
| Letting <math>\boldsymbol{\rho}</math> be the [[Position vector|radius vector]] of the [[osculating orbit]], <math>\mathbf{r}</math> the radius vector of the perturbed orbit, and <math>\delta \mathbf{r}</math> the variation from the osculating orbit,
| |
| | |
| {{NumBlk|:|<math>\delta \mathbf{r} = \mathbf{r} - \boldsymbol{\rho},</math> and the [[Equations of motion|equation of motion]] of <math>\delta \mathbf{r}</math> is simply|{{EquationRef|1}}}}
| |
| | |
| {{NumBlk|:|<math>\ddot{\delta \mathbf{r}} = \mathbf{\ddot{r}} - \boldsymbol{\ddot{\rho}}.</math>|{{EquationRef|2}}}}
| |
| | |
| <math>\mathbf{\ddot{r}}</math> and <math>\boldsymbol{\ddot{\rho}}</math> are just the equations of motion of <math>\mathbf{r}</math> and <math>\boldsymbol{\rho}</math>,
| |
| | |
| {{NumBlk|:|<math>\mathbf{\ddot{r}} = \mathbf{a}_{\text{per}} - {\mu \over r^3} \mathbf{r}</math> for the perturbed orbit and |{{EquationRef|3}}}}
| |
| | |
| {{NumBlk|:|<math>\boldsymbol{\ddot{\rho}} = - {\mu \over \rho^3} \boldsymbol{\rho}</math> for the unperturbed orbit,|{{EquationRef|4}}}}
| |
| | |
| where <math>\mu = G(M+m)</math> is the [[Standard gravitational parameter|gravitational parameter]] with <math>M</math> and <math>m</math> the [[mass]]es of the central body and the perturbed body, <math>\mathbf{a}_{\text{per}}</math> is the perturbing [[acceleration]], and <math>r</math> and <math>\rho</math> are the magnitudes of <math>\mathbf{r}</math> and <math>\boldsymbol{\rho}</math>.
| |
| | |
| Substituting from equations ({{EquationNote|3}}) and ({{EquationNote|4}}) into equation ({{EquationNote|2}}),
| |
| | |
| {{NumBlk|:|<math>\ddot{\delta \mathbf{r}} = \mathbf{a}_{\text{per}} + \mu \left( {\boldsymbol{\rho} \over \rho^3} - {\mathbf{r} \over r^3} \right),</math> |{{EquationRef|5}}}}
| |
| | |
| which, in theory, could be integrated twice to find <math>\delta \mathbf{r}</math>. Since the osculating orbit is easily calculated by two-body methods, <math>\boldsymbol{\rho}</math> and <math>\delta \mathbf{r}</math> are accounted for and <math>\mathbf{r}</math> can be solved. In practice, the quantity in the brackets, <math> {\boldsymbol{\rho} \over \rho^3} - {\mathbf{r} \over r^3} </math>, is the difference of two nearly equal vectors, and further manipulation is necessary to avoid the need for extra [[Significant figures|significant digits]].<ref>
| |
| Bate, Mueller, White (1971), sec. 9.3.</ref><ref>
| |
| Roy (1988), sec. 7.4.</ref>
| |
| | |
| === Sperling–Burdet method ===
| |
| In 1991 Victor R. Bond and Michael F. Fraietta created an efficient and highly accurate method for solving the two-body perturbed problem.<ref name=Peláez>{{cite journal|last=Peláez|first=Jesús|coauthors=José Manuel Hedo, Pedro Rodríguez de Andrés|title=A special perturbation method in orbital dynamics|journal=Celestial Mech Dyn Astr|date=13|year=2006|month=October|volume=|DOI =10.1007/s10569-006-9056-3}}</ref> This method uses the linearized and regularized differential equations of motion derived by Hans Sperling and a perturbation theory based on these equations developed by C.A. Burdet. In 1973, Bond and Hanssen improved Burdet's set of differential equations by using the total energy of the perturbed system as a parameter instead of the two-body energy and by reducing the number of elements to 13. In 1989 Bond and Gottlieb embedded the Jacobian integral, which is a constant when the potential function is explicitly dependent upon time as well as position in the Newtonian equations. The Jacobian constant was used as an element to replace the total energy in a reformulation of the differential equations of motion. In this process, another element which is proportional to a component of the angular momentum is introduced. This brought the total number of elements back to 14. In 1991, Bond and Fraietta made further revisions by replacing the Laplace vector with another vector integral as well as another scalar integral which removed small secular terms which appeared in the differential equations for some of the elements.<ref name=Bond>{{cite journal|last=Bond|first=Victor|coauthors=Michael F. Fraietta|title=Elimination Of Secular Terms From The Differential Equations For The Elements of Perturbed Two-Body Motion|journal=Flight Mechanics and Estimation Theory Symposium|year=1991}}</ref>
| |
| | |
| The Sperling–Burdet method is executed in a 5 step process as follows:<ref name="Bond"/>
| |
| :'''Step 1: Initialization'''
| |
| ::Given an initial position, <math>\mathbf{r}_0</math>, an initial velocity, <math>\mathbf{v}_0</math>, and an initial time, <math>t_0</math>, the following variables are initialized:
| |
| ::<math>s=0</math>
| |
| ::<math>r_0=(\mathbf{r}_0\cdot\mathbf{r}_0)^{1/2}</math>
| |
| ::<math>a=r_0</math>
| |
| ::<math>b=\mathbf{r}_0\cdot\mathbf{v}_0</math>
| |
| ::<math>\tau=t_0</math>
| |
| ::<math>\boldsymbol{\alpha}=\mathbf{r}_0</math>
| |
| ::<math>\boldsymbol{\beta}=a\mathbf{v}_0</math>
| |
| ::Perturbations due to perturbing masses, defined as <math>V_0</math> and <math>\Bigg[{\partial{V}\over{\partial{\mathbf{r}}}}\Bigg]_0</math>, are evaluated
| |
| ::Perturbations due to other accelerations, defined as <math>\mathbf{P}_0</math>, are evaluated
| |
| ::<math>\alpha_J=\frac{2\mu}{r_0}-\mathbf{v}_0\cdot\mathbf{v}_0-2V_0</math>
| |
| ::<math>\gamma=\mu-\alpha_J\alpha</math>
| |
| ::<math>\boldsymbol{\delta}=-(\mathbf{v}_0\cdot\mathbf{v}_0)\mathbf{r}_0+(\mathbf{r}_0\cdot\mathbf{v}_0)\mathbf{v}_0+\frac{\mu}{r_0}\mathbf{r}_0-\alpha_J\mathbf{r}_0</math>
| |
| ::<math>\sigma=0</math>
| |
| :'''Step 2: Transform Elements to Coordinates'''
| |
| ::<math>\mathbf{r}=\boldsymbol{\alpha}+\boldsymbol{\beta}sc_1+\boldsymbol{\delta}s^2c_2</math>
| |
| ::<math>\mathbf{r'}=\boldsymbol{\beta}c_0+\boldsymbol{\delta}sc_1</math>
| |
| ::<math>\mathbf{x}_3=\alpha_j(\boldsymbol{\alpha}-\mathbf{r})+\boldsymbol{\delta}</math>
| |
| ::<math>\gamma=\mu-\alpha_ja</math>
| |
| ::<math>r=a+bsc_1+\gamma s^2c_2</math>
| |
| ::<math>\mathbf{v}=\mathbf{r'}/r</math>
| |
| ::<math>r'=bc_0+\gamma sc_1</math>
| |
| ::<math>t=\tau+as+bs^2c_2+\gamma s^3c_3</math>
| |
| ::where <math>c_0, c_1, c_2, c_3</math> are [[Stumpff function]]s
| |
| :'''Step 3: Evaluate Differential Equations for the Elements'''
| |
| ::<math>\mathbf{F}=\mathbf{P}-{\partial{V}\over\partial{\mathbf{r}}}</math>
| |
| ::<math>\mathbf{Q}=r^2\mathbf{F}+2\mathbf{r}(-V+\sigma)</math>
| |
| ::<math>\alpha'_j=2(-\mathbf{r'}+r\boldsymbol{\omega}\times\mathbf{r})\cdot\mathbf{P}</math>
| |
| ::<math>\mu\boldsymbol{\epsilon}'=2(\mathbf{r'}\cdot\mathbf{F})\mathbf{r}-(\mathbf{r}\cdot\mathbf{F})\mathbf{r'}-(\mathbf{r}\cdot\mathbf{r'})\mathbf{F}</math>
| |
| ::<math>\boldsymbol{\alpha}'=-\mathbf{Q}sc_1-\mu\boldsymbol{\epsilon}'s^2c_2-\alpha'_j\big[\boldsymbol{\alpha}s^2c_2+2\boldsymbol{\beta}s^3\bar{c}_3+\frac{1}{2}\boldsymbol{\delta}s^4c^2_2\big]</math>
| |
| ::<math>\boldsymbol{\beta}'=\mathbf{Q}c_0+\mu\boldsymbol{\epsilon}'sc_1+\alpha'_j\big[\boldsymbol{\alpha}sc_1+\boldsymbol{\beta}s^2\bar{c}_2-\boldsymbol{\delta}s^3(2\bar{c}_3-c_1c_2)\big]</math>
| |
| ::<math>\boldsymbol{\delta}'=\mathbf{Q}\alpha_jsc_1-\mu\boldsymbol{\epsilon}'c_0+\alpha'_j\big[-\boldsymbol{\alpha}c_0+2\alpha_j\boldsymbol{\beta}s^3\bar{c}_3+\frac{1}{2}\boldsymbol{\delta}\alpha_js^4c^2_2\big]</math>
| |
| ::<math>\sigma'=r\boldsymbol{\omega}\cdot\mathbf{r}\times\mathbf{F}</math>
| |
| ::<math>a'=-\frac{1}{r}\mathbf{r}\cdot\mathbf{Q}sc_1-\alpha_j'\big[as^2c_2+2bs^3\bar{c}_3+\frac{1}{2}\gamma s^4c^2_2\big]</math>
| |
| ::<math>b'=\frac{1}{r}\mathbf{r}\cdot\mathbf{Q}c_0+\alpha_j'\big[asc_1+bs^2\bar{c}_2-\gamma s^3(2\bar{c}_3-c_1c_2)\big]</math>
| |
| ::<math>\gamma'=-\frac{1}{r}\mathbf{r}\cdot\mathbf{Q}\alpha_jsc_1+\alpha_j'\big[-ac_0+2b\alpha_js^3\bar{c}_3+\frac{1}{2}\gamma \alpha_j s^4c^2_2\big]</math>
| |
| ::<math>\tau'=\frac{1}{r}\mathbf{r}\cdot\mathbf{Q}s^2c_2+\alpha_j'\big[as^3c_3+\frac{1}{2}bs^4c^2_2-2\gamma s^5(c_5-4\bar{c}_5)\big]</math>
| |
| :'''Step 4: Integration'''
| |
| ::Here the differential equations are integrated over a period <math>\Delta s</math> to obtain the element value at <math>s+\Delta s</math>
| |
| :'''Step 5: Advance'''
| |
| ::Set <math>s=s+\Delta s</math> and return to step 2 until simulation stopping conditions are met.
| |
| | |
| == Models of perturbing forces ==
| |
| Perturbing forces cause orbits to become perturbed from a perfect Keplerian orbit. Models for each of these forces are created and executed during the orbit simulation so their effects on the orbit can be determined.
| |
| | |
| === Non-spherical gravity ===
| |
| The Earth is not a perfect sphere nor is mass evenly distributed within the Earth. This results in the point-mass gravity model being inaccurate for orbits around the Earth, particularly [[Low Earth Orbit|Low Earth orbits]]. To account for variations in gravitational potential around the surface of the Earth, the gravitational field of the Earth is modeled with spherical harmonics<ref name=Roithmayr>{{cite journal|last=Roithmayr|first=Carlos|title=Contributions of Spherical Harmonics to Magnetic and Gravitational Fields|journal=NASA/TM–2004–213007|year=2004|month=March}}</ref> which are expressed through the equation:
| |
| :<math>{\mathbf{f}}=-\frac{\mu}{R^2}\mathbf{\hat{r}}+\sum_{n=2}^\infty \sum_{m=0}^n {\mathbf{f}}_{n,m}</math>
| |
| where
| |
| :<math>{\mu}</math> is the gravitational parameter defined as the product of G, the [[universal gravitational constant]], and the mass of the primary body.
| |
| :<math>\mathbf{\hat{r}}</math> is the unit vector defining the distance between the primary and secondary bodies, with <math>{R}</math> being the magnitude of the distance.
| |
| :<math>{\mathbf{f}}_{n,m}</math> represents the contribution to <math>{\mathbf{f}}</math> of the spherical harmonic of degree ''n'' and order ''m'', which is defined as:<ref name="Roithmayr"/>
| |
| | |
| : <math>
| |
| \begin{align}
| |
| \mathbf{f}_{n,m} & = \frac{\mu R_O^2}{R^{n+m+1}} \left(\frac{C_{n,m}\mathcal{C}_m+S_{n,m}\mathcal{S}_m}{R}(A_{n,m+1}\mathbf{\hat{e}}_3 - \left(s_{\lambda} A_{n,m+1}+(n+m+1)A_{n,m}\right)\mathbf{\hat{r}}\right) \\[10pt]
| |
| & {}\quad {}+ mA_{n,m} ((C_{n,m}\mathcal{C}_{m-1} + S_{n,m}\mathcal{S}_{m-1})\mathbf{\hat{e}}_1+(S_{n,m}\mathcal{C}_{m-1}-C_{n,m}\mathcal{S}_{m-1})\mathbf{\hat{e}}_2))
| |
| \end{align}
| |
| </math>
| |
| | |
| where:
| |
| :<math>R_O</math> is the mean equatorial radius of the primary body.
| |
| :<math>R</math> is the magnitude of the position vector from the center of the primary body to the center of the secondary body.
| |
| :<math> C_{n,m}</math> and <math>S_{n,m}</math> are gravitational coefficients of degree ''n'' and order ''m''. These are typically found through [[gravimetry]] measurements.
| |
| :The unit vectors <math>\mathbf{\hat{e}}_1, \mathbf{\hat{e}}_2, \mathbf{\hat{e}}_3</math> define a coordinate system fixed on the primary body. For the Earth, <math>\mathbf{\hat{e}}_1</math> lies in the equatorial plane parallel to a line intersecting Earth's geometric center and the [[Greenwich meridian]],<math>\mathbf{\hat{e}}_3</math> points in the direction of the North polar axis, and <math>\mathbf{\hat{e}}_2=\mathbf{\hat{e}}_3\times\mathbf{\hat{e}}_1</math>
| |
| :<math>A_{n,m}</math> is referred to as a derived [[Legendre polynomial]] of degree ''n'' and order ''m''. They are solved through the recurrence relation: <math>A_{n,m}(u)=\frac{1}{n-m}((2n-1)uA_{n-1,m}(u)-(n+m-1)A_{n-2,m}(u))</math>
| |
| :<math>s_\lambda</math> is sine of the geographic latitude of the secondary body, which is <math>\mathbf{\hat{r}}\cdot\mathbf{\hat{e}}_3</math>.
| |
| :<math>\mathcal{C}_m, \mathcal{S}_{m}</math> are defined with the following recurrence relation and initial conditions:<math>\mathcal{C}_m=\mathcal{C}_1\mathcal{C}_{m-1}-\mathcal{S}_1\mathcal{S}_{m-1}, \mathcal{S}_m=\mathcal{S}_1\mathcal{C}_{m-1}+\mathcal{C}_1\mathcal{S}_{m-1}, \mathcal{S}_0=0, \mathcal{S}_1=\mathbf{R}\cdot\mathbf{\hat{e}}_2, \mathcal{C}_{0}=1, \mathcal{C}_1=\mathbf{R}\cdot\mathbf{\hat{e}}_1</math>
| |
| When modeling perturbations of an orbit around a primary body only the sum of the <math>{\mathbf{f}}_{n,m}</math> terms need to be included in the perturbation since the point-mass gravity model is accounted for in the <math>-\frac{\mu}{R^2}\mathbf{\hat{r}}</math> term
| |
| | |
| === Third-body perturbations ===
| |
| Gravitational forces from third bodies can cause perturbations to an orbit. For example, the [[Sun]] and [[Moon]] cause perturbations to Orbits around the Earth.<ref name=SMAD>{{cite book|last=Larson|first=Wiley|title=Space Mission Analysis and Design|year=1999|publisher=Microcosm Press|location=California|isbn=1-881883-10-8}}</ref> These forces are modeled in the same way that gravity is modeled for the primary body by means of [[N-body simulation|Direct gravitational N-body simulations]]. Typically, only a spherical point-mass gravity model is used for modeling effects from these third bodies.<ref name=Delgado>{{cite web|last=Delgado|first=Manuel|title=Third Body Perturbation Modeling the Space Environment|url=http://ocw.upm.es/ingenieria-aeroespacial/modeling-the-space-environment/contenidos/material-de-clase/mse05_3rdbody.pdf|work=European Masters in Aeronautics and Space|publisher=Universidad Polit ´ecnica de Madrid|accessdate=27 November 2012}}</ref>
| |
| Some special cases of third-body perturbations have approximate analytic solutions. For example, perturbations for the right ascension of the ascending node and argument of perigee for a circular Earth orbit are:<ref name="SMAD"/>
| |
| :<math>\dot{\Omega}_\mathrm{MOON}=-0.00338(\cos(i))/n</math>
| |
| :<math>\dot{\omega}_\mathrm{MOON}=-0.00169(4-5\sin^2(i))/n</math>
| |
| :where:
| |
| :<math>\dot{\Omega}</math> is the change to the right ascension of the ascending node in degrees per day.
| |
| :<math>\dot{\omega}</math> is the change to the argument of perigee in degrees per day.
| |
| :<math>i</math> is the orbital inclination.
| |
| :<math>n</math> is the number of orbital revolutions per day.
| |
| | |
| === Solar radiation ===
| |
| {{main|Radiation pressure}}
| |
| Solar radiation pressure causes perturbations to orbits. The magnitude of acceleration it imparts to a spacecraft in Earth orbit is modeled using the equation below:<ref name="SMAD"/>
| |
| :<math>a_R\approx -4.5 \times 10^{-6}(1+r)A/m</math>
| |
| where:
| |
| :<math>a_R</math> is the magnitude of acceleration in meters per second-squared.
| |
| :<math>A</math> is the cross-sectional area exposed to the [[Sun]] in meters-squared.
| |
| :<math>m</math> is the spacecraft mass in [[kilograms]].
| |
| :<math>r</math> is the reflection factor which depends on material properties. <math>r=0</math> for absorption, <math>r=1</math> for specular reflection, and <math>r\approx0.4</math> for diffuse reflection.
| |
| | |
| For orbits around the Earth, solar radiation pressure becomes a stronger force than drag above 800 km altitude.<ref name="SMAD"/>
| |
| | |
| === Propulsion ===
| |
| {{main|Spacecraft propulsion}}
| |
| There are many different types of spacecraft propulsion. Rocket engines are one of the most widely used. The force of a rocket engine is modeled by the equation:<ref>{{cite book|author=George P. Sutton and Oscar Biblarz|title=Rocket Propulsion Elements|edition=7th|publisher=Wiley Interscience|year=2001|ISBN=0-471-32642-9}} See Equation 2-14.</ref>
| |
| :<math>F_n = \dot{m}\;v_\text{e} = \dot{m}\;v_\text{e-act} + A_\text{e}(p_\text{e} - p_\text{amb})</math>
| |
| | |
| :{| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right|<math>\dot{m}</math>
| |
| |align=left|= exhaust gas mass flow
| |
| |-
| |
| !align=right|<math>v_\text{e}</math>
| |
| |align=left|= effective exhaust velocity
| |
| |-
| |
| !align=right|<math>v_\text{e-act}</math>
| |
| |align=left|= actual jet velocity at nozzle exit plane
| |
| |-
| |
| !align=right|<math>A_\text{e}</math>
| |
| |align=left|= flow area at nozzle exit plane (or the plane where the jet leaves the nozzle if separated flow)
| |
| |-
| |
| !align=right|<math>p_\text{e}</math>
| |
| |align=left|= static pressure at nozzle exit plane
| |
| |-
| |
| !align=right|<math>p_\text{amb}</math>
| |
| |align=left|= ambient (or atmospheric) pressure
| |
| |}
| |
| | |
| Another possible method is a [[solar sail]]. Solar sails use [[radiation pressure]] in a way to achieve a desired propulsive force.<ref>{{cite web|url=http://messenger.jhuapl.edu/news_room/details.php?id=102 |title=MESSENGER Sails on Sun's Fire for Second Flyby of Mercury |date=2008-09-05 |quote=On September 4, the MESSENGER team announced that it would not need to implement a scheduled maneuver to adjust the probe's trajectory. This is the fourth time this year that such a maneuver has been called off. The reason? A recently implemented navigational technique that makes use of solar-radiation pressure (SRP) to guide the probe has been extremely successful at maintaining MESSENGER on a trajectory that will carry it over the cratered surface of Mercury for a second time on October 6. |postscript=<!--None-->}}</ref> The perturbation model due to the solar wind can be used as a model of propulsive force from a solar sail.
| |
| | |
| === Drag ===
| |
| {{main|Drag (physics)}}
| |
| The primary non-gravitational force acting on satellites in low Earth orbit is atmospheric drag.<ref name="SMAD"/> Drag will act in opposition to the direction of velocity and remove energy from an orbit. The force due to drag is modeled by the following equation:
| |
| | |
| :<math>F_D\, =\, \tfrac12\, \rho\, v^2\, C_d\, A,</math>
| |
| | |
| where
| |
| :<math> \mathbf{F}_D </math> is the [[force]] of drag,
| |
| :<math> \mathbf{} \rho </math> is the [[density]] of the fluid,<ref>Note that for the [[Earth's atmosphere]], the air density can be found using the [[barometric formula]]. It is 1.293 kg/m<sup>3</sup> at 0 °C and 1 [[atmosphere (unit)|atmosphere]].</ref>
| |
| :<math> \mathbf{} v </math> is the [[velocity]] of the object relative to the fluid,
| |
| :<math> \mathbf{} C_d </math> is the [[drag coefficient]] (a [[dimensionless number|dimensionless]] [[parameter]], e.g. 2 to 4 for most satellites<ref name="SMAD"/>)
| |
| :<math> \mathbf{} A </math> is the reference [[area]].
| |
| | |
| Orbits with an altitude below 120 km generally have such high drag that the orbits decay too rapidly to give a satellite a sufficient lifetime to accomplish any practical mission. On the other hand, orbits with an altitude above 600 km have relatively small drag so that the orbit decays slow enough that it has no real impact on the satellite over its useful life.<ref name="SMAD"/> [[Density of air]] can vary significantly in the [[thermosphere]] where most low Earth orbiting satellites reside. The variation is primarily due to solar activity, and thus solar activity can greatly influence the force of drag on a spacecraft and complicate long-term orbit simulation.<ref name="SMAD"/>
| |
| | |
| === Magnetic fields ===
| |
| Magnetic fields can play a significant role as a source of orbit perturbation as was seen in the [[Long Duration Exposure Facility]].<ref name="Roithmayr"/> Like gravity, the magnetic field of the Earth can be expressed through spherical harmonics as shown below:<ref name="Roithmayr"/>
| |
| | |
| :<math>{\mathbf{B}}=\sum_{n=1}^\infty \sum_{m=0}^n {\mathbf{B}}_{n,m}</math>
| |
| where
| |
| :<math>{\mathbf{B}}</math> is the magnetic field vector at a point above the Earth's surface.
| |
| :<math>{\mathbf{B}}_{n,m}</math> represents the contribution to <math>{\mathbf{B}}</math> of the spherical harmonic of degree ''n'' and order ''m'', defined as:<ref name="Roithmayr"/>
| |
| | |
| :<math>
| |
| \begin{align}
| |
| \mathbf{B}_{n,m} & = \frac{K_{n,m}a^{n+2}}{R^{n+m+1}}\left[\frac{g_{n,m}\mathcal{C}_m+h_{n,m}\mathcal{S}_m}{R}((s_{\lambda} A_{n,m+1}+(n+m+1)A_{n,m})\mathbf{\hat{r}})-A_{n,m+1}\mathbf{\hat{e}}_3\right] \\[10pt]
| |
| & {}\quad {}-mA_{n,m}((g_{n,m}\mathcal{C}_{m-1}+h_{n,m}\mathcal{S}_{m-1})\mathbf{\hat{e}}_1+(h_{n,m}\mathcal{C}_{m-1}-g_{n,m}\mathcal{S}_{m-1})\mathbf{\hat{e}}_2))
| |
| \end{align}
| |
| </math>
| |
| | |
| where:
| |
| :<math>a</math> is the mean equatorial radius of the primary body.
| |
| :<math>R</math> is the magnitude of the position vector from the center of the primary body to the center of the secondary body.
| |
| :<math>\mathbf{\hat{r}}</math> is a unit vector in the direction of the secondary body with its origin at the center of the primary body.
| |
| :<math> g_{n,m}</math> and <math>h_{n,m}</math> are Gauss coefficients of degree ''n'' and order ''m''. These are typically found through [[magnetic field]] measurements.
| |
| :The unit vectors <math>\mathbf{\hat{e}}_1, \mathbf{\hat{e}}_2, \mathbf{\hat{e}}_3</math> define a coordinate system fixed on the primary body. For the Earth, <math>\mathbf{\hat{e}}_1</math> lies in the equatorial plane parallel to a line intersecting Earth's geometric center and the [[Greenwich meridian]],<math>\mathbf{\hat{e}}_3</math> points in the direction of the North polar axis, and <math>\mathbf{\hat{e}}_2=\mathbf{\hat{e}}_3\times\mathbf{\hat{e}}_1</math>
| |
| :<math>A_{n,m}</math> is referred to as a derived [[Legendre polynomial]] of degree ''n'' and order ''m''. They are solved through the recurrence relation: <math>A_{n,m}(u)=\frac{1}{n-m}((2n-1)uA_{n-1,m}(u)-(n+m-1)A_{n-2,m}(u))</math>
| |
| :<math>K_{n,m}</math> is defined as: 1 if ''m'' = 0, <math>\big[\frac{n-m}{n+m}\big]^{0.5}K_{n-1,m}</math> for <math> n\ge (m+1)</math> and <math> m=[1\ldots\infty]</math> , and <math>[(n+m)(n-m+1)]^{-0.5}K_{n,m-1}</math> for <math> n\ge m</math> and <math>m=[2\ldots\infty]</math>
| |
| :<math>s_{\lambda}</math> is sine of the geographic latitude of the secondary body, which is <math>\mathbf{\hat{r}}\cdot\mathbf{\hat{e}}_3</math>.
| |
| :<math>\mathcal{C}_{m}, \mathcal{S}_{m}</math> are defined with the following recurrence relation and initial conditions:<math>\mathcal{C}_{m}=\mathcal{C}_{1}\mathcal{C}_{m-1}-\mathcal{S}_{1}\mathcal{S}_{m-1}, \mathcal{S}_{m}=\mathcal{S}_{1}\mathcal{C}_{m-1}+\mathcal{C}_{1}\mathcal{S}_{m-1}, \mathcal{S}_{0}=0, \mathcal{S}_{1}=\mathbf{R}\cdot\mathbf{\hat{e}}_2, \mathcal{C}_{0}=1, \mathcal{C}_{1}=\mathbf{R}\cdot\mathbf{\hat{e}}_1</math>
| |
| | |
| == See also ==
| |
| * [[Perturbation (astronomy)]]
| |
| * [[Osculating orbit]]
| |
| * [[Orbital resonance]]
| |
| * [[n-body problem]]
| |
| * [[two-body problem]]
| |
| * [[sphere of influence (astrodynamics)]]
| |
| | |
| == Notes and references ==
| |
| <references/>
| |
| | |
| ==External links==
| |
| * [http://www.fas.org/irp/imint/docs/rst/Intro/Part2_1b.html] Gravity maps of the Earth
| |
| | |
| [[Category:Orbital perturbations]]
| |
| [[Category:Dynamical systems]]
| |
| [[Category:Dynamics of the Solar System]]
| |