9.2
general documentation
cs_cdofb_navsto.h
Go to the documentation of this file.
1#ifndef CS_CDOFB_NAVSTO_H
2#define CS_CDOFB_NAVSTO_H
3
4/*============================================================================
5 * Functions shared among all face-based schemes for the discretization of the
6 * Navier--Stokes system
7 *============================================================================*/
8
9/*
10 This file is part of code_saturne, a general-purpose CFD tool.
11
12 Copyright (C) 1998-2026 EDF S.A.
13
14 This program is free software; you can redistribute it and/or modify it under
15 the terms of the GNU General Public License as published by the Free Software
16 Foundation; either version 2 of the License, or (at your option) any later
17 version.
18
19 This program is distributed in the hope that it will be useful, but WITHOUT
20 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
21 FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
22 details.
23
24 You should have received a copy of the GNU General Public License along with
25 this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
26 Street, Fifth Floor, Boston, MA 02110-1301, USA.
27*/
28
29/*----------------------------------------------------------------------------*/
30
31#include "base/cs_defs.h"
32
33/*----------------------------------------------------------------------------
34 * Standard C library headers
35 *----------------------------------------------------------------------------*/
36
37/*----------------------------------------------------------------------------
38 * Local headers
39 *----------------------------------------------------------------------------*/
40
41#include "alge/cs_matrix.h"
42#include "base/cs_base.h"
43#include "base/cs_field.h"
44#include "base/cs_math.h"
45#include "base/cs_time_plot.h"
46#include "base/cs_time_step.h"
47#include "cdo/cs_cdo_connect.h"
48#include "cdo/cs_cdo_quantities.h"
49#include "cdo/cs_cdo_turbulence.h"
50#include "cdo/cs_iter_algo.h"
51#include "cdo/cs_navsto_context.h"
52#include "cdo/cs_navsto_param.h"
53#include "cdo/cs_sdm.h"
54#include "mesh/cs_mesh.h"
55
56/*============================================================================
57 * Macro definitions
58 *============================================================================*/
59
60/*============================================================================
61 * Type definitions
62 *============================================================================*/
63
69typedef enum {
70
80 CS_CDOFB_NAVSTO_BOUSSINESQ_CELL_DOF,
81
91 CS_CDOFB_NAVSTO_BOUSSINESQ_FACE_DOF
92
93} cs_cdofb_navsto_boussinesq_type_t;
94
95
105typedef struct {
106
115 cs_real_t rho_c;
116
128 cs_real_t *div_op;
129
139 cs_real_t mass_rhs;
140
159 cs_boundary_type_t *bf_type; /* Size: n_fc */
160 cs_real_t *pressure_bc_val; /* Size: n_fc */
161
166} cs_cdofb_navsto_builder_t;
167
168/*----------------------------------------------------------------------------*/
181/*----------------------------------------------------------------------------*/
182
183typedef void
184(cs_cdofb_navsto_source_t)(const cs_navsto_param_t *nsp,
185 const cs_cell_mesh_t *cm,
186 const cs_cdofb_navsto_builder_t *nsb,
187 cs_cell_sys_t *csys);
188
189/*============================================================================
190 * Public function prototypes
191 *============================================================================*/
192
193/*----------------------------------------------------------------------------*/
199/*----------------------------------------------------------------------------*/
200
201void
202cs_cdofb_navsto_set_boussinesq_algo(cs_cdofb_navsto_boussinesq_type_t type);
203
204/*----------------------------------------------------------------------------*/
213/*----------------------------------------------------------------------------*/
214
215cs_cdofb_navsto_builder_t
216cs_cdofb_navsto_create_builder(const cs_navsto_param_t *nsp,
217 const cs_cdo_connect_t *connect);
218
219/*----------------------------------------------------------------------------*/
225/*----------------------------------------------------------------------------*/
226
227void
228cs_cdofb_navsto_free_builder(cs_cdofb_navsto_builder_t *nsb);
229
230/*----------------------------------------------------------------------------*/
242/*----------------------------------------------------------------------------*/
243
244void
245cs_cdofb_navsto_define_builder(cs_real_t t_eval,
246 const cs_navsto_param_t *nsp,
247 const cs_cell_mesh_t *cm,
248 const cs_cell_sys_t *csys,
249 const cs_cdo_bc_face_t *pr_bc,
250 const cs_boundary_type_t *bf_type,
251 cs_cdofb_navsto_builder_t *nsb);
252
253/*----------------------------------------------------------------------------*/
265/*----------------------------------------------------------------------------*/
266
267void
268cs_cdofb_navsto_mass_flux(const cs_navsto_param_t *nsp,
269 const cs_cdo_quantities_t *quant,
270 const cs_real_t *face_vel,
271 cs_real_t *mass_flux);
272
273/*----------------------------------------------------------------------------*/
286/*----------------------------------------------------------------------------*/
287
288double
289cs_cdofb_navsto_cell_divergence(const cs_lnum_t c_id,
290 const cs_cdo_quantities_t *quant,
291 const cs_adjacency_t *c2f,
292 const cs_real_t *f_vals);
293
294/*----------------------------------------------------------------------------*/
306/*----------------------------------------------------------------------------*/
307
308void
309cs_cdofb_navsto_compute_face_pressure(const cs_mesh_t *mesh,
310 const cs_cdo_connect_t *connect,
311 const cs_cdo_quantities_t *quant,
312 const cs_time_step_t *ts,
313 const cs_navsto_param_t *nsp,
314 const cs_real_t *p_cell,
315 cs_real_t *p_face);
316
317/*----------------------------------------------------------------------------*/
327/*----------------------------------------------------------------------------*/
328
329void
330cs_cdofb_navsto_add_grad_div(short int n_fc,
331 const cs_real_t zeta,
332 const cs_real_t div[],
333 cs_sdm_t *mat);
334
335/*----------------------------------------------------------------------------*/
344/*----------------------------------------------------------------------------*/
345
346void
347cs_cdofb_navsto_init_pressure(const cs_navsto_param_t *nsp,
348 const cs_cdo_quantities_t *quant,
349 const cs_time_step_t *ts,
350 cs_field_t *pr);
351
352/*----------------------------------------------------------------------------*/
364/*----------------------------------------------------------------------------*/
365
366void
367cs_cdofb_navsto_check_init(const cs_navsto_param_t *nsp,
368 const cs_cdo_connect_t *connect,
369 const cs_cdo_quantities_t *quant,
370 const cs_time_step_t *ts,
371 const cs_field_t *velocity,
372 const cs_real_t *face_vel,
373 const cs_field_t *pressure);
374
375/*----------------------------------------------------------------------------*/
385/*----------------------------------------------------------------------------*/
386
387void
388cs_cdofb_navsto_init_face_pressure(const cs_navsto_param_t *nsp,
389 const cs_cdo_connect_t *connect,
390 const cs_time_step_t *ts,
391 cs_real_t *pr_f);
392
393/*----------------------------------------------------------------------------*/
402/*----------------------------------------------------------------------------*/
403
404void
405cs_cdofb_navsto_rescale_pressure_to_ref(const cs_navsto_param_t *nsp,
406 const cs_cdo_quantities_t *quant,
407 cs_real_t values[]);
408
409/*----------------------------------------------------------------------------*/
417/*----------------------------------------------------------------------------*/
418
419void
420cs_cdofb_navsto_set_zero_mean_pressure(const cs_cdo_quantities_t *quant,
421 cs_real_t values[]);
422
423/*----------------------------------------------------------------------------*/
455/*----------------------------------------------------------------------------*/
456
457void
458cs_cdofb_navsto_extra_op(const cs_navsto_param_t *nsp,
459 const cs_mesh_t *mesh,
460 const cs_cdo_quantities_t *quant,
461 const cs_cdo_connect_t *connect,
462 const cs_time_step_t *ts,
463 cs_time_plot_t *time_plotter,
464 const cs_turbulence_t *turb,
465 const cs_adv_field_t *adv_field,
466 const cs_real_t *mass_flux,
467 const cs_real_t *p_cell,
468 const cs_real_t *u_cell,
469 const cs_real_t *u_face,
470 const cs_cdo_navsto_psteady_cvg_t &ps_cvg);
471
472/*----------------------------------------------------------------------------*/
487/*----------------------------------------------------------------------------*/
488
489void
490cs_cdofb_block_dirichlet_alge(short int f,
491 const cs_equation_param_t *eqp,
492 const cs_cell_mesh_t *cm,
493 const cs_property_data_t *pty,
494 cs_cell_builder_t *cb,
495 cs_cell_sys_t *csys);
496
497/*----------------------------------------------------------------------------*/
514/*----------------------------------------------------------------------------*/
515
516void
517cs_cdofb_block_dirichlet_pena(short int f,
518 const cs_equation_param_t *eqp,
519 const cs_cell_mesh_t *cm,
520 const cs_property_data_t *pty,
521 cs_cell_builder_t *cb,
522 cs_cell_sys_t *csys);
523
524/*----------------------------------------------------------------------------*/
541/*----------------------------------------------------------------------------*/
542
543void
544cs_cdofb_block_dirichlet_weak(short int fb,
545 const cs_equation_param_t *eqp,
546 const cs_cell_mesh_t *cm,
547 const cs_property_data_t *pty,
548 cs_cell_builder_t *cb,
549 cs_cell_sys_t *csys);
550
551/*----------------------------------------------------------------------------*/
568/*----------------------------------------------------------------------------*/
569
570void
571cs_cdofb_block_dirichlet_wsym(short int fb,
572 const cs_equation_param_t *eqp,
573 const cs_cell_mesh_t *cm,
574 const cs_property_data_t *pty,
575 cs_cell_builder_t *cb,
576 cs_cell_sys_t *csys);
577
578/*----------------------------------------------------------------------------*/
594/*----------------------------------------------------------------------------*/
595
596void
597cs_cdofb_symmetry(short int fb,
598 const cs_equation_param_t *eqp,
599 const cs_cell_mesh_t *cm,
600 const cs_property_data_t *pty,
601 cs_cell_builder_t *cb,
602 cs_cell_sys_t *csys);
603
604/*----------------------------------------------------------------------------*/
617/*----------------------------------------------------------------------------*/
618
619void
620cs_cdofb_prescribed_smooth_wall_n_pena_t_robin(
621 short int fb,
622 const cs_equation_param_t *eqp,
623 const cs_cell_mesh_t *cm,
624 const cs_property_data_t *pty,
625 cs_cell_builder_t *cb,
626 cs_cell_sys_t *csys);
627
628/*----------------------------------------------------------------------------*/
641/*----------------------------------------------------------------------------*/
642
643void
644cs_cdofb_prescribed_smooth_wall_n_alge_t_robin(
645 short int fb,
646 const cs_equation_param_t *eqp,
647 const cs_cell_mesh_t *cm,
648 const cs_property_data_t *pty,
649 cs_cell_builder_t *cb,
650 cs_cell_sys_t *csys);
651
652/*----------------------------------------------------------------------------*/
665/*----------------------------------------------------------------------------*/
666
667void
668cs_cdofb_prescribed_smooth_wall_n_alge_t_neumann (
669 short int fb,
670 const cs_equation_param_t *eqp,
671 const cs_cell_mesh_t *cm,
672 const cs_property_data_t *pty,
673 cs_cell_builder_t *cb,
674 cs_cell_sys_t *csys);
675
676/*----------------------------------------------------------------------------*/
689/*----------------------------------------------------------------------------*/
690
691cs_sles_convergence_state_t
692cs_cdofb_navsto_nl_algo_cvg(const cs_navsto_param_t *nsp,
693 const cs_real_t *pre_iterate,
694 cs_real_t *cur_iterate,
695 cs_iter_algo_t *algo);
696
697/*----------------------------------------------------------------------------*/
706/*----------------------------------------------------------------------------*/
707
708void
709cs_cdofb_navsto_set_gravity_func(const cs_navsto_param_t *nsp,
710 cs_cdofb_navsto_source_t **p_func);
711
712/*----------------------------------------------------------------------------*/
724/*----------------------------------------------------------------------------*/
725
726void
727cs_cdofb_navsto_gravity_term(const cs_navsto_param_t *nsp,
728 const cs_cell_mesh_t *cm,
729 const cs_cdofb_navsto_builder_t *nsb,
730 cs_cell_sys_t *csys);
731
732/*----------------------------------------------------------------------------*/
745/*----------------------------------------------------------------------------*/
746
747void
748cs_cdofb_navsto_boussinesq_at_cell(const cs_navsto_param_t *nsp,
749 const cs_cell_mesh_t *cm,
750 const cs_cdofb_navsto_builder_t *nsb,
751 cs_cell_sys_t *csys);
752
753/*----------------------------------------------------------------------------*/
766/*----------------------------------------------------------------------------*/
767
768void
769cs_cdofb_navsto_boussinesq_at_face(const cs_navsto_param_t *nsp,
770 const cs_cell_mesh_t *cm,
771 const cs_cdofb_navsto_builder_t *nsb,
772 cs_cell_sys_t *csys);
773
774/*----------------------------------------------------------------------------*/
787/*----------------------------------------------------------------------------*/
788
789void
790cs_cdofb_navsto_boussinesq_by_part(const cs_navsto_param_t *nsp,
791 const cs_cell_mesh_t *cm,
792 const cs_cdofb_navsto_builder_t *nsb,
793 cs_cell_sys_t *csys);
794
795/*----------------------------------------------------------------------------*/
807/*----------------------------------------------------------------------------*/
808
809void
810cs_cdofb_navsto_stream_source_term(cs_lnum_t n_elts,
811 const cs_lnum_t *elt_ids,
812 bool dense_output,
813 void *input,
814 cs_real_t *retval);
815
816/*----------------------------------------------------------------------------*/
832/*----------------------------------------------------------------------------*/
833
834cs_cdo_balance_t *
835cs_cdofb_navsto_balance(const cs_navsto_param_t *nsp,
836 const cs_cdo_connect_t *connect,
837 const cs_cdo_quantities_t *quant,
838 const cs_time_step_t *ts,
839 const cs_cdo_bc_face_t *pr_bc,
840 const cs_boundary_type_t *bf_type,
841 const cs_real_t *pr_c);
842
843/*----------------------------------------------------------------------------*/
857/*----------------------------------------------------------------------------*/
858
859void
860cs_cdofb_navsto_check_convergence(const cs_navsto_param_t *nsp,
861 const cs_cdo_quantities_t *quant,
862 const cs_time_step_t *ts,
863 const cs::cdo_navsto_ctx_t *sc,
864 const cs_turbulence_t *tbs,
865 cs_cdo_navsto_psteady_cvg_t &ps_cvg);
866
867#endif /* CS_CDOFB_NAVSTO_H */
cs_field_t
Field descriptor.
Definition: cs_field.h:275
cs_time_step_t
time step descriptor
Definition: cs_time_step.h:60
cs_base.h
cs_boundary_type_t
int cs_boundary_type_t
Definition: cs_boundary.h:65
cs_cdo_connect.h
cs_cdo_quantities.h
cs_cdo_turbulence.h
cs_cdofb_navsto_gravity_term
void cs_cdofb_navsto_gravity_term(const cs_navsto_param_t *nsp, const cs_cell_mesh_t *cm, const cs_cdofb_navsto_builder_t *nsb, cs_cell_sys_t *csys)
Take into account the gravity effects. Compute and add the source term to the local RHS....
Definition: cs_cdofb_navsto.cpp:2609
cs_cdofb_navsto_nl_algo_cvg
cs_sles_convergence_state_t cs_cdofb_navsto_nl_algo_cvg(const cs_navsto_param_t *nsp, const cs_real_t *pre_iterate, cs_real_t *cur_iterate, cs_iter_algo_t *algo)
Test if one has to do one more non-linear iteration. Test if performed on the relative norm on the in...
Definition: cs_cdofb_navsto.cpp:2514
cs_cdofb_block_dirichlet_pena
void cs_cdofb_block_dirichlet_pena(short int f, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a Dirichlet BCs on the three velocity components. For instance for a velocity inlet...
Definition: cs_cdofb_navsto.cpp:1665
cs_cdofb_prescribed_smooth_wall_n_pena_t_robin
void cs_cdofb_prescribed_smooth_wall_n_pena_t_robin(short int fb, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a wall BCs by Robin mix BC on vectors'tangential component, and strong penalization...
Definition: cs_cdofb_navsto.cpp:2042
cs_cdofb_navsto_add_grad_div
void cs_cdofb_navsto_add_grad_div(short int n_fc, const cs_real_t zeta, const cs_real_t div[], cs_sdm_t *mat)
Add the grad-div part to the local matrix (i.e. for the current cell)
Definition: cs_cdofb_navsto.cpp:618
cs_cdofb_navsto_boussinesq_at_face
void cs_cdofb_navsto_boussinesq_at_face(const cs_navsto_param_t *nsp, const cs_cell_mesh_t *cm, const cs_cdofb_navsto_builder_t *nsb, cs_cell_sys_t *csys)
Take into account the buoyancy force with the Boussinesq approx. Compute and add the source term to t...
Definition: cs_cdofb_navsto.cpp:2710
cs_cdofb_navsto_compute_face_pressure
void cs_cdofb_navsto_compute_face_pressure(const cs_mesh_t *mesh, const cs_cdo_connect_t *connect, const cs_cdo_quantities_t *quant, const cs_time_step_t *ts, const cs_navsto_param_t *nsp, const cs_real_t *p_cell, cs_real_t *p_face)
Compute an estimation of the pressure at faces.
Definition: cs_cdofb_navsto.cpp:546
cs_cdofb_block_dirichlet_wsym
void cs_cdofb_block_dirichlet_wsym(short int fb, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a Dirichlet BCs on the three velocity components. For instance for a velocity inlet...
Definition: cs_cdofb_navsto.cpp:1820
cs_cdofb_symmetry
void cs_cdofb_symmetry(short int fb, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a boundary defined as 'symmetry' (treated as a sliding BCs on the three velocity co...
Definition: cs_cdofb_navsto.cpp:1927
cs_cdofb_navsto_init_pressure
void cs_cdofb_navsto_init_pressure(const cs_navsto_param_t *nsp, const cs_cdo_quantities_t *quant, const cs_time_step_t *ts, cs_field_t *pr)
Initialize the pressure values.
Definition: cs_cdofb_navsto.cpp:691
cs_cdofb_block_dirichlet_alge
void cs_cdofb_block_dirichlet_alge(short int f, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a Dirichlet BCs on the three velocity components. For instance for a velocity inlet...
Definition: cs_cdofb_navsto.cpp:1563
cs_cdofb_navsto_check_init
void cs_cdofb_navsto_check_init(const cs_navsto_param_t *nsp, const cs_cdo_connect_t *connect, const cs_cdo_quantities_t *quant, const cs_time_step_t *ts, const cs_field_t *velocity, const cs_real_t *face_vel, const cs_field_t *pressure)
Check initial state.
Definition: cs_cdofb_navsto.cpp:894
cs_cdofb_navsto_boussinesq_type_t
cs_cdofb_navsto_boussinesq_type_t
Type of algorithm to compute the Boussinesq approximation.
Definition: cs_cdofb_navsto.h:69
CS_CDOFB_NAVSTO_BOUSSINESQ_CELL_DOF
@ CS_CDOFB_NAVSTO_BOUSSINESQ_CELL_DOF
Boussinesq approximation relyong on a cell contribution.
Definition: cs_cdofb_navsto.h:80
CS_CDOFB_NAVSTO_BOUSSINESQ_FACE_DOF
@ CS_CDOFB_NAVSTO_BOUSSINESQ_FACE_DOF
Boussinesq approximation relyong on face contributions.
Definition: cs_cdofb_navsto.h:91
cs_cdofb_navsto_rescale_pressure_to_ref
void cs_cdofb_navsto_rescale_pressure_to_ref(const cs_navsto_param_t *nsp, const cs_cdo_quantities_t *quant, cs_real_t values[])
Update the pressure field in order to get a field with a mean-value equal to the reference value.
Definition: cs_cdofb_navsto.cpp:947
cs_cdofb_navsto_free_builder
void cs_cdofb_navsto_free_builder(cs_cdofb_navsto_builder_t *nsb)
Destroy the given cs_cdofb_navsto_builder_t structure.
Definition: cs_cdofb_navsto.cpp:239
cs_cdofb_navsto_cell_divergence
double cs_cdofb_navsto_cell_divergence(const cs_lnum_t c_id, const cs_cdo_quantities_t *quant, const cs_adjacency_t *c2f, const cs_real_t *f_vals)
Compute the divergence in a cell of a vector-valued array defined at faces (values are defined both a...
Definition: cs_cdofb_navsto.cpp:487
cs_cdofb_navsto_set_boussinesq_algo
void cs_cdofb_navsto_set_boussinesq_algo(cs_cdofb_navsto_boussinesq_type_t type)
Set the way to compute the Boussinesq approximation.
Definition: cs_cdofb_navsto.cpp:191
cs_cdofb_navsto_boussinesq_at_cell
void cs_cdofb_navsto_boussinesq_at_cell(const cs_navsto_param_t *nsp, const cs_cell_mesh_t *cm, const cs_cdofb_navsto_builder_t *nsb, cs_cell_sys_t *csys)
Take into account the buoyancy force with the Boussinesq approx. Compute and add the source term to t...
Definition: cs_cdofb_navsto.cpp:2650
cs_cdofb_navsto_check_convergence
void cs_cdofb_navsto_check_convergence(const cs_navsto_param_t *nsp, const cs_cdo_quantities_t *quant, const cs_time_step_t *ts, const cs::cdo_navsto_ctx_t *sc, const cs_turbulence_t *tbs, cs_cdo_navsto_psteady_cvg_t &ps_cvg)
Check the convergence of the pseudo-steady algorithm when the unsteady Navier-Stokes system with a CD...
Definition: cs_cdofb_navsto.cpp:3234
cs_cdofb_navsto_stream_source_term
void cs_cdofb_navsto_stream_source_term(cs_lnum_t n_elts, const cs_lnum_t *elt_ids, bool dense_output, void *input, cs_real_t *retval)
Get the source term for computing the stream function. This relies on the prototype associated to the...
Definition: cs_cdofb_navsto.cpp:2767
cs_cdofb_navsto_init_face_pressure
void cs_cdofb_navsto_init_face_pressure(const cs_navsto_param_t *nsp, const cs_cdo_connect_t *connect, const cs_time_step_t *ts, cs_real_t *pr_f)
Initialize the pressure values when the pressure is defined at faces.
Definition: cs_cdofb_navsto.cpp:787
cs_cdofb_navsto_define_builder
void cs_cdofb_navsto_define_builder(cs_real_t t_eval, const cs_navsto_param_t *nsp, const cs_cell_mesh_t *cm, const cs_cell_sys_t *csys, const cs_cdo_bc_face_t *pr_bc, const cs_boundary_type_t *bf_type, cs_cdofb_navsto_builder_t *nsb)
Set the members of the cs_cdofb_navsto_builder_t structure.
Definition: cs_cdofb_navsto.cpp:263
cs_cdofb_navsto_set_gravity_func
void cs_cdofb_navsto_set_gravity_func(const cs_navsto_param_t *nsp, cs_cdofb_navsto_source_t **p_func)
Set the function pointer computing the source term in the momentum equation related to the gravity ef...
Definition: cs_cdofb_navsto.cpp:2565
cs_cdofb_navsto_extra_op
void cs_cdofb_navsto_extra_op(const cs_navsto_param_t *nsp, const cs_mesh_t *mesh, const cs_cdo_quantities_t *quant, const cs_cdo_connect_t *connect, const cs_time_step_t *ts, cs_time_plot_t *time_plotter, const cs_turbulence_t *turb, const cs_adv_field_t *adv_field, const cs_real_t *mass_flux, const cs_real_t *p_cell, const cs_real_t *u_cell, const cs_real_t *u_face, const cs_cdo_navsto_psteady_cvg_t &ps_cvg)
Perform extra-operation related to Fb schemes when solving Navier-Stokes. Computation of the followin...
Definition: cs_cdofb_navsto.cpp:1053
cs_cdofb_navsto_create_builder
cs_cdofb_navsto_builder_t cs_cdofb_navsto_create_builder(const cs_navsto_param_t *nsp, const cs_cdo_connect_t *connect)
Create and allocate a local NavSto builder when Fb schemes are used.
Definition: cs_cdofb_navsto.cpp:208
cs_cdofb_block_dirichlet_weak
void cs_cdofb_block_dirichlet_weak(short int fb, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a Dirichlet BCs on the three velocity components. For instance for a velocity inlet...
Definition: cs_cdofb_navsto.cpp:1731
cs_cdofb_navsto_set_zero_mean_pressure
void cs_cdofb_navsto_set_zero_mean_pressure(const cs_cdo_quantities_t *quant, cs_real_t values[])
Update the pressure field in order to get a field with a zero-mean average.
Definition: cs_cdofb_navsto.cpp:985
cs_cdofb_prescribed_smooth_wall_n_alge_t_neumann
void cs_cdofb_prescribed_smooth_wall_n_alge_t_neumann(short int fb, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a wall BCs by a algebraic enforcement for the normal component, and neumann for the...
Definition: cs_cdofb_navsto.cpp:2321
cs_cdofb_navsto_source_t
void() cs_cdofb_navsto_source_t(const cs_navsto_param_t *nsp, const cs_cell_mesh_t *cm, const cs_cdofb_navsto_builder_t *nsb, cs_cell_sys_t *csys)
Compute and add a source term to the local RHS. This is a special treatment to enable source involvin...
Definition: cs_cdofb_navsto.h:184
cs_cdofb_navsto_mass_flux
void cs_cdofb_navsto_mass_flux(const cs_navsto_param_t *nsp, const cs_cdo_quantities_t *quant, const cs_real_t *face_vel, cs_real_t *mass_flux)
Compute the mass flux playing the role of the advection field in the Navier-Stokes equations One cons...
Definition: cs_cdofb_navsto.cpp:432
cs_cdofb_prescribed_smooth_wall_n_alge_t_robin
void cs_cdofb_prescribed_smooth_wall_n_alge_t_robin(short int fb, const cs_equation_param_t *eqp, const cs_cell_mesh_t *cm, const cs_property_data_t *pty, cs_cell_builder_t *cb, cs_cell_sys_t *csys)
Take into account a wall BCs by an algebraic enforcement for the normal component and robin on tangen...
Definition: cs_cdofb_navsto.cpp:2129
cs_cdofb_navsto_balance
cs_cdo_balance_t * cs_cdofb_navsto_balance(const cs_navsto_param_t *nsp, const cs_cdo_connect_t *connect, const cs_cdo_quantities_t *quant, const cs_time_step_t *ts, const cs_cdo_bc_face_t *pr_bc, const cs_boundary_type_t *bf_type, const cs_real_t *pr_c)
Compute the balance for an NS system over the full computational domain between time t_cur and t_cur ...
Definition: cs_cdofb_navsto.cpp:2808
cs_cdofb_navsto_boussinesq_by_part
void cs_cdofb_navsto_boussinesq_by_part(const cs_navsto_param_t *nsp, const cs_cell_mesh_t *cm, const cs_cdofb_navsto_builder_t *nsb, cs_cell_sys_t *csys)
Take into account the buoyancy force with the Boussinesq approx. Compute and add the source term to t...
cs_defs.h
cs_real_t
double cs_real_t
Floating-point value.
Definition: cs_defs.h:332
cs_lnum_t
int cs_lnum_t
local mesh entity id
Definition: cs_defs.h:325
cs_field.h
cs_iter_algo.h
cs_math.h
cs_matrix.h
cs_mesh.h
cs_navsto_context.h
cs_navsto_param.h
cs_sdm.h
cs_sles_convergence_state_t
cs_sles_convergence_state_t
Definition: cs_sles.h:57
cs_time_plot.h
cs_time_plot_t
struct _cs_time_plot_t cs_time_plot_t
Definition: cs_time_plot.h:44
cs_time_step.h
algo
char * algo
Definition: field_names.h:106
mesh
Definition: mesh.f90:26
cs::cdo_navsto_ctx_t
Context related to CDO face-based discretization when dealing with Navier-Stokes equations and vector...
Definition: cs_navsto_context.h:48
cs_adjacency_t
Definition: cs_mesh_adjacencies.h:68
cs_adv_field_t
Definition: cs_advection_field.h:151
cs_cdo_balance_t
Definition: cs_cdo_toolbox.h:64
cs_cdo_bc_face_t
Definition: cs_cdo_bc.h:105
cs_cdo_connect_t
Definition: cs_cdo_connect.h:57
cs_cdo_navsto_psteady_cvg_t
Definition: cs_param_types.h:660
cs_cdo_quantities_t
Definition: cs_cdo_quantities.h:142
cs_cdofb_navsto_builder_t
Structure storing additional arrays related to the building of the Navier-Stokes system.
Definition: cs_cdofb_navsto.h:105
cs_cdofb_navsto_builder_t::mass_rhs
cs_real_t mass_rhs
Definition: cs_cdofb_navsto.h:139
cs_cdofb_navsto_builder_t::div_op
cs_real_t * div_op
Definition: cs_cdofb_navsto.h:128
cs_cdofb_navsto_builder_t::rho_c
cs_real_t rho_c
Definition: cs_cdofb_navsto.h:115
cs_cdofb_navsto_builder_t::bf_type
cs_boundary_type_t * bf_type
Definition: cs_cdofb_navsto.h:159
cs_cdofb_navsto_builder_t::pressure_bc_val
cs_real_t * pressure_bc_val
Definition: cs_cdofb_navsto.h:160
cs_cell_builder_t
Set of local and temporary buffers.
Definition: cs_cdo_local.h:56
cs_cell_mesh_t
Set of local quantities and connectivities related to a mesh cell.
Definition: cs_cdo_local.h:242
cs_cell_sys_t
Set of arrays and local (small) dense matrices related to a mesh cell This is a key structure for bui...
Definition: cs_cdo_local.h:163
cs_equation_param_t
Set of parameters to handle an unsteady convection-diffusion-reaction equation with term sources.
Definition: cs_equation_param.h:190
cs_iter_algo_t
Structure to handle the convergence of an iterative algorithm.
Definition: cs_iter_algo.h:285
cs_mesh_t
Definition: cs_mesh.h:85
cs_navsto_param_t
Structure storing the parameters related to the resolution of the Navier-Stokes system.
Definition: cs_navsto_param.h:296
cs_property_data_t
Structure storing the evaluation of a property and its related data.
Definition: cs_property.h:219
cs_turbulence_t
Structure storing the parameters related to the resolution of the turbulence modelling....