Hello all, I'm trying to model viscous flow due to self-gravity in the interior of a partially molten asteroid, which existed in the early solar system. I'm using an altered version of step-22, and everything works fine until I try to use an axisymmetric domain- specifically, a 2-D domain representing a solid of rotation around the y axis.
Somehow, the pressure solution shows anomalies near r=0 and the div(u) = 0 condition is not well-satisfied. I think I have altered the system setup properly, modifying the divergence and the integration to reflect cylindrical coordinates, but something is still not right! Have other people had experience in implementing the Stokes equations in axisymmetric coordinates? Is there an example somewhere you can point me to? Thanks for any and all help! Roger _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
