Extrusion v3

.gitignore

@@ -54,6 +54,7 @@ docs/simulation/readme.md
 docs/post_process/readme.md
 docs/api/readme.md
 
+examples/**/IC/*.dat
 examples/*batch/*/
 examples/**/D/*
 examples/**/p*

docs/documentation/case.md

@@ -262,8 +262,8 @@ The code provides three pre-built patches for dimensional extrusion of initial c
 - `case(270)`: Extrude 1D data to 2D domain
 - `case(370)`: Extrude 2D data to 3D domain
 
-Setup: Only requires specifying `init_dir` and filename pattern via `zeros_default`. Grid dimensions are automatically detected from the data files.
-Implementation: All variables and file handling are managed in `src/common/include/ExtrusionHardcodedIC.fpp` with no manual grid configuration needed.
+Setup: Only requires specifying `files_dir` and filename pattern via `file_extension`. The files are located, for example, at `examples/1D_flamelet/IC`, and their format is `prim.XX.YY.file_extension.dat`.
+Implementation: All variables and file handling are managed in the `case.py` file of the simulation.
 Usage: Ideal for initializing simulations from lower-dimensional solutions, enabling users to add perturbations or modifications to the base extruded fields for flow instability studies.
 
 #### Parameter Descriptions

examples/1D_flamelet/case.py

@@ -0,0 +1,73 @@
+#!/usr/bin/env python3
+import argparse
+import json
+import math
+import os
+
+import cantera as ct
+
+current_dir = os.path.dirname(os.path.abspath(__file__))
+ctfile = "h2o2.yaml"
+sol_L = ct.Solution(ctfile)
+sol_L.TPX = 300, 8000, "O2:2,N2:2,H2O:5"
+L = 0.016
+Nx = 1199
+dx = L / Nx
+dt = 1e-8
+Tend = 0.60e-3
+NT = int(Tend / dt)
+SAVE_COUNT = 1000
+NS = 1000
+case = {
+    "run_time_info": "T",
+    "x_domain%beg": -L / 2,
+    "x_domain%end": +L / 2,
+    "m": Nx,
+    "n": 0,
+    "p": 0,
+    "dt": float(dt),
+    "t_step_start": 0,
+    "t_step_stop": NT,
+    "t_step_save": NS,
+    "t_step_print": 100,
+    "parallel_io": "F",
+    "model_eqns": 2,
+    "num_fluids": 1,
+    "num_patches": 1,
+    "mpp_lim": "F",
+    "mixture_err": "F",
+    "time_stepper": 3,
+    "weno_order": 5,
+    "weno_eps": 1e-16,
+    "weno_avg": "F",
+    "mapped_weno": "T",
+    "mp_weno": "T",
+    "riemann_solver": 2,
+    "wave_speeds": 2,
+    "avg_state": 1,
+    "bc_x%beg": -8,
+    "bc_x%end": -8,
+    "viscous": "F",
+    "files_dir": os.path.join(current_dir, "IC"),
+    "file_extension": "000000",
+    "chemistry": "T",
+    "chem_params%diffusion": "T",
+    "chem_params%reactions": "T",
+    "chem_params%transport_model": 2,
+    "format": 1,
+    "precision": 2,
+    "prim_vars_wrt": "T",
+    "patch_icpp(1)%geometry": 1,
+    "patch_icpp(1)%hcid": 170,
+    "patch_icpp(1)%x_centroid": 0,
+    "patch_icpp(1)%length_x": L,
+    "patch_icpp(1)%vel(1)": "0",
+    "patch_icpp(1)%pres": 1.01325e5,
+    "patch_icpp(1)%alpha(1)": 1,
+    "patch_icpp(1)%alpha_rho(1)": "1",
+    "fluid_pp(1)%gamma": 1.0e00 / (1.5e00 - 1.0e00),
+    "fluid_pp(1)%pi_inf": 0,
+    "cantera_file": ctfile,
+}
+if __name__ == "__main__":
+    print(json.dumps(case))