Initial commit

Intermediate_Canyon/cs_meg_boundary_function.c

+/*----------------------------------------------------------------------------*/
+/*
+  This file is generated by Code_Saturne, a general-purpose CFD tool.
+*/
+/*----------------------------------------------------------------------------*/
+
+#include "cs_defs.h"
+
+/*----------------------------------------------------------------------------
+ * Standard C library headers
+ *----------------------------------------------------------------------------*/
+
+#include <assert.h>
+#include <math.h>
+
+#if defined(HAVE_MPI)
+#include <mpi.h>
+#endif
+
+/*----------------------------------------------------------------------------
+ *  Local headers
+ *----------------------------------------------------------------------------*/
+
+#include "cs_headers.h"
+/*----------------------------------------------------------------------------*/
+
+BEGIN_C_DECLS
+
+/*----------------------------------------------------------------------------*/
+
+cs_real_t *
+cs_meg_boundary_function(const cs_zone_t *zone,
+                         const char      *field_name,
+                         const char      *condition)
+{
+  cs_real_t *new_vals = NULL;
+
+  /*--------------------------------------------------*/
+  /* User defined formula for "velocity" over BC=BC_2 */
+  if (strcmp(field_name, "velocity") == 0 &&
+      strcmp(condition, "norm_formula") == 0 &&
+      strcmp(zone->name, "BC_2") == 0) {
+
+    const cs_real_3_t *xyz = (cs_real_3_t *)cs_glob_mesh_quantities->b_face_cog;
+
+    const int vals_size = zone->n_elts * 1;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+    const cs_real_t H = cs_notebook_parameter_value_by_name("H");
+    const cs_real_t zref = cs_notebook_parameter_value_by_name("zref");
+    const cs_real_t Uref = cs_notebook_parameter_value_by_name("Uref");
+    const cs_real_t z0 = cs_notebook_parameter_value_by_name("z0");
+    const cs_real_t k = cs_notebook_parameter_value_by_name("k");
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t f_id = zone->elt_ids[e_id];
+      const cs_real_t z = xyz[f_id][2];
+
+      cs_real_t ustar = k * Uref / log((zref - H + z0) / z0) ;
+      new_vals[0 * zone->n_elts + e_id] =  ustar / k * log((z - H + z0) / z0) ;
+    }
+  }
+  /*--------------------------------------------------*/
+
+  /*---------------------------------------------------*/
+  /* User defined formula for "direction" over BC=BC_2 */
+  if (strcmp(field_name, "direction") == 0 &&
+      strcmp(condition, "formula") == 0 &&
+      strcmp(zone->name, "BC_2") == 0) {
+
+    const int vals_size = zone->n_elts * 3;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+    const cs_real_t phi = cs_notebook_parameter_value_by_name("phi");
+    const cs_real_t pi = cs_math_pi;
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t f_id = zone->elt_ids[e_id];
+
+      new_vals[0 * zone->n_elts + e_id] =  sin(phi * pi / 180) ;
+      new_vals[1 * zone->n_elts + e_id] =  cos(phi * pi / 180) ;
+      new_vals[2 * zone->n_elts + e_id] =  0 ;
+    }
+  }
+  /*---------------------------------------------------*/
+
+  /*-------------------------------------------------------*/
+  /* User defined formula for "turbulence_ke" over BC=BC_2 */
+  if (strcmp(field_name, "turbulence_ke") == 0 &&
+      strcmp(condition, "formula") == 0 &&
+      strcmp(zone->name, "BC_2") == 0) {
+
+    const cs_real_3_t *xyz = (cs_real_3_t *)cs_glob_mesh_quantities->b_face_cog;
+
+    const int vals_size = zone->n_elts * 2;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+    const cs_real_t H = cs_notebook_parameter_value_by_name("H");
+    const cs_real_t zref = cs_notebook_parameter_value_by_name("zref");
+    const cs_real_t Uref = cs_notebook_parameter_value_by_name("Uref");
+    const cs_real_t z0 = cs_notebook_parameter_value_by_name("z0");
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t f_id = zone->elt_ids[e_id];
+      const cs_real_t z = xyz[f_id][2];
+
+      cs_real_t ustar = 0.41 * Uref / log((zref - H + z0) / z0) ;
+      new_vals[0 * zone->n_elts + e_id] =  (cs_math_pow2((ustar))) / (pow(0.09, 0.5)) ;
+      new_vals[1 * zone->n_elts + e_id] =  (cs_math_pow3((ustar))) / (0.41 * (z - H + z0)) ;
+    }
+  }
+  /*-------------------------------------------------------*/
+
+  return new_vals;
+}
+
+/*----------------------------------------------------------------------------*/
+
+END_C_DECLS
+

