Face nonlinear advection term discretization

Functions
subroutine	mod_add_u_div_u_term::add_u_div_u_term (matrix, density_face, divergence, navier_ls_map)
	Add ρ·u·div(u) to A.
subroutine	mod_compute_div_u_u_term_centered_o2_scheme::compute_div_u_u_term_centered_o2_scheme (matrix, density_face, theta, velocity, navier_stencil, navier_ls_map)
	Add ρ·div(uu) to A with a second-order centered scheme. More...
subroutine	mod_compute_div_u_u_term_theta_hybrid_centered_o2_upwind_scheme::compute_div_u_u_term_theta_hybrid_centered_o2_upwind_scheme (theta, velocity)
	Compute theta for hybrid centered/upwind scheme.
subroutine	mod_compute_div_u_u_term_upwind_o1_scheme::compute_div_u_u_term_upwind_o1_scheme (matrix, density_face, theta, velocity, navier_stencil, navier_ls_map)
	Add ρ·div(uu) to A with a first-order upwind scheme.
subroutine	mod_compute_div_u_u_term_upwind_o2_scheme::compute_div_u_u_term_upwind_o2_scheme (matrix, density_face, theta, velocity, navier_stencil, navier_ls_map)
	Add ρ·div(uu) to A with a second-order upwind scheme.

Detailed Description

Function Documentation

compute_div_u_u_term_centered_o2_scheme()

subroutine mod_compute_div_u_u_term_centered_o2_scheme::compute_div_u_u_term_centered_o2_scheme	(	double precision, dimension(:), intent(inout)	matrix,
		type(t_face_field), intent(in)	density_face,
		type(t_face_field), intent(in)	theta,
		type(t_face_field), intent(in)	velocity,
		type(t_face_stencil), intent(in)	navier_stencil,
		type(t_face_ls_map), intent(in)	navier_ls_map
	)

The nonlinear momentum term \( \rho * \nabla \cdot (u_{n+1} \otimes u_{n+1}) \) with a second-order centered scheme.