Build an algebraic CDO face-based system for the Navier-Stokes equations and solved it with a monolithic approach. More...

#include "cs_defs.h"
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>
#include <assert.h>
#include <string.h>
#include <bft_mem.h>
#include "cs_blas.h"
#include "cs_cdo_bc.h"
#include "cs_cdofb_priv.h"
#include "cs_cdofb_scaleq.h"
#include "cs_cdofb_vecteq.h"
#include "cs_cdofb_navsto.h"
#include "cs_dbg.h"
#include "cs_equation_bc.h"
#include "cs_equation_common.h"
#include "cs_equation_priv.h"
#include "cs_evaluate.h"
#include "cs_fp_exception.h"
#include "cs_log.h"
#include "cs_math.h"
#include "cs_navsto_coupling.h"
#include "cs_param.h"
#include "cs_post.h"
#include "cs_sles.h"
#include "cs_source_term.h"
#include "cs_static_condensation.h"
#include "cs_timer.h"
#include "cs_cdofb_monolithic.h"

Include dependency graph for cs_cdofb_monolithic.c:

Functions
static void	_build_shared_structures (void)
	Define the cs_range_set_t, cs_interface_set_t, cs_matrix_assembler_t and cs_matrix_structure_t structures. More...

static void	_apply_bc_partly (const cs_cdofb_monolithic_t sc, const cs_equation_param_t eqp, const cs_cell_mesh_t cm, const cs_boundary_type_t bf_type, cs_cell_sys_t csys, cs_cell_builder_t cb)
	Apply the boundary conditions to the local system when this should be done before the static condensation. More...

static void	_apply_remaining_bc (const cs_cdofb_monolithic_t sc, const cs_equation_param_t eqp, const cs_cell_mesh_t cm, cs_cell_sys_t csys, cs_cell_builder_t cb, cs_cdofb_navsto_builder_t nsb, cs_real_t *mass_rhs)
	Apply the boundary conditions to the local system when this should be done after the static condensation. More...

static void	_assemble (const cs_cell_sys_t csys, const cs_cell_mesh_t cm, const cs_real_t div_op, const bool has_sourceterm, cs_matrix_assembler_values_t mav, cs_real_t rhs[], cs_real_t eqc_st[])
	Perform the assembly stage for a vector-valued system obtained with CDO-Fb schemes Shares similarities with cs_equation_assemble_block_matrix() More...

static void	_solve_system (cs_sles_t sles, const cs_matrix_t matrix, const cs_equation_param_t eqp, cs_real_t vel_f, cs_real_t pre_c, cs_real_t b_f, cs_real_t *b_c)
	Solve a linear system arising from a scalar-valued CDO-Fb scheme. More...

static void	_compute_and_update_fields (const cs_matrix_t matrix, cs_cdofb_monolithic_t sc, cs_equation_t mom_eq, cs_real_t mom_rhs, cs_real_t *mass_rhs)
	Solve the linear system and update the velocity fields and the pressure field. More...

void	cs_cdofb_monolithic_init_common (const cs_mesh_t mesh, const cs_cdo_quantities_t quant, const cs_cdo_connect_t connect, const cs_time_step_t time_step)
	Set shared pointers from the main domain members. More...

void *	cs_cdofb_monolithic_init_scheme_context (const cs_navsto_param_t nsp, cs_boundary_type_t bf_type, void *nsc_input)
	Initialize a cs_cdofb_monolithic_t structure. More...

void *	cs_cdofb_monolithic_free_scheme_context (void *scheme_context)
	Destroy a cs_cdofb_monolithic_t structure. More...

void	cs_cdofb_monolithic_set_sles (const cs_navsto_param_t nsp, void context)
	Start setting-up the Navier-Stokes equations when a monolithic algorithm is used to couple the system. No mesh information is available at this stage. More...

void	cs_cdofb_monolithic_compute_steady (const cs_mesh_t mesh, const cs_navsto_param_t nsp, void *scheme_context)
	Solve the steady Navier-Stokes system with a CDO face-based scheme using a monolithic approach. More...

void	cs_cdofb_monolithic_compute_implicit (const cs_mesh_t mesh, const cs_navsto_param_t nsp, void *scheme_context)
	Solve the unsteady Navier-Stokes system with a CDO face-based scheme using a monolithic approach and an Euler time scheme. More...

