Preguntas etiquetadas con finite-difference

Estoy tratando de resolver la ecuación de advección, pero aparece una extraña oscilación en la solución cuando la onda se refleja desde los límites. ¡Si alguien ha visto este artefacto antes, me interesaría saber la causa y cómo evitarlo! Este es un gif animado, se abre en una ventana separada...

Descripción del experimento: En la interpolación de Lagrange, la ecuación exacta se muestrea en puntos (orden polinomial ) y se interpola en 101 puntos. Aquí es variado de 2 a 64. Cada vez , y parcelas de error se preparan. Se ve que, cuando la función se muestrea en puntos equidistantes, gotas...

Cuál sería una buena discretización de diferencia finita para la siguiente ecuación: ∂ρ∂t+∇⋅(ρu)=0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? Podemos tomar el caso 1D: ∂ρ∂t+ddx(ρu)=0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho...

Las fórmulas de diferencia finita estándar se pueden usar para calcular numéricamente una derivada con la expectativa de que tenga valores de función en puntos espaciados uniformemente, de modo que es una constante. ¿Qué sucede si tengo puntos desigualmente espaciados, de modo que ahora varía de un...

Para un proyecto en el que estoy trabajando (en PDE hiperbólicas), me gustaría obtener una idea aproximada del comportamiento observando algunos números. Sin embargo, no soy un muy buen programador. ¿Puede recomendarme algunos recursos para aprender a codificar efectivamente esquemas de...

Probablemente sea una pregunta a nivel de estudiante, pero no puedo dejarla clara para mí. ¿Por qué es más preciso usar cuadrículas no uniformes en los métodos numéricos? Estoy pensando en el contexto de algún método de diferencia finita para el PDE de la forma . Y suponga que estoy interesado en...

Estoy interesado en resolver la ecuación de Poisson utilizando el enfoque de diferencias finitas. Me gustaría entender mejor cómo escribir la ecuación matricial con las condiciones de contorno de Neumann. ¿Alguien revisaría lo siguiente, es correcto? La matriz de diferencia finita La ecuación de...

Estoy tratando de resolver una ecuación de Poisson en 2D por diferencias finitas. En el proceso, obtengo una matriz dispersa con solo variables en cada ecuación. Por ejemplo, si las variables fueran , entonces la discretización produciría:5 55 55UUU Ui - 1 , j+ Ui + 1 , j- 4 Ui , j+ Ui , j - 1+ Ui...

Existen numerosos esquemas FD para la ecuación de advección discutir en la web. Por ejemplo aquí: http://farside.ph.utexas.edu/teaching/329/lectures/node89.html∂T∂t+ u ∂T∂X= 0∂T∂t+tu∂T∂X=0 0\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=0 Pero no he visto a nadie proponer un esquema...

Tengo un problema cuando quiero usar la aproximación de diferencia de centro de orden superior: (−ui+2,j+16ui+1,j−30ui,j+16ui−1,j−ui−2,j12)(−ui+2,j+16ui+1,j−30ui,j+16ui−1,j−ui−2,j12)\left(\frac{-u_{i+2,j}+16u_{i+1,j}-30u_{i,j}+16u_{i-1,j}-u_{i-2,j}}{12}\right) para la ecuación de Poisson en un...

Por mucho que trato de encontrar una explicación concisa en Internet, parece que no puedo entender el concepto de una diferencia finita mimética, o cómo incluso se relaciona con las diferencias finitas estándar. Sería realmente útil ver algunos ejemplos simples de cómo se implementan para PDE...

Estoy tratando de encontrar algunos recursos para ayudar a explicar cómo elegir las condiciones de contorno cuando se utilizan métodos de diferencias finitas para resolver PDE. Todos los libros y notas a los que tengo acceso actualmente dicen cosas similares: Las reglas generales que gobiernan...

Estoy interesado en implementar una malla móvil para un problema de advección-difusión. Adaptive Moving Mesh Methods da un buen ejemplo de cómo hacer esto para la ecuación de Burger en 1D usando diferencia finita. ¿Alguien podría ofrecer un ejemplo práctico para resolver la ecuación de...

Estoy tratando de aprender sobre la resolución numérica de PDE por mí mismo. He estado comenzando con el método de diferencia finita (FDM) durante algún tiempo porque escuché que FDM es el fundamento de numerosos métodos numéricos para PDE. Hasta ahora, tengo algunos conocimientos básicos sobre...

Estoy trabajando en resolver las ecuaciones de poroelasticidad unidimensionales acopladas (modelo de biot), dadas como: ∂−(λ+2μ)∂2u∂x2+∂p∂x=0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} =

Cuando se quiere calcular derivadas numéricas, el método presentado por Bengt Fornberg aquí (y reportado aquí ) es muy conveniente (preciso y simple de implementar). Como el documento original data de 1988, me gustaría saber si hay una mejor alternativa hoy (como (o casi como) simple y más...

Estoy intentando resolver una ecuación del tipo: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Donde f(x)f(x)f(x) tiene un polo simple en 000 , para los valores propios y vectores propios más pequeños...

Estoy implementando el documento " Transporte masivo óptimo para el registro y la deformación ", mi objetivo es ponerlo en línea ya que simplemente no puedo encontrar ningún código de transporte masivo euleriano en línea y esto sería interesante al menos para la comunidad de investigadores en el...

Tengo una pregunta sobre la codificación de las condiciones de contorno para la mecánica de sólidos (elasticidad lineal). En el caso especial, tengo que usar diferencias finitas (3D). Soy muy nuevo en este tema, por lo que quizás algunas de las siguientes preguntas pueden ser muy básicas. Para...

¿Alguien puede ayudarme a encontrar los libros sobre soluciones numéricas (diferencia finita y métodos de Crank-Nicolson) de Poisson y ecuaciones de difusión que incluyen ejemplos sobre geometría irregular, como un dominio que consiste en el área entre un rectángulo y un círculo (especialmente...