手机浏览器扫描二维码访问
作者君在作品相关中其实已经解释过这个问题。
不过仍然有人质疑。
那么作者君在此列出相关参考文献中的一篇开源论文。
以下是文章内容:
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).
凌天神帝 农门悍妻太嚣张 天价妈咪:爹地闪开宝宝来 宠妻成狂:闪婚总裁太霸道 极品全能狂少 我家太子妃要上天 修真聊天群 绝品神相 都市逍遥兵王 穿越六十年代农家女 医鸣惊人:残王独宠废材妃 狂武战尊 三国之蜀汉中兴 都市最强高手 都市逍遥邪医 美女总裁的神级兵王 大周王侯 龙王之我是至尊 我从天界归来 医手遮天
穿越到了大明万历皇帝的身上,你们将会听到一个声音万历驾到!简介无能,就是穿越附身万历做皇帝的故事。...
谁说穿越后宫就是圣宠不衰六宫无妃的幸福生涯?她保证不打死他!过来,朕不打你。放屁!渺渺,过来朕抱。谁信!苏渺,再不过来打断你的狗腿!皇上~~~苏渺一睁眼就是冷宫的四堵灰墙,简直凄凄惨惨戚戚!为了保住她的腿,只好狗腿的腻在皇帝陛下身边,惹得众人纷纷嫉恨白眼,直骂妖妃。可惜皇帝陛下不知怎的就宠...
工作不顺心,女朋友闹分手,房租交不上霉运连连的方少俊想死的心都有哎?这是什么?蓝药水,喷了方少俊一脸从此,人生,逆转!...
我只是一名普通的商人,我的职责,就是将两个时空的商品相互贩卖。末世的食粮,如同金子一般珍贵,这个当然不能放手。现代里稀缺的金属?末世里随便你去捡去拉,拉到...
医生,居全球最受尊重的职业排行首位,高学历高智商高收入。故事的主人公曾经也这么认为,而真正做了医生才知道,这是一个高投入高强度高压力高风险的职业。当下,医患纠纷呈不断上升之势,医患关系已经紧张到了白热化的程度。曾经被赞美称医者仁心的医生们,处于人人自危又人人自保的两难中,这种尴尬的环境又加剧了患者看病难看病贵的现状。在这种情况规培医生政策运应而生。2015年之前,是否规培规培年限规培方式在不同地区同一地区不同医院各有差异,不少医院借此机会让医学生成为了廉价劳动力...
书海阁小说网免费提供作者落鱼的经典小说天武僵神最新章节全文阅读服务本站更新及时无弹窗广告欢迎光临wwwshgtw观看小说夜洛重生天武大陆,成为了曾经叱咤整个大陆的超级种族僵尸一族的普通僵尸。时过万载,僵尸一族再次出世,作为普通僵尸的夜洛,在人族妖族围剿之下一步步成长,逆天地铸神位,成为一代僵尸之神!!!...