Symmetry BC with PyFR

I am trying to run the following test case from NASA: 3D ONERA M6 Wing Validation for Turbulence Model Numerical Analysis. How do I correctly setup the bc? Right now, my .ini file looks like:

 [soln-bcs-farfield]
type = char-riem-inv
rho = 0.648
u = 290.97
v = 15.57
w = 0.0
p = 101325

[soln-bcs-symmetry]
type = slp-adia-wall

[soln-bcs-wing3d]
type = no-slp-adia-wall

[soln-ics]
rho = 0.648
u = 290.97
v = 15.57
w = 0.0
p = 101325

because I did not know what to do about the symmetry plane. Is this the correct way to go about this? I cannot find any documentation on how to setup a symmetry plane with PyFR. Thank you!

Yes, the slp-adia-wall condition should be the one you want for a symmetry plane.

Regards, Freddie.

Thank you Freddie, I really appreciate it. This is my first time running a sim with a mesh that is not from one of PyFR’s test cases. I am running into some issues with the minimum time step rejected. I have read through others issues with this and I see that it is usually an issue with the mesh (this mesh is from NASA so I am leaving it be), bc, initial conditions, etc. I have adjusted the dt and the both the boundary and initial conditions. I pasted my .ini file and dtstats.csv. If you could point me in any direction I would greatly appreciate it. Thanks again.

.ini file:

[backend]
precision = single

[backend-cuda]
device-id = round-robin

[constants]
gamma = 1.4
mu = 0.0003333333333333333
Pr = 0.71

[solver]
system = navier-stokes
order = 3

[solver-time-integrator]
formulation = std
scheme = rk45
controller = pi
tstart = 0.0
tend = 10.0
dt = 0.0001
atol = 0.00001
rtol = 0.00001
errest-norm = l2
safety-fact = 0.9
min-fact = 0.3
max-fact = 2.5

[solver-interfaces]
riemann-solver = rusanov
ldg-beta = 0.5
ldg-tau = 0.1

[solver-interfaces-quad]
flux-pts = gauss-legendre

[solver-interfaces-tri]
flux-pts = williams-shunn

[solver-elements-hex]
soln-pts = gauss-legendre

[solver-elements-pri]
soln-pts = williams-shunn~gauss-legendre

[solver-interfaces-line]
flux-pts = gauss-legendre

[soln-plugin-nancheck]
nsteps = 50

[soln-plugin-writer]
dt-out = 0.001
basedir = .
basename = nasa-case-{t:.2f}

[soln-bcs-farfield]
type = char-riem-inv
rho = 0.648
u = 90.97
v = 15.57
w = 0.0
p = 101325

[soln-bcs-symmetry]
type = slp-adia-wall

[soln-bcs-wing3d]
type = no-slp-adia-wall

[soln-ics]
rho = 0.648
u = 90.97
v = 15.57
w = 0.0
p = 101325

[soln-plugin-dtstats]
flushsteps = 100
file = dtstats.csv
header = true

