Resumen
Resumo
Apresentamos neste trabalho a solução numérica de algumas equações de Poisson, uma equação diferencial parcial elíptica de segunda ordem, em malhas estruturadas bidimensionais e tridimensionais. Na determinação da solução numérica, empregamos o método iterativo SOR para solucionar o sistema de equações lineares proveniente da discretização da equação de Poisson por intermédio do método de diferenças finitas. Além disso, construímos algumas soluções manufaturadas 2D e 3D para a equação de Poisson, testamos valores ótimos para o parâmetro de sobrerrelaxação no método SOR e analisamos o comportamento dos métodos empregados na solução numérica de problemas 2D com singularidades. Na visualização das soluções manufaturadas e numéricas 2D e 3D, utilizamos, respectivamente, o Matlab e o Tecplot 360. Concluímos que a convergência do método SOR é lenta em problemas com condições de contorno de Neumann e em problemas com singularidades fortes.
Palavras-chave: Método de diferenças finitas. Método SOR. Soluções manufaturadas.
Abstract
Numerical solution of 2d and 3d poisson equation in structured meshes
We present in this work the numerical solution of some Poisson equations, an elliptic partial differential equation of second order, in two-dimensional and three-dimensional structured meshes. In determining the numerical solution, we used the iterative SOR method to solve the system of linear equations arising from the discretization of the Poisson equation using the finite difference method. Furthermore, we build some 2D and 3D manufactured solutions for the Poisson equation, and test optimal values for the over-relaxation parameter in the SOR method and analyze the behavior of the methods used in the numerical solution of 2D problems with singularities. In the visualization of the 2D and 3D manufactured and numerical solutions, we used, respectively, Matlab and Tecplot 360. We concluded that the convergence of the SOR method is slow in problems with Neumann boundary conditions and in problems with strong singularities.
Keywords: Finite difference method. SOR method. Manufactured solutions.
Citas
BURDEN, R. L.; FAIRES, D. J.; BURDEN, A. M. Análise numérica. 3. ed. São Paulo: Cengage Learning, 2016.
CENICEROS, H. D.; NÓS, R. L.; ROMA, A. M. Three-dimensional fully adaptive simulations of phase-field fluid models. Journal of Computational Physics, Amsterdam, v. 229, p. 6135-6155, 2010.
CHENG, L. Finite difference methods for Poisson equation. 2020. Disponível em: https://www.math.uci.edu/~chenlong/ 226/FDM.pdf. Acesso em: 16 out. 2021.
FERZIGER, J. H. Numerical methods for engineering application. New York: John Wiley & Sons, 1981.
FORTUNA, A. de O. Técnicas computacionais para dinâmica dos fluidos: conceitos básicos e aplicações. São Paulo: Edusp, 2000.
GARCIA, M. V. P.; HUMES, C.; STERN, J. M. Generalized line criterion for Gauss-Seidel method. Computational and Applied Mathematics, New York, v. 22, n. 1, p. 91-97, 2003.
IÓRIO, V. EDP: um curso de graduação. Rio de Janeiro: IMPA, 1989.
KALABA, R. E.; SPINGARN, K. A criterion for the convergence of the Gauss-Seidel method. Applied Mathematics and Computation, Amsterdam, v. 4, n. 4, p. 359-367, 1978.
LAPLACE, C. Dev-C++ official website. 2020. Disponível em: https://www.bloodshed.net/. Acesso em: 16 out. 2021.
LOGAN, J. D. Applied partial differential equations. 3. ed. New York: Springer, 2015.
MATLAB. Math Works. 2021. Disponível em: https://www.mathworks.com/products/matlab.html. Acesso em: 16 out. 2021.
MALISKA, C. R. Transferência de calor e mecânica dos fluidos computacional. Rio de Janeiro: LTC, 2004.
MICHELETTI, J. P. S. B. Solução numérica da equação de Poisson em malhas estruturadas bidimensionais e tridimensionais. 2021. 111 f. Trabalho de Conclusão de Curso (Licenciatura em Matemática) – Universidade Tecnológica Federal do Paraná, Curitiba, 2021.
MITCHELL, W. F. A collection of 2D elliptic problems for testing adaptive grid refinement algorithms. Preprint submitted to Elsevier, Amsterdam, p. 1-19, 2017.
NÓS, R. L. et al. Three-dimensional coarsening dynamics of a conserved, nematic liquid crystal-isotropic fluid mixture. Journal of Non-Newtonian Fluid Mechanics, Amsterdam, v. 248, p. 62-73, 2017.
NÓS, R. L.; CENICEROS, H. D.; ROMA, A. M. Simulação tridimensional adaptativa da separação das fases de uma mistura bifásica usando a equação de Cahn-Hilliard. TEMA Tendências em Matemática Aplicada e Computacional, São Carlos, v. 13, p. 37-50, 2012.
NÓS, R. L.; ROMA, A. M.; CENICEROS, H. D. Solução de equações diferenciais parciais elípticas por técnicas multinível-multigrid em malhas tridimensionais bloco-estruturadas com refinamento localizado. In: CONGRESSO NACIONAL DE MATEMÁTICA APLICADA E COMPUTACIONAL, 28, 2005, Santo Amaro. Anais[...]. Santo Amaro: SBMAC, 2005.
ONABID, M. A. Solving three-dimensional (3D) Laplace equations by successive over-relaxation method. African Journal of Mathematics and Computer Science Research, Lagos, v. 5, n. 13, p. 204-208, 2012. Disponível em: https://academicjournals.org/article/article1379673859_Onabid.pdf. Acesso em: 22 out. 2021.
SCHILDT, H. C. Completo e total. 3. ed. São Paulo: Makron Books, 1997.
STRIKWERDA, J. C. Finite difference schemes and partial differential equations. New York: Chapman & Hall, 1989.
TECPLOT. Tecplot 360. 2021. Disponível em: https://www.tecplot.com/products/tecplot-360/. Acesso em: 16 out. 2021.
UM, S. Y. Dispersion and local-error analysis of compact LFE-27 formula for obtaining sixth-order accurate numerical solutions of 3D Helmholtz equation. 2013. Disponível em: https://www.researchgate.net/publication/275876371_Dispersion_and_local-error_analysis_of_compact_LFE-27_formula_for_obtaining_sixth-order_accurate_numerical_solutions_of_3D_Helmholtz_equation/figures?lo=1. Acesso em: 16 out. 2021.
WILLIAMS, G. Five-point stencil. 2011. Disponível em: https://source.ggy.bris.ac.uk/wiki/File:Five-point-stencil.jpg. Acesso em: 16 out. 2021.
KALABA, R. E.; SPINGARN, K. A criterion for the convergence of the Gauss-Seidel method. Applied Mathematics and Computation, v. 4, n. 4, p. 359-367, 1978.
LAPLACE, C. Dev-C++ official website. 2020. Disponível em: https://www.bloodshed.net/. Acesso em: 16 out. 2021.
LOGAN, J. D. Applied partial differential equations. 3. ed. New York: Springer, 2015.
MATLAB. MathWorks. 2021. Disponível em: https://www.mathworks.com/products/matlab.html. Acesso em: 16 out. 2021.
MALISKA, C. R. Transferência de calor e mecânica dos fluidos computacional. Rio de Janeiro: LTC, 2004.
MICHELETTI, J. P. S. B.. Solução numérica da equação de Poisson em malhas estruturadas bidimensionais e tridimensionais. Orientador: Rudimar Luiz Nós. 2021. 111 f. Trabalho de Conclusão de Curso (Licenciatura em Matemática) – Departamento Acadêmico de Matemática, UTFPR, Curitiba, 2021.
MITCHELL, W. F. A collection of 2D elliptic problems for testing adaptive grid refinement algorithms. Preprint submitted to Elsevier, p. 1-19, 2017.
NÓS, R. L. et al. Three-dimensional coarsening dynamics of a conserved, nematic liquid crystal-isotropic fluid mixture. Journal of Non-Newtonian Fluid Mechanics, v. 248, p. 62-73, 2017.
NÓS, R. L.; CENICEROS, H. D.; ROMA, A. M. Simulação tridimensional adaptativa da separação das fases de uma mistura bifásica usando a equação de Cahn-Hilliard. TEMA Tendências em Matemática Aplicada e Computacional, v. 13, p. 37-50, 2012.
NÓS, R. L.; ROMA, A. M.; CENICEROS, H. D. Solução de equações diferenciais parciais elípticas por técnicas multinível-multigrid em malhas tridimensionais bloco-estruturadas com refinamento localizado. In: XXVIII Congresso Nacional de Matemática Aplicada e Computacional. Santo Amaro: SBMAC, 2005.
ONABID, M. A. Solving three-dimensional (3D) Laplace equations by successive over-relaxation method. African Journal of Mathematics and Computer Science Research, v. 5, n. 13, p. 204-208, 2012. Disponível em: https://academicjournals.org/article/article1379673859_Onabid.pdf. Acesso em: 22 out. 2021.
SCHILDT, H. C – Completo e total. 3. ed. São Paulo: Makron Books, 1997.
STRIKWERDA, J. C. Finite difference schemes and partial differential equations. New York: Chapman & Hall, 1989.
TECPLOT. Tecplot 360. 2021. Disponível em: https://www.tecplot.com/products/tecplot-360/. Acesso em: 16 out. 2021.
UM, S.-Y. Dispersion and local-error analysis of compact LFE-27 formula for obtaining sixth-order accurate numerical solutions of 3D Helmholtz equation. 2013. Disponível em: https://www.researchgate.net/publication/275876371_Dispersion_and_local-error_analysis_of_compact_LFE-27_formula_for_obtaining_sixth-order_accurate_numerical_solutions_of_3D_Helmholtz_equation/figures?lo=1. Acesso em: 16 out. 2021.
WILLIAMS, G. Five-point stencil. 2011. Disponível em: https://source.ggy.bris.ac.uk/wiki/File:Five-point-stencil.jpg. Acesso em: 16 out. 2021.
Esta obra está bajo licencia Creative Commons Attribution-NonCommercial 4.0 International License.
Derechos de autor 2022 Array