Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 6 vistas

    Ejercicio 1 Jesus Saucedo Lopez

    Cargado por

    e.jsaucedolopez
    Solución de ejercicio de modelos numéricos

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Ejercicio 1 Jesus Saucedo Lopez

    Cargado por

    e.jsaucedolopez
    6 vistas11 páginas
    Solución de ejercicio de modelos numéricos

    Título original

    EJERCICIO_1_JESUS_SAUCEDO_LOPEZ

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Solución de ejercicio de modelos numéricos

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    6 vistas11 páginas

    Ejercicio 1 Jesus Saucedo Lopez

    Cargado por

    e.jsaucedolopez
    Solución de ejercicio de modelos numéricos

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 11

    Resolución del Problema de Dirichlet

    Máster en Física: Radiaciones, Nanotecnología, Partículas y Astrofísica.

    Asignatura: Complementos Matemáticos y Numéricos

    Catedráticos: PhD. María Gómez Rocha

    Estudiante: Jesús Saucedo López.

    Granada, España a 07 de febrero del 2024.

    1
    2
    Resolución del ejercicio
    El problema de Dirichlet representa un tipo de desafío en la resolución de ecuaciones
    diferenciales parciales, donde el objetivo radica en descubrir la solución de una ecuación dentro
    de un área específica, considerando valores definidos en su perímetro. Este problema se
    manifiesta en múltiples campos, abarcando la física, ingeniería y matemática aplicada, y su
    resolución resulta fundamental para comprender cómo se desarrollan fenómenos físicos dentro
    de un espacio delimitado.
    Método de ejecución:
    Lo que he realizado fue aplicar el método de diferencias finitas que segmenta el dominio Ω
    mediante una rejilla, separando el espacio en nodos distribuidos de manera uniforme. Luego,
    aproximamos las derivadas parciales en la ecuación de Laplace utilizando diferencias finitas, lo
    que nos permitirá construir un sistema de ecuaciones lineales.
    Primera parte correspondiente a la primera sugerencia del enunciado:
    De esta forma , examinaremos el problema de Dirichlet asociado a la ecuación de Laplace en un
    dominio bidimensional Ω, el cual se define como {(x, y), 0 < x < 1, 0 < y < 1}. La ecuación a
    abordar es ∇2u(x, y) = 0 dentro de Ω, con condiciones de contorno u(x, y) = g(x, y) en el límite
    ∂Ω. Otro aspecto, importante es que las iteraciones del modelo son de 2000, su tolerancia es
    1x10^-6.

    Entonces, al resolver el sistema de ecuaciones lineales, utilizaremos el método de Gauss-Seidel


    con control de convergencia. Este método se distingue por su habilidad para ajustar de manera
    iterativa los valores de las incógnitas hasta lograr una solución convergente.

    3
    En este sentido el algoritmo para resolver numéricante las 5 funciones del primer apartado del
    ejercicio:

    Por esto utilizaré tres niveles de simulación: uno con N=2, otro con N=3, y finalmente N=100,
    porque con N=1000 mi computadora ha colapsado, y porque consideré que es suficiente por las
    razones que defenderé a continuación.
    En una primera simulación, mediante N=2, aplicando este método se obtiene una malla muy
    liviana que solo sirve para evaluar si el código que me hemos implementado es capaz de evaluar
    o no el método en las 5 funciones de partida. Es decir si el método ha detectado problemas de
    frontera o no. El resultado del script que pondré después, son las 5 figuras de colores arcoiris
    mostradas de a continuación (figura 1).

    4
    Resolución numérica del problema con N=2.

    Figura 1: notese lo angulosa y sesgada representación, tal como se puede apreciar que el
    tener una malla poco densa no favorece para nada la interpretación de lo que realmente
    se esta intentando observar. No obstante al igual que con el caso de N=3, el resultado
    nos permitío verificar la funcionalidad del método.

    Dicho esto, presento el script de la parte correspondiente para N=2.

    5
    Resolución numérica del problema con N=3.

    La simulación con N=3, nos arroja el siguiente resultado reflejado en las 5 gráficas de la
    figura 2.

    Figura 2: Se nota una mejora en la representación de las funciones, de una malla con
    N=3, con esto podemos tener un primer visto de los ajustes para aumentar precison y
    estabilidad. En la etapa siguiente, nos acercaremos en teoría mucho más a la solución
    exacta de cada una de las funciones, por lo tanto al ser menores las diferencias entre el
    modelo numérico y el exacta la precisión en la solución del problema de Diritlech será
    mayor.

    Dicho, esto le proporciono la copia de script correspondiente con N=3.

    6
    Resolución numérica del problema con malla de N=100.

    Tal como se puede apreciar en la figura 3, las funciones contorno se han suavizado
    como era lo esperado, no obstante, para resolver el primer paso real de esta etapa hay
    que solucionar el cálculo para las 3 funciones frontera definidas en el enunciado.

    Figura 3: en este caso la solución del problema de Diritlech es mejor, ya que con las
    malla de 100 la representación es más parecida a la función exacta correspondiente a
    cada figura.

    El script de ejecución de este ejercicio es:

    7
    Segunda parte correspondiente con la segunda sugerencia descrita en el
    enunciado:

    Conociendo que el método funciona para mallas pequeñas de N=2, y N=3, y N=100,
    entonces procedemos a trabajar para resolver de forma numérica y exacta las tres
    funciones de contorno del enunciado (ver figura 4).

    Figura 4: resultados de las tres funciones 𝑥, 𝑥𝑦, 4𝑥𝑦(𝑥 − 𝑦) ∗ (𝑥 + 𝑦). Se puede


    distinguir que el mallado tiene una relación directa con la resolución espacial del
    problema, esto se debatirá mas adelante.

    El script de ejecución de esta parte es:

    8
    Realizado el ploteo de las solución aproximada, pasamos a probar si las soluciones
    numéricas cumplen o no con la formulación de Laplace, entonces luego de programar el
    siguiente script, se obtiene debajo de este una salida que muestra algo que es interesante
    saber:

    Lo interesante de ver el pequeño código siguiente, es que la salida del anterior script, en
    donde se nos dice que las tres soluciones numéricas del contorno, cumple la función de
    Laplace, porque las sumas de sus derivadas parciales segundas respecto de sus
    respectivas variables x e y dan 0.

    Realizado esto podemos calcular de forma exacta las soluciones de los tres contornos.

    La ejecución del código corrobora que únicamente las funciones g1(x, y) y g3(x, y)
    cumplen con la ecuación de Laplace. Por consiguiente, la solución u(x, y) en todo el
    dominio también cumple con dicha ecuación. En este escenario, el resultado derivado
    mediante métodos numéricos, como el de Gauss-Seidel, converge hacia la solución
    exacta del problema de Dirichlet, equiparable a la función de contorno en el límite.

    9
    3. Resolución del problema de Dirichlet para encontrar la convergencia de una
    solución numérica con respecto de una solución exacta.

    Para empezar con la resolución primero hay que plotear las 3 funciones exactas para lo
    cual tenemos el siguiente script , y probamos gráficamente si el resultado es una
    solución del problema.

    Cuyo resultado es el ploteo de la figur a 5.

    10
    Figura 5: podemos observar un aumento de la resolución espacial con respecto a la
    resolución de la figura 4, no obstante, muestran cierta relación espacial. Pero, para
    evitar susceptibilidades tenemos que determinar si aplicando la ecuación de Laplace a
    las soluciones exactas de las 5 funciones del problema.

    Entonces luego de haber programado un función de Laplace aplicada a las funciones del
    enunciada por medio del siguiente script:

    Se tiene que:

    Lo que implica que las 3 funciones exactas 𝑔(𝑥, 𝑦) = 𝑥; 𝑔(𝑥, 𝑦) = 4𝑥𝑦(𝑥 − 𝑦)(𝑥 +
    𝑦); y 𝑔(𝑥, 𝑦) = 𝑥 3 − 3𝑥𝑦 2 , cumplen la ecuación de Laplace, por lo tanto la solución
    numérica 𝑢(𝑥, 𝑦) de estas tres funciones dentro del dominio Ω, satisface la ecuación de
    Laplace. Por lo que debido Gauss Seidel convergerá en sus soluciones numéricas a una
    solución exacta del problema de Dirichlet.

    11

    También podría gustarte