2002
DOI: 10.1002/num.10006
|View full text |Cite
|
Sign up to set email alerts
|

Analyzing and visualizing a discretized semilinear elliptic problem with Neumann boundary conditions

Abstract: A semilinear elliptic equation d∆u − u + u p = 0 over the unit ball in R 2 with positive solution and the homogeneous Neumann boundary condition is considered. This equation models applications like chemotactic aggregation and biological pattern formation. Focusing on solving the discretized version of the equation, this work proposes an efficient algorithm that combines a newly developed discretization scheme on polar coordinates with a fast Fourier solver. An analysis of the induced matrix structures proves … Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1

Citation Types

0
2
0

Year Published

2003
2003
2008
2008

Publication Types

Select...
4

Relationship

1
3

Authors

Journals

citations
Cited by 4 publications
(2 citation statements)
references
References 16 publications
0
2
0
Order By: Relevance
“…If the components of q 0 are nonnegative, this property is preserved by each iteration q j , and hence also by the limit vector if it exists (see [15,Theorem 3.1]). The convergence of a subsequence of this iteration method to a nonzero vector is proved in [15,Theorem 2.1]. Although the convergence of the entire sequence is not proved, it is observed numerically to be very robust.…”
Section: Goto Step 1 Endmentioning
confidence: 99%
“…If the components of q 0 are nonnegative, this property is preserved by each iteration q j , and hence also by the limit vector if it exists (see [15,Theorem 3.1]). The convergence of a subsequence of this iteration method to a nonzero vector is proved in [15,Theorem 2.1]. Although the convergence of the entire sequence is not proved, it is observed numerically to be very robust.…”
Section: Goto Step 1 Endmentioning
confidence: 99%
“…Note that fine meshes are used in the heterojunctions and a half of the mesh length is shined in the radial coordinate• coordinate r and the natural axial coordinate z with the following two special treatments• First, in the heterojunction area, fine meshes are used to capture the rapid change in the wave functions. Secondly, half of the mesh length is shifted in the radial coordinate to avoid incorporating the pole condition [38,39]. See Figure 2 for details• Based on the grid points, equation (1) is discretized by the 3D centered seven-point finite difference method …”
Section: Numerical Schemesmentioning
confidence: 99%