Changes between Version 31 and Version 32 of 1DPulsedJets

Timestamp:

03/06/12 14:52:07 (13 years ago)

Author:

ehansen

Comment:

—

1DPulsedJets

@@ -141 +141 @@
 ==== Results ====
 Here we look at what post-shock temperature the simulation gives as a function of the effective resolution.  The effective resolution = (# cells) x 2^(# AMR levels)^.  Each run used a base grid of 100 cells, so the only thing that changed from run to run was the number of AMR levels.  For some reason, I could not get it to run with more than 7 levels, so the 25600 run is actually a base grid of 200 cells with 7 levels.  For Tps I will use the highest temperature that I get from the simulation.  The highest temperature always occurs early on in the simulation and then the post-shock temperature slowly decreases over time, but stays relatively constant at some lower temperature.
+[[BR]]
+Z
 ||= Effective Resolution =||= cells/lcool =||= Tps (10^3^ K) =||= Relative Error (%) =||
 ||= 100 =||= 10.68 =||= 34.37 =||= 41.06 =||
@@ -158 +160 @@
 
 My next idea was to revert back to just DM cooling instead of the new Z cooling.  Perhaps, the Z cooling routine is not implemented quite right, and there is some "double counting" for cooling strength or something like this.  So here is a data table where I used DM instead of Z.  The improvement column is the decrease in relative error from the comparable Z cooling run.
-
+[[BR]]
+DM
 ||= Effective Resolution =||= cells/lcool =||= Tps (10^3^ K) =||= Relative Error (%) =||
-||= 100 =||= 10.68 =||=  =||=  =||
-||= 200 =||= 21.36 =||=  =||=  =||
-||= 400 =||= 42.71 =||=  =||=  =||
+||= 100 =||= 10.68 =||= 25.71 =||= 55.91 =||
+||= 200 =||= 21.36 =||= 31.63 =||= 45.76 =||
+||= 400 =||= 42.71 =||= 41.07 =||= 29.57 =||
 ||= 800 =||= 85.43 =||= 47.65 =||= 18.28 =||
 ||= 1600 =||= 170.85 =||= 51.37 =||= 11.91 =||
@@ -168 +171 @@
 ||= 6400 =||= 683.40 =||= 53.98 =||= 7.43 =||
 ||= 12800 =||= 1366.81 =||= 54.14 =||= 7.16 =||
-||= 25600 =||= 2733.61 =||= =||= =||
+||= 25600 =||= 2733.61 =||= 58.16 =||= 0.26 =||
 
-The DM cooling gets closer to what I expect but the Tps is still a little low.  In both cases, Tps seems to be converging to some other value different from my expected value of 58.3123.  However, when I do a run with no cooling, I get post-shock temperatures reaching 58.3123, and some even higher...closer to 60.75 with just 100 cells and no AMR.  So I decided to take a closer look at the adiabatic case to see if I could figure out the discrepancy.  Again, here is the data for no cooling:
-
-
+The DM cooling gets closer to what I expect but the Tps is still a little low.  In both cases, Tps seems to be converging to some other value different from my expected value of 58.3123.  However, when I do a run with no cooling, I get post-shock temperatures reaching 58.3123, and some even higher...closer to 60.97 with just 100 cells and no AMR.  So I decided to take a closer look at the adiabatic case to see if I could figure out the discrepancy.  Again, here is the data for no cooling:
+[[BR]]
+None
 ||= Effective Resolution =||= cells/lcool =||= Tps (10^3^ K) =||= Relative Error (%) =||
-||= 100 =||= 10.68 =||=  =||=  =||
-||= 200 =||= 21.36 =||=  =||=  =||
-||= 400 =||= 42.71 =||=  =||=  =||
-||= 800 =||= 85.43 =||=  =||=  =||
-||= 1600 =||= 170.85 =||=  =||=  =||
-||= 3200 =||= 341.70 =||=  =||=  =||
-||= 6400 =||= 683.40 =||=  =||=  =||
-||= 12800 =||= 1366.81 =||=  =||=  =||
-||= 25600 =||= 2733.61 =||=  =||=  =||
+||= 100 =||= 10.68 =||= 60.97 =||= 4.56 =||
+||= 200 =||= 21.36 =||= 61.09 =||= 4.76 =||
+||= 400 =||= 42.71 =||= 61.00 =||= 4.61 =||
+||= 800 =||= 85.43 =||= 61.06 =||= 4.71 =||
+||= 1600 =||= 170.85 =||= 61.76 =||= 5.91 =||
+||= 3200 =||= 341.70 =||= 63.45 =||= 8.81 =||
+||= 6400 =||= 683.40 =||= 64.17 =||= 10.04 =||
+||= 12800 =||= 1366.81 =||= 65.20 =||= 11.81 =||
+||= 25600 =||= 2733.61 =||= 64.16 =||= 10.02 =||
 
 You can see that these simulations actually give higher than expected post-shock temperatures.  If I used these values as the "correct" values then the post-shock temperatures from the cooling simulations are even further off.