1 // Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors. 2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details. 3 // 4 // SPDX-License-Identifier: BSD-2-Clause 5 // 6 // This file is part of CEED: http://github.com/ceed 7 8 #ifndef stg_shur14_type_h 9 #define stg_shur14_type_h 10 11 #include <ceed.h> 12 #include "newtonian_types.h" 13 14 /* Access data arrays via: 15 * CeedScalar (*sigma)[ctx->nmodes] = (CeedScalar (*)[ctx->nmodes])&ctx->data[ctx->offsets.sigma]; 16 * CeedScalar *eps = &ctx->data[ctx->offsets.eps]; */ 17 typedef struct STGShur14Context_ *STGShur14Context; 18 struct STGShur14Context_ { 19 CeedInt nmodes; // !< Number of wavemodes 20 CeedInt nprofs; // !< Number of profile points in STGInflow.dat 21 CeedInt nynodes; // !< Number of mesh nodes in the y direction 22 CeedScalar alpha; // !< Geometric growth rate of kappa 23 CeedScalar u0; // !< Convective velocity 24 CeedScalar time; // !< Solution time 25 CeedScalar P0; // !< Inlet pressure 26 CeedScalar theta0; // !< Inlet temperature 27 bool is_implicit; // !< Whether using implicit time integration 28 bool mean_only; // !< Only apply the mean profile 29 CeedScalar dx; // !< dx used for h calculation 30 CeedScalar dz; // !< dz used for h calculation 31 bool prescribe_T; // !< Prescribe temperature weakly 32 struct NewtonianIdealGasContext_ newtonian_ctx; 33 34 struct { 35 size_t sigma, d, phi; // !< Random number set, [nmodes,3], [nmodes,3], [nmodes] 36 size_t kappa; // !< Wavemode frequencies in increasing order, [nmodes] 37 size_t prof_dw; // !< Distance to wall for Inflow Profie, [nprof] 38 size_t ubar; // !< Mean velocity, [nprof, 3] 39 size_t cij; // !< Cholesky decomposition [nprof, 6] 40 size_t eps; // !< Turbulent Disspation [nprof, 6] 41 size_t lt; // !< Tubulent Length Scale [nprof, 6] 42 size_t ynodes; // !< Locations of nodes in y direction [nynodes] 43 } offsets; // !< Holds offsets for each array in data 44 size_t total_bytes; // !< Total size of struct plus array 45 CeedScalar data[1]; // !< Holds concatenated scalar array data 46 }; 47 48 #endif 49