void	cs_cdofb_monolithic_compute_theta (const cs_mesh_t mesh, const cs_navsto_param_t nsp, void *scheme_context)
	Solve the unsteady Navier-Stokes system with a CDO face-based scheme using a monolithic approach and a theta time scheme. More...

Detailed Description

Build an algebraic CDO face-based system for the Navier-Stokes equations and solved it with a monolithic approach.

Function Documentation

◆ _apply_bc_partly()

static void _apply_bc_partly	(	const cs_cdofb_monolithic_t *	sc,
		const cs_equation_param_t *	eqp,
		const cs_cell_mesh_t *	cm,
		const cs_boundary_type_t *	bf_type,
		cs_cell_sys_t *	csys,
		cs_cell_builder_t *	cb
	)

static

Apply the boundary conditions to the local system when this should be done before the static condensation.

Parameters

[in]	sc	pointer to a cs_cdofb_monolithic_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in]	cm	pointer to a cellwise view of the mesh
[in]	bf_type	type of boundary for the boundary face
[in,out]	csys	pointer to a cellwise view of the system
[in,out]	cb	pointer to a cellwise builder

◆ _apply_remaining_bc()

static void _apply_remaining_bc	(	const cs_cdofb_monolithic_t *	sc,
		const cs_equation_param_t *	eqp,
		const cs_cell_mesh_t *	cm,
		cs_cell_sys_t *	csys,
		cs_cell_builder_t *	cb,
		cs_cdofb_navsto_builder_t *	nsb,
		cs_real_t *	mass_rhs
	)

static

Apply the boundary conditions to the local system when this should be done after the static condensation.

Parameters

[in]	sc	pointer to a cs_cdofb_monolithic_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in]	cm	pointer to a cellwise view of the mesh
[in,out]	csys	pointer to a cellwise view of the system
[in,out]	cb	pointer to a cellwise builder
[in,out]	nsb	builder structure for the NavSto system
[in,out]	mass_rhs	pointer to the rhs for the mass eq. in this cell

◆ _assemble()

static void _assemble	(	const cs_cell_sys_t *	csys,
		const cs_cell_mesh_t *	cm,
		const cs_real_t *	div_op,
		const bool	has_sourceterm,
		cs_matrix_assembler_values_t *	mav,
		cs_real_t	rhs[],
		cs_real_t	eqc_st[]
	)

static

Perform the assembly stage for a vector-valued system obtained with CDO-Fb schemes Shares similarities with cs_equation_assemble_block_matrix()

Parameters

[in]	csys	pointer to a cs_cell_sys_t structure
[in]	cm	pointer to a cs_cell_mesh_t structure
[in]	div_op	array with the divergence op. values
[in]	has_sourceterm	has the equation a source term?
[in,out]	mav	pointer to cs_matrix_assembler_values_t
[in,out]	rhs	right-end side of the system
[in,out]	eqc_st	source term from the context view

◆ _build_shared_structures()

static void _build_shared_structures ( void )

static

Define the cs_range_set_t, cs_interface_set_t, cs_matrix_assembler_t and cs_matrix_structure_t structures.

◆ _compute_and_update_fields()

static void _compute_and_update_fields	(	const cs_matrix_t *	matrix,
		cs_cdofb_monolithic_t *	sc,
		cs_equation_t *	mom_eq,
		cs_real_t *	mom_rhs,
		cs_real_t *	mass_rhs
	)

static

Solve the linear system and update the velocity fields and the pressure field.

Parameters

[in]	sc	scheme context
[in]	mom_eq	equation structure related to the momentum
[in]	matrix	pointer to a cs_matrix_t structure
[in,out]	mom_rhs	right-hand side for the momentum equation
[in,out]	mass_rhs	right_hand side for the mass equation

◆ _solve_system()

static void _solve_system	(	cs_sles_t *	sles,
		const cs_matrix_t *	matrix,
		const cs_equation_param_t *	eqp,
		cs_real_t *	vel_f,
		cs_real_t *	pre_c,
		cs_real_t *	b_f,
		cs_real_t *	b_c
	)

static

