Build an algebraic CDO face-based system for unsteady convection-diffusion-reaction of scalar-valued equations with source terms. 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_cdo_advection.h"
#include "cs_cdo_bc.h"
#include "cs_cdo_diffusion.h"
#include "cs_equation_bc.h"
#include "cs_equation_common.h"
#include "cs_evaluate.h"
#include "cs_hodge.h"
#include "cs_log.h"
#include "cs_math.h"
#include "cs_mesh_location.h"
#include "cs_post.h"
#include "cs_quadrature.h"
#include "cs_reco.h"
#include "cs_search.h"
#include "cs_sles.h"
#include "cs_source_term.h"
#include "cs_static_condensation.h"
#include "cs_cdofb_scaleq.h"

Include dependency graph for cs_cdofb_scaleq.c:

Functions
bool	cs_cdofb_scaleq_is_initialized (void)
	Check if the generic structures for building a CDO-Fb scheme are allocated. More...

void	cs_cdofb_scaleq_init_common (const cs_cdo_quantities_t quant, const cs_cdo_connect_t connect, const cs_time_step_t time_step, const cs_matrix_structure_t ms)
	Allocate work buffer and general structures related to CDO scalar-valued face-based schemes. Set shared pointers from the main domain members. More...

const cs_matrix_structure_t *	cs_cdofb_scaleq_matrix_structure (void)
	Get the pointer to the related cs_matrix_structure_t. More...

void	cs_cdofb_scaleq_get (cs_cell_sys_t csys, cs_cell_builder_t cb)
	Retrieve work buffers used for building a CDO system cellwise. More...

void	cs_cdofb_scaleq_finalize_common (void)
	Free work buffer and general structure related to CDO face-based schemes. More...

void *	cs_cdofb_scaleq_init_context (const cs_equation_param_t eqp, int var_id, int bflux_id, cs_equation_builder_t eqb)
	Initialize a cs_cdofb_scaleq_t structure storing data useful for building and managing such a scheme. More...

void *	cs_cdofb_scaleq_free_context (void *data)
	Destroy a cs_cdofb_scaleq_t structure. More...

void	cs_cdofb_scaleq_init_values (cs_real_t t_eval, const int field_id, const cs_mesh_t mesh, const cs_equation_param_t eqp, cs_equation_builder_t eqb, void context)
	Set the initial values of the variable field taking into account the boundary conditions. Case of scalar-valued CDO-Fb schemes. More...

void	cs_cdofb_scaleq_solve_steady_state (const cs_mesh_t mesh, const int field_id, const cs_equation_param_t eqp, cs_equation_builder_t eqb, void context)
	Build and solve the linear system arising from a scalar steady-state convection/diffusion/reaction equation with a CDO-Fb scheme One works cellwise and then process to the assembly. More...

void	cs_cdofb_scaleq_solve_implicit (const cs_mesh_t mesh, const int field_id, const cs_equation_param_t eqp, cs_equation_builder_t eqb, void context)
	Build and solve the linear system arising from a scalar convection/diffusion/reaction equation with a CDO-Fb scheme and an implicit Euler scheme. One works cellwise and then process to the assembly. More...

void	cs_cdofb_scaleq_solve_theta (const cs_mesh_t mesh, const int field_id, const cs_equation_param_t eqp, cs_equation_builder_t eqb, void context)
	Build and solve the linear system arising from a scalar convection/diffusion/reaction equation with a CDO-Fb scheme and an implicit/explicit theta scheme. One works cellwise and then process to the assembly. More...

cs_equation_balance_t *	cs_cdofb_scaleq_balance (const cs_equation_param_t eqp, cs_equation_builder_t eqb, void *context)
	Compute the balance for an equation over the full computational domain between time t_cur and t_cur + dt_cur Case of scalar-valued CDO face-based scheme. More...

void	cs_cdofb_scaleq_boundary_diff_flux (const cs_real_t t_eval, const cs_equation_param_t eqp, const cs_real_t pot_f, const cs_real_t pot_c, cs_equation_builder_t eqb, cs_real_t *bflux)
	Compute the diffusive flux across each boundary face. Case of scalar-valued CDO-Fb schemes. More...

void	cs_cdofb_scaleq_extra_op (const char eqname, const cs_field_t field, const cs_equation_param_t eqp, cs_equation_builder_t eqb, void *data)
	Predefined extra-operations related to this equation. More...

cs_real_t *	cs_cdofb_scaleq_get_cell_values (void *context)
	Get the computed values at mesh cells from the inverse operation w.r.t. the static condensation (DoF used in the linear system are located at primal faces) The lifecycle of this array is managed by the code. So one does not have to free the return pointer. More...

cs_real_t *	cs_cdofb_scaleq_get_face_values (void *context)
	Retrieve an array of values at mesh faces for the current context. The lifecycle of this array is managed by the code. So one does not have to free the return pointer. More...

