Module

for

Sources and Sinks

 

11.11  Sources and Sinks

    If the two-dimensional motion of an ideal fluid consists of an outward radial flow from a point and is symmetrical in all directions, then the point is called a simple source.  A source at the origin can be considered as a line perpendicular to the z plane along which fluid is being emitted.  If the rate of emission of volume of fluid per unit length  [Graphics:Images/SourceSinkMod_gr_1.gif],  then the origin is said to be a source of strength  m,  the complex potential for the flow is  

            [Graphics:Images/SourceSinkMod_gr_2.gif],  

and the velocity  V  at the point [Graphics:Images/SourceSinkMod_gr_3.gif] is given by  

            [Graphics:Images/SourceSinkMod_gr_4.gif].  

For fluid flows, a sink is a negative source and is a point of inward radial flow at which the fluid is considered to be absorbed or annihilated.  Sources and sinks for flows are illustrated in Figure 11.93.

 

            Figure 11.93  Sources and sinks for an ideal fluid.

 

11.11.1  Source: A Charged Line

    In the case of electrostatics a source will correspond to a uniformly charged line perpendicular to the z-plane at the point [Graphics:Images/SourceSinkMod_gr_5.gif].  We will show that if the line L is located at  [Graphics:Images/SourceSinkMod_gr_6.gif]  and carries a charge density of  [Graphics:Images/SourceSinkMod_gr_7.gif]  coulombs per unit length, then the magnitude of the electric field is  [Graphics:Images/SourceSinkMod_gr_8.gif].  Hence  E  is given by  

(11-41)        [Graphics:Images/SourceSinkMod_gr_9.gif],  

and the complex potential is  

            [Graphics:Images/SourceSinkMod_gr_10.gif],  and  
            
            [Graphics:Images/SourceSinkMod_gr_11.gif].  

A sink for electrostatics is a negatively charged line perpendicular to the z-plane.  The electrostatic field for electrostatic problems corresponds to the velocity field for fluid flow problems, except that their corresponding potentials differ by a sign change.

    To establish Equation (11-41), we start with Coulomb's law, which states that two particles with charges q and Q exert a force on one another with magnitude   [Graphics:Images/SourceSinkMod_gr_12.gif],  where r is the distance between particles and C is a constant that depends on the scientific units.  For simplicity, we assume that [Graphics:Images/SourceSinkMod_gr_13.gif] and the test particle at the point z has charge [Graphics:Images/SourceSinkMod_gr_14.gif].  

    The contribution [Graphics:Images/SourceSinkMod_gr_15.gif] induced by the element of charge  [Graphics:Images/SourceSinkMod_gr_16.gif]  along the segment of length [Graphics:Images/SourceSinkMod_gr_17.gif] situated at a height h above the plane has magnitude  [Graphics:Images/SourceSinkMod_gr_18.gif]  given by  

            [Graphics:Images/SourceSinkMod_gr_19.gif].  

    It has the same magnitude as [Graphics:Images/SourceSinkMod_gr_20.gif] induced by the element  [Graphics:Images/SourceSinkMod_gr_21.gif]  located a distance -h  below the plane.  From the vertical symmetry involved their sum,  [Graphics:Images/SourceSinkMod_gr_22.gif],  lies parallel to the plane along the ray from the origin, as shown in Figure 11.94.

 

            Figure 11.94  Contributions to [Graphics:Images/SourceSinkMod_gr_23.gif] from the elements of charge [Graphics:Images/SourceSinkMod_gr_24.gif] situated at [Graphics:Images/SourceSinkMod_gr_25.gif], above and below the z plane.

    By the principal of superposition, we add all contributions from the elements of charge along L to obtain  [Graphics:Images/SourceSinkMod_gr_26.gif].  By vertical symmetry,  [Graphics:Images/SourceSinkMod_gr_27.gif] lies parallel to the complex plane along the ray from the origin through the point z.  Hence the magnitude of  [Graphics:Images/SourceSinkMod_gr_28.gif] is the sum of all components  [Graphics:Images/SourceSinkMod_gr_29.gif]  that are parallel to the complex plane, where t is the angle between  [Graphics:Images/SourceSinkMod_gr_30.gif] and the plane.  Letting  [Graphics:Images/SourceSinkMod_gr_31.gif]  in this summation process produces the definite integral   

            [Graphics:Images/SourceSinkMod_gr_32.gif].  

    Next, we use the change of variable  [Graphics:Images/SourceSinkMod_gr_33.gif]  and the trigonometric identity  [Graphics:Images/SourceSinkMod_gr_34.gif]  to obtain the equivalent integral:  

            [Graphics:Images/SourceSinkMod_gr_35.gif].  

    Multiplying this magnitude  [Graphics:Images/SourceSinkMod_gr_36.gif]  by the unit vector  [Graphics:Images/SourceSinkMod_gr_37.gif]  establishes Formula (11-41).  If [Graphics:Images/SourceSinkMod_gr_38.gif], then the field is directed away from [Graphics:Images/SourceSinkMod_gr_39.gif] and, if [Graphics:Images/SourceSinkMod_gr_40.gif], then it is directed toward [Graphics:Images/SourceSinkMod_gr_41.gif].  An electrical field located at [Graphics:Images/SourceSinkMod_gr_42.gif] is given by   

            [Graphics:Images/SourceSinkMod_gr_43.gif],  

