Brinkman Namespace Reference

Classes
class	Matrix
class	Parameters

Functions
template<int dim>
void	tangential_friction (FullMatrix< double > &M, const FEValuesBase< dim > &int_fe, double friction_coefficient)

Detailed Description

Classes and functions pertaining the discretization of coupled Brinkman-Darcy-Stokes problems.

The bilinear form considered is

\[ \begin{array}{ccccl} a(u,v) &-& b(v,p) &=& (f,v) \\ b(u,q) && &=& (g,q) \end{array} \]

The operators are

\begin{eqnarray*} a(u,v) &=& \bigl(\nu \nabla u, \nabla v\bigr) + \bigl(\rho u,v\bigr) + \bigl<\gamma\sqrt\rho u_{S,\tau},v_{S,\tau}\bigr>_{\Gamma_{SD}} + {\rm IP} + \bigl<\sigma \nabla\!\cdot\! u, \nabla\!\cdot\! v \bigr> \\ b(u,q) &=& \bigl(\nabla\!\cdot\! u,q\bigr) \end{eqnarray*}

Here, IP refers to the interior penalty face terms. The last term in the form \(a(.,.)\) is an optional grad-div stabilization.

Function Documentation

template<int dim>

void Brinkman::tangential_friction	(	FullMatrix< double > &	M,
		const FEValuesBase< dim > &	int_fe,
		double	friction_coefficient
	)

One-sided friction interface term as it shows up in the Beavers-Joseph-Saffman condition.