Academy of Mathematics and Systems Science, CAS Colloquia & Seminars
Prof. Shengtai Li,Los Alamos National Lab., NM，USA
Numerical simulations for planet-disk interactions
We develop a simulation code for interaction between the proto-planetary disk and embedded proto-planets. The disk evolution is subject to viscous (or inviscid) hydrodynamics (Euler) equations in a 2-D polar grid with external gravity force from central star and embedded planets. The planets are advanced according to N-body Newtonian mechanics. This science problem poses exceptional challenges to the numerical simulations. It requires millions time steps. We proposed a hydro algorithm that is at least an order of magnitude faster than standard Godunov solver. This problem also requires a self-gravity solver that must be updated every time step. Our parallelized 2D self-gravity solver costs less than 10% of the total computational cost. This problem also requires enormous resolution around the planet to accurately represent the dust around planet. We have developed an embedded Lagrangian adaptive mesh refinement method that speeds up the computation by another order of magnitude.
Copyright@2008, All Rights Reserved, Academy of Mathematics and Systems Science, CAS
Tel: 86-10-62553063 Fax: 86-10-62541829 E-mail: email@example.com