kohyun / contur

CONTUR code for contoured nozzle design by J.C. Sivells ADEC TR 78 63 December 1978

# contur
CONTUR code for contoured nozzle design by J.C. Sivells ADEC TR 78 63 December 1978

This respository contains the CONTUR code in 
Appendix D of J.C. Sivells, A computer program for the aerodynamic
design of axisymmetric and planar nozzles for supersonic and hypersonic
wind tunnels, ARO Inc., a Sverdrup Corporation Company, ADEC TR 78 63,
December 1978.

The code uses a combination of analytical solutions, the method of 
characteristics, and centerline distributions in order to calculate the
divergent section of a convergent-divergent de Laval nozzle.
 
Folder sivells contains a FORTRAN77 source code which runs through 
a series of 16 subroutines and uses 7 user-defined input cards
describing the flow conditions of the desired nozzle profile.

Folder src contains the FORTRAN90 source code version of folder sivells

Sivells reports the input cards and output file of a MACH 4 
axisymmetric nozzle. The version of the code in this repository was
run in May 2019 on this test case and the same output was obtained.

Florentina-Luiza Zavalan, under the guide of project
supervisor Aldo Rona, typeset Appendix D from the freely available
source code listing at 
https://apps.dtic.mil/dtic/tr/fulltext/u2/a062944.pdf
Significant effort was put in to interpret the low-resolution scan and
to disambigue similar typographical symbols (e.g. * from +, 5 from S).

Aldo Rona translated the code in FORTRAN90. If you find any bug
while testing, please report it to the authors: 
[e-mail], [e-mail]

License: See LICENSE.txt. Code users must acknowledge the
provenance of the listing by using the following acknowledgement in
their published work: "This work
used the CONTUR source code by L.F. Zavalan and A. Rona, based on the
computer program by J.C. Sivells".

Compilation instructions

Linux/Unix users:
1. Create a new directory, suggested name: sivells
2. Download in the SAME directory the bundle of 20 *.f source files,
   one input.txt file, and makefile.
3. You are advised to also download Sivells.pdf, which is article by
   Sivells, and output.txt, which is the sample output from the MACH 4 
   test case. These files are in folders sivells/ and docs/.
4. Compile the code by just typing: make. This will create the object
   files *.o and the executable MAIN.exe
5. Run the executable MAIN.exe

Windows users:
1. Download and install Microsoft Visual Studio
2. Open Microsoft Visual Studio and create a new Project, suggested name: contur
3. Move input.txt and input2d.txt from the /src directory into the
   'resources' directory of the Project.
4. Move all .f files from the /src directory to the 'source' directory of the Project.
   Make sure you only move the .f files from /src, which have all small case names, do
   not add any of the .f files from the /sivells directory, which have CAPITAL CASE names,
   as this will over-define symbols in the compilation.
5. Right click on the makefile you have downloaded and open it with a text editor.
   Do not double click on it. Check the list of files you now have in 'source' 
   of your Project. You should only have files with root names <root>.f matching 
   the object list root names <root>.o in makefile, e.g. axial.f -> axial.o. You should
   also have main.f in your 'source'. You should not 
   have any .txt, makefile, or .pdf files in your 'source'.
6. Remove makefile. You do not need this file for compiling under Windows.
7. Select Compile and Run from the taskbar of Microsoft Visual Studio

Frequently Asked Questions
1. You may need to let makefile know what compiler you are using. To
do so, use a text editor to edit makefile and change FC=ifort to your
own compiler, e.g. gfort. Run the command "man -k fortran" to find out
what compiler is installed on your system.

2. You may need to change the compiler flags according to what is
available from your compiler. In unix/linux, you can type
"man mycompiler" to find out what flags are available.

08 May 2019 - A. Rona, F.L. Zavalan

Card notes:

Card 5: Use an integer value for XJ as JX=XJ and JX is integer

Code updates:
14May2019: AR moved label 5 from line BOU 72 to BOU 69 to avoid jumping
into the DO IV=1,IW loop.
16May2019: AR variables ns and nc removed from AXIAL.f as not used
16May2019: AR removed M from common block of NEO.f
11Jun2019: AR code changed to f90. Added input2d.txt as sample input
file to generate a two-dimensional contoured nozzle. To use rename
input2d.txt as input.txt and rund the code

-----
end of notes