Weak-Prescribed BCs

Hi, I would like to explore some weak-prescribed boundary conditions mentioned in “A Guide to the Implementation of Boundary Conditions in Compact High-Order Methods for Compressible Aerodynamics”, Mengaldo et al. Most of the Mako templates in PyFR for the advective fluxes at boundaries seem to be for the Weak-Riemann approach. I understand that unlike the Weak-Riemann approach, for the Weak-prescribed approach, there isn’t two set of states and there is no need for a Riemann Solver at the boundaries. Would it be possible to test the Weak-prescribed approach in PyFR, and if so, may I get some pointers on how I could go about doing that?

Edit: Perhaps it involves something like setting cflux_state = None in inters.py, and creating <%pyfr:macro name=‘bc_common_flux_state’..> similar to that in ghost.mako/slp-adia-wall.mako?

This form of enforcement is supported. In particular your BC needs to provide:

This is controlled by the BC class setting cflux_state = None which ensures that no default implementations are included here:

Then, everything is on your BC to provide a complete bc_common_flux_state along with the appropriate common solution states needed for LDG.

Regards, Freddie.

2 Likes