and the corresponding complex potential is  

            [Graphics:Images/SourceSinkMod_gr_44.gif].

 

Example 11.30.  (Source and sink of equal strength)  Let a source and sink of unit strength be located at the points [Graphics:Images/SourceSinkMod_gr_45.gif] and [Graphics:Images/SourceSinkMod_gr_46.gif], respectively.  The complex potential for a fluid flowing from the source at  [Graphics:Images/SourceSinkMod_gr_47.gif]  to the sink at  [Graphics:Images/SourceSinkMod_gr_48.gif]  is  

            [Graphics:Images/SourceSinkMod_gr_49.gif].  

[Graphics:Images/SourceSinkMod_gr_50.gif]     [Graphics:Images/SourceSinkMod_gr_51.gif]     

            Figure 11.95 (a)  Source and sink of unit strength.

The velocity potential and stream function are  

            [Graphics:Images/SourceSinkMod_gr_52.gif],   and  

            [Graphics:Images/SourceSinkMod_gr_53.gif],  

respectively.  Solving for the streamline  [Graphics:Images/SourceSinkMod_gr_54.gif],  we start with   

            [Graphics:Images/SourceSinkMod_gr_55.gif]  

and obtain the equation

            [Graphics:Images/SourceSinkMod_gr_56.gif].

A straightforward calculation shows that points on the streamline must satisfy the equation  

            [Graphics:Images/SourceSinkMod_gr_57.gif],  

which is the equation of a circle with center at  [Graphics:Images/SourceSinkMod_gr_58.gif]  that passes through the points  [Graphics:Images/SourceSinkMod_gr_59.gif].  Several streamlines are indicated in Figure 11.95(a).

Explore Solution 11.30.

 

Example 11.31.  (Two sources of equal strength)  Let two sources of unit strength be located at the points  [Graphics:Images/SourceSinkMod_gr_84.gif].  The resulting complex potential for a fluid flow is  

            .  

 

[Graphics:Images/SourceSinkMod_gr_86.gif]     [Graphics:Images/SourceSinkMod_gr_87.gif]     

                 Figure 11.95 (b)  Two sources of unit strength.

The velocity potential and stream function are


            [Graphics:Images/SourceSinkMod_gr_88.gif],   and  

            [Graphics:Images/SourceSinkMod_gr_89.gif],  

respectively.  Solving for the streamline  [Graphics:Images/SourceSinkMod_gr_90.gif],  we start with  

            [Graphics:Images/SourceSinkMod_gr_91.gif]

and obtain the equation  

            [Graphics:Images/SourceSinkMod_gr_92.gif].  

     If we express this equation in the form  [Graphics:Images/SourceSinkMod_gr_93.gif]  then we can get  

            [Graphics:Images/SourceSinkMod_gr_94.gif],   or  
            
            [Graphics:Images/SourceSinkMod_gr_95.gif].

Now use the rotation of axes  

            [Graphics:Images/SourceSinkMod_gr_96.gif]  

then the streamlines must satisfy the equation   

            [Graphics:Images/SourceSinkMod_gr_97.gif]  

and are rectangular hyperbolas with centers at the origin that pass through the points  [Graphics:Images/SourceSinkMod_gr_98.gif].  Several streamlines are indicated in Figure 11.95(b).

Explore Solution 11.31.

 

 

Ideal Fluid Flow

    Let an ideal fluid flow in a domain in the z plane be effected by a source located at the point [Graphics:Images/SourceSinkMod_gr_121.gif].  Then the flow at points z, which lie in a small neighborhood of the point [Graphics:Images/SourceSinkMod_gr_122.gif], is approximated by that of a source with the complex potential  

            [Graphics:Images/SourceSinkMod_gr_123.gif].  