Intermediate_Canyon/cs_meg_initialization.c

+/*----------------------------------------------------------------------------*/
+/*
+  This file is generated by Code_Saturne, a general-purpose CFD tool.
+*/
+/*----------------------------------------------------------------------------*/
+
+#include "cs_defs.h"
+
+/*----------------------------------------------------------------------------
+ * Standard C library headers
+ *----------------------------------------------------------------------------*/
+
+#include <assert.h>
+#include <math.h>
+
+#if defined(HAVE_MPI)
+#include <mpi.h>
+#endif
+
+/*----------------------------------------------------------------------------
+ *  Local headers
+ *----------------------------------------------------------------------------*/
+
+#include "cs_headers.h"
+/*----------------------------------------------------------------------------*/
+
+BEGIN_C_DECLS
+
+/*----------------------------------------------------------------------------*/
+
+cs_real_t *
+cs_meg_initialization(const cs_zone_t *zone,
+                      const char      *field_name)
+{
+  cs_real_t *new_vals = NULL;
+
+  /*-----------------------------------------------------------------------*/
+  /* User defined initialization for variable velocity over zone all_cells */
+  if (strcmp(zone->name, "all_cells") == 0 &&
+      strcmp(field_name, "velocity") == 0 ) {
+
+    const int vals_size = zone->n_elts * 3;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+
+      new_vals[3 * e_id + 0] =  0. ;
+      new_vals[3 * e_id + 1] =  0. ;
+      new_vals[3 * e_id + 2] =  0. ;
+    }
+  }
+  /*-----------------------------------------------------------------------*/
+
+  /*----------------------------------------------------------------------*/
+  /* User defined initialization for variable thermal over zone all_cells */
+  if (strcmp(zone->name, "all_cells") == 0 &&
+      strcmp(field_name, "thermal") == 0 ) {
+
+    const int vals_size = zone->n_elts * 1;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+
+      new_vals[e_id] =  293.15 ;
+    }
+  }
+  /*----------------------------------------------------------------------*/
+
+  /*----------------------------------------------------------------------*/
+  /* User defined initialization for variable traceur over zone all_cells */
+  if (strcmp(zone->name, "all_cells") == 0 &&
+      strcmp(field_name, "traceur") == 0 ) {
+
+    const int vals_size = zone->n_elts * 1;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+
+      new_vals[e_id] =  0 ;
+    }
+  }
+  /*----------------------------------------------------------------------*/
+
+  /*------------------------------------------------------------------------*/
+  /* User defined initialization for variable velocity over zone dessus_rue */
+  if (strcmp(zone->name, "dessus_rue") == 0 &&
+      strcmp(field_name, "velocity") == 0 ) {
+
+    const cs_real_3_t *xyz = (cs_real_3_t *)cs_glob_mesh_quantities->cell_cen;
+
+    const int vals_size = zone->n_elts * 3;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+    const cs_real_t H = cs_notebook_parameter_value_by_name("H");
+    const cs_real_t zref = cs_notebook_parameter_value_by_name("zref");
+    const cs_real_t Uref = cs_notebook_parameter_value_by_name("Uref");
+    const cs_real_t z0 = cs_notebook_parameter_value_by_name("z0");
+    const cs_real_t phi = cs_notebook_parameter_value_by_name("phi");
+    const cs_real_t pi = cs_math_pi;
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+      const cs_real_t z = xyz[c_id][2];
+
+      new_vals[3 * e_id + 0] =  sin(phi * pi / 180) * Uref * log((z - H + z0) / z0) / log((zref - H + z0) / z0) ;
+      new_vals[3 * e_id + 1] =  cos(phi * pi / 180) * Uref * log((z - H + z0) / z0) / log((zref - H + z0) / z0) ;
+      new_vals[3 * e_id + 2] =  0. ;
+    }
+  }
+  /*------------------------------------------------------------------------*/
+
+  /*--------------------------------------------------------------------------*/
+  /* User defined initialization for variable turbulence over zone dessus_rue */
+  if (strcmp(zone->name, "dessus_rue") == 0 &&
+      strcmp(field_name, "turbulence") == 0 ) {
+
+    const cs_real_3_t *xyz = (cs_real_3_t *)cs_glob_mesh_quantities->cell_cen;
+
+    const int vals_size = zone->n_elts * 2;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+    const cs_real_t H = cs_notebook_parameter_value_by_name("H");
+    const cs_real_t zref = cs_notebook_parameter_value_by_name("zref");
+    const cs_real_t Uref = cs_notebook_parameter_value_by_name("Uref");
+    const cs_real_t z0 = cs_notebook_parameter_value_by_name("z0");
+    const cs_real_t k = cs_notebook_parameter_value_by_name("k");
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+      const cs_real_t z = xyz[c_id][2];
+
+      cs_real_t ustar = k * Uref / (log((zref - H + z0) / z0)) ;
+      new_vals[2 * e_id + 0] =  (cs_math_pow2((ustar))) / (pow(0.09, 0.5)) ;
+      new_vals[2 * e_id + 1] =  (cs_math_pow3((ustar))) / (k * (z - H + z0)) ;
+    }
+  }
+  /*--------------------------------------------------------------------------*/
+
+  /*-----------------------------------------------------------------------*/
+  /* User defined initialization for variable thermal over zone dessus_rue */
+  if (strcmp(zone->name, "dessus_rue") == 0 &&
+      strcmp(field_name, "thermal") == 0 ) {
+
+    const int vals_size = zone->n_elts * 1;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+
+      new_vals[e_id] =  293.15 ;
+    }
+  }
+  /*-----------------------------------------------------------------------*/
+
+  /*-----------------------------------------------------------------------*/
+  /* User defined initialization for variable traceur over zone dessus_rue */
+  if (strcmp(zone->name, "dessus_rue") == 0 &&
+      strcmp(field_name, "traceur") == 0 ) {
+
+    const int vals_size = zone->n_elts * 1;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+
+      new_vals[e_id] =  0 ;
+    }
+  }
+  /*-----------------------------------------------------------------------*/
+
+  return new_vals;
+}
+
+/*----------------------------------------------------------------------------*/
+
+END_C_DECLS
+

