Talk:BA Fynn Wawrzyniak

From Arbeitsgruppe Kuiper
Revision as of 08:51, 26 April 2024 by Fynn.W (talk | contribs) (logs für dt=1e-10,dt=1e-15,dt=1e-20)
Jump to navigation Jump to search

ConsToPrim(): p(E) < 0

Hilf ein Absenken des initialen Zeitschritts? --Lothar (talk) 07:53, 26 April 2024 (CEST)


Es scheint, dass das Absenken des initialen Zeitschritts nicht hilft. Ich habe einfach mal die logs für verschiedene Einstellung des initialen Zeitschritts beigefügt.

Für \( dt=10^{-10} \)

> Starting computation... 

step:0; t = 0.0000e+00; dt = 1.0000e-10; 0.0 %
        [Mach = 49.994517]
! ConsToPrim(): p(E) < 0 (-2.34e-05),   @step = 66 (stage = 1); [i,j = 2, 73], [x1,x2 =  0.000001, 2.246239]
! ConsToPrim(): p(E) < 0 (-1.27e-05),   @step = 66 (stage = 1); [i,j = 2, 74], [x1,x2 =  0.000001, 2.277655]

Für \( dt=10^{-15} \)

> Starting computation... 

step:0; t = 0.0000e+00; dt = 1.0000e-15; 0.0 %
        [Mach = 49.994517]
step:100; t = 1.3780e-10; dt = 1.3781e-11; 0.0 %
          [Mach = 49.995660]
! ConsToPrim(): p(E) < 0 (-3.49e-06),   @step = 186 (stage = 2); [i,j = 2, 73], [x1,x2 =  0.000001, 2.246239]
! ConsToPrim(): p(E) < 0 (-1.84e-06),   @step = 187 (stage = 1); [i,j = 2, 72], [x1,x2 =  0.000001, 2.214823]


Für \( dt=10^{-20} \)

> Starting computation... 

step:0; t = 0.0000e+00; dt = 1.0000e-20; 0.0 %
        [Mach = 49.994517]
step:100; t = 1.3780e-15; dt = 1.3781e-16; 0.0 %
          [Mach = 49.994517]
step:200; t = 1.8991e-11; dt = 1.8991e-12; 0.0 %
          [Mach = 49.994517]
step:300; t = 1.3811e-07; dt = 5.1143e-09; 0.0 %
          [Mach = 57.016462]
! ConsToPrim(): p(E) < 0 (-7.80e-06),   @step = 307 (stage = 1); [i,j = 2, 73], [x1,x2 =  0.000001, 2.246239]
! ConsToPrim(): p(E) < 0 (-8.17e-06),   @step = 307 (stage = 2); [i,j = 2, 73], [x1,x2 =  0.000001, 2.246239]

--Fynn.W (talk) 09:51, 26 April 2024 (CEST)