idaholab / moose

Multiphysics Object Oriented Simulation Environment
https://www.mooseframework.org
GNU Lesser General Public License v2.1
1.66k stars 1.03k forks source link

Adaptivity (coarsening) violates solute conservation #8349

Open dschwen opened 7 years ago

dschwen commented 7 years ago

Description of the enhancement or error report

If coarsening happens in a phase field simulation conserved order parameters do not have their value perfectly conserved, as the integral over the coarse element is likely to be different than the integral over the refined children. I would like to attempt to fix this by adding source terms that exactly compensate the local non-conservation.

Rationale for the enhancement or information for reproducing the error

Solute conservation is essential to produce reproducable results between runs with and without adaptivity. Solute drift affects the system's free energy considerably.

Identified impact

Eliminate the unexpected solute drift. Existing simulations are unaffected.

friedmud commented 7 years ago

There are other other options. Instead of doing interpolation we could use an L2 projection to get the new values. That would perfectly preserve the concentration. On Wed, Jan 18, 2017 at 11:16 PM Daniel Schwen notifications@github.com wrote:

Description of the enhancement or error report

If coarsening happens in a phase field simulation conserved order parameters do not have their value perfectly conserved, as the integral over the coarse element is likely to be different than the integral over the refined children. I would like to attempt to fix this by adding source terms that exactly compensate the local non-conservation. Rationale for the enhancement or information for reproducing the error

Solute conservation is essential to produce reproducable results between runs with and without adaptivity. Solute drift affects the system's free energy considerably. Identified impact

Eliminate the unexpected solute drift. Existing simulations are unaffected.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/idaholab/moose/issues/8349, or mute the thread https://github.com/notifications/unsubscribe-auth/AA1JMVzNnOeuHUFYCc-chQ6KNNsLvKnuks5rTuOKgaJpZM4Lnr7P .

dschwen commented 7 years ago

Yeah, I know. Let me run this with an actual problem. I have it done and working.

dschwen commented 7 years ago

Wait, there is no interpolation going on when coarsening, is there? I thought nodes are just dropped. Conversely I don't see where L2 projection would come into play.

YaqiWang commented 7 years ago

With DG, projection is local thus can be easily done. Just saying.

friedmud commented 7 years ago

Just dropping the nodes is interpolation for Lagrange...

On Thu, Jan 19, 2017 at 11:03 AM Yaqi notifications@github.com wrote:

With DG, projection is local thus can be easily done. Just saying.

— You are receiving this because you commented.

Reply to this email directly, view it on GitHub https://github.com/idaholab/moose/issues/8349#issuecomment-273816523, or mute the thread https://github.com/notifications/unsubscribe-auth/AA1JMbA8C2jy2BP_n7zQ4nUi1YzMbwTBks5rT4k2gaJpZM4Lnr7P .