Intermediate_Canyon/cs_meg_source_terms.c

+/*----------------------------------------------------------------------------*/
+/*
+  This file is generated by Code_Saturne, a general-purpose CFD tool.
+*/
+/*----------------------------------------------------------------------------*/
+
+#include "cs_defs.h"
+
+/*----------------------------------------------------------------------------
+ * Standard C library headers
+ *----------------------------------------------------------------------------*/
+
+#include <assert.h>
+#include <math.h>
+
+#if defined(HAVE_MPI)
+#include <mpi.h>
+#endif
+
+/*----------------------------------------------------------------------------
+ *  Local headers
+ *----------------------------------------------------------------------------*/
+
+#include "cs_headers.h"
+/*----------------------------------------------------------------------------*/
+
+BEGIN_C_DECLS
+
+/*----------------------------------------------------------------------------*/
+
+cs_real_t *
+cs_meg_source_terms(const cs_zone_t *zone,
+                    const char      *name,
+                    const char      *source_type)
+{
+  cs_real_t *new_vals = NULL;
+
+  /*----------------------------------------------------*/
+  /* User defined source term for traceur over zone rue */
+  if (strcmp(zone->name, "rue") == 0 &&
+      strcmp(name, "traceur") == 0 && 
+      strcmp(source_type, "scalar_source_term") == 0 ) {
+
+    const int vals_size = zone->n_elts * 2;
+    BFT_MALLOC(new_vals, vals_size, cs_real_t);
+
+    const cs_real_t volume = zone->measure;
+
+    for (cs_lnum_t e_id = 0; e_id < zone->n_elts; e_id++) {
+      cs_lnum_t c_id = zone->elt_ids[e_id];
+
+      new_vals[2 * e_id + 0] =  1 / volume ;
+      new_vals[2 * e_id + 1] =  0 ;
+    }
+  }
+  /*----------------------------------------------------*/
+
+  return new_vals;
+}
+
+/*----------------------------------------------------------------------------*/
+
+END_C_DECLS
+

Intermediate_Canyon/cs_user_extra_operations_IC.c

