# Low Reynolds, Transonic, High AoA condition

**URL:** <https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840>\
**Category:** Cases\
**Created:** [14 June 2023 05:24 UTC](https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840 "2023-06-14T05:24:14Z")\
**Posts on this page:** 3\
**Page:** 2

<div class="post-metadata">

**Author:** ![YoonSuKyung](https://avatars.discourse-cdn.com/v4/letter/y/5e9695/32.png) [@YoonSuKyung](https://pyfr.discourse.group/u/YoonSuKyung)\
**Post date:** [19 June 2023 06:04 UTC](https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840/21 "2023-06-19T06:04:45Z")

</div>

This is the configuration file I tried.

```auto
[backend]
precision = single ;double
rank-allocator = linear

[backend-cuda]
device-id = local-rank ;round-robin
;mpi-type = cuda-aware ;standard

[constants]
gamma = 1.2941
Pr = 0.72

; aoa = 10 ; 
Uin = 185.416386 ; 
Vin = 32.693911 ; 

; aoa = 20 ; 
; Uin = 176.922256 ; 
; Vin = 64.394435 ; 
w = 0.0

Pc = 64.21955525 ;pa
rhoc= 0.001500451 ;kg/m^3

mu = 1.370e-5
cpTref =230401 ;j~ * K
cpTs = 187359.12

[solver]
system = navier-stokes ;compressible n-s
order = 4 ;order of polynomial solution basis
;anti-alias = flux
viscosity-correction = sutherland ;none
shock-capturing = entropy-filter ;artificial-viscosity ;none

[solver-entropy-filter]
d-min = 1e-6
p-min = 1e-6
e-tol = 1e-4

;[solver-artificial-viscosity]
;max-artvisc = 0.01
;s0 = 0.01
;kappa = 5.0

[solver-time-integrator]
formulation = std
scheme = rk45
controller = none
tstart = 0.0
tend = 1.001
dt = 1e-6

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

[solver-interfaces-line]
flux-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-interfaces-tri]
flux-pts = alpha-opt ;williams-shunn
quad-deg = 10
quad-pts = williams-shunn

[solver-interfaces-quad]
flux-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-elements-tri]
flux-pts = alpha-opt ;williams-shunn
soln-pts = williams-shunn
quad-deg = 10
quad-pts = williams-shunn

[solver-elements-quad]
flux-pts = gauss-legendre-lobatto
soln-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-elements-hex]
flux-pts = gauss-legendre-lobatto
soln-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-elements-tet]
flux-pts = alpha-opt ;shunn-ham
soln-pts = shunn-ham
quad-deg = 10
quad-pts = shunn-ham

[solver-elements-pri]
flux-pts = alpha-opt ;williams-shunn~gauss-legendre-lobatto
soln-pts = williams-shunn~gauss-legendre-lobatto
quad-deg = 10
quad-pts = williams-shunn~gauss-legendre-lobatto

[solver-elements-pyr]
flux-pts = gauss-legendre-lobatto
soln-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = witherden-vincent

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

[soln-plugin-tavg-1]
nsteps = 10
dt-out = 0.5
tstart = 2.6
mode = continuous
basedir = .
basename = tavg-{t:04.2f}

avg-p = p
avg-u = u
avg-v = v
;avg-w = w
avg-pp = p*p
avg-uu = u*u
avg-uv = u*v
;avg-uw = u*w
avg-vv = v*v
;avg-vw = v*w
;avg-ww = w*w

fun-avg-upup = uu - u*u
fun-avg-upvp = uv - u*v
;fun-avg-upwp = uw - u*w
fun-avg-vpvp = vv - v*v
;fun-avg-vpwp = vw - v*w
;fun-avg-wpwp = ww - w*w
fun-avg-urms = sqrt(uu - u*u + vv - v*v)
;fun-avg-urms = sqrt(uu - u*u + vv - v*v +ww - w*w)

[soln-plugin-nancheck]
nsteps = 10

[soln-plugin-residual]
nsteps = 10
file = residual.csv
header = true

;[soln-plugin-pseudostats]
;flushsteps = 10
;file = residual.csv
;header = true

[soln-plugin-fluidforce-upper]
nsteps = 10
file = upper_forces.csv
header = true
morigin = (0.0, 0.0, 0.0)

[soln-plugin-fluidforce-lower]
nsteps = 10
file = lower_forces.csv
header = true
morigin = (0.0, 0.0, 0.0)

[soln-bcs-inlet]
type = char-riem-inv
;type =sub-in-frv
rho = rhoc
u = Uin
v = Vin
w = 0
p = Pc

[soln-bcs-outlet]
type = char-riem-inv
rho = rhoc
u = Uin
v = Vin
w = 0
p = Pc

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

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

[soln-ics]
rho = rhoc
u = Uin
v = Vin
w = 0
p = Pc

```

---

<div class="post-metadata">

**Author:** ![WillT](https://avatars.discourse-cdn.com/v4/letter/w/bc79bd/32.png) [@WillT](https://pyfr.discourse.group/u/WillT)\
**Post date:** [19 June 2023 07:34 UTC](https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840/22 "2023-06-19T07:34:50Z")

</div>

> [@YoonSuKyung](#):
>
> ```auto
> [solver-elements-pri]
> flux-pts = alpha-opt ;williams-shunn~gauss-legendre-lobatto
> soln-pts = williams-shunn~gauss-legendre-lobatto
> 
> ```

For prism elements, you need to use:

```auto
soln-pts = alpha-opt~gauss-legendre-lobatto

```

---

<div class="post-metadata">

**Author:** ![YoonSuKyung](https://avatars.discourse-cdn.com/v4/letter/y/5e9695/32.png) [@YoonSuKyung](https://pyfr.discourse.group/u/YoonSuKyung)\
**Post date:** [19 June 2023 07:49 UTC](https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840/23 "2023-06-19T07:49:17Z")

</div>

I judged that I only need to change the flux point. Thanks to you, my calculation started now.  
My final setting is this.

```auto
[backend]
precision = single ;double
rank-allocator = linear

[backend-cuda]
device-id = local-rank ;round-robin
;mpi-type = cuda-aware ;standard

[constants]
gamma = 1.2941
Pr = 0.72

; aoa = 10 ; 
Uin = 185.416386 ; 
Vin = 32.693911 ; 

; aoa = 20 ; 
; Uin = 176.922256 ; 
; Vin = 64.394435 ; 
w = 0.0

Pc = 64.21955525 ;pa
rhoc= 0.001500451 ;kg/m^3

mu = 1.370e-5
cpTref =230401 ;j~ * K
cpTs = 187359.12

[solver]
system = navier-stokes ;compressible n-s
order = 4 ;order of polynomial solution basis
;anti-alias = flux
viscosity-correction = sutherland ;none
shock-capturing = entropy-filter ;artificial-viscosity ;none

[solver-entropy-filter]
d-min = 1e-6
p-min = 1e-6
e-tol = 1e-4

;[solver-artificial-viscosity]
;max-artvisc = 0.01
;s0 = 0.01
;kappa = 5.0

[solver-time-integrator]
formulation = std
scheme = euler ;rk4
controller = none
tstart = 0.0
tend = 1.001
dt = 5e-8

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

[solver-interfaces-line]
flux-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-interfaces-tri]
flux-pts = alpha-opt ;williams-shunn
quad-deg = 10
quad-pts = williams-shunn

[solver-interfaces-quad]
flux-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-elements-tri]
flux-pts = alpha-opt ;williams-shunn
soln-pts = alpha-opt ;williams-shunn
quad-deg = 10
quad-pts = williams-shunn

[solver-elements-quad]
flux-pts = gauss-legendre-lobatto
soln-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-elements-hex]
flux-pts = gauss-legendre-lobatto
soln-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = gauss-legendre-lobatto

[solver-elements-tet]
flux-pts = alpha-opt ;shunn-ham
soln-pts = alpha-opt ;shunn-ham
quad-deg = 10
quad-pts = shunn-ham

[solver-elements-pri]
flux-pts = alpha-opt ;williams-shunn~gauss-legendre-lobatto
soln-pts = alpha-opt~gauss-legendre-lobatto ;williams-shunn~gauss-legendre-lobatto
quad-deg = 10
quad-pts = williams-shunn~gauss-legendre-lobatto

[solver-elements-pyr]
flux-pts = gauss-legendre-lobatto
soln-pts = gauss-legendre-lobatto
quad-deg = 10
quad-pts = witherden-vincent

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

[soln-plugin-tavg-1]
nsteps = 10
dt-out = 0.5
tstart = 2.6
mode = continuous
basedir = .
basename = tavg-{t:04.2f}

avg-p = p
avg-u = u
avg-v = v
;avg-w = w
avg-pp = p*p
avg-uu = u*u
avg-uv = u*v
;avg-uw = u*w
avg-vv = v*v
;avg-vw = v*w
;avg-ww = w*w

fun-avg-upup = uu - u*u
fun-avg-upvp = uv - u*v
;fun-avg-upwp = uw - u*w
fun-avg-vpvp = vv - v*v
;fun-avg-vpwp = vw - v*w
;fun-avg-wpwp = ww - w*w
fun-avg-urms = sqrt(uu - u*u + vv - v*v)
;fun-avg-urms = sqrt(uu - u*u + vv - v*v +ww - w*w)

[soln-plugin-nancheck]
nsteps = 10

[soln-plugin-residual]
nsteps = 10
file = residual.csv
header = true

;[soln-plugin-pseudostats]
;flushsteps = 10
;file = residual.csv
;header = true

[soln-plugin-fluidforce-upper]
nsteps = 10
file = upper_forces.csv
header = true
morigin = (0.0, 0.0, 0.0)

[soln-plugin-fluidforce-lower]
nsteps = 10
file = lower_forces.csv
header = true
morigin = (0.0, 0.0, 0.0)

[soln-bcs-inlet]
type = char-riem-inv
;type =sub-in-frv
rho = rhoc
u = Uin
v = Vin
w = 0
p = Pc

[soln-bcs-outlet]
type = char-riem-inv
rho = rhoc
u = Uin
v = Vin
w = 0
p = Pc

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

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

[soln-ics]
rho = rhoc
u = Uin
v = Vin
w = 0
p = Pc

```

Let me check the result of applying the entry filter. Thanks a lot !!

[Previous page](https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840.md?page=1)