Hi,
I ran into a problem (non convergence in the Riemann solver)
using the new 'quadratic_cartesian' interpolation scheme in a
3D simulation including selfgravity.
Looking into the details after a refinement step I figured
that the data in the newly refined blocks develop stripe-like
structures. This leads, for instance, to overshooting in
velocities, resp. mach numbers.
Using the old interpolation scheme (second_order_old)
everything is fine and the variables are smoothly prolongated
onto the newly refined blocks.
Is this a known problem, or is there a problem with the
new (quadratic_cartesian) interpolation scheme?
If you need more information, please don't hesitate to
ask for (e.g. 2D slices of the simulation data with the
two different interpolations schemes)
Thank you,
Robi
====================================
Robi Banerjee
Department of Physics and Astronomy
ABB-320, McMaster University
Hamilton, ON L8S 4M1
CANADA
e-mail: banerjee@physics.mcmaster.ca
phone : (905) 525-9140 x 23189
fax : (905) 546-1252
This archive was generated by hypermail 2b30 : Wed Nov 05 2003 - 11:05:22 CST