PeleC icon indicating copy to clipboard operation
PeleC copied to clipboard

Clarification Needed on LES Extrapolation Warning and SGS Model Implementation

Open JhonCordova opened this issue 9 months ago • 4 comments

Hello,

I am working with LES and I would like to know specifically what this message means: 'WARNING -- Need to implement LES extrap here, see diffusion,' which appears in the getLESTerm method defined in LES.cpp. Could you please provide me with more information?

void
PeleC::getLESTerm(
  amrex::Real time,
  amrex::Real dt,
  amrex::MultiFab& LESTerm,
  amrex::Real reflux_factor)
{
  BL_PROFILE("PeleC::getLESTerm()");

  if (!do_les) {
    LESTerm.setVal(0, 0, NVAR, LESTerm.nGrow());
    return;
  }

  if (verbose != 0) {
    amrex::Print() << "... Computing LES term at time " << time << std::endl;
  }

  switch (les_model) {

  case 0:
    getSmagorinskyLESTerm(time, dt, LESTerm, reflux_factor);
    break;

  case 1:
    getDynamicSmagorinskyLESTerm(time, dt, LESTerm, reflux_factor);
    break;

  default:
    amrex::Error("Invalid les_model number.");
    break;
  }

  amrex::Print()
    << "WARNING -- Need to implement LES extrap here, see diffusion"
    << std::endl;

  // Filter the SGS source term
  if (use_explicit_filter) {
    les_filter.apply_filter(LESTerm, filtered_les_source);
    LESTerm.define(
      grids, dmap, NVAR, filtered_les_source.nGrow(), amrex::MFInfo(),
      Factory());
    amrex::MultiFab::Copy(
      LESTerm, filtered_les_source, 0, 0, NVAR, filtered_les_source.nGrow());
  }
}

Additionally, I want to implement a one equation SGS model in LES. I have seen in the PeleC documentation, specifically in the equations section, that there are two advected quantities, Ak and Bk, that can be used. How can I enable their use in PeleC? How do I define them?

Thank you.

JhonCordova avatar May 26 '24 19:05 JhonCordova