Hello PyFR Community, I'm currently trying to recreate the simulations from
the paper of Park, Vincent and Witherden (2017) on a NACA0021 Aerofoil in
deep stall.
I am using the provided mesh and configuration file for the p4.ini file and
p4_mesh and started with a warmup phase until t=100 with p=1 and Re=27000.
After the warmup phase I changed the parameters to p=4 and Re=270000.
However, I am getting a 'minimum sized time step rejected' error at around
t=226.
Could you help me and share any information and tips, how I can solve this
error?
I attached the configuration file after the warmup phase.
Best regards, Vu
--
You received this message because you are subscribed to the Google Groups "PyFR
Mailing List" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web, visit
https://groups.google.com/d/msgid/pyfrmailinglist/45754f4d-85e2-4715-bacc-c9c008942cd5n%40googlegroups.com.
{\rtf1\ansi\ansicpg1252\cocoartf1671\cocoasubrtf600
{\fonttbl\f0\fnil\fcharset0 Menlo-Regular;}
{\colortbl;\red255\green255\blue255;\red0\green0\blue0;}
{\*\expandedcolortbl;;\csgray\c0;}
\paperw11900\paperh16840\margl1440\margr1440\vieww10800\viewh8400\viewkind0
\pard\tx560\tx1120\tx1680\tx2240\tx2800\tx3360\tx3920\tx4480\tx5040\tx5600\tx6160\tx6720\pardirnatural\partightenfactor0
\f0\fs22 \cf2 \CocoaLigature0 precision = double\
rank-allocator = linear\
\
[backend-openmp]\
;cblas-st = \
;cblas-mt =\
\
[backend-cuda]\
device-id = local-rank\
mpi-type = cuda-aware\
\
[constants]\
gamma = 1.4\
mu = 0.0000037\
Pr = 0.72\
\
M = 0.1\
uc = 1.0\
rhoc = 1.0\
\
[solver]\
system = navier-stokes\
order = 4\
anti-alias = flux, surf-flux\
\
[solver-time-integrator]\
scheme = rk45\
controller = pi\
tstart = 0.0\
tend = 250.0\
dt = 2.5e-6\
atol = 1e-5\
rtol = 1e-5\
\
[solver-interfaces]\
riemann-solver = roem\
ldg-beta = 0.5\
ldg-tau = 0.1\
\
[solver-interfaces-quad]\
flux-pts = gauss-legendre\
quad-deg = 11\
quad-pts = gauss-legendre\
\
[solver-interfaces-line]\
flux-pts = gauss-legendre\
quad-deg = 11\
\
[solver-elements-tri]\
soln-pts = williams-shunn\
quad-deg = 9\
\
[solver-elements-quad]\
soln-pts = gauss-legendre\
quad-deg = 11\
quad-pts = gauss-legendre\
\
[solver-elements-hex]\
soln-pts = gauss-legendre\
quad-deg = 11\
quad-pts = gauss-legendre\
\
\
[soln-plugin-writer]\
dt-out = 5.0\
basedir = .\
basename = naca_p4_1c_\{t:.1f\}\
\
\
[soln-bcs-wall]\
type = no-slp-adia-wall\
\
[soln-bcs-inflow]\
type = char-riem-inv\
rho = rhoc\
u = uc\
v = 0\
w = 0\
p = rhoc*uc*uc/(M*M*gamma)\
\
[soln-bcs-outflow]\
type = char-riem-inv\
p = rhoc*uc*uc/(M*M*gamma)\
u = uc\
v = 0\
w = 0\
rho = rhoc\
\
[soln-ics]\
rho = rhoc\
u = uc\
v = 0\
w = 0\
p = rhoc*uc*uc/(M*M*gamma)\
\
[soln-filter]\
freq = 0\
cutoff = 0\
order = 16\
alpha = 36\
\
[soln-plugin-fluidforce-wall]\
nsteps = 500\
file = wall-forces.csv\
header = true \
\
[soln-plugin-dtstats]\
flushsteps = 500\
file = dtstats.csv\
header = true\
\
[soln-plugin-residual]\
nsteps = 500\
file = residual.csv\
header = true\
\
\
[soln-plugin-tavg]\
nsteps = 500\
dt-out = 5.0\
basedir = .\
basename = tavg_\{t:.1f\}\
\
avg-p = p\
avg-p2 = p*p\
avg-vel = sqrt(u*u + v*v)}
