**Fundamental Solutions of Co-dimension 1, 2, and 3**

**Method of Images**

**Finite Domains: Separation of Variables**

__Periodic or Finite Domains: Separation of Variables__

Consider the diffusion equation on the finite one-dimensional domain, , with Dirichlet BCs:

(11-1) |

and ICs

(11-2) |

Assuming that and inserting into the diffusion equation resulted in a pair of separated ordinary differential equations each of which was equal to the same constant, .

The possible solutions are:

(11-3) |

for and for ,

(11-4) |

The BCs cannot be solved simultaneously with a non-trivial solution unless .

Applying the BCs to Eq. , , but there are an infinite number of other values that will satisfy the boundary conditions:

(11-5) |

where is any integer.

The general solutions (for the BCs and ICs given above) are:

(11-6) |

However, to satisfy the initial conditions the general solution must be taken as a linear superposition of each possible solution--and with the weighting factor at chosen to satisfy the initial conditions. Therefore, seek such that:

(11-7) |

which is a Fourier sine series representation of .

__Remarks on the Fourier Method of Solution__

**(A)**

From the solution above, a periodic solution in can be constructed for and for all integer (use the in , etc).

**(B)**

Note how each term in the Fourier series decays much faster than the previous term--the
rate at which amplitudes get smaller is exponential with smaller wavelengths.

This dependence on wavelengths is a useful intuitive device for understanding the diffusion equation.

For example, consider two different initial conditions on the infinite domain :

The solutions and are obviously,

(11-9) |

Therefore shorter wavelengths, , decay times faster than longer wavelengths .

__Time Dependent Boundary Conditions, Semi-Infinite Domains__

The final analytic method for obtaining an analytic solution to to the diffusion equation is the method of Laplace transforms. This method is applicable to semi-infinite domains and is especially effective for problems with time-dependent boundary conditions.

The method of Laplace transforms is an example of operator calculus. In this case, the operator is a transformation of time derivatives to algebraic expressions in a transformation variable. The transformed diffusion equation becomes an inhomogeneous ordinary differential equation in the spatial variable. The ordinary differential equation is solved for the transformed boundary conditions and then the transformation is reversed--usually through a table of Laplace transform pairs.

The Laplace Transform is defined as the linear operator:

(11-10) |

that maps to another function .

The essential property of the Laplace transform is its operation on time derivatives:

(11-11) |

which can be integrated by parts,

(11-12) |

so the initial condition gets integrated into the the transformed derivative.

The transform has no effect on spatial derivatives, therefore the transformed diffusion equation becomes:

which is an ordinary differential equations since there are only derivatives with respect to .

__Two Examples of the Laplace Transform Method__

For this initial condition, the transformed diffusion equation, Eq. 11-13, becomes

The general solution, , to this differential equation is the sum of two different kinds of solutions. The first kind of solution is the one that satisfies Eq. 11-14 for a finite value of the right-hand-side ( , called the

For this problem, the homogeneous solution is:

(11-15) |

where

(11-16) |

The particular solution is:

(11-17) |

so the general solution is given by:

(11-18) |

The transformed boundary conditions become

(11-19) |

and

(11-20) |

This equation immediately above implies that and the other determines ,

(11-21) |

This equation requires transformation back to the variable through the inverse of a Laplace transform. The inverses are usually much more difficult to find than the Laplace transforms. Fortunately, Laplace transforms and their inverses are usually tabulated in math handbooks. For example,

Selected Laplace Transform Pairs | |

erfc | |

erfc |

So, the solution can be obtained through the use of the above table:

(11-22) |

which could have been guessed from the symmetry of the infinite domain solution.

For a second example with time-dependent solutions,

For this case, the transformed diffusion equation is the same as that calculated in the previous example:

(11-23) |

and the zero-flux boundary condition at again implies that . The transformed constant flux boundary condition is:

(11-24) |

Solving for the constant ,

(11-25) |

Using the table of Laplace transform pairs,

(11-26) |

Therefore, the surface concentration changes as :

(11-27) |

where is the constant surface flux.

__Anisotropic Diffusion Coefficients__

The methods of solution that have been treated do not account for the general case of an anisotropic diffusion coefficient. However, a simple strategem can be used to reconstruct the diffusion equation from its tensor form into the scalar form that is treated above.

Recall form of the the flux relationship in a particular frame of reference:

(11-28) |

The diffusion equation becomes:

(11-29) |

This would be difficult to solve in the general case. However, there is a trick for the case where the anisotropic diffusivity does not depend on .

First, find the rotation which diagonalizes :

(11-30) |

The are the eigenvalues of and the coordinate frame of reference of have axes parallel to the eigenvalues of .

The diffusion equation is much simpler in the eigenframe:

(11-31) |

Define an effective diffusivity :

det det | (11-32) |

and rescale the lengths of the coordinate axes with:

(11-33) |

Using the and the rescaled lengths,

(11-34) |

which is the simple scalar isotropic diffusion equation to which the solution methods (scaling, method of images, Fourier sum, and Laplace transform) were applied.