Solve a linear system arising from a scalar-valued CDO-Fb scheme.

Parameters

[in,out]	sles	pointer to a cs_sles_t structure
[in]	matrix	pointer to a cs_matrix_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	vel_f	initial velocity on faces
[in,out]	pre_c	initial pressure in cells
[in,out]	b_f	right-hand side (scatter/gather if needed) on faces
[in,out]	b_c	right_hand side on cells (mass equation)

◆ cs_cdofb_monolithic_compute_implicit()

void cs_cdofb_monolithic_compute_implicit	(	const cs_mesh_t *	mesh,
		const cs_navsto_param_t *	nsp,
		void *	scheme_context
	)

Solve the unsteady Navier-Stokes system with a CDO face-based scheme using a monolithic approach and an Euler time scheme.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	scheme_context	pointer to a structure cast on-the-fly

◆ cs_cdofb_monolithic_compute_steady()

void cs_cdofb_monolithic_compute_steady	(	const cs_mesh_t *	mesh,
		const cs_navsto_param_t *	nsp,
		void *	scheme_context
	)

Solve the steady Navier-Stokes system with a CDO face-based scheme using a monolithic approach.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	nsp	pointer to a cs_navsto_param_t structure
[in,out]	scheme_context	pointer to a structure cast on-the-fly

◆ cs_cdofb_monolithic_compute_theta()

void cs_cdofb_monolithic_compute_theta	(	const cs_mesh_t *	mesh,
		const cs_navsto_param_t *	nsp,
		void *	scheme_context
	)

Solve the unsteady Navier-Stokes system with a CDO face-based scheme using a monolithic approach and a theta time scheme.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	scheme_context	pointer to a structure cast on-the-fly

◆ cs_cdofb_monolithic_free_scheme_context()

void* cs_cdofb_monolithic_free_scheme_context ( void * scheme_context )

Destroy a cs_cdofb_monolithic_t structure.

Parameters

[in] scheme_context pointer to a scheme context structure to free

Returns: a NULL pointer

◆ cs_cdofb_monolithic_init_common()

void cs_cdofb_monolithic_init_common	(	const cs_mesh_t *	mesh,
		const cs_cdo_quantities_t *	quant,
		const cs_cdo_connect_t *	connect,
		const cs_time_step_t *	time_step
	)

Set shared pointers from the main domain members.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	quant	additional mesh quantities struct.
[in]	connect	pointer to a cs_cdo_connect_t struct.
[in]	time_step	pointer to a cs_time_step_t structure

◆ cs_cdofb_monolithic_init_scheme_context()

void* cs_cdofb_monolithic_init_scheme_context	(	const cs_navsto_param_t *	nsp,
		cs_boundary_type_t *	bf_type,
		void *	nsc_input
	)

Initialize a cs_cdofb_monolithic_t structure.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in]	bf_type	type of boundary for each boundary face
[in]	nsc_input	pointer to a cs_navsto_ac_t structure

Returns: a pointer to a new allocated cs_cdofb_monolithic_t structure

◆ cs_cdofb_monolithic_set_sles()

void cs_cdofb_monolithic_set_sles	(	const cs_navsto_param_t *	nsp,
		void *	context
	)

Start setting-up the Navier-Stokes equations when a monolithic algorithm is used to couple the system. No mesh information is available at this stage.

Start setting-up the Navier-Stokes equations when a monolithic algorithm is used to coupled the system. No mesh information is available at this stage.

Parameters

[in]	nsp	pointer to a cs_navsto_param_t structure
[in,out]	context	pointer to a context structure cast on-the-fly

Functions

Detailed Description

Function Documentation

◆ _apply_bc_partly()

◆ _apply_remaining_bc()

◆ _assemble()

◆ _build_shared_structures()

◆ _compute_and_update_fields()

◆ _solve_system()

◆ cs_cdofb_monolithic_compute_implicit()

◆ cs_cdofb_monolithic_compute_steady()

◆ cs_cdofb_monolithic_compute_theta()

◆ cs_cdofb_monolithic_free_scheme_context()

◆ cs_cdofb_monolithic_init_common()

◆ cs_cdofb_monolithic_init_scheme_context()

◆ cs_cdofb_monolithic_set_sles()