/*============================================================================ * 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) *============================================================================*/ /* VERS */ /* This file is part of Code_Saturne, a general-purpose CFD tool. Copyright (C) 1998-2021 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" /*---------------------------------------------------------------------------- * Standard C library headers *----------------------------------------------------------------------------*/ #include #include #include #if defined(HAVE_MPI) #include #endif /*---------------------------------------------------------------------------- * PLE library headers *----------------------------------------------------------------------------*/ #include /*---------------------------------------------------------------------------- * Local headers *----------------------------------------------------------------------------*/ #include "cs_headers.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) { CS_UNUSED(domain); FILE *file = NULL; const cs_lnum_t n_b_faces = cs_glob_mesh->n_b_faces; const cs_real_3_t *b_face_cog = (const cs_real_3_t *)cs_glob_mesh_quantities->b_face_cog; const cs_real_3_t *surfbo = (const cs_real_3_t *)cs_glob_mesh_quantities->b_face_normal; const cs_real_t *surfbn = cs_glob_mesh_quantities->b_face_surf; //const cs_real_t *f_b_temp // = (const cs_real_t *)cs_field_by_name("boundary_temperature")->val; const cs_real_t *f_b_temp = (const cs_real_t *)cs_field_by_name("T_Mean")->val; const cs_real_3_t *forbr = (const cs_real_3_t *)cs_field_by_name("boundary_forces")->val; //cs_field_t *f = cs_thermal_model_field(); cs_field_t *f = (cs_real_t *)cs_field_by_name("T_Mean"); double visls_0 = 0.027316647887; //const double visls_0 = // = cs_field_get_key_double(f, cs_field_key_id("diffusivity_ref")); if (cs_glob_time_step->nt_cur == cs_glob_time_step->nt_max) { /* Compute the thermal fluxes at selected boundary faces */ const int location_id = CS_MESH_LOCATION_BOUNDARY_FACES; const cs_lnum_t n_elts = cs_mesh_location_get_n_elts(location_id)[0]; cs_real_t *boundary_flux = NULL; BFT_MALLOC(boundary_flux, n_elts, cs_real_t); cs_post_boundary_flux(f->name, n_elts, NULL, boundary_flux); /* Some declarations */ cs_real_t srfbn, srfnor[3], fornor, stresses[3]; cs_real_t *loc_nusselt = NULL, *loc_friction = NULL, *loc_coords1 = NULL, *loc_coords2 = NULL, *loc_coords3 = NULL; cs_real_t *glo_nusselt = NULL, *glo_friction = NULL, *glo_coords1 = NULL, *glo_coords2 = NULL, *glo_coords3 = NULL; /* Print Nusselt and friction coeff file header*/ if (cs_glob_rank_id <= 0) { file = fopen("Surface_values.dat","w"); fprintf(file, "# This routine writes values at walls\n"); fprintf(file, "# 1:Coords, 2:Cf, 3:Nu \n"); } /* Select boundary faces to print */ /* -------- TO MODIFY ---------- */ const char criteria[] = "bc_cold_2-wall, bc_cold_1-natural_convection-67, bc_cold_1-natural_convection-68, bc_cold_1-natural_convection-69, bc_cold_1-natural_convection-70, bc_cold_1-natural_convection-71, bc_cold_1-natural_convection-72, bc_hot-natural_convection_ca-57, bc_hot-natural_convection_ca-58, bc_hot-natural_convection_ca-59, bc_hot-natural_convection_ca-60, bc_hot-natural_convection_ca-61, bc_hot-natural_convection_ca-62"; /* -------- TO MODIFY ---------- */ cs_lnum_t n_selected_faces = 0; cs_gnum_t n_selected_faces_g = 0; cs_lnum_t *selected_faces = NULL; BFT_MALLOC(selected_faces, n_b_faces, cs_lnum_t); cs_selector_get_b_face_list(criteria, &n_selected_faces, selected_faces); /* Get the total number of faces shared on all ranks */ n_selected_faces_g = n_selected_faces; cs_parall_sum(1 , CS_GNUM_TYPE, &n_selected_faces_g); /*Allocate local and global arrays to store the desired data */ BFT_MALLOC(loc_nusselt, n_selected_faces, cs_real_t); BFT_MALLOC(loc_friction, n_selected_faces, cs_real_t); BFT_MALLOC(loc_coords1, n_selected_faces, cs_real_t); BFT_MALLOC(loc_coords2, n_selected_faces, cs_real_t); BFT_MALLOC(loc_coords3, n_selected_faces, cs_real_t); if (cs_glob_n_ranks > 1) { BFT_MALLOC(glo_nusselt , n_selected_faces_g, cs_real_t); BFT_MALLOC(glo_friction, n_selected_faces_g, cs_real_t); BFT_MALLOC(glo_coords1, n_selected_faces_g, cs_real_t); BFT_MALLOC(glo_coords2, n_selected_faces_g, cs_real_t); BFT_MALLOC(glo_coords3, n_selected_faces_g, cs_real_t); } /* Add some reference values */ /* -------- TO MODIFY ---------- */ cs_real_t length_ref = 0.8; cs_real_t temp_ref = 48.0; /* -------- TO MODIFY ---------- */ for (cs_lnum_t ielt = 0; ielt < n_selected_faces; ielt++) { cs_lnum_t f_id = selected_faces[ielt]; /* Compute the Nusselt number */ cs_real_t tfac = f_b_temp[f_id]; loc_nusselt[ielt] = boundary_flux[f_id] * length_ref / (visls_0 * (tfac-temp_ref)); //loc_nusselt[ielt] = boundary_flux[f_id]/(visls_0); /* Compute the friction coefficient */ srfbn = surfbn[f_id]; srfnor[0] = surfbo[f_id][0] / srfbn; srfnor[1] = surfbo[f_id][1] / srfbn; srfnor[2] = surfbo[f_id][2] / srfbn; fornor = forbr[f_id][0]*srfnor[0] + forbr[f_id][1]*srfnor[1] + forbr[f_id][2]*srfnor[2]; stresses[0] = (forbr[f_id][0] - fornor*srfnor[0]) / srfbn; stresses[1] = (forbr[f_id][1] - fornor*srfnor[1]) / srfbn; stresses[2] = (forbr[f_id][2] - fornor*srfnor[2]) / srfbn; loc_friction[ielt] = cs_math_3_norm(stresses); /* Here we plot the results with respect to X */ loc_coords1[ielt] = b_face_cog[f_id][0]; loc_coords2[ielt] = b_face_cog[f_id][1]; loc_coords3[ielt] = b_face_cog[f_id][2]; } /* Gather the data of all ranks */ if (cs_glob_n_ranks > 1) { cs_parall_allgather_r(n_selected_faces, n_selected_faces_g, loc_nusselt, glo_nusselt); cs_parall_allgather_r(n_selected_faces, n_selected_faces_g, loc_friction, glo_friction); cs_parall_allgather_r(n_selected_faces, n_selected_faces_g, loc_coords1, glo_coords1); cs_parall_allgather_r(n_selected_faces, n_selected_faces_g, loc_coords2, glo_coords2); cs_parall_allgather_r(n_selected_faces, n_selected_faces_g, loc_coords3, glo_coords3); } /* Print in file */ for (cs_gnum_t ielt = 0; ielt < n_selected_faces_g; ielt++) { if (cs_glob_rank_id == -1) fprintf(file,"%17.9e %17.9e %17.9e %17.9e %17.9e\n", loc_coords1[ielt],loc_coords2[ielt],loc_coords3[ielt],loc_friction[ielt], loc_nusselt[ielt]); if (cs_glob_rank_id == 0) fprintf(file,"%17.9e %17.9e %17.9e %17.9e %17.9e\n", glo_coords1[ielt],glo_coords2[ielt],glo_coords3[ielt],glo_friction[ielt], glo_nusselt[ielt]); } /* Free allocated memory */ BFT_FREE(boundary_flux); BFT_FREE(selected_faces); BFT_FREE(loc_nusselt); BFT_FREE(loc_friction); BFT_FREE(loc_coords1); BFT_FREE(loc_coords2); BFT_FREE(loc_coords3); BFT_FREE(glo_nusselt); BFT_FREE(glo_friction); BFT_FREE(glo_coords1); BFT_FREE(glo_coords2); BFT_FREE(glo_coords3); } } /*----------------------------------------------------------------------------*/ END_C_DECLS