#include "cs_defs.h"
#include <assert.h>
#include <stdarg.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include "bft_mem.h"
#include "bft_error.h"
#include "bft_printf.h"
#include "cs_array.h"
#include "cs_combustion_model.h"
#include "cs_cogz_headers.h"
#include "cs_field.h"
#include "cs_field_default.h"
#include "cs_field_pointer.h"
#include "cs_log.h"
#include "cs_math.h"
#include "cs_mesh.h"
#include "cs_mesh_quantities.h"
#include "cs_parall.h"
#include "cs_soot_model.h"

Include dependency graph for cs_soot_model.c:

Functions
void	cs_soot_production (int f_id, cs_real_t smbrs[], cs_real_t rovsdt[])
	Specific physical models: soot production model. More...

Detailed Description

Soot related models.

Function Documentation

◆ cs_soot_production()

void cs_soot_production	(	int	f_id,
		cs_real_t	smbrs[],
		cs_real_t	rovsdt[]
	)

Specific physical models: soot production model.

Specific physic subroutine : soot production model.

Define the source terms for the soot mass fraction and the precursor number for soot model. The equations read: $rovsdt \delta a = smbrs$

et could already contain source term and don't have to be erased but incremented.

For stability sake, only positive terms should be add in . There is no constrain for . 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 are set and containning $\rho \Omega$

in $kg_a.s^{-1}$ (ex: for velocity: $kg.m.s^{-2}$ , for temperature: $kg.°C.s^{-1}$ , for enthalpy: $J.s^{-1}$ )
en $kg.s^{-1}$

Parameters

[in]	f_id	scalar field id
[out]	smbrs	explicit right hand side
[out]	rovsdt	implicit terms

Functions

Detailed Description

Function Documentation

◆ cs_soot_production()