dtstats.csv:
n,t,dt,action,error
0,0.0,0.0001,reject,100
1,0.0,4.615752455922284e-05,reject,100
2,0.0,2.1305170734352593e-05,reject,100
3,0.0,9.833939414093153e-06,reject,100
4,0.0,4.539103000199141e-06,reject,100
5,0.0,2.095137582085339e-06,reject,100
6,0.0,9.670636440005477e-07,reject,100
7,0.0,4.4637263898286805e-07,reject,100
8,0.0,2.0603456046416841e-07,reject,100
9,0.0,9.510045284673536e-08,reject,100
10,0.0,4.3896014878664006e-08,reject,100
11,0.0,2.026131384813945e-08,reject,100
12,0.0,9.352120915476184e-09,reject,100
13,0.0,4.316707508369135e-09,reject,100
14,0.0,1.9924853283252993e-09,reject,100
15,0.0,9.196819047606618e-10,reject,28.938164234507518
16,0.0,5.081206668177727e-10,reject,15.84871021437981
17,0.0,3.0634461437824554e-10,reject,2.9362822702902642
18,0.0,2.358425042910731e-10,reject,1.1382910246964804
19,0.0,2.0830891508316063e-10,accept,0.7190350855388457
20,2.0830891508316063e-10,1.966625221175001e-10,accept,0.16431863286567536
21,4.0497143720066075e-10,2.2215100624074355e-10,accept,0.07799477528473359
22,6.271224434414043e-10,2.3943542762597077e-10,accept,0.028473276084589227
23,8.66557871067375e-10,2.761877231771012e-10,accept,0.019717862066721194
24,1.1427455942444763e-09,3.022791107137764e-10,accept,0.017840527667066865
25,1.4450247049582528e-09,3.2296600470771967e-10,accept,0.016948978160258376
26,1.7679907096659725e-09,3.440101987047923e-10,accept,0.01707195429604582
27,2.1120009083707647e-09,3.64076650835509e-10,accept,0.016961272362129732
28,2.4760775592062738e-09,3.8597004889844304e-10,accept,0.01679514946424965
29,2.862047608104717e-09,4.0948461354300185e-10,accept,0.016504030543987483
30,3.2715322216477187e-09,4.350847398337429e-10,accept,0.016251306172763207
31,3.7066169614814617e-09,4.624710087615291e-10,accept,0.016079099595960156
32,4.1690879702429905e-09,4.915439339092066e-10,accept,0.01613388128587755
33,4.660631904152197e-09,5.216031400839975e-10,accept,0.02443463131839532
34,5.182235044236195e-09,5.213560852897643e-10,accept,0.20261291055711889
35,5.703591129525959e-09,4.0054455539284236e-10,reject,1.2173134192919912
36,5.703591129525959e-09,2.9628525573637306e-10,accept,0.9099031078211285
37,5.9998763852623326e-09,2.2861161933766772e-10,accept,0.06943504850525511
38,6.228488004600001e-09,2.9992065935340944e-10,accept,0.027480441761855438
39,6.52840866395341e-09,3.4352331789545816e-10,accept,0.005970321830075839
40,6.871931981848868e-09,4.454270096036279e-10,accept,0.005485237203176594
41,7.317358991452496e-09,4.980987886692866e-10,accept,0.007196219341477931
42,7.815457780121782e-09,5.307541229503111e-10,accept,0.05403133148521194
43,8.346211903072093e-09,4.3440841052064004e-10,accept,0.2033311669924494
44,8.780620313592734e-09,3.6255922224054683e-10,accept,0.799820778324288
45,9.143179535833281e-09,2.8513395765744114e-10,accept,0.2297738012846753
46,9.428313493490723e-09,3.1025379995407556e-10,accept,0.04859607038926159
47,9.738567293444799e-09,3.709701837709702e-10,accept,0.00676790125487262
48,1.0109537477215769e-08,5.014859888458266e-10,accept,0.0059377095391961055
49,1.0611023466061596e-08,5.617256260281819e-10,accept,0.13913323363022348
50,1.1172749092089777e-08,3.928251260909575e-10,accept,0.5722901395907516
51,1.1565574218180734e-08,3.116360220995868e-10,accept,0.5382569616687836
52,1.1877210240280321e-08,2.893652458993946e-10,accept,0.039590762011300705
53,1.2166575486179716e-08,3.897549886307243e-10,accept,0.010244575614795911
54,1.2556330474810439e-08,4.855790280093976e-10,accept,0.006048926506055424
55,1.3041909502819837e-08,5.665796128733548e-10,accept,0.3630252741621902
56,1.3608489115693192e-08,3.4545288986313194e-10,reject,1.3939190809167303
57,1.3608489115693192e-08,2.663853049209589e-10,accept,0.8078983096635685
58,1.3874874420614151e-08,2.2232045466771408e-10,accept,0.0974995872237953
...

I believe that mesh is meant for steady-state RANS solvers which would make it very ill-suited for high-order LES/DNS codes and would likely result in the instabilities that you’re seeing.

Thanks for the response! Just curious, would make a mesh more suitable for RANS solvers vs high order codes like PyFR?

RANS meshes generally have very high aspect ratio/stretched cells in the boundary layer with high growth rates (which is usually fine for low-order steady state solvers) whereas ILES/DNS meshes (particularly for high-order explicit schemes) generally need lower aspect ratio cells in regions of high gradients and lower growth rates. This can be mitigated slightly with anti-aliasing but is primarily a stability/accuracy issue.