There are all sorts of resonances around us, in the world, in our culture, and in our technology. A tidal resonance causes the 55 foot tides in the Bay of Fundy. Mechanical and acoustical resonances and their control are at the center of practically every musical instrument that ever existed. Even our voices and speech are based on controlling the resonances in our throat and mouth. Technology is also a heavy user of resonance. All clocks, radios, televisions, and gps navigating systems use electronic resonators at their very core. Doctors use magnetic resonance imaging or MRI to sense the resonances in atomic nuclei to map the insides of their patients. In spite of the great diversity of resonators, they all share many common properties. In this blog, we will delve into their various aspects. It is hoped that this will serve both the students and professionals who would like to understand more about resonators. I hope all will enjoy the animations.

For a list of all topics discussed, scroll down to the very bottom of the blog, or click here.

Origins of Newton's laws of motion

Non-mathematical introduction to relativity

Three types of waves: traveling waves, standing waves and rotating waves new

History of mechanical clocks with animations
Understanding a mechanical clock with animations
includes pendulum, balance wheel, and quartz clocks

Water waves, Fourier analysis



Saturday, May 7, 2011

Conformal mapping I

Figure 1. Some common conformal mappings as applied to images, fields, and a grid. The four mapping functions are shown in the horizontal set of buttons (they require pushing "transform" to activate). The first function, zn has a user adjustable n value (type in a new value between 0 and 2 in the n = box). The various images are shown as the vertical set of buttons on the right (give each image a little time to load).

This flash animation is intended to illustrate that conformal mapping is very much like morphing an image in photo shop, but in this case we morph one valid fluid flow or electric field pattern into another. This allows us to start with a very simple flow or electric field and from it to "graphically" obtain all sorts of other flow and electric field patterns.

The initial coordinates, i.e. z-space or (x,y), are shown as a green grid and are warped by the transformations. The final coordinates, i.e. w-space or (u,v), are shown in light gray. Note that these are complex spaces, i.e.:
 z = x + iy  and  w = u + iv .

To aid understanding, we show the transformation process as a series of small steps. Normally, this is done a single step. The transformations also require scaling factors, e.g. w = zn is really more like w = a(bz)n

where a and b are real scaling factors.

Conformal mappings use complex functions that are "analytic". One of their properties is that they are orthonormal. This means that even though a mapped grid is warped, the originally perpendicular grid lines will remain perpendicular at the intersection points. Another property is that if we map a field that satisfies Laplace's equation, such as simple fluid flow or magnetic fields, the mapped fields will also satisfy Laplace's equations and are therefore valid fields.

We show the transformations w = zn, w = ez, w = sinz, w = z + 1/z. There are many more transformations, such as the rest of the complex trigonometric functions, complex hyperbolic functions, and polynomials in z.

Other postings by the author
up: Conformal mapping contents
next: Conformal mapping II
This posting includes a flash animation showing the physics discussed. Most computers have a flash player already installed, but if yours does not, download the free Adobe flash player here.


Conformal mapping to solve Laplace's Equation

Laplace's equation is a defining equation of electrostatics, low speed fluid flow, and gravitational fields. In two dimensions Laplace's Equation can be written as

two dimensional Laplace's equation       (1)

where Φ is the potential as appropriate for the problem at hand. In the usual problem, we know the values of Φ on the boundaries of a particular region of space and we want to solve for Φ in the interior of that region. Conformal mapping is one of the standard methods for producing such a solution. For another solution method, i.e. separation of variables, see my postings on water waves.

The mappings

Conformal mapping makes use of a particular type of function of complex numbers called analytic functions. These functions map the points in the complex plane into a distorted second complex space. Figure 1 at the right demonstrates a few such mappings, where the complex planes shown there also include images to aid in the understanding of this mapping process. This mapping process is similar to the projections of the Earth's spherical surface onto a flat sheet of paper to make the standard maps of the world. (Histories of map projections are here and here.) What makes conformal mappings so useful is that they map a solution of Laplace's equation into another solution of Laplace's equation.

There are many complex analytic functions and each one defines a particular mapping (or distortion) of the complex plane as shown in Figure 1. Additional mappings can be found at Wikipedia, Wolfram, Mathews-Howell/Fullerton or Rotorbrain/Foote-interactive mapping, to name a few places.

