#include <stdio.h>
#include <stdlib.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"
Go to the source code of this file.
Functions | |
| int | s_function (struct block_data_3D *A) |
| Subroutines for relaxation regularisationSecond order structure function. | |
| int | s_function_chi (struct block_data_3D *A, int chi) |
| Second order structure function. | |
| int | relaxation (struct block_data_3D *A, int step_count, int global_step, int time_level) |
| Calculate the relaxation coefficient. | |
| int | filter_relax (struct block_data_3D *A) |
|
||||||||||||||||||||
|
Calculate the relaxation coefficient.
Use the time derivatives to advance the deconvolved quantities forward by time step dt. Results stored in _sf variables. Filter deconvolved values Recompute advancement of deconvolved quantities with chi_0 Estimate Relaxation coefficient (chi) ==================================================================== NB: This function requires the fluxes (and hence net flux changes) to already be calculated from the deconvolved values. It may be necessary to filter the relaxation coefficients as they can be noisy. Using the same filtering functions 3 - 5 times achieves the same results as using a Pade filter - Stolz, 2004, private communication. ==================================================================== |
|
|
Subroutines for relaxation regularisationSecond order structure function. #################################################################### ###################################################################### This file containes the subroutines used for relaxation regularisation described in 1) Stolz S, Adams NA, Kleiser L The approximate deconvolution model for large-eddy simulations of compressible flows and its applicatino to shock-turbulent boundary layer interation. PHYS FLUIDS 13 (10): 2985-3001 OCT 2001 The implementation assumes, that the y-direction is perpendicular to the x-z plane => xrechen and zrechen does not depend on j (index in y-direction). Andrew Denman, 2004 ###################################################################### ######################################################################
|
|
||||||||||||
|
Second order structure function.
This function uses the _sf variables as this is where the Euler time step results are temporarily stored so as not to overwrite the deconvolved values. |
1.4.4