/* This file is part of Code_Saturne, a general-purpose CFD tool. Copyright (C) 1998-2020 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_les_inflow-base.c * * \brief Generation of synthetic turbulence at LES inlets initialization. * * See \ref les_inflow for examples. */ /*----------------------------------------------------------------------------*/ /*============================================================================ * Local (user defined) function definitions *============================================================================*/ /*----------------------------------------------------------------------------*/ /*! * \brief Define parameters of synthetic turbulence at LES inflow. */ /*----------------------------------------------------------------------------*/ void cs_user_les_inflow_define(void) { /*! [set_restart] */ cs_les_inflow_set_restart(false, /* allow_read */ false); /* allow_write */ /*! [set_restart] */ /* Only one synthetic turbulence inlet: the Synthetic Eddy Method Batten Method is used for boundary faces of zone "Inlet" */ /*! [init_2] */ { /* Number of turbulence structures */ int n_entities = 50; /* Velocity, turbulent kinetic energy and dissipation rate */ cs_real_t vel_r[3] = {0, 0, 0}; cs_real_t k_r = 0.0; cs_real_t eps_r = 0.0; cs_les_inflow_add_inlet(CS_INFLOW_SEM, false, /* volume mode */ cs_boundary_zone_by_name("Inlet"), n_entities, 0, /* verbosity */ vel_r, k_r, eps_r); } /*! [init_2] */ } /*----------------------------------------------------------------------------*/ /*! * \brief Update of the characteristics of a given synthetic turbulence inlet. * * \param[in] zone pointer to associated boundary zone * \param[out] vel_r reference mean velocity * \param[out] k_r reference turbulent kinetic energy * \param[out] eps_r reference turbulent dissipation */ /*----------------------------------------------------------------------------*/ void cs_user_les_inflow_update(const cs_zone_t *zone, cs_real_t vel_r[3], cs_real_t *k_r, cs_real_t *eps_r) { /*! [update_1] */ if (strcmp(zone->name, "Inlet") == 0) { /* Velocity, turbulent kinetic energy and dissipation rate */ vel_r[0] = 0.1335; vel_r[1] = 0.0; vel_r[2] = 0.0; *k_r = 1e-12; *eps_r = 1e-12; } /*! [update_1] */ } /*----------------------------------------------------------------------------*/ /*! * \brief Definition of mean velocity, Reynolds stresses and dissipation rate * for each boundary face of the given synthetic turbulence inlet. * * Rij components are ordered as usual: XX, YY, ZZ, XY, YZ, XZ * * Arrays are pre-initialized before this function is called * (see \ref cs_user_les_inflow_define). * * \param[in] zone pointer to associated boundary zone * \param[in, out] vel_l velocity a zone faces * \param[in, out] rij_l reynods stresses at zone faces * \param[in, out] eps_l reference turbulent dissipation */ /*----------------------------------------------------------------------------*/ void cs_user_les_inflow_advanced(const cs_zone_t *zone, cs_real_3_t vel_l[], cs_real_6_t rij_l[], cs_real_t eps_l[]) { /* Example 2: * - Reynolds stresses and dissipation at the inlet are computed * using the turbulence intensity and standard laws for * a circular pipe for the first synthetic turbulence inlet, * - no modification of the statistics of the flow is provided for the * second synthetic turbulence inlet */ /*! [example_2] */ if (strcmp(zone->name, "Inlet") == 0) { const cs_lnum_t n_faces = zone->n_elts; const cs_real_t d2_s3 = 2.0/3.0; for (cs_lnum_t i = 0; i < n_faces; i++) { vel_l[i][0] = 0.1335; vel_l[i][1] = 0; vel_l[i][2] = 0; cs_real_t uref2 = cs_math_3_square_norm(vel_l[i]); uref2 = cs_math_fmax(uref2, 1e-12); /* Hydraulic diameter */ const cs_real_t x_hd = 2.4; /* Turbulence intensity */ const cs_real_t x_ti = 0.049; cs_real_t x_k = 1e-12; cs_real_t x_eps = 1e-12; cs_turbulence_bc_ke_turb_intensity(uref2, x_ti, x_hd, &x_k, &x_eps); rij_l[i][0] = d2_s3 * x_k; rij_l[i][1] = d2_s3 * x_k; rij_l[i][2] = d2_s3 * x_k; rij_l[i][3] = 0; rij_l[i][4] = 0; rij_l[i][5] = 0; eps_l[i] = x_eps; } } /*! [example_2] */ } /*----------------------------------------------------------------------------*/ END_C_DECLS