[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. ], i got error like below.
(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.