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

Include dependency graph for cs_soot_model.cpp:

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

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)$