If  [Graphics:Images/SourceSinkMod_gr_124.gif]  is a conformal mapping and  [Graphics:Images/SourceSinkMod_gr_125.gif],  then S(z) has a nonzero derivative at [Graphics:Images/SourceSinkMod_gr_126.gif] and   

            [Graphics:Images/SourceSinkMod_gr_127.gif]

where  [Graphics:Images/SourceSinkMod_gr_128.gif]  as  [Graphics:Images/SourceSinkMod_gr_129.gif].  Taking logarithms yields

            [Graphics:Images/SourceSinkMod_gr_130.gif].

Because  [Graphics:Images/SourceSinkMod_gr_131.gif],  the term  [Graphics:Images/SourceSinkMod_gr_132.gif]  approaches the constant value  [Graphics:Images/SourceSinkMod_gr_133.gif]  as  [Graphics:Images/SourceSinkMod_gr_134.gif].  As  [Graphics:Images/SourceSinkMod_gr_135.gif]  is the complex potential for a source located at the point [Graphics:Images/SourceSinkMod_gr_136.gif], the image of a source under a conformal mapping is a source.

    We can use the technique of conformal mapping to determine the fluid flow in a domain D in the z plane that is produced by sources and sinks.  If we can construct a conformal mapping  [Graphics:Images/SourceSinkMod_gr_137.gif]  so that the image of sources, sinks, and boundary curves for the flow in D are mapped onto sources, sinks, and boundary curves in a domain G where the complex potential is known to be [Graphics:Images/SourceSinkMod_gr_138.gif], then the complex potential in D is given by [Graphics:Images/SourceSinkMod_gr_139.gif].

 

Example 11.32.  Suppose that the lines  [Graphics:Images/SourceSinkMod_gr_140.gif]  are considered as walls of a containing vessel for a fluid flow produced by a single source of unit strength located at the origin.  The conformal mapping  [Graphics:Images/SourceSinkMod_gr_141.gif]   maps the infinite strip bounded by the lines  [Graphics:Images/SourceSinkMod_gr_142.gif]  onto the w plane slit along the boundary rays  [Graphics:Images/SourceSinkMod_gr_143.gif]  and  [Graphics:Images/SourceSinkMod_gr_144.gif],  and the image of the source at  [Graphics:Images/SourceSinkMod_gr_145.gif]  is a source located at  [Graphics:Images/SourceSinkMod_gr_146.gif].  The complex potential  

            [Graphics:Images/SourceSinkMod_gr_147.gif]  

determines a fluid flow in the w plane past the boundary curves  [Graphics:Images/SourceSinkMod_gr_148.gif]  and  [Graphics:Images/SourceSinkMod_gr_149.gif], which lie along streamlines of the flow.  Therefore the complex potential for the fluid flow in the infinite strip in the z plane is  

            [Graphics:Images/SourceSinkMod_gr_150.gif].  

Several streamlines for the flow are illustrated in Figure 11.96. 

[Graphics:Images/SourceSinkMod_gr_151.gif]     [Graphics:Images/SourceSinkMod_gr_152.gif]     

                 Figure 11.96  A source in the center of a strip.

Explore Solution 11.32.

 

Example 11.33.  Suppose that the lines  [Graphics:Images/SourceSinkMod_gr_170.gif]  are considered as walls of a containing vessel for the fluid flow produced by a single source of unit strength located at the point  [Graphics:Images/SourceSinkMod_gr_171.gif]  and a sink of unit strength located at the point  [Graphics:Images/SourceSinkMod_gr_172.gif].  The conformal mapping  [Graphics:Images/SourceSinkMod_gr_173.gif]  maps the infinite strip bounded by the lines  [Graphics:Images/SourceSinkMod_gr_174.gif]  onto the w plane slit along the boundary rays  [Graphics:Images/SourceSinkMod_gr_175.gif]  and  [Graphics:Images/SourceSinkMod_gr_176.gif].  The image of the source at [Graphics:Images/SourceSinkMod_gr_177.gif] is a source at [Graphics:Images/SourceSinkMod_gr_178.gif], and the image of the sink at [Graphics:Images/SourceSinkMod_gr_179.gif] is a sink at [Graphics:Images/SourceSinkMod_gr_180.gif].  The potential  

            [Graphics:Images/SourceSinkMod_gr_181.gif]  

determines a fluid flow in the w plane past the boundary curves [Graphics:Images/SourceSinkMod_gr_182.gif], which lie along streamlines of the flow.  Therefore the complex potential for the fluid flow in the infinite strip in the z plane is

            [Graphics:Images/SourceSinkMod_gr_183.gif].  

