Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Make meshblock-ic's compatible with AMR #557

Open
adityakpandare opened this issue Aug 1, 2022 · 1 comment
Open

Make meshblock-ic's compatible with AMR #557

adityakpandare opened this issue Aug 1, 2022 · 1 comment
Assignees
Labels
amr Adaptive Mesh Refinement inciter

Comments

@adityakpandare
Copy link
Member

adityakpandare commented Aug 1, 2022

Currently, the element-meshblock data that assigns element-ids to mesh blocks does not get updated after AMR step. To have meshblocks compatible with AMR, the following needs to be done:

  • Use Refiner::coarseBlkElems to update elemblockid in Refiner after AMR. This will use element-ancestry from AMR::tetstore.
  • ALECG needs to update m_nodeblockid after AMR. This will require a setup similar to ALECG::lhs() (which updates normals via norm()) which is called every time AMR happens.
@adityakpandare adityakpandare added inciter amr Adaptive Mesh Refinement labels Aug 1, 2022
@adityakpandare adityakpandare self-assigned this Aug 1, 2022
@adityakpandare
Copy link
Member Author

adityakpandare commented Aug 9, 2022

This will also involve using the refined elemblockid to re-initialize solutions in Refiner::writeMesh() and Refiner::nodeinit(). Currently these functions are not using the updated elemblockid for post-AMR initialization.

  • nodeinit() is particularly important since errors are computed based on the solutions from this function; and these errors are used to drive initial refinement t0ref. However, if initial refinement based on mesh-block ICs is not required, then this can be dropped; since it will not affect solution.
  • writeMesh() is not too important since it is used only for writing t0ref results to file, and does not affect AMR or solution.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
amr Adaptive Mesh Refinement inciter
Projects
None yet
Development

No branches or pull requests

1 participant