-
Notifications
You must be signed in to change notification settings - Fork 98
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge branch 'development' into exact_riemann_merge
- Loading branch information
Showing
45 changed files
with
3,352 additions
and
119 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,3 +1,4 @@ | ||
| ALLOW_GPU_PRINTF | ||
| AMREX_DEBUG | ||
| AMREX_PARTICLES | ||
| AMREX_SPACEDIM | ||
|
|
||
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
132 changes: 66 additions & 66 deletions
132
Exec/science/Detonation/ci-benchmarks/sdc_det_plt00040_extrema.out
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,47 @@ | ||
| # Shock Burning Experiments | ||
|
|
||
| This directory is meant to explore shock burning with detonation. Compile as: | ||
|
|
||
| ``` | ||
| make USE_SIMPLIFIED_SDC=TRUE USE_SHOCK_VAR=TRUE NETWORK_DIR=aprox13 -j 4 | ||
| ``` | ||
|
|
||
| Then the script `setup_runs.py` will setup a suite of simulations with | ||
| the following resolutions into separate directories (using the | ||
| `inputs-shock-burn.template`): | ||
|
|
||
|
|
||
| | resolution | base grid | levels (4x jumps) | | ||
| | ------------ | ----------- | ------------------- | | ||
| | 24 km | 48 | 1 | | ||
| | 12 km | 96 | 1 | | ||
| | 6 km | 192 | 1 | | ||
| | 3 km | 384 | 1 | | ||
| | 1.5 km | 768 | 1 | | ||
| | 0.1875 km | 6144 | 1 | | ||
| | 2343.74 cm | 12288 | 2 | | ||
|
|
||
| you can set the value of the shock detection threshold there | ||
| and the directory names will reflect that setting. | ||
|
|
||
| ## plotting | ||
|
|
||
| The following scripts can make useful plots (some use the | ||
| `detonation.py` module as support): | ||
|
|
||
| * `det_speed_comp.py` : plot detonation speed vs. resolution, using | ||
| simple differencing to estimate the detonation speed from the last 3 | ||
| plotfiles. | ||
|
|
||
| * `profile_compare.py` : given a list of pruns (from different | ||
| resolutions), make a plot showing all of their profiles together. | ||
|
|
||
| * `profiles.py` : for a single run, plot profiles of T and enuc for | ||
| several times. | ||
|
|
||
| * `show_shock_flag.py` : simply plot the shock variable on top of T | ||
| and enuc profiles. | ||
|
|
||
| * `zoom_summary.py` : given a list of runs (from different | ||
| resolutions), plot the last plotfile from each run, zoomed in on | ||
| where the peak energy generation is. |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,74 @@ | ||
| #!/usr/bin/env python | ||
|
|
||
| import matplotlib.pyplot as plt | ||
|
|
||
| import detonation | ||
|
|
||
| runs = [("res24.0km", 24), | ||
| ("res12.0km", 12), | ||
| ("res6.0km", 6), | ||
| ("res3.0km", 3), | ||
| ("res1.5km", 1.5), | ||
| ("res0.1875km", 0.1875)] #, | ||
| #("res0.024km", 0.024)] #, | ||
| #("res0.003km", 0.003)] | ||
|
|
||
| nsb1_runs = [("res24.0km_noshockburn_1", 24), | ||
| ("res12.0km_noshockburn_1", 12), | ||
| ("res6.0km_noshockburn_1", 6), | ||
| ("res3.0km_noshockburn_1", 3), | ||
| ("res1.5km_noshockburn_1", 1.5), | ||
| ("res0.1875km_noshockburn_1", 0.1875)] #, | ||
|
|
||
| nsb23_runs = [("res24.0km_noshockburn_0.666", 24), | ||
| ("res12.0km_noshockburn_0.666", 12), | ||
| ("res6.0km_noshockburn_0.666", 6), | ||
| ("res3.0km_noshockburn_0.666", 3), | ||
| ("res1.5km_noshockburn_0.666", 1.5), | ||
| ("res0.1875km_noshockburn_0.666", 0.1875)] #, | ||
|
|
||
| res = [] | ||
| v = [] | ||
| dv = [] | ||
|
|
||
| for ddir, dx in runs: | ||
| res.append(dx) | ||
| d = detonation.Detonation(ddir) | ||
| v.append(d.v) | ||
| dv.append(d.v_sigma) | ||
|
|
||
| nsb23_res = [] | ||
| nsb23_v = [] | ||
| nsb23_dv = [] | ||
|
|
||
| for ddir, dx in nsb23_runs: | ||
| nsb23_res.append(dx) | ||
| d = detonation.Detonation(ddir) | ||
| nsb23_v.append(d.v) | ||
| nsb23_dv.append(d.v_sigma) | ||
|
|
||
| nsb1_res = [] | ||
| nsb1_v = [] | ||
| nsb1_dv = [] | ||
|
|
||
| for ddir, dx in nsb1_runs: | ||
| nsb1_res.append(dx) | ||
| d = detonation.Detonation(ddir) | ||
| nsb1_v.append(d.v) | ||
| nsb1_dv.append(d.v_sigma) | ||
|
|
||
|
|
||
| fig, ax = plt.subplots() | ||
|
|
||
| ax.errorbar(res, v, yerr=dv, fmt="o", label="burning in shocks allowed") | ||
| ax.errorbar(nsb23_res, nsb23_v, yerr=nsb23_dv, fmt="d", label="no shock burning (f=2/3)") | ||
| ax.errorbar(nsb1_res, nsb1_v, yerr=nsb1_dv, fmt="d", label="no shock burning (f=1)") | ||
|
|
||
| ax.set_xlabel(r"$\Delta x$ (km)") | ||
| ax.set_ylabel(r"$v_\mathrm{det}$ (cm / s)") | ||
|
|
||
| ax.legend() | ||
|
|
||
| ax.set_xscale("log") | ||
|
|
||
| fig.savefig("det_speeds.png") |
Oops, something went wrong.