Elliptic equations background

Elliptic equations can be thought of as being the steady state limit of the diffusion equation,

u_t = (\kappa u_x)_x + (\kappa u_y)_y + \psi

when both the 1) boundary conditions and 2) the source (or forcing term) \psi is time-independent. Under these conditions, the time-dependent terms vanish, and we are led to the elliptic equation (here in 2D),

(\kappa u_x)_x + (\kappa u_y)_y = f

where u is the dependent variable we are solving for, and here f is the forcing term.

This equation needs to 1) be satisfied by all points within a bounding region and 2) satisfy the boundary conditions on that region. This then can be interpreted as an instantaneous constraint on the system, much different than the wave-like solutions of hyperbolic equations which travel with finite speed.

Equation Discretization

The Poisson equation can be thought of as being the steady state limit of the diffusion equation. Can be thought of the steady state, static form of the diffusion equaiton - begin with equaiotn, expect relaxation to steady state (i.e. as t→ inf, d/dt → 0), get the following Poisson equation:

Matrix form, relaxation form

• Discretizatiion leads to a syustenm of eqns. There are 2 ways to solve the equations. Matrix direct methods and relaxtiuon iteration mrethods.

Matrix form (not followed)

I did not write a code that uses a direct matrix method for solving the system of equations, but include this section here for completeness.

Relaxation form (followed)

Error, accuracy, convergence

Laplace Equation (source term = 0) Physical Meaning

Note that the discretization equation with f=0 describes each point as being the average of its neighbors. This is in fact one of the key properties of the harmonic functions, which are the class of solution to the Laplace Equation. Since this is such an ubiquitious equation in physics, and a subset of the Poisson equation, I felt it was important to review some of the properties behind it.

• Boundary conditions
• Minimizing distance b/w boundaries (in 1d this is a line, in 2d this is soap bubble)
• Physical interprestaion - no charge within the domain, but elsewhere.

Poisson Equation Physical Meaning

While the interpreation of the solution propertoies are not as clear cut and intuitive as Laplace, the solution can be thoguht of as composing green's functions (cite references).

The code

Explain what my thing does, post it for download

Tests and Results

Test 1 - Laplace's Equation

Test 2 - Poisson's Equation with simple forcing

Test 3 - Poisson's Equation with complicated forcing

Properties - aka only an approximate solution now . Verify 2nd order (how?)

Cite References

References