next up previous
Next: Examples of Gradient Flows Up: Variational Methods for Microstructural Previous: Inner Products and Norms

Inner Products on Fields

The same properties for norms and inner products on vectors apply to norms on fields. What does the magnitude of a function mean? What does it mean for a function to be near zero?

There are an infinite number of inner products and norms on fields. A common choice is to sum the point-wise multiplication of two fields over a domain tex2html_wrap_inline1044 :

  equation687

this is the definition of the tex2html_wrap_inline1016 inner product and norm (functions defined on a domain with an inner product make up a Hilbert space). The tex2html_wrap_inline1016 norm is what is minimized in a parametric least squares fit of a function to data. If f and g have integral zero (which can be produced by subtracting of their averages), then the norm is a measure of `how far the function is from zero.' For the rest of this section, all functions assumed to have integral zero.

However, our notion of the distance of a function from zero may depend on the physical context. For example, consider the function defined on the periodic domain tex2html_wrap_inline1054 :

  equation689

which is plotted in figure 2.1

tex2html_wrap1176
* Figure tex2html_wrap1428 . tex2html_wrap1425

With the tex2html_wrap_inline1016 norm, tex2html_wrap_inline1058 is independent of tex2html_wrap_inline1060 in Fig. 2.1; this is reasonable if ``stuff'' need not flow from one bump to another. However, if flow is required then the tex2html_wrap_inline1016 norm is not a useful measure of the distance of the function from zero.

The tex2html_wrap_inline1064 (called `H minus one') norm is particularly useful for diffusive type processes in which there is a conservation law.

equation692

where tex2html_wrap_inline1066 is the solution to the Poisson PDE: tex2html_wrap_inline1068 . As a result the set of functions in tex2html_wrap_inline1064 inner products is limited to those with integral 0, and thus tex2html_wrap_inline1064 is a prime candidate for variables that are conserved. It may be helpful to visualize tex2html_wrap_inline1066 as the steady state concentration profile for a mass source density given by tex2html_wrap_inline960 ; then the tex2html_wrap_inline1064 norm is the tex2html_wrap_inline1016 norm of the flux densitygif.

The natural inner product for tex2html_wrap_inline1064 is:

  equation700

We examine how the inner product behaves on the functions defined in Figure 2.1, by integrating the function twice to get tex2html_wrap_inline1066 :

tex2html_wrap1178
* Figure tex2html_wrap1429 . tex2html_wrap1426 tex2html_wrap_inline1066 and tex2html_wrap_inline1098 for the two functions in Figure 2.1 in the limit of small tex2html_wrap_inline1100 .

If the tex2html_wrap_inline1064 norm of f is now plotted as a function of tex2html_wrap_inline1060 in the limit of small tex2html_wrap_inline1100 , we get Figure 2.3.

tex2html_wrap1180
* Figure tex2html_wrap1430 . tex2html_wrap1427 tex2html_wrap_inline1064 norm of the function in Figure 2.1 as a function of the bump separation in the limit of small tex2html_wrap_inline1100 : tex2html_wrap_inline1114 . This illustrates how the magnitude of the function increases with how much the system must diffuse to reach a uniform value.

The tex2html_wrap_inline1064 inner product can be written in a more convenient form, by noting that tex2html_wrap_inline1118 and using the divergence theorem:

  equation708

where the integral at the boundary is assumed to be zero and the operator tex2html_wrap_inline1120 is the `inverse Laplacian.' The last form is directly usable for finding the gradient flow for f.

To see how the inner product is related to the functional norm, let F be the total free energy of a system with a free energy density which can be written a function of position, concentration and concentration gradients, tex2html_wrap_inline1126 at time t=0: tex2html_wrap_inline1130 . With a flow field, tex2html_wrap_inline1132 , the total free energy as a function of time is approximately F(f + vt) for small t>0. The initial rate of change of F is:

  equation711

where tex2html_wrap_inline1142 is the variational derivative of F-the left-hand-side of Euler's Equation from variational calculus[5] from using the divergence theorem to remove derivatives of v in the volume integral, assuming the boundary integral is zerogif. Again, tex2html_wrap_inline1156 will depend on the choice of the inner product. The integral over the boundary tex2html_wrap_inline1148 usually vanishes identically, either because admissible flows v vanish on tex2html_wrap_inline1148 (Dirichlet conditions), or because the function g, which depends on derivatives of f with respect to its gradient and higher order-derivatives, vanishes (Neumann conditions).

Since the inner product depends on the rate of change, it is natural to locally `weight' these inner products with some constant or function tex2html_wrap_inline1168 representing the mobility. For example in tex2html_wrap_inline1016 mobility weighting, a function with a `bump' in a region of high mobility is closer to zero than it would be if the bump is a region of low mobility. Thus the tex2html_wrap_inline1016 inner product is usually modified to

equation715

and the tex2html_wrap_inline1064 inner product modified to

equation717

We now show how standard evolutions result from doing gradient flows with these inner products.


next up previous
Next: Examples of Gradient Flows Up: Variational Methods for Microstructural Previous: Inner Products and Norms

W. Craig Carter
Tue Sep 30 16:07:27 EDT 1997