Add new setup for 3D flame propagation through a tube (#2551)

Exec/science/flame_tube/GNUmakefile

@@ -0,0 +1,38 @@
+PRECISION = DOUBLE
+PROFILE = FALSE
+
+DEBUG = FALSE
+
+DIM = 3
+
+COMP = gnu
+
+USE_MPI = TRUE
+
+USE_GRAV = TRUE
+USE_REACT = TRUE
+
+USE_ROTATION = FALSE
+USE_DIFFUSION = TRUE
+
+# define the location of the CASTRO top directory
+CASTRO_HOME := ../../..
+
+USE_JACOBIAN_CACHING = TRUE
+USE_CXX_MODEL_PARSER = TRUE
+NUM_MODELS := 2
+
+# This sets the EOS directory in $(MICROPHYSICS_HOME)/eos
+EOS_DIR := helmholtz
+
+# This sets the network directory in $(MICROPHYSICS_HOME)/networks
+NETWORK_DIR := aprox13
+
+INTEGRATOR_DIR := VODE
+
+CONDUCTIVITY_DIR := stellar
+
+Bpack := ./Make.package
+Blocs := .
+
+include $(CASTRO_HOME)/Exec/Make.Castro

Exec/science/flame_tube/Make.package

@@ -0,0 +1,2 @@
+CEXE_headers += initial_model.H
+

Exec/science/flame_tube/README.md

@@ -0,0 +1,5 @@
+# flame_tube
+
+This is a slightly modified version of flame_wave that models a 3-d flame
+propagating through a tube shaped domain, with a periodic boundary along the
+transverse horizontal direction (y).

Exec/science/flame_tube/_prob_params

@@ -0,0 +1,70 @@
+
+dtemp real 3.81e8_rt y
+
+x_half_max real 1.2e5_rt y
+
+x_half_width real 3.6e4_rt y
+
+# cutoff mass fraction of the first species for refinement
+X_min real 1.e-4_rt y
+
+# do we dynamically refine based on density? or based on height?
+tag_by_density integer 1 y
+
+# used for tagging if tag_by_density = 1
+cutoff_density real 500.e0_rt y
+
+# used if we are refining based on height rather than density
+refine_height real 3600 y
+
+dx_model real 10.0_rt y
+
+T_hi real 5.e8_rt y
+
+T_star real 1.e8_rt y
+
+T_lo real 5.e7_rt y
+
+dens_base real 2.e6_rt y
+
+H_star real 500.e0_rt y
+
+atm_delta real 25.e0_rt y
+
+fuel1_name character "helium-4" y
+
+fuel2_name character "" y
+
+fuel3_name character "" y
+
+fuel4_name character "" y
+
+ash1_name character "iron-56" y
+
+ash2_name character "" y
+
+ash3_name character "" y
+
+fuel1_frac real 1.0_rt y
+
+fuel2_frac real 0.0_rt y
+
+fuel3_frac real 0.0_rt y
+
+fuel4_frac real 0.0_rt y
+
+ash1_frac real 1.0_rt y
+
+ash2_frac real 0.0_rt y
+
+ash3_frac real 0.0_rt y
+
+low_density_cutoff real 1.e-4_rt y
+
+smallx real 1.e-10_rt y
+
+x_refine_distance real 1.e30_rt y
+
+max_hse_tagging_level integer 2 y
+
+max_base_tagging_level integer 1 y

Exec/science/flame_tube/analysis/vis_3d/andes-single.submit

@@ -0,0 +1,34 @@
+#!/bin/bash
+#SBATCH -A ast106
+#SBATCH -J flame_tube_vr
+#SBATCH -N 1
+#SBATCH -t 4:00:00
+
+set -u
+
+cd $SLURM_SUBMIT_DIR
+
+source "/ccs/proj/ast106/$USER/mambaforge_$(uname -m)/etc/profile.d/conda.sh"
+conda activate andes_yt_dev
+
+plotfiles=(run_256/flame_tube_25cm_smallplt*)
+plotfiles=(run_256_extra_plotfiles/flame_tube_25cm_smallplt0039200)
+#plotfile=flame_wave_1000Hz_25cm_smallplt207887
+#plotfile=flame_wave_1000Hz_25cm_smallplt40842
+
+#for i in flame_wave_1000Hz_25cm_smallplt*[0-9]
+todo=()
+for f in "${plotfiles[@]}"; do
+ dest=${f/run_/analysis_}
+ # check the last image generated for each plotfile
+ #if ! [[ -f "${dest}_enuc_annotated_top.png" ]]; then
+ todo+=("$f")
+ #fi
+done
+
+export OMP_NUM_THREADS=16
+if [[ ${#todo[@]} -gt 0 ]]; then
+ srun --cpu-bind=no python ~/dev/Castro/Exec/science/flame_tube/analysis/vis_3d/vol-xrb.py "${todo[@]}"
+fi
+
+echo "done!"

Exec/science/flame_tube/analysis/vis_3d/andes-slice.submit

@@ -0,0 +1,17 @@
+#!/bin/bash
+#SBATCH -A ast106
+#SBATCH -J plots
+#SBATCH -N 1
+#SBATCH -t 8:00:00
+#SBATCH -p gpu
+
+cd $SLURM_SUBMIT_DIR
+
+#plotfile=flame_wave_1000Hz_25cm_smallplt207887
+#plotfile=flame_wave_1000Hz_25cm_smallplt40842
+
+source activate andes_env
+
+srun python slice_vertical.py flame_wave_1000Hz_25cm_smallplt174562
+
+

Exec/science/flame_tube/analysis/vis_3d/andes.submit

@@ -0,0 +1,26 @@
+#!/bin/bash
+#SBATCH -A ast106
+#SBATCH -J plots
+#SBATCH -N 1
+#SBATCH -t 8:00:00
+
+cd $SLURM_SUBMIT_DIR
+
+source "/ccs/proj/ast106/$USER/mambaforge_$(uname -m)/etc/profile.d/conda.sh"
+conda activate andes_yt_dev
+
+plotfile=run_256_extra_plotfiles/flame_tube_25cm_smallplt0039200
+#plotfile=flame_wave_1000Hz_25cm_smallplt207887
+#plotfile=flame_wave_1000Hz_25cm_smallplt40842
+
+#for i in flame_wave_1000Hz_25cm_smallplt*[0-9]
+export OMP_NUM_THREADS=16
+for i in $plotfile; do
+ dest=${i/run_/analysis_}
+ if ! [[ -f "${dest}_abar_annotated_top.png" ]]; then
+ srun --cpu-bind=no python ~/dev/Castro/Exec/science/flame_tube/analysis/vis_3d/vol-xrb-abar.py "${i}"
+ srun --cpu-bind=no python ~/dev/Castro/Exec/science/flame_tube/analysis/vis_3d/vol-xrb-enuc.py "${i}"
+ fi
+done
+
+echo "done!"

Exec/science/flame_tube/analysis/vis_3d/slice_vertical.py

@@ -0,0 +1,97 @@
+import argparse
+import os
+import sys
+
+import matplotlib
+import matplotlib.pyplot as plt
+import numpy as np
+from mpl_toolkits.axes_grid1 import ImageGrid
+
+import yt
+from yt.frontends.boxlib.api import CastroDataset
+# assume that our data is in CGS
+from yt.units import amu, cm
+
+matplotlib.use('agg')
+
+
+
+
+def doit(plotfiles):
+
+ plotfile = sys.argv[1]
+ ds = CastroDataset(plotfile)
+
+ xmin = ds.domain_left_edge[0]
+ xmax = ds.domain_right_edge[0]
+ xctr = 0.5 * (xmin + xmax)
+ L_x = (2./3.) * (xmax - xmin)
+
+ ymin = ds.domain_left_edge[1]
+ ymax = ds.domain_right_edge[1]
+ yctr = 0.5*(ymin + ymax)
+
+ zmin = 0.0*cm
+ zmax = 1.e4*cm
+
+ zctr = 0.5*(zmin + zmax)
+ L_z = zmax - zmin
+
+ fig = plt.figure()
+
+ grid = ImageGrid(fig, 111, nrows_ncols=(len(plotfiles), 1),
+ axes_pad=0.25, label_mode="L",
+ cbar_mode="single", cbar_size="0.5%")
+
+ for i, pf in enumerate(plotfiles):
+
+ ds = CastroDataset(pf)
+
+ f = "abar"
+
+ sp = yt.SlicePlot(ds, "y", f, origin="native", center=[xctr, yctr, zctr],
+ width=[L_z, L_x], fontsize="9")
+ sp.set_buff_size((4800,4800))
+ sp.swap_axes()
+
+ sp.set_zlim(f, 4, 5)
+ sp.set_log(f, False)
+ sp.set_cmap(f, "plasma_r")
+
+ sp.set_axes_unit("cm")
+
+ sp.annotate_text((0.8, 0.9), "{:5.2f} ms".format(1000.0*float(ds.current_time.in_cgs())),
+ coord_system="axis", text_args={"color": "black", "size": 9})
+
+ plot = sp.plots["abar"]
+ plot.figure = fig
+ plot.axes = grid[i].axes
+ plot.cax = grid.cbar_axes[i]
+ if i < len(plotfiles)-1:
+ grid[i].axes.xaxis.offsetText.set_visible(False)
+
+ sp._setup_plots()
+
+ fig.set_size_inches(10.0, 3.5)
+ fig.subplots_adjust(left=0.05, right=0.95, top=0.95, bottom=0.05)
+ fig.savefig("time_series_3D.pdf")
+
+if __name__ == "__main__":
+
+ p = argparse.ArgumentParser()
+
+ p.add_argument("--skip", type=int, default=1,
+ help="interval between plotfiles")
+ p.add_argument("plotfiles", type=str, nargs="+",
+ help="list of plotfiles to plot")
+
+ args = p.parse_args()
+
+ plot_prefix = args.plotfiles[0].split("plt")[0] + "plt"
+ plot_nums = sorted([p.split("plt")[1] for p in args.plotfiles], key=int)
+
+ plotfiles = []
+ for n in range(0, len(plot_nums), args.skip):
+ plotfiles.append("{}{}".format(plot_prefix, plot_nums[n]))
+
+ doit(plotfiles)