Long-Time-Step Methods for Oscillatory Differential Equations
Considered are numerical integration schemes for nondissipative dynamical systems in which multiple time scales are present. It is assumed that one can do an explicit separation of the RHS "forces" into fast forces and slow forces such that (i) the fast forces contain the high frequency pa...
Gespeichert in:
Veröffentlicht in: | SIAM journal on scientific computing 1998, Vol.20 (3), p.930-963 |
---|---|
Hauptverfasser: | , , |
Format: | Artikel |
Sprache: | eng |
Schlagworte: | |
Online-Zugang: | Volltext |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
container_end_page | 963 |
---|---|
container_issue | 3 |
container_start_page | 930 |
container_title | SIAM journal on scientific computing |
container_volume | 20 |
creator | García-Archilla, B. Sanz-Serna, J. M. Skeel, R. D. |
description | Considered are numerical integration schemes for nondissipative dynamical systems in which multiple time scales are present. It is assumed that one can do an explicit separation of the RHS "forces" into fast forces and slow forces such that (i) the fast forces contain the high frequency part of the solution, (ii) the fast forces are conservative, and (iii) the reduced problem consisting only of the fast forces can be integrated much more cheaply than the full problem. The fast forces are allowed to have low frequency components. Particular applications for which the schemes are intended include N-body problems (for which most of the forces are slow) and nonlinear wave phenomena (for which the fast forces can be propagated by spectral methods). The assumption of cheap integration of fast forces implies that the overall cost of integration is primarily determined by the step size used to sample the slow forces. A long-time-step method is one in which this step size exceeds half the period of the fastest normal mode present in the full system. An existing method that comes close to qualifying is the "impulse" method, also known as Verlet-I and r-RESPA. It is shown that it might fail, though, for a couple of reasons. First, it suffers a serious loss of accuracy if the step size is near a multiple of the period of a normal mode, and, second, it is unstable if the step size is near a multiple of half the period of a normal mode. Proposed in this paper is a "mollified" impulse method having an error bound that is independent of the frequency of the fast forces. It is also shown to possess superior stability properties. Theoretical results are supplemented by numerical experiments. The method is efficient and reasonably easy to implement. |
doi_str_mv | 10.1137/S1064827596313851 |
format | Article |
fullrecord | <record><control><sourceid>proquest_cross</sourceid><recordid>TN_cdi_proquest_miscellaneous_25821952</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><sourcerecordid>2586040781</sourcerecordid><originalsourceid>FETCH-LOGICAL-c367t-124d40b063190d9bfe40dbe065075923ba8f4c12876ed14316be80c1abc0085b3</originalsourceid><addsrcrecordid>eNplkDtPwzAUhS0EEqXwA9giBrbAvX7FHhgQlIdU1KFltuzEgVRp3NrJ0H9PqjLBdK50Pl19OoRcI9whsuJ-iSC5ooXQkiFTAk_IBEGLvEBdnB5uyfNDf04uUloDoOSaTsjDPHRf-arZ-HzZ-2324fvvUKWsDjFbpLJpW9uHuM-em7r20Xd9Y9tsthts34QuXZKz2rbJX_3mlHy-zFZPb_l88fr-9DjPSyaLPkfKKw4ORjUNlXa151A5D1LAKEyZs6rmJVJVSF8hZyidV1CidSWAEo5Nye3x7zaG3eBTbzZNKv0o1_kwJEOFoqgFHcGbP-A6DLEb3YymyKGQgo8QHqEyhpSir802Nhsb9wbBHNY0_9ZkP-oFZaA</addsrcrecordid><sourcetype>Aggregation Database</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype><pqid>921407654</pqid></control><display><type>article</type><title>Long-Time-Step Methods for Oscillatory Differential Equations</title><source>SIAM Journals Online</source><creator>García-Archilla, B. ; Sanz-Serna, J. M. ; Skeel, R. D.</creator><creatorcontrib>García-Archilla, B. ; Sanz-Serna, J. M. ; Skeel, R. D.</creatorcontrib><description>Considered are numerical integration schemes for nondissipative dynamical systems in which multiple time scales are present. It is assumed that one can do an explicit separation of the RHS "forces" into fast forces and slow forces such that (i) the fast forces contain the high frequency part of the solution, (ii) the fast forces are conservative, and (iii) the reduced problem consisting only of the fast forces can be integrated much more cheaply than the full problem. The fast forces are allowed to have low frequency components. Particular applications for which the schemes are intended include N-body problems (for which most of the forces are slow) and nonlinear wave phenomena (for which the fast forces can be propagated by spectral methods). The assumption of cheap integration of fast forces implies that the overall cost of integration is primarily determined by the step size used to sample the slow forces. A long-time-step method is one in which this step size exceeds half the period of the fastest normal mode present in the full system. An existing method that comes close to qualifying is the "impulse" method, also known as Verlet-I and r-RESPA. It is shown that it might fail, though, for a couple of reasons. First, it suffers a serious loss of accuracy if the step size is near a multiple of the period of a normal mode, and, second, it is unstable if the step size is near a multiple of half the period of a normal mode. Proposed in this paper is a "mollified" impulse method having an error bound that is independent of the frequency of the fast forces. It is also shown to possess superior stability properties. Theoretical results are supplemented by numerical experiments. The method is efficient and reasonably easy to implement.</description><identifier>ISSN: 1064-8275</identifier><identifier>EISSN: 1095-7197</identifier><identifier>DOI: 10.1137/S1064827596313851</identifier><language>eng</language><publisher>Philadelphia: Society for Industrial and Applied Mathematics</publisher><subject>Accuracy ; Approximation ; Eigenvalues ; Experiments ; Grants ; Methods</subject><ispartof>SIAM journal on scientific computing, 1998, Vol.20 (3), p.930-963</ispartof><rights>[Copyright] © 1998 Society for Industrial and Applied Mathematics</rights><lds50>peer_reviewed</lds50><woscitedreferencessubscribed>false</woscitedreferencessubscribed><citedby>FETCH-LOGICAL-c367t-124d40b063190d9bfe40dbe065075923ba8f4c12876ed14316be80c1abc0085b3</citedby><cites>FETCH-LOGICAL-c367t-124d40b063190d9bfe40dbe065075923ba8f4c12876ed14316be80c1abc0085b3</cites></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><link.rule.ids>314,780,784,3184,4024,27923,27924,27925</link.rule.ids></links><search><creatorcontrib>García-Archilla, B.</creatorcontrib><creatorcontrib>Sanz-Serna, J. M.</creatorcontrib><creatorcontrib>Skeel, R. D.</creatorcontrib><title>Long-Time-Step Methods for Oscillatory Differential Equations</title><title>SIAM journal on scientific computing</title><description>Considered are numerical integration schemes for nondissipative dynamical systems in which multiple time scales are present. It is assumed that one can do an explicit separation of the RHS "forces" into fast forces and slow forces such that (i) the fast forces contain the high frequency part of the solution, (ii) the fast forces are conservative, and (iii) the reduced problem consisting only of the fast forces can be integrated much more cheaply than the full problem. The fast forces are allowed to have low frequency components. Particular applications for which the schemes are intended include N-body problems (for which most of the forces are slow) and nonlinear wave phenomena (for which the fast forces can be propagated by spectral methods). The assumption of cheap integration of fast forces implies that the overall cost of integration is primarily determined by the step size used to sample the slow forces. A long-time-step method is one in which this step size exceeds half the period of the fastest normal mode present in the full system. An existing method that comes close to qualifying is the "impulse" method, also known as Verlet-I and r-RESPA. It is shown that it might fail, though, for a couple of reasons. First, it suffers a serious loss of accuracy if the step size is near a multiple of the period of a normal mode, and, second, it is unstable if the step size is near a multiple of half the period of a normal mode. Proposed in this paper is a "mollified" impulse method having an error bound that is independent of the frequency of the fast forces. It is also shown to possess superior stability properties. Theoretical results are supplemented by numerical experiments. The method is efficient and reasonably easy to implement.</description><subject>Accuracy</subject><subject>Approximation</subject><subject>Eigenvalues</subject><subject>Experiments</subject><subject>Grants</subject><subject>Methods</subject><issn>1064-8275</issn><issn>1095-7197</issn><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>1998</creationdate><recordtype>article</recordtype><sourceid>8G5</sourceid><sourceid>ABUWG</sourceid><sourceid>AFKRA</sourceid><sourceid>AZQEC</sourceid><sourceid>BENPR</sourceid><sourceid>CCPQU</sourceid><sourceid>DWQXO</sourceid><sourceid>GNUQQ</sourceid><sourceid>GUQSH</sourceid><sourceid>M2O</sourceid><recordid>eNplkDtPwzAUhS0EEqXwA9giBrbAvX7FHhgQlIdU1KFltuzEgVRp3NrJ0H9PqjLBdK50Pl19OoRcI9whsuJ-iSC5ooXQkiFTAk_IBEGLvEBdnB5uyfNDf04uUloDoOSaTsjDPHRf-arZ-HzZ-2324fvvUKWsDjFbpLJpW9uHuM-em7r20Xd9Y9tsthts34QuXZKz2rbJX_3mlHy-zFZPb_l88fr-9DjPSyaLPkfKKw4ORjUNlXa151A5D1LAKEyZs6rmJVJVSF8hZyidV1CidSWAEo5Nye3x7zaG3eBTbzZNKv0o1_kwJEOFoqgFHcGbP-A6DLEb3YymyKGQgo8QHqEyhpSir802Nhsb9wbBHNY0_9ZkP-oFZaA</recordid><startdate>1998</startdate><enddate>1998</enddate><creator>García-Archilla, B.</creator><creator>Sanz-Serna, J. M.</creator><creator>Skeel, R. D.</creator><general>Society for Industrial and Applied Mathematics</general><scope>AAYXX</scope><scope>CITATION</scope><scope>3V.</scope><scope>7WY</scope><scope>7WZ</scope><scope>7X2</scope><scope>7XB</scope><scope>87Z</scope><scope>88A</scope><scope>88F</scope><scope>88I</scope><scope>88K</scope><scope>8AL</scope><scope>8FE</scope><scope>8FG</scope><scope>8FH</scope><scope>8FK</scope><scope>8FL</scope><scope>8G5</scope><scope>ABJCF</scope><scope>ABUWG</scope><scope>AFKRA</scope><scope>ARAPS</scope><scope>ATCPS</scope><scope>AZQEC</scope><scope>BBNVY</scope><scope>BENPR</scope><scope>BEZIV</scope><scope>BGLVJ</scope><scope>BHPHI</scope><scope>CCPQU</scope><scope>D1I</scope><scope>DWQXO</scope><scope>FRNLG</scope><scope>F~G</scope><scope>GNUQQ</scope><scope>GUQSH</scope><scope>HCIFZ</scope><scope>JQ2</scope><scope>K60</scope><scope>K6~</scope><scope>K7-</scope><scope>KB.</scope><scope>L.-</scope><scope>L6V</scope><scope>LK8</scope><scope>M0C</scope><scope>M0K</scope><scope>M0N</scope><scope>M1Q</scope><scope>M2O</scope><scope>M2P</scope><scope>M2T</scope><scope>M7P</scope><scope>M7S</scope><scope>MBDVC</scope><scope>P5Z</scope><scope>P62</scope><scope>PATMY</scope><scope>PDBOC</scope><scope>PQBIZ</scope><scope>PQBZA</scope><scope>PQEST</scope><scope>PQQKQ</scope><scope>PQUKI</scope><scope>PTHSS</scope><scope>PYCSY</scope><scope>Q9U</scope><scope>7SC</scope><scope>8FD</scope><scope>L7M</scope><scope>L~C</scope><scope>L~D</scope></search><sort><creationdate>1998</creationdate><title>Long-Time-Step Methods for Oscillatory Differential Equations</title><author>García-Archilla, B. ; Sanz-Serna, J. M. ; Skeel, R. D.</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-c367t-124d40b063190d9bfe40dbe065075923ba8f4c12876ed14316be80c1abc0085b3</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>1998</creationdate><topic>Accuracy</topic><topic>Approximation</topic><topic>Eigenvalues</topic><topic>Experiments</topic><topic>Grants</topic><topic>Methods</topic><toplevel>peer_reviewed</toplevel><toplevel>online_resources</toplevel><creatorcontrib>García-Archilla, B.</creatorcontrib><creatorcontrib>Sanz-Serna, J. M.</creatorcontrib><creatorcontrib>Skeel, R. D.</creatorcontrib><collection>CrossRef</collection><collection>ProQuest Central (Corporate)</collection><collection>Access via ABI/INFORM (ProQuest)</collection><collection>ABI/INFORM Global (PDF only)</collection><collection>Agricultural Science Collection</collection><collection>ProQuest Central (purchase pre-March 2016)</collection><collection>ABI/INFORM Global (Alumni Edition)</collection><collection>Biology Database (Alumni Edition)</collection><collection>Military Database (Alumni Edition)</collection><collection>Science Database (Alumni Edition)</collection><collection>Telecommunications (Alumni Edition)</collection><collection>Computing Database (Alumni Edition)</collection><collection>ProQuest SciTech Collection</collection><collection>ProQuest Technology Collection</collection><collection>ProQuest Natural Science Collection</collection><collection>ProQuest Central (Alumni) (purchase pre-March 2016)</collection><collection>ABI/INFORM Collection (Alumni Edition)</collection><collection>Research Library (Alumni Edition)</collection><collection>Materials Science & Engineering Collection</collection><collection>ProQuest Central (Alumni Edition)</collection><collection>ProQuest Central UK/Ireland</collection><collection>Advanced Technologies & Aerospace Collection</collection><collection>Agricultural & Environmental Science Collection</collection><collection>ProQuest Central Essentials</collection><collection>Biological Science Collection</collection><collection>ProQuest Central</collection><collection>Business Premium Collection</collection><collection>Technology Collection</collection><collection>Natural Science Collection</collection><collection>ProQuest One Community College</collection><collection>ProQuest Materials Science Collection</collection><collection>ProQuest Central Korea</collection><collection>Business Premium Collection (Alumni)</collection><collection>ABI/INFORM Global (Corporate)</collection><collection>ProQuest Central Student</collection><collection>Research Library Prep</collection><collection>SciTech Premium Collection</collection><collection>ProQuest Computer Science Collection</collection><collection>ProQuest Business Collection (Alumni Edition)</collection><collection>ProQuest Business Collection</collection><collection>Computer Science Database</collection><collection>Materials Science Database</collection><collection>ABI/INFORM Professional Advanced</collection><collection>ProQuest Engineering Collection</collection><collection>ProQuest Biological Science Collection</collection><collection>ABI/INFORM Global</collection><collection>Agricultural Science Database</collection><collection>Computing Database</collection><collection>Military Database</collection><collection>Research Library</collection><collection>Science Database</collection><collection>Telecommunications Database</collection><collection>Biological Science Database</collection><collection>Engineering Database</collection><collection>Research Library (Corporate)</collection><collection>Advanced Technologies & Aerospace Database</collection><collection>ProQuest Advanced Technologies & Aerospace Collection</collection><collection>Environmental Science Database</collection><collection>Materials Science Collection</collection><collection>ProQuest One Business</collection><collection>ProQuest One Business (Alumni)</collection><collection>ProQuest One Academic Eastern Edition (DO NOT USE)</collection><collection>ProQuest One Academic</collection><collection>ProQuest One Academic UKI Edition</collection><collection>Engineering Collection</collection><collection>Environmental Science Collection</collection><collection>ProQuest Central Basic</collection><collection>Computer and Information Systems Abstracts</collection><collection>Technology Research Database</collection><collection>Advanced Technologies Database with Aerospace</collection><collection>Computer and Information Systems Abstracts Academic</collection><collection>Computer and Information Systems Abstracts Professional</collection><jtitle>SIAM journal on scientific computing</jtitle></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext</fulltext></delivery><addata><au>García-Archilla, B.</au><au>Sanz-Serna, J. M.</au><au>Skeel, R. D.</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>Long-Time-Step Methods for Oscillatory Differential Equations</atitle><jtitle>SIAM journal on scientific computing</jtitle><date>1998</date><risdate>1998</risdate><volume>20</volume><issue>3</issue><spage>930</spage><epage>963</epage><pages>930-963</pages><issn>1064-8275</issn><eissn>1095-7197</eissn><abstract>Considered are numerical integration schemes for nondissipative dynamical systems in which multiple time scales are present. It is assumed that one can do an explicit separation of the RHS "forces" into fast forces and slow forces such that (i) the fast forces contain the high frequency part of the solution, (ii) the fast forces are conservative, and (iii) the reduced problem consisting only of the fast forces can be integrated much more cheaply than the full problem. The fast forces are allowed to have low frequency components. Particular applications for which the schemes are intended include N-body problems (for which most of the forces are slow) and nonlinear wave phenomena (for which the fast forces can be propagated by spectral methods). The assumption of cheap integration of fast forces implies that the overall cost of integration is primarily determined by the step size used to sample the slow forces. A long-time-step method is one in which this step size exceeds half the period of the fastest normal mode present in the full system. An existing method that comes close to qualifying is the "impulse" method, also known as Verlet-I and r-RESPA. It is shown that it might fail, though, for a couple of reasons. First, it suffers a serious loss of accuracy if the step size is near a multiple of the period of a normal mode, and, second, it is unstable if the step size is near a multiple of half the period of a normal mode. Proposed in this paper is a "mollified" impulse method having an error bound that is independent of the frequency of the fast forces. It is also shown to possess superior stability properties. Theoretical results are supplemented by numerical experiments. The method is efficient and reasonably easy to implement.</abstract><cop>Philadelphia</cop><pub>Society for Industrial and Applied Mathematics</pub><doi>10.1137/S1064827596313851</doi><tpages>34</tpages></addata></record> |
fulltext | fulltext |
identifier | ISSN: 1064-8275 |
ispartof | SIAM journal on scientific computing, 1998, Vol.20 (3), p.930-963 |
issn | 1064-8275 1095-7197 |
language | eng |
recordid | cdi_proquest_miscellaneous_25821952 |
source | SIAM Journals Online |
subjects | Accuracy Approximation Eigenvalues Experiments Grants Methods |
title | Long-Time-Step Methods for Oscillatory Differential Equations |
url | https://sfx.bib-bvb.de/sfx_tum?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2024-12-25T21%3A19%3A01IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-proquest_cross&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=Long-Time-Step%20Methods%20for%20Oscillatory%20Differential%20Equations&rft.jtitle=SIAM%20journal%20on%20scientific%20computing&rft.au=Garc%C3%ADa-Archilla,%20B.&rft.date=1998&rft.volume=20&rft.issue=3&rft.spage=930&rft.epage=963&rft.pages=930-963&rft.issn=1064-8275&rft.eissn=1095-7197&rft_id=info:doi/10.1137/S1064827596313851&rft_dat=%3Cproquest_cross%3E2586040781%3C/proquest_cross%3E%3Curl%3E%3C/url%3E&disable_directlink=true&sfx.directlink=off&sfx.report_link=0&rft_id=info:oai/&rft_pqid=921407654&rft_id=info:pmid/&rfr_iscdi=true |