runellipductdisk.m

clear all
format compact
close all

% Freestream velocity
W = 1.0;
Minf = 0.0;

Sref = 1.0;

xlim = [-2 1];
rlim = [0 2];
nsurvey = 51;

nsl = 21;

drawplots = true;

ntstep = 6;

% Construct an ellipse
ncirc = 136;
rada = 1.0;
radb = .25;
xcen = 0.0;

% Construct a duct
naf = 151;

chord = 0.6;
alpha = 15;       % Degrees about LE
xoff = 0.7;      % X-Offset of LE
roff = 0.4;      % R-Offset of LE

% NACA 4-Digit airfoil parameters
dig1 = 4;
dig2 = 4;
dig34 = 12;

flipaf = false;

% Airfoil point spacing
% 1 -- Cosine cluster LE, TE
% 2 -- Cosine cluster LE
% 3 -- Uniform
spacing = 1;


% Actually generate bodies
[xep, rep, rad, Vex] = setupellipsoid( ncirc, rada, radb, xcen );

jtelow = -1;
jteup = - 1;
kutta = false;

names{1} = 'Ellipsoid';
xepts{1} = xep;
repts{1} = rep;
kuttas{1} = kutta;
props{1} = false;
deltaCP{1} = 0.0;
jtels{1} = jtelow;
jteus{1} = jteup;
rads{1}=rad;
Vexs{1}=Vex;


[xep, rep, rad, Vex] = setupNACAduct( naf, chord, alpha, xoff, roff, dig1, dig2, dig34, flipaf, spacing );

% Index for Kutta condition
jtelow = 1;
jteup = naf - 1;
kutta = true;


names{2} = ['NACA ' num2str(dig1) num2str(dig2) num2str(dig34) ' Duct'];
xepts{2} = xep;
repts{2} = rep;
kuttas{2} = kutta;
props{2} = false;
deltaCP{2} = 0.0;
jtels{2} = jtelow;
jteus{2} = jteup;
rads{2}=rad;
Vexs{2}=Vex;

% Actuator disk
dCP = .75;

% Disk radius
rdisk = 1.0;

rdisk = repts{2}(1);
xstart = xepts{2}(1);

% End of contracting streamtube
xend = 5;

% Vortex ring spacing
dxring = rdisk * 0.1;

% Number of vortex ring panels
npan = ( xend - xstart ) / dxring;

% Set up initial geometry
xpts = linspace( xstart, xend, npan + 1 );
rpts = rdisk * ones( size( xpts ) );

names{3} = 'Disk';
xepts{3} = xpts;
repts{3} = rpts;
kuttas{3} = false;
props{3} = true;
deltaCP{3} = dCP;
xdisk{3} = 0.85;
jtels{3} = 0;
jteus{3} = 0;
rads{3} = 0;
Vexs{3} = nan;

% Execute script
run('bor')