¿Cuál es la tasa de convergencia teórica para un solucionador de veneno FFT? Estoy resolviendo una ecuación de Poisson: con en el dominio con periódico condición de frontera. Esta densidad de carga es neta neutra. La solución viene dada por: donde . En el espacio recíproco donde son los...
Recientemente he estado comparando diferentes solucionadores no lineales de scipy y quedé particularmente impresionado con el ejemplo de Newton-Krylov en el Scipy Cookbook en el que resuelven una ecuación de ecuación diferencial de segundo orden con un término de reacción no lineal 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...
Considere una simétrica definida positiva tridiagonal sistema lineal , donde A ∈ R n × n y b ∈ R n . Dados tres índices 0 ≤ i < j < k < n , si asumimos que solo las filas de ecuaciones estrictamente entre i y k se mantienen, podemos eliminar variables intermedias para obtener una...
Me gustaría saber cómo se aplican normalmente las condiciones de Dirichlet cuando se utiliza el método de volumen finito en una cuadrícula no uniforme centrada en celdas, Mi implementación actual simplemente impone la condición de límite al fijar el valor de la primera celda, ϕ1= gre(...
He oído que se puede usar una transformada rápida de Fourier para resolver el problema de Poisson cuando las condiciones de contorno son todas de un solo tipo ... Serie seno para dirichlet, coseno para neumann y ambas para periódico. Considerando un dominio rectangular 2D, suponga que dos lados...
He estado tratando de depurar este error los últimos días, me preguntaba si alguien tiene consejos sobre cómo proceder. Estoy resolviendo la ecuación de Poisson para una distribución de carga escalonada (un problema común en la electrostática / física de semiconductores) en una malla de volumen...
Planteamiento del problema Implementé multirredes geométricas para donde f = 3 π 2−∇2=f−∇2=f-\nabla^{2}=f enΩ∈[0,1]en uncubounitario. Los límites de Dirichlet en la cara izquierda, la cara inferior y la cara frontal son0. Los límites de Neumann en la cara superior, derecha y posterior...