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

TimeSeriesCallback for dimensions other than 2D #1869

Open
andrewwinters5000 opened this issue Mar 9, 2024 · 1 comment
Open

TimeSeriesCallback for dimensions other than 2D #1869

andrewwinters5000 opened this issue Mar 9, 2024 · 1 comment
Labels

Comments

@andrewwinters5000
Copy link
Member

andrewwinters5000 commented Mar 9, 2024

Further extensions of the functionality from #1855 are possible. Adding the ability of having gauge points on, e.g., a 1D TreeMesh should be straightforward. For instance, one does not need a Newton iteration as the affine mapping between computational coordinates and physical coordinates is easily inverted exactly.

Instead of using the barycenters, it could be possible to use the available normal directions instead. See this discussion where this idea is partially implemented for the UnstructuredMesh2D.

For 3D meshes, like P4estMesh, a bit more is involved but it should still be straightforward. The search to locate which element a given (x,y,z) points lives is likely easily extended from the strategy in #1855. However, it might become more performance critical, especially if AMR is used. This is because the search would need to be redone after each time the AMRCallback is applied. The Newton iteration to compute the reference coordinates could reuse the existing straight_side_quad_map, straight_side_quad_map_metrics, transfinite_quad_map, and transfinite_quad_map_metrics functions that are working in 3D. However, an efficient Newton iteration may need implemented that requires inverting a 3x3 matrix.

@patrickersing
Copy link
Contributor

I have opened PR #1873 to extend the TimeSeriesCallback to other dimensions specifically for TreeMesh

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

No branches or pull requests

2 participants