void	cs_cdofb_scaleq_read_restart (cs_restart_t restart, const char eqname, void *scheme_context)
	Read additional arrays (not defined as fields) but useful for the checkpoint/restart process. More...

void	cs_cdofb_scaleq_write_restart (cs_restart_t restart, const char eqname, void *scheme_context)
	Write additional arrays (not defined as fields) but useful for the checkpoint/restart process. More...

Detailed Description

Build an algebraic CDO face-based system for unsteady convection-diffusion-reaction of scalar-valued equations with source terms.

Function Documentation

◆ cs_cdofb_scaleq_balance()

cs_equation_balance_t* cs_cdofb_scaleq_balance	(	const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		void *	context
	)

Compute the balance for an equation over the full computational domain between time t_cur and t_cur + dt_cur Case of scalar-valued CDO face-based scheme.

Parameters

[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in,out]	context	pointer to a scheme builder structure

Returns: a pointer to a cs_equation_balance_t structure

◆ cs_cdofb_scaleq_boundary_diff_flux()

void cs_cdofb_scaleq_boundary_diff_flux	(	const cs_real_t	t_eval,
		const cs_equation_param_t *	eqp,
		const cs_real_t *	pot_f,
		const cs_real_t *	pot_c,
		cs_equation_builder_t *	eqb,
		cs_real_t *	bflux
	)

Compute the diffusive flux across each boundary face. Case of scalar-valued CDO-Fb schemes.

Parameters

[in]	t_eval	time at which one performs the evaluation
[in]	eqp	pointer to a cs_equation_param_t structure
[in]	pot_f	array of values at faces
[in]	pot_c	array of values at cells
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in,out]	bflux	pointer to the values of the diffusive flux

◆ cs_cdofb_scaleq_extra_op()

void cs_cdofb_scaleq_extra_op	(	const char *	eqname,
		const cs_field_t *	field,
		const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		void *	data
	)

Predefined extra-operations related to this equation.

Parameters

[in]	eqname	name of the equation
[in]	field	pointer to a field structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in,out]	data	pointer to cs_cdofb_scaleq_t structure

◆ cs_cdofb_scaleq_finalize_common()

void cs_cdofb_scaleq_finalize_common ( void )

Free work buffer and general structure related to CDO face-based schemes.

◆ cs_cdofb_scaleq_free_context()

void* cs_cdofb_scaleq_free_context ( void * data )

Destroy a cs_cdofb_scaleq_t structure.

Parameters

[in,out] data pointer to a cs_cdofb_scaleq_t structure

Returns: a NULL pointer

◆ cs_cdofb_scaleq_get()

void cs_cdofb_scaleq_get	(	cs_cell_sys_t **	csys,
		cs_cell_builder_t **	cb
	)

Retrieve work buffers used for building a CDO system cellwise.

Parameters

[out]	csys	pointer to a pointer on a cs_cell_sys_t structure
[out]	cb	pointer to a pointer on a cs_cell_builder_t structure

◆ cs_cdofb_scaleq_get_cell_values()

cs_real_t* cs_cdofb_scaleq_get_cell_values ( void * context )

Get the computed values at mesh cells from the inverse operation w.r.t. the static condensation (DoF used in the linear system are located at primal faces) The lifecycle of this array is managed by the code. So one does not have to free the return pointer.

Parameters

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

Returns: a pointer to an array of cs_real_t

◆ cs_cdofb_scaleq_get_face_values()

cs_real_t* cs_cdofb_scaleq_get_face_values ( void * context )

Retrieve an array of values at mesh faces for the current context. The lifecycle of this array is managed by the code. So one does not have to free the return pointer.

Parameters

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

Returns: a pointer to an array of cs_real_t (size n_faces)

◆ cs_cdofb_scaleq_init_common()

void cs_cdofb_scaleq_init_common	(	const cs_cdo_quantities_t *	quant,
		const cs_cdo_connect_t *	connect,
		const cs_time_step_t *	time_step,
		const cs_matrix_structure_t *	ms
	)

Allocate work buffer and general structures related to CDO scalar-valued face-based schemes. Set shared pointers from the main domain members.

Parameters

[in]	quant	additional mesh quantities struct.
[in]	connect	pointer to a cs_cdo_connect_t struct.
[in]	time_step	pointer to a time step structure
[in]	ms	pointer to a cs_matrix_structure_t structure

◆ cs_cdofb_scaleq_init_context()

void* cs_cdofb_scaleq_init_context	(	const cs_equation_param_t *	eqp,
		int	var_id,
		int	bflux_id,
		cs_equation_builder_t *	eqb
	)

Initialize a cs_cdofb_scaleq_t structure storing data useful for building and managing such a scheme.

Parameters

[in]	eqp	pointer to a cs_equation_param_t structure
[in]	var_id	id of the variable field
[in]	bflux_id	id of the boundary flux field
[in,out]	eqb	pointer to a cs_equation_builder_t structure

