Fuel initializartion example

Fuel example

Initialization

The following initialization block needs to be added for the following examples:

  /* If this is restarted computation, do not reinitialize values */

  if (domain->time_step->nt_prev > 0)
    return;

  const cs_real_t n_cells = domain->mesh->n_cells;

  const cs_combustion_model_t  *cm = cs_glob_combustion_model;
  const int io2 = cm->io2 - 1;
  const int in2 = cm->in2 - 1;
  const int ico2 = cm->ico2 - 1;
  const int ih2o = cm->ih2o - 1;

  const cs_real_t *wmole = cm->wmole;
  const cs_real_t *oxyo2 = cm->oxyo2;
  const cs_real_t *oxyn2 = cm->oxyn2;
  const cs_real_t *oxyh2o = cm->oxyh2o;
  const cs_real_t *oxyco2 = cm->oxyco2;

  /* The domain is filled with air at t1init;

     Computation of h1init */

  cs_real_t t1init = 1000.0;
  cs_real_t t2init = 1000.0;

  /* Transported variables for droplets */

  const cs_real_t trefth = 25.0 + cs_physical_constants_celsius_to_kelvin;

  cs_real_t h2init = cm->fuel.h02fol + cm->fuel.cp2fol*(t2init-trefth);

  /* Transported variables for the mix (droplets and carrying gases) */

  cs_real_t coefe[CS_COMBUSTION_GAS_MAX_GLOBAL_SPECIES];
  for (int ige = 0; ige < CS_COMBUSTION_GAS_MAX_GLOBAL_SPECIES; ige++)
    coefe[ige] = 0.;

  /* Consider oxydant 1 */

  coefe[io2] =    wmole[io2]*oxyo2[1]
               / (  wmole[io2]*oxyo2[1]  + wmole[in2]*oxyn2[1]
                  + wmole[ih2o]*oxyh2o[1] + wmole[ico2]*oxyco2[1]);
  coefe[ih2o] =   wmole[ih2o]*oxyh2o[1]
                / (  wmole[io2]*oxyo2[1] + wmole[in2] * oxyn2[1]
                   + wmole[ih2o]*oxyh2o[1] + wmole[ico2]*oxyco2[1]);
  coefe[ico2] =   wmole[ico2]*oxyco2[1]
                / (  wmole[io2]*oxyo2[1] +wmole[in2]*oxyn2[1]
                   + wmole[ih2o]*oxyh2o[1]+wmole[ico2]*oxyco2[1]);
  coefe[in2] = 1. - coefe[io2] - coefe[ih2o] - coefe[ico2];

  cs_real_t h1init = cs_fuel_t2h_gas(coefe, t1init);

  cs_array_real_set_scalar(n_cells, h1init, CS_F_(h)->val);

  /* Transported variables for the mix (passive scalars, variance)
   * Variables not present here are initialized to 0. */

  if (CS_F_(yco2) != NULL) {

    /* Consider oxydant 1 */

    const int ioxy = 0;
    const cs_real_t wmo2 = wmole[io2];
    const cs_real_t wmco2 = wmole[ico2];
    const cs_real_t wmh2o = wmole[ih2o];
    const cs_real_t wmn2 = wmole[in2];
    const cs_real_t dmas = ( oxyo2 [ioxy]*wmo2  + oxyn2[ioxy]*wmn2
                           + oxyh2o[ioxy]*wmh2o + oxyco2[ioxy]*wmco2);
    const cs_real_t co2 = oxyco2[ioxy]*wmco2/dmas;

    cs_array_real_set_scalar(n_cells, co2, CS_F_(yco2)->val);

  }