Blog Entries

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...

Direct Numerical Simulation on a wing profile

1 billion points DNS (Direct Numerical Simulation) on a NACA4412 profile at 5 degrees angle of attack. Reynolds number is 350000 per airfoil chord and Mach number is 0.117. Both upper and lower turbulent boundary layers are tripped respectively at 15% and 50% by roughness elements evenly spaced in the boundary layer created by a zonal immersed boundary condition (Journal of Computational Physics, Volume 363, 15 June 2018, Pages 231-255,...

https://github.com/mcavallerin/airFoil_tools


generate 3D model for foils

Hello,
I am new to OF and gnuplot. I am working on the T3A flat plate case in tutorials of OpenFOAM. I was struggling a lot to plot graph of friction coefficient from simulation and experimental data using the default plot_file (creatGraphs.plt) provided with the tutorial. I looked on the internet for a solution but remained unsuccessful. So, after trying for some time, I got the graph correct so I decided to share it here for use of anyone else.

The trick is we have to edit...

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...