nextuppreviouscontents
Next:4 Fourier Series for boundary Up:3 Construction of Green's function Previous:3.1 One Dimension

3.2 Two Dimensions

Our development of the Green's function method for the partial differential equation, such as problem (2), will follow essentially the same lines as for the case of ordinary differential equations. Our starting point is equation (1.12). We briefly choose "v" to be our Green's function G, and use x, y to denote the fixed point, with tex2html_wrap_inline6093tex2html_wrap_inline6113 as the corresponding dummy integration variables. Considering the problem (2) as an example, requiring that G satisfy

eqnarray802

over tex2html_wrap_inline6025 , the last term in (1.12) reduces to u(x,y). Replacing tex2html_wrap_inline6233 in the left-hand side by the known function tex2html_wrap_inline6235 , and subjecting G to boundary conditions on tex2html_wrap_inline6239 which result in the elimination of any normal derivative terms in the boundary integral. We can solve for the desired u(x,y) in terms of quantities which are known, and eventually, the Green's function can be determined. As we mentioned before, the domain tex2html_wrap_inline6025 that we considered in this report is a unit square, lets say tex2html_wrap_inline6245 , and according to equation (3.18) imposed the boundary condition on G, therefore, the corresponding boundary value problem on the Green's function will become

eqnarray804

and the eigenfunction method can be used to determine G.

Now, the eigenvalue problem corresponding to (3.19) is

eqnarray806

Solving this by separation of variables, we seek

eqnarray808

Inserting this into (3.20), and dividing through by XY we have

eqnarray810

Since the left-hand side is a function of tex2html_wrap_inline6093 alone, and the right-hand side is a function tex2html_wrap_inline6113 alone, we observe that (3.22) can hold for all tex2html_wrap_inline6093 's and tex2html_wrap_inline6113 's in tex2html_wrap_inline6025 only if both sides are constant, say tex2html_wrap_inline6271 . As for the boundary conditions on X and Y, the conditions

eqnarray812

imply that X(0)=X(1)=Y(0)=Y(1)=0. Thus

eqnarray814

For the first ordinary differential equation,

eqnarray816

Either A and/or tex2html_wrap_inline6281 must be zero; A=0 is unacceptable since it leads to the trivial solution tex2html_wrap_inline6285 and hence tex2html_wrap_inline6287 , whereas eigenfunctions are, by definition, nontrivial solutions of (3.20). To have a nontrivial solution, the constant k must be chosen such that k coincides with a zero of the sine function, that is, tex2html_wrap_inline6293 , where tex2html_wrap_inline6295 . The constant A remains arbitrary, and

eqnarray818

Inserting tex2html_wrap_inline6293 into the second equation of (3.23), and proceeding as above, we have

eqnarray820

so that, for nontrivial solutions, we need

displaymath6299

or, the eigenvalues of (3.23) are given by

eqnarray822

and the corresponding eigenfunctions are

eqnarray824

Next, we expand the quantities in (3.19) in terms of these eigenfunctions

eqnarray826

eqnarray828

Unknowns tex2html_wrap_inline6301 are the Fourier coefficient and tex2html_wrap_inline6303 are given by

eqnarray830

where we define the two-dimensional inner product by

eqnarray832

Accordingly,

eqnarray834

so that

eqnarray836

Now, inserting (3.29) and (3.30) into (3.19), and noting from (3.20) that tex2html_wrap_inline6305 , we have

eqnarray838

so that, equating coefficients of tex2html_wrap_inline6307 ,

eqnarray840

and

eqnarray842

With tex2html_wrap_inline6309 given by (3.28), tex2html_wrap_inline6311 by (3.27), and

eqnarray844

we have, finally,

eqnarray846


nextuppreviouscontents
Next:4 Fourier Series for boundary Up:3 Construction of Green's function Previous:3.1 One Dimension
Cheung Sau Hung

Wed Sep 15 10:03:39 HKT 1999