#include <opm/core/grid.h>

Include dependency graph for trans_tpfa.h:

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

Functions
void	tpfa_htrans_compute (struct UnstructuredGrid G, const double perm, double *htrans)

void	tpfa_trans_compute (struct UnstructuredGrid G, const double htrans, double *trans)

void	tpfa_eff_trans_compute (struct UnstructuredGrid G, const double totmob, const double htrans, double trans)

Detailed Description

Routines to assist in the calculation of two-point transmissibilities.

Function Documentation

void tpfa_eff_trans_compute	(	struct UnstructuredGrid *	G,
		const double *	totmob,
		const double *	htrans,
		double *	trans
	)

Calculate effective two-point transmissibilities from one-sided, total mobility weighted, transmissibilities.

Specifically, compute the following product

$\mathsf{T}_f = \big(\frac{1}{\lambda_1t_1} + \frac{1}{\lambda_2t_2}\big)^{-1} = \lambda_1\lambda_2 \frac{t_1t_2}{t_1 + t_2}$

in which $t_1$ and $t_2$ are the one-sided, static transmissibility values connecting the cells of face f and $\lambda_1$ and $\lambda_2$ denote the total mobilities of the respective cells.

Parameters

[in]	G	Grid.
[in]	totmob	Total mobilities. One positive scalar value for each cell.
[in]	htrans	One-sided transmissibilities as defined by function tpfa_htrans_compute().
[out]	trans	Effective, two-point transmissibilities. Array of size at least `G->number_of_faces`.

void tpfa_htrans_compute	(	struct UnstructuredGrid *	G,
		const double *	perm,
		double *	htrans
	)

Calculate static, one-sided transmissibilities for use in the two-point flux approximation method.

The one-sided transmissibilities are defined by the formula

$t_i = \frac{\vec{n}_f \mathsf{K}_c \vec{c}_c}{\lVert \vec{c}_c \rVert^2}$

in which i is the half-face index corresponding to the cell-face index pair (c,f) and $\vec{c}_{cf} = \Bar{x}_f - \Bar{x}_c$ is the centroid difference vector.

Parameters

[in]	G	Grid.
[in]	perm	Permeability. One symmetric, positive definite tensor per grid cell.
[out]	htrans	One-sided transmissibilities. Array of size at least `G->cell_facepos[ G->number_of_cells ]`.

void tpfa_trans_compute	(	struct UnstructuredGrid *	G,
		const double *	htrans,
		double *	trans
	)

Compute two-point transmissibilities from one-sided transmissibilities.

The two-point transmissibilities are given by the simple, near-harmonic average (save a factor of two)

$\mathsf{T}_f = \big(\frac{1}{t_1} + \frac{1}{t_2}\big)^{-1} = \frac{t_1t_2}{t_1 + t_2}$

in which $t_1$ and $t_2$ are the one-sided transmissibilities that connect the neighbouring cells of face f.

Parameters

[in]	G	Grid.
[in]	htrans	One-sided transmissibilities as defined by function tpfa_htrans_compute().
[out]	trans	Interface, two-point transmissibilities. Array of size at least `G->number_of_faces`.

Functions

Detailed Description

Function Documentation