Hello everyone,
I recently read the book about high-order flux reconstruction schemes and found the vector correction function for triangle element. I would like to ask a question about how to get these vectors in PyFR? The references include the Ph.D. thesis of P. Castonguay and the chapter 10 of F.D. Witherden, P.E. Vincent and A. Jameson.
Regards
Kiny
Hi,
During an FR step one only needs the divergence of these functions. Moreover, it turns out to be simpler to directly compute the divergence of these functions. This is the approach taken by PyFR.
Regards, Freddie.
Hi Kiny,
As Freddie says, PyFR only deals with the divergence of these correction functions, and not the vector correction functions themselves.
The DG variants are the default coded into PyFR.
To try others you could adopt an approach similar to that which Brian adopted here for tensor product elements:
https://www.sciencedirect.com/science/article/pii/S0021999116304557
in which I believe he modified the code for calculating the m3 matrix to read from a file instead of calling gbasis_at (there should be a patch in the Supplementary Material).
One warning about this approach, is that the m3 coefficients you read from a file will have to have been generated in the context of a specific set of solution points, so if you run with anything but those solutions points it will lead to problems! I guess another approach would be to generalise gbasis_at.
Thanks
Peter
Thank you Freddie. I finally understand what you mean. It is more convenient to use the divergence of correction function as they are fixed form. Regards Kiny
Hi Peter,
Thank you very much. I found the ‘patch.patch’ in the supplementary material as follow. It replaces gbasis_at with the python operation of loadtxt for specific degree of polynomial. I will study the theory in this paper.
Regards
Kiny
From 44d8c84475255aed777c4fd728fcca2feaa35620 Mon Sep 17 00:00:00 2001
