Scribd Logo
Upload

Welcome to Scribd!

  • 17 views

    All All: %nodos

    Uploaded by

    Ronaldo Brioso
    This document contains code to simulate fluid flow through a channel with two inlets of different flow rates (Q1 and Q2) using the finite difference method. It defines parameters like channel length and width, discretizes the geometry into a grid, sets boundary conditions, and iterates a solution over time steps to calculate flow velocities throughout the channel. The results are plotted with contour and quiver plots to visualize the velocity profile.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd

    All All: %nodos

    Uploaded by

    Ronaldo Brioso
    17 views5 pages
    This document contains code to simulate fluid flow through a channel with two inlets of different flow rates (Q1 and Q2) using the finite difference method. It defines parameters like channel length and width, discretizes the geometry into a grid, sets boundary conditions, and iterates a solution over time steps to calculate flow velocities throughout the channel. The results are plotted with contour and quiver plots to visualize the velocity profile.

    Original Title

    Brioso

    Copyright

    © © All Rights Reserved

    Available Formats

    DOCX, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    This document contains code to simulate fluid flow through a channel with two inlets of different flow rates (Q1 and Q2) using the finite difference method. It defines parameters like channel length and width, discretizes the geometry into a grid, sets boundary conditions, and iterates a solution over time steps to calculate flow velocities throughout the channel. The results are plotted with contour and quiver plots to visualize the velocity profile.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    17 views5 pages

    All All: %nodos

    Uploaded by

    Ronaldo Brioso
    This document contains code to simulate fluid flow through a channel with two inlets of different flow rates (Q1 and Q2) using the finite difference method. It defines parameters like channel length and width, discretizes the geometry into a grid, sets boundary conditions, and iterates a solution over time steps to calculate flow velocities throughout the channel. The results are plotted with contour and quiver plots to visualize the velocity profile.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    You are on page 1of 5

    Programa 2

    clear all
    close all
    L=300;
    W1=50;
    W2=10;
    h=5;

    Q1=10;
    Uoo1=Q1/(W1*h);
    q1=Q1/W1;

    Q2=2;
    Uoo2=Q2/(W2*h);
    q2=Q2/W2;

    %Nodos
    dx=1; nx=L/dx+1;
    dy=1; ny=W1/dy+1;
    alpha=(dx/dy)^2;
    Lei=100; ie1=(Lei/dx)+1;
    Lef=110; ie2=(Lef/dx)+1;

    %Malla
    for j=1:ny
    for i=1:nx
    x(i,j)=(i-1)*dx;
    y(i,j)=(j-1)*dy;
    end
    end

    for j=1:ny
    for i=1:nx
    x(i,j)=(i-1)*dx;
    y(i,j)=(j-1)*dy;
    end
    end

    %Funcion Corriente A

    A(1:nx,1:ny,1)=0;

    kmax=1000;

    for k=1:kmax
    for j=2:ny-1
    for i=2:nx-1
    A(i,j,k+1)=((A(i+1,j,k)+A(i-1,j,k)+alpha*(A(i,j+1,k)+A(i,j-1,k)))/
    (2*(1+alpha)));
    end
    end

    % condiciones de frontera
    % Arriba y Abajo
    A(1:nx,ny,k+1)=q1;
    A(1:ie1,1,k+1)=0;
    A(ie2:nx,1,k+1)=q2;

    % Entrada y Salida
    for j=1:ny
    A(1,j,k+1)=(j-1)*q1/(ny-1);
    A(nx,j,k+1)=((j-1)*(q1-q2))/(ny-1);
    end
    % Entrada del segundo caudal
    for i=ie1:ie2
    A(i,1,k+1)=(i-ie1)*(q2)/(ie2-ie1);
    end

    end

    %Velocidades
    u(1:nx,1:ny)=0;
    v(1:nx,1:ny)=0;

    for j=2:ny
    for i=2:nx
    u(i,j)=(A(i,j,k)-A(i,j-1,k))/dy;
    v(i,j)=-(A(i,j,k)-A(i-1,j,k))/dx;
    end
    end

    figure
    contour(x,y,A(:,:,kmax-1),20)
    hold on
    quiver(x,y,u,v,1)
    axis equal
    Problema 2
    clear all
    close all
    L=300;
    W1=50;
    W2=10;
    h=5;

    Q1=10;
    Uoo1=Q1/(W1*h);
    q1=Q1/W1;

    Q2=2;
    Uoo2=Q2/(W2*h);
    q2=Q2/W2;

    %
    dx=1; nx=L/dx+1;
    dy=1; ny=W1/dy+1;
    alpha=(dx/dy)^2;
    Lei=100; ie1=(Lei/dx)+1;
    Lef=110; ie2=(Lef/dx)+1;

    % definiendo la geometria
    for j=1:ny
    for i=1:nx
    x(i,j)=(i-1)*dx;
    y(i,j)=(j-1)*dy;
    end
    end

    for j=1:ny
    for i=1:nx
    x(i,j)=(i-1)*dx;
    y(i,j)=(j-1)*dy;
    end
    end

    A(1:nx,1:ny,1)=0;

    kmax=1500;

    for k=1:kmax
    for j=2:ny-1
    for i=2:nx-1
    A(i,j,k+1)=((A(i+1,j,k)+A(i-1,j,k)+alpha*(A(i,j+1,k)+A(i,j-1,k)))/
    (2*(1+alpha)));
    end
    end

    % condiciones de frontera

    %arriba y abajo
    A(1:nx,ny,k+1)=q1;
    A(1:ie1,1,k+1)=0;
    A(ie2:nx,1,k+1)=q2;

    % entrada y salida
    for j=1:ny
    A(1,j,k+1)=(j-1)*q1/(ny-1);
    A(nx,j,k+1)=q1+((j-1)*(-q2))/(ny-1);
    end
    % en la entrada del segundo caudal
    for i=ie1:ie2
    A(i,1,k+1)=(i-ie1)*(q2)/(ie2-ie1);
    end

    end

    % calculando las velocidades


    u(1:nx,1:ny)=0;
    v(1:nx,1:ny)=0;

    for j=2:ny
    for i=2:nx
    u(i,j)=(A(i,j,k)-A(i,j-1,k))/dy;
    v(i,j)=-(A(i,j,k)-A(i-1,j,k))/dx;
    end
    end

    figure
    contour(x,y,A(:,:,kmax-1),20)
    hold on
    quiver(x,y,u,v,1)
    axis equal

    You might also like