Anuncio

Colapsar
No hay ningún anuncio todavía.

Kutta JOUKOWSKI

Colapsar
X
 
  • Filtro
  • Hora
  • Mostrar
Borrar todo
nuevos mensajes

  • Kutta JOUKOWSKI

    Hola!!!!!!!!
    Espero y esten super bien les quisiera pedir un poquito de ayuda en un programa que estoy haciendo en python estoy mapeando en el plano complejo , para hacer la simulacion del ala de un avion miren aqui encontre el codigo en matlab veanlo
    % FLUID PARAMETER
    rho = 1.225;
    % TRANSFORMATION PARAMETER
    lambda = r-s;
    % CIRCULATION
    beta = (theta);
    k = 2*r*v*sin(beta);
    Gamma = k/(2*pi); %CIRCULATION

    %COMPLEX ASYMPTOTIC SPEED
    w = v * exp(i*theta);

    %TOLLERANCE
    toll = +5e-2;

    % GENERATING MESH
    x = meshgrid(-5:.1:5);
    y = x';

    % COMPLEX PLANE
    z = x + i*y;

    % Inside-circle points are Excluded!
    for a = 1:length(x)
    for b = 1:length(y)
    if abs(z(a,b)-s) <= r - toll
    z(a,b) = NaN;
    end
    end
    end

    % AERODYNAMIC POTENTIAL
    f = w*(z) + (v*exp(-i*theta)*r^2)./(z-s) + i*k*log(z);

    % JOUKOWSKI TRANSFORMATION,
    J = z+lambda^2./z;

    %GRAPHIC - Circle and Joukowski Airfoil
    angle = 0:.1:2*pi;
    z_circle = r*(cos(angle)+i*sin(angle)) + s;
    z_airfoil = z_circle+lambda^2./z_circle;

    % KUTTA JOUKOWSKI THEOREM
    L = v_inf*rho*Gamma;
    L_str = num2str(L);


    Pero en esta parte donde usan la transpuesta segun matlab no se como hacerlo en python vean % AERODYNAMIC POTENTIAL
    f = w*(z) + (v*exp(-i*theta)*r^2)./(z-s) + i*k*log(z);
    se supone q es la transpuesta no ?
    Espero y alguien me pueda ayudar
    Gracias
    Saludos
    Naoko Konichiiwa!!!!!!!!!

Contenido relacionado

Colapsar

Trabajando...
X