Main Page | Modules | Data Structures | Directories | File List | Data Fields | Globals | Related Pages
elmer » source

e_dynamic.c File Reference

Subroutines for filtering and dynamic eddy viscosity model. More...

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include "../../util/source/useful.h"
#include "../../util/source/compiler.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 "../../mb_cns/source/mb_cfg.h"
#include "./e_block.h"
#include "../../mb_cns/source/interpolate.h"
#include "./e_turb.h"
#include "./e_filtered_deriv.h"
#include "./e_dynamic.h"

Defines

#define COPY_TO_FILTER_ARRAY(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define COPY_PRODUCT_TO_FILTER_ARRAY(WORK, QUANTITY1, QUANTITY2, ILAST, JLAST, KLAST)
#define COPY_TO_BLOCK(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define COPY_PRODUCT_TO_IFI_FILTER_ARRAY(WORK, QUANTITY1, QUANTITY2, ILAST, JLAST, KLAST)
#define COPY_TO_IFI_FILTER_ARRAY(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define COPY_IFI_TO_BLOCK(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define COPY_PRODUCT_TO_IFJ_FILTER_ARRAY(WORK, QUANTITY1, QUANTITY2, ILAST, JLAST, KLAST)
#define COPY_TO_IFJ_FILTER_ARRAY(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define COPY_IFJ_TO_BLOCK(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define COPY_PRODUCT_TO_IFK_FILTER_ARRAY(WORK, QUANTITY1, QUANTITY2, ILAST, JLAST, KLAST)
#define COPY_TO_IFK_FILTER_ARRAY(WORK, QUANTITY, ILAST, JLAST, KLAST)
#define COPY_IFK_TO_BLOCK(WORK, QUANTITY, ILAST, JLAST, KLAST)

Functions

int dynamic_eddy_viscosity (struct block_data_3D *A)
int evaluate_C (struct block_data_3D *A)
int evaluate_C_frac (struct block_data_3D *A)
int filter_centered (struct block_data_3D *A)
int filter_iface (struct block_data_3D *A)
int filter_x (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int filter_y (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int filter_z (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int filter_iface_x (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int filter_iface_y (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int filter_iface_z (struct block_data_3D *A, double ***work_array, double ***work_array_2)
int filter_coefficients (struct block_data_3D *A)
int A2work (struct block_data_3D *A, double ***work_array)
int initialise_dev_filter_memory (int nnx, int nny, int nnz)
double *** allocate_filter_work (int nnx, int nny, int nnz)
int test_filter (struct block_data_3D *A)

Variables

double *** work_deconv
int dev_filter_memory_alloc = 0
double *** work_array
double *** work_array_2

Detailed Description

Subroutines for filtering and dynamic eddy viscosity model.

Author:
AD
Version:
March 2005

Define Documentation

#define COPY_IFI_TO_BLOCK (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            A->ifi[i][j][k].QUANTITY = WORK[i][j][k] ;                      \
         }                                                                  \
      }                                                                     \
   }

#define COPY_IFJ_TO_BLOCK (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            A->ifj[i][j][k].QUANTITY = WORK[i][j][k] ;                      \
         }                                                                  \
      }                                                                     \
   }

#define COPY_IFK_TO_BLOCK (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            A->ifk[i][j][k].QUANTITY = WORK[i][j][k] ;                      \
         }                                                                  \
      }                                                                     \
   }

#define COPY_PRODUCT_TO_FILTER_ARRAY (  WORK,
QUANTITY1,
QUANTITY2,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ctr[i][j][k].QUANTITY1 * A->ctr[i][j][k].QUANTITY2; \
         }                                                                  \
      }                                                                     \
   }

#define COPY_PRODUCT_TO_IFI_FILTER_ARRAY (  WORK,
QUANTITY1,
QUANTITY2,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ifi[i][j][k].QUANTITY1 * A->ifi[i][j][k].QUANTITY2 ; \
         }                                                                  \
      }                                                                     \
   }

#define COPY_PRODUCT_TO_IFJ_FILTER_ARRAY (  WORK,
QUANTITY1,
QUANTITY2,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ifj[i][j][k].QUANTITY1 * A->ifj[i][j][k].QUANTITY2 ; \
         }                                                                  \
      }                                                                     \
   }

#define COPY_PRODUCT_TO_IFK_FILTER_ARRAY (  WORK,
QUANTITY1,
QUANTITY2,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ifk[i][j][k].QUANTITY1 * A->ifk[i][j][k].QUANTITY2 ; \
         }                                                                  \
      }                                                                     \
   }

#define COPY_TO_BLOCK (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            A->ctr[i][j][k].QUANTITY = WORK[i][j][k] ;                      \
         }                                                                  \
      }                                                                     \
   }

#define COPY_TO_FILTER_ARRAY (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ctr[i][j][k].QUANTITY ;                      \
         }                                                                  \
      }                                                                     \
   }

#define COPY_TO_IFI_FILTER_ARRAY (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ifi[i][j][k].QUANTITY ;                      \
         }                                                                  \
      }                                                                     \
   }

#define COPY_TO_IFJ_FILTER_ARRAY (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ifj[i][j][k].QUANTITY ;                      \
         }                                                                  \
      }                                                                     \
   }

#define COPY_TO_IFK_FILTER_ARRAY (  WORK,
QUANTITY,
ILAST,
JLAST,
KLAST   ) 
 

Value:

for ( i = 0; i <= ILAST; ++i ) {                                         \
      for ( j = 0; j <= JLAST; ++j ) {                                      \
         for ( k = 0; k <= KLAST; ++k) {                                    \
            WORK[i][j][k] = A->ifk[i][j][k].QUANTITY ;                      \
         }                                                                  \
      }                                                                     \
   }

Generated on Tue Oct 25 10:24:12 2005 for CFCFD by  doxygen 1.4.4