Numerical long-term integration of geodesic equations of motion

More documents

Recommendations

Info

Numerical experimentsρ(0) = 1.7: still regular motion1.51.2510.750.50.25z0-0.25-0.5-0.75-1-1.250 0.25 0.5 0.75 1 1.25 1.5 1.75 2ρ0.250.150.05ρ .-0.05-0.151.5 1.6ρ1.7 1.8Figure: TrajectoryFigure: Sections at z = 018/26 | Jonathan Seyrich (Numerical long-term integration) Tübingen | January 21, 2013
Numerical experimentsρ(0) = 1.7:still regular motion∆H(τ)1e-0071e-0081e-0091e-0101e-0111e-0121e-013CM2, ε=0.1, s=4CM1, ε=0.1, s=4RK5con, h=1e-4RK5var, h=1e-4IGEM, ε=0.1, s=40 100000 200000 300000 400000 500000τFigure: Energy error ∆H = | [H(yn)−H 0]/H 0 |.Integrator T calc [s]IGEM, ɛ = 0.1, s = 4 2202.1CK5con, h = 10 −4 23813.2CK5var, h max = 10 −4 19207.7Astro1, ɛ = 0.1, s = 4 266.9Astro2, ɛ = 0.1, s = 4 2100.9Table: CPU calculation times T calc19/26 | Jonathan Seyrich (Numerical long-term integration) Tübingen | January 21, 2013
Page 1 and 2: Numerical long-term integration of
Page 3 and 4: Problem and motivationGeodesic equa
Page 5 and 6: Numerical properties of the equatio
Page 11 and 12: Structure preserving integrator for
Page 13 and 14: Numerical experimentsTestcase: Mank
Page 15 and 16: Numerical experiments◮ We choose:
Page 17: Numerical experimentsρ(0) = 30.7:
Page 21 and 22: Numerical experimentsρ(0) = 0.7: c
Page 23 and 24: Numerical experiments◮ IGEM calcu
Page 25 and 26: Summary◮ Energy conservation impo
Page 27: SummaryReferences II26/26 | Jonatha

Numerical long-term integration of geodesic equations of motion

Create successful ePaper yourself

Delete template?

Save as template?