Filter.h File Reference

#include "variables.h"
#include "logging.h"

Include dependency graph for Filter.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
double	ApplyLESTestFilter (const SimCtx *simCtx, double values[3][3][3], double weights[3][3][3])
	Applies a numerical "test filter" to a 3x3x3 stencil of data points.

Function Documentation

ApplyLESTestFilter()

double ApplyLESTestFilter	(	const SimCtx *	simCtx,
		double	values[3][3][3],
		double	weights[3][3][3]
	)

Applies a numerical "test filter" to a 3x3x3 stencil of data points.

EXPLANATION of "Test Filter": In the dynamic Smagorinsky LES model, the model coefficient (Cs) is not constant. It is calculated "dynamically" by sampling the turbulent energy cascade at two different scales. The first scale is the grid filter itself (implicit in the discretization). The second, larger scale is this explicit "test filter". By comparing the behavior of the resolved flow at these two scales, the model can determine the appropriate amount of dissipation needed, and thus calculate the Smagorinsky coefficient.

This function is the primary interface for this process. It selects the appropriate filtering algorithm based on simulation settings.

Parameters

simCtx	The global simulation context, containing flags like testfilter_ik.
values	The 3x3x3 array of scalar values at the stencil points.
weights	The 3x3x3 array of weights (typically 1/cell volume) for averaging.

Returns

The resulting scalar value after applying the test filter.

Applies a numerical "test filter" to a 3x3x3 stencil of data points.

Definition at line 145 of file Filter.c.

{
    // This function acts as a dispatcher. It reads the simulation configuration
    // and calls the appropriate internal filtering routine.
    if (simCtx->testfilter_ik) {
        // If the "-testfilter_ik" flag is set, it indicates that the simulation
        // is homogeneous in the i and k directions, so we use the specialized,
        // more accurate Simpson's rule filter. The `weights` are ignored in this case.
        return ApplySimpsonRuleHomogeneousFilter(values);
    } else {
        // This is the default case for general, non-uniform, curvilinear grids.
        // The volume-weighted box filter is used to ensure correct averaging.
        return ApplyVolumeWeightedBoxFilter(values, weights);
    }
}

Here is the call graph for this function:

Here is the caller graph for this function: