Local definitions and initialization
Mesh quantities
It may be useful to access some mesh adjacencies and quantities, in which case local pointers allow for more readable code:
#define restrict
Definition: cs_defs.h:124
double cs_real_t
Floating-point value.
Definition: cs_defs.h:304
cs_lnum_t cs_lnum_2_t[2]
vector of 2 local mesh-entity ids
Definition: cs_defs.h:310
cs_real_t cs_real_3_t[3]
vector of 3 floating-point values
Definition: cs_defs.h:317
Definition: cs_mesh_quantities.h:89
cs_real_t * b_face_cog
Definition: cs_mesh_quantities.h:104
cs_real_t * i_face_surf
Definition: cs_mesh_quantities.h:106
cs_real_t * i_f_face_surf
Definition: cs_mesh_quantities.h:109
cs_real_t * i_face_cog
Definition: cs_mesh_quantities.h:103
cs_real_t * b_face_normal
Definition: cs_mesh_quantities.h:97
cs_real_t * cell_cen
Definition: cs_mesh_quantities.h:91
cs_real_t * i_f_face_normal
Definition: cs_mesh_quantities.h:99
cs_real_t * b_f_face_normal
Definition: cs_mesh_quantities.h:101
cs_real_t * i_face_normal
Definition: cs_mesh_quantities.h:95
cs_lnum_2_t * i_face_cells
Definition: cs_mesh.h:110
Associated properties
Accessing cell porosity property values is required so values may be set:
cs_field_t * cs_field_by_name(const char *name)
Return a pointer to a field based on its name.
Definition: cs_field.c:2338
cs_real_t * val
Definition: cs_field.h:146
Example: define porosity by geometric zones
Individual cell porosity values can be assigned to each cell, so they may be based on groups, geometric criteria, or any other time-independent functions:
if (x < 20.)
cpro_porosi[cell_id] = 1.;
else
cpro_porosi[cell_id] = 0.5;
}
halo_type);
int cs_lnum_t
local mesh entity id
Definition: cs_defs.h:298
void cs_field_synchronize(cs_field_t *f, cs_halo_type_t halo_type)
Synchronize current parallel and periodic field values.
Definition: cs_field_operator.c:1067
cs_halo_type_t
Definition: cs_halo.h:50
@ CS_HALO_STANDARD
Definition: cs_halo.h:52
cs_lnum_t n_cells
Definition: cs_mesh.h:96
Matching face equivalent surfaces should also be assigned in a corresponding manner, for interior faces:
if (x > 19.9999)
face_porosity = 0.5;
for (int i = 0; i < 3; i++)
i_f_face_normal[face_id][i] = face_porosity * i_face_normal[face_id][i];
mq->
i_f_face_surf[face_id] = cs_math_3_norm(i_f_face_normal[face_id]);
}
cs_lnum_t n_i_faces
Definition: cs_mesh.h:97
and for boundary faces:
if (x > 19.9999)
face_porosity = 0.5;
for (int i = 0; i < 3; i++)
b_f_face_normal[face_id][i] = face_porosity * b_face_normal[face_id][i];
mq->
b_f_face_surf[face_id] = cs_math_3_norm(b_f_face_normal[face_id]);
}
cs_real_t * b_f_face_surf
Definition: cs_mesh_quantities.h:110
cs_lnum_t n_b_faces
Definition: cs_mesh.h:98