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 12: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 95 Comments 0 sbaffini is offline Edit Tags
Old

A generalized thermal/dynamic wall function: Part 4

Posted February 7, 2019 at 10: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, 100 views)
File Type: txt musker.txt (930 Bytes, 112 views)
sbaffini's Avatar
Senior Member
Posted in Uncategorized
Views 79 Comments 0 sbaffini is offline Edit Tags
Old

A generalized thermal/dynamic wall function: Part 3

Posted October 17, 2016 at 11:25 by sbaffini (NuTBox)
Updated November 18, 2018 at 05: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, 170 views)
sbaffini's Avatar
Senior Member
Views 751 Comments 0 sbaffini is offline Edit Tags
Old

A generalized thermal/dynamic wall function: Part 2

Posted October 17, 2016 at 08:24 by sbaffini (NuTBox)
Updated December 21, 2016 at 09: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 743 Comments 0 sbaffini is offline Edit Tags
Old

A generalized thermal/dynamic wall function: Part 1

Posted October 14, 2016 at 12:27 by sbaffini (NuTBox)
Updated December 21, 2016 at 09: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 1019 Comments 0 sbaffini is offline Edit Tags

All times are GMT -4. The time now is 18:53.