  Young planets embedded in their protoplanetary disk interact gravitationallywith it leading to energy and angular momentum exchange. This interactiondetermines the evolution of the planet through changes to the orbitalparameters. We investigate changes in the orbital elements of a 20 Earth--massplanet due to the torques from the disk. We focus on the non-linear evolutionof initially non-vanishing eccentricity $e$ and/or inclination $i$. We treatthe disk as a two- or three-dimensional viscous fluid and performhydrodynamical simulations with an embedded planet. We find rapid exponentialdecay of the planet orbital eccentricity and inclination for small initialvalues of $e$ and $i$, in agreement with linear theory. For larger values of $e> 0.1$ the decay time increases and the decay rate scales as $\dot{e} \proptoe^{-2}$, consistent with existing theoretical models. For large inclinations($i$ > 6 deg) the inclination decay rate shows an identical scaling $di/dt\propto i^{-2}$. We find an interesting dependence of the migration on theeccentricity. In a disk with aspect ratio $H/r=0.05$ the migration rate isenhanced for small non-zero eccentricities ($e < 0.1$), while for larger valueswe see a significant reduction by a factor of $\sim 4$. We find no indicationfor a reversal of the migration for large $e$, although the torque experiencedby the planet becomes positive when $e \simeq 0.3$. This inward migration iscaused by the persisting energy loss of the planet.  For non gap forming planets, eccentricity and inclination damping occurs on atime scale that is very much shorter than the migration time scale. The resultsof non linear hydrodynamic simulations are in very good agreement with lineartheory for small $e$ and $i$.
We treat the disk as a two- or three-dimensional viscous fluid and perform hydrodynamical simulations with an embedded planet.