Pulverized coal

The examples provided here are only concern inlet boundary conditions. Other boundary condition types are not specific to the pulverized coal model, so generic examples apply (and as usual, the GUI is recommended).

Note that turbulence variables BC's are usually computed automatically, based on definitions using cs_boundary_conditions_inlet_set_turbulence_hyd_diam or cs_boundary_conditions_inlet_set_turbulence_internsity.

This following examples show how user code con be defined in cs_user_boundary_conditions_setup, so as to define inlet properties using the high-level for non-coal-specific variables, and access the coal inlet structure definition to set the required properties.

Example 1

Secondary or tertiary air inlet, for zone 'inlet', using either the base pulverized coal model or the pulverized coal with moisture model (with Lagrangian particle tracking).

  {
    const cs_zone_t *z = cs_boundary_zone_by_name("inlet");
 
    cs_coal_bc_inlet_t *ci = cs_coal_boundary_conditions_get_inlet(z);
 
    /* Inlet velocity direction based on constant vector;
       value based on oxidizer mass flow rate (in kg/s);*/
 
    cs_real_t u[] = {0, 0, 0.5};
    cs_boundary_conditions_open_set_velocity_by_value(z, u);
 
    cs_boundary_conditions_open_set_mass_flow_rate_by_value(z, 1.46e-03);
 
    ci->ientat = 1;
 
    /* Oxidizer's number (1 to 3) */
    ci->inmoxy = 1;
 
    /* Oxidizer Temperature in K */
    ci->t_air = 400.0 + cs_physical_constants_celsius_to_kelvin;
 
    /* Turbulence variables BC's calculated based on hydraulic diameter */
 
    cs_boundary_conditions_inlet_set_turbulence_hyd_diam(z, 0.032);
 
    /* Automatic treatment of non-user scalars */
 
    /* Treatment of user-defined scalars */
 
    const char *names[] = {"scalar1", "scalar2"};
 
    for (int i = 0; i < 2; i++) {
 
      cs_equation_param_t  *eqp = cs_equation_param_by_name(names[i]);
      if (eqp == nullptr)
        continue;
 
      cs_real_t val[1] = {1.0};
      cs_equation_add_bc_by_value(eqp,
                                  CS_BC_DIRICHLET,
                                  z->name,
                                  val);
 
    }
  }

Example 2

Primary air and coal inlet, for zone 'primary_inlet'.

  {
    const cs_zone_t *z = cs_boundary_zone_by_name("primary_inlet");
 
    cs_coal_bc_inlet_t *ci = cs_coal_boundary_conditions_get_inlet(z);
 
    /* Inlet velocity direction based on constant vector;
       value based on oxidizer mass flow rate (in kg/s);*/
 
    cs_real_t u[] = {0, 0, 0.5};
    cs_boundary_conditions_open_set_velocity_by_value(z, u);
 
    cs_coal_boundary_conditions_inlet_set_air_mass_flow_rate_by_value
      (z, 1.46e-03);
 
    /* Turbulence variables BC's calculated based on hydraulic diameter */
 
    cs_boundary_conditions_inlet_set_turbulence_hyd_diam(z, 0.1);
 
    /* Air inlet with pulverized coal */
    ci->ientcp = 1;
 
    /* Oxidizer's number (1 to 3) */
    ci->inmoxy = 1;
 
    /* Oxidizer Temperature in K */
    ci->t_air = 800.0 + cs_physical_constants_celsius_to_kelvin;
 
    /* code_saturne deals with ncha different coals (component of blend)
       every coal is described by n_classes_per_coal[icha] class of particles
       (each of them described by an inlet diameter) */
 
    cs_coal_model_t *coal = cs_glob_coal_model;
 
    /* Treatment for the first coal */
 
    int icha = 0;
 
    /* Coal mass flow rate in kg/s */
    ci->qimpcp[icha] = 1.46e-4;
 
    /* Percentage mass fraction of each granulometric class */
    for (int iclapc = 0; iclapc < coal->n_classes_per_coal[icha]; iclapc++) {
      ci->distch[icha][iclapc] = 100.0/(double)coal->n_classes_per_coal[icha];
    }
 
    /* Inlet temperature for coal & primary air */
    ci->timpcp[icha] = 800.0 + cs_physical_constants_celsius_to_kelvin;
 
    /* Automatic treatment of non-user scalars */
  }