e_turb.h

Go to the documentation of this file.

#ifndef E_TURB_HEADER_ALREADY_INCLUDED

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "../../util/source/useful.h"
#include "../../util/source/logfile.h"
#include "../../gas_models/source/gas.h"
#include "../../flux_calc/source/flux_calc.h"
#include "../../mb_cns/source/cns_case_id.h"
#include "../../mb_cns/source/cns_work_vector.h"
#include "../../mb_cns/source/cns_debug_level.h"
#include "../../mb_cns/source/cns_cell.h"
#include "../../mb_cns/source/cns_echo.h"
#include "../../mb_cns/source/cns_bc_defs.h"
#include "./e_block.h"

/*---------------------------------------------------------------------*/

int turbulent_flux(struct block_data_3D *A);
int calc_alpha(struct block_data_3D *A);
double wall_damping(double z);
int eddy_viscosity_constants(struct block_data_3D *A);
int strain_rate(struct block_data_3D *A);
int filtered_strain_rate(struct block_data_3D *A);
int baldwin_lomax_3D(struct block_data_3D *bp);
int smagorinsky_3D(struct block_data_3D *bp);
int mixed_scale_3D(struct block_data_3D *bp);
int mixed_scale_sgs_energy(struct block_data_3D *bp);
int smagorinsky_strain_rate(struct block_data_3D *A);

/*---------------------------------------------------------------------*/

#define E_TURB_HEADER_ALREADY_INCLUDED
#endif

/* Overall characteristics of the boundary layer model */
#define COMPRESSIBILITY 1
#define TRANSITIONAL 1

/* Parameters that are not often varied */
#define C_mutm  14.0
#define C_WK    1.0
#define K       0.0168

/* A_plus should be varied for boundary layers in pressure gradients */
#define A_plus  26


/* Parameters that I (Richard) think should be varied in response 
 * to changes in the gas and Mach number.
 * These parameters are grouped according to publications 
 * or cases in which they seem to work well.
 *
 * note 1: For Richard's simulations of the entire Drummond tunnel facility, 
 *         the original coefficients described in Baldwin and Lomax (1978) 
 *         worked badly.
 * note 2: k_inner is not varied in any other study?
 *
 * note 3: One and one only set of coefficients should be active.
 */

#if 1
/* Baldwin and Lomax (1978)
 * (Used successfully in Peter's models of a Mach 8 nozzle using air)
 */
 #define C_CP    1.6
 #define C_KLEB  0.3
 #define k_inner 0.4
#endif


#if 0
/* Kim, Harloff and Sverdrup (1988) 
 * These parameters were obtained from the graphs for the conditions corresponding to Richard's 
 * simulations of the Drummond Tunnel incident shock. Other conditions can be looked up on
 * these graphs
 */
 #define C_CP    1.74
 #define C_KLEB  0.4716
 #define k_inner 0.4
#endif

#if 0
/* Other attempts by Richard at calibrating the inner layer thermal conductivity of the 
 * turbulent boundary layer */
 #define k_inner 0.24
 #define k_inner 0.20 (too high)
 #define k_inner 0.14 (shown to be too high for both single and double res simulations)
 #define k_inner 0.08 (also too low, did show a dip in the pressure right behind the reflected shock?)
 #define k_inner 0.06
 #define k_inner 0.026 (much too low - approximated laminar boundary layers)
#endif


#if 0
 /* Notes for Drummond Tunnel modelling:
  * Modifying the outer layer coefficients C_cp and C_kleb did not affect the shock speeds or
  * post shock pressures. They do affect the outer layer profile and, therefore, the contact 
  * surface shape. 
  * with k_inner = 0.18 working well for Nitrogen driving Nitrogen, the values of these outer
  * layer coefficients were increased. This increased the curvature of the contact surface,
  * to be closer to that of the experimental results. 
  * These values should be close to that of Kim, Harloff and Sverdrup (1988) above  
  */
  /* (high Ccp 8x) */
 #define C_CP    12.8
 /* (mid Ccp 4x) */
 #define C_CP    6.4
 /* (original Ccp) */
 #define C_CP    1.6
 #define C_KLEB  0.3
 #define k_inner 0.18
#endif


#if 0
 /* Optimised coefficients for Nitrogen driving Nitrogen in the Drummond Tunnel
  * Provides a good model of the turbulent boundary layer behind the Mach 2.4 incident shock
  * (used with the aug98 and may99 experiments by David Buttsworth) */
 #define C_CP    6.4
 #define C_KLEB  0.3
 #define k_inner 0.18
#endif


#if 0
 /* Optimised coefficients for Helium driving Nitrogen in the Drummond Tunnel
  * Provides a good model of the turbulent boundary layer behind the Mach 2.4 incident shock
  * (used with the feb00 experiments by David Buttsworth) */
 #define C_CP    1.6
 #define C_KLEB  0.3
 #define k_inner 0.12
#endif

---