Chislennoe modelirovanie - vazhneishii instrument v kosmologicheskih issledovaniyah. Sushestvuet mnogo razlichnyh podhodov k "sozdaniyu vselennoi v probirke". Sootvetstvenno, neobhodimo sravnivat' raznye podhody monopolizm i otsutstvie svobodnogo kompetentnogo obsuzhdeniya nauke protivopokazano).

Avtory (izvestneishie specialisty v svoei oblasti) sravnivayut dva chislennyh koda, osnovannyh na raznyh metodikah. Obsuzhdayutsya prichiny neskol'ko razlichayushihsya rezul'tatov raboty dvuh programm.

Daetsya obzor novoi versii parallel'nyh superkomp'yuterov serii GRAPE (GRAvity piPE) dlya modelirovaniya metodom mnogih tel, aktivno primenyayushimsya v astrofizike (struktura i vzaimodeistvie galaktik, zvezdnye skopleniya, krupnomasshtabnaya struktura, ustoichivost' i formirovanie planetnyh sistem i t.p.). Seriya nachalas' v nachale 90-h godov i, razumeetsya, moshnost' komp'yuterov bystro rosla. Komp'yutery GRAPE-3 i -5 proizvodilis' dazhe malymi seriyami, chto pozvolilo sdelat' ih ne ochen' dorogimi (t.e. dostupnymi bol'shinstvu krupnyh universitetov, gde byli sootvetstvuyushie gruppy issledovatelei). GRAPE-6 imeet 2048 processorov. Maksimal'naya skorost' - 64 Teraflopa.

Stat'ya budet bezuslovno interesna vsem, kto interesuetsya chislennymi metodami, komp'yuternym modelirovanie i t.p.

Korotkii otchet o hode rabot nad interesnym proektom. Avotry postavili svoei zadachei s pomosh'yu chislennogo modelirovaniya poluchit' biblioteku slivayushihsya galaktik. Kak horosho izvestno, vzaimodeistvie galaktik mozhet porozhdat' ves'ma prichudlivye konfiguracii. Bylo by interesno zaranee imet' nabor simulyacii. Dlya togo, chtoby eto davalo otdachu, nuzhno imet' okolo 1000 smodelirovannyh sliyanii. Bol'shaya rabota!

Sushestvuet dostatochno mnogo otkrytyh kodov dlya samyh razlichnyh celei. V tom chisle i dlya kosmologicheskogo modelirovaniya. V dannoi stat'e avtory predstavlyayut opisanie programmy FLY. Rabota, razumeetsya, ochen' tehnicheskaya, no budet interesna tem, kto interesuetsya chislennym modelirovaniem.

Dve stat'i (vtoraya "Smoothed Particle Magnetohydrodynamics II. Variational principles and variable smoothing length terms"), posvyashennye modelirovaniyu metodom sglazhennyh chastic (Smoothed Particle Hydrodynamics (SPH)).

Predlozhen novyi algoritm, kotoryi proveren mnozhestvom odnomernyh testov, predstavlennymi v etih dvuh rabotah. Predpolagaetsya, chto algoritm budet kraine polezen v astrofizicheskih prilozheniyah (vklyuchaya neideal'nuyu MGD).

Dissertaciya posvyashena chislennomu modelirovaniya gravitacionnogo kollapsa. Tema arhislozhnaya i arhiaktual'naya.

Eta stat'ya - prazdnik trehmernyh raschetov, kotorye byli vypolneny dlya serii modelei s razlichnoi skorost'yu vrasheniya i magnitnym polem. Samoe bystroe iz rassmotrennyh vrashenii ne privodit k fragmentacii zvezdy i ne vyzyvaet sushestvennyh izmenenii v processah konvekcii v centre ee yadra. Analogichno, rassmatrivalis' ne slishkom sil'nye magnitnye polya, kotorye ne dominiruyut v processe vzryva sverhnovoi (bolee sil'nye magnitnye polya, nablyudaemye u neitronnyh zvezd, mogut voznikat' pri ih posleduyushem ohlazhdenii i szhatii). V modeli s samym bystrym vrasheniem mozhet obrazovat'sya pul'sar, energii vrasheniya kotorogo budet dostatochno dlya sbrosa obolochki sverhnovoi. No sovremennye evolyucionnye teorii predskazyvayut slishkom medlennoe vrashenie u predsverhnovyh, chtoby takoi (magnitorotacionnyi) mehanizm mog srabotat'.

Pravil'no razdelit' nablyudatel'nye dannye na intervaly, chtoby zatem ih udachno approksimirovat' - pochti iskusstvo. V dannoi zametke opisan prostoi, no moshnyi algoritm, kotoryi prosmatrivaet vse vozmozhnye varianty razbieniya N tochek dannyh za vremya poryadka N² (obrashaem vnimanie, chto obshee chislo variantov razbieniya proporcional'no e^N). Algoritm garantirovanno nahodit global'nyi optimum i avtomaticheski opredelyaet chislo intervalov na kotoroe nado razbit' ishodnye dannye. Algoritm legko obobshaetsya na mnogomernye dannye. Vozmozhna, takzhe, ego rabota v real'nom vremeni.

Primer razbieniya nabora dannyh na intervaly

Bazy dannyh bol'shih obzorov neba (naprimer, SuperCOSMOS Sky Survey, SSS), takzhe kak i otdel'nye nablyudeniya, "zasoryayutsya" proletami iskusstvennyh sputnikov i samoletov, fizicheskimi defektami priemnikov (segodnya eto, v osnovnom, PZS) i t.d. Odnako, esli poyavlenie podobnyh defektov v odinochnyh nablyudeniyah mozhno obnaruzhit' i, po vozmozhnosti, ustranit' v ruchnuyu, to v avtomaticheski provodimyh obzorah bol'shogo ob'ema takoe reshenie nevozmozhno.

Dannaya stat'ya posvyashena obnaruzheniyu i udaleniyu iz dannyh obzorov "trekov" iskusstvennyh sputnikov zemli i samoletov i drugih lineinyh defektov.

Polnuyu versiyu stat'i (s risunkami) mozhno vzyat' na saite avtorov http://www.anc.ed.ac.uk/~amos/publications.html, a detali proekta uznat' zdes'.

Baza dannyh o galaktikah "GOLDmine" ("Zolotye kopi") byla sushestvenno obnovlena k sentyabryu 2003 g. Adres bazy dannyh: http://goldmine.mib.infn.it/.

Skromnoe nazvanie sistemy MODEST ne dolzhno vvodit' v zabluzhdenie. Eto krupnyi proekt, kotoryi daet vozmozhnost' provodit' kompleksnoe modelirovanie plotnyh zvezdnyh sistem, naprimer, sharovyh skoplenii. Nekotorye kuski koda dostupny v Internet. Proekt ob'edinyaet neskol'ko issledovatel'skih grupp po vsemu miru. V stat'e daetsya obzor pervogo goda raboty.

Avtor stat'i rabotaet v Department of Energy (chto-to tipa nashego MinSredMasha) uzhe davno. Yasno, chto imenno takie vedomstva imeyut vozmozhnost' provodit' ochen' slozhnye chislennye eksperimenty.

Stat'ya ochen' i ochen' interesnaya. Vsem sovetuem prochest'. Opisano mnogo konkretnyh rabot (osobenno v oblasti izucheniya povedeniya haoticheskih sistem).

Interesno, "a kak u nas?".....

Bol'shoe kolichestvo astronomicheskih ob'ektov vybrasyvaet strui (ili, kak ih eshe nazyvayut, dzhety). Eto i obrazuyushiesya zvezdy, i kompaktnye ob'ekty v tesnyh dvoinyh sistemah, i sverhmassivnye chernye dyry v yadrah galaktik.... Esli skorost' dzheta velika, to ih nazyvayut relyativistskimi, t.k. nuzhno uchityvat' effekty special'noi teorii otnositel'nosti. Krome togo, vo mnogih dzhetah nablyudayut zakruchennye (spiral'nye) struktury. Ob'yasneniyu etogo fenomena, a takzhe popytkam vosproizvesti podobnye struktury v chislennyh eksperimentah, posvyashena dannaya rabota.

Eto ne pervyi raschet sliyaniya dvoinoi neitronnoi zvezdy (i uzh tem bolee ne poslednii). Osnovnaya problema zdes' sostoit v sleduyushem provodit' etot raschet mozhno tol'ko v treh izmereniyah (poskol'ku zadacha ne obladaet osevoi simmetriei), a esli, kak v dannom sluchae, uchityvayutsya magnitnye polya, to trebovaniya k ob'emam komp'yuternoi pamyati stanovyatsya ochen' bol'shimi. Bol'shinstvo predydushih raschetov velos' na malyh setkah (s nizkim razresheniem), imenno poetomu kazhdoe sushestvennoe uvelichenie razresheniya v takih raschetah vyzyvaet bol'shoi interes.

Chislennym raschetam sliyaniya dvoinyh neitronnyh zvezd i izluchayushihsya pri etom gravitacionnyh voln posvyashen takzhe preprint astro-ph/0306481 k kotoromu prilagaetsya fil'm.

V etoi stat'e rassmatrivaetsya i kolichestvenno ocenivaetsya, s pomosh'yu chislennogo modelirovaniya metodom N-tel, nablyudayushiisya vklad dvuh effektov s ochen' poeticheskimi nazvaniyami "Pal'cy Boga" (Fingers of God) i "effekt Bych'ego Glaza" (Bull's Eye Effect). Oba etih effekta svyazany s otkloneniyami dvizhenii otdel'nyh ob'ektov (galaktik) ot zakona Habbla. Eti otkloneniya svyazany s tem, chto krome uchastiya v obshem Habblovskom rasshirenii imeyut sobstvennye pekulyarnye skorosti (otnositel'no dvizheniya Habblovskogo potoka v srednem).

Pervyi effekt svyazan s dvizheniyami galaktik v kompaktnyh stacionarnyh skopleniyah. Virial'nye skorosti ob'ektov tam dostigayut tysyach kilometrov v sekundu, polovina galaktik dvizhetsya k nam (ih krasnoe smeshenie z umen'shaetsya), vtoraya polovina - ot nac (z stanovitsya bol'she). Togda pri postroenii trehmernogo raspredeleniya ob'ektov, kogda rasstoyanie vychislyaetsya cherez luchevuyu skorost' (krasnoe smeshenie) po zakonu Habbla, v napravlenii virializovannyh skoplenii my uvidim vytyanutye vdol' lucha zreniya "sigary" - "'Pal'cy Boga' so vseh storon ukazyvayushie na nas".

Vtoroi effekt ("effekt Bych'ego Glaza") rabotaet na bol'shih masshtabah i deistvuet v protivopolozhnuyu storonu (t.e. ne rastyagivaet, a szhimaet raspredelenie ob'ektov po z). On deistvuet v bolee krupnyh strukturah (skopleniyah i sverh skopleniyah), kotorye ne dostigli stacionarnogo sostoyaniya i vse eshe prodolzhayut kollapsirovat'. V etom sluchae ob'ekty raspolozhennye na blizhnei k nas storone takogo ob'ekta udalyayutsya ot nas v srednem bystree, chem po zakonu Habbla, a na dal'nei storone - v srednem medlennee. V itoge takoi ob'ekt po krasnomu smesheniyu zanimaet men'shii interval, chem sootvetstvuyushii ego geometricheskim razmeram. K sozhaleniyu, proishozhdenie nazvaniya etogo effekta ne ponyatno.

Illyustraciya vozniknoveniya effektov "Pal'cev Boga" i "Bych'ego Glaza".

Metody raschetov nuzhno postoyanno uluchshat'. V stat'e avtory obsuzhdayut problemu postanovki pravil'nyh nachal'nyh uslovii pri modelirovanii v kosmologii. Delo v tom, chto struktura poyavlyaetsya iz uzhe sushestvovavshih fluktuacii plotnosti. Vazhno pravil'no ih zadat'. Avtory delayut neskol'ko "progonov" s raznymi parametrami, i dayut svoi rekomendacii po vyboru nachal'nyh uslovii i rabote s nimi.

^{(Arhiv Chislennoe modelirovanie:} v.2, 2004, v.1, 2002-2003)