Mathematical method for simulating evolution of solar system improved

Apr 11, 2013
From left to right, researchers from UPV/EHU Joseba Makazaga, Ander Murua and Mikel Antoñana, researchers at UPV/EHU Credit: UPV/EHU

In order to improve a simulation designed to study the evolution of the solar system through time, numerical mathematical methods have been developed at the Computing Faculty of the University of the Basque Country (UPV/EHU). Specifically, the methods proposed enable the simulation calculations to be done faster and more accurately.

The methodology developed at the UPV/EHU's Computing Faculty is a clear example of interdisciplinarity and collaboration. Indeed, mathematicians, computer scientists, physicists and astronomers have been working together on this task, and even though a large proportion of the work was done at the UPV/EHU, the Universities of Valencia and Castellon and the were also involved.

The mathematician Ander Murua has explained how the collaboration took place. "At the Paris Observatory there is a well-known astronomer by the name of Jacques Laskar who is doing research into the evolution of the solar system. Among other things, Laskar has developed precise mathematical models of the solar system and by using numerical methods in , he has made calculations to find out how the solar system has evolved over millions of years. When the astronomical information obtained by Laskar by means of these calculations and simulations is compared with geological data, it is possible to know how changes in the Earth's orbit are related to glaciations and warming, and this can help to predict what could happen in the future. Both the mathematical model of the and the numerical methods used are important for carrying out these simulations."

Laskar's team did their latest simulation about three years ago and went back 250 million years. It took the computers a whole year to do it. Nevertheless, according to Laskar, even though the results obtained for the last 50 million years are reliable, the ones further back lose reliability because of the chaotic behaviour of the system. Murua says, "Apparently, the next time Laskar does a simulation, he would like to achieve a reliability of the results for a period of 70 million years, by fine-tuning the and improving the numerical methods to do the calculations."

"That being the case," adds Murua, "Laskar issued a challenge during a sectoral conference that we held in Castellon". We had had no contact with Laskar, but through the physicists Fernando Casas and Sergio Blanes of the Universities of Castellon and Valencia, an opportunity presented itself to get in touch with the Catalan Ariadna Farres who works with Laskar." That is how the collaboration started.

Ten times faster

The Joseba Makazaga and Mikel Antoñana are working together with the mathematician Ander Murua at the UPV/EHU's Computing Faculty. "We responded to the challenge issued by Laskar and we have improved the numerical methods used in the simulation. Our team has largely been responsible for developing the numerical methods that are more effective than the ones known so far. Firstly, we have achieved greater accuracy, and secondly we have considerably reduced the time needed to do the calculations," stressed Murua.

In fact, they have carried out many experiments to verify the validity of the numerical methods they have developed, and they have seen that the simulation can be done ten times faster than with the previous methods. "We don't know when Laskar is planning to do the again, but when he does, we won't have to wait a year to find out the results, the work will be completed in a matter of weeks," he reckons.

In the meantime, an explanation of the developed has been published in the scientific journal Applied Numerical Mathematics in an article entitled New families of symplectic splitting methods for numerical integration in dynamical astronomy. In addition, Murua has announced that they will also be publishing the results of the comparison between the new methods and the previous ones in the journal Celestial Mechanics and Dynamical Astronomy.

Explore further: Magnetar discovered close to supernova remnant Kesteven 79

add to favorites email to friend print save as pdf

Related Stories

A glacier's life

Oct 29, 2008

EPFL researchers have developed a numerical model that can re-create the state of Switzerland's Rhône Glacier as it was in 1874 and predict its evolution until the year 2100. This is the longest period of time ever modeled ...

New method for solving differential equations

Jan 24, 2008

Dutch-sponsored mathematician Valeriu Savcenco has developed new methods for the numerical solution of ordinary differential equations. These so-called multirate methods are highly efficient for large systems, where some ...

A dynamical quantum simulator

Mar 16, 2012

(PhysOrg.com) -- An international collaboration demonstrates the superiority of a dynamical quantum simulator over state-of-the-art numerical calculations.

Recommended for you

How can we find tiny particles in exoplanet atmospheres?

Aug 29, 2014

It may seem like magic, but astronomers have worked out a scheme that will allow them to detect and measure particles ten times smaller than the width of a human hair, even at many light-years distance.  ...

Spitzer telescope witnesses asteroid smashup

Aug 28, 2014

(Phys.org) —NASA's Spitzer Space Telescope has spotted an eruption of dust around a young star, possibly the result of a smashup between large asteroids. This type of collision can eventually lead to the ...

User comments : 0