9.1
general documentation
Functions
cs_combustion_lw.h File Reference
#include "base/cs_defs.h"
#include "base/cs_field.h"
#include "pprt/cs_physical_model.h"
+ Include dependency graph for cs_combustion_lw.h:

Go to the source code of this file.

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...
 

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