Hi John,
I'm following your saga! I was planning on looking into this problem,
but I haven't had the time to devise a test problem that could better
isolate it. But I can give advise on trying to troubleshoot this
problem. I would use two tests that need to be modified to follow a
color field:
1. 1D shock tube (ProblemType 1)
2. 2D Kelvin-Helmholtz Instability (ProblemType 8)
You can use the CollapseTest or something else that uses a color or
metallicity field as an example. These two tests are quick to run, and
the shock (boundary in the KH test) are well-defined and are easily
inspected in a debugger. Also, you can setup static refine regions
around the boundary layer to test whether AMR affects your solution.
About your changes to the Riemann solver: I'm not completely certain
multiplying the color left/right states by "dls(i,j)/dlX(i)" is the
thing to do, and I'd have to think about a better solution later.
About the CFRF warnings: this would indicate that a field is negative
that's not supposed be negative. Your warnings are saying that density
and total energy are going to be negative with the flux correction. Did
you see these warnings before you changes to the Riemann solver?
Good luck with this problem, and let us know if you have any other
questions.
Cheers,
John
> --
> You received this message because you are subscribed to the Google
> Groups "enzo-dev" group.
> To unsubscribe from this group and stop receiving emails from it, send
> an email to
enzo-dev+u...@googlegroups.com.
> To post to this group, send email to
enzo...@googlegroups.com.
> Visit this group at
http://groups.google.com/group/enzo-dev.
> For more options, visit
https://groups.google.com/groups/opt_out.
--
John Wise
Assistant Professor of Physics
Center for Relativistic Astrophysics, Georgia Tech
http://cosmo.gatech.edu