dd6abb687427149075d0cf7c608b7b64fc551115,trunk/SUAVE/Methods/Aerodynamics/Common/Fidelity_Zero/Lift/weissinger_vortex_lattice.py,,weissinger_vortex_lattice,#,19

Before Change

        
        ya = np.atleast_2d((i)*deltax)
        yb = np.atleast_2d((i+1)*deltax)
        xa = np.atleast_2d(((i+1)*deltax-deltax/2)*np.tan(sweep) + 0.25*sl)
        x  = np.atleast_2d(((i+1)*deltax-deltax/2)*np.tan(sweep) + 0.75*sl)
        y  = np.atleast_2d(((i+1)*deltax-deltax/2))
        
        xloc_leading  = ((i+1)*deltax)*np.tan(sweep)
        xloc_trailing = ((i+1)*deltax)*np.tan(sweep) + sl        
                
        RHS  = np.atleast_2d(np.sin(twist_distri+aoa)).T
        A = (whav(x,y,xa.T,ya.T)-whav(x,y,xa.T,yb.T)\
            -whav(x,y,xa.T,-ya.T)+whav(x,y,xa.T,-yb.T))*0.25/np.pi
    
        // Vortex strength computation by matrix inversion
        T = np.linalg.solve(A.T,RHS)
        
        // Calculating the effective velocty         
        A_v = A*0.25/np.pi*T
        v   = np.sum(A_v,axis=0)
        
        Lfi = -T.T * (np.sin(twist_tc)-v)
        Lfk =  T.T * np.cos(twist_tc)   
        Lft = -Lfi*np.sin(twist_tc)+Lfk*np.cos(twist_tc)
        Dg  = Lfi*np.cos(twist_tc)+Lfk*np.sin(twist_tc)
            
        L  = deltax * Lft
        D  = deltax * Dg
        
        // Total lift
        LT = np.sum(L)
        DT = np.sum(D)
    
        Cl = 2*LT/(0.5*Sref)

After Change

        ya = np.atleast_2d((i)*deltax)
        yb = np.atleast_2d((i+1)*deltax)
        xa = np.atleast_2d(((i+1)*deltax-deltax/2)*np.tan(sweep) + 0.25*sl)
        x  = np.atleast_2d(((i+1)*deltax-deltax/2)*np.tan(sweep) + 0.75*sl)
        y  = np.atleast_2d(((i+1)*deltax-deltax/2))
        
        xloc_leading  = ((i+1)*deltax)*np.tan(sweep)
        xloc_trailing = ((i+1)*deltax)*np.tan(sweep) + sl        
                
        RHS  = np.atleast_2d(np.sin(twist_distri+aoa)).T
        A = (whav(x,y,xa.T,ya.T)-whav(x,y,xa.T,yb.T)\
            -whav(x,y,xa.T,-ya.T)+whav(x,y,xa.T,-yb.T))*0.25/np.pi
    
        // Vortex strength computation by matrix inversion
        T = np.linalg.solve(A.T,RHS)
        
        // Calculating the effective velocty         
        A_v = A*0.25/np.pi*T
        v   = np.sum(A_v,axis=0)
        
        Lfi = -T.T * (np.sin(twist_tc)-v)
        Lfk =  T.T * np.cos(twist_tc)   
        Lft = -Lfi*np.sin(twist_tc)+Lfk*np.cos(twist_tc)
        Dg  = Lfi*np.cos(twist_tc)+Lfk*np.sin(twist_tc)
            
        L  = deltax * Lft
        D  = deltax * Dg
        
        // Total lift
        LT = np.sum(L)
        DT = np.sum(D)
    
        Cl = 2*LT/(0.5*Sref)

Instances