Hi Niki,
Excuse me for interrupting you.
I'm afraid that I had doubts about the solution about 2D cylinder with
Re=200. There were frequently errors like "NaNs detected at t = xxx" in
PyFR run when I increased the Reynolds number by only modify the "Uin"
(Re>140, and Re=140 is fine without errors).
The mailing list encouraged me about the section named "shock-capturing"
is not suitable for the AC method.
My question is :
a. What's wrong with the inc_cylinder_2d case when the Re > 140;
b. Are there easy way changing the Reynolds number without failure
of NaN check
I attached the diff file ($diff inc_cylinder_2d_my_edit.ini
inc_cylinder_2d.ini) and residual.csv below. By any chance, could you
help check the configuration file ?
Regards,
Lin
在 2018年04月20日 18:16, Niki Loppi 写道:
>
> Hi Will,
>
> Thanks for your message.
>
> I have only run external flows with a range of Reynolds numbers,
> including 2D cylinder Re=200, turbulent jet Re=10,000, SD7003
> Re=60,000. All these cases converged ok.
>
> In your previous email, you mentioned that you are using the
> forward-euler pseudo-time scheme and "the pseudo time step is 1 over
> 10 of the physical time step. 10 iteration per physical time step."
> You cannot expect any level of convergence with these numbers! I would
> encourage you to use RK4 / tvd-RK3 and make the iteration count 10-20
> times bigger. Using P-multigrid always helps.
>
> Would you be able to share your case or at least the .ini file? I can
> try to quickly run it or at least have a look.
>
> Regards,
> Niki
> --
> 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 pyfrmailinglist+unsubscr...@googlegroups.com
> <mailto:pyfrmailinglist+unsubscr...@googlegroups.com>.
> To post to this group, send email to pyfrmailinglist@googlegroups.com
> <mailto:pyfrmailinglist@googlegroups.com>.
> Visit this group at https://groups.google.com/group/pyfrmailinglist.
> For more options, visit https://groups.google.com/d/optout.
--
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 pyfrmailinglist+unsubscr...@googlegroups.com.
To post to this group, send an email to pyfrmailinglist@googlegroups.com.
Visit this group at https://groups.google.com/group/pyfrmailinglist.
For more options, visit https://groups.google.com/d/optout.
7d6
< cblas= /usr/lib/libcblas.so
10d8
< ;43x16 R=0.5*1 8 35 8 8
12c10
< Uin = 1.6
---
> Uin = 1.0
20,25d17
< ;shock-capturing = artificial-viscosity
<
< ;[solver-artificial-viscosity]
< ;max-artvisc = 0.01
< ;s0 = 0.01
< ;kappa = 5.0
61d52
< ;nsteps = 1
71c62
< basename = inc_cylinder_2d_Re{t:.2f}
---
> basename = cylinder_2d_{t:.2f}
n,t,i,p,u,v
1,0.0,1,-,-,-
2,0.0,2,-,-,-
3,0.0,3,2.9001386958319424,2.517386527261391,1.3455911043919124
4,0.0088,1,-,-,-
5,0.0088,2,-,-,-
6,0.0088,3,nan,nan,nan
7,0.0176,1,-,-,-
8,0.0176,2,-,-,-
9,0.0176,3,nan,nan,nan
10,0.0264,1,-,-,-
11,0.0264,2,-,-,-
12,0.0264,3,nan,nan,nan
13,0.0352,1,-,-,-
14,0.0352,2,-,-,-
15,0.0352,3,nan,nan,nan
16,0.044000000000000004,1,-,-,-
17,0.044000000000000004,2,-,-,-
18,0.044000000000000004,3,nan,nan,nan
19,0.05280000000000001,1,-,-,-
20,0.05280000000000001,2,-,-,-
21,0.05280000000000001,3,nan,nan,nan
22,0.06160000000000001,1,-,-,-
23,0.06160000000000001,2,-,-,-
24,0.06160000000000001,3,nan,nan,nan
25,0.0704,1,-,-,-
26,0.0704,2,-,-,-
27,0.0704,3,nan,nan,nan
28,0.0792,1,-,-,-
29,0.0792,2,-,-,-
30,0.0792,3,nan,nan,nan
31,0.08800000000000001,1,-,-,-
32,0.08800000000000001,2,-,-,-
33,0.08800000000000001,3,nan,nan,nan
34,0.09680000000000001,1,-,-,-
35,0.09680000000000001,2,-,-,-
36,0.09680000000000001,3,nan,nan,nan
37,0.10560000000000001,1,-,-,-
38,0.10560000000000001,2,-,-,-
39,0.10560000000000001,3,nan,nan,nan
40,0.11440000000000002,1,-,-,-
41,0.11440000000000002,2,-,-,-
42,0.11440000000000002,3,nan,nan,nan
43,0.12320000000000002,1,-,-,-
44,0.12320000000000002,2,-,-,-
45,0.12320000000000002,3,nan,nan,nan
46,0.132,1,-,-,-
47,0.132,2,-,-,-
48,0.132,3,nan,nan,nan
49,0.1408,1,-,-,-
50,0.1408,2,-,-,-
51,0.1408,3,nan,nan,nan
52,0.1496,1,-,-,-
53,0.1496,2,-,-,-
54,0.1496,3,nan,nan,nan
55,0.1584,1,-,-,-
56,0.1584,2,-,-,-
57,0.1584,3,nan,nan,nan
58,0.16720000000000002,1,-,-,-
59,0.16720000000000002,2,-,-,-
60,0.16720000000000002,3,nan,nan,nan
61,0.17600000000000002,1,-,-,-
62,0.17600000000000002,2,-,-,-
63,0.17600000000000002,3,nan,nan,nan
64,0.18480000000000002,1,-,-,-
65,0.18480000000000002,2,-,-,-
66,0.18480000000000002,3,nan,nan,nan
67,0.19360000000000002,1,-,-,-
68,0.19360000000000002,2,-,-,-
69,0.19360000000000002,3,nan,nan,nan
70,0.20240000000000002,1,-,-,-
71,0.20240000000000002,2,-,-,-
72,0.20240000000000002,3,nan,nan,nan
73,0.21120000000000003,1,-,-,-
74,0.21120000000000003,2,-,-,-
75,0.21120000000000003,3,nan,nan,nan
76,0.22000000000000003,1,-,-,-
77,0.22000000000000003,2,-,-,-
78,0.22000000000000003,3,nan,nan,nan
79,0.22880000000000003,1,-,-,-
80,0.22880000000000003,2,-,-,-
81,0.22880000000000003,3,nan,nan,nan
82,0.23760000000000003,1,-,-,-
83,0.23760000000000003,2,-,-,-
84,0.23760000000000003,3,nan,nan,nan
85,0.24640000000000004,1,-,-,-
86,0.24640000000000004,2,-,-,-
87,0.24640000000000004,3,nan,nan,nan
88,0.25520000000000004,1,-,-,-
89,0.25520000000000004,2,-,-,-
90,0.25520000000000004,3,nan,nan,nan
91,0.264,1,-,-,-
92,0.264,2,-,-,-
93,0.264,3,nan,nan,nan
94,0.2728,1,-,-,-
95,0.2728,2,-,-,-
96,0.2728,3,nan,nan,nan
97,0.28159999999999996,1,-,-,-
98,0.28159999999999996,2,-,-,-
99,0.28159999999999996,3,nan,nan,nan
100,0.29039999999999994,1,-,-,-
101,0.29039999999999994,2,-,-,-
102,0.29039999999999994,3,nan,nan,nan
103,0.2991999999999999,1,-,-,-
104,0.2991999999999999,2,-,-,-
105,0.2991999999999999,3,nan,nan,nan
106,0.3079999999999999,1,-,-,-
107,0.3079999999999999,2,-,-,-
108,0.3079999999999999,3,nan,nan,nan
109,0.31679999999999986,1,-,-,-
110,0.31679999999999986,2,-,-,-
111,0.31679999999999986,3,nan,nan,nan
112,0.32559999999999983,1,-,-,-
113,0.32559999999999983,2,-,-,-
114,0.32559999999999983,3,nan,nan,nan
115,0.3343999999999998,1,-,-,-
116,0.3343999999999998,2,-,-,-
117,0.3343999999999998,3,nan,nan,nan
118,0.3431999999999998,1,-,-,-
119,0.3431999999999998,2,-,-,-
120,0.3431999999999998,3,nan,nan,nan
121,0.35199999999999976,1,-,-,-
122,0.35199999999999976,2,-,-,-
123,0.35199999999999976,3,nan,nan,nan
124,0.36079999999999973,1,-,-,-
125,0.36079999999999973,2,-,-,-
126,0.36079999999999973,3,nan,nan,nan
127,0.3695999999999997,1,-,-,-
128,0.3695999999999997,2,-,-,-
129,0.3695999999999997,3,nan,nan,nan
130,0.3783999999999997,1,-,-,-
131,0.3783999999999997,2,-,-,-
132,0.3783999999999997,3,nan,nan,nan
133,0.38719999999999966,1,-,-,-
134,0.38719999999999966,2,-,-,-
135,0.38719999999999966,3,nan,nan,nan
136,0.39599999999999963,1,-,-,-
137,0.39599999999999963,2,-,-,-
138,0.39599999999999963,3,nan,nan,nan
139,0.4047999999999996,1,-,-,-
140,0.4047999999999996,2,-,-,-
141,0.4047999999999996,3,nan,nan,nan
142,0.4135999999999996,1,-,-,-
143,0.4135999999999996,2,-,-,-
144,0.4135999999999996,3,nan,nan,nan
145,0.42239999999999955,1,-,-,-
146,0.42239999999999955,2,-,-,-
147,0.42239999999999955,3,nan,nan,nan
