CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums

NuTBox

Register Blogs Members List Search Today's Posts Mark Forums Read

Old

Connecting Fortran with VTK - the MPI way

Posted May 24, 2019 at 13:12 by sbaffini (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...
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 271 Comments 0 sbaffini is online now Edit Tags
Old

A generalized thermal/dynamic wall function: Part 4

Posted February 7, 2019 at 11:29 by sbaffini (NuTBox)

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...
Attached Files
File Type: txt comparewf.txt (2.0 KB, 155 views)
File Type: txt musker.txt (930 Bytes, 165 views)
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 147 Comments 0 sbaffini is online now Edit Tags
Old

A generalized thermal/dynamic wall function: Part 3

Posted October 17, 2016 at 12:25 by sbaffini (NuTBox)
Updated November 18, 2018 at 06:57 by sbaffini

In this post i summarize the initial problem and the procedure to determine the wall function value (i.e., the solution) for given y^+,F_T^+,Pr and Pr_t.

We looked for a solution T^+\left(y^+,F_T^+,Pr,Pr_t\right) to the problem:

\frac{dT^+}{dy^+}=\frac{Pr\left(1+F_T^+y^+\right)}{\left[1+\left(\frac{Pr}{Pr_t}\right)\left(\frac{\mu_t}{\mu}\right)\right]}

with:

\frac{\mu_t}{\mu}=\frac{\left(ky^+\right)^3}{\left(ky^+\right)^2+\left(ka_0\right)^3-\left(ka_0\right)^2}...
Attached Files
File Type: txt wallfn.txt (2.2 KB, 179 views)
sbaffini's Avatar
Senior Member
Views 804 Comments 0 sbaffini is online now Edit Tags
Old

A generalized thermal/dynamic wall function: Part 2

Posted October 17, 2016 at 09:24 by sbaffini (NuTBox)
Updated December 21, 2016 at 10:07 by sbaffini

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

f^+\left(y^+,\frac{Pr}{Pr_t}\right)=\int_0^{y+}{\frac{1}{\left[1+\left(\frac{Pr}{Pr_t}\right)\left(\frac{\mu_t}{\mu}\right)\right]}dz^+}

with:

\frac{\mu_t}{\mu}\left(y^+,a,k\right)=\frac{\left(ky^+\right)^3}{\left(ky^+\right)^2+\left(ka\right)^3-\left(ka\right)^2}

I added all the explicit functional dependencies here because we know f^+\left(y^+,1\right)...
sbaffini's Avatar
Senior Member
Views 799 Comments 0 sbaffini is online now Edit Tags
Old

A generalized thermal/dynamic wall function: Part 1

Posted October 14, 2016 at 13:27 by sbaffini (NuTBox)
Updated December 21, 2016 at 10:07 by sbaffini

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 (k is the Von Karman constant):

f^+\left(y^+\right) = \frac{1}{k}\log\left(1+ky^+\right) +A\left(1-e^{-\frac{y^+}{B}}-\frac{y^+}{B}e^{-\frac{y^+}{C}}\right)

with:
...
sbaffini's Avatar
Senior Member
Views 1070 Comments 0 sbaffini is online now Edit Tags

All times are GMT -4. The time now is 10:01.