+/*============================================================================
+ * This function is called at the end of each time step, and has a very
+ *  general purpose
+ *  (i.e. anything that does not have another dedicated user function)
+ *============================================================================*/
+
+/* Code_Saturne version 6.0-beta */
+
+/*
+  This file is part of Code_Saturne, a general-purpose CFD tool.
+
+  Copyright (C) 1998-2019 EDF S.A.
+
+  This program is free software; you can redistribute it and/or modify it under
+  the terms of the GNU General Public License as published by the Free Software
+  Foundation; either version 2 of the License, or (at your option) any later
+  version.
+
+  This program is distributed in the hope that it will be useful, but WITHOUT
+  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
+  details.
+
+  You should have received a copy of the GNU General Public License along with
+  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
+  Street, Fifth Floor, Boston, MA 02110-1301, USA.
+*/
+
+/*----------------------------------------------------------------------------*/
+
+#include "cs_defs.h"
+#include "cs_math.h"
+
+/*----------------------------------------------------------------------------
+ * Standard C library headers
+ *----------------------------------------------------------------------------*/
+
+#include <assert.h>
+#include <math.h>
+
+#if defined(HAVE_MPI)
+#include <mpi.h>
+#endif
+
+/*----------------------------------------------------------------------------
+ * PLE library headers
+ *----------------------------------------------------------------------------*/
+
+#include <ple_coupling.h>
+
+/*----------------------------------------------------------------------------
+ *  Local headers
+ *----------------------------------------------------------------------------*/
+
+#include "cs_headers.h"
+
+/*----------------------------------------------------------------------------
+ *  Header for the current file
+ *----------------------------------------------------------------------------*/
+
+#include "cs_prototypes.h"
+
+/*----------------------------------------------------------------------------*/
+
+BEGIN_C_DECLS
+
+/*----------------------------------------------------------------------------*/
+/*!
+ * \file cs_user_extra_operations.c
+ *
+ * \brief This function is called at the end of each time step, and has a very
+ * general purpose (i.e. anything that does not have another dedicated
+ * user function)
+ */
+/*----------------------------------------------------------------------------*/
+
+/*============================================================================
+ * User function definitions
+ *============================================================================*/
+
+/*----------------------------------------------------------------------------*/
+/*!
+ * \brief This function is called at the end of each time step.
+ *
+ * It has a very general purpose, although it is recommended to handle
+ * mainly postprocessing or data-extraction type operations.
+ *
+ * \param[in, out]  domain   pointer to a cs_domain_t structure
+ */
+/*----------------------------------------------------------------------------*/
+
+void
+cs_user_extra_operations(cs_domain_t     *domain)
+{
+  /* Get number of cells */
+  const cs_mesh_t *m = domain->mesh;
+  const cs_lnum_t n_cells = m->n_cells;
+
+  /* Get physical fields */
+  const cs_real_t *prop = cs_field_by_name("traceur")->val;
+  const cs_real_t *densite = cs_field_by_name("density")->val;
+  const cs_real_3_t *vel = cs_field_by_name("velocity")->val;
+  const cs_real_t *k = cs_field_by_name("k")->val;
+  const cs_real_t *epsilon = cs_field_by_name("epsilon")->val;
+
+  /* Name of the first zone and first output file */
+  char zone[100] = "x>27.5 and x<55 and z>0 and z<14";  
+  char fname[100] = "concentration_traceur_rue.dat";
+  /* Name of the second zone and output file */
+  char zone2[100] = "x>27.5 and x<55 and z>14 and z<28";
+  char fname2[100] = "concentration_haut_rue.dat";
+
+  /* Write header at the first time step */
+  FILE *f = NULL;
+  FILE *f2 = NULL;
+  if (cs_glob_time_step->nt_cur == cs_glob_time_step->nt_ini && cs_glob_rank_id <= 0) {
+    f = fopen(fname,"w");
+    fprintf(f, "#%17s%100s\n", "zone : ", zone);
+    fprintf(f, "#%s\n", "temps (sec), traceur (microg/m3), Vx, Vy, Vz, densité (kg/m3), k, epsilon");
+    fclose(f);
+    f2 = fopen(fname2,"w");                                                                      
+    fprintf(f2, "#%17s%100s\n", "zone : ", zone2);                                               
+    fprintf(f2, "#%s\n", "temps (sec), traceur (microg/m3), Vx, Vy, Vz, densité (kg/m3), k, epsilon");
+    fclose(f2);
+  }
+
+  /* Open file and print time */
+  if (cs_glob_rank_id <= 0) {
+    f = fopen(fname,"a");
+    fprintf(f, "%17.9e ", cs_glob_time_step->t_cur);
+    f2 = fopen(fname2,"a");
+    fprintf(f2, "%17.9e ", cs_glob_time_step->t_cur);
+  }
+
+  /* Allocate memory */
+  cs_lnum_t n_cells_sel = 0;
+  cs_lnum_t *cells_sel_ids = NULL;
+  BFT_MALLOC(cells_sel_ids, m->n_cells, cs_lnum_t);
+
+  /* Temporary variables */
+  cs_real_t sum, rrho, sumx, sumy, sumz, sumk, sumepsilon, sumvolume;
+
+  /* Get cells for the first zone */
+  cs_selector_get_cell_list(zone, &n_cells_sel, cells_sel_ids);
+
+  /* Average various quantities over the first zone */
+  {
+    sum = 0;
+    rrho = 0;
+    sumx = 0;
+    sumy = 0;
+    sumz = 0;
+    sumk = 0;
+    sumepsilon = 0;
+    sumvolume = 0;
+    const cs_real_t *restrict volume = domain->mesh_quantities->cell_f_vol;
+    for (cs_lnum_t iel = 0; iel < n_cells_sel; iel++) {
+      const cs_lnum_t cell_id = cells_sel_ids[iel];
+      sum += 1000. * prop[cell_id] * densite[cell_id] * volume[cell_id];
+      rrho += densite[cell_id] * volume[cell_id];
+      sumx += vel[cell_id][0] * volume[cell_id];
+      sumy += vel[cell_id][1] * volume[cell_id];
+      sumz += vel[cell_id][2] * volume[cell_id];
+      sumk += k[cell_id] * volume[cell_id];
+      sumepsilon += epsilon[cell_id] * volume[cell_id];
+      sumvolume += volume[cell_id];
+    }
+  }
+
+  /* Parallel sum then write to file */
+  cs_parall_sum(1, CS_DOUBLE, &sum);
+  cs_parall_sum(1, CS_DOUBLE, &rrho);
+  cs_parall_sum(1, CS_DOUBLE, &sumx);
+  cs_parall_sum(1, CS_DOUBLE, &sumy);
+  cs_parall_sum(1, CS_DOUBLE, &sumz);
+  cs_parall_sum(1, CS_DOUBLE, &sumvolume);
+  cs_parall_sum(1, CS_DOUBLE, &sumk);
+  cs_parall_sum(1, CS_DOUBLE, &sumepsilon);
+  if (cs_glob_rank_id <= 0) fprintf(f," %17.9e", sum/sumvolume);
+  if (cs_glob_rank_id <= 0) fprintf(f," %17.9e", sumx/sumvolume);
+  if (cs_glob_rank_id <= 0) fprintf(f," %17.9e", sumy/sumvolume);
+  if (cs_glob_rank_id <= 0) fprintf(f," %17.9e", sumz/sumvolume);
+  if (cs_glob_rank_id <= 0) fprintf(f," %17.9e", rrho/sumvolume);
+  if (cs_glob_rank_id <= 0) fprintf(f," %17.9e", sumk/sumvolume);
+  if (cs_glob_rank_id <= 0) fprintf(f," %17.9e", sumepsilon/sumvolume);
+
+  /* Print end of line and close the file */
+  if (cs_glob_rank_id <= 0) {
+    fprintf(f, "\n");
+    fclose(f);
+  }
+
+  /* Get cells for the second zone */
+  cs_selector_get_cell_list(zone2, &n_cells_sel, cells_sel_ids);
+
+  /* Average various quantities over the second zone */
+  {
+    sum = 0;
+    rrho = 0;
+    sumx = 0;
+    sumy = 0;
+    sumz = 0;
+    sumk = 0;
+    sumepsilon = 0;
+    sumvolume = 0;
+    const cs_real_t *restrict volume = domain->mesh_quantities->cell_f_vol;
+    for (cs_lnum_t iel = 0; iel < n_cells_sel; iel++) {
+      const cs_lnum_t cell_id = cells_sel_ids[iel];
+      sum += 1000. * prop[cell_id] * densite[cell_id] * volume[cell_id];
+      rrho += densite[cell_id] * volume[cell_id];
+      sumx += vel[cell_id][0] * volume[cell_id];
+      sumy += vel[cell_id][1] * volume[cell_id];
+      sumz += vel[cell_id][2] * volume[cell_id];
+      sumk += k[cell_id] * volume[cell_id];
+      sumepsilon += epsilon[cell_id] * volume[cell_id];
+      sumvolume += volume[cell_id];
+    }
+  }
+