Non-Boussinesq Convection
Non-Boussinesq Convection
I want to simulate non-Boussinesq convection and I am having a lot of difficultly getting convergence. I went back to tutorial #8 Transient flow and heat equations - Rayleigh-Benard instability, setup the model, and ran it with the Boussinesq approximation enabled. After this went well, I went through the model and started making the changes I thought I needed in order to do this without the Boissinesq approximation. I want to do this so I can use this approach on another model with large temperature differences where the Boussinesq approximation is invalid.
Some of the changes made (among others):
-Changed to transient model
-Added in thermal compressibility model and reference temperature
-Added in the buoyancy force: body force, force 2 = -9.81
-Tripled the the density of the mesh
From what I understand, I should be able to get similar results to the model using the Boussinesq approximation. The simulation runs up to a point where I can start to see the convection currents form and then the run diverges. I'm going to try making my time step sizes smaller and see if that does it. However, I am hoping somebody who has done this before might be able to take a look at what I am doing and see if anything stands out as an error.
Project attached.
Cheers,
Clayton
Some of the changes made (among others):
-Changed to transient model
-Added in thermal compressibility model and reference temperature
-Added in the buoyancy force: body force, force 2 = -9.81
-Tripled the the density of the mesh
From what I understand, I should be able to get similar results to the model using the Boussinesq approximation. The simulation runs up to a point where I can start to see the convection currents form and then the run diverges. I'm going to try making my time step sizes smaller and see if that does it. However, I am hoping somebody who has done this before might be able to take a look at what I am doing and see if anything stands out as an error.
Project attached.
Cheers,
Clayton
- Attachments
-
- non-boussinesq convection.rar
- (91.36 KiB) Downloaded 329 times
Last edited by MrChips on 09 Dec 2016, 18:33, edited 1 time in total.
Re: Non-Boussinesq Convection
With the smaller time step size (400 at 0.5) I was able to get the simulation to run a lot longer, but it eventually diverged coming to NaN. The image below is a comparison between what I got with the the full Navier Stokes model and the one from the tutorial that uses the Boussinesq approximation. If anybody has any ideas on how to make this simulation more stable it would be appreciated. I've listed out the new sif file below.
Code: Select all
Header
CHECK KEYWORDS Warn
Mesh DB "." "."
Include Path ""
Results Directory ""
End
Simulation
Max Output Level = 5
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
Simulation Type = Transient
Steady State Max Iterations = 20
Output Intervals = 1
Timestepping Method = BDF
BDF Order = 2
Timestep intervals = 400
Timestep Sizes = 0.5
Solver Input File = case.sif
Post File = t.ep
End
Constants
Gravity(4) = 0 -1 0 9.82
Stefan Boltzmann = 5.67e-08
Permittivity of Vacuum = 8.8542e-12
Boltzmann Constant = 1.3807e-23
Unit Charge = 1.602e-19
End
Body 1
Target Bodies(1) = 1
Name = "Body Property 1"
Equation = 1
Material = 1
Body Force = 1
Initial condition = 1
End
Solver 2
Equation = Heat Equation
Procedure = "HeatSolve" "HeatSolver"
Variable = Temperature
Exec Solver = Always
Stabilize = False
Bubbles = True
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 1
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-10
BiCGstabl polynomial degree = 2
Linear System Preconditioning = Diagonal
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
End
Solver 1
Equation = Navier-Stokes
Procedure = "FlowSolve" "FlowSolver"
Variable = Flow Solution[Velocity:2 Pressure:1]
Exec Solver = Always
Stabilize = False
Bubbles = True
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 1
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-7
BiCGstabl polynomial degree = 2
Linear System Preconditioning = Diagonal
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
End
Equation 1
Name = "Equation 1"
Convection = Computed
Active Solvers(2) = 2 1
End
Material 1
Name = "Water (room temperature)"
Reference Temperature = 293
Viscosity = 1.002e-3
Heat expansion Coefficient = 0.207e-3
Compressibility Model = Thermal
Heat Conductivity = 0.58
Relative Permittivity = 80.1
Sound speed = 1497.0
Heat Capacity = 4183.0
Density = 998.3
End
Body Force 1
Name = "BodyForce 1"
Flow Bodyforce 2 = -9.81
End
Initial Condition 1
Name = "InitialCondition 1"
Velocity 2 = 1e-9
Velocity 1 = 0
Temperature = 293
End
Boundary Condition 1
Target Boundaries(1) = 1
Name = "Top"
Velocity 1 = 0
Velocity 2 = 0
Temperature = 293
End
Boundary Condition 2
Target Boundaries(1) = 3
Name = "Bottom"
Velocity 1 = 0
Velocity 2 = 0
Temperature = 293.5
End
Boundary Condition 3
Target Boundaries(1) = 2
Name = "Sides"
Velocity 1 = 0
Velocity 2 = 0
Last edited by MrChips on 09 Dec 2016, 21:14, edited 1 time in total.
Re: Non-Boussinesq Convection
Can you try to set for Navier-Stokes solver:
Stabilize = True
Stabilize = True
-
- Site Admin
- Posts: 4828
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Non-Boussinesq Convection
Hi
You could try to set in both solver relaxation factor to ~0.7 and perhaps try with for both solvers:
There is a simple test case "CompressibleIdealGas" using these.
Are the two cases the same for small temperature differences? For ideal gas the linear heat expansion coefficient is 1/T.
-Peter
You could try to set in both solver relaxation factor to ~0.7 and perhaps try with for both solvers:
Code: Select all
Nonlinear Timestepping = Logical True
Are the two cases the same for small temperature differences? For ideal gas the linear heat expansion coefficient is 1/T.
-Peter
Re: Non-Boussinesq Convection
Thanks Peter, I will try this out. Where would I be able to download all the test cases? Is there a central repository for these?
Vencel, the simulation ended up diverging at exactly the same time, but thank you for the suggestion.
Vencel, the simulation ended up diverging at exactly the same time, but thank you for the suggestion.
-
- Site Admin
- Posts: 4828
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Non-Boussinesq Convection
Hi
Is it the linear or nonlinear system that diverges. You are using pretty crappy preconditioner (diagonal). More expensive preconditioner should give your better convergence on linear system level.
-Peter
Is it the linear or nonlinear system that diverges. You are using pretty crappy preconditioner (diagonal). More expensive preconditioner should give your better convergence on linear system level.
-Peter
Re: Non-Boussinesq Convection
Hi Clayton,MrChips wrote:Where would I be able to download all the test cases? Is there a central repository for these?
The test example for Compressible Ideal Gas can be found in the Github repository of Elmer as:
https://github.com/ElmerCSC/elmerfem/tr ... ssIdealGas
Yours Sincerely,
Anil Kunwar
Last edited by annier on 13 Dec 2016, 06:33, edited 1 time in total.
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
Re: Non-Boussinesq Convection
Embarrassingly I'm not entirely sure how to tell, I'm new to this and still have a lot to learn. I did change the preconditioner to ILU1 for both solvers, but the run doesn't get passed the first time step in this case. Here are the last number of lines from the log file (with diagonal preconditioning):raback wrote:Hi
Is it the linear or nonlinear system that diverges. You are using pretty crappy preconditioner (diagonal). More expensive preconditioner should give your better convergence on linear system level.
-Peter
Code: Select all
ERROR:: IterSolve: Failed convergence tolerances.
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.47989792E+74 2.0000000 ) :: heat equation
HeatSolve: iter: 1 Assembly: (s) 0.45 0.45
HeatSolve: iter: 1 Solve: (s) 0.25 0.25
HeatSolve: Result Norm : 4.7989792055454271E+073
HeatSolve: Relative Change : 2.0000000000000000
ComputeChange: SS (ITER=3) (NRM,RELC): ( 0.47989792E+74 2.0000000 ) :: heat equation
SolveEquations: -------------------------------------
SolveEquations: Coupled system iteration: 4
SolveEquations: -------------------------------------
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: navier-stokes
FlowSolve:
FlowSolve:
FlowSolve: -------------------------------------
FlowSolve: NAVIER-STOKES ITERATION 1
FlowSolve: -------------------------------------
FlowSolve:
FlowSolve: Starting Assembly...
FlowSolve: Assembly:
FlowSolve: Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
FlowSolve: Dirichlet conditions done
1 0.3869E+03
2 0.3229E+03
3 0.4516E+03
4 0.8149E+03
5 0.4841E+02
6 0.3914E+02
7 0.3395E+02
8 0.2861E+02
9 0.1054E+03
10 0.2101E+03
11 0.1394E+03
12 0.5918E+02
13 0.1001E+02
14 0.9994E+01
15 0.2639E+02
16 0.3533E+03
17 0.6608E+02
18 0.3472E+04
19 0.1027E+05
20 0.5832E+04
21 0.5818E+04
22 0.2839E+04
23 0.2117E+04
24 0.2128E+04
25 0.8592E+03
26 0.6754E+03
27 0.6626E+03
28 0.6548E+03
29 0.6548E+03
30 0.6439E+03
31 0.5953E+03
32 0.2797E+04
33 0.3843E+04
34 0.5409E+04
35 0.1557E+05
36 0.1745E+05
37 0.1721E+05
38 0.1718E+05
39 0.1913E+05
40 0.2254E+05
41 0.2911E+05
42 0.2783E+05
43 0.2987E+05
44 0.1646E+06
45 0.1092E+07
46 0.3392E+06
47 0.2016E+06
48 0.1235E+07
49 0.1107E+08
50 0.1128E+08
51 0.2505E+09
52 0.2329E+09
53 0.3711E+09
54 0.9509E+11
55 0.7596E+11
56 0.1282E+12
57 0.7190E+11
58 0.4221E+11
59 0.4593E+11
60 0.4766E+11
61 0.6924E+11
62 0.7112E+11
63 0.6596E+11
64 0.7494E+11
65 0.7465E+11
66 0.7395E+11
67 0.7300E+11
68 0.8723E+11
69 0.7451E+11
70 0.7355E+11
71 0.2872E+12
72 0.1569E+12
73 0.1700E+12
74 0.2032E+12
75 0.1790E+12
76 0.1535E+12
77 0.1486E+12
78 0.1740E+12
79 0.1403E+12
80 0.4196E+12
81 0.2678E+14
82 0.1135E+13
83 0.1507E+14
84 0.4672E+15
85 0.3749E+14
86 0.3606E+14
87 0.1201E+16
88 0.5837E+15
89 0.4647E+15
90 0.5023E+17
91 0.3005E+16
92 0.5697E+16
93 0.3441E+17
94 0.1527E+17
95 0.2788E+17
96 0.4114E+17
97 0.4116E+17
98 0.4142E+17
99 0.4084E+17
100 0.4132E+17
101 0.5692E+17
102 0.3996E+17
103 0.3902E+17
104 0.2449E+18
105 0.2353E+18
106 0.2361E+18
107 0.2338E+18
108 0.2271E+18
109 0.2416E+18
110 0.3025E+18
111 0.1824E+18
112 0.5750E+18
113 0.5669E+18
114 0.5433E+18
115 0.5483E+18
116 0.5499E+18
117 0.5779E+18
118 0.4954E+18
119 0.1294E+19
120 0.3608E+19
121 0.1031E+20
122 0.6594E+19
123 0.7181E+19
124 0.7931E+20
125 0.2989E+20
126 0.5778E+20
127 0.6063E+20
128 0.7922E+20
129 0.8373E+20
130 0.8367E+20
131 0.8364E+20
132 0.1967E+21
133 0.1430E+21
134 0.1365E+21
135 0.1481E+22
136 0.1237E+22
137 0.1964E+22
138 0.7245E+22
139 0.1450E+24
140 0.1777E+23
141 0.1742E+23
142 0.6035E+24
143 0.2437E+23
144 0.3662E+23
145 0.1297E+24
146 0.4331E+23
147 0.1688E+24
148 0.1695E+24
149 0.1797E+24
150 0.1972E+24
151 0.1914E+24
152 0.1773E+24
153 0.1810E+24
154 0.1787E+24
155 0.1805E+24
156 0.1853E+24
157 0.1808E+24
158 0.1664E+24
159 0.1654E+24
160 0.1696E+24
161 0.1743E+24
162 0.1689E+24
163 0.1852E+24
164 0.1851E+24
165 0.1856E+24
166 0.1868E+24
167 0.1889E+24
168 0.1889E+24
169 0.1947E+24
170 0.1980E+24
171 0.2421E+24
172 0.2274E+24
173 0.2130E+24
174 0.2477E+24
175 0.2615E+24
176 0.2603E+24
177 0.2598E+24
178 0.2649E+24
179 0.3739E+24
180 0.4345E+24
181 0.3956E+24
182 0.3973E+24
183 0.3989E+24
184 0.4480E+26
185 0.7157E+25
186 0.1027E+26
187 0.4977E+25
188 0.4833E+25
189 0.5107E+25
190 0.5784E+25
191 0.4947E+25
192 0.3065E+25
193 0.4011E+25
194 0.3791E+25
195 0.3128E+25
196 0.5732E+26
197 0.1811E+26
198 0.3001E+26
199 0.3602E+26
200 0.8397E+25
201 0.2245E+25
202 0.2528E+25
203 0.2776E+25
204 0.2689E+26
205 0.3199E+26
206 0.4921E+27
207 0.2729E+27
208 0.8629E+29
209 0.1008E+28
210 0.4227E+29
211 0.5417E+28
212 0.2607E+29
213 0.9517E+29
214 0.1723E+30
215 0.4818E+29
216 0.2958E+30
217 0.1239E+31
218 0.4307E+32
219 0.3235E+31
220 0.3338E+32
221 0.3338E+32
222 0.3340E+32
223 0.3340E+32
224 0.3340E+32
225 0.3340E+32
226 0.3341E+32
227 0.3341E+32
228 0.3341E+32
229 0.3343E+32
230 0.3343E+32
231 0.3341E+32
232 0.3346E+32
233 0.3351E+32
234 0.3312E+32
235 0.3176E+32
236 0.3166E+32
237 0.3154E+32
238 0.3207E+32
239 0.4065E+32
240 0.3338E+33
241 0.2643E+33
242 0.1670E+33
243 0.1067E+35
244 0.2440E+34
245 0.2974E+33
246 0.1610E+34
247 0.1251E+34
248 0.2430E+34
249 0.5367E+34
250 0.1653E+34
251 0.1190E+34
252 0.9299E+33
253 0.1243E+34
254 0.1260E+34
255 0.4369E+33
256 0.4604E+33
257 0.4711E+33
258 0.6165E+33
259 0.6407E+33
260 0.2137E+34
261 0.5191E+34
262 0.2219E+35
263 0.2531E+34
264 0.1943E+34
265 0.1834E+34
266 0.2493E+33
267 0.1845E+33
268 0.1879E+33
269 0.1855E+33
270 0.3002E+33
271 0.2771E+33
272 0.2271E+33
273 0.2857E+33
274 0.1435E+34
275 0.3133E+32
276 0.2550E+34
277 0.4228E+32
278 0.4218E+32
279 0.4218E+32
280 0.4219E+32
281 0.4239E+32
282 0.4239E+32
283 0.4240E+32
284 0.4239E+32
285 0.4237E+32
286 0.4240E+32
287 0.4204E+32
288 0.4233E+32
289 0.4231E+32
290 0.4233E+32
291 0.4232E+32
292 0.4235E+32
293 0.4234E+32
294 0.4234E+32
295 0.4243E+32
296 0.4243E+32
297 0.4250E+32
298 0.4250E+32
299 0.4270E+32
300 0.4244E+32
301 0.4174E+32
302 0.1434E+33
303 0.1202E+33
304 0.1298E+33
305 0.1163E+33
306 0.1781E+34
307 0.5018E+33
308 0.3213E+33
309 0.2067E+33
310 0.2148E+33
311 0.1102E+33
312 0.1615E+33
313 0.1367E+33
314 0.1317E+33
315 0.1321E+33
316 0.1328E+33
317 0.1286E+33
318 0.1285E+33
319 0.4729E+33
320 0.1557E+34
321 0.1254E+33
322 0.2204E+33
323 0.3765E+32
324 0.1323E+33
325 0.3991E+33
326 0.3251E+33
327 0.1429E+34
328 0.1158E+34
329 0.1515E+34
330 0.3820E+34
331 0.6888E+33
332 0.6051E+33
333 0.5931E+33
334 0.5937E+33
335 0.6796E+33
336 0.4805E+33
337 0.4609E+33
338 0.4607E+33
339 0.4605E+33
340 0.4614E+33
341 0.4641E+33
342 0.4663E+33
343 0.4657E+33
344 0.4657E+33
345 0.4661E+33
346 0.4650E+33
347 0.4714E+33
348 0.2230E+34
349 0.1245E+34
350 0.1414E+34
351 0.9135E+33
352 0.9067E+33
353 0.3486E+34
354 0.2728E+34
355 0.5584E+34
356 0.2411E+34
357 0.2053E+34
358 0.3551E+33
359 0.3534E+33
360 0.3283E+33
361 0.3322E+33
362 0.3341E+33
363 0.3313E+33
364 0.3294E+33
365 0.3344E+33
366 0.3340E+33
367 0.3327E+33
368 0.3320E+33
369 0.3319E+33
370 0.3307E+33
371 0.3294E+33
372 0.3330E+33
373 0.3332E+33
374 0.3344E+33
375 0.3252E+33
376 0.5773E+33
377 0.6800E+33
378 0.3903E+33
379 0.1011E+34
380 0.9596E+33
381 0.2176E+34
382 0.1545E+34
383 0.1671E+34
384 0.8867E+33
385 0.5806E+33
386 0.1114E+34
387 0.1586E+34
388 0.1413E+34
389 0.7029E+34
390 0.1886E+34
391 0.4841E+34
392 0.9673E+33
393 0.9335E+33
394 0.1502E+34
395 0.4288E+33
396 0.4961E+33
397 0.6194E+33
398 0.6324E+33
399 0.6404E+33
400 0.6505E+33
401 0.5038E+32
402 0.1109E+34
403 0.3037E+33
404 0.3080E+33
405 0.3540E+33
406 0.4642E+33
407 0.4653E+33
408 0.4705E+33
409 0.5054E+33
410 0.3335E+33
411 0.3421E+33
412 0.3420E+33
413 0.3427E+33
414 0.3426E+33
415 0.3427E+33
416 0.3427E+33
417 0.3434E+33
418 0.3425E+33
419 0.3431E+33
420 0.3435E+33
421 0.3425E+33
422 0.3435E+33
423 0.3439E+33
424 0.3438E+33
425 0.3435E+33
426 0.3440E+33
427 0.3444E+33
428 0.3450E+33
429 0.3457E+33
430 0.3451E+33
431 0.3461E+33
432 0.3463E+33
433 0.3457E+33
434 0.3466E+33
435 0.3464E+33
436 0.3466E+33
437 0.3456E+33
438 0.3462E+33
439 0.3470E+33
440 0.3444E+33
441 0.3684E+33
442 0.3580E+33
443 0.2683E+33
444 0.1900E+33
445 0.1872E+33
446 0.2068E+33
447 0.1963E+33
448 0.2638E+33
449 0.1150E+34
450 0.6129E+33
451 0.5969E+33
452 0.6022E+33
453 0.7884E+33
454 0.4374E+34
455 0.1739E+34
456 0.1105E+35
457 0.5117E+34
458 0.8620E+34
459 0.6370E+35
460 0.4444E+34
461 0.3932E+34
462 0.2076E+34
463 0.1627E+34
464 0.1401E+35
465 0.5818E+34
466 0.1031E+34
467 0.1959E+34
468 0.1816E+34
469 0.1797E+34
470 0.1881E+34
471 0.2717E+34
472 0.2605E+34
473 0.2101E+34
474 0.2709E+34
475 0.4640E+33
476 0.4828E+33
477 0.5365E+33
478 0.6214E+33
479 0.6155E+33
480 0.3524E+33
481 0.5317E+34
482 0.1234E+35
483 0.3776E+33
484 0.3640E+33
485 0.5351E+33
486 0.1886E+36
487 0.1495E+34
488 0.1094E+34
489 0.7059E+34
490 0.2551E+34
491 0.2598E+34
492 0.3918E+34
493 0.3186E+34
494 0.1829E+34
495 0.9378E+33
496 0.4490E+34
497 0.1012E+34
498 0.1009E+34
499 0.1920E+34
500 0.2112E+34
501 0.2112E+34
ERROR:: IterSolve: Failed convergence tolerances.
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.70159443+110 2.0000000 ) :: navier-stokes
FlowSolve: iter: 1 Assembly: (s) 0.94 0.94
FlowSolve: iter: 1 Solve: (s) 1.50 1.50
FlowSolve: Result Norm : 7.0159442990262395E+109
FlowSolve: Relative Change : 2.0000000000000000
ComputeChange: SS (ITER=4) (NRM,RELC): ( 0.70159443+110 2.0000000 ) :: navier-stokes
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: heat equation
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 1
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:
HeatSolve: Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
1 NaN
2 NaN
3 NaN
4 NaN
5 NaN
6 NaN
7 NaN
8 NaN
9 NaN
10 NaN
11 NaN
12 NaN
13 NaN
14 NaN
15 NaN
16 NaN
17 NaN
18 NaN
19 NaN
20 NaN
21 NaN
22 NaN
23 NaN
24 NaN
25 NaN
26 NaN
27 NaN
28 NaN
29 NaN
30 NaN
31 NaN
32 NaN
33 NaN
34 NaN
35 NaN
36 NaN
37 NaN
38 NaN
39 NaN
40 NaN
41 NaN
42 NaN
43 NaN
44 NaN
45 NaN
46 NaN
47 NaN
48 NaN
49 NaN
50 NaN
51 NaN
52 NaN
53 NaN
54 NaN
55 NaN
56 NaN
57 NaN
58 NaN
59 NaN
60 NaN
61 NaN
62 NaN
63 NaN
64 NaN
65 NaN
66 NaN
67 NaN
68 NaN
69 NaN
70 NaN
71 NaN
72 NaN
73 NaN
74 NaN
75 NaN
76 NaN
77 NaN
78 NaN
79 NaN
80 NaN
81 NaN
82 NaN
83 NaN
84 NaN
85 NaN
86 NaN
87 NaN
88 NaN
89 NaN
90 NaN
91 NaN
92 NaN
93 NaN
94 NaN
95 NaN
96 NaN
97 NaN
98 NaN
99 NaN
100 NaN
101 NaN
102 NaN
103 NaN
104 NaN
105 NaN
106 NaN
107 NaN
108 NaN
109 NaN
110 NaN
111 NaN
112 NaN
113 NaN
114 NaN
115 NaN
116 NaN
117 NaN
118 NaN
119 NaN
120 NaN
121 NaN
122 NaN
123 NaN
124 NaN
125 NaN
126 NaN
127 NaN
128 NaN
129 NaN
130 NaN
131 NaN
132 NaN
133 NaN
134 NaN
135 NaN
136 NaN
137 NaN
138 NaN
139 NaN
140 NaN
141 NaN
142 NaN
143 NaN
144 NaN
145 NaN
146 NaN
147 NaN
148 NaN
149 NaN
150 NaN
151 NaN
152 NaN
153 NaN
154 NaN
155 NaN
156 NaN
157 NaN
158 NaN
159 NaN
160 NaN
161 NaN
162 NaN
163 NaN
164 NaN
165 NaN
166 NaN
167 NaN
168 NaN
169 NaN
170 NaN
171 NaN
172 NaN
173 NaN
174 NaN
175 NaN
176 NaN
177 NaN
178 NaN
179 NaN
180 NaN
181 NaN
182 NaN
183 NaN
184 NaN
185 NaN
186 NaN
187 NaN
188 NaN
189 NaN
190 NaN
191 NaN
192 NaN
193 NaN
194 NaN
195 NaN
196 NaN
197 NaN
198 NaN
199 NaN
200 NaN
201 NaN
202 NaN
203 NaN
204 NaN
205 NaN
206 NaN
207 NaN
208 NaN
209 NaN
210 NaN
211 NaN
212 NaN
213 NaN
214 NaN
215 NaN
216 NaN
217 NaN
218 NaN
219 NaN
220 NaN
221 NaN
222 NaN
223 NaN
224 NaN
225 NaN
226 NaN
227 NaN
228 NaN
229 NaN
230 NaN
231 NaN
232 NaN
233 NaN
234 NaN
235 NaN
236 NaN
237 NaN
238 NaN
239 NaN
240 NaN
241 NaN
242 NaN
243 NaN
244 NaN
245 NaN
246 NaN
247 NaN
248 NaN
249 NaN
250 NaN
251 NaN
252 NaN
253 NaN
254 NaN
255 NaN
256 NaN
257 NaN
258 NaN
259 NaN
260 NaN
261 NaN
262 NaN
263 NaN
264 NaN
265 NaN
266 NaN
267 NaN
268 NaN
269 NaN
270 NaN
271 NaN
272 NaN
273 NaN
274 NaN
275 NaN
276 NaN
277 NaN
278 NaN
279 NaN
280 NaN
281 NaN
282 NaN
283 NaN
284 NaN
285 NaN
286 NaN
287 NaN
288 NaN
289 NaN
290 NaN
291 NaN
292 NaN
293 NaN
294 NaN
295 NaN
296 NaN
297 NaN
298 NaN
299 NaN
300 NaN
301 NaN
302 NaN
303 NaN
304 NaN
305 NaN
306 NaN
307 NaN
308 NaN
309 NaN
310 NaN
311 NaN
312 NaN
313 NaN
314 NaN
315 NaN
316 NaN
317 NaN
318 NaN
319 NaN
320 NaN
321 NaN
322 NaN
323 NaN
324 NaN
325 NaN
326 NaN
327 NaN
328 NaN
329 NaN
330 NaN
331 NaN
332 NaN
333 NaN
334 NaN
335 NaN
336 NaN
337 NaN
338 NaN
339 NaN
340 NaN
341 NaN
342 NaN
343 NaN
344 NaN
345 NaN
346 NaN
347 NaN
348 NaN
349 NaN
350 NaN
351 NaN
352 NaN
353 NaN
354 NaN
355 NaN
356 NaN
357 NaN
358 NaN
359 NaN
360 NaN
361 NaN
362 NaN
363 NaN
364 NaN
365 NaN
366 NaN
367 NaN
368 NaN
369 NaN
370 NaN
371 NaN
372 NaN
373 NaN
374 NaN
375 NaN
376 NaN
377 NaN
378 NaN
379 NaN
380 NaN
381 NaN
382 NaN
383 NaN
384 NaN
385 NaN
386 NaN
387 NaN
388 NaN
389 NaN
390 NaN
391 NaN
392 NaN
393 NaN
394 NaN
395 NaN
396 NaN
397 NaN
398 NaN
399 NaN
400 NaN
401 NaN
402 NaN
403 NaN
404 NaN
405 NaN
406 NaN
407 NaN
408 NaN
409 NaN
410 NaN
411 NaN
412 NaN
413 NaN
414 NaN
415 NaN
416 NaN
417 NaN
418 NaN
419 NaN
420 NaN
421 NaN
422 NaN
423 NaN
424 NaN
425 NaN
426 NaN
427 NaN
428 NaN
429 NaN
430 NaN
431 NaN
432 NaN
433 NaN
434 NaN
435 NaN
436 NaN
437 NaN
438 NaN
439 NaN
440 NaN
441 NaN
442 NaN
443 NaN
444 NaN
445 NaN
446 NaN
447 NaN
448 NaN
449 NaN
450 NaN
451 NaN
452 NaN
453 NaN
454 NaN
455 NaN
456 NaN
457 NaN
458 NaN
459 NaN
460 NaN
461 NaN
462 NaN
463 NaN
464 NaN
465 NaN
466 NaN
467 NaN
468 NaN
469 NaN
470 NaN
471 NaN
472 NaN
473 NaN
474 NaN
475 NaN
476 NaN
477 NaN
478 NaN
479 NaN
480 NaN
481 NaN
482 NaN
483 NaN
484 NaN
485 NaN
486 NaN
487 NaN
488 NaN
489 NaN
490 NaN
491 NaN
492 NaN
493 NaN
494 NaN
495 NaN
496 NaN
497 NaN
498 NaN
499 NaN
500 NaN
501 NaN
ERROR:: IterSolve: Failed convergence tolerances.
ERROR:: ComputeChange: Norm of solution appears to be NaN
Re: Non-Boussinesq Convection
Well, I wasn't able to get the ILU preconditioning to work. Anytime I had it enabled the simulation would stop after the first few timesteps. I played around with some of the tolerances, but in the end I set both solvers back to diagonal preconditioning. I also wasn't able to get non-linear timestepping working, but it was more successful than the ILU preconditioning. With more practice I could probably get this working.
However, I was able to get the simulation to run just fine. I ended up creating a restart file from this model where I only had the heat equation active. I then used this restart file to run the model with the Navier-Stokes equations enabled. I thought it might make things more stable. I also adjusted the maximum number of non-linear iterations. Honestly, it was probably the last change that stabilized the model. See below.
I'm going to try to refine this a bit and see if I can get the model to work without using the restart file for my initial conditions.
Peter, you asked me whether the linear or non-linear system diverges first. Are you able to determine this from the few lines of the solver log that I attached? Would you be able to describe how to determine this from it? I only have a half a clue about how to read the solver log. I assume it was the non-linear convergence that was the problem just based on the changes I made.
However, I was able to get the simulation to run just fine. I ended up creating a restart file from this model where I only had the heat equation active. I then used this restart file to run the model with the Navier-Stokes equations enabled. I thought it might make things more stable. I also adjusted the maximum number of non-linear iterations. Honestly, it was probably the last change that stabilized the model. See below.
I'm going to try to refine this a bit and see if I can get the model to work without using the restart file for my initial conditions.
Peter, you asked me whether the linear or non-linear system diverges first. Are you able to determine this from the few lines of the solver log that I attached? Would you be able to describe how to determine this from it? I only have a half a clue about how to read the solver log. I assume it was the non-linear convergence that was the problem just based on the changes I made.
-
- Site Admin
- Posts: 4828
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Non-Boussinesq Convection
Hi
You could perhaps sent the case for inspection. It should not be so much more difficult from the Boussinesq case.
-Peter
You could perhaps sent the case for inspection. It should not be so much more difficult from the Boussinesq case.
-Peter