7.0
general documentation
cs_post_util.h
Go to the documentation of this file.
1 #ifndef __CS_POST_UTIL_H__
2 #define __CS_POST_UTIL_H__
3 
4 /*============================================================================
5  * Postprocessing utility functions.
6  *============================================================================*/
7 
8 /*
9  This file is part of Code_Saturne, a general-purpose CFD tool.
10 
11  Copyright (C) 1998-2021 EDF S.A.
12 
13  This program is free software; you can redistribute it and/or modify it under
14  the terms of the GNU General Public License as published by the Free Software
15  Foundation; either version 2 of the License, or (at your option) any later
16  version.
17 
18  This program is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
20  FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
21  details.
22 
23  You should have received a copy of the GNU General Public License along with
24  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
25  Street, Fifth Floor, Boston, MA 02110-1301, USA.
26 */
27 
28 /*----------------------------------------------------------------------------*/
29 
30 /*----------------------------------------------------------------------------
31  * Local headers
32  *----------------------------------------------------------------------------*/
33 
34 #include "cs_defs.h"
35 #include "cs_mesh_location.h"
36 #include "cs_field_operator.h"
37 
38 /*----------------------------------------------------------------------------*/
39 
40 BEGIN_C_DECLS
41 
42 /*=============================================================================
43  * Type Definitions
44  *============================================================================*/
45 
46 typedef enum {
47 
48  CS_POST_UTIL_Q_CRITERION,
50  CS_POST_UTIL_BOUNDARY_CLASS_ID,
52  CS_POST_UTIL_N_TYPES
54 } cs_post_util_type_t;
55 
56 /*============================================================================
57  * Global variables
58  *============================================================================*/
59 
62 extern int cs_glob_post_util_flag[];
63 
64 /*============================================================================
65  * Public function prototypes
66  *============================================================================*/
67 
68 /*----------------------------------------------------------------------------*/
93 /*----------------------------------------------------------------------------*/
94 
95 void
96 cs_cell_segment_intersect_select(void *input,
97  cs_lnum_t *n_cells,
98  cs_lnum_t **cell_ids);
99 
100 /*----------------------------------------------------------------------------*/
127 /*----------------------------------------------------------------------------*/
128 
129 void
130 cs_cell_polyline_intersect_select(void *input,
131  cs_lnum_t n_points,
132  cs_lnum_t *n_cells,
133  cs_lnum_t **cell_ids,
134  cs_real_t **seg_c_len);
135 
136 /*----------------------------------------------------------------------------*/
163 /*----------------------------------------------------------------------------*/
164 
165 void
166 cs_cell_segment_intersect_probes_define(void *input,
167  cs_lnum_t *n_elts,
168  cs_real_3_t **coords,
169  cs_real_t **s);
170 
171 /*----------------------------------------------------------------------------*/
183 /*----------------------------------------------------------------------------*/
184 
185 void
186 cs_b_face_criterion_probes_define(void *input,
187  cs_lnum_t *n_elts,
188  cs_real_3_t **coords,
189  cs_real_t **s);
190 
191 /*----------------------------------------------------------------------------*/
202 /*----------------------------------------------------------------------------*/
203 
204 cs_real_t
205 cs_post_turbomachinery_head(const char *criteria_in,
206  cs_mesh_location_type_t location_in,
207  const char *criteria_out,
208  cs_mesh_location_type_t location_out);
209 
210 /*----------------------------------------------------------------------------*/
221 /*----------------------------------------------------------------------------*/
222 
223 cs_real_t
224 cs_post_moment_of_force(cs_lnum_t n_b_faces,
225  const cs_lnum_t b_face_ids[],
226  cs_real_t axis[3]);
227 
228 /*----------------------------------------------------------------------------*/
237 /*----------------------------------------------------------------------------*/
238 
239 void
240 cs_post_stress_tangential(cs_lnum_t n_b_faces,
241  const cs_lnum_t b_face_ids[],
242  cs_real_3_t stress[]);
243 
244 /*----------------------------------------------------------------------------*/
252 /*----------------------------------------------------------------------------*/
253 
254 void
255 cs_post_b_pressure(cs_lnum_t n_b_faces,
256  const cs_lnum_t b_face_ids[],
257  cs_real_t pres[]);
258 
259 /*----------------------------------------------------------------------------*/
272 /*----------------------------------------------------------------------------*/
273 
274 void
275 cs_post_evm_reynolds_stresses(cs_field_interpolate_t interpolation_type,
276  cs_lnum_t n_cells,
277  const cs_lnum_t cell_ids[],
278  const cs_real_3_t *coords,
279  cs_real_6_t *rst);
280 
281 /*----------------------------------------------------------------------------*/
292 /*----------------------------------------------------------------------------*/
293 
294 void
295 cs_post_anisotropy_invariant(cs_lnum_t n_cells,
296  const cs_lnum_t cell_ids[],
297  const cs_real_t coords[][3],
298  cs_real_t inv[][2]);
299 
300 /*----------------------------------------------------------------------------*/
317 /*----------------------------------------------------------------------------*/
318 
319 void
320 cs_post_q_criterion(const cs_lnum_t n_loc_cells,
321  const cs_lnum_t cell_ids[],
322  cs_real_t q_crit[]);
323 
324 /*----------------------------------------------------------------------------*/
335 /*----------------------------------------------------------------------------*/
336 
337 void
338 cs_post_boundary_flux(const char *scalar_name,
339  cs_lnum_t n_loc_b_faces,
340  const cs_lnum_t b_face_ids[],
341  cs_real_t b_face_flux[]);
342 
343 /*----------------------------------------------------------------------------*/
356 /*----------------------------------------------------------------------------*/
357 
358 void
359 cs_post_field_cell_to_b_face_values(cs_field_t *f,
360  cs_lnum_t n_loc_b_faces,
361  const cs_lnum_t b_face_ids[],
362  cs_real_t *b_val);
363 
364 /*----------------------------------------------------------------------------*/
365 
366 END_C_DECLS
367 
368 #endif /* __CS_POST_UTIL_H__ */
cs_post_turbomachinery_head
cs_real_t cs_post_turbomachinery_head(const char *criteria_in, cs_mesh_location_type_t location_in, const char *criteria_out, cs_mesh_location_type_t location_out)
Compute the head of a turbomachinery (total pressure increase)
Definition: cs_post_util.c:746
cs_post_moment_of_force
cs_real_t cs_post_moment_of_force(cs_lnum_t n_b_faces, const cs_lnum_t b_face_ids[], cs_real_t axis[3])
Compute the magnitude of a moment of force torque) given an axis and the stress on a specific boundar...
Definition: cs_post_util.c:880
cs_cell_segment_intersect_select
void cs_cell_segment_intersect_select(void *input, cs_lnum_t *n_cells, cs_lnum_t **cell_ids)
Select cells cut by a given segment.
Definition: cs_post_util.c:144
cs_post_boundary_flux
void cs_post_boundary_flux(const char *scalar_name, cs_lnum_t n_loc_b_faces, const cs_lnum_t b_face_ids[], cs_real_t b_face_flux[])
Compute scalar flux on a specific boundary region.
Definition: cs_post_util.c:1269
cs_real_6_t
cs_real_t cs_real_6_t[6]
vector of 6 floating-point values
Definition: cs_defs.h:322
cs_post_anisotropy_invariant
void cs_post_anisotropy_invariant(cs_lnum_t n_cells, const cs_lnum_t cell_ids[], const cs_real_t coords[][3], cs_real_t inv[][2])
Compute the invariant of the anisotropy tensor.
cs_field_t
Field descriptor.
Definition: cs_field.h:125
cs_mesh_location_type_t
cs_mesh_location_type_t
Definition: cs_mesh_location.h:60
cs_defs.h
BEGIN_C_DECLS
#define BEGIN_C_DECLS
Definition: cs_defs.h:495
cs_post_q_criterion
void cs_post_q_criterion(const cs_lnum_t n_loc_cells, const cs_lnum_t cell_ids[], cs_real_t q_crit[])
Compute the Q-criterion from Hunt et. al over each cell of a specified volume region.
Definition: cs_post_util.c:1225
cs_field_interpolate_t
cs_field_interpolate_t
Field interpolation modes.
Definition: cs_field_operator.h:54
cs_field_operator.h
cs_real_t
double cs_real_t
Floating-point value.
Definition: cs_defs.h:307
cs_mesh_location.h
cs_post_evm_reynolds_stresses
void cs_post_evm_reynolds_stresses(cs_field_interpolate_t interpolation_type, cs_lnum_t n_cells, const cs_lnum_t cell_ids[], const cs_real_3_t *coords, cs_real_6_t *rst)
Compute Reynolds stresses in case of Eddy Viscosity Models.
Definition: cs_post_util.c:1044
cs_glob_post_util_flag
int cs_glob_post_util_flag[]
cs_post_stress_tangential
void cs_post_stress_tangential(cs_lnum_t n_b_faces, const cs_lnum_t b_face_ids[], cs_real_3_t stress[])
Compute tangential stress on a specific boundary.
Definition: cs_post_util.c:918
cs_post_b_pressure
void cs_post_b_pressure(cs_lnum_t n_b_faces, const cs_lnum_t b_face_ids[], cs_real_t pres[])
Compute pressure on a specific boundary region.
Definition: cs_post_util.c:956
cs_real_3_t
cs_real_t cs_real_3_t[3]
vector of 3 floating-point values
Definition: cs_defs.h:320
cs_lnum_t
int cs_lnum_t
local mesh entity id
Definition: cs_defs.h:301
cs_cell_polyline_intersect_select
void cs_cell_polyline_intersect_select(void *input, cs_lnum_t n_points, cs_lnum_t *n_cells, cs_lnum_t **cell_ids, cs_real_t **seg_c_len)
Select cells cut by a line composed of segments.
Definition: cs_post_util.c:314
END_C_DECLS
#define END_C_DECLS
Definition: cs_defs.h:496
CS_POST_UTIL_N_TYPES
Definition: cs_post_util.h:52
cs_cell_segment_intersect_probes_define
void cs_cell_segment_intersect_probes_define(void *input, cs_lnum_t *n_elts, cs_real_3_t **coords, cs_real_t **s)
Define probes based on the centers of cells intersected by a given segment.
Definition: cs_post_util.c:635
CS_POST_UTIL_Q_CRITERION
Definition: cs_post_util.h:48
cs_post_util_type_t
cs_post_util_type_t
Definition: cs_post_util.h:46
CS_POST_UTIL_BOUNDARY_CLASS_ID
Definition: cs_post_util.h:50
cs_post_field_cell_to_b_face_values
void cs_post_field_cell_to_b_face_values(cs_field_t *f, cs_lnum_t n_loc_b_faces, const cs_lnum_t b_face_ids[], cs_real_t *b_val)
Compute values at a selection of boundary faces of a given field located on cells.
Definition: cs_post_util.c:1318
cs_b_face_criterion_probes_define
void cs_b_face_criterion_probes_define(void *input, cs_lnum_t *n_elts, cs_real_3_t **coords, cs_real_t **s)
Define a profile based on centers of faces defined by a given criterion.
Definition: cs_post_util.c:696