programmer's documentation
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
Functions
cs_gwf_physical_properties.c File Reference
#include "cs_defs.h"
#include <assert.h>
#include <stdio.h>
#include <stdarg.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include "bft_mem.h"
#include "bft_error.h"
#include "bft_printf.h"
#include "cs_gwf_parameters.h"
#include "cs_math.h"
#include "cs_field.h"
#include "cs_field_pointer.h"
#include "cs_mesh.h"
#include "cs_gwf_physical_properties.h"
Include dependency graph for cs_gwf_physical_properties.c:

Functions

void cs_gwf_physical_properties (void)
 Update physical properties of the Ground Water Flow module. More...
 
void cs_gwf_sorbed_concentration_update (cs_real_t *c_scal, cs_real_t *kp, cs_real_t *km, cs_real_t *sorb)
 Update sorbed concentration for scalars with kinetic sorption. More...
 

Function Documentation

void cs_gwf_physical_properties ( void  )

Update physical properties of the Ground Water Flow module.

Species transport is delayed by retention in solid phase. This delay is computed as follows:

\[ R = 1 + \rho \dfrac{K_d}{\theta} ; \]

where $ R $ is the delay factor, $ \rho $ the soil (bulk) density, $ K_d $ the contaminant distribution coefficient and $ \theta $ the moisture content.

Please refer to the dedicated section of the theory guide for more informations.

void cs_gwf_sorbed_concentration_update ( cs_real_t c_scal,
cs_real_t kp,
cs_real_t km,
cs_real_t sorb 
)

Update sorbed concentration for scalars with kinetic sorption.

It is estimated by the following analytical expression :

\[ S^{n+1} = S^n \exp(- k^{-} \Delta t) - C^n \dfrac{k^{+}}{k^{-}} \left(\exp(- k^{-} \Delta t) - 1 \right) \]

Parameters
[in]c_scalconcentration field
[in]kpkplus field
[in]kmkminus field
[in]sorbsorbed concentration field