Example #1
0
def main():
    # Input
    epsilon = 0.025
    tau = 30
    alpha = 10
    
    # Generate coordinates and cp distribution
    x, y, cp = van_de_vooren(tau, epsilon, alpha)
    
    # Plot the pressure distribution
    plot_cp(x, cp, x, y)
Example #2
0
def main():
    # This function runs the vortex panel method on a NACA 2412
    
    # Generate airfoil coordinates
    airfoil_coordinates = NACA4('2412', 50)
        
    x_airfoil = airfoil_coordinates[:, 0]
    y_airfoil = airfoil_coordinates[:, 1]
    
    # Set the angle of attack
    alpha = 10

    # Run the vortex panel method
    x, cp, cl, cm = vortex_panel_method(alpha, airfoil_coordinates)
    
    # Print the values of cl and cm in the terminal 
    print('C_l = ', cl)
    print('C_m = ', cm)
    
    # Plot the pressure distribution
    plot_cp(x, cp, x_airfoil, y_airfoil)
    
    # Plot the lift coefficient and moment coefficient versus angle of attack
    plot_polar(airfoil_coordinates)
Example #3
0
rtstruct = sb.build_structure_set()
rtstruct.add_external_box()
sph = rtstruct.add_sphere(radius=70,
                          center=[-50, 0, -100],
                          name='Sph-Organ',
                          interpreted_type='CAVITY')
sph2 = rtstruct.add_sphere(radius=25,
                           center=[0, 0, 0],
                           name='Sph-PTV',
                           interpreted_type='PTV')

rtplan = sb.build_static_plan()
b1 = rtplan.build_beam(gantry_angle=180, collimator_angle=15, meterset=100)
b1.conform_to_rectangle(1, 1, [0, 0])

rtplan.build(finalize_mlc=False)

modules.conform_mlc_to_roi(b1.rtbeam, sph.roi_contour, sb.current_study)
modules.conform_mlc_to_roi(b1.rtbeam, sph2.roi_contour, sb.current_study)
b1.finalize_mlc()

sb.write(tmpdir)
print tmpdir

import plotting as p

p.plot_cp(rtplan.datasets[0].Beams[0], rtplan.datasets[0].Beams[0].CPs[0])
p.plot_roi_in_cp(rtplan.datasets[0].Beams[0],
                 rtplan.datasets[0].Beams[0].CPs[0], sph.roi_contour,
                 sb.current_study)
pp.axis('image')
Example #4
0
print ct.pixel_array.max(),ct.pixel_array.min()
ct.add_sphere(radius = 25, center = [0,0,0], real_value = -1000, mode = 'set')
print ct.pixel_array.max(),ct.pixel_array.min()
ct.add_box(size = [25,50,25], center = [0,0,0], stored_value = 300, mode = 'add')
print ct.pixel_array.max(),ct.pixel_array.min()

assert sb.seriesbuilders['CT'] == [ct]

rtstruct = sb.build_structure_set()
rtstruct.add_external_box()
sph = rtstruct.add_sphere(radius=70, center = [-50,0,-100], name='Sph-Organ', interpreted_type='CAVITY')
sph2 = rtstruct.add_sphere(radius = 25, center = [0,0,0], name='Sph-PTV', interpreted_type='PTV')

rtplan = sb.build_static_plan()
b1 = rtplan.build_beam(gantry_angle = 180, collimator_angle = 15, meterset = 100)
b1.conform_to_rectangle(1,1,[0,0])

rtplan.build(finalize_mlc = False)

modules.conform_mlc_to_roi(b1.rtbeam, sph.roi_contour, sb.current_study)
modules.conform_mlc_to_roi(b1.rtbeam, sph2.roi_contour, sb.current_study)
b1.finalize_mlc()

sb.write(tmpdir)
print tmpdir

import plotting as p

p.plot_cp(rtplan.datasets[0].Beams[0], rtplan.datasets[0].Beams[0].CPs[0])
p.plot_roi_in_cp(rtplan.datasets[0].Beams[0], rtplan.datasets[0].Beams[0].CPs[0], sph.roi_contour, sb.current_study)
pp.axis('image')