Abstract
When we use penalty method to solve the pure Neumann problem, the associated stiffness matrix we get is ill-conditioned, which leads to the unstable computation. In this paper, we design an iterative penalty method for the problem and prove some error estimates. In our algorithm, we can use a not very small penalty parameter to avoid the unstable computation. Numerical examples are given to show the algorithm is very effective and powerful.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: 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.
