Cell advection term

Functions
subroutine	mod_add_cell_advection_div_u_term::add_cell_advection_div_u_term (matrix, coefficient, divergence, stencil_size, equation_ls_map)
	Add to the matrix diagonal the second part of the cell advection term. More...

Detailed Description

Function Documentation

add_cell_advection_div_u_term()

subroutine mod_add_cell_advection_div_u_term::add_cell_advection_div_u_term	(	double precision, dimension(:), intent(inout)	matrix,
		double precision, dimension(:,:,:), intent(in)	coefficient,
		double precision, dimension(:,:,:), intent(in)	divergence,
		integer, intent(in)	stencil_size,
		type(t_ls_map), intent(in)	equation_ls_map
	)

The term writes \( \varrho S^{n+1} \nabla \cdot \mathbf{u^{n+1}} \). It is usually neglected since velocity field is solenoidal. But it may not be the case up to computer precision in some cases, so this term is kept.