Returns: a pointer to a new allocated cs_cdofb_scaleq_t structure

◆ cs_cdofb_scaleq_init_values()

void cs_cdofb_scaleq_init_values	(	cs_real_t	t_eval,
		const int	field_id,
		const cs_mesh_t *	mesh,
		const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		void *	context
	)

Set the initial values of the variable field taking into account the boundary conditions. Case of scalar-valued CDO-Fb schemes.

Parameters

[in]	t_eval	time at which one evaluates BCs
[in]	field_id	id related to the variable field of this equation
[in]	mesh	pointer to a cs_mesh_t structure
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in,out]	context	pointer to the scheme context (cast on-the-fly)

◆ cs_cdofb_scaleq_is_initialized()

bool cs_cdofb_scaleq_is_initialized ( void )

Check if the generic structures for building a CDO-Fb scheme are allocated.

Returns: true or false

◆ cs_cdofb_scaleq_matrix_structure()

const cs_matrix_structure_t* cs_cdofb_scaleq_matrix_structure ( void )

Get the pointer to the related cs_matrix_structure_t.

Returns: a pointer to a cs_matrix_structure_t structure

◆ cs_cdofb_scaleq_read_restart()

void cs_cdofb_scaleq_read_restart	(	cs_restart_t *	restart,
		const char *	eqname,
		void *	scheme_context
	)

Read additional arrays (not defined as fields) but useful for the checkpoint/restart process.

Parameters

[in,out]	restart	pointer to cs_restart_t structure
[in]	eqname	name of the related equation
[in]	scheme_context	pointer to a data structure cast on-the-fly

◆ cs_cdofb_scaleq_solve_implicit()

void cs_cdofb_scaleq_solve_implicit	(	const cs_mesh_t *	mesh,
		const int	field_id,
		const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		void *	context
	)

Build and solve the linear system arising from a scalar convection/diffusion/reaction equation with a CDO-Fb scheme and an implicit Euler scheme. One works cellwise and then process to the assembly.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	field_id	id of the variable field related to this equation
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in,out]	context	pointer to cs_cdofb_scaleq_t structure

◆ cs_cdofb_scaleq_solve_steady_state()

void cs_cdofb_scaleq_solve_steady_state	(	const cs_mesh_t *	mesh,
		const int	field_id,
		const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		void *	context
	)

Build and solve the linear system arising from a scalar steady-state convection/diffusion/reaction equation with a CDO-Fb scheme One works cellwise and then process to the assembly.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	field_id	id of the variable field related to this equation
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in,out]	context	pointer to cs_cdofb_scaleq_t structure

◆ cs_cdofb_scaleq_solve_theta()

void cs_cdofb_scaleq_solve_theta	(	const cs_mesh_t *	mesh,
		const int	field_id,
		const cs_equation_param_t *	eqp,
		cs_equation_builder_t *	eqb,
		void *	context
	)

Build and solve the linear system arising from a scalar convection/diffusion/reaction equation with a CDO-Fb scheme and an implicit/explicit theta scheme. One works cellwise and then process to the assembly.

Parameters

[in]	mesh	pointer to a cs_mesh_t structure
[in]	field_id	id of the variable field related to this equation
[in]	eqp	pointer to a cs_equation_param_t structure
[in,out]	eqb	pointer to a cs_equation_builder_t structure
[in,out]	context	pointer to cs_cdofb_scaleq_t structure

◆ cs_cdofb_scaleq_write_restart()

void cs_cdofb_scaleq_write_restart	(	cs_restart_t *	restart,
		const char *	eqname,
		void *	scheme_context
	)

Write additional arrays (not defined as fields) but useful for the checkpoint/restart process.

Parameters

[in,out]	restart	pointer to cs_restart_t structure
[in]	eqname	name of the related equation
[in]	scheme_context	pointer to a data structure cast on-the-fly

Functions

Detailed Description

Function Documentation

◆ cs_cdofb_scaleq_balance()

◆ cs_cdofb_scaleq_boundary_diff_flux()

◆ cs_cdofb_scaleq_extra_op()

◆ cs_cdofb_scaleq_finalize_common()

◆ cs_cdofb_scaleq_free_context()

◆ cs_cdofb_scaleq_get()

◆ cs_cdofb_scaleq_get_cell_values()

◆ cs_cdofb_scaleq_get_face_values()

◆ cs_cdofb_scaleq_init_common()

◆ cs_cdofb_scaleq_init_context()

◆ cs_cdofb_scaleq_init_values()

◆ cs_cdofb_scaleq_is_initialized()

◆ cs_cdofb_scaleq_matrix_structure()

◆ cs_cdofb_scaleq_read_restart()

◆ cs_cdofb_scaleq_solve_implicit()

◆ cs_cdofb_scaleq_solve_steady_state()

◆ cs_cdofb_scaleq_solve_theta()

◆ cs_cdofb_scaleq_write_restart()