Request for Guidance on Shock- Flame Interaction Setup in PeleC

[Alex-jian522]

Hello, everyone!
I’m working on setting up a case for shock-flame interaction in PeleC, with the shock wave defined as an initial condition. Is there a tutorial or example case I can refer to for guidance on how to do this?
Thank you very much for your help!

[marchdf]

There's the Sod regression test and the PMF test. Combining those should get you to where you want to be.

• https://github.com/AMReX-Combustion/PeleC/tree/development/Exec/RegTests/Sod
• https://github.com/AMReX-Combustion/PeleC/tree/development/Exec/RegTests/PMF

[Alex-jian522]

Thank you very much!
I have another question: how can I set the lower boundary as NSCBC to prevent the shock from reflecting at this boundary? Currently, I use "hard" and do nothing in the bcnormal of prob.H, as shown in the following.

pelec.lo_bc = "Interior" "Hard"
pelec.hi_bc = "Interior" "NoSlipWall"

void
bcnormal(
const amrex::Real* /x[AMREX_SPACEDIM]/,
const amrex::Real* /s_int[NVAR]/,
amrex::Real s_ext[NVAR],
const int idir,
const int sgn,
const amrex::Real /time/,
amrex::GeometryData const& geomdata,
ProbParmDevice const& prob_parm,
const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& /turb_fluc/)
{
const amrex::Real* prob_lo = geomdata.ProbLo();
const amrex::Real* prob_hi = geomdata.ProbHi();
amrex::GpuArray<amrex::Real, NUM_SPECIES + 4> pmf_vals = {{0.0}};
amrex::Real u[3] = {0.0};
amrex::Real molefrac[NUM_SPECIES] = {0.0};
amrex::Real massfrac[NUM_SPECIES] = {0.0};
amrex::Real /* pert, */ rho, energy, T, pres;

auto eos = pele::physics::PhysicsType::eos();
if (sgn == -1) {
} else {
}
}

[baperry2]

See documentation on BCs here: https://amrex-combustion.github.io/PeleC/BoundaryConditions.html#boundary-conditions, NSCBC is not currently supported