I have very recently picked up PyFR. I have been trying to compute the Taylor-Green vortex solutions reported here: On the utility of GPU accelerated high-order methods for unsteady flow simulations: A comparison with industry-standard tools - ScienceDirect
I used the .ini file provided with the MMC1 extras, via the Science Direct link, and ran the p5 case on the 43^3 mesh.
This is what I was expecting:
This is what I got:
This is the published ini file (to which I have added tend=169.0):
;cblas-mt = Enter path to local multi-threaded BLAS library for OpenMP backend [constants] gamma = 1.4 mu = 7.395099728874521e-05 Pr = 0.71 M=0.1 [solver] system = navier-stokes order = 5 [solver-time-integrator] scheme = rk45 controller = pi t0 = 0.0 dt = 0.0001 atol = 0.000001 rtol = 0.000001 safety-fact = 0.9 min-fact = 0.3 max-fact = 2.5 tend = 169.0 [solver-interfaces] riemann-solver = rusanov ldg-beta = 0.5 ldg-tau = 0.1 [solver-interfaces-quad] flux-pts = gauss-legendre [solver-elements-hex] soln-pts = gauss-legendre [soln-output] format = pyfrs-dir basedir = . basename = taylor_green-1-%(t).3f times = range(0, 169.03, 1001) [soln-ics] u = +0.118321595661992*sin(x)*cos(y)*cos(z) v = -0.118321595661992*cos(x)*sin(y)*cos(z) w = 0.0 p = 1.0+1.0*0.118321595661992*0.118321595661992/16*(cos(2*x)+cos(2*y))*(cos(2*z)+2) rho = (1.0+1.0*0.118321595661992*0.118321595661992/16*(cos(2*x)+cos(2*y))*(cos(2*z)+2))/1.0
I would be grateful for any advice.