Hi Eduardo,
I think determining first layer of elements from the wall of interest is
the best place to start for this. Native PyFR mesh file (.pyfrm) uses hdf5
format. Datasets starting with 'bcon_' stores the connectivity of a
boundary to the nearby cells. If you check a bcon_ dataset in a pyfrm file
you will see something like this;
DATASET "bcon_wall_p0" {
...
DATA {
...
(3): {
"quad",
185,
2,
0
},
...
It tells you that boundary face 3 shares the face 2 of element number 185
of the quad elements. Based on this you may want to store the element-local
indices of solution points that you are interested in. Once you construct
this you'll have a corresponding data set for the boundary you are working
on. solvers/euler/kernels/bccflux.mako is a good example of how boundary
conditions are implemented in PyFR.
quad element local solution point indices for p=2, face indices, and mesh
point indices
2________2________3
| |
| 6 7 8 |
| |
3| 3 4 5 |1
| |
| 0 1 2 |
|________________ |
0 0 1
and for hex elements this repeats itself in z direction
Best,
Semih
On Tuesday, February 25, 2020 at 1:18:21 AM UTC, Eduardo Molina wrote:
>
> Hi all,
>
> I would like to implement a basic Algebraic Wall Model in PyFR. The model
> consists of 2 simple equations for the turbulent boundary layer: one for
> the shear stress and another for the heat flux (in case of an isothermal
> BC). As a 1st attempt, we can restrict the Wall Model to probe the LES
> quantities within the first element off the wall. Another modification is
> that the we might need a slip condition that will let the velocity to be
> driven by the imposed shear stress.
>
> Regarding elements types at the boundary, I don't know how PyFR is
> structured but I was thinking that quad/hexa(boundary layer) is a good
> start.
>
> This will be my first interaction with PyFR and I think that the best way
> to learn a new code is by coding/adding things....
>
> Thus, I really appreciate if someone could help pointing me where I can
> start the implementation and where I can find LES test cases other than the
> ones already available in GitHub.
>
> Best,
> Eduardo
>
--
You received this message because you are subscribed to the Google Groups "PyFR
Mailing List" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web, visit
https://groups.google.com/d/msgid/pyfrmailinglist/29918e8a-52e4-403d-8ec1-69f8093669e9%40googlegroups.com.
