9.1
general documentation
Functions
cs_combustion_lw.cpp File Reference
#include "base/cs_defs.h"
#include <algorithm>
#include <assert.h>
#include <stdio.h>
#include <stdarg.h>
#include <stdlib.h>
#include <string.h>
#include "bft/bft_error.h"
#include "bft/bft_printf.h"
#include "base/cs_array.h"
#include "base/cs_array_reduce.h"
#include "base/cs_base.h"
#include "base/cs_dispatch.h"
#include "cdo/cs_equation_param.h"
#include "base/cs_field.h"
#include "base/cs_field_default.h"
#include "base/cs_field_operator.h"
#include "base/cs_field_pointer.h"
#include "base/cs_physical_constants.h"
#include "base/cs_log.h"
#include "base/cs_math.h"
#include "base/cs_parall.h"
#include "base/cs_physical_properties.h"
#include "base/cs_restart.h"
#include "base/cs_restart_default.h"
#include "mesh/cs_mesh.h"
#include "rayt/cs_rad_transfer.h"
#include "turb/cs_turbulence_model.h"
#include "pprt/cs_combustion_model.h"
#include "pprt/cs_physical_model.h"
#include "cogz/cs_combustion_gas.h"
#include "cogz/cs_combustion_boundary_conditions.h"
#include "cogz/cs_combustion_ht_convert.h"
#include "cogz/cs_combustion_lw.h"
+ Include dependency graph for cs_combustion_lw.cpp:

Functions

void cs_combustion_lw_fields_init (void)
 Initialize specific fields for Libby-Williams gas combustion model. More...
 
void cs_combustion_lw_physical_prop (void)
 
void cs_combustion_lw_source_terms (cs_field_t *f_sc, cs_real_t smbrs[], cs_real_t rovsdt[])
 Compute source terms for premixed flame Libby-Williams combustion model. More...
 

Detailed Description

Libby-Williams gas combustion model.

Function Documentation

◆ cs_combustion_lw_fields_init()

void cs_combustion_lw_fields_init ( void  )

Initialize specific fields for Libby-Williams gas combustion model.

◆ cs_combustion_lw_physical_prop()

void cs_combustion_lw_physical_prop ( void  )

◆ cs_combustion_lw_source_terms()

void cs_combustion_lw_source_terms ( cs_field_t *  f_sc,
cs_real_t  smbrs[],
cs_real_t  rovsdt[] 
)

Compute source terms for premixed flame Libby-Williams combustion model.

Define the source terms for a given scalar over one time step.

The equations read: $ rovsdt \delta a = smbrs $

$ rovsdt $ et $ smbrs $ may already contain source term so must not be overwritten, but incremented.

For stability, only positive terms should be add in $ rovsdt $. There is no constraint for $ smbrs $. For a source term written $ S_{exp} + S_{imp} a $, source terms are: $ smbrs = smbrs + S_{exp} + S_{imp} a $ $ rovsdt = rovsdt + \max(-S_{imp},0) $

Here we set $ rovsdt $ and $ smbrs $ containing $ \rho \Omega $

  • $ smbrs $ in $ kg_a.s^{-1} $ (ex: for velocity: $ kg.m.s^{-2} $, for temperature: $ kg.C.s^{-1} $, for enthalpy: $ J.s^{-1} $)
  • $ rovsdt $ in $ kg.s^{-1} $
Parameters
[in]f_scpointer to scalar field
[in,out]smbrsexplicit right hand side
[in,out]rovsdtimplicit terms