This routine adds Brinkman (porous-media) damping to the velocity stored in the auxField for the incompressible fluid model. Reference: Dominique d’Humières, Irina Ginzburg, (2009), "Viscosity independent numerical errors for Lattice Boltzmann models: From recurrence equations to “magic” collision numbers", Computers & Mathematics with Applications, Volume 58, Issue 5, Pages 823-840
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(mus_source_op_type), | intent(inout) | :: | fun |
Description of method to update source |
||
| real(kind=rk), | intent(inout) | :: | auxField(:) |
output auxField array |
||
| integer, | intent(in) | :: | iLevel |
current level |
||
| type(tem_time_type), | intent(in) | :: | time |
current timing information |
||
| type(tem_varSys_type), | intent(in) | :: | varSys |
variable system definition |
||
| type(mus_convertFac_type), | intent(in) | :: | phyConvFac |
Physics conversion factor for current level |
||
| type(mus_derVarPos_type), | intent(in) | :: | derVarPos(:) |
position of derived quantities in varsys |