T-I-1: Difference between revisions

From ESPCI Wiki
Jump to navigation Jump to search
Rosso (talk | contribs)
Rosso (talk | contribs)
Line 85: Line 85:


* Consider a constant a uniform flow, parallel to the real axis and with velocity $V_0$. Show that the complex potential writes <math>g_0(z)=V_0 z</math>.
* Consider a constant a uniform flow, parallel to the real axis and with velocity $V_0$. Show that the complex potential writes <math>g_0(z)=V_0 z</math>.
* Consider a fluid in presence of an obstacle. The obstacle is a circle with <math>R=1</math>. Far from the circle the velocity is $V_0$. Use the Joukovski's transformation to show that  the complex potential writes
* Consider a fluid in presence of an obstacle. The obstacle is a circle with <math>R=1</math>. Far from the circle the velocity is <math>V_0</math>. Use the Joukovski's transformation to show that  the complex potential writes
<center><math>g(z)=V_0(z+\frac{1}{z})</math></center>
<center><math>g(z)=V_0(z+\frac{1}{z})</math></center>
* Compute the velocity along the real and the imaginary axis. Draw the streamlines (<math>\psi(z)=\text{const.}</math>) of the flow.
* Compute the velocity along the real and the imaginary axis. Draw the streamlines (<math>\psi(z)=\text{const.}</math>) of the flow.
* Explain (without calculation) how you can use the oukovski's transformation to study the flow in presence of the airfoil.
* Explain (without calculation) how you can use the Joukovski's transformation to study the flow if the circle is replaced by the airfoil.

Revision as of 23:04, 15 October 2011


Analytical functions: conformal map and applications to hydrodynamics

This homework deals with the application of conformal maps to the study of two-dimensional hydrodynamics. A conformal map is a geometrical transformation which preserves all (oriented) crossing angles between lines. In dimension d3 a conformal map is necessarily composed from the following limited number of transformations: translations, rotations, homothetic transformation and special conformal transformation (which is the composition of a reflection and an inversion in a sphere). However in two dimensions, d=2, the space of conformal mappings is much larger and one can show that, given an open set Ω, any holomorphic function f:Ω such that f(z)0, zΩ defines a conformal map from Ω to f(Ω). The aim of this HW is to exploit this property to study some hydrodynamic flows in two spatial dimensions.

Joukovski's transformation

The Joukovski's transformation is defined by the following application

J:{0}zz+1z
  • Compute J(z) and deduce from it the maximal ensemble on which J is a conformal map. Show that J is always injective. Under which condition on the set Ω the application J on Ω is surjective ? Give some examples of such (maximal) ensembles.
  • Give the image by J of the following sub-sets: (a) the half-line passing through the origin O and making an angle α with the x-axis, (b) the circle centered at the origin of radius R (analyse in particular the case R=1). What is the image, by J, of the outside of the unit circle |z|>1.

Hint: it might be useful to use polar coordinates, writing z=reiθ. Get a better idea of this Joukowski's transformation using the following code in Mathematica:

- for the half-line passing through the origin:

 Jouk[z_] := z + 1/z
 Jouk[R Cos[u] + I  R Sin[u]];
 ParametricPlot[{{Re[%], Im[%]} /. {u -> 0.5}, {R Cos[u], R Sin[u]} /. {u -> 0.5}}, {R, .01, 10}]

- for the circle centered at the origin of radius R:

 Jouk[R Cos[u] + I  R Sin[u]];
 ParametricPlot[{{Re[%], Im[%]} /. {R -> 0.79}, {R Cos[u], R Sin[u]} /. {R -> 0.79}}, {u, 0, 2 \[Pi]}, 
 PlotRange -> {{-3, 3}, {-1.5, 1.5}}]


  • Study the conformal map J in the vicinity of z=1: we consider a "smooth" curve γ passing through z=1, with a well defined tangent. Show that the image of γ exhibits a cusp in J(1). In this purpose, we parametrize this curve γ by z(t) with z(0)=1 and z(0)0. Write then the Taylor expansion of z in t=0 up to first order and the expansion of J close to 1 up to second order.


Joukowski showed that the image of a circle passing through z=1 and containing the point z=1 is mapped onto a curve shaped like the cross section of an airplane wing. We call this curve the Joukowski airfoil.

  • Convince yourself that the parametric curve
1R(cos(u)+sin(α))+iR(cos(α)+sin(u))with0<u<2π

identifies a circle of radius R, passing through z=1. α being the angle between the real axis and the tangent at z=1. You can now visualize the Joukowski airfoil using the following code:

Jouk[z_] := z + 1/z
Jouk[1 - R Sin[\[Alpha]] + R Cos[u] +   I (R Cos[\[Alpha]] + R Sin[u])];
ParametricPlot[{{Re[%], Im[%]} /. {R -> 1.15, \[Alpha] -> 1.3}, {{1 - R Sin[\[Alpha]] + R Cos[u]}, {R Cos[\[Alpha]] + R Sin[u]}} /. {R -> 1.15, \[Alpha] -> 1.3}}, {u, 0, 2 \[Pi]}]

Harmonic functions and hydrodynamics in the plane

We recall that a function φ:Ω or (Ω being an open set of ) is called a "harmonic function" if it satisfies the Laplace equation

Δφ=0whereΔφ2φx2+2φy2

in all point z=x+iyΩ. Similarly to conformal maps, harmonic functions in two dimensions, are closely related to holomorphic functions.

  • Let us consider g:Ω a holomorphic function. Show that g,φ=Reg,ψ=Img are harmonic functions.
  • Geometric interpration of φ and ψ: show that the streamlines of φ are the level curves of ψ.
  • Show that, if φ:Ω is a harmonic function and f:ΩΩ a conformal map, then Φ=φf is also a harmonic function.

We now consider the irrotational flow of a non-viscous and incompressible fluid in some region of the plane. We denote by v=(vx,vy) its velocity field.

  • Show that v=(vx,vy) is the gradient of a scalar potential φ(x,y) which satisfies φ=0
  • Show that you can construct ψ(x,y) such that g=φ+iψ is holomorphic and vx+ivy=g(z). g(z) is the complex potential associated to the 2-dimensional fluid flow.

Back to the Joukovski's transformation

  • Consider a constant a uniform flow, parallel to the real axis and with velocity $V_0$. Show that the complex potential writes g0(z)=V0z.
  • Consider a fluid in presence of an obstacle. The obstacle is a circle with R=1. Far from the circle the velocity is V0. Use the Joukovski's transformation to show that the complex potential writes
g(z)=V0(z+1z)
  • Compute the velocity along the real and the imaginary axis. Draw the streamlines (ψ(z)=const.) of the flow.
  • Explain (without calculation) how you can use the Joukovski's transformation to study the flow if the circle is replaced by the airfoil.