myprojects-hg-reborn

Fork

I updated the file VELA/CVS.mph by fixing a problem with the temperature
close to the layer.

I corrected a small bug in the code R-codes/read_two_sets.R (before
lendat was appearing twice unnecessarily).

I added the file VELA/CVS.mph which tries modelling the CVS by using a
mesh with bricks and a boundary layer.

I added the file VELA/anaconda.mph which contains the modelling with
Comsol of the anaconda (mesh with bricks and the use of wall functions
for the velocity)
)

I added the file VELA/test-grid-square2-bis+anaconda.mph, which contains
the results of the latest vela simulations. It has both the anaconda and
the CVS. It contains a boundary layer both along the anaconda and along
the CVS. I maninly used bricks mesh the geometry, but I also used a mesh
with triangular elements for some domains.

I added the script R-codes/read_two_sets.R, which automatically reads
and plots 2 datasets of data to ease the comparison.
I also updated the readme file

I modified the code R-codes/Housiadas_T_profile.R by simply adding a
legend to the plot of the two theoretical temperature profiles (the one
by Christos and the one obtained using the Diffus-Boelter correlation)

I added a new directory which contains a simple calculation of the
mean axial temperature profile in the anaconda

I updated the code R-codes/Housiadas_T_profile.R, which now also
computes the temperature profile using the Dittus-Boelter correlation

I updated the file R-codes/Housiadas_T_profile.R correcting also a few
bugs.
Now, it calculates the Re from the physics of the flow in a pipe and it
also selects the appropriate power law for the velocity profile
according to the Re.

I added a few lines to the code R-codes/Housiadas_T_profile.R.
It is now able to plot read the radial velocity profiles and compare
them to the one-seventh law or similar empirical expressions.

I added some parts to R-codes/Housiadas_T_profile.R: now the code is
able to read & integrate some radial temperature profile, integrate them
and get the numerical dimensionless temperature

I added the small code R-codes/Housiadas_T_profile.R, which calculates
& prints out the temperature profile in a tube using the turbulent
expression in the paper by Housiadas & Yannis & I updated the readme
file

I updated the code R-codes/readme after adding the file R-vortex-inlet.R

I added the code R-codes/R-vortex-inlet.R. It is a very short and simple
script which plots the initial velocity profile at the inlet given by
eq. 7.9 in the book by Orlandi.

As usual, I created my own Makefile to compile the project and I made
some tiny modifications to the PIP/pipma.f file (I simply fixed the
paths leading to some input files which have to be read by the code).
The code now compiles and seems to work, but I did not wait till the
completion of the simulation ---> careful! There may be problems later
on!

I created a new Makefile to be able to compile the code and run it.
I had to slightly modify the codes TR3D/tr3dnuma.f and
TR3D/tr3dnuinri.f. I basically had to trivially change the path leading
to some files which had to be read and the format to save some output
(or better: to print it out), since the compiler used to complain.
The code now seems to be running but I have not completed a single run
yet---> careful, there may be problems!

I added the R code RK/read_burgers.R which reads and plots some of the
output files from the burgts.f code and calculates numerically the first
derivative of the velocity profile wrt the physical coordinate according
to eq. 2.14 in the book by Orlandi

I modified the file RK/burgts.f by adding some comments and slighting
changing the output format in order to make it easier to do
post-processing with an R code.

I added the R file ELLIP/read_stream.R which reads the analytical and
numerical solution to Poisson equation and plots then for some selected
values.

I modified files ELLIP/param.f and ELLIP/capro.f in order to be able to
save the numerical results in a convenient way for post-processing them
with R. I add to modify the file param.f in order to be able to add some
extra arrays to save the data.
NB: Orlandi uses implicit declarations for his arrays!!!!!
Be very careful with that!!!!!

I modified the file ELLIP/capro.f by adding some comments which are
supposed to make the code more readable

I modified the files ELLIP/capro.f and ELLIP/Makefile essentially by
removing the parts dealing with the p2fftess.f file which uses a library
which is not available on the system.

Trivial modification to Makefile (basically I am now compiling with
gfortran-4.1 instead of f95). It goes without saying that I can create
the .o files in directory NCARFFT by inputing gfortran-4.1 -c filename

I did a few things, most of which are described in the readme file I
added.
1) I actually added a readme file. This is a modification of one of the
two makefiles initially provided by Orlandi in this directory.
2)The Makefile now requires that the .o files are preliminarly created
in the NCARFFT folder. (just run f95 -c filename for both files and you
are done)
3)The file p2fftess needed some minor modification (fundamentally I had
to take away the lines invoking some specific library available on
IBM computers)
4)On can use [so far] both f95 and gfortran-4.1 commands to compile the
code. Perhaps in the future it will be better to use only the gfortran
compiler [the other one does not seem to be supported any longer]

Some modifications again to SECDBL/r_code_secdbls_f.R and
SECDBL/secdbls.f. Nothing substantial but it helps polishing a few
details.

I added the R code SECDBL/r_code_secdbls_f.R to be read and plot some of
the output files of secdbls.f (this is still preliminary)

I modified the file SECDBL/secdbls.f basically by changing the way in
which the output files are saved (to make it easier to post-process them
with an R code)

I modified SECD/r_code_secdml_f.R by adding a few lines to plot also the
velocity profile saved in file u_profile

I modified the code SECD/secdml.f by changing the names and formats of
the output files. Furthermore, I added an R code,
SECD/r_code_secdml_f.R, which reads and plots out the quantities of
interest