NuTBox

I wrote a little couple of programs, respectively in Fortran and C++, as a proof of concept for connecting a Fortran program to a sort of visualization server based on VTK. The nice thing is that it uses MPI for the connection, so on the Fortran side nothing new and scary.

The code (you can find it at https://github.com/plampite/vtkForMPI) and the idea strongly predate a similar example in Using Advanced MPI by W. Gropp et al., but makes it more concrete by adding actual visualization...

In previous posts of this series I presented an elaboration of the Musker-Monkewitz analytical wall function that allowed extensions to non equilibrium cases and to thermal (scalar) cases with, in theory, arbitrary Pr/Pr_t (Sc/Sc_t) ratios.

In the meanwhile, I worked on a rationalization and generalization of the framework, derivation of averaged production term for the TKE equation, etc.

While the new material is presented in a substantially different manner and will...

In this post i summarize the initial problem and the procedure to determine the wall function value (i.e., the solution) for given ,, and .

We looked for a solution to the problem:



with:

...

In the first part of this post we left with the problem of computing the following integral:



with:



I added all the explicit functional dependencies here because we know ...

In a previous post i wrote about an extension of the Reichardt law of the wall to pressure gradient effects. That was derived by assuming the Reichardt profile as a solution for the case without pressure gradient and using integration by parts. In particular, given the the Reichardt function ( is the Von Karman constant):



with:
...