Recently, the class of Hamiltonian Boundary Value Methods (HBVMs) has been introduced with the aim of preserving the energy associated with polynomial Hamiltonian systems (and, more in general, with all suitably regular Hamiltonian systems). However, many interesting problems admit other invariants besides the Hamiltonian function. It would be therefore useful to have methods able to preserve any number of independent invariants. This goal is achieved by generalizing the line-integral approach which HBVMs rely on, thus obtaining a number of generalizations which we collectively name Line Integral Methods. In fact, it turns out that this approach is quite general, so that it can be applied to any numerical method whose discrete solution can be suitably associated with a polynomial, such as a collocation method, as well as to any conservative problem. In particular, a completely conservative variant of both HBVMs and Gauss collocation methods is presented. Numerical experiments confirm the effectiveness of the proposed methods.
Line Integral Methods which preserve all invariants of conservative problems / Luigi Brugnano; Felice Iavernaro. - In: JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS. - ISSN 0377-0427. - STAMPA. - 236:(2012), pp. 3905-3919. [10.1016/j.cam.2012.03.026]
Line Integral Methods which preserve all invariants of conservative problems
BRUGNANO, LUIGI;
2012
Abstract
Recently, the class of Hamiltonian Boundary Value Methods (HBVMs) has been introduced with the aim of preserving the energy associated with polynomial Hamiltonian systems (and, more in general, with all suitably regular Hamiltonian systems). However, many interesting problems admit other invariants besides the Hamiltonian function. It would be therefore useful to have methods able to preserve any number of independent invariants. This goal is achieved by generalizing the line-integral approach which HBVMs rely on, thus obtaining a number of generalizations which we collectively name Line Integral Methods. In fact, it turns out that this approach is quite general, so that it can be applied to any numerical method whose discrete solution can be suitably associated with a polynomial, such as a collocation method, as well as to any conservative problem. In particular, a completely conservative variant of both HBVMs and Gauss collocation methods is presented. Numerical experiments confirm the effectiveness of the proposed methods.File | Dimensione | Formato | |
---|---|---|---|
jcam 236 (2012) 3905–3919.pdf
Accesso chiuso
Tipologia:
Versione finale referata (Postprint, Accepted manuscript)
Licenza:
Tutti i diritti riservati
Dimensione
789.99 kB
Formato
Adobe PDF
|
789.99 kB | Adobe PDF | Richiedi una copia |
I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.