Our first example

Below, in Fig. 2 we start with a very simple solution to Laplace's equation, a solution representing constant horizontal flow of a fluid above a flat plate. We then use various conformal mappings to distort this solution into other solutions of Laplace's equation, representing flow around various distortions of that plate. We see that our one simple starting solution can give rise to a whole family of solutions.

We go on to use conformal mapping to get algebraic solutions to all the various distorted solutions.

Example 1 - Fluid flow transformed by zn
starting fluid flow pattern Fig. 2a. Very simple, constant horizontal flow with the associated velocity potential shown. The associated potential is given by the equation:
    Φ = x       (2.1)   
and results in a fluid velocity as given by:
    v = −∇Φ = −∂Φ/∂x ax − ∂Φ/∂y ay = −ax    ,    (2.2)   
where ax and ay are the unit vectors in the x and y directions. (For an explanation of velocity potentials see my postings on water waves.)

This flow pattern can be mapped into more complex geometries using conformal mapping as we show next.
flow around a corner Fig. 2b. This figure shows the result of applying the conformal transform,
    w = z½       (2.3)   
to the solution shown in Fig. 2a. This transform is also shown in Figure 1, as applied to various images.

The flow shown in Fig. 2b is that of fluid circulating around a corner made of rigid material. Because the original flow (that in Fig. 2a) satisfied Laplaces equation and because the transformation was conformal, we can be assured that the transformed flow will also satisfy Laplace's equation. It therefore will be the correct flow solution for irrotational, incompressible flow in a corner (see my earlier posting for a discussion of irrotational incompressible flow.)

Basically the square root transform used here bends the x axis at the origin and rotates it by 90° as shown by the red arrow. The y axis is rotated by 45°.

A detailed, more accurate mapping of this function is shown in a later posting.

Mathematics of z½

Sometimes we only need a picture of the flow in a given situation, in which case we simply need to do the image transformation as shown in Fig. 1 and we are done. Other times, we need an algebraic solution. We next use the transformation to derive an algebraic solution for this corner flow pattern.

We are given the initial potential in Equation (2.1) in the form of Φ as a function of (x,y) (actually, in this case it is a function of only x). We now need Φ as a function of (u,v) in Fig. 2b, where (u,v) are the rectilinear coordinates of the warped image. Thus we need x as a function of (u,v) so we can substitute this into Equation 2.1 to convert this equation into one that varies with u and v.

We have the relation between (x,y) and (u,v) in Equation (2.3) above. We might explain that the complex variable z is assumed to be made up of a real part, x, and an imaginary part, y, i.e.
    z = x+iy  .  (2.4) 
Likewise, w is assumed to be made up of real and imaginary parts u and v:
   w = u+iv  .  (2.5)

We need to invert Equation (2.3):
   z = w2  .  (2.6)  
Remembering that both z and w are complex, we need to solve this for x in terms of u and v.

There are several ways to do the solution, but perhaps the easiest is just to square Equation (2.5):
   z = w2  = (u + iv)2 = u2v2 + 2iuv  .  (2.7)  
Since Equation (2.7) equals z and the real part of z is defined to be x, the real part of Equation (2.7) therefore equals x, i.e.
   x = u2v2  .  (2.8)  
Substituting this into our equation for Φ, i.e. Equation (2.1), we have:
   Φ = u2v2  .  (2.9)  

In the (u,v) space, the flow velocity will be negative the gradient of this flow potential, i.e.
    velocity = −∇Φ = −∂Φ/∂u au − ∂Φ/∂v av  .  (2.10)  
Substituting (2.9) into (2.10), we have:
   velocity = −2u au + 2v av  .  (2.11)  
Notice that this result indicates the fluid velocity will be in the negative horizontal (i.e. u) direction and in the positive vertical (v) direction, consistent with what we see in Figure 2b above. The equation also indicates that the velocity will be zero in the very corner where both u and v are zero.

fluid flow around the end of a plate Fig. 2c. This figure shows the result of applying the conformal transform:
    w = z2   (2.12)  
