Mapping relative to ALE
 
  const cs_lnum_t n_b_faces = domain->mesh->n_b_faces;
 
  const cs_lnum_t *b_face_vtx_idx = domain->mesh->b_face_vtx_idx;
 
  const cs_lnum_t *b_face_vtx_lst = domain->mesh->b_face_vtx_lst;
 
  const cs_lnum_t *b_face_cells = domain->mesh->b_face_cells;
 
 
  const int nt_cur = domain->time_step->nt_cur;
 
 
  
 
 
double cs_real_t
Floating-point value.
Definition: cs_defs.h:332
cs_real_t cs_real_3_t[3]
vector of 3 floating-point values
Definition: cs_defs.h:347
int cs_lnum_t
local mesh entity id
Definition: cs_defs.h:325
cs_field_t * cs_field_by_name(const char *name)
Return a pointer to a field based on its name.
Definition: cs_field.c:2489
@ dt
Definition: cs_field_pointer.h:65
#define CS_F_(e)
Macro used to return a field pointer by its enumerated value.
Definition: cs_field_pointer.h:51
cs_real_t * val
Definition: cs_field.h:152
Assign boundary conditions to boundary faces
One may use selection criteria to filter boundary case subsets.
Loop on faces from a subset.
Set the boundary condition for each face.
Calculation of displacement at current time step
 
Example 1
Example : For boundary faces of color 4 assign a fixed velocity
 
  
 
 
  
  const cs_real_t deltaa = sin(3.141596*(nt_cur-1)/50);
 
  const cs_real_t delta  = sin(3.141596*nt_cur/50.);
 
 
 
    const cs_lnum_t c_id = b_face_cells[face_id];
 
 
    mesh_u->bc_coeffs->rcodcl1[n_b_faces*0 + face_id] = 0;
 
    mesh_u->bc_coeffs->rcodcl1[n_b_faces*1 + face_id] = 0;
 
    mesh_u->bc_coeffs->rcodcl1[n_b_faces*2 + face_id] = (delta-deltaa)/
dt[c_id];
 
 
  }
@ CS_BOUNDARY_ALE_IMPOSED_VEL
Definition: cs_boundary.h:142
const cs_zone_t * cs_boundary_zone_by_name(const char *name)
Return a pointer to a boundary zone based on its name if present.
Definition: cs_boundary_zone.c:711
@ mesh_u
Definition: cs_field_pointer.h:86
Field descriptor.
Definition: cs_field.h:131
const cs_lnum_t * elt_ids
Definition: cs_zone.h:65
cs_lnum_t n_elts
Definition: cs_zone.h:64
Example 2
Example: For boundary face of color 5 assign a fixed displacement on nodes
 
 
 
         ii < b_face_vtx_idx[face_id+1];
         ii++) {
      if (impale[vtx_id] == 0) {
        disale[vtx_id][0] = 0;
        disale[vtx_id][1] = 0;
        disale[vtx_id][2] = delta;
        impale[vtx_id] = 1;
      }
    }
  }
 
  
Example 3
Example : For boundary faces of color 6 assign a sliding boundary
 
 
 
 
  }
@ CS_BOUNDARY_ALE_SLIDING
Definition: cs_boundary.h:141
  
Example 4
Example : Prescribe elsewhere a fixed boundary
 
 
 
 
 
  }
@ CS_BOUNDARY_ALE_FIXED
Definition: cs_boundary.h:140