Several streamlines for the flow are illustrated in Figure 11.97. 

     [Graphics:Images/SourceSinkMod_gr_185.gif]     

                 Figure 11.97  A source and a sink on the edges of a strip.

Explore Solution 11.33.

 

    We can use the technique of transformation of a source to determine the effluence from a channel extending from infinity.  In this case, we construct a conformal mapping  [Graphics:Images/SourceSinkMod_gr_205.gif]  from the upper half-plane [Graphics:Images/SourceSinkMod_gr_206.gif] so that the single source located at [Graphics:Images/SourceSinkMod_gr_207.gif] is mapped to the point [Graphics:Images/SourceSinkMod_gr_208.gif] at infinity that lies along the channel.  The streamlines emanating from [Graphics:Images/SourceSinkMod_gr_209.gif] in the upper half-plane are mapped onto streamlines issuing from the channel.

 

Example 11.34.  Consider the conformal mapping   

            [Graphics:Images/SourceSinkMod_gr_210.gif],  

which maps the upper half-plane [Graphics:Images/SourceSinkMod_gr_211.gif] onto the domain consisting of the upper half-plane [Graphics:Images/SourceSinkMod_gr_212.gif] joined to the channel  [Graphics:Images/SourceSinkMod_gr_213.gif].  The point [Graphics:Images/SourceSinkMod_gr_214.gif] is mapped onto the point [Graphics:Images/SourceSinkMod_gr_215.gif] and [Graphics:Images/SourceSinkMod_gr_216.gif] is mapped onto [Graphics:Images/SourceSinkMod_gr_217.gif] along the channel.  Images of the rays  [Graphics:Images/SourceSinkMod_gr_218.gif], where [Graphics:Images/SourceSinkMod_gr_219.gif] is a constant, are streamlines issuing from the channel as indicated in Figure 11.98. 

[Graphics:Images/SourceSinkMod_gr_220.gif]          

                 Figure 11.98  Effluence from a channel into a half-plane.

Explore Solution 11.34.

 

Extra Example 1.  Use a Schwarz-Christoffel transformation, (see Section 11.9), to show that the conformal mapping  

            [Graphics:Images/SourceSinkMod_gr_229.gif]  

will map the flow in the upper half-plane from a source at  [Graphics:Images/SourceSinkMod_gr_230.gif]  onto the flow from a channel with  [Graphics:Images/SourceSinkMod_gr_231.gif]  into the first quadrant, as shown in Figure 11.105.
Hint: Set  [Graphics:Images/SourceSinkMod_gr_232.gif] and [Graphics:Images/SourceSinkMod_gr_233.gif].   

[Graphics:Images/SourceSinkMod_gr_234.gif]     

                 Figure 11.105  Flow from a channel into a quadrant.

Solution.  Using the Schwarz-Christoffel transformation method, (see Section 11.9), the exterior angles are  [Graphics:Images/SourceSinkMod_gr_235.gif],  and the formula for the derivative [Graphics:Images/SourceSinkMod_gr_236.gif] is  

            [Graphics:Images/SourceSinkMod_gr_237.gif]   
integrate

            [Graphics:Images/SourceSinkMod_gr_238.gif]
            
            [Graphics:Images/SourceSinkMod_gr_239.gif]
            
            [Graphics:Images/SourceSinkMod_gr_240.gif]

    
use the conditions  [Graphics:Images/SourceSinkMod_gr_241.gif],  [Graphics:Images/SourceSinkMod_gr_242.gif]  and solve the resulting system for A and B.  The desired result is

            [Graphics:Images/SourceSinkMod_gr_243.gif].

Explore Extra Solution 1.

 

Extra Example 2.  Use a Schwarz-Christoffel transformation, (see Section 11.9), to show that the conformal mapping  

            [Graphics:Images/SourceSinkMod_gr_259.gif]
            

will map the flow in the upper half-plane from a source at  [Graphics:Images/SourceSinkMod_gr_260.gif]  onto the flow from a channel with  [Graphics:Images/SourceSinkMod_gr_261.gif]  into a sector, as shown in Figure 11.106.
Hint: Set  [Graphics:Images/SourceSinkMod_gr_262.gif] and [Graphics:Images/SourceSinkMod_gr_263.gif].  Use the change of variable  [Graphics:Images/SourceSinkMod_gr_264.gif]  in the resulting integral. 

[Graphics:Images/SourceSinkMod_gr_265.gif]     

                 Figure 11.106  Flow from a channel into a sector in the first quadrant.

Solution.  Using the Schwarz-Christoffel transformation method, (see Section 11.9), the exterior angles are  [Graphics:Images/SourceSinkMod_gr_266.gif],  and the formula for the derivative [Graphics:Images/SourceSinkMod_gr_267.gif] is  

            [Graphics:Images/SourceSinkMod_gr_268.gif]   