to the solution of Laplace's equation shown in Fig. 2a. This transform can also be shown in the animation in Fig. 1 if n in that animation is set to the value 2.

The flow shown in Fig. 2.c is that of a fluid circulating around the end of a horizontal plate that lays along the positive x-axis. The flow starts in the upper right hand quadrant and flows in the negative x-direction until it passes the y-axis at which point it sweeps downward and to the right, exiting in the lower right hand quadrant.

The transformation has rotated the negative x-axis in the counter-clockwise direction all the way so it coincides with the positive x-axis. This is indicated by the red curved arrow. The grid and flow potential have been warped consistent with this. A more accurate detailed mapping of the mapping is shown in a future posting.

Using the animation in Fig. 1, you can play with the more general transformation
    w = zn   (2.13)  
where n varies between 0 and 2. This class of transformations rotates the negative x axis so that it makes an angle with the positive x-axis that is proportional to the power n (actually the angle is given by , i.e. n×3.14 in radians or n×180° in degrees.) Using n larger than 2 usually produces confusing results, since the function is warped around on top of itself.

Mathematics of z2

We now repeat the mathematics we did in the z½ case only now for the z2 case. Like we did in the earlier case, we look for the relation between the w and z coordinate spaces, i.e. between (u,v) and (x,y). This we have in the form of Equation (2.12). We expand this equation out using Equations (2.4) and (2.5):
    u + iv = (x + iy)2  .  (2.14) 
Expanding the right hand side of this yields:
    u + iv = x2y2 + 2ixy   .  (2.15) 
Equating the real part of the left side with the real part of the right side yields:
    u  = x2y2   .  (2.16) 
Repeating with the imaginary parts yields:
    v  = 2xy   .  (2.17) 

Equations (2.16) and (2.17) can be solved for x and y using the quadratic formula to yield:
    x  = 0.707 × [ u + (u2+ v2)½ ]½     (2.18) 

    y  = v/2x   .  (2.19) 

Because our original solution to Laplace's equation was  Φ = x  as given by Equation (2.1), we now just substitute Equation (2.18) in for x to get:
    Φ = 0.707 × [ u + (u2+ v2)½ ]½   .  (2.20) 

To get an equation for the fluid velocity as a function of u and v we need to take minus the gradient of Φ as prescribed by Equation (2.10). This is very straightforward to do, but it yields a rather complicated formula. To make the result simpler to write, we express it in terms of complex polar coordinates rw = (u2+ v2)½ and φw = arctan2(v,u):
    velocityu  = −∂Φ/∂u = (0.354/rw½) (1 + cosφw)    (2.21) 

    velocityv  = −∂Φ/∂v = (0.354/rw½) sinφw / (1 + cosφw)½   .  (2.22) 

Fig. 2d - A sample of other flow patterns from the zn transform
fluid flow in a tight corner fluid flow in a relaxed corner fluid flow around an outside corner
fluid flow diverted by a flat plate fluid flow diverted by an angle 'keep the air clean' cartoon

More on incompressible, irrotational flow

Conformal mapping applies to incompressible, irrotational flow because this type of fluid flow obeys Laplace's equation. Conformal mapping will map one solution of Laplace's equation into another solution. Often we apply a series of conformal mappings to fit complex boundary conditions.

Fluids behave incompressibly and irrotationally if viscosity is not a significant force. Many commonly occurring air and water flow patterns can be approximated as incompressible and irrotational. The air flow around wings of an aircraft traveling considerably less than the speed of sound, fits this criteria, as does the flow around fan and windmill blades. The same is true of airflow around the front 60% of an aerodynamically clean automobile. A previous posting in this series dealt with water waves which involve incompressible, irrotational water flow having a gravitationally controlled free surface.

skydiver cartoon

In the last series of images above we show some other flow patterns obtainable by the zn transform. Note that the last two simply involve starting with the whole space filled with flow, instead of only the upper half space as shown in Fig. 1. Also note that the bending point of the x-axis can be moved by using the (z − x1)n transform in place of zn, where x1 is the point on the x-axis where the bend is to occur.



Copyright, P. Ceperley, September 2010.

Other postings by author up: Conformal mapping contents next: Conformal mapping II