# 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:** 1\
**Showing post:** 7

<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:** [14 June 2023 09:00 UTC](https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840/7 "2023-06-14T09:00:35Z")

</div>

```auto
[solver]
system = navier-stokes ;compressible n-s
order = 1 ;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-6

```

When applying entropy filter with reference to [[Gradients without writing to disk - #8 by tdzanic](https://pyfr.discourse.group/t/gradients-with-writing-to-disk/436/8). ], i got error like below.

```auto
(pyfr) sk@sk-comsik:~/python3.10/pyfr/airfoil/clf5605_DB/Re25000_M08/aoa10/p1_euler$ mpirun -np 4 pyfr run -b cuda -p grid1.pyfrm Re25000_M08_euler.ini
Traceback (most recent call last):
  File "/home/sk/python3.10/pyfr/bin/pyfr", line 8, in <module>
    sys.exit(main())
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/ __main__.py", line 118, in main
    args.process(args)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/ __main__.py", line 251, in process_run
    _process_common(
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/ __main__.py", line 236, in _process_common
    solver = get_solver(backend, rallocs, mesh, soln, cfg)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/solvers/ __init__.py", line 16, in get_solver
    return get_integrator(backend, systemcls, rallocs, mesh, initsoln, cfg)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/integrators/ __init__.py", line 36, in get_integrator
    return integrator(backend, systemcls, rallocs, mesh, initsoln, cfg)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/integrators/dual/phys/controllers.py", line 8, in __init__
    super(). __init__ (*args, **kwargs)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/integrators/dual/phys/steppers.py", line 16, in __init__
    super(). __init__ (*args, **kwargs)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/integrators/dual/phys/base.py", line 16, in __init__
    self.pseudointegrator = get_pseudo_integrator(
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/integrators/dual/pseudo/ __init__.py", line 76, in get_pseudo_integrator
    return pseudointegrator(backend, systemcls, rallocs, mesh,
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/integrators/dual/pseudo/pseudocontrollers.py", line 11, in __init__
    super(). __init__ (*args, **kwargs)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/integrators/dual/pseudo/base.py", line 51, in __init__
    self.system = systemcls(backend, rallocs, mesh, initsoln,
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/solvers/base/system.py", line 35, in __init__
    eles, elemap = self._load_eles(rallocs, mesh, initsoln, nregs, nonce)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/solvers/base/system.py", line 122, in _load_eles
    ele.set_backend(self.backend, nregs, nonce, linoff)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/solvers/navstokes/elements.py", line 34, in set_backend
    super().set_backend(*args, **kwargs)
  File "/home/sk/python3.10/pyfr/lib/python3.10/site-packages/pyfr/solvers/baseadvecdiff/elements.py", line 147, in set_backend
    raise ValueError('Invalid shock capturing scheme')
ValueError: Invalid shock capturing scheme
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

```

Do I need a version or action to use entropy filters? I’ve only used the artificial viscosity method.

---

_[View the full topic](https://pyfr.discourse.group/t/low-reynolds-transonic-high-aoa-condition/840)._