integrate

            [Graphics:Images/SourceSinkMod_gr_269.gif]
            
use the substitutions

            [Graphics:Images/SourceSinkMod_gr_270.gif]  
and get

            [Graphics:Images/SourceSinkMod_gr_271.gif]  

Use the identity  [Graphics:Images/SourceSinkMod_gr_272.gif]  and write the integral as

            [Graphics:Images/SourceSinkMod_gr_273.gif]

Now use the substitution [Graphics:Images/SourceSinkMod_gr_274.gif] and get

            [Graphics:Images/SourceSinkMod_gr_275.gif]
            
use the conditions  [Graphics:Images/SourceSinkMod_gr_276.gif],  [Graphics:Images/SourceSinkMod_gr_277.gif]  and solve the resulting system for A and B.  The desired result is

            [Graphics:Images/SourceSinkMod_gr_278.gif].

Explore Extra Solution 2.

 

Extra Example 3.  Use a Schwarz-Christoffel transformation, (see Section 11.9), to show that the conformal mapping  

            [Graphics:Images/SourceSinkMod_gr_292.gif]  

will map the flow in the upper half-plane from a source at  [Graphics:Images/SourceSinkMod_gr_293.gif]  onto the flow in a right-angled channel with  [Graphics:Images/SourceSinkMod_gr_294.gif], as shown in Figure 11.107. 

[Graphics:Images/SourceSinkMod_gr_295.gif]     

                 Figure 11.107  Flow in a right-angled channel.

Solution.  Using the Schwarz-Christoffel transformation method, (see Section 11.9), if we choose  [Graphics:Images/SourceSinkMod_gr_296.gif]  and  [Graphics:Images/SourceSinkMod_gr_297.gif],  then the formula  

            [Graphics:Images/SourceSinkMod_gr_298.gif]

will determine a mapping [Graphics:Images/SourceSinkMod_gr_299.gif] of the upper half-plane onto the domain indicated in Figure 11.74(a).  With  [Graphics:Images/SourceSinkMod_gr_300.gif],  we let [Graphics:Images/SourceSinkMod_gr_301.gif],  then  [Graphics:Images/SourceSinkMod_gr_302.gif]  and  [Graphics:Images/SourceSinkMod_gr_303.gif].  

    The limiting formula for the derivative [Graphics:Images/SourceSinkMod_gr_304.gif] becomes   

            [Graphics:Images/SourceSinkMod_gr_305.gif]

where  [Graphics:Images/SourceSinkMod_gr_306.gif], which will determine a mapping [Graphics:Images/SourceSinkMod_gr_307.gif] from the upper half plane onto the channel as indicated in Figure 11.74(b) and Figure 11.107.

    Using integrals in Table 11.2, we obtain

        [Graphics:Images/SourceSinkMod_gr_308.gif]   

    If we use the principal branch of the inverse sine function, then the boundary values  [Graphics:Images/SourceSinkMod_gr_309.gif]  lead to the system  

            [Graphics:Images/SourceSinkMod_gr_310.gif]   and   [Graphics:Images/SourceSinkMod_gr_311.gif],  

which we can solve to obtain [Graphics:Images/SourceSinkMod_gr_312.gif].  Hence the required solution is

            [Graphics:Images/SourceSinkMod_gr_313.gif]  

Explore Extra Solution 3.

 

Extra Example 4.  Use a Schwarz-Christoffel transformation, (see Section 11.9), to show that the conformal mapping  

            [Graphics:Images/SourceSinkMod_gr_322.gif]  

will map the flow in the upper half-plane onto the flow from a channel back into a quadrant, as shown in Figure 11.108,
where  [Graphics:Images/SourceSinkMod_gr_323.gif]  and  [Graphics:Images/SourceSinkMod_gr_324.gif]

[Graphics:Images/SourceSinkMod_gr_325.gif]     

                 Figure 11.108  Flow from a channel  back into the first quadrant.

Explore Extra Solution 4.

 

Library Research Experience for Undergraduates

Electrostatics

Ideal Fluid Flow

Steady State Temperature

Joukowski Transformation and Airfoils

Schwarz-Christoffel Transformation

Complex Potential

Green's Function

 

Download This Mathematica Notebook

  Download The Maple Worksheet

   

The Next Module is

 Fourier Series

 

 

Return to the Complex Analysis Modules  

 

 

Return to the Complex Analysis Project

 

 

(c) 2006 John H. Mathews, Russell W. Howell