手机浏览器扫描二维码访问
作者君在作品相关中其实已经解释过这个问题。
不过仍然有人质疑。
那么作者君在此列出相关参考文献中的一篇开源论文。
以下是文章内容:
Long-termintegrationsandstabilityofplanetaryorbitsinourSolarsystem
Abstract
Wepresenttheresultsofverylong-termnumericalintegrationsofplanetaryorbitalmotionsover109-yrtime-spansincludingallnineplanets.Aquickinspectionofournumericaldatashowsthattheplanetarymotion,atleastinoursimpledynamicalmodel,seemstobequitestableevenoverthisverylongtime-span.Acloserlookatthelowest-frequencyoscillationsusingalow-passfiltershowsusthepotentiallydiffusivecharacterofterrestrialplanetarymotion,especiallythatofMercury.ThebehaviouroftheeccentricityofMercuryinourintegrationsisqualitativelysimilartotheresultsfromJacquesLaskarssecularperturbationtheory(e.g.emax~0.35over~±4Gyr).However,therearenoapparentsecularincreasesofeccentricityorinclinationinanyorbitalelementsoftheplanets,whichmayberevealedbystilllonger-termnumericalintegrations.Wehavealsoperformedacoupleoftrialintegrationsincludingmotionsoftheouterfiveplanetsoverthedurationof±5×1010yr.TheresultindicatesthatthethreemajorresonancesintheNeptune–Plutosystemhavebeenmaintainedoverthe1011-yrtime-span.
1Introduction
1.1Definitionoftheproblem
ThequestionofthestabilityofourSolarsystemhasbeendebatedoverseveralhundredyears,sincetheeraofNewton.Theproblemhasattractedmanyfamousmathematiciansovertheyearsandhasplayedacentralroleinthedevelopmentofnon-lineardynamicsandchaostheory.However,wedonotyethaveadefiniteanswertothequestionofwhetherourSolarsystemisstableornot.Thisispartlyaresultofthefactthatthedefinitionoftheterm‘stability’isvaguewhenitisusedinrelationtotheproblemofplanetarymotionintheSolarsysteActuallyitisnoteasytogiveaclear,rigorousandphysicallymeaningfuldefinitionofthestabilityofourSolarsyste
Amongmanydefinitionsofstability,hereweadopttheHilldefinition(Gladman1993):actuallythisisnotadefinitionofstability,butofinstability.Wedefineasystemasbecomingunstablewhenacloseencounteroccurssomewhereinthesystem,startingfromacertaininitialconfiguration(Chambers,Wetherill&Boss1996;Ito&Tanikawa1999).AsystemisdefinedasexperiencingacloseencounterwhentwobodiesapproachoneanotherwithinanareaofthelargerHillradius.Otherwisethesystemisdefinedasbeingstable.HenceforwardwestatethatourplanetarysystemisdynamicallystableifnocloseencounterhappensduringtheageofourSolarsystem,about±5Gyr.Incidentally,thisdefinitionmaybereplacedbyoneinwhichanoccurrenceofanyorbitalcrossingbetweeneitherofapairofplanetstakesplace.Thisisbecauseweknowfromexperiencethatanorbitalcrossingisverylikelytoleadtoacloseencounterinplanetaryandprotoplanetarysystems(Yoshinaga,Kokubo&Makino1999).OfcoursethisstatementcannotbesimplyappliedtosystemswithstableorbitalresonancessuchastheNeptune–Plutosyste
1.2Previousstudiesandaimsofthisresearch
Inadditiontothevaguenessoftheconceptofstability,theplanetsinourSolarsystemshowacharactertypicalofdynamicalchaos(Sussman&Wisdom1988,1992).Thecauseofthischaoticbehaviourisnowpartlyunderstoodasbeingaresultofresonanceoverlapping(Murray&Holman1999;Lecar,Franklin&Holman2001).However,itwouldrequireintegratingoveranensembleofplanetarysystemsincludingallnineplanetsforaperiodcoveringseveral10Gyrtothoroughlyunderstandthelong-termevolutionofplanetaryorbits,sincechaoticdynamicalsystemsarecharacterizedbytheirstrongdependenceoninitialconditions.
Fromthatpointofview,manyofthepreviouslong-termnumericalintegrationsincludedonlytheouterfiveplanets(Sussman&Wisdom1988;Kinoshita&Nakai1996).Thisisbecausetheorbitalperiodsoftheouterplanetsaresomuchlongerthanthoseoftheinnerfourplanetsthatitismucheasiertofollowthesystemforagivenintegrationperiod.Atpresent,thelongestnumericalintegrationspublishedinjournalsarethoseofDuncan&Lissauer(1998).Althoughtheirmaintargetwastheeffectofpost-main-sequencesolarmasslossonthestabilityofplanetaryorbits,theyperformedmanyintegrationscoveringupto~1011yroftheorbitalmotionsofthefourjovianplanets.TheinitialorbitalelementsandmassesofplanetsarethesameasthoseofourSolarsysteminDuncan&Lissauerspaper,buttheydecreasethemassoftheSungraduallyintheirnumericalexperiments.Thisisbecausetheyconsidertheeffectofpost-main-sequencesolarmasslossinthepaper.Consequently,theyfoundthatthecrossingtime-scaleofplanetaryorbits,whichcanbeatypicalindicatoroftheinstabilitytime-scale,isquitesensitivetotherateofmassdecreaseoftheSun.WhenthemassoftheSunisclosetoitspresentvalue,thejovianplanetsremainstableover1010yr,orperhapslonger.Duncan&Lissaueralsoperformedfoursimilarexperimentsontheorbitalmotionofsevenplanets(VenustoNeptune),whichcoveraspanof~109yr.Theirexperimentsonthesevenplanetsarenotyetcomprehensive,butitseemsthattheterrestrialplanetsalsoremainstableduringtheintegrationperiod,maintainingalmostregularoscillations.
Ontheotherhand,inhisaccuratesemi-analyticalsecularperturbationtheory(Laskar1988),Laskarfindsthatlargeandirregularvariationscanappearintheeccentricitiesandinclinationsoftheterrestrialplanets,especiallyofMercuryandMarsonatime-scaleofseveral109yr(Laskar1996).TheresultsofLaskarssecularperturbationtheoryshouldbeconfirmedandinvestigatedbyfullynumericalintegrations.
Inthispaperwepresentpreliminaryresultsofsixlong-termnumericalintegrationsonallnineplanetaryorbits,coveringaspanofseveral109yr,andoftwootherintegrationscoveringaspanof±5×1010yr.Thetotalelapsedtimeforallintegrationsismorethan5yr,usingseveraldedicatedPCsandworkstations.Oneofthefundamentalconclusionsofourlong-termintegrationsisthatSolarsystemplanetarymotionseemstobestableintermsoftheHillstabilitymentionedabove,atleastoveratime-spanof±4Gyr.Actually,inournumericalintegrationsthesystemwasfarmorestablethanwhatisdefinedbytheHillstabilitycriterion:notonlydidnocloseencounterhappenduringtheintegrationperiod,butalsoalltheplanetaryorbitalelementshavebeenconfinedinanarrowregionbothintimeandfrequencydomain,thoughplanetarymotionsarestochastic.Sincethepurposeofthispaperistoexhibitandoverviewtheresultsofourlong-termnumericalintegrations,weshowtypicalexamplefiguresasevidenceoftheverylong-termstabilityofSolarsystemplanetarymotion.Forreaderswhohavemorespecificanddeeperinterestsinournumericalresults,wehavepreparedawebpage(access),whereweshowraworbitalelements,theirlow-passfilteredresults,variationofDelaunayelementsandangularmomentumdeficit,andresultsofoursimpletime–frequencyanalysisonallofourintegrations.
InSection2webrieflyexplainourdynamicalmodel,numericalmethodandinitialconditionsusedinourintegrations.Section3isdevotedtoadescriptionofthequickresultsofthenumericalintegrations.Verylong-termstabilityofSolarsystemplanetarymotionisapparentbothinplanetarypositionsandorbitalelements.Aroughestimationofnumericalerrorsisalsogiven.Section4goesontoadiscussionofthelongest-termvariationofplanetaryorbitsusingalow-passfilterandincludesadiscussionofangularmomentumdeficit.InSection5,wepresentasetofnumericalintegrationsfortheouterfiveplanetsthatspans±5×1010yr.InSection6wealsodiscussthelong-termstabilityoftheplanetarymotionanditspossiblecause.
2Descriptionofthenumericalintegrations
(本部分涉及比较复杂的积分计算,作者君就不贴上来了,贴上来了起点也不一定能成功显示。)
2.3Numericalmethod
Weutilizeasecond-orderWisdom–Holmansymplecticmapasourmainintegrationmethod(Wisdom&Holman1991;Kinoshita,Yoshida&Nakai1991)withaspecialstart-upproceduretoreducethetruncationerrorofanglevariables,‘warmstart’(Saha&Tremaine1992,1994).
Thestepsizeforthenumericalintegrationsis8dthroughoutallintegrationsofthenineplanets(N±1,2,3),whichisabout111oftheorbitalperiodoftheinnermostplanet(Mercury).Asforthedeterminationofstepsize,wepartlyfollowthepreviousnumericalintegrationofallnineplanetsinSussman&Wisdom(1988,7.2d)andSaha&Tremaine(1994,22532d).Weroundedthedecimalpartofthetheirstepsizesto8tomakethestepsizeamultipleof2inordertoreducetheaccumulationofround-offerrorinthecomputationprocesses.Inrelationtothis,Wisdom&Holman(1991)performednumericalintegrationsoftheouterfiveplanetaryorbitsusingthesymplecticmapwithastepsizeof400d,110.83oftheorbitalperiodofJupiter.Theirresultseemstobeaccurateenough,whichpartlyjustifiesourmethodofdeterminingthestepsize.However,sincetheeccentricityofJupiter(~0.05)ismuchsmallerthanthatofMercury(~0.2),weneedsomecarewhenwecomparetheseintegrationssimplyintermsofstepsizes.
Intheintegrationoftheouterfiveplanets(F±),wefixedthestepsizeat400d.
WeadoptGaussfandgfunctionsinthesymplecticmaptogetherwiththethird-orderHalleymethod(Danby1992)asasolverforKeplerequations.ThenumberofmaximumiterationswesetinHalleysmethodis15,buttheyneverreachedthemaximuminanyofourintegrations.
Theintervalofthedataoutputis200000d(~547yr)forthecalculationsofallnineplanets(N±1,2,3),andabout8000000d(~21903yr)fortheintegrationoftheouterfiveplanets(F±).
Althoughnooutputfilteringwasdonewhenthenumericalintegrationswereinprocess,weappliedalow-passfiltertotheraworbitaldataafterwehadcompletedallthecalculations.SeeSection4.1formoredetail.
2.4Errorestimation
2.4.1Relativeerrorsintotalenergyandangularmomentum
Accordingtooneofthebasicpropertiesofsymplecticintegrators,whichconservethephysicallyconservativequantitieswell(totalorbitalenergyandangularmomentum),ourlong-termnumericalintegrationsseemtohavebeenperformedwithverysmallerrors.Theaveragedrelativeerrorsoftotalenergy(~10?9)andoftotalangularmomentum(~10?11)haveremainednearlyconstantthroughouttheintegrationperiod(Fig.1).Thespecialstartupprocedure,warmstart,wouldhavereducedtheaveragedrelativeerrorintotalenergybyaboutoneorderofmagnitudeormore.
RelativenumericalerrorofthetotalangularmomentumδAA0andthetotalenergyδEE0inournumericalintegrationsN±1,2,3,whereδEandδAaretheabsolutechangeofthetotalenergyandtotalangularmomentum,respectively,andE0andA0aretheirinitialvalues.ThehorizontalunitisGyr.
Notethatdifferentoperatingsystems,differentmathematicallibraries,anddifferenthardwarearchitecturesresultindifferentnumericalerrors,throughthevariationsinround-offerrorhandlingandnumericalalgorithms.IntheupperpanelofFig.1,wecanrecognizethissituationinthesecularnumericalerrorinthetotalangularmomentum,whichshouldberigorouslypreserveduptomachine-εprecision.
2.4.2Errorinplanetarylongitudes
SincethesymplecticmapspreservetotalenergyandtotalangularmomentumofN-bodydynamicalsystemsinherentlywell,thedegreeoftheirpreservationmaynotbeagoodmeasureoftheaccuracyofnumericalintegrations,especiallyasameasureofthepositionalerrorofplanets,i.e.theerrorinplanetarylongitudes.Toestimatethenumericalerrorintheplanetarylongitudes,weperformedthefollowingprocedures.Wecomparedtheresultofourmainlong-termintegrationswithsometestintegrations,whichspanmuchshorterperiodsbutwithmuchhigheraccuracythanthemainintegrations.Forthispurpose,weperformedamuchmoreaccurateintegrationwithastepsizeof0.125d(164ofthemainintegrations)spanning3×105yr,startingwiththesameinitialconditionsasintheN?1integration.Weconsiderthatthistestintegrationprovidesuswitha‘pseudo-true’solutionofplanetaryorbitalevolution.Next,wecomparethetestintegrationwiththemainintegration,N?1.Fortheperiodof3×105yr,weseeadifferenceinmeananomaliesoftheEarthbetweenthetwointegrationsof~0.52°(inthecaseoftheN?1integration).Thisdifferencecanbeextrapolatedtothevalue~8700°,about25rotationsofEarthafter5Gyr,sincetheerroroflongitudesincreaseslinearlywithtimeinthesymplecticmap.Similarly,thelongitudeerrorofPlutocanbeestimatedas~12°.ThisvalueforPlutoismuchbetterthantheresultinKinoshita&Nakai(1996)wherethedifferenceisestimatedas~60°.
3Numericalresults–I.Glanceattherawdata
Inthissectionwebrieflyreviewthelong-termstabilityofplanetaryorbitalmotionthroughsomesnapshotsofrawnumericaldata.Theorbitalmotionofplanetsindicateslong-termstabilityinallofournumericalintegrations:noorbitalcrossingsnorcloseencountersbetweenanypairofplanetstookplace.
3.1Generaldescriptionofthestabilityofplanetaryorbits
First,webrieflylookatthegeneralcharacterofthelong-termstabilityofplanetaryorbits.Ourinterestherefocusesparticularlyontheinnerfourterrestrialplanetsforwhichtheorbitaltime-scalesaremuchshorterthanthoseoftheouterfiveplanets.AswecanseeclearlyfromtheplanarorbitalconfigurationsshowninFigs2and3,orbitalpositionsoftheterrestrialplanetsdifferlittlebetweentheinitialandfinalpartofeachnumericalintegration,whichspansseveralGyr.Thesolidlinesdenotingthepresentorbitsoftheplanetsliealmostwithintheswarmofdotseveninthefinalpartofintegrations(b)and(d).Thisindicatesthatthroughouttheentireintegrationperiodthealmostregularvariationsofplanetaryorbitalmotionremainnearlythesameastheyareatpresent.
Verticalviewofthefourinnerplanetaryorbits(fromthez-axisdirection)attheinitialandfinalpartsoftheintegrationsN±1.Theaxesunitsareau.Thexy-planeissettotheinvariantplaneofSolarsystemtotalangularmomentu(a)TheinitialpartofN+1(t=0to0.0547×109yr).(b)ThefinalpartofN+1(t=4.9339×108to4.9886×109yr).(c)TheinitialpartofN?1(t=0to?0.0547×109yr).(d)ThefinalpartofN?1(t=?3.9180×109to?3.9727×109yr).Ineachpanel,atotalof23684pointsareplottedwithanintervalofabout2190yrover5.47×107yr.Solidlinesineachpaneldenotethepresentorbitsofthefourterrestrialplanets(takenfromDE245).
极品全能狂少 天价妈咪:爹地闪开宝宝来 三国之蜀汉中兴 都市最强高手 穿越六十年代农家女 宠妻成狂:闪婚总裁太霸道 医鸣惊人:残王独宠废材妃 凌天神帝 修真聊天群 大周王侯 狂武战尊 龙王之我是至尊 我从天界归来 都市逍遥邪医 农门悍妻太嚣张 美女总裁的神级兵王 绝品神相 我家太子妃要上天 医手遮天 都市逍遥兵王
作者暖暖的经典小说甜宠霸爱重生娇妻吻上瘾最新章节全文阅读服务本站更新及时无弹窗广告小说甜宠霸爱重生娇妻吻上瘾前世因为识人不清,一生落魄惨死。重活一世,她只想逆转翻身,不再重蹈覆辙。重回18岁生日宴的前夕,她无奈向他求助,却不料,他单手掐着她的小脸儿。真是个可爱又甜美的小女人。...
上一世,瞎了眼!真情错付,死不瞑目!这一世,再重来,仇恨滔天!绝不手软!惩庶母,虐渣男,斗小三!专治各种不服和小人!相府嫡女,涅槃重生!当最善良的人学会狠毒,当最软弱的人硬起心肠,凤眸半挑,素手倾城,欺她害她者,定要千百倍偿还!翻云覆雨中,将自己的命运掌握在自己的手上!且看她锋芒喋血,风华惊天下!...
爽文当道,欲逆世而行!圣纹世界所流传的,是名扬天下流传千古的传奇佳话,而非天选之子气运皇帝背景板!强者,非朝夕!天赋,非恩赐!运气,非资本!世事所不容者,皆有挣扎化朽为奇迹之能!欲由心生,此为贪婪!展开收起...
作者今歌的经典小说总裁蜜宠小萌妻最新章节全文阅读服务本站更新及时无弹窗广告小说总裁蜜宠小萌妻那男的刚说什么?宗继泽黑着脸问。助理呃,他他说,再看,再看我就吃了你?宗继泽的脸更黑了,这男人,找死!陆丁宁不退反进,一把墙咚,嘿,帅哥,听说你对女人不敢兴趣,要不咋两聊聊滚!宗继泽俊脸黑了又白,白了又红。某天,宗继泽被吃干抹净后宗继泽咬牙,这女人,死定了!...
麦成雄满意的站在媒体镜头前大声说道我并不是你们所说的神,更不是香港梦的缔造者,我能做的无非就是指引了一个方向。我是80后,看着港片长大的,不知道什么时候,我离它越来越远了,已经好久没有进电影院,就像工作之后,远离故乡,但这份情难以割舍,虽然现在是好莱坞大片当道,不过港片一直在我心中,熟悉港片历史的人都知道港片曾经的辉煌,即使好莱坞影人也会为港片喝彩,假如环境改变,谁敢说香港梦不为人乐道?谢谢马桶清洁工提供的书友群119609144...
男神他又套路我前世苏倩倩一心学习,被母亲和妹妹玩弄于股掌之中,到最后连命都送了!意外重生。重活一世,虐渣渣,谈恋爱,成为人生最大的赢家!唯独婚姻是最不想要的,可是傅少,您别缠着我啊,我哪里好,我跑不掉改还不行吗?...