A spectral method for elliptic equations: The Neumann problem

Boundary value problems for second-order elliptic equations (35J25) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Spectral, collocation and related methods for boundary value problems involving PDEs (65N35)

Abstract: Let

O m e g a

be an open, simply connected, and bounded region in

m a t h b b R^{d}

,

d g e q 2

, and assume its boundary

p a r t i a l O m e g a

is smooth. Consider solving an elliptic partial differential equation

- D e l t a u + g a m m a u = f

over

O m e g a

with a Neumann boundary condition. The problem is converted to an equivalent elliptic problem over the unit ball

B

, and then a spectral Galerkin method is used to create a convergent sequence of multivariate polynomials

u_{n}

of degree

l e q n

that is convergent to

u

. The transformation from

O m e g a

to

B

requires a special analytical calculation for its implementation. With sufficiently smooth problem parameters, the method is shown to be rapidly convergent. For

u i n C^{i n f t y} (o v e r l i n e O m e g a)

and assuming

p a r t i a l O m e g a

is a

C^{i n f t y}

boundary, the convergence of

V e r t u - u_{n} V e r t_{H^{1}}

to zero is faster than any power of

1 / n

. Numerical examples in

m a t h b b R^{2}

and

m a t h b b R^{3}

show experimentally an exponential rate of convergence.

Recommendations

Cites work

Cited in

(21)

This page was built for publication: A spectral method for elliptic equations: The Neumann problem

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2431114)