Skip to content
Snippets Groups Projects
Commit 127f33a9 authored by Ned Coltman's avatar Ned Coltman
Browse files

Update README.md

parent 50c7f8a0
No related branches found
No related tags found
No related merge requests found
......@@ -54,13 +54,22 @@ struct LocalResidual<TypeTag, TTag::TwoPIncompressible>
using type = MyLocalResidual<TypeTag>;
};
```
We are hoping to evaluate the following function:
```math
\frac{\partial (\phi S_\alpha)}{\partial t} + \nabla \cdot \textbf{v}_\alpha = \frac{q_\alpha}{\rho_\alpha}
```
```math
\textbf{v}_\alpha = -\lambda_\alpha \textbf{K} \left( \nabla p_\alpha - \rho_\alpha \textbf{g}\right)
```
#### Task 2.1 computeFlux Method
First, look at the `computeFlux` method in `mylocalresidual.hh`. Within the lambda function we need to return the upwinded term. As this is an incompressible model, no density term is required, only the mobility term. This value is stored in the volume variables object (volVars). Complete the lambda function.
First, look at the `computeFlux` method in `mylocalresidual.hh`.
Within the lambda function we need to return the upwinded term. As this is an incompressible model, no density term is required, only the mobility term. This value is stored in the volume variables object (volVars). Complete the lambda function.
* complete the upwind term lambda in the `computeFlux` method.
#### Task 2.2 computeStorage Method
Second, complete `computeStorage` method. This term, when incompressible, looks like this:
```math
\frac{\partial (\phi S_\alpha)}{\partial t}
```
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment