Abstract
AbstractIn problems such as the computation of incompressible flows with moving boundaries, it may be necessary to solve Poisson's equation on a large sequence of related grids. In this paper the LU decomposition of the matrix A0 representing Poisson's equation discretized on one grid is used to efficiently obtain an approximate solution on a perturbation of that grid. Instead of doing an LU decomposition of the new matrix A, the RHS is perturbed by a Taylor expansion of A−1 about A0. Each term in the resulting series requires one ‘backsolve’ using the original LU.Tests using Laplace's equation on a square/rectangle deformation look promising; three and seven correction terms for deformations of 20% and 40% respectively yielded better than 1% accuracy.As another test, Poisson's equation was solved in an ellipse (fully developed flow in a duct) of aspect ratio 2/3 by perturbing about a circle; one correction term yielded better than 1% accuracy.Envisioned applications other than the computation of unsteady incompressible flow include: three‐dimensional parabolic problems in tubes of varying cross‐section, use of ‘elimination’ techniques other than LU decomposition, and the solution of PDEs other than Poisson's equation.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
More From: International Journal for Numerical Methods in Fluids
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.