{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "direct-toner",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "# Recap of Fundamentals in Aerodynamics - 2D and 3D"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "exact-dividend",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "```{epigraph}\n",
    "... if we all designed bodies to have stagnation points at the back, then D'Alembert would rest happy, and we would have no drag. Presumably that would be a Good Thing generally, except for parachutists!\n",
    "```\n",
    "<p style=\"text-align:right\" > <a href=\"https://arc.aiaa.org/doi/abs/10.2514/6.1996-161\"> PBS Lissaman</a> </p>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "periodic-humanity",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "```{epigraph}\n",
    "It’s easy to explain how a rocket works, but explaining how a wing works takes a rocket scientist.\n",
    "```\n",
    "<p style=\"text-align:right\" > <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\"> Philippe Spalart</a> </p>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "advance-slovenia",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "```{epigraph}\n",
    "Aerodynamic theory is now (1930) rather like the physical theory of light; Sir William Bragg recently said that physicists use the electron theory on Mondays, Wednesdays and Fridays, and the wave theory on alternate days. Both have uses but reconciliation of the two ideas has not yet been achieved. So it is in aeronautics. In our experimental work we assume that viscosity is an essential property of air and the building of a compressed-air tunnel is the latest expression of that belief. The practically useful theory of Prandtl comes from considering air as frictionless or inviscid.\n",
    "```\n",
    "<p style=\"text-align:right\" > <a href=\"https://press.uchicago.edu/ucp/books/book/chicago/E/bo12024403.html\"> Leonard Bairstow</a> </p>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "strange-matter",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "```{epigraph}\n",
    "To those who fear flying, it is probably disconcerting that physicists and aeronautical engineers still passionately debate the fundamental issue underlying this endeavor: what keeps planes in the air?\n",
    "```\n",
    "<p style=\"text-align:right\" > <a href=\"https://www.nytimes.com/2003/12/09/news/staying-aloft-what-does-keep-them-up-there.html\"> STAYING ALOFT; What Does Keep Them Up There? </a></p>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "floating-unknown",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "```{epigraph}\n",
    "On a strictly mathematical level, engineers know how to design planes that will stay aloft. But equations don't explain why aerodynamic lift occurs. There are two competing theories that illuminate the forces and factors of lift. Both are incomplete explanations. Aerodynamicists have recently tried to close the gaps in understanding. Still, no consensus exists.\n",
    "```\n",
    "<p style=\"text-align:right\" > <a href=\"https://www.scientificamerican.com/article/no-one-can-explain-why-planes-stay-in-the-air/\"> Scientific American</a> </p>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "devoted-consumer",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "This chapter is a whirlwind tour of some of the basics of airfoil aerodynamics that is relevant for this course. Since airfoils can be imagined as the basic building blocks that make up the rotor, the fundamental understanding of airfoil aerodynamics laid down here would help develop a more intuitive understanding of rotor aerodynamics. As some of the quotes above imply, the mathematical theory of fluid behavior, and by extension the mathematiacl theory of flight, are well established. However, given the complex equations involved a straightforward cause and effect relation is difficult to establish and hence a simple intuitive explanation for something even as common and fundamental as the process of lift generation on an airfoil is still missing. Nevertheless, below we try to make as much progress as possible in the improving/refreshing the understanding of the fundamental fluid dynamics concepts involved. It is hoped that the discussions in this chapter will make the fundamental concepts, underlying much of fluid dynamics, more intuitive than the involved mathematical equations themselves.\n",
    "\n",
    "This chapter is predominantly based on the content in {cite}`BA3`, in particular Chapter 7 from this reference. Additionally, {cite}`A23` and {cite}`A23_1` provide a good overview of the historical development of the field of fluid dynamics. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "canadian-semiconductor",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Bodies in fluids"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "prompt-paraguay",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Understanding the nature of forces acting on bodies of arbitrary shape is in our interest if the objective is to minimise the amount of work done in order to overcome the resistive forces due to friction caused due to relative motion between the body and the surrounding fluid. Optimum hull design of a cruise ship, minimum drag profile of a competitive road cyclist+cycle etc are a few examples where an understanding of fluid behavior would help achieve the objectives. Of course, we need not look too far for examples because a helicopter has various entities that might benefit from design based on a sound knowhow in fluid dynamics. A streamlined fuselage design such that drag (rotorcraft fuselage drag is referred to as parasite drag) is minimised would go a long way in reducing operating costs.  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "adjacent-working",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Some fundamental laws are followed by real fluids -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "prepared-virus",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "- conservation laws\n",
    "    - energy, linear momentum and mass"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "removed-exception",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "- they cannot penetrate the body surface ($F$)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "overall-packet",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "$$\\mathbf{v}. \\nabla F = 0$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8255182f",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Fluid cannot penetrate the surface of the body in consideration. This seems like a trivial concept to state explicitly but remember that when a problem involving fluids is mathematically stated, the condition at the boundaries of the fluid need to be clearly defined. In the case of an airfoil moving about in *infinite* fluid, the surface of the airfoil forms one of the boundaries (the other boundary is at infinite distance away and not a source of any disturbance). What happens to the fluid at this boundary, i.e. the surface of the airfoil, dictates what happens in the rest of the fluid. This behavior of the fluid, specifically the velocity in the case of an incompressible flow, is the quantity of interest because this would then eventually lead to the amount of lift generated from the subsequent behavior of the fluid due to the airfoil motion. Since we know from physical observations that fluid parcels do not permeate through object surfaces, the motion of the fluid at this boundary is exactly matched with the motion of the airfoil surface. This is refered to as the no penetration boundary condition."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "formal-beverage",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Any model of (Newtonian) fluids would have to satisfy the above laws and we do have a model that does that - the Navier-Stokes equations. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "mounted-edinburgh",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Navier-Stokes equations (NSE)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "viral-parks",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The NSE are essential to so many analyses you encounter involving fluid flows yet they can be solved without the use of a computer for only a handful of idealised scenarios. Even today, the fact that these equations are challenging to solve is what keeps many aerospace and mechanical engineers employed. Which is why in the early days of the invention of the airplane, and essentially the beginning of the field of aerodynamics, the exact mathematical theory to analyse the flow already existed (NSE had already been around for atleast half a century by the time of the first Wright brothers' flight) but was of little help since the equations couldn't be solved for the body of interest - a wing! Clever simplifications had to be devised based on meticulous experimental observations and astute engineering judgement so that the mathematical problem of lift generation could be posed in a way that it was atleast solvable, while of course simultaneously leading to results of acceptable accuracy. The schematic below shows the simplifications that are often adopted in fluid analyses based on the problem at hand.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "subtle-coordinator",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/FluidEqsChart2.svg\" class=\"large\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Fluid models and equations with different approximations  </p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "productive-hormone",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "The assumptions that lead to simplifications are the following-"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "circular-terrorism",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "- no viscosity"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "scenic-recipient",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Viscosity is a property of the fluid. Given the low kinematic viscosity of air, the viscous terms can be assumed to be much smaller than the other terms in certain scenarios leading to the Euler equation. However, even the Euler equation cannot be solved analytically. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "adjacent-swiss",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "- incompressibility"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "pursuant-winner",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "assuming that the air in incompressible registers an error in the density of only about 3% for flow velocity upto Mach number 0.3. Clearly, for low speed applications this assumption is a reasonable one."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "meaningful-field",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "- irrotationality"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "timely-strand",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Irrotationality is a property of the flow (and not the fluid itself). The influence of this assumption is mathematically very apparent in that it leads to the fluid flow behavior to be represented using the Laplace's equation. Which is an incredible simplification because of the properties that the function $\\phi$, that statisfies the Laplace equation, is bound to exhibit. This is a linear equation which means it is amenable to superposition - i.e. flowfield around two bodies placed close to each other is the sum of the flow fields of the two bodies in isolation. The function $\\phi$ is continuous and differentiable (referred to as holomorphic in mathematical parlance) which means **any differentiable function represents the flow over some body!** We just need to identify the function that represents the flowfield over the body of interest to us - an airfoil or a flat plate (airfoil with zero thickness). Lastly, and from the point of view of early practioners in the field of aerodynamics, the ubiquity with which this function appears in different areas of physics - electrostatics, gravitation etc. - means that aerodynamiscists would not be alone in their struggle to identify solution strategies for the Laplace's equation."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "lightweight-receipt",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<div style=\"float:right; width:45%;\" >\n",
    "    <img src=\"../assets/MT_heirarchy_noannot.png\" class=\"big\" />\n",
    "    <p style=\"text-align:center; color:gray; font-size:60%;\"> Heirarchy of physical theories [taken from \"Understanding Aerodynamics\" by Doug McLean] <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\">[source] </a>\n",
    "  </p>\n",
    "</div>\n",
    "\n",
    "Principle aim of a model\n",
    "- include *just* relevant behavior of physical systems\n",
    "- superfluousness (generally) increases computational complexity\n",
    "- (too much) parsimony can make it *useless*\n",
    "\n",
    "NSE are an *approximate* model of fluids\n",
    "- rarely do these *approximations* become limitations \n",
    "- fluids modelled as a continuous medium \n",
    "    - *granularity* (individual molecules) is ignored\n",
    "- basic conservation laws are satisfied \n",
    "    - mass, momentum, energy\n",
    "- constitutive relations between involved variables\n",
    "    - to make as many equations as there are variables\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "local-geneva",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Simple governing laws, complex results"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "plastic-broadway",
   "metadata": {
    "cell_style": "center",
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The NSE are generally all that an aerospace engineer needs for all their aerodynamics analyses needs - five field partial differential equations and three constitutive relations to solve for three space coordinates (Lagrangian)/velocity components (Eulerian), pressure, density, temperature, viscosity and thermal conductivity. fuselage [^RotorWake] Being field equations, the NSE can be solved if the boundary conditions of the domain of investigation are known. In case of boundaries that form the interface with another material - for e.g. the surface of a wing - additional conditions need to be imposed since NSE do not cover the physics at such interfaces. Examples include \"no slip\" boundary condition typically employed for flow over solid boundaries and \"no temperature jump\" in case of a temperature differential.\n",
    "  \n",
    "Since the momentum conservation equations are part of the NSE, Newton's 2nd is definitely involved when describing flow behavior. However, simply attributing lift generated by an airfoil, for e.g., to just Newton's 2nd law is not entirely accurate. There is an obvious cause and effect relation when one talks about Newton's 2nd law. Within a fluid, however, the motion of continuum parcels is influenced by the forces exerted and the forces exerted are in turn influenced by the motion of these parcels. For e.g., the momentum equation above can be viewed as a form of Newton's second law of motion F=mass x acceleration, except that the forcing terms on the RHS of the equation - the pressure gradient and the shear stresses - are themselves a function of the fluid velocity. The complex interaction between the three sets of forces - inertial, pressure and shear - is what gives fluid dynamics its richness.\n",
    "\n",
    "Such implicit relations of relatively simple local physics (conservation laws) buried within the NSE make an intuitive foresight of what would happen in the global flowfield difficult. Below is an informal discussion of the mechanism by which airfoils generate lift to exemplify this. \n",
    "\n",
    "[^RotorWake]: When modelling rotor wake there is not much to be gained from solving the most complex form of the problem. For e.g., most rotorcraft (atleast conventional designs) don't fly faster than M=0.3 so right away density can be assumed to be constant without incurring much error. So can temperature and viscosity, and thermal conductivity is just zero. From the continuity and the momentum conservation equations, the problem then boils down to solving for the three velocity components and the pressure distribution. More on that later.  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "gorgeous-rover",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "<div class=\"image-grid2\">\n",
    "    <img  src=\"../assets/airfoil_flow_newton.png\" class=\"small\"/>\n",
    "    <img  src=\"../assets/airfoil_flow_real.png\" class=\"small\"/>\n",
    "</div>\n",
    "\n",
    "<div class=\"image-grid2\">\n",
    "  <p > Streamlines representing hypothetical fluid flow idealised as independent fluid 'particles' (Newton's approach)\n",
    "  </p> \n",
    "  <p > Streamlines representing real fluid flow\n",
    "  </p>       \n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "judicial-artist",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Some of the salient features of fluid behavior have been summarised below -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "minor-happiness",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "- fluid (infinitely deformable) follows airfoil contour\n",
    "    - a result of rapid random motion of molecules\n",
    "- this influences fluid velocity in a large region in the vicinity of airfoil\n",
    "    - natural by-product of fluids behaving as continuous media     \n",
    "- flow sped up on upper surface and slowed down on lower surface\n",
    "    - so Bernoulli-based explanation of lift is correct (but not sufficient!)     \n",
    "- most lift occurs due to suction pressure on upper side   \n",
    "- airfoil shape has a significant effect on lift generated\n",
    "    - this is completely ignored by Newton's approach\n",
    "- mutual/reciprocal relation between pressure and velocity\n",
    "    - pressure difference causes fluid acceleration\n",
    "    - existence of pressure difference due velocity difference \n",
    "    - inertia plays central role\n",
    "        - air is *light* but at fast enough speeds, large enough mass can be influenced\n",
    "- magnitude of lift generated can be controlled \n",
    "    - linear dependence in $\\alpha$\n",
    "    - square dependence in fluid velocity\n",
    "        - amount of fluid influenced \n",
    "        - downward momentum imparted\n",
    "- the above effects die down with distance from the airfoil surfaces  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "unusual-replacement",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "```{admonition} Aside\n",
    "If the picture on the left above is assumed to be the *true* fluid model, then it is easy to see that the lift generated would be proportional to square of the angle of attack. However, we now know that airfoil lift is proportional to angle of attack (at least when the flow is not significantly separated). This is also evident from the pressure contours around an airfoil shown below - doubling the anlge of attack doubles the pressure difference, the amount of fluid influenced stays the same.\n",
    "<div class=\"image-grid2\">\n",
    "    <img  src=\"../assets/BA3_fig7_3_24a.png\" class=\"medium\"/>\n",
    "    <img  src=\"../assets/BA3_fig7_3_24b.png\" class=\"medium\"/>\n",
    "</div>\n",
    "\n",
    "<div class=\"image-grid2\"> \n",
    "  <p > Pressure field around NACA0010 at angle of attack of 4° <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\">[source] </a>\n",
    "  </p> \n",
    "  <p > Pressure field around NACA0010 at angle of attack of 8° <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\">[source] </a>\n",
    "  </p>       \n",
    "</div>\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "professional-bennett",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Notice that viscosity did not come up in explanation above. Viscosity plays a key role in the ability of airfoil sections to produce lift in that it ensures that the fluid doesn't turn around the trailing edge. However, it is not necessary in order to quantify the amount of lift which is why simple potential flow theory models, such as thin airfoil theory, are capable to predicting airfoil lift to a good degree of accuracy. Even when the flow conditions are unsteady, potential flow theory-based unsteady models are just as capable in predicting the unsteady lift and moment variation.  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "polar-single",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Generation of lift"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "rising-worcester",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "\n",
    "\"Lift is an interaction in which the airfoil and the fluid exchange equal and opposite forces (as per Newton's third law)\" - {cite}`BA3`."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "spread-touch",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "<div style=\"float:right; width:25%;\" >\n",
    "    <img src=\"../assets/AA_C1_CP_contour.png\" class=\"big\" />\n",
    "    <p style=\"text-align:center; color:gray; font-size:60%;\"> Coefficient of pressure contour plot of NACA23012 airfoil at $M=0.4$ and $\\alpha = 2^°$ <a href=\"https://arc.aiaa.org/doi/pdf/10.2514/6.2019-1101#:~:text=The%20aerodynamic%20efficiency%20was%20evaluated,%3D%200.6%20at%20Cl%20%3D%200.85.\">[source] </a>\n",
    "  </p>\n",
    "</div>\n",
    "\n",
    "- mathematically\n",
    "    - its because of the NSE\n",
    "- intuitively\n",
    "    - ?\n",
    "- always accompanied by circulation and vorticity\n",
    "- manifests as pressure difference \n",
    "- envisaged as flow turning \n",
    "    - below the airfoil AND above it!!\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "formed-check",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA8_fig5_2.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Pressure distribution around an airfoil that contributes to drag (left) and lift (right). Red arrows show positive contribution and blue arrows negative contribution.  <a href=\"https://mitpress.mit.edu/9780262526449/flight-vehicle-aerodynamics/\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "deluxe-transsexual",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Vorticity and circulation"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "marine-malpractice",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Vorticity $\\mathbf{\\omega}$ is a measure of the local rotation of fluid elements i.e. the element rotation about its own axis. If the velocity flow field $\\mathbf{V}$ is known at every point then mathematically this 'rotation' can be quantified using the following equation -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "functional-appraisal",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "$$\\mathbf{\\omega} = \\nabla \\times \\mathbf{V}$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "sonic-theory",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA3_fig3_3_5.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Vorticity of a continuum particle/element (modified from <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\">[source]) </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "controlled-richards",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Fluid dynamics, as it concerns aerospace engineering, is largely about keeping track of vorticity generation and vorticity transport. It is introduced into the flow due to no-slip condition at the airfoil surface (or the surface of any body for that matter). This vorticity doesn‘t just immediately dissipate/dissapear and keeping track of it is imperative to get overall aerodynamics right. Note that viscosity is not necessary in order to generate vorticity, based on the flow profiles shown above if *somehow* a velocity shear exists between two layers of fluid then we know from the definition that a sheet of vorticity exists between those two layers and viscosity does not even come into the picture. \n",
    "\n",
    "Recall that the velocity profile within the boundary layer on the upper surface looks like that shown in figure (c) above, where the lower fluid layers always have a lower velocity than the ones above them leading to a shearing effect and consequently generate vorticity in the fluid. The figure below shows varying degrees of modelling fidelity where the vorticity distributed within the fluid volume in the boundary layer is approximated as lumped as a sheet over the surface of the airfoil itself or lumped over a line (usually the camber line of the airfoil). Alternately, the entire distributed circulatory effect of the airfoil can also be lumped at just one point on the airfoil (usually the quarter-chord). With each simplification the accuracy of the velocity field prediction is affected, especially close (within a few chord lengths) to the airfoil. However, for most applications (including rotary-wing analysis) representing the circulation around the airfoil as lumped at the quarter chord location suffices and modelling the wake of the blade (finite-wing) using Prandtl's lifting-line theory becomes more relevant for accurate analysis.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "confidential-maple",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA8_fig2_12.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Schematic of the actual vorticity distribution around an airfoil due to the boundary layer and in the wake compared in real flow versus a simplified model  <a href=\"https://mitpress.mit.edu/9780262526449/flight-vehicle-aerodynamics/\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fixed-adult",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA3_BA1_mix.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Schematic showing the origins of flow vorticity in the boundary layer and further levels of modelling fidelity (modified from <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\">[source]) </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "distinguished-version",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The circulation of the flow is a relevant quantity because it quanitifies the amount of lift being generated per the Kutta-Joukowski theorem."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "historic-exploration",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "$$\\Gamma = \\oint \\mathbf{V}\\cdot \\mathbf{dl}$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "governing-ivory",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The following video presents some interesting experimental demonstrations that you might find useful to further help clarify the two concepts. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "alien-bench",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEABALDBgXFRcXFxcdHRceHx0fHR0gHSUdHR0fLicxMC0nLS01PVBCNThLOS0tRWFFS1NWW1xbMkFlbWRYbVBZW1cBERISDAkMIw4OI1c2LDZXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV1dXV//AABEIAWgB4AMBIgACEQEDEQH/xAAaAAACAwEBAAAAAAAAAAAAAAAAAQIDBAUG/8QAOxAAAgEDAgQEAggEBgMBAAAAAAECAwQREiEFMUFREzJhcSKBFCMzQlJikbEVcqHBBjRDU4KSJKLhY//EABQBAQAAAAAAAAAAAAAAAAAAAAD/xAAUEQEAAAAAAAAAAAAAAAAAAAAA/9oADAMBAAIRAxEAPwDuaSSiTwSSArwPBZpDAFeA0lukNIFekNJbpDAFWkNJZgMAV6RaS3SLAFWkNJbgMAVaQ0lmAwBVpDSW4FgCvSGkswGkCrSGktwGkCrSLSW4DAFWkNJbgWAKtI9JZgMAV6RaS3AYAq0hpLMBgCrAaS3SGAKtIaS3AYAp0hpLdIaQKtItJdpFgCpRHpLFEekCnSLQXaQwBVpFpLsCwBVpHpLMBgCvSPSWYDAFWkNJbgNIFWkeks0hgCvSGks0jwBVpDBbgMAV4DBZpDAF2CSAaAAwSwGAI4DBLAwI4DBIAI4DBIAIYAngWAI4FgmLAEcBgkAEcBgkIBYFgkAEcBgYARwGCQgFgWCWBYAWAJCAQsEsBgCOAwSwJPICwGB4ACOAwSDAEcCJAwEIg7imuc4/qRld0l/qR/UCxDIU60JLVGSa7lgCDApTSko9Xlorr3MaeM5bfJJZAtwLBkfE49Kc3/xNNGrrjq0tejWGBLAYKrm6hSWZv2XNsjQvYVHhZT7SWMgX4AouL2nSeJt59myv+JU3yU3/AMWBsAy0r+EpacSTfLMWsk7m7p0ceJLGeWzYF+BGL+L0Pxv/AKs1U68JQ1qXw9+QFgYBbrK5DAWAK3cR8VU08yw2/QtAWAwSwAFo0AwAUqkVzaXzMlWtKcnCm8Y80u3ohrh8H5syfq2BsTzyAqpUVCOI8iqwrSqU25c8yWfmBqyLUu6/U5NCjrc/ErPaTWM42I1aUG9FFSc/xZeIgdkUppc2l8zPXrOlR1PeWyXqzPS4dr+KtJyk+mcJegHQTzyGU21sqSai3p7N5wRrSarU1nZpgaBSaSyxlF99lP2AtAVLeMfZFHD56oPPNSkgNADACutUUIuT6ELWuqsFNfNdn2Mt/cwVWnTnJKPmfr2IW9aEbhqEk6dTfbpMDoVFJxlpeJY2ObON3GLcqkEvY6hn4j9hU9v7gZVb3T/1o/oabSnVjnxJqXbCwX0nmEX6IkBRbzcnUz0lhFrkkm3slzZRaeasvz/2I3y1OlDpKW/sgK3WrVfskox/HLfPshKVWlh1JKcM4bxhx9TclhYRXcxTpzT5aWBVeXapKEukpJN9k+poTzujk3b1WlFy3+Kn+5fmVs996D5d4f8AwDfJ7P2ZnsPsovvll0mnBtPKcXh/Ip4a80IewGkWBjAjgCQsARDHclgAM7tKW78OP6Gaxt6cqUZOnHLz0N8llMy8NeaEPTKAjfRjCk8JRSafbqaVUjhfEuXcqv6EZ0pqSzhNr0Zmt+G0ZU4Sw94p+ZgTnWjK5pKMk/hlnDyTqzaqvEdT0cih20KVxR0RwmpJ/oacf+R/w/uBCnc1HJLwJL1fJFl1WcUlHeUniKNBkvU4uFVLOhvK64fUCVK1jH4pfFPq3/YL6lqpt/ej8SfbBZSrwmtUZpr3Kbqq5Lw6W8pbN9IruwIXVRStZTx93JG2vnVgnSpZxhZbSwWXdPRaziukSl2LSjVoPRPSsr7s9uoF0bepKSlVktt1FLZMpv6+itS+rc9pbJZLrW+U5eHNOFVc4vr7BXmo3FJtpLTLcBW1yqkseBKPrKKSL61NaJJpYw9ug3c01znH9RTmpU5OLysPdAYbOtKkoQq+SS+CfRejL6teVSTpUXy89TpH0XqWUqUalCMZrKcUWW9vGnFRjyAz0aEadfEV/p7t8289TYZp7XMfzQa/qasAAAGALUFR/C36MaHgDFwvDp56tvP6m45kZO3qSTT8OW6fZmifEaMV5032W7A1sw8L2VRdpyNdGpripYaz0fMos6EoOpnlKWV7Acp2UpTqVYbtTacekkdSxqwlHEVpkvNHGGmSs6Dg6mrrJtewrmwjN64txn0kv7gV8U2jTb5KcW/Y3IzU6U503Csk3yyuvqUL6RR+FRVSK5POGB0DHePTUoy6Zab9yrNxVnDMfDgnl75b9DZcUFUg4y5fsBYZr2ScfDb3nsl1KnRuVtGrFru1uStLGUZupUnrqPk+i9gNdOOFFdlgyWDw6y7VGbcGStw6E5OWqUW+eHjIGrIpSSTb5JZZjXC4pp655Tz5mbNOVhgYbS3jVi6lSKbm8rK5LoF5ZR8J+HFKcfijhdUboxSWFyADPaV1Vpxl16+jC+i3RqJLL0vYnSt4wcnHbU8v3LcAU2zzSp7Y+GP7Fo8BgDLQptVara+F4aY7uk5RUo+eL1Jd/Q0MAMC4pBbVIyhLs0yNW6lVzClB77OT2SR0HFPmgwBiu7Nyowpw+7KD+SZrccrD3RIAMEbWpSeKUk6T5xl932LbChKnT0S6N49jUACAYAIBgAgAAEQp0lBYisLmWABCUU00+TWGKnTUYqK5LkTQAZ69vrlTlnDg8+4V7SNR5bafdPBoADFCxcZJqtPC6N5RsGAGKrwujJ50tPrh4yaKNCNNYgsItACLQYJBgDPc2lOqsTjuuT6olKhGSSlFSx3LsCAz/QaP+3H9C2nSUViKwuxYACAYgE4ptNrdchgAAAwwBaMBgJxzzK1awznQs+xcACSHgYALADABAMAEIkIBYDAwAQDACIEhMBCJCAQDABAMQCwGBiAQIYAIBhgBCGwAQDEAAAAIYAACGAEUMIgAAAAAAAAAAAAAAAiQgAQwAEAAAgAYCAYAWjBDABgAAMQwAAABDwMQCAYgAAABAMAEIYAIQwAQDwIAEMAEJDFEBiGIAAU5KMXJvCW7ZRZ3kK8XKGcJ43QGgQwAQDwIAAAAQwAAABAJIYkMAAAARGc1FZk0l6vBMwxgqleo5LKglFJ8k+oGyMk1lPKHnp1KJ6aNObisJJvHqUWdnhqtJt1ZLL32x2wBuGV1XJRbgsy6JvGTBWnduLXhxjnqpbgdIDm8Jq3EnLxN6a2Taw8nSAAATmsqOd3yQDIqpHLjqWpdM7lVe18R58Scf5XgwVeEOVRfG9POUs/G/QDrARjHCS7bEkAAMALUMEMAGIYAMAAAGACAYmAgGIAEMQAAAAgGIAEMAEAAAgAAASGIAENlNzcRpQc5dOS7vsBwuMXlSTdNrSt8R6v1Opwah4dvBd9yu1svEzVrxzOe6T+6uh0YxwklyAYAMBCJCAAAQAAABXXk1Tm1zSbRz7ejcVIRn9IxqWcJHSqLMZL0Zl4U829P2wBTDxKVanGVRzjPVzXJo6JXUoRk4yfOLyiVaTUW4x1NclyyBIRzXxCsv9Bf90dCnLVFPuuXPAEjHZfaXC/P/Y2GO32uK676WBLiX2E/l+5ph5Y+yM/ElmhP2yX0/LH2QBVqRgtUnhbIlqMvEIZjDspxz7GXiVdSnRUMtRmpTa3UUBv+kLxPDaw8ZXZotyc2daNW6peG8qCk5Ppv0NF9Z+MorW44ae3Ngam8LL5GGjWhVuZOLyoQxn1bLL9PwuTkk1qS5tdjFZt0as3KnLFXDjpXl9GB1KtXRCUn0Rh4XKanVhU8zamvaRbfQdRwpJtJ/E32wZqdvUo3NOUqjmppxe3LHIDqgAAMQxAXjQhoBgAwBDEhgAAAAIYAAhiAQAACABgIQwAQAAAIYALAhsQAJDEApyUU23suZzKEXc1FVl9jF/BHu+47iTuanhR+yi/jl3fY6MIKKSSwlyQEgAYCAYgAQxAAAMCIwABGHhW1HD5KUl/U3mL6PFTjSXkeqbXdgaYSTWU8kiNOMUsRwkTwBjfC6Lbbhz9WaKVKMIqMVhIsEAGJbXbX4qef0ZtMVb/N0f5JAXXqzRqfyslayzSg/wAqC5X1VT+WX7FfD39RT/lA0kYwS2SS+RIAIxpxXJJfIkAAAAAAJoYAAAAAAABcMQwGMQwAYgAYAAAAAAAAAIQxAAAACAAAQAAAAAAhDYAI59/XlKSoUvM/PL8ES6+uvDSjHepLaK/uFla+Gm3vOW8n6gWW1vGlBRj8/VloDAQwABAAADEMAEMAAAAAAqrUFPGea5NbMtBgVUaSgml8/UmCGACGIAMsbeXjOrKWcLEVjkagAz3ds6sdOtx746rsW04KMVFcksEwAQwAAAAAAAAEAxAADAAEMQFwxEkADQhoBiGAARjBJtrrzJAAAAAAhiYAIAAAATeOYABiuOK0YPGrVLtHdlH064qfZUNK7zeAOmBzPo93PzVlFdooP4ZUfmuagHSA5q4ZNcrmf7jVjXXK5fzQHQKLy6jSjnnLpHq2cy9rXFFJeMpzeyio7iocOuNXiznFzfRrOkDdY20k3Vq71Zf+q7G05ztbib+KvpXaKwH0GsuVxL5rIHRA5+q6hzUai9NmTo8RhKWmScJ9pbZ+YG0QAAAAAAAACQwAAYCGAxAAEYjCIAADEAYDAhgIYCAAAAGIAAAAAAQwAQwAAEMALSQhgAxDAYCGAAAAAgYIAYmMQAJvCy3hGe8vqVBZqSS9OrOfLxbmPiOL8LnGmnhz9wNFTieqThQg6kusuUF8yuPD6tXe4qt/khtEu4fKfldDw4Jbbmi4rSglppuft0AVG0p01iMEv3LiqhVco5lBxfZ8yr+JUs4baecbpgagATeAAxXd1NS8OlTbm/vPyo0wnKT8uI/1ZYBhsrDw25zeuq+cn09EbQAAEMQAV17eFSOmcU1/VFoAc1VJ2zxUblRfKXWHozoQkpJNPKfJhOCknFrKfNHOoydtVVKT+ql9m/wv8IHSAAAAK69eFNZnJJGJ8RnPajSb/NL4YgdE5vEeJ+C9oqXz3D6JVqb1arS/DHZEp8IoyUVpxjn6+4Ge1vritNKMYqHV8zsEIU1FJRWEiYAAxAKIxIYAJjEAAAAAMAAQAAAAAAAAAAAAAIBgAAAFwCGAMYAAwBAAwI5ABgIMgMyX914ccR3qPaKLrisqcHN8l/UxWNBybr1PNLkvwoCmy4Rv4lw/Eqvf4t1E6sYpbJYQwAAAAAi4J80n8iRT42qWIp46voBC8TcdKUsvk08YK7ejWg1qqa4+q3RsACE5xisyaS9SSeeRGdNSWJLK7Mqt7VUm9Lel9G84AvApldRjLTLK7N8mXJgIBgwEAAAii8tlWpuD2fNPqmO5uoUlmckuy6sxupXr+T6ql3a+NgFvxOMaeKz+si9LS3cvXA1WuK3kiqUH96W8sexmnaRtq1Or5oy+Gbe7y+p2ljpyAw0OGwi9U81Jd5bm1RSGACQwAAABgAmMQCiMUQAYmMQAAAAAAwIgPAYAQAAAAAACGACAYgGIAAuGIYDGRGAxkQyAGK5vZa/Dox1TXmfSJtZh4ZjFR/e1yz35gW2tSq8qrFJrk1yZoAUnhNgc69fi1oUV5V8UzopYRz+Exz4lR85Sf6HRAAAAAjUqRisyeBOqtWnm/ToLwouWprL9QHTlqWcYJYKbq7hSWZPfpFbt+xlUritusUo/rIDoC1LujF/Dm/PWnL54GuF0+8/+zA2iMf0Scd6dWS9Jbpmqm3j4ufUAqQUlhpNepVQoKGdOcPo3yLxAMRnuaVTKnTlulvFvaSK63EIwgnJPW+UFvJvsBrk0llvC7nPq386j0W0c96j8q9u4lb1a7Uq/w0+lJP8Adm+FNRWIpJeiwBktuHKD1zeuo+bf9jbgBgZ76h4lGcO629yHDa2ujBvmvhfujUc7h3w1rmn+bUvmB0gAAEhgDAAEMAABAKIxRGAgwMQAAAADEMAEMAEIkIAYhiAAAAAQDAQAAFwxDABiABgIYAYK/Dczc4VJQb545M3iYGaxp6Yv43Pd5b7llxtTn/KzHb3UKScJ5UlJ9HvuYrviFapUdGlHCktsrfHUDpcJ/wAvA2HP4ImreKb3Taf6nQACFOopZxyXUTxNNZ264JJJLCAIwSzgjWqqnCU5ckTMHFfi8KHSU1n2AVjQc5OvV8z8i/DE6AkgAx1b5qcowpSljm1yJWl74kpwcHGUcZT9SHDnl1Zd5slRtpRr1qjeVPTj0wBrABAMQFHiylLEV8K5t9fYCm4u5a/CorNTrJ+WHqwt7FQbqTbqVebk/wCyNaWBVKihFyk8RXNgQt7iNRPHNbNPZplxgr1FplWoJTqJdHzQuE3069PVOKW+NuoHQAQABzV8N8/zQ/Y6Rzbva8oPupIDpIYhgAMAAqqVYw3lJJepW76l/uR/UtqUozWJJNepV9Co/wC3H9AM0ONUW8Syt8ZxsdHJz6HD6aqylzw9o9InQASGRiSABDIsBhkQAAxDAYAIBiAAAGAgAQwAQDEAAAAXDEMAAAABiGACGAEXFPmkR8OOrVhZ5Z64JgBy7afg3E6T2jN6oZ79Ua7mNSTUYvEX5pdTBxSXjVI0aazNPLn+Attr5wapXG01yl92QHQhBRWFyJCTyMAMXE4/DCovuSTft1NopxTTT5PZgEWmk1yYMwUav0eXhVH8H3JdMdmW2d9Gq5pc4vHuu4FdvBwzCGFGL3b5vJtRVK3g5a2viLQGJsZBSUltugI06rlnC+Ho+5PAJEa1WMIuUuQEZtPMM4bXfDx3MFfhiacfGnqfJOWf6Bb3Cq3UpJNKEEnnY59S8Trwr6+VTQo/k5ZA2aYWFCU0tUnhPpqZptqkdMakMaJrLS6Mwf4jpTnBSX2cFl+rIcPXhRUdTcJJPD6P0A78XlDMNO40NRb58jVGqmBYcy+/zdt7SOlk5t1ve0P5ZAdMBA2A8gcDiHF61Oo4x0ad8PmRsLi7rzUtaVNPfbZ+gHVrSbuKcE2kouT9R1rzTXp0lh6k2+6HXtlNqSk4ySxldir+GwwmpNVE86/vZAlVem6hj78XqXt1NjMlCzUZa5Sc58svovQ1AESRGIwGJgIBgAwEMQwAAEAyIxAAAAAACAYgGAgAALhiAAGIYANCAAAAADn8SvXDFKnvVly9PUt4heqhDPOT2iu7KOGWbWatXerL/wBV2AusLNUY45ye8n1bLq9vCotM45RaAHNjY1af2NT4fwy3RN38qf2tNr8y3RuYNJrDWUBVQuYVF8Mky0w1+Gxb10n4c/TkypXdej9rDVH8UewE7zhFOtJylKXtnZEbbhXhTjKFR4W2GuaNtvcwqx1QkmWAAARqR1LGQIzSqRxnb0JpJLCCMUlhLCGAGO/p1JaHTxlPO/I1gwOOuGVl4jVVfWefbf5F1bhNJ0XCMEpY2fr3OixAcnjGqNi1PGrEU/crrUtMItc0kX/4hf8A4z/miXVKeYr2QHNmvEjHfEk9idKvOlLTP9TNcvwpZzhEa3FqWjS1qf7AehoV01zJzpRlJTx8SzhnmbG+a3edJ3ra6UktwM3EOKyoy0Rpty6PozFGheXHnloi/kd/CfRP5GO6nc6sUowx+JsDIv8AD1PCzOTlnd9zp0oQpxUVhJfI5krS8l5q6ivyrkQ/hTxqqXEpR6vOwHWdxBc5x/Ug76iv9WP6mGHBqGE2216s0wsbaks6YJd3gCUuJ0V9/PsmyD4tS6an7RZdSVB50KDa7JFidPpp/oBz7ji/1ctEJp9G47GOHHqyaUqal7ZydupCFSEoSxvtzSYqFpTp+WK9+bALO58WGrRKD7MvAAHkBDQAMQAMQAAMQxAMQAACGIAABAMQxAXZGIYAMQAMBAAyuvXjTg5SeEiU5JJtvCXU4jze1sb+BB/9mBOxpO5q/SKnk+5F/udlIUIpLC2SGAAAAAAAAJoYAc244e4t1KD0T5tfdkW2d74mYyWmoucWbDk8W3nTjT+3zt7eoG+q5NqMdu8i2PIxWF46icJrFWOzXf1NoDAQwEAAAjPc3lOkszkvbqXVIaotZxlc1zRwFZ1qVTOiNWT5OT3Ao4zxKdWnhU3Gk5L4nzbQqvH5umowgtWMNnSqTrSWmtaxlH0a2KVbV5rTCnClT9ssDg6K1w22849cEaVjieJvCPR0rDwU922+bOXfLFRpLdLIG2zsqKW+X7tlVv4ii6sF9Xqe3VLJjldtU5Lq9l8ztcNahShD0QGm0u1JLc2J5OXc2TX1lL5xHZXuXh7MDpyjlNHn3PQp2ab86w/yPdnoIyTM8rCDrqv9/GPQCiVsqtXRL7Omksd2yi6wqMvh1eHUxFehvq2sZT1ZafJ4eMkqdtCMdONs5+YGXh9JSfjbbrCS6B9DpfSPL91S5vmaqdvGM5Tjtq5rpnuS0xc9eVqxgDl0aEZ1qmqnN/G/jUtkXXV27eWPNGXLvF+voaFbxVRzU3Ft7rOzLnbQbbcctrDyBC0pOKzKWqUt2+hoIU6ajFRXJciYAMQAMYhASECAAABADAGIAAAABDI5AkIYgLRiGAwAQDDIjk8RvJTl9Ho+d+Z9EgK7u4lc1fApbU155I6trbxpQUIrCX9Suws40YaVz6vq2agAAAAAAAAAAEAEZzUU29kt2BVdXCpQcnz6Lu+xnsLdrNWpvUnz/KuxVbxdxU8af2cfs1/c6QGG/otYrQ88f6o1W9ZVIKS6osa6HNtZeFXnRflfxQA6QxAAEZySWW9iNWtGCy+uy7sbSkt0BClV1rKTS6Z6k8EkJgAsDACmrHY5VzbLU5Y3Z1LxtU5aV8T2Rzb6ThCEFvVlsvT1A89Vj9bleWL27ZOzaV843K61mow0pcjHQk4SwwPS0K2diq8sFP46e0/6MzW9U6VKplAc61vHGWme0kdWnVUii7s41VnlNcmc+lWnRloqfJ9GB16ucfDz2M1W8cKuJr6t4+L8L9S+jWUkSlSjJttZysP1QGK4qtRrNb7xxhhw6jHLn4coS/My2HD4qM4JvTL15Erez8OWdcmuzArnRcXKWlTWctZ+JexqozUoqS5MqqWcJSb3WeeHsy+EFFJJYSAkIYAIYAAwwAAAAACAAABDEACGIBkNO+SQZAAAALMjIkgAYgA53FLyUcUqSbqS69l3J8NsFRjl7ze8mbNKznG5MAGIAGIAyAAAADEAAByrmbuKqpRf1UfO11fYs4jdtNUaf2kub/Cu5ps7ZUoKK+b7sC6MUkkuSGAABz+LU3ojUj5oPPyN4pR1JxfJrAEaNVThGS5NE3I53DJaJVKEvuvMfZmilCUpa57PdRXZAWeCtet7vp6FohgMQCABiM17ewoxy3mT8sVzkwHcXHh5cl8KXPq32Rlt6Tk3Vn5pcl+GPYdtQnUxUr8+cYdIm3SBz60Dl3NHfJ3a0DDXpAYbapjY6dvVOVKnhmmhUwB2oTyhV7eNWOJL59jNQqmyEgOO1Ut5Ye8ejOlbXSkluW1KanFqSyjk17aVCWqOXD9gO2BitLxSSNiYAhiQwAQxZAYCyAEgAQDEAAAAIBiAAAAABC6jFgBgAgJkiI0AwAADI8kRgMBAAwEADAWQyAzJxC9VGHeb2iu7J3l3GlByfPou7MNjaynPx63nflj0igLeG2bhmpU3qz3b7ehvBAAwEGQFKWBiMPErpwSp096k9ku3qBhvbl/SlOlFtUtqjXVHZpVIzipReU+pRY2qo01Hm+cn3ZXUsXGWuhLS+sXvFsDaM49TidenJRnbt+sd0Wfxf/8ACpn2A6gm8bt4Rz1cXE/LTjD1byL+HObzWqOf5VtFARueJylJ07eOuX4vuodlw1qXi15a6vTtH2N1KlGCSikl6FgAAABCccmSrA3GerEDl16ZnjsdGrAx1IAW0ahvo1DlReDXRqAdNMGs7MppTLkwOXc2TpvXS5dUW2V6pczoHPvOH5euntLqujA6MJJobORY3cs6Xs08HU15AU5bbFSbLUiSSAgppDU13M93PfCMfjNMDq6hazArkauAN6qIlkweMThVA2AUxqFikBIBZAAAAAAEwAYCABkiIwHkMiGAwyIAGAhgAALIDK61WMIuUnhIlOelOT5Lc5MVK7qapbUV5Y92A7SjKvU8aovhXkj/AHOrgEsJJAAwQCAeQEAFdzXjSg5y5IwcMpOcpXFTnLyp9EV1ZO7r6F9jB/F+ZnWSwsIBgAAGAwAAAAAAAhgACABkJIkDAyVYmSpA6FRGSpEDG0OLwWSiQaA1UahspzycuDNdKYG5MCnxFgqlcdgCvTjrzhZ7lsZGXLbLYAakwk9iMRtAZqkMs53EFpaZ12jncRp6oS9gOargkrk5Lr4bQvpIHaV0aaVwefp18tHRo1gO1Tq5L4zOXSrGmNYDdrDWY1WJqqBrU0PJnUh6gNAFcJ5JgMAAAJIgSTAYyORgMBAADEIBgIYAJLAxZAYCABgIADJzuJ3j2o096kv6Iu4hfKjHvN7Rj1bKOG2clmrV3qy3fouwGqytVRpqC+b7s0ZFkAHkBZABgIAGAshkBgwyAAAmxZAkIMiAUkZ6kTSVTQGGa3INGicSmSArwThIiIAu6+MFELn1K+IJ6VJdDmq4A9NZy1RbNcInN4ZP6tHRjIC1DIZDIBIorrKwWyZRUYHkuKWsozlJL4efscxyPbKlFyakspnF4vwRxbnTW3bsByKFTDN1OucrkyyFQDu0q5qhXRwIXGC5XYHcdcnTrNvY49Gs5vCOzb09KA1wkWpmaLLoAWReGjUnsZC+k9gLQEGQIkkwABhkAAExgABkAAAAAABAAAGQAAyZ767VGm54b6JLuAAYeH2spy8etvN7xX4UdQAAAGAALIwAWR5AAFkEwAAyGRgACGAAIAAMkWAAUziZpoAArDAABXWilF55HmrjaTxyyAAd/hdX6qJ0YVgAC1V0DroYAQdUjJgAGerPDTNsJKcQADg8Y4KnmdNb9UeanFxeGAAEZEk9wADucMpJLPXY6sWAAWxL4gAEhwniQABoyGQAD//Z\n",
      "text/html": [
       "\n",
       "        <iframe\n",
       "            width=\"800\"\n",
       "            height=\"400\"\n",
       "            src=\"https://www.youtube.com/embed/pOA3VJHCnWs\"\n",
       "            frameborder=\"0\"\n",
       "            allowfullscreen\n",
       "            \n",
       "        ></iframe>\n",
       "        "
      ],
      "text/plain": [
       "<IPython.lib.display.YouTubeVideo at 0x7f885c3d5ac0>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import YouTubeVideo\n",
    "YouTubeVideo('pOA3VJHCnWs', width=800, height=400)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "hidden-crystal",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "## Potential flow theory"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "artistic-austin",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "We are fixated quite a bit on the rotationality of the flow $\\nabla \\times \\mathbf{V}$, or its lack thereof. This is for a good reason. The irrotationality condition $\\nabla \\times \\mathbf{V} = 0$ guarantees the existence of a $\\phi$ such that $\\mathbf{V} = \\nabla \\phi$. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "153dfae8",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "$$\\mathbf{V} = u\\hat{i} + v\\hat{j}$$\n",
    "$$u = \\frac{\\partial \\phi }{\\partial x }, v = \\frac{\\partial \\phi }{\\partial y } $$\n",
    "$$\\Rightarrow \\nabla \\times \\mathbf{V} = 0$$\n",
    "$$\\nabla . \\mathbf{V} = \\frac{\\partial^2 \\phi }{\\partial x^2 } +\\frac{\\partial^2 \\phi }{\\partial y^2 } = \\nabla^2\\phi = 0$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "visible-combine",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "This, together with the continuity equation for an incompressible flow gives the Laplace's equation $ \\nabla^2 \\phi = 0$. This is valid irrespective of whether the flow is unsteady and is true for both 2D and 3D representations. Alternately, if we begin with the incompressibility assumption $\\nabla.\\mathbf{V} = 0$, then for the special case of 2D flow scenarios the irrotationality assumption allows the possibility of some function $\\psi$ such that it automatically satisfies the irrotationality condition.  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "miniature-slide",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$\\mathbf{V} = u\\hat{i} + v\\hat{j}$$\n",
    "$$u = \\frac{\\partial \\psi }{\\partial y }, v = -\\frac{\\partial \\psi }{\\partial x } $$\n",
    "$$\\Rightarrow \\nabla.\\mathbf{V} = 0$$\n",
    "$$\\nabla \\times \\mathbf{V} = \\frac{\\partial^2 \\psi }{\\partial x^2 } -\\frac{\\partial^2 \\psi }{\\partial y^2 } = \\nabla^2\\psi = 0$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "wound-party",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Physically, $\\psi$ denotes a streamline in a flow and is called the stream function."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "polyphonic-surfing",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "\\begin{equation*}\n",
    "\\begin{split}\n",
    "d \\psi & =\\frac{\\partial \\psi}{\\partial x} d x+\\frac{\\partial \\psi}{\\partial y} d y \\\\\n",
    "& =u d y-v d x\n",
    "\\end{split}\n",
    "\\end{equation*}\n",
    "$$\\text {streamline equation} \\rightarrow \\frac{d x}{u}=\\frac{d y}{v} \\text { or } u d y-v d x=0$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "rural-genius",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Not just that, the difference in the magnitude of $\\psi$ between two streamlines is related to the rate of mass flow of fluid between the two streamlines (for $\\rho$=const). "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "prescribed-danger",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$ \\rho V = \\rho\\Delta \\psi = c_2 - c_1$$\n",
    "<center>\n",
    "    <img src=\"../assets/B37_fig2_40.png\" class=\"large\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Flow between two streamlines. <a href=\"https://www.mheducation.com/highered/product/fundamentals-aerodynamics-anderson/M9781259129919.html\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "portuguese-hierarchy",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Note that $\\psi$ and $\\phi$ are not equivalent here. $\\psi$ can only be used for 2D flow representations and its form is coordinate dependent while $\\phi$ can be used for modelling 3D flows and is coordinate-free. Of course all of this is contingent on the continuity equation, which is a fundamental law of nature, and the irrotationality assumption. Its the latter condition that one needs to be mindfull of, since not all fluid flows satisfy that condition and potential flow theory based solutions would not be accurate in such scenarios. Fortunately, the assumption applies for attached flow around an airfoil (the experiment of Bryant and Williams {cite}`Add54` discussed below) and airfoils are the fundamental building blocks of most aerospace applications potential flow theory has developed into a viable theory for analysing flows."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "valuable-accessory",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Since $\\psi$ and $\\phi$ are different representation of the the velocity field itself, from the above discussion it follows that - "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "golden-climb",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "$$u = \\frac{\\partial \\psi }{\\partial y } = \\frac{\\partial \\phi }{\\partial x }, v = -\\frac{\\partial \\psi }{\\partial x } = \\frac{\\partial \\phi }{\\partial y } $$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "possible-leonard",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The above equations present a pair of <a href=\"https://en.wikipedia.org/wiki/Cauchy%E2%80%93Riemann_equations\"> Cauchy-Reimann Equations </a> such that a complex function given by F(z = x+iy) is complex differentiable (*infinitely* differentiable) if it is given as -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "lesbian-bachelor",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$F(z) = \\phi + i\\psi $$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "selective-number",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The derivative is given as"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "numeric-parker",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$F'(z) = \\frac{dF}{dz} =  \\frac{1}{2}(\\phi_x + \\psi_y) - \\frac{i}{2}(\\phi_y + \\psi_x) = \\frac{1}{2}(u + u) -\\frac{i}{2}(v + v) = u -i v = \\bar{w}$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "artificial-guyana",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "where $\\bar{w}$ is the conjugate of the velocity in the complex domain."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "impressive-layout",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "### Elementary flows"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "solved-population",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "So posing the problem of a *real-world* fluid flow in abstract terms of complex analysis provides us the insight that the functions $\\psi$ and $\\phi$, which are used to represent *real* fluid flow velocities, are indeed continuous and infinitely differentiable everywhere in the domain. Additionally, $F(z)$ also satisfies the Laplace's equation since both $\\psi$ and $\\phi$ satisfy it.  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "integral-translation",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$\\nabla^2 F(z) = 0$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "outstanding-average",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "This is a useful because the Laplace's equation is linear. If the function $F(z)$ can simply be evaluated for simple elementary flows and these elementary flows can be used to create the flow field around any complicated geometry then $F(z)$ can be obtained essentially for any complicated geometry. Knowing $F(z)$, the flow velocity around any geometry can be evaluated and thereafter using the Bernoulli's equation to get the pressure distribution is a straightforward excercise. The question arises, what sort elementary flows are useful?  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "designed-quebec",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "- Uniform flow\n",
    "$$F(z) = U_\\infty z $$\n",
    "- Uniform flow at an angle of attack $\\alpha$ \n",
    "$$F(z) = U_\\infty e^{-i\\alpha} z $$\n",
    "- Source/Sink\n",
    "$$F(z) = \\frac{m}{2\\pi}\\text{ln} z$$\n",
    "- Line vortex\n",
    "$$F(z) = \\frac{\\Gamma}{2\\pi i}\\text{ln} z$$\n",
    "- Doublet\n",
    "$$F(z) = \\frac{\\mu}{2\\pi z}$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "rural-extension",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The elementary flow fields represented by the complex function $F(z)$ satisfy the continuity equation as well as the irrotationality condition everywhere in the domain except at the point where the singularity exists. That is okay because the flow filed right at the point of singularity, or close to it, is not of interest to us. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "third-saying",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "### Practical application from elementary solutions"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bridal-harvey",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "- Source+ image Source pair\n",
    "$$F(z) = \\frac{m}{2\\pi}\\text{ln} \\frac{z+a}{z-a}$$\n",
    "<center>\n",
    "    <img src=\"../assets/BA1_fig3_12b.png\" class=\"large\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Solid wall modelled using the method of images. <a href=\"https://shop.elsevier.com/books/aerodynamics-for-engineering-students/merken/978-0-08-096632-8\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "descending-charity",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "- Source+Sink pair in horizontal flow\n",
    "$$F(z) = \\frac{-m}{2\\pi}\\text{ln} (z-a) + \\frac{m}{2\\pi}\\text{ln} (z+a) + + U_\\infty z$$\n",
    "<center>\n",
    "    <img src=\"../assets/BA9_fig8_13.png\" class=\"large\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Flow around a Rankine body. <a href=\"https://www.mheducation.com/highered/product/fluid-mechanics-white-xue/M9781260258318.html\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "hundred-whale",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "- Doublet in horizontal flow\n",
    "$$F(z) = \\frac{\\mu}{2\\pi z} + U_\\infty z$$\n",
    "<center>\n",
    "    <img src=\"../assets/B24_fig3_12.png\" class=\"large\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Flow around a circular cylinder. <a href=\"https://www.cambridge.org/core/books/lowspeed-aerodynamics/077FAF851C4582F1B7593809752C44AE\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "traditional-wallpaper",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "It is already shown above that $F'(z)$ gives the complex conjugate of the complex velocity. So, once $F(z)$ for any arbitrary geometry is constructed the complex derivative would give the velocity distribution around the body."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "acting-liverpool",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "\\begin{equation*} \n",
    "\\begin{split} \n",
    "F^{\\prime}&=U_{\\infty}-\\frac{\\mu}{2 \\pi z^2} \\\\\n",
    "&=U_{\\infty}-\\frac{\\mu}{2 \\pi r^2} e^{-2 i \\theta};\\quad \\text{where  } z=re^{i \\theta} \\\\\n",
    "&=e^{-i \\theta}\\left[U_{\\infty} e^{i \\theta}-\\frac{\\mu}{2 \\pi r^2} e^{-i \\theta}\\right] \\\\\n",
    "&=\\left[\\left(U_{\\infty}-\\frac{\\mu}{2 \\pi r^2}\\right) \\cos \\theta\\right] - i \\left[\\left(U_{\\infty}+\\frac{\\mu}{2 \\pi r^2}\\right) \\sin \\theta\\right] \\\\\n",
    "&=V_r-iV_\\theta\n",
    "\\end{split}\n",
    "\\end{equation*}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "moderate-regular",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "where $V_r$ and $V_\\theta$ represent the radial and the tangential components of the flow velocity around the cylinder. In order to confirm that a doublet in a free stream indeed represents flow around a cyclinder, simply the no-penetration boundary condition should be satisfied, i.e. $V_r(R)=0$. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "military-experience",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "\\begin{align} \n",
    "V_r(R)&=0 \\\\\n",
    "\\Rightarrow R &= \\sqrt{\\frac{\\mu }{2\\pi U_\\infty}} = \\text{const.}\n",
    "\\end{align}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ranking-blake",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The radius of the cylinder depends on the strength of the doublet as well as the free stream velocity combination used to construct the flow field. Usually, you'd want to model flow around a specific cylinder in a flow so $R$ and $U_\\infty$ are known and $\\mu$ can be obtained by satisfying the no-penetration condition above. This gives-  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "comparable-arthur",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$V_\\theta(R) = 2 U_\\infty \\sin \\theta$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "loaded-mortgage",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Applying the Bernoulli's equation, the pressure distribution around the cylinder can be obtained -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "equipped-allowance",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$p = \\frac{1}{2}\\rho U_\\infty^2 \\left ( 1-4\\sin^2 \\theta \\right )$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "deluxe-pantyhose",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "This pressure distribution results in a zero net lift and zero drag. The zero net lift and zero drag come from perfect fore-aft and top-down symmetry of the flow and the fact that viscosity is ignored in potential flow theory. Additionally, the stream function is given by the imaginary part of $F(z)$."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "hearing-decision",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "$$\\psi=\\operatorname{Im}(F)=\\left(U_\\infty-\\frac{\\mu}{2\\pi r}\\right)\\sin \\theta$$\n",
    "$$\\Rightarrow \\psi |_A = \\psi |_B = 0$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "collect-psychology",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The streamline $\\psi=0$ defines the boundary of the cylinder."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "random-greene",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "One of the early experiments by Fuhrmann and Prandtl involving flow around a slender airship-like body in a wind tunnel showed that indeed the measured pressure distribution around the body matched that predicted by theory {cite}`BA13`. The theoretical pressure was obtained using potential flow theory where the geometry was modelled by placing multiple sinks and sources in a uniform flow as shown below."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "solid-campbell",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA13_Fig53.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Fuhrmann and Prandtl study of a uniform flow around an airship-like solid body. <a href=\"https://press.uchicago.edu/ucp/books/book/chicago/E/bo12024403.html\">[source]) </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "controlled-backing",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Due to the low magnitude of kinematic viscosity ($\\nu/\\rho$) of air, most flows of interest to an aerospace engineer occur at high Reynolds numbers (>$10^6$) and low angles of attack $\\alpha$. In such conditions the flow can be modelled reasonably well by dropping the viscous terms (i.e. assuming $\\nu=0$, limiting case of a viscous fluid).  Of course, one runs the risk of not being able to model effects of flow separation (which are viscous in nature), viscous drag etc. But the speed and ease of solving the fluid flow equations assuming $\\nu=0$ as well as the physicality of the results has led to the popularity of the potential flow theory. Even extensions to this theory have been proposed that account for the viscous effects in some adhoc way but without having to solve the full NSE equations. In fact for low Mach number flows, the fluid equations can be further simplified and yet be physically representative of the real flow field.  \n",
    "\n",
    "When the irrotationality and the incompressibilty assumptions are made along with zero viscosity, the NSE reduce to two eqs - the Laplacian and the Bernoulli's eq."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "precise-modem",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "$$ \\nabla^2 \\phi = 0$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "excited-break",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "$$\\frac{\\rho V^2}{2} + p = const$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "supreme-patrick",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Here $\\phi$ is some unknown function that needs to be found out and all the quantities of interest in the entire flowfield would be known. $\\rho$ is a constant quantity that is known apriori, if the Laplacian can be solved i.e. $\\phi$ is obtained over the entire fluid domain then the fluid velocity can be obtained using the following relation -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "physical-chorus",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "$$\\mathbf{V} = \\nabla \\phi$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "alpine-kidney",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "This is highly convenient because $\\phi$ is a scalar function, only one known quantity at a given coordinate location, but velocity has three components at each point that need to be found out. Additionally, $\\phi$ is independent of the coordinate system and, based on the form of the $\\nabla$ operator, the velocity vector $\\mathbf{V}$ can be evaluated in any coordinate system.[^nablecoordinates]  \n",
    "\n",
    "[^nablecoordinates]: For most applications the gradient of the potential function $\\phi$ is evaluated in the cartesian coordinates. However, finite-state inflow models exist (also called dynamic inflow models) where the induced flow of a rotor is evaluated using a potential flow theory-based solution of the rotor. Here, the Laplacian is solved in ellipsoidal coordinates because the problem is amenable to an analytical solution in such a system. One of the most common such inflow model used extensively in real-time rotor simulations is the Peter-He inflow model. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "brave-enlargement",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "- the early go-to model for obtaining pressure distribution of slender bodies\n",
    "    - airfoils, airships etc.\n",
    "- bodies of any shape could be analysed\n",
    "- conformal theory used to obtain analytic results of pressure distribution  \n",
    "- also used in unsteady aerodynamic theories - for e.g. Theodorsen theory"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "configured-twist",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The potential flow theory discussion above is partially based on the details presented in {cite}`BA5`(Chapter 3). Students interested in the derivation of some of the results presented above can refer this source."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "foreign-craft",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Kutta-Joukowski theorem"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "devoted-posting",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Addition of a line vortex doesn't disturb the cylindrical geometry established by the doublet+uniform flow combination. However, the flow around the cylinder is changed dramatically so much so that now when the integrated effect of the pressure distribution is evaluated a non-zero lift exists. The drag force is still zero. Intuitively this makes sense because the flow is no longer has a top-down symmetry. Instead, a clockwise circulation leads to higher flow velocity at the upper surface of the cylinder and flow slows down at the lower surface. We know from the Bernoulli's relation that this should lead to lower pressure at the top and higher pressure at the bottom. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "saved-cabin",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "$$F(z) = \\frac{\\Gamma}{2\\pi i}\\text{ln} z + \\frac{\\mu}{2\\pi z} + U_\\infty z$$\n",
    "<center>\n",
    "    <img src=\"../assets/B37_fig3_32.png\" class=\"large\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Flow around a circular cylinder with circulation. <a href=\"https://www.mheducation.com/highered/product/fundamentals-aerodynamics-anderson/M9781259129919.html\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "awful-commissioner",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Clearly, if indeed there is lift generated in the above case it has to be related to the circulation strength $\\Gamma$. In fact, if carry out the same excercise as the non-lifting cylinder case discussed earlier, the amount of lift generated (per unit length since its 2D) turns out to be -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "extreme-celebrity",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    " $$L´ = \\rho V_{\\infty}\\Gamma$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "moving-pendant",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "This result leads to the more general case that the lifting force generated by any entity in an incompressible, inviscid, irrotational flow in an unbounded fluid is proportional to the circulation around it and perpendicular to the oncoming flow velocity direction. This relation is a consequence of the Kutta-Joukowski theorem. A more general notation for identifying the direction of resultant force using the Kutta-Joukowski theorem can be obtained by defining the direction of $\\Gamma$ as a vector quantity based on the right hand rule and using the following relation to give the resultant lift force per unit span ($\\mathbf{Q}_{\\infty}=\\mathbf{V}_{\\infty}$) -"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "medical-customs",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    " $$ \\mathbf{L´}= \\rho \\mathbf{Q}_{\\infty} \\times  \\mathbf{\\Gamma}$$\n",
    "<center>\n",
    "    <img src=\"../assets/B24_Fig316.png\" class=\"large\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Vector notation for Kutta-Joukowski theorem (<a href=\"https://www.cambridge.org/core/books/lowspeed-aerodynamics/077FAF851C4582F1B7593809752C44AE\">[source]) </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "corresponding-stuart",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "This relation tells us the simple fact that lift is related to the circulation around the body in question. So if the circulation were zero, there would be no lift generated. On the other hand, if you knew precisely the circulation around the cross-section of cylinder, or more practical entities of interest such as a wing or a blade, you could get the lift per unit span and the total lift would simply be an integral over all such sections over its span. Precise experiments carried out in the wind tunnel by Bryant and Williams {cite}`Add54` with wing sections, where the flow velocities around these section have been measured, show that indeed the effect of an airfoil is to introduce a circulation around itself and that circulation everywhere else in the fluid domain is zero, or rather close to zero, such that the potential flow theory assumption of irrotational flow is not misplaced. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cathedral-cowboy",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/Add54_Fig8.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Schematic of the experimentally measured circulation around an airfoil section <a href=\"https://royalsocietypublishing.org/doi/epdf/10.1098/rsta.1926.0005\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "impaired-graphics",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Note that the contours used in the experiments of Bryant and Williams were perpendicular to the wake of the airfoil. This is to ensure that that the effect of the vorticity in the wake of the airfoil is not included in the circulation measurement around the airfoil. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "incoming-raise",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Kutta condition"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "defensive-celebration",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Based on the potential flow theory, lift generation can be predicted using a simple relation given by the Kutta-Joukowski theorem. This raises the natural question, how does one obtain $\\Gamma$!? Potential flow theory, on its own cannot uniquely predict the amount of lift an airfoil is capable of generating. Since the amount of circulation in the flow is undecided, some more physics needs to be modelled in order to account for this. We could either go to the origins of the potential flow theory and see if some of the assumptions made can be relaxed to include more physics for this purpose. This is bound to make the equations more involved that the Laplace equation. Alternately, this is done using the Kutta condition. To illustrate the situation schematically, consider the two flow representations shown below for the same $\\alpha$ and $V_{\\infty}$ but different circulations around the airfoils. Both are valid/possible within the framework of potential flow theory, which one is correct?"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "invisible-affiliation",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "<div class=\"image-grid2\">\n",
    "    <img  src=\"../assets/BA3_fig7_1_3b.png\" class=\"medium\"/>\n",
    "    <img  src=\"../assets/BA3_fig7_1_3c.png\" class=\"medium\"/>\n",
    "</div>\n",
    "\n",
    "<div class=\"image-grid2\">\n",
    "  <p > Streamlines around an airfoil (real flow) <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\">[source] </a>\n",
    "  </p> \n",
    "  <p >  Streamlines representation around an airfoil based on a mathematically valid potential flow solution (not a real flow) <a href=\"https://www.wiley.com/en-us/Understanding+Aerodynamics%3A+Arguing+from+the+Real+Physics-p-9781119967514\">[source]) </a>\n",
    "  </p>       \n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "tutorial-daughter",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "In reality, only the one on the left is observed. This is a consequence of the viscosity of fluids in nature and the fact that they cannot turn around sharp corners (airfoil trailing edges are sharp enough). Since potential flow theory doesn't model viscosity, a condition needs to be mathematically enforced so that the resulting flow field obtained is physically correct and unique[^uniqueness]. This is achieved by enforcing what is called the Kutta condition. Based on what your favorite reference textbook for fluid dynamics is, this could also be referred to as the Joukowski condition {cite}`BDI2`. \n",
    "\n",
    "[^uniqueness]: Recall that we carried out a similar excercise when we ignored the negative root of the quadratic induced flow velocity equation."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "employed-soccer",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<div class=\"image-grid2\">\n",
    "    <img  src=\"../assets/B37_fig4_18a.png\" class=\"medium\"/>\n",
    "    <img  src=\"../assets/B37_fig4_18c.png\" class=\"medium\"/>\n",
    "</div>\n",
    "\n",
    "<div class=\"image-grid2\">\n",
    "  <p > Instantaneous flow field just from impulsive start of airfoil <a href=\"https://www.mheducation.com/highered/product/fundamentals-aerodynamics-anderson/M9781259129919.html\">[source] </a>\n",
    "  </p> \n",
    "  <p >  Steady state flow field around an airfoil <a href=\"https://www.mheducation.com/highered/product/fundamentals-aerodynamics-anderson/M9781259129919.html\">[source]) </a>\n",
    "  </p>       \n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0cd0341e",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Potential flow theory alone cannot determine the amount of circulation around a body. This is true in the case of steady as well as unsteady aerodynamics. Even the widely used thin airfoil theory results of the lift curve slope=$2\\pi$, $C_{m_{c/4}}=0$ etc. are obtained after imposing an additional condition based on the flow physics. Purely following the mathematical construct of potential flow theory would not give us a unique solution (and hence *no* solution). The amount of circulation decides where the stagnation points would be. The Kutta condition tells us where the stagnation points ought to be in the case of flow around an airfoil. Therefore, by extension it helps identify a unique circulation that exists around an airfoil at a given angle of attack."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "similar-cotton",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The Kutta condition is an engineering approximation that is based on experimental observations rather than mathematically derived from first principles. Images of flow field development taken at the instant flow is started around an airfoil suggest that there some flow turning occurs at the airfoil trailing edge but large shear forces encountered quickly make the flow pass smoothly at the trailing edge rather than turn around it. The Kutta condition is an indispensable part of introductory aerodynamics textbooks owing to the success of potential flow theory together with this empirically observed assumption."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "studied-delta",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Kelvin's circulation theorem"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "vanilla-savings",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "For an inviscid, incompressible fluid of constant density, the circulation around any closed loop within the fluid remains time-invariant. For an airfoil initially at rest within a fluid, the circulation around any contour within the fluid is going to be zero. When the fluid is set in motion there is vorticity due to the uniform fluid velocity, yet we now know (from the discussion aboce) that there would be a circulation set around the airfoil. The Kelvin's circulation theorem implies that a vorticity of opposite sign must be getting introduced in the fluid so that the net circulation around any closed contour, that includes the airfoil in it, must be zero. The video below demonstrates this perfectly. A series of stillshots from the video above are also available in {cite}`BA14`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "1e90d731",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEABALDA4MChAODQ4SERATGCgaGBYWGDEjJR0oOjM9PDkzODdASFxOQERXRTc4UG1RV19iZ2hnPk1xeXBkeFxlZ2MBERISGBUYLxoaL2NCOEJjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY//AABEIAWgB4AMBIgACEQEDEQH/xAAbAAACAwEBAQAAAAAAAAAAAAAEBQIDBgEAB//EAD0QAAIBAwMDAwMCBAUDBAEFAAECAwAEEQUSIRMxQQYiURQyYXGBFSNCkTNSocHwJLHxBzRi0eEWJUNygv/EABgBAQEBAQEAAAAAAAAAAAAAAAABAgME/8QAIREBAQEAAwEBAQEBAAMAAAAAAAERAiExQRJRYSIDMkL/2gAMAwEAAhEDEQA/AMZGPa1W26RMyb2x7uapjY9vmrlXAD7CV8nHFNDDV47VZIDA24Y5oO9ljlC7VxiuRMrMFbBUnk1C5EfVIhPtpqOxyiOPAFVMARkV1c4xXDmg4K8MYNcNWxpmMk96KrHapdJwgYjg1wLzgdzR8yy/RqrDgecdqiF5OeKtitzLE7hh7R2qMkZTBxwfNRRnQNtYjNUewa8VGM1IMcVwg/NBzvxXAOK8Ac0TdWptQoY5JoQOgXPuNEXwh6i9BsjHNDqCzBf8xxVwtnEhTIyO9RVGwgZrmKuYjGzyKqINB4rXKkOK4aDx5r1eHNdwBVR0DNeOKshgaTnGF+a5JFtPByPmioeKjUu1e4qiOK7XTXKzR7NWJIVXFVYr371YJHBqWwdPcDUBU027Gz3oIivRsFfJr2a5QSZwzZAr2agK6O/NB0HBqwDeOTjFRcJkbePmrBGhA91UUnHiu8bfzXdqg4zXCAG+ag62dozXO9dZ93FRwRyKDo4rzAiuZNTZi3epRWKIhidvdjgVSq5NEJcSLCUVsD9BUg9KwcfpVS5auZ780fYiAK3VUH966QLtxru6pyctkDionA8c0E94KEYrkUYfPuxioZ4wKvOFg/JrIi1pKqb8e2qGyKJjuZAoQtlR4qE5DPx2qCCPyM9qIV0TJHc0LiuN2BpQU7M6OQOMUB85ppFIh09lUAPjuaVnnOasEowd+1eTWk05Vg06WG9ULG6llPekFrC0ivKh+znvRH17uFWc7kXsKxRfYW8EguJGyAPtH96WNx2+aJkn25MRwH7gUMfzVR4EivZ5rsYG8b/t80wltIp0eS1ICL3zQLjzVm1405U4Pmr7OKGR9j/d80druyF4UgwU29xTVLV2pHuxlq411I6bGb2nvXOjOY2fpttHnHFcgt5ZlYpGzBe5A7VpBd9dQywRxQx7dnnPegt2FIxUcAHg5r2Mnyag5mpDOM1ZFETk4yP1qcMSyuQTtUCqLNN0641OVkt1B2DLE+KonmeVtzjtWl9NanY6ZZXgdttxJwBjuMfNJkt/5kMM3tSU8uagHW3QW3WZvdngVAnaAyvlj3plr1lDYPFDBKJV25yKUsMHjnNCduZOc+a8ST3rq4zg14j4NFRIrldFeaqPA13II5FRrtVFq3DrGYwfaa51ONviq8V3OPFZVfIItg2kE1R+lcHJrp9vbmg8QRXq9vJXBrmTQSJz4rmBXq5xSI72rma7j4qNFSroUtwBUQcdqsikMbFh8YoICvEV5jnJrmeOaCYcBcbefmvA1DNd7d6okBnsMmvdjzxXYZem+4c12V+q5bGKg5kfFeB3D4qIr36VQUIdkJ3LnPINC5waLimyNsh4xxQzgBqg6kgHirbW3e4ZtniqAAT3oywuRb9QE4B81ZgGnjMTlW70Zplm151D1Qm0f3qm7G5g4O4nmqEkZC20sOPBpBIPjIIziq2OSatX2kgjmucM321UVDvVzFpFBCnA81FkJPtFEWl10FMRQNn5FMUMP81eXvu70ZPHHJEWjx35oHJUYpReCmzOOaqYgntgVyPlwPmrJVVScVBPfGLcgLyR80DRCoWDZ7AUMRgmrguty+WVTwRzUR3Kn5wa0Poj6UalL9ZGrI0ZAZhkA5FD+qLOCz1R3syvSY5AXxxXPe8Al9aQQQwmFyxYc58UJEnVlCfNGS6a8GnJdu6kP2FC7SsQMikZ7VpEZQFYqPFMNGilkaZwu6NU5FaTTPSNreaMkyzHrSLnPgH4oGMv6fln01QplnH3Z7U9Z3+ANGiQ3wnuI9sGfd+BQeqBX1BxAxZM+2qbppYpGgd8sDzTDR79I51SeLqqOcGk49tKzq0yaabB41B7Bsc0PHeXNpEUQhVYc8d6L9Q31tfXSm2h6W3g8YpYwZ5QpbJNUw6k9PNBo0WpF94bGUXxkUsguYobkOYwVHipxahc28LWxfKfFEWD2D28ou1w4B2mpIF8k++Z2UbVc9qLubRUtUuISSuAG/WrdL0+O+Ztxwc+0VRedWy6lkRwTVULBBJdS7I15p1dXdtEI4c9TppjP5pZY6lLp6SLH3fiu3um3FtbRXcigRz8qQaiYEkkeQksxPPk1wZfAHejI7VUshdEhgeNvmhC219y8UVKON3l2Ac1YkadOTecMKut7d9v1khGwc4zyaElcNKzAcVBAkV4H5qS7c5NRYc8VodGAcmvMRnioivHigka8hCvlhkV1TxxUSMeagk3uJK1HPzU4ZAjZIqJOWyKDleAzXia9kAcVR3OO9cODzRMAgeN+qdrD7aqlgaJVZwCrdqgiDxzXDg1Z75UAROFqoD3YPFB1cbsVeHiVSNuTVBG098iuNQSIySa6qZH6VFACDk1wEgjzQd7HtXWkLHBqczlyMrtrkkQVQd2c1RWR8V0CpGUmMLjtXMMUyF4Hmoq2B403dVc57VxipB6a1AkFQMVwcUR5eWAJwKLvDblE6P3DvVMcLOMjHNRwMgeBQRA5rykDOaOgtkvZVRR0xjljQ0ts0TMmQcGiitFW2mvkF6+Ih3/ALVy9FqmpOLY7oAe9VWzxxxOHQ7/AAa5D0em/UHv8VcQTKYFugynKAYoaRN8pMfY9qpUZHbmrYJOg25hmrpjiQysTwRipW42zHqLmrJbosvtOKv0/UUhUxTR5DcFqKDG55Gij8muS2sgbaRzVt30ob3fb8p3FeN/Ox3HseKVFVnHm5AbxXbj/wBw2O1cuImjCvn7+TUYSc471B1klWMueFNB5yaY3F2zw9EjhaXd81QVaXLRwPGON39Q70RcxAWay9VnY/uBS9AckCmEN8ItPktWjDB/JPasYND6d0kapaRrfQSLBGMhiSA1J/UrQrqDW1qymCLgY5/1rS6Hqc91oS21uE63CLmstr+kT6XdATlWaTkEGqz/APSGm67qWmjbbXDLH/lIDD/Wp2+oG612G71Jy4DZY4x/2oMWrrbiVwdh7EVKxtvq76K13bQ7Y3YrUOjr1Jc6Vf3qNZ4Df1sARml1rssbgT9QN+O9GeodMtNKkjiiYu5XLH4pH3OOc1asWX0qz3Tyou0McgCitIjtZJHN2+3A9vOKpubC4tokllXCOPaaqgYA4IAP5NZorZcynBzz3op4Q1uCXGR4qgbQ4J5we1WKplmwvCk9qij9EllivI7hvdFEeR2px6i0szWja3Gw2k52jwCcUo1HSp9NskkWYPHL3AGMf8zUl11v4B/DmQnn7t3jOe2Kz73AlyZGyfNF3Go3NxBHBLKWiiGEXA4qiNVwSTj4rqQSO6IF5fsTWhyMs4KbsL3xniosuB+Km0LxzFGHK/diuzmPf/Kzj4NB23WScrArHk+0ZrlzbSWlx0pgVPeuRmS3ZJ04IORVl9dyX0vVlxuxiiBjXsmvba9iqqUhUt7BgUZbaZLcWpuNwEY70D+lEwSbYGUyso+B5qwQuIVhk2pJuGO9UGpsdx45FRCnOMftUo8q7jxXiCDipbf25ou8t4EWI28vUYj3DHaoBGTaBznNRyMV7JGQa4KD3epmR2AVmJC9qiBUl2nOTiqJQzSRZEZxnvxUOXbnivZ28jvXOTUHWGPNdHuHbtXPFdIK/vVEfNTT2MD3qHaug5GRQFXRclTJFjI9vFCtxVjTyOoDHOOBUFQuSBQcKnaD80ZFcqlg0GwEt5quMCOUGT3qPAqqTb1CRwuaKlBH1ZdvYea7Mq9QrH48VFtoYGNsZqSGSJiRgn5oKyzr7SSMeKiCTxnipM5kYs3epdLC7twNAXaSk4hkk2xn+rt/rVN2u2YhJN6jsaP0XThqLGNWyw/px4oe5jjtruSIqfZxitZ0B94KgYx+agEY+/B2/NTjX6glVOCewoiY/T2nQcfzDzUAJJ8Gikt98ed2MDNUQRGaVUB70Zc2zpKsS+4nAGKhgUBAjdQ8jtVcY3MB4NEXtjPaqrTLgN2NWQ6czWZuNwAAqiT6bcRxbgoZTzQrOelsYcjtRUF1MsJV2O0eaY2tlZ/StNISz44pggdBu2sY7meQCMrkZpRE628/HIBp1LJPc6ao+owi8BaSzwhZCN4P5xUK5dMksrMgABoTtTOx0ue+gnliICxDJzSs8Eg+Dimoeen9EuNRie4gYHa20rkZxTPVdDtrW+t7bjfIAMDyTVXobVY7C6nimGI3Gd3xVuu63AfU1td2xEiRFQc9vFZ71Poezsb7QdatpZ4JOiWO0A53fsK76w1FNQ1ABVdRGMAMpFaDVvV1iIIHtlWaTdnBH20Lot7b65rkrS2q7dhxWp4zbl1mNK1FIGMd0u+E9gRnFPNV0KytbCO7s7jEzMCvu/Fd0vRrGT1LfWlwh2rgxjPbNCa7pC6fqkcH1UjQMNyhj2qU3arsbazuGUX9yesXAckk+3Pz+lGep9M0+xSOTTmjHl1BzWeukDXey2YkkhRVd3Fc279K43q/wT3p21D/AFpnv9KsWEi7QFGAaE0b03Pqs0qLIqpH3c0vmSWKKPLNzghQe1fSdOWGD04JIvYDGWLeaW1uR8weNra6MRAZkP7Gm2nPLfXAtbO1UzsOTgCkskjGUsCSc9zTb0/Pc2F0dQSMvGoIamMh7/6xJJLW53/y+SucgVfPpkCaFBfJNmRzhk54qi41Bri8uJ3H+L/pVdxqDTWUVqUASPsRTAGeBVxnkeNUZvs7Y8VZa232TyqegGwxphr9hbW8wms2JgftRS2Ni3LS7WJxz5onU9Lm09IpXIZJRkEGl4Uk4FEXF1cyQpDMxKp2BqoqjzI+wk4+KiABJg9q8jFZMiuM2TRHWI3cdq4RXO9G2Mse7ZOAV7Cigu1TwGHJxRF/a9BwV+1u34qsQYg6pcY+KorVijZHipNKTJvGM1xY2kGUBxmuEYOKgi2d2T5rynblgcGvD3Gmmk2VnOxe7n2KvjI5pArUNIwUDJJ4qU0LwPslXa3xVt30o7xzak9MN7TUJOrMdzgsfmmCpWxRMNk8tu8wICr8mqY0DtjODVsyTW2ELe1vigHOR3rueKJt5IBE6TJkns3xVXSaQM0Yyq+aCteSAO9F2SxSXW24faiqe9De3YCPuq6JUNvIzqc+D8UFc+wTuU5QHii7qyjXT4rqJxlu657UAH9u3uPmul2K7Sx2/FBHNTTbtOTz4qBrlA5vLW0FjA1k2+f/APkHwKUsWJOQe/NWQb+ptjbbnjnzRV1GbeFUbG8jNVoJN0yF6YIPmrnDR26HGAwrqQwtaGVpMPntUbiVDGioc471EDqVDZapRugb3fbUMZOB3rzJtbB70Qws9QOmXX1FmTnGMGqJZ/qZZJZj725ocjAr3iqqUUhiYMvBHapyztL7pMl/mq4kLuBir3RdvuOMUFduD1RhtufNOEshAizyz9+Qc5pMjnGFHNOrH6U2TG9Zt/8ASuasWK9U1DrRCFvfjsSKViR+kEDtt8jPFGX08cyqkcfPigWjZJArjGaJTRbiGa02MoB7YojTNQggsGt5IwzHgGlMaxCYKzHHmjt1mt6Gx7QOasIo2h7goxZV/Haqo4ouqS7ZCnt80VJdRtI/STIFUQ280iPIqHaOSaxVXNqJtoZY7VigkGGA4pMTk9v1p1a6S97ZS3Stwg7UlIwWHkGolPPSssENzL1gTlcD+4p9r/p+yTSHvYR05Au7jzWa0GNpLg7UB8k1ub64tLvSXt5JFG1OVz8V14zYzj52sDoitMhAb7c+a3/pu0ttG0Y6hMwUyLuY/jwKxep6m11GkbRALHkKQMVrfS0Q1HRDHqLb7ZPtUnH75rlyh+bSK0ki1L1VLJK7GNuVYd6DmnW01SbqFpAAQMitB6Vs7FvUV20W0iM4jXORj5pT60Vhrb/yOkuAOB93fmpOXeJ+Se1ZjfJJGMkSBgP3pr6nvjf3cT7drou3BqXpa0g+tFxeMFiXlQxwDXtXgiv9SupbZwsaABQfPeumdGdlmn3X0eowTyIJFjbODWyvl1GC3e/tpQbaQbmi+OPFYXaVJEp/StTpvqPdpX8PmjP2ld34pFZyyeH68PKvsJ5q+4upbdpoYDi3c9vxQcgEbsIzuweGNekkL4J/cVkRdkJ9oIFXafAJ7+GNhlWbBqCQHK7gVB/zDArkTtBKJEOGU8GitlP6WuejMkc6x2pG4KSe9ZFpSIBGeabX3qa9utPFtvKgjDYxk/vS1Y7Vmj3Sf/25otxTEWEivEhYpzxXbq6luX3S+KLs74adcymJFmUjGCKEmzOrzKu0Z5A8UZxTGdsikfNTnjIbJ81WpIPHjzTHTrV9Rn6TOM/PamBbiugZ7UyTT4pdRMAl2IByTTnWfTlra6H9ZZOZGUAs27P61fyMqcg9+1cJ3HNcB4/NExGB+muNrZ5NRVJDxLj/ADVxQx7fvRF7GwlP9QAq2ygi2l55Nn4+aoqurYwrEx5DjORVMrI2OmpHg05mltLiyMIfaV7UniAEyknIBoIvFJGFLoVzyM14TMOBTTVZnvpkSBMiNQMKPxStQFkIkXBHcUovguI0tnj6WZGPDVQzM3Dk8fNThdUuBIB7VPau3kqT3LOihQfipiKshR2q23u3gDKn2t3FdhgcsCF3A8VZf2ZtSoYFWIziqodx1G9qmuwyCGTLruHxXI5Whbcpq8wNcQvPkAjxUA7YZiyjAJohbCc2Ru9o6QPzQq8A5plLdB7ERQkhB3XNAFhrglkQDA5qod+a8sjDhSRmpxwtJKqDGTTBEHY2VPNWGd3XDHNTuIFtm2sdxoY4zVVLliAOxNF31g9nDEzsD1ACMUGu7HFWyTzTgCRi20cZ8VBSODUiSeTUyi7N2Rn4r0MXVdUBwCe9EWCSNoghXBz3o3VLOG2toWgO7cPcaBvLZrSXaXDDuMV6OV5iI2c7fzQQSQIc1OVi+Grk+wHamOPNW2TQBHWfJPigs0m9SxvRLJFvTsas1W4ivLsyQLtVuwpbkg8dqM023F1PteQIo571uX4p5pM1paaZIk9v/wBR3DYHNJnFzqV4zKnIHFO57q2jsntUXrOezDn/AFpeX+hj3WkhMp7jGcVb4g7Sbmyt9OljuYN9wCfANJLlElkZ4HyD3HxR4mtLe3+oB3Tt94zSl5CrFoxtV+cVmqM0aSK0vA9xH1IiMEVZd3L27SrGoEMhyF+KEhUyzIjZTf5pnrmnxWUUASbqyOPnNPgrtIboWMkhOIXXNIG+5s881s9RSe19PQIqAJsBbj8VjP6ialmI2HoeCFrO6kcjf2A/FU6hF0CR1PfIT48UF6UtJrq5dIpiiheRuxmtDqXp+4uEjitn6kyjLe7t+9X9STD9fCew9MXV9A1xGyhB4I70DDqktmjWyMemWKstbX+NppltHZTQPFcBduMZB/cVl9QsLboySqR1t27AHmsfuam1fpIuLC8GoWtm7w7ffhu9Q1HUZfUlw0UVvgjkH/KBWu0i8t4fTsc1wBGiL7uKz+h3E8l1dS6ZYgxtn38DA/elnHNh+tZy6uJbfZbuQekefGa1OrwafL6eiudyrIRxs7msfqBll1CcP9+85H5ryF2tOnlmKnt3Aq7YLJolLhYFLAf1ntUY7owXayOobZ4B4Nab05LZ3OgtYyqqTFtvUI57/NJtX0gWmo/RwP1mbBGOTzU/S2PTyWlweoybd3x81Vpukm+vhAkij+rn/tV+vaU2lmCDdncM1zQYI5tQhjd3SYyDGMjArUnbO9C/VGoR3S21ksHTa3zkg5+P/qs/NJuwNuABgVpvUTR6d6gjJQSBo8sMfmkEMTXupH6ePO5i2z8VasgYZZl4J/FTuo+lNgrtz4pg7R2l4ZRFtKcbSPNBmczXKvN/MA71hTL0rZwXmoPDdMFQr3NUXlmI9Yls7SYPFn7scGqLu7DyD6ZOkuMccE0w0/SLm+tgbdVV0PJ3AGrJar1hZ2ctvcx3J6bQ5wf81L4IJhI0to32nvmmE9tMlrOlwyDaf3PNJ4bh7dj0yQD3rSD57mK7Yb4zE6/B7mjr6+urLQ47NgphlHfzyaSiUSB1OC7djV0MhLrHduXiA9uTnFZ/VAwKoNxTIPbmrYowtt1dmT45qmZFSUhTlfFEWsoEUqN2x/arBFM3TEM+wAdsZq2BooHK3YZvxjGKI0JrL6kC7IwDnkVHWJkvdSEix9OAnG7GBWs61VLzxOxjgiLA/b80AwKvjBBHg0yMsWm6hBcWeJFXBOeQahqlwl5dmYIIiw5AGKwPabe/SyscA7h5riNDNeOZhgP5HigexyDV1tbz3kvSgUu/wKI5HEzNIsalgB3qKKnRbOd3ijbC6bTXnjeLcSu0+cUK0wMPSRR34PmqG0N39BojRSwbnlOVf4pXd3cl3tMh9yjGavu3uYYkiu1JO325oj039AbthqABBX25Hmr7TSdF3uFz3oq5t5rMhGbKsPFSvDDDfzdNAY8+0VW9wXTHJH5qCengMzhlDBhjk9q5eWz2coQNncM1SrFeUOKsN2zMOoA5HYmgirDpFWj5P9VWSRxJsZJOcc1Gd3kwwj2oe2BxVOwqRuBAPkiiprlpNzkkCoSDEhpkY7e2kiiZg8cmNzDxUteh0+J4jp8hfcvu796uBaGZBjFSjJXll9p4zXupvlXdgCpdYgSR4BUnjNQWXVmIlRo5A+/x8VW0zJB0duCDnNEWNncyq80CB+lyRkUdNZ2tzo5unnCXanlKYEZLuMnJx814dgRxRNtkssLYCseTVV1GIpiinIqI9FA8iM+DgealBKsEmSM/qKssppjG1tEhcvwAKi9lOkgSZChPzVVXGA83/wAe9aG20aG9aL6SfpuR7gRT7TvTthb2TMhWaXb92c4rHXFxNa6k0i716bY4OK1mNZnpktrNbawbJQpYLksfNSn0iaa4Z7NxKyA7+MYoSxebUtS6jzMh2/cT3ozTr66sp7pYkZ0GQWAJ/erh0F07TXFyZryAtADg4oG5VRqRCRt0gchcc0ztddvIpOkoDjcSBV41RIdYFzcW68ptwADirkxkVfxx3WkrcRW3TEa/eaycl1M8o3+4r2rfJM+sWotxB0rWT7m7ZGfFZr1Lo0dhdJ9EG24Gcc4qcpfVplp+qW13oEkFw381Uxg1iHAEjfAbimkFpKZ2jz79pOF88UrkBEjAjBBrFtqNF6Nib6iWfdtSNea3WgQOhmnaQuspyoPgVk/RAsxpl+9ww6g7E9wP+YpS2q36xboHcIhPux+a5e1l9GvLWBtUinlUMGjK8+OR/wDms02kfVapcx2kv8otye+KzketT9PEmWfP3U/sLu90O3/iFzbs8E2DweRmt8OMnpm9mevWn0Xph7NGLHGMnvWf0LXjoCmGWLqRyDJx3BrTWHqCw1aGTrjp7TjbJjkVT6ptdPvdLM3UTqKPYR3NdLJiZjKXf0d/Hc3sT7Lh2JCUHavNaqf5YIPzV2i6dcS3auqBlQ5YZ7in2rpZ3e1ETphBzx5rnFLbCzGoRtLvNu0fIx5NVI19Z6rDcuokdTwxFeuNQmWz/hscWdrewn4FRm1q7ePoyqpcLtB8itdGtZd3KGL6vUbRXVRhSvjNLZtOsbuexl02ZhK8m5tpHAHeq7/UbK69K28SykXMW3IYc9uab+nruxW1jtrJcyBNz8Y5rcus8eJDdW6XXq0Ws0hlDJt3f5TUvUOnQ+nWt5rGVhMRg7sHNMfTaQXF/e3d2qrPvwpbx3oT1T1YLpJ/qElSQ7do7qKlnWuknTKvcC4ZpLhmLk5wKbyaIw9MRagg97PyPxzS/Vba2gZGt5N5cZYfFGWfqW5hhgt5IxJFCOF/bFc7LDCpRkbcYkXmirbUryGRUhk6ZbiiY9Rgku7q4mj2NKCAo8cYpQ2C+R2zVlMMNQju4pSkspfcM0LbxfUyLGv3Gr7a3l1GUqhG4DABNSuRNaILdk+w5JFa9RRawdSfpsOx9xHfFaR9BtUjjljnwpGRv+as9MNEL63S1BclSZT8D/zip+uM3Cgwj2wfcan51rGY1CMC76UfJ/Hmro4B9PIoBjcDkHzQcUgMsbyMSBRl3MJ7YrGmSD91IyXxxFm2ltvHemR1CGbTorGRQNrfeO9LpGkdQWDbRwOOKncwIgRkbcCOaurr0iJDKdrb1HauRyJJcZn+34FUYq+yljt7kPKm5cdqyg+G0sriCSUymNl7LmhdPuJbW7MsRwR5qEbAShmQiBnyR+M0w169sLwwLYW/SCLhjgDJpoCEjz3btj3OCSKEYEOfnNEw3ht89MDGMV2ykkJmVUV2cHk+K1exG4vpLkDqgMVGAaHL8gqMH8V3pnDHyDyKivcVFSbcxJOTV9u0Bt2WThsgj9KuEj2LlQFkV1zQpVZGRY19x8UErlYAQbdiQRyDVBGcGiobYxXyRTrsyRmq7tBHcOo+3xRBdpqQhtHtpIldT2PkUyMceq6Wkm0RiIkHFIYLaSeQrENxAzXY2aNXTeVOcEUiqnDKxXduA7GrFh3W/VDZI7ipW6pKwiwFY/1V64t5LdipOU8EUFum2D384RThfLfFe1CyazuGizuAP3VK0a8S0ZoUfog+5gOKM1DVvrdOhgMPEQHvoArS6ntGb6eXAYYYEZqwKEuY3un3xMMkLQka9aQRocA/NWNAqLljTRy8eNrktbZCZ4q3S7f67UEt35L9qCbAfK1oPTtlE91HcyXGzZyae0gyx0S50S/N1KAYYyfcR88ULJqcMWqyT3S/URN2+BTX1drlvLY/S2jhyx9xHgVjOo/QKYBUnvVrVv8ADW21O8tZJZ7ZiIGOSp54qm7v4tQudzoFUjnHzV8GqQPpq2EsOMD7xSmRUSUhDuUGlTTGS1uoLNbvO2FTgEVp4dRstP0eNbVOvNOgyO53Eeaz1tZahdae2XIhQ8CtH6c086Zay3TIJjt3Ar3FGppZYwxDTrvrRbLoksM+KHsrJb3Ry8ZL3QPIJqlrttS1WYyN09/z4FGaJZzWl4HtmEknbHyKu/GR9nc3l9aNYlej0BhivBqWqXMNrpRgjfrO45ZuSKJj1G7/AIp0prfoI2A2496aJYafabp5AuW5O7mnKf6sY7R4rjTVW8uot0boQpI/BrMTv1biV8YDuSK3OvX5vomSz5iTIP5rBnO87u+6s2YzWv8AQ+lpfQXjSMQp9mB/3/0p3cppWj6JPYzFGcIxG77iTnH+1Zb0hqF1bXxgtxvWTkrnvWkutN+p1YSajCelNhRjkA9qZxTKR+mdMguJZLi9G2IcqGOBWw1lYLjQ+kGTpuUA58bh/tSPW/Tlzaw7tOcmLuynxWVF/dPGIRI7JHyMeMVesXwf6nsbSG6jWwwPYNyhs816K2jh0lTcxvHNK4EbMSAR/wBqFs4k1KR2uJijqMD80z9RXjvp0FrNCMx/a6ms2s8j610aPS7EXEc7PcOowCRhiaKi0+1lEbzJsdDukDZANfPo725OGaZwAMDJzWs068uLy3tLWQtszlpD5qTjdTLi7XYIG1uyNnEDMSdwA4IxWW1WGR9clhZOnISAFxjFPtYuTD6rtWZysa/bgfiuT29s+tyXUkwkPtOfit2E40svrX+FqsNxGsjuPvFG6bqFitvJ0t0V0ylEC/NV+qJrW9u4PpzsUcM1JLiL6O+X6cs4GDkd803Gpv0+WO6v9KliaPotb8s+MFv+YrPS2d+LIXTpIYWPDHmmmn6lLYLdJfIzrKvtycVGbXhLoUdkU2tGR7s98U2VosS2HUVCHYuM481RIvTY/Px8VrPTmmQ6hE948wWdTtA74FI10i8uZW6QDnuW7CpYii1tRedNYz/MZsGrr7S1srpI3mBVu/4oi706TSLaGdZAZSfcAO1DS2N/PELqWNumx4bNZzPRG1uV07UC4HUA7EUxuL6znzIsTOzjBBz3pQG+nDxsm5zwD8URpLD61JHyBEckYzmtS4G2i3B0UT3bWzAsNoDAj/nao6hqdvcaVNIxPXuCTtzwKK1r1FHfwNZRW4yT33f/AIpBYactwXaeQpHGeV71dW/4Et3hQ7Z13LRMqsLHqxsoQnsO+Khfi1E2LfgY71dbZFszKm+MD3Zoy9eagHsIrZI0GOSw7ml0auxO1WIHxRtxprJYpd9RdshOFHcVdHqEFvpnRij/AJzfcxqAWz+l6hFzkj4qmaMNORCCy+MVfZabc6jNstVLHux+KewHT7FVgiX/AKogiRm8VMVmuq/T6RPCnjiqjkHnijIzbi+lNwSV3HtV94bSVR0s58DFL0hdj25rsUrwvuU4NTkhkijBccGq8DvmqCkzMjPIRj+1DlVJG01DcRkckGuqu4gDuaqusW6gBJOO1Eta3UGJjA6Y5DFTUJECOnTyXHJzR0+v3k8LQyhcdqgXXFxLcTCSRsuPNEQWt1qcgWNGZwMkgeKCbDNkDFMrXV5bXaYFCuBgkeaCCSy6ZeFf6sYYYoOUh5Gb5OaunlErNI5JkfkmqtjbQRyTQHfw5YYo5ZZMKw7g9qokk6q9Ld/LXs3zXbZS8wgnkITFQnhaPcMfyweD80BcWrTQadJYpgxN34FFWt9YLoL2zRA3TcA4/NJosKCWGa6mSeOSORSU1XsaNsg4YVfbhGfFw5K4riL15EUfcTg0ZM6Wl2faHXGM0Ry9s44IEkhDMrDvjtQizSIMREjNM21iW8tVslt1Az9wrkP015J0Z/5TR+fmilJ3vknOPJpulzbHQugI1M2e/mjha2j6ZMqRkMnZ/mhNK0W4mhe7G0wJnjPerFwu0+wuLy5Kwwu4XlsDtRLaXLO0jW6ZEX3fg0fp3qBdMaTp24YSD/N2P9qXpq9xbmdU25uGLH96dHRoNfltNM+hMALuMb8Uw0yWbSrKBJpuolxgbT35/wBaSalpjw2UV2JtzNyRjGK0mgWvXt4L1iJGjTG0/pRuUj9RWcFrLush3GSRzg1DRtdFlLGAmSTtOa0TvGouri/t9kbcKe9Zaw0ma/kl+mUFe4Y8Yp96Lm9HXqO6kmu7RwdsYwW2HmjZNOl1qWJ4JWW3ReWJzk/pSmBEDm2viVdOOa1ujmFLcR2r7lXvWrNOmXvNLudKuGkOWgIIJ+axkrB5nPjccV9bvJ4bmK4t5Gwyqc5FfJboBbiQKcgOefmpymOdh16LD/x6LYQvBzmvp011BHcR2zHdM4yqAePmsB6Ps45NNv7pmKtEQFYeKMaK+W5tdR06Y3DSkIWfkDxiuNvYferri8ttL3Wg4OQ9YjQNUtLCO6S6iDNIpAyuea+kXkbS2OJMM23njjNfJtUG3UJRtC+49q1J/wAlNdB0SPV4rudpekUPtwe1KJ5rh4x1mZl+Sc0Tp13cpGbO3ZV65wx81o7r0ZJ0lWG53gckMMVrNS2MqIriO3ineE9Ddwx7HmtVZPf3WlylIjBGg9pxil2uSbNOhtVYYiOGA8c1q7TU7GKyht5Zky64H5q8CM7Y2l9IVnurZ7iAA+5u5HzzTKD0/p99Cs9lJsbP8xd3+lVatqd7ZLFpcACiXhJMf0/+Kt0hUtbBzYFmu956m7kE/muXLlW4Uz2kd3qf/wC2xE/T8yRvx/3orRboG/urqax3RR8HamdpA5pdeX+oQa3LLDGizbSH2DiqNO9R3WnwXNuI42SYsWyOQT3q91NPdR6Wv6dcXq25jSH7CRyfmsWqh5Iw3tV8DNP09RG10V7OGNT1M5J8ZpLNMstrFEsYVk/q+a0Voda0mHRdOhurGc7nADDd91VaTra2mmPFKGR2+18d6RtJdTIke6RwPtBOQKd6zYrFpNmkZDSDG7H6VuX7EoT6oXs6JcSkW27lsE1prvUNKbQHt7W4LPEvtGG70qu9MhsvT0UoBFxLgBD5yaSR28mXtVUrIeTXPlARa3STW6W4h6lzI2FOOaOhil0OGX620O6Q5Ukcf3pMly1pfRSxKOpA3Y9iRWr1bW47/Qmtb6IJePyqr/oaDMmdDe/UGPYuOP1qBYSXJ/nFY3+7uKYaPpMmsSNasTGYULf9qWXlhLBLLn7Y22k/pREJ4wsxjUZ/Tmmui/TCxuVmLByMKhB5qn01C9xqg2x7xtPeibeedNedOkgMbN7WHHGa1Kqm30i9uY1jfdHEze3ccYph6j9PpZ2EUluN204cjvRNzqsd5av9YrW+w4QpwKk0czaeqafNJfIDuYvzW/8AlCi0ujpsbCzaTqMnv4PFAXM0Uqb+RN5Y9zWj024USz9a2G7b7h+1QuE0mWPrBMtzuXPanvisiQOSefzRcVpcCFJ8BYnOFYkd/wDatjDoml6vYwNa/wAt4wMgeeOxrO6r17aae1KgW8bAbcecd6xymGALiKeS3FywzGTjvQ0sTRqN3BPOKm08gi6O/Kd8VFizoGY5x2oj0DqpO8eKl0WEXWXgA/NU45q2F9rLn7M8iiuK7huoOa4T1GPyaa6nbWjxLLYvkY9w/NC6bbx3Cyh3CyAe0fNCgiu04NdVV755qy4tpYJdsqkVUePFRF9uxWcMydRB3wM4qd10S5eDj8UXpd5b2sTJIgbqd8+KKtxp4ZkxkP8Acc9q16pEMyNktRcb9eMQzPtRR3ojULCCDUY4oZN8bc5+KGu7VUnJhbK1MRA20iwmQAGMec0Z6fhhnvtk4JUqcYGeaDtZykqrJlogfctH8219HNphyX7KecVYqzUmtbSN7ZIdlwGyDjnFAKyfTKJAWZj3+K0lhpVrqcry6lc7bntsUgUjudO6UUjpLlY2+2lMWy/T6U8E1u3WLKCwI7UtupluLl5wm0HuBT54zqWkQ9OzKBRzJSVdOlmkkEXuWP7jShpDq0K6SIWPfgjaad6C897ps9rblI4ewbzyPistcPDLZxQwrh1ODTXTbi0h0CVhcNFdAkYDYzxVlaiVha2ltcSQSxLcJEfcwWuXWn2cV+t90StkSCM+f2q30pDdSNNIoDiRed1D6ra6jHlLwgwqcqqjsKz9T4I1nVdN1S2ihhUxuDheMf8AarhNN6c0wJ1dwk5H70u0KbSreSX6wdx7c0NrExv2xE5KIfaPxWv1D/RcOrS3QEEzPIHOcPzinMsMnp8reRn+XIQjR/r/AOKzWiWl2bkTpHuSHk5rVSzz63NEqRAQJ7jnya1LqaEmaw1C7e7vMRnb/tTj0tEkFhI6tlWbI/ShV0GzmEquzdTscHtTC20qT6IWwlZEHle5q5/Vi76S21EPLgHcpXcK+W6zZrp+pzW6NuAbNfV9NtDYWxt1cuozgnv3r5X6gWRNZuhP9xYkfp4rHJeQz06Lt47qOBwqSj3Emn2lx3kLQ6Qk6iM+4uvcc5oH0fYG60q+kRirjgY80w9EWMkskl5NISyMUIJ7EVy1wrURS3EVyttJGZY9met8fgist6ug09tptiiPn+ZsH/etY+o20c3RMi5PB57Uv1HQrG9snhgVUZ/duBzz3q/jbqzxnvRGn9G6kvp9nSAwhPim8eiy3WozXbXBWBuY1Q9qBl0s6Vbrbw3LsznkE8CtFaRtZW8ce5pQeO3auk4Z26xnZrawjjvLe4QuzN/iuv48Gl/qLT7XTbWzks8iUtwfntWp1WCPUF+jXCycHdjtVGtaU89lBEuHkRuCRWrOlLb0yXVtYfxCIopK/wA3PnFBXE1xpGpyJZ/zIZACzU51MTwaQF1CITW4UBiowVPg8fms7FrVqYmiDYC59zDk1ztpekH1KMXD3EHBKlTSN2aWcsgxuPmiI2tZJETqMqufcaO1yxiso7f6Yk7hnce5qMoWmm2wtTNdTBXXnb80NdJA0LTQvwf6fzXL+2vIRDLdxuEb7WK4U1V0fqFaSMcg5KD4oNJ6feCz9Oz3ElvvYMfdx8CibS2s7VxfCU3azjDDHKUg+omaA2is0KFc7CODR1tbHRr+zcz7kn4ZOMVZadGvrFYxpNtJEQWRgy1kri8uNQuuuwwUXFOPUrtJqHShcrGq5Cnkf60qCXMlpI0MDbR9zquaW36XBGhaVdXmpW87RH6Yvln8VofUV6LfWrWFYciKPfu/c/8A1RPou/S404WZGDCvf5rNepL+4k1qcwvkRjaDgcVZcaySGPpqaNbfUNQf/wBx4pEzvKxvHbqdUkug/Wo6XqU1uZIVXcso5Xb5o7Qru0trx0vkCQ8na2e9XdR3SZEtYjJaT7bmU7WX4FNdb0jTYtH66sqXRCkkd2J71nbuezVpTZ+1mfK4zV1rGL65hhnmZnbuT4/FRioi5W60j6IITcGQgN4wMVbb3+o2ln0I4ysR9pPzV62RtL25V2VdoGztyKnFqX18RtOmIxHyDRviBj1CSOQQs5Cnip3t2qx9ONsn+o/NW6ctpLqTnUkKW6Dg4wG/tVWp6ZHPKr6e+Ucnav4oXwJp19Lp92t7bj/D7gnvnj/emupyXGp6e1+0i7W5ZM8ik0kRjk6Dr02AwxNcVJYo2KSkRtwRnINZztg70+PT9ZiMcmLWWMZB+aVXNlJY3T70M0R7MvmqXikgbp7va39XzWj0HUoLQfTaqAyN/hu44qNztmEfph8xkbu34qD5Za0etW1tMrfSMrEHPt+KQXYWNlCDaAORVSnUc+nt6ZMWP+oBPcVn0EiuHAPFdgXrSYDbVptfJFFaR9AgEHDeaHr2oWU9xp8V4V9mOTRGkWFvqthN9TcBHiHtzVuk6jcTaZJbzxLLarkcjGP7UPe6Ybe16i/yy54UMRxVX/SJlKTFScgHGauNtNsMqKen80faW1jNaut3J0ps+3nvRc84trGOzRRtByT80IRByqYFF9WBYlwDkjn9a87WqvN00C4Ht5zQSk5BPK5yRQXSWskY3SKVDDK580RFa3arBMGADfac1PVdTbUkhVVwIxjhcV4XrTRwW6nYI+KgHneWOdxK2ZPkVKxu0tbnqyAyJ/lo76SC7kEaNhz3OaWSxdFnSZcOOwFawN5NenMTRWiskT91A7VRpt0tvDMkyMpbs+Ku0V5YLKUfSNITk5KZqQWXV41hhtmLA4LKv2/rRSVmMZfY5IJp1p+k2suhyXk8g3KTx57URY6FHbaytpq20IyZQ5wG4+aW6jBFb6nNawTkQDkYbIpDMW2HqC4sokS3G1Qe1MbbXDq2pxw3fCOQox2pJp9lPeI5hXdsByAO9UwGW3YtGPerYxjsabU1ofU+hQWcsNxG3sfIZaN9OafpwcvIquV+3cKQm8utQhZZHaR1HIPOKosrpoEJEjZHzRrZrYX1+mmXMsCWf8mRc7lHBrmngXtvLJb3DQgVRoWojVmK3m1gi4GcCvW0lubm5tSwitGB2kHGT+tdJWLw3uHGkRuluxkkErsxyw80TcatBbAxuSCO4xST06BG0ztckxROVTLcEU0lS3ubpWMRbdwTg4q3uLOGDbS4FxCJIz7GGea+Ua87nXLot3396+pXQFhYSNCu0KpwPivkd3K0t1NLIclmJNc+RW29CBl0S+aMZbdwP704toX0zQXlASK4cFmY88k96T+gSf4NfBf83B/Y1W+nazqMJ60jva+F39/9a4/jla59Hmj2UL2geRupK59zY70NfSXmlalCkYBglO1fyf8AarLjUjo8VlbW8Gd2QV8+KhqmpLcIv1ETxLE28Mynv+K68f3J3UFXds31sDzDKSLgt/lNOIYkhix3VR5rNxXVzqhjIkAjB4HbNNdTinTTNlu7sSwBwcHFZvPlvbW2F91dg+o7UKGCEbQR2JNH31zs1SztgcdQNk/pil2o3Nrp8FqFDNcQsrdPByeefxTUpaazbLIp5U5VhwVNa20nOlWtajNpk6Q3CiSzkON2Of0rN6zbwaxcg6LASI0/mALgZ5rRPKl6wsb+VWMcuwcHLEHFHPDY+nLSWW1hO+TnYCTuIrMu+pedfMLU9G6iYrnawyDTzXrxLi8ikkwEVBtQfpSaO3mdTdlP5av7ua0M1tYeoJYhZMIGii/mFhjJAqruKNV9Q/xPS47X6bCxkbn3f7YqvRdPeZmuYJQYYiC+4YJGe1KQzW7yKVBTdtJ8GnmlLCqSSpKTE3HSU9z+lRr0bqcVtqWp28lt7YYkAkcDPOTSRr5Hv1aRSyp9o3dqMhN5pIlnmQx28j/4R54pDIUM52NkH8YpCtZpVlH6hurlWl6XTQeN3f8AtQ+iakumPc6fKgljZtpbOCP2pFBPc2u57eWSM9mZCRx+oo/TLi3s4Fu7iIzSh89s5pN3s0x0CBF1i6IdgUBKjGMj4pFq80lxqcxC7TnGAa09pINZvpb63RoFROM8FqSeoYre2njmgkDOw9w/NdLxmGg7eboL9u6YHgYomz059SWeeT2MMnBqnTZ5I9TSURdUn7lAzTDWL0FHe3ieEtwRtK81nED2CwxljNFjacA5oib6eOdHWTB+8KF80HGdsCrL9+ckH4o7SntGvGma36qjOOKA2y0iPW4GuXuenKGwRjPFZ2/PRneOIFShwXB70Rq16iXkn0W6LcMFVOKVPI5blmYHvSgoXLiEEknPfNE2hkWI3Ec3vXkIRS7cDGFGc/AFFWUwRgMZYf0/NZUyvLX+J6cdSSQB1GHjx/vSj6lDaiFoyvP35/2prFA8EwmDBYZPvi3Uq1KNUvGCqVU8quKi2LJLwPbLE6fb2Oasn1CKexWGSLLJ9rA0vVS5xRUUERgZS387+kVrETeOS2iWWKXdv7jHao23SmmzOf2ocs6exyRt8ZqBB3frQFzoiuegmR+tW6eY572CK4YJHuy7Hn9qDgnMRIPNQJ43fmg115pkMjvJps2LZSC6AcfmmXqPT7a60hLgTbOkvGOc1jrC+mgtZo4X+/grnvV0sl01qIMyMM+4ckCtbGtDW5+rQKSI2iGVOM7qqu52aQZyGxg80dc2a28SSwNllxuUUtnk6rbtm01lHoSqSZdN+eO+KKvLFrVoi5BEwyAPFCDJZQoyc8AU3ksZbWOObUWYMwzGCc4qoWSAxERrwaksEsZildDtY1GUtPdjp+4/2p/rs0cenWcceOoijeB4OKYFcsotrtZk44+2rdRigDfVRTdRj4x2oB1adg7Z2fNegLIzxbdwHfNA9sPVRhsmt5YFIIK7gcd6Ms7640PTDIkKOkp3Bycbf/us5p+nNqMskcXccgZrt7c3ccIsp2YCM/aamrtG3eoT6nqFs00g5OMquMZ/ehbm0xqgg3FwfIHeqfrz9PHH0lV1/rHeneiarY294kl9HtIH+IVz/wBqRPQszz6RLvsjtVwQyMO1VaJLDNqm26HEjc0z9WywXskVzYuHRvvwMfp3pPozR2mpLNcqCoOfBqL9ONant9KvgLCNSZFw1Z53kBbcoO7xWgnvbTUb8OsYCjscV2G00+Wzu5JhiRCQODW8LQXpy9gs5ZDJH7iMAZzmmthC+oQdRkChT2+aVenJ7CzmmlvI97/0grnit7bR213CksSBVOGAAxWuKyk+k6OPqpHm3BONqeK0cSKvPxVcYEbtgE5Pmr8BF3Ecmrf4Wq51+pgkTHBUivkOs2b2WqTQOcgHII+DX2RSojJ8kV8o9WTLLr03THbAP61m+MjPSWvJpYuIHj3JId2fitV6f9QpfTtCyBAWOz8isT6buDFcvElr1ZJRsB+Oa0F3Yvo+oR3KuqOBv2+D+KzLfHLcrVnT4pNRFy3JVcbT2oua1hnQrIgKkYwa+YahrNze6gsjybeR9vjFba31R77ToktVdpsoGLDABBGTU5c8bk1PVvo9PtY4iRChOAR4qenXoe8MLN1Y2XdGw/H/AJpH6sttQu54oZSjpjKKnk0X6d0y+toZJpCkLNgEN8ZrF/8AJbCwp9Y3GNTk6bbdqgE/tTSx1K3svT6NaMryAc8+fzRPqjTvro4pMAwL97DuB81kzo1xp7RXF0oex5KuDnFXjziYH1Fru2uVu3kJ6jbwR8mu/wAYvVk6l5KZBj2g9hVN9JCwd4GyAeDQAglkjMwUlfJqrgudrj6fdKuyCQ+KHZyjr0HZFx3BxmmWqvcrZQW8/vXuuKjqV+LuztbVLdllQAZI70XA6XcTaUYWhy4b7qosblrWXqKD7W4HjNW6hai2kgEXtdx7s0Vf6C9np0N006MrYyoPbNMXEX1W6vbtXndZFHaMDih5bRp7qdwNqx0HG5jkDr4oxLyWS86ioNo7igshmh+mWAthXOCfIqccBl320UuYl5Jr019bTXIYwlSBihotkqlAcAnNSDSem3iud0Ky7di4UA96zGqKy6hPGcsA3miPppoJlubZCUXnIpvfW8GuaX9baLtlg9siHua3eUswEaZe6Xa6bC1sqtegfae5PnNC63cy3sCXEluqbG7J/vQnplDPNPCqj3LnJ8U2srS5U3UKokw2EcH8U/XwxyzbTr7SWYbPrgv255/tQ2lSRwWkkLoBKQaUQ+67SD/BIY7qdX2nr9Dvt2BKDlj5qxCHUoelJub7mNBCm+k2y3twwnk6ar3Pg16awP1rYIMfYNUwC6fcpas5eMPkcZqESPeXRxhDRsFzDpdy4VGl8ZYcVAR3GoSMsMXtX3Hb2ANSmuXaiJVjWQmRTk5NGXNjcX0cEkMeSBycUpdMXASNC2w80XPqU8ZWOPK4Hms5WtLWVo5tvY1L3LKPxUM/NFWVx0BJj7mGBWtZev7YxCN2P3jNeilWaEQMoDHs1Rmg6MQaUEF+RUbB4472GVhwrU1pB49jEfFejQyceO9FanIJNQaVRgEcUGo3du5NRKa6fBZnS3lkcrJvxWj0Rrdry7jDKwEa9/3rI2Vuby66JIBPGT4pu1h9G0TLIHcHDY+Ksqz3Q1s1wuozw9PqSsWAHjFA3aiJTEU2vuJNFam30+qiZDlc5qOqXCXTxyKMNjmgX2shtrhJSMlTkZo+/wBYl1C4SW4AKp2UcUDLGR7vBqEaM7bVHIpqCY0a4nZohsI5xV6mW8V4zyyjNCIrr39pPGTVm+VI9pQnJ4PyaajsVy8SNGyA4NctJke8HU9scje41U8MiNmZCpPYGuwiMA9TvQaq0g0u01MywXJSMJy2R3xSHV5Prb6SRXDqOAw81SpdIG2/bmpi6ie3MJQqT/VUVyy097q5WPPejdQ0K4tFBl7HtQ9nbXd1JFaxEMxPtOe1Nba8ey1R7XUuRgAEVrjBbZaEyaed7538/pQml6Ab3V2indhEDyV70NPyZnilxHu9opx6VvZ3lWAgFcn3ZrU4y00L6k0RNHniNnIyxScNnnFPJLe00304safzXlC9+7MSKA9TXksc7pKMqB7R+KBsLcW/0+plg6jtH5WtZIoq29PSzMZ2QJ8JRUOqS6MpiaEufAz2p7Y6hbXUJMcgDf5T3FVCG2t7FpbwrcEckkZ/tVyfFnGqtHvZNSD3G7awODH4FNbS8EzOkgwVOCKyMOpwjUIjYnoRythkIq0fUza662sqg45anpjS6jcpaIXZvbjivkuozm4v7iY/1SGvqNxApsmgunDts5Jr5VcKq3MqL9okIFZ5eJWh9IXYsTPMbYyso+4AnFa2zey9Sp1buBvZxtJKj+4pR/6eRb7G7ViMFxxjNbCO2itoNiYVfwK4XTOLHW2m2Vp6kEIG6Ld2bnite6pBi3t49hcEgqOBWaKSafqL316m6NjhAKMvPUNxHCtzHaFYl+4MfFTlw5XvE2PZlt79J78mR04XPGB+KH1fVZ7rZbQIwDnkock/2rmk3qa5qrvclSioNif3pp/B7Sxke6twQwGOT4q8OO3Kze0NFfq6M9vIzbgGQ7u4zn5rOJqTRWdxpmofzY0J27uMDx2qh7+abW2S0YqZiM88fFEXmkzz37wsyh5Uzg9v711vCfFjKwQma4EYJWNz38VtYNNtbXSUtZGB65wZM9qzF9eG3tW01oQJIm27wfg0fd3cI0K0ERYzlvcM96SyA31SLdJLOKwAdkb3BTu4q71TZCHT7O+hVUliCkjHekt1dyx30BiQKWI9vfBqz1Nd3lzcW9tLICrAYVR5qVon1C8e9dZXxvAxxUHupZLdY3kZgvYE1yW1kguBFOpjz81BkVZigfIB74rKCrO0jmX3TBP1qzT9GvLuOS5iUiKPu3zQvRJcLG+4/wBsVpLC8utEijhuwv0s3BA7rQhNcWZMKSqQCTih7WMq75BOPimGqMt/dOunMTFGM4xQem3ZtbrLxdTwRmrMBVlqkloVtZiDbMff7RnFaPTbvQ4Lpo7X2pOoBJLYJ/eslfvDd3ubZOnu7gnOKpTer4LfYeCDirOg41jT5vT971YGIt5cgc80Z6f1WKxgm6kmJZBkMxpX9deav/09yTIVHtAHal8oaIPbsp3g9zTZqaZ3+lXcLC8lUqkr/wCL4Apsj2yQR2E0ytGwwJVYUpl9Qz3ljBZzovTjPu29yKL1FdKMUUti+xmGCO4FWUtgo2diZHs7NhINuThuSaRra3zOschZCDwTxRGkLPDcmaNNwQnJz3qN1fT3l6IlxnP9qWxE7uJLS0NveoDLIcrIPAq6w1SPS4ili3VaX2v7Qcf8zSrUDOLgQzPvKcZFHtd2llpKRwR75Zc5c8Y/5ms7/DEtUgl01kurcrtm5bgHvSeWYXUpZ8BvntXOpNOVRmYjPk1O6tlhVcHJI70aDL35q61jE0wjztZiADUDsWLaT7qN0iGF90s7lQnIOKiOau8u9IZ190Y4OMZFAJ/iKOwzR2qzm6kWRTuVeN2MUCn+ImOcmkUckC3V2kXUxuruoWYsrnpx8Mvap6i8HUiNqNkqgZOfNCXTztL1JuW+aLXo3kjkMxYq45+KYafO/wDMnkYnI80sX+e4DNivLI0e5RyKJBV7GrJ1klDFjkLntXLUQvC/XbDDtQpGMMD+1dRWmlCqOTVF0GyZxFJISviiVWK11FArDpke7zQ11aPYzhHIJ78VV1MSFiN36mgvv5BI+FPtHbFW280v06ADe0bBhxyKEQG4mVEXGfzV00r205VRggY/WgaXd4mpXdsHVV2jDcYpXqNmILkhGDJ4INcILR71fMnkDxTP0wthJeSC/PG32huRV9WdldtKTKkMj+3cPA7Vqrqz0WO+t4QyRoyEt7s88Y7/AL0n1q1tILqY2xABGU4oTT7C81V2K8lR3puGfB9pDNa64yWko2AnY/GAKEv1vL66mdxvCHBlA71dYwyRm4+oZx08jIFOdEeyh0+drqRfp3PGe/5qaSaF0HQI7vTJZZpiDjhAe35oLTTdWFz1baGaVUJGQpINT0e6gttXk3uzW+CFG7ArXWl2JtFf6KFQ/uCgn8nmkt+E4s6JY9T1QS6tiFQu0ZOKp1myiidPoLsSQlvsDZwad2FrYX6PHcxhbuM+7JpGrvDrc8cAR4kY7iR2q/8AX1T239LEaeHS6f6gjO4Y4/FERWF3LoZsZFAlAxvPn8060u5S4s0dSDxg4oKaa6uLiSC2jMew8yN2/arx5xe2TT07LZxSXDsxlh/pXnNN/S9hKryXNyrKWxjPGaZ319FZGG2k980zBcDzXdRMzWTC0wkgrpn8MBazDdXkki2ynaqEbvk/FfNJYmilkjkJDqxyPzX0WPVntbNon29RUyWz5/Svnt1L1rqWQ8szE5rHJmxoPS+oy2Wn3PSQgnkNjNN19U3bWqK0Dk490m3gUl9LalHBHNatGGLDIzWr0ZRdQGG5iVUzx+aRHp9dtLqCIFdwB5yvai9RtpL7SXWDGXXgfNI/UlxZWNwsMUYRguSR2P8AzFMNE1Ce/wBNHTwhUYDEZFbnmM1h45b3Rb5gA0cqjkfitv6Xv7jUtOY3cgOR+9S1zSreTRp5pSDdKhIfyTjgUl0nUNLisrcyTNFKOHAOOa5T/wBsa8AXMclhqxuUjYQwzcNjg81qpp11ayLxo8M0a5ViMEfoRQGvanay6PJHFgo/9eO1Bx397FBHHblSki4zit+BBZ3cUF/JLexfUrzyRk5p/pekQS2cuosRtbJVCfspafT99HM+8AqQXJ+R3phZ2kc+jSzxXMi9IkGMHjIrM/0hIbgWt4JpsuQ3AI8U01WD6iCHWbbGUIbZkeD8UvlSXV5gIYxvReQKqtluwstspKhAdynxUR7WtU/itxHKy7Cq7TS3jdgc5PFeOMnnODV9i8KXAM4ynkVmhyPTsjyWotZA8kq5cKw9tNb70xPHpbXV9eF5UXOxm4FQ9L6jYW9xKsbHqN9u81VreoHUIZIrq7AdX/w4zgU05fMA6KPp7O5vFGQAQfxSmOQi6EgwNxoqW4WGI2lsxKSDB/Wl5UgAnsKBhPFbxXTbW37hxgeaBZWjYq2VJ8VdFdhWViAcVdNLFcxmXABFaE9F1ZtJnZxGJCwxyO1UahcPf3r3GApc9hxQZPNSD7SG74OaxnYY2BTT98tzGHLjABGa88ELaO0qjbMXyBg9s1ZqGpWk8sDxw+1Vww/NXG/trizW2CYkfAQVtCu1vJrc7Vc7W7jNRdi1xvjbHnPajZtKa3ukS4O1mGePioXlpHA52saiwJPOXkyCc/PzUJZjKUyMEd6rbvUowGbb5NQTjkZJfaSfgUXYqtzcHrsAo+TVmnaY01/CjNhWPeo6xYjTtWNuCXBAPFaEbC3iluGluP8ADFcvp4us4tCFiIxgDvU7jMKFHG34xQEZCzK3cZqUSWSTomMg7fjFXaaIhfQ9c7UDc5GaLjM13cGK3iXJHHFMn9Ny2mnNe3BG5Rkr8UxZCXUdn1rvD9meCOKEeRn+4lq88hYn4rsLiN92M1BwQybC+0gDyeK7FJsB4yT81dc3T3AA2hQPAoWrgs3ZA7A1IAgqUOCar2nGa7GxWVM9s81EN9ZuLaS3iSCM9UD3MaTLzWjmTT7hJjCcyMPaM0lNpMBnpGt8mrFQPuQhthHc0fbWX1coMko6ZPMh8VCCwurxDsjwEHPFE/wa+Swa5jJIHDAfFSJhhbvo+mWs0boLiU/a23vxWehuNkicAIDyPxV4srtrQ3BiPQHc1C1jglfEmRx2FKp5rRsNSsY5LIjqom0gDFIrS+vNPY9CUp4PNNOlEluY7Vs+SaU3BWSQbAeO9SrTy61uMaIIo9sk8g/mHGKEtLNLuwijVGMkjY3eBQEkWIf18U8l1oR6ZBDa24SVT3xQ99A6zo8mizxbm6iuD7q1uihJtCjnt9ouVVtq5884zSa+1Ma0kNpIFEx7Y+ats/Teq2ZScXAG05KA8YrXGdr5ehMeiSX9zNNcT9G4AwTGcUhk0nUbe4eO1SSQO5Xevn85rTQNFqErQrMyYHvKnBzQ8GvSWTGzSFZRG2xG8kDyavNcT0/UX9PwJZ3kDqz8huCKYf8A6ntxcJbFWEz+ccD96TzQXnqSd1deiIsYOO1Vv6VmjAkmviHHauE9Dq401JNThuJJd8jHOR4x8V03rxak8LkNb7R7vg0q0m/W2vRBO7P/AEBmPammp6bttneyBLtztznNeqXo+lepWMTmeSAfy2XuOeawbgJK4HPOK+m6VH9JZsLtdi453V85viv8QuTEPY0h21jlGaM0SwuL6SQ2jBWQZJJxTiz1GR1WyuZdgV8M9B+l7yGzgujOSuRww8URqcem3dtDFpzmS5lOWxknNcuNuuenWlaTo+ozSrcTdeRT7RuxgUz1ezhstHmj0+2XqBe6jkD5rNQyWOj20QmhYXiHk+4VGx9TbdcM0zubZiQVySMYpz421Fel6abmznl1OaWEAezeCM1n4yfqEWPkhsCtl6r120vtOFtZjqBuSQO1Y+N+jcxyDurAmp+bFMtQjezaNbhcbgDmp2GsdCymsunuRiGUjvmpanqSa1fWqEhVUqpYjFMLLQbW39QI/WV7dE38t5q/jSeHF0txqfowPMOgygHB7kA+f2rEXN5bLKj2yuBgBs/NbzXdVtpNBuFilQ5G0AEfNfPtOsvq5SkknTTuCfNJwvFda70xqsN1dOi2+HCE54rOa1evJf3LFDE5fGPxUNOE9lrIjtZcSZ2g9waM1fSrszTXl5hR+PJrUnTXxVe+nngs4bmG4V+oPcM80onhMJ2Hv+KZ3lrfW4gaTrdN/tJztFMG0ZNWKLp4XqRpmT3ZycVnlcqJaB6VvLz+fK/RiZeOeTSa9t/pZ5LaYZmDfd4rbemtWur6KayYpBJa+1j80nh0c6xplxdvKDOh5b5qfv8AphTo+jy6rePBBIqgD3MTRerxPpjppc+3pkZLjzQmi3FzpTfXw8hWw6/5h5pr6tWLVbW31SGQfZtdQe3/ADmr9WeM9dxRW5HQbcp70OuJAwUHPzUVJHsP20Xa27xv/MiIU881WTS09PTXUscauoZk3Ak9u3/3SnUraSxvpbWQAsh5YeaMXXLy3uw0U2Cg2r7RwP7UJdzy3dy08rb5H7n5q3F1RbFPqEMwzGDzRt2IJLpZLRSqZ79qC2cMSO1SinZYumT7PikDGWe5RNm/co5HNQ1AOEiZDlm7gUAZSG9pOPzzVn1sodWDcr24FNSqemXcqRhhUSCPaRyKtkuHmnMpwGPfArx3Mnjj/Wpotgubi0AljOMds1Ge+muLn6iTlsYqdxazwwQSTxuI35GRQx4YlOB4qzkqVzcSXL7n5NRKjpg55z2qGTXhuPFZR9A9J6pp0lqUSEQzRL72IHP71nNa1K+vpJpTNi33FQoPjNJ4pXjUqjFQe+Kttg9wfpw2Fb5qYvzFMUzRHirorrE6O4yoPIqm5i6UxQ+KrPJ9vYVqIYzhdTvwtuu3OO9XQaVdTzOje1Yx3PmgbKZ47uN1bb7hk4rQXt+0s6/TPsUL7vzW4rOyoyMw8IcVBXLMGPijXeMpJnGSTmgM8YHasgyyuhBeLKecU51L1G7xqoiGaRWdsspJaQAAZ70QmnSXcbSI2EUcCpVlqz+KXNxIFifYTwRVslxqGnxtHvISTkj5oKO2ltWSbgnPH7Uz1PUzqkaNHBgwrhjjzXT4ApdVuhaGBZT0T3HmmEGsWS6S1vLZ5uOxfApf1Lee0IIVJh5oEJJjncA54c9jXNNbCx0CK/0sz2U7JMVKnPAzVNr6auLedI3iBOdzMDwaYaXfRadorQyyqkvTLKw8nHFU6b6gIsH+pugZyfbux2pPWsiPqOfTY7i2W3jVrmGQFto8DuKVatdy63cJJBDgINm3NWGwN3dtcW6PJIxLcinekaJDpkb3t8xExOQu7gft81v0y1mYLZbSYSXKtHJGdygU+l1candQWkdyUjkXa4x3z3ou40JtRnkuZ5ChK8IvakF3ZWVvcWr2sv8ANDhXTdz3q9xPBNpa2Gm6jPFdSBo1GVzR1pFp4tDfRynIbKD457UovoDqF0/0o3ooG4Dk09NrFYWFq7xggAb9wx4rUahhbi5ubiGTARMcnPLVbc6TDPOXbIl/zA0lkn1In6ywLfTr/RgEGjNL11b64kSdSjp8dqqhbi0KTyQzwb27xuKY6O5W1IuD71Pb4ovImdZI2Vhnv3qqe2YThhgqe4rchXZ5oBYztKMDac5r5NMVN3IwH8vccf3r63LDHNaSxzjAKkV8mulEd5MinKK5Az+tcv8AyOdEafBNOkkcXx2+aeWen3FrZ/xG3KB4AdyEc5pX6cufpNQEpUsqgg8Uz1S7n2zuiOkM39q5yRnWi07UdP1jSmF0Y1lI96mk/p/RGlgurkIqqCyqp5/1pVpOnrcFd77VbvzWqTU7fSrA2Cj3bSAfnNa21m2sparLJcPZpFmWRyF8ChdV0u70qUJdqAWHBByDWx0m102e6hCSsskfuwcgk17/ANQrFHtYbwy+5DtCZ75/8VjlbMal1i1miax6aRESDneKgl1NG4/mNt/WnPo9VuLqa0NqsvWTBY49gx3ofUdDNrrn0CuCrDdnPYVP12qqO2ub23lli4hjHu5zTawsIdYhjtLSYRzIm4uR3/ak7Ncafdy2NvNlHbbweDT8em77S7i2vILtUZxhiWxtzS8qYzd7Z3WlX7Qt7poznenIphp12+rX6Wt/PsRuAx55p6twtleSwXaw3jyfdMQCR/vWan055daNtagnc25SOMCrLcaw61jWGkhfTIoxMIsDqDjNJNH1i60e6d7dAxYYZTzXJriS1LwsuJV4LVPQ3gS6607jI/pYd6lm+onuuVS6vJkaMz59ynGaFgkuoI3Lu6RzeQe9H61eXccpVwBbyD2L4qvTre81ieO1jj6iwLuPIxjt3qZAuhleNW+4oOK8xlMI2FtjdxnNa3S7bTP+rtLhUzgggj7Tjms/bzW+malIpYTQcjt2qmAYIm9zsvCc1OW8llXJYDHAAq2PUEhNxGq9RJc447UvJ7iqixYmkBcY4qKBmkABwa8NxXg4rikrzQERsFZkkXdn84qc8EYi3xNnHcVTBBNdORGCzAZNQDNHuU5yOMUHkG4hfmvPDsPJyK4p5znmrJAu1Sj7m8igjhNoA71bbrHCzNMMnHt5qjHIPmrkVbhv5r7MDj80B0uuy3GliyuI1cJwjZwRSrPA8V5lwxx4rmKSYOjipxjeQqj3VOVo90ezHH3VeDCY3CjaT2YUAvTcswUZxUVLRnKnBFFQlY4WfPv7YoZ2BOfNQRd2dtzHJqQIVTxya8oXHJr0UbSyhFGSewqjsDhZASO1XSXDHlV2/NUzQPBIVcYYV0uCv5pogXJyPmojtivea7iqCNPCtPtdtuaYSXR0+IxA7t3AFLLeIy3CIDtLHANN9Z0ebTY4Z5JBIGx3PamaofUdMu7WziupW/lvnA+KAiuHjGEJw3cfNGPdXOoRCO4nzHH9qE0CrgTqcAAMKtBMlk/UVUOXYZxjFENerdQJbTLsEfZgfNFavdwStC1u22QJg4FcFnaS6K8iyD6hWzj8VjQBCGmk2u5MaHjmmWgWVpeTXHVJLJgotKImxcAK3t4zTPR9Rt9KvXmwXyNuMVqEPk11g21IFEkY2qoPei9Ev7y+upjfRpFHjCqTWcsb4trovOjlBJlsDO0Hzx+tGep5PrL1JrGR8KgyEyMnJ5qz+t/qnEGoW+l3s9vdXO4E5SgPU8VhcWf1Ng6/U5BwvmsrEclPqXbfuBy2TmnN/wBZNoggClkySMdsVP1WdD6LfmwlaSYe7PPNH6trR1oR28UW2Pdk+/7v9KRpbXNyW/lM+PirdL2wXSwy+1t4BJ4K1ZfibWzg1a2sNLWLb7kH2Z70lvdVgclrePpu49xB7URrVhDaCO7aXqDH2mg1W0kcvBH1F2+4he1dTacem9UgTTT15FBDHvTE65ZGVYxJknsQMisK08PMZiwd/A+RWvjGmx6cVt1i3Y745qew08VVZGdnBGK+Ra2ytrF0ygFS5AArXkanPpszRXgaJUIyDisK+SWBPIPJrHIaL0zcQRafdo0G+QkbTTh7uLU7IWPR2SIM/vXvQFjBNYzyvGGYtjkVoU0S2trs3cCbJD3x5pHOx3S9KtLWIMIl3gck0PqejW2qXlvcRY9jYcqe4ptt3wuOVyO9I/SmogpcWch/wJGCsexGaxy5ZcMLPVOmJpKpf20km4cYzWbN1fau6m5kZoFOC3ha2PrmZhp6DaDGfOazXpbVrSwt7qG9BMcinbgZ5rP63triu9MvZ2WvS5uRFGOFYkDdVvqHqyauZ7FllRlw7nkCsq5UTMYx7dx2/pU+tJsxu4q41vWHF/YRwbFt2e4uXUMceD3onVL65/gtulxMzTKftPGKv9BakqXr2MyjbMDgn9O3+lT9cwxwzQxQ24VEXJYD8ms73gUSreW0VvfumVc+fNGwXN1k3sCqJX7Jil0upyXdpBDJjowMGC+eKY6nqlneW9t9HHtuUxyRjFdJNigNUlhlhUFenc5Jc/PalvV9igDse9H6ksrgO+C7fdig7e0ea5SI+1XOCx7CmYg251MXtiLeZQHj7NVnpzVH06/dt3TR0KkeK5r2mWundFLe4ErMMtjxS+acSwIhTlPNY9BlzcdHUZZon6iyEsWHyaGlMbox2EMfNSt0FtLFNMQV7gCrmnt89Yw7kY8jFblAcBjMLJj+YexrxtnhKGXs3aouRLcMYE2jwKi0kn2ljxUvaCdQtRbSDYfa1BhuMGpRyYlVm5AOcUzuJdOvG3IhhcLg8d+KFoG2mlgctEeSMVWwZnJbuTzUo52hDoh4aqQfNBPGPFe70ws72I2LWtwPaWyD8VTHYSXe4Wih9vPeihR+akQSoyML81fBEn1UUd0+xAcN+KImu47cy28B6kLdiauAGeJYgpD7twzVea53J+Kukt2hRXOCG7VEQRVZc5wRUNxAODXK9V0d3HFe71yulcc1B2uqzI4YHBFcRGcEgdq6ZCUC47VB15HkbLkkmoV6vVR2vecVGrbc4mU4zg9qBjpMaJqUS3WRg5FGeqrwXF6Y9/8ALjUYAP4o5o0imiuJxs3cKPzSLWox/EG2nc5POPNa3pVFvEZImZDVlmsRjkV49zAZojS5I7GZ1vYGG4e0EUPcyKt05jQxq3agPstLjW1N5cSlFB4HxQsFrLKQ6q3Sk44plotraP7r+UMp8GmupXFtpM1tb9HMK+5SPFX8darK30AspenjDEZFWaUtqyyG6PFe1m8/impNJEh21VaxCINIzDHxWc7Q59LpPHeymOMPanOWI7Cndpc2enC4kulADNlSR4pQ2uQW2jJFZAhmxv8A96U6jKt00UaScN3J8VrxdTuLTqW63Qk3CRsKAPtya0WsBLDSLbODKVAz81K802GDQ447aQGRSGBz3NI743lwVFzG3VVfaKXjh4qsdZm0y73vErI3ginltoy61eHUpT04X9yovms5HZy3cxiRD9QR2/FPLq6TTdGiiil2XyhUcfHzU4zCQrvrqT6uaEAypG+zk5wKcLHFBZOLU8N92azjTqlytwME92/NajSFGrOZSyIjd18mtyphRawR3DunSO4dqK07R7m/MkSTGPb3NPb64i0np20CHqS8KfApNa6jc6Rqf/UQ5M/G7xTl4uRE506O6sZ5T7Vxn5rGsQHcKcgmvrV1ZWuq2Rcqpdl4YdxXyq/hEF9NCP6GIrn39LMb/wD9Odx0ucnHMnH+tbCsX/6bjNhdkZC9QYH960uqSz2tjLLApkkA9q1NjGVPU7lbOzkkPwcV8nk1GaNHhBZG3E7l4NbFP4vqto6XbRxDxxzSyKxsJ7KaKeRRcwk5OPuxW/ynhBJqE91bCG5mkkC/aC3ahCpC1bdQlW3ohCHtVgt2hXNyCgce3jNYrQTO3jxVhQiLJjOD55qc9nLHCsjDhjxTG8tJU02NwwZQORjmgVQyPFIskTlXQ5BHcVpP4tca3YtHPgyQry5Uc0lbTGXSxfFsKx4XH5qmJ2ijOyTG7g1nNJVP2OQwqdvFLPLtiBLfiobstyMk1p/TCxWO+7nAK47/ABW+M2qUGb6ddjhv/wDXNWyXUMdky7fe3aoa/PDdX/VgG1G8UA+XcKOcCrQTddJrSHpA7z9xJzVcdqzSGKVunxnJ816OMsPb28/iuXNw8zhm4HYEVMFL5X27iQOBTLTmhuY/p52AAGQe2KVHNTCuqh9p2nzUEpB0pnCNwCQCKIsEinlbrsAo+eM0FzT1ZtKOmxwbSJj9z0gSygB2K8p4qIBIz4q67iWGUpG+9fmpSGMwKF4bzSoHxg9q9x8VakigYKZ/Oav6tu8JBXaw7eaATNW291NaSGS3kZGIxxVeA2SDirI7Z5Iy6c47iiq5JGkcs/LE5JqOfirAg3YY4qDcHAoj2096kXZwFLHioZOKkiM7hF+41R2NVMihz7T3PxUrlIklIhbevzU7mznsZNky4JofzUHD/ar7VUmk2yvtX5qkDJwBmrGb2BNuCPNAx0vTJr2adLeQBUHP5pW6lHZT4NW211LauWhcqTwagVZgWYd6iq6mqbga5jipoT2HFVEMe7Bq61kEFzHIyk7GyRioISrqxAbBzijbm9jmKvHBt2gcfNFPPU93BfCNrIdVVHuI/prNoJUkjuF95Uhuee1ExXNysUsscX8g8NQ8ccs7BYgSh8DxTC3RmtaqdUuIpjEI3RcHAquJTebpZ8bVGB4oO4ieGXY42kVZBcFYWiKlgfApJIC5em0aJbBup/8AHJq2aee7tGW4yxh4BbvVmi6kmnTBzCZUx7sD7ar+vjlW5xGP5zZC/wCWt7BRptzDbJKrR72ZTzQMjs3IGMntRd1bx2/TMbglhyPihpY2JwMmpQwitIZbAukoM+PtqiWwkjhSV2581RHFNE6ZyoY011HpNZRrHLlgORigJ9NpJd3cSO0hQZbJzjijNT1K4/jX08UYZlYKDgdqaWMsNn6fja2UGQqMAfNZlbTUJ9VJUAzMd5ycba3mRRc9/JpuvpObXa4Tayj+oE9/9KCvsavq7H2xdUdvNeu5rgaqGu06kqLgnsKAk/8AfFzlZM5XBrNDKb0teWlo1xIVIX45yKM0yS0l+nt7TMNwO5DHmu2eo6tq1u9uETpx/ccc4pAlw8F6ZYhtlU8fimwrTepZzbNErPlowGGe9More39RaUksrEFR/Twc1htQ1Ca/ud1xgtgCtGzvolpZPFJ/iH+YPxxWt9JVuna/Do9lJZTq6SJu2bgeeaxFzMZ7maU/dIxbNbvXYLXXbZJrRCzouS48CsE42My+Qcc1zq1tP/T/AFWKBbi1l9ufeG/5+tbNdQt7iJjDIGxx+9fH7G4e3Z2jOCRijLfULiFwxmZQxyMHtTjJrm3WrxPJaPJLK8OBxtJ/2r59K4jEio5LZ5NbOHWgLZU1BkdWGMqKU616fggsjfW8jHqEEL45rrzmrii5vraT0ylvtHXB74pfd6i1zYRwSRjK8h6XncvtNeLHaAK43sFLfuXh6nuSMg4NXb5b+5xGxjiPyeBQSQSOhcLwO5qUMjIwBOBUBlyk8VubYzGRFP2g8UvxjvRErlGJRiQfmqt4wSeTVwdeBo0V2BAPamljHcwWH1jnNuxK4JyD+1LZ7p51VGwAtER6jI2m/RHGzOasHbhIpbdpAw357Cg4Zmibco5xiuHA+0kijINPluk6kY9oFAXoVzbC6KXirsfucUBqXTS4aKA5jU8H5qkxsJ9nYioSqyOVbuKarrRERq+e/ijDqJbTRaGNeD91X3SwSaXE6HEgGMUtkiaNRu80wQwT25xXlGTluwr2cE17OKIsJjL1WwyeO1crpqDnfvRUEMMlux34cfNC10HHY0BE6wiNDG3u81VDM8LEoTg96hViwSGHqge3OKqq2YuSfmuHPmu+a6x3GojyDJxnFMNPsWu5wLdtrL3YnFLyMGr4Ll4ImETEE+aqi9RlM4Ill3yqcFqWbTUwCSCT3qyRAoUCsooUlGBHep78ksw5NEw6fJNDvWqxayEE4xiqKdhKFwOKtWOV4gcezPeoINwKE4GatUzqpiQZFBGdUT2kYIH96qVsMv60fOkTaegJ/wCoz2oN7eSN0Qjlu1AVqFskccc8bAq4GQPFNLmfSrfRYks0Rp2HvJU5FIrhZIm6UjEkePiqyQU7nPxV1R9hOz200DOQpGcDzWi9Dx2y29xPLtyp5yO2KyMEnTJ+DTGIvawDoSHZMff+lWInqTR6jqN1JGBhT7QKA064FnerLNHuVTypHejbmK2gvIjay5THvz80NdxqbkuT7G5qYp2NRE1lcXFlp6rF/UxUccUrmsuhawzyDmXkcYzTvRJoX9L3luzqhO7B88ihNS1SG9sLRdgSaLGF+eK0tIJwVkBAOKd6THFLYTTzLyOBkVfcenkWGKaa56ck/OzildxLNYxvZq+VHc08XFdzdgwlFGeeGqExIt1dlOT24xQ8bDcqsO5plfTtKY4iFwoGMVGDtdWt7MWxtosED3Aio3uuLNfxzQARuo5IH3fiqNKU6qhXohTGvf5qlITEsoeL3gkCt7Vcv9WW4uw5iAZxtJ+KovNIltijxyiTeODntXBpVzelnjTG0iitStruxtkMo2dhuNZrSWgawdNeWCZcu+Bn/n60skgK3heT27jR+o6WLFba5Llml5NKbq5Lv+lYS/w+ttBtJbAXpnB2klh+hoHWbuG+lRYwVSAYAP8Az8VPQNPu9RSfoylYV7jPBOKnbaDNendGPcGwTXT4Gnp/Uo4dHeOKP3qpyAKxFwxlnldhgs5NbhbST07FK0wEglUgY8cVhpX6kzyD+picU5WYlRUc8VYcCPjv8VGM96sjheXdsHauXiHekaTc6jbssSA7Rnk1RdXF9p0X0M7jap+2qtP1m+00n6eYqTwRtB/70NeXk19KZZzvc9zim3R2cR3cymEBTjkdqYaPbaebK5kvWAkUYQUvsrOSZ1fJRQcE0z1vSbeyulhhm3swyADmiwNazv8Awm5jB9hJxSxUOM+M1fNa3FmMSo6B/njNSt1lkiKLGSPnFVEJEDYweMV4WrtE0igFR3rzKUYxnIJ8UXNbXdlZ4nhljEnYspAxRQlsI+sDIePNGrYNeh2s0GI+5NQtvoUtC0il5s9s1KwvzDcMiP0oX4IohfkpJ7hyDzRjTTmIzRe2M+3vQ11sNw5Q7h80zntoo9DglWYFyQWjBqhcyl0Lufd4o7TdGbUYZJI5VXpjPuPel8028jAwAO1eS5lhyYpGUHjANBxy0b7R9yHGa66SzKZW5A7nNSdW6IlP9feqgzY2hjtPjNBAmvdxXT3rynANQTV1ClSM1CrktwbdpCf0qEUZkbHgeaCsjBroGaKtYFnmMIwSftJNRuLSe2kZJI2BHnFNA34oyJZvoXIcCLPIoVtuOO9T3EQ7A5xntVFe3j8V1RkjIzXWfKgeBXY32+O1Ba8XXn2Wyk8dqi69AlSAW8j4r0VzJBOZYzhj+KjI5kcu3c96Kgp5o2OGKVDmTBx2oLYSu4dq8oJxg85xiogyK5EMTRjvnvVMkhcHL5qNxbSQbeoPuGRVWCasUw0mRILtJZV6ijxTuC7t49de4MYMZU+yssGcDAJAFS60m8sXOTWtyAu+k+s1SaWEbU3ZA+K7ZTy/xFZSNxQUCJWDcHFTVypyGwfmsT1BeoSi8vXn243eKH+lPUCZGSOKbWGkxvpkl5dSbFHKc96UtJI8xCkkjsQKombCZH2tHn/5V6XdDGyA5FELqk6RGBu54GQKDcSRvskzz80Wq8jZ3oorGkSSj3HjiqTbkt7fNdCr0fvwwOKI2egwaVb263UrpvfuG8UJew6N07xsjqbsxkVn4xgxor+0/dzVupQxfUbYXyR35rWt705OLl7aOSS4DMOy5q3UbR49LgnkZSz54FLoo5ZJumhJokxy7DE7s+3sD4qIDV8jbjmpCLJyzd+9eQiN+RmnUGjXF1p7X6gLGAfaR3A80xAjTm3VRbP+DV0esTpC0DgNxnNRjexW0KRgvKe+M1VFYsHEgBbyc+KpB2mT3N7cCCA7FflieO1M9V1m3FlLp92vWnU7QQOMg96zMd7JaSuEbGfgVG1YXF0WlO7nPNNblMIb5XPRuSTAo4BHaglhWW/VxGehu+PFSvp1uJ1VF2dh2p0L+ztoViMQD7eSaiWbXrvUV0RmXTCrLKvuH+U4qn0nrEVrqchumISQd/g0BZ2seoaqqI5WIt7j8UbqXp429w4t2zHgGpy8TPrValdabr1lcWkcwLhdwJ7ZFfLZEMczxk52sRTZNPvIWLe+MAEgjjxShshyWOTnmsRlKEF32gcmrj1IH28qT4qqCTpyq45xRF1dG4lEhGCPitCy1tWuZjGTtb81TOjQzMh7qcZFO/TumW+ru5luDHKvYA4Joy99NTWVtcSKwkI7bj3GaZVwgtGuLkmKHjAyaa+nBatcSXmpS7mjGADUNL0W5VTcXRNtGexJ70ibiZo1yecd+9MwaX1Xq9pfxQxWqZKnO6lEE1xE+0jaCOARUptO+ighuXcMxIJTv5phq91Y6hBDNCBFJGuNoGM0CNi6XYkI9wbdTzW/ULarYxQG227By27PP9qSQvtl6so3D4NSluRh1RQA34rKxelnCdP6u8iWhCqwvj7jR2l/TPBMLmTBA9o/NDWlstxeiGR9gY/ca0yhPCYwrgZVhVakswTJweBzTLUiUH0+PanZvmlmdv60Bt3pF1ZwrNKn8thw1L/0ppc63dXVgtrMwKL481VcaZLbWkVy+ArnjBq3/AOwkjRRIDtPIqfQYYkZCIz5NMorGe5hFyYWkgRe9D6pfvOscOAsaDsKYBVgWR3Kt7QMioKyCNlK5PzXA2xSF81FMMcE4qDiM3254Pii7iP6W2VVcOX5OPFDOmxsGvRIZJApOP1oIxyGORXXuDmjLnVJ7k5cAcYoWeMRSbe/5qs/rQerw5NWiBjb9bPGcd6roOGvCpYphotkl/qMcLsAvnNUABCcHHFSmi6YGDnNaD1Vp0GmyQpBhVf4rPyKB9rbqWYtQSQrkY4Pio+c1bEoOcjtUGxk4qI88skgAY5A7VwZXnNSCFU3H9qhnJ5oNX6f9O22oaabu5uBGMkYxWbvVSO5dIzlB2NWWzOYiizsq/5NxwaikBuZSiDsM1NqoLCVjErL7T5qMSCSTBOBXd7f4TP7Ae1XBFguEkzuRSDViJ9SdSsTlgi9lq/TDHNfGNsRjyTRGt6pb37RtaQ9MoOW24qqzhs/4fJcTyYuM+1cGrTQuoAfXMsZyE81G3gn1CcqnLKKiVkiUuykFu2aZ6BcRwiVXTDODtcVKoEWUyo5kbpsvYHzUxDE9oq4xMDzVuqS3E8StIvsQkA0uV2J5z+opEE28fRcySjKiutD1HxE25j2NQt+vPugiUsT4q6JRC8UaKUuQcYPk1poKoktLjMikEdxRsWpxxxODHuLec9qL1fTbqJUmviAXHHNIGADEeKeJ4bXtgIdPt7wOHEhBK/FMdT1zOmxW1odilcMKz9s8krpFuZlB4XxV9wh3Mpi2nxTaupWt2lqjYUFm8/FX21+xLqo3Fh/alIGc570TYz/AEdwsuzf8ipqSoSKWdnZf2pxpGmfUabJNGD1AO2KAWbfcySFPbJ2XFMoNVnsbbbGpXjBGKsahE7MtxhuCDg0e5Wdd6LnYPmgZ361wX4O40YD04nVfbURXBcPBK/T9uR2pt/FGawIbd1cd8d6RW7oLqN35UMMr8it5NqWg/RhgsRbHCiPBz/arO/T4z9hdyX0UiSnZ04z+/FZdv8AFbzgmjZJpGnmaMlFOfaPigO2az0ykvFSBHzXa9QFxSnT5Y57a4Bf4HimN56nuLyNUcgfOK9XquijU9dmvYI4GYBU/wAvGaVKwVsjuK7Xqvova5DkdViR8V23MBmBkbC16vUF1+1rIQYGwB35oJNhfDvgV2vUsiuMUDnY3t+a6suGDZ5ByDXq9WWRlxfx3EY3cMO1QCW7Q+1iXPf8V6vUVc2nowjjtpRJI/8ASO9e1C21GyjSG7Vwn9INer1LWs61bp2s3cVjJZIV2EHuOaXjpuSZjjHx5r1eprIqxayWSQTH2lTtoBwvUbZ9ueK7Xqo8jlHDEZxVscg+pErplc8gV6vVBdeiKWQPEQq47UF3rteq0eyQNuTj4rucrjH716vVB7aVIzRulTpDdFpG2YHBrtepBG9vJLqQtJKZMcLmhEEgbKDNer1UcLHJHYmpGIjb5yM16vVBEtxtJ/SrpTbi0QIcy+a9XqAdM5yOaOsJJQXeNBgDmvV6qoJihmJbgE81KUgkhGJX812vVBdbyLt6brgHzXZIZrZg7qRGxyPzXq9VZGXtxHeNGHXo4HFVJI1tZ52A+7hq7XqVUbnUhJbiM4NUGaIIuDivV6oOQXcttc9W3cBvkiu3Fw8s6TBhvBzn816vUNFXOrTXwWK9fciDAwMUtkI34ByK9XqqiLR2gmDquanc3UjzEkAZr1eqgVsYyTzR9pJYiwkWYn6jOVr1erNDS7vNMGjxmIA3g70ueWeW3VpAoRvIFer1VBE9vpI09WSf+fnkZFLurlmjDZz2r1eosXWtp07tOrjaRmirqOwW2ZUcm4Y8LmvV6s1aVjfbM4cZJGKDKnOTXq9UjL//2Q==\n",
      "text/html": [
       "\n",
       "        <iframe\n",
       "            width=\"800\"\n",
       "            height=\"400\"\n",
       "            src=\"https://www.youtube.com/embed/VcggiVSf5F8\"\n",
       "            frameborder=\"0\"\n",
       "            allowfullscreen\n",
       "            \n",
       "        ></iframe>\n",
       "        "
      ],
      "text/plain": [
       "<IPython.lib.display.YouTubeVideo at 0x7fd38b4df400>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import YouTubeVideo\n",
    "YouTubeVideo('VcggiVSf5F8', width=800, height=400)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "extreme-bobby",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<div class=\"image-grid2\">\n",
    "    <img  src=\"../assets/B37_fig421_a.png\" class=\"medium\"/>\n",
    "    <img  src=\"../assets/B37_fig421_b.png\" class=\"medium\"/>\n",
    "</div>\n",
    "\n",
    "<div class=\"image-grid2\">\n",
    "  <p > Zero circulation in a fluid at rest (<a href=\"https://www.mheducation.com/highered/product/fundamentals-aerodynamics-anderson/M9781259129919.html\">[source]) </a>\n",
    "  </p> \n",
    "  <p >  Starting vortex observed behind an airfoil in a fluid (<a href=\"https://www.mheducation.com/highered/product/fundamentals-aerodynamics-anderson/M9781259129919.html\">[source]) </a>\n",
    "  </p>       \n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "handy-uganda",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "Section 9.4 in {cite}`B30_1` presents a concise and clear explanation of the development of the bound and starting vortices and the interested students are encouraged to read the relevant couple of pages for this reference."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "improving-tucson",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "# 3D aerodynamics"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "requested-canyon",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The validity of Prandtl's theory of circulation was established for the 2D airfoil (infinite wing) case by experiments conducted by Bryant and Williams {cite}`Add54`. The corresponding verification tests for the finite-wing case were conducted by Fage and Simmons {cite}`Add54_1`. This is shown schematically below where flowfield velocity measurements were carried out in 4 different planes (A, B, C and D) before and after the wing. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "defined-remove",
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA13_Fig95.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\">  Schematic representation of flow profiles behind a finite wing \n",
    "    <a href=\"https://press.uchicago.edu/ucp/books/book/chicago/E/bo12024403.html\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "extended-announcement",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA13_fig9_8.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Schematic of the experimentally measured circulation around an airfoil section <a href=\"https://press.uchicago.edu/ucp/books/book/chicago/E/bo12024403.html\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "stuffed-paragraph",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The contours $C$ and $C_o$ are seemingly unrelated as far as the flow field is concerned. However, if Prandtl's lifting line theory is to be believed then the circulation about the two contours should be equal. Additionally, since the circulation about $C_o$ represents bound circulation, it is related to the lifting force generated by the wing - "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fallen-cowboy",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "$$\\Gamma_b = \\frac{2VSC_L}{b\\pi}$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "adjustable-communications",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "where $\\Gamma_b$ is the circulation due to vorticity bound to the wing, $V$ is the oncoming flow velocity, $S$ is the area of the wing and $b$ is its span. Based on their meticulous flow field measurements using a hot wire probe, Fage and Simmons showed all the above assertions to be true. The two experiments of {cite}`Add54` and {cite}`Add54` had helped establish that indeed the *theory* of circulation conformed with reality."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "productive-acrylic",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<div class=\"image-grid2\">\n",
    "    <img  src=\"../assets/BA14_Fig147.png\" class=\"small\"/>\n",
    "    <img  src=\"../assets/BA14_Fig148.png\" class=\"small\"/>   \n",
    "</div>\n",
    "\n",
    "<div class=\"image-grid2\">  \n",
    "  <p > Tip vortices generated at the wing tips of a finite wing <a href=\"https://books.google.de/books/about/Applied_Hydro_and_Aeromechanics.html?id=Ds-bd0zAwIYC&redir_esc=y\">[source] </a>\n",
    "  </p>\n",
    "  <p>  Modelling a finite wing using bound and trailing vortices <a href=\"https://books.google.de/books/about/Applied_Hydro_and_Aeromechanics.html?id=Ds-bd0zAwIYC&redir_esc=y\">[source] </a>\n",
    "  </p>\n",
    "\n",
    "</div>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "harmful-corpus",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "<center>\n",
    "    <img src=\"../assets/BA8_fig5_6.png\" class=\"big\">\n",
    "    <p  style=\"font-size:16px; color:gray; text-align:center;\"> Source of induced velocity on a finite-wing using Prandtl's lifting-line theory  <a href=\"https://mitpress.mit.edu/9780262526449/flight-vehicle-aerodynamics/\">[source] </a></p>\n",
    "\n",
    "</center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "91d5fffd",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "As {cite}`B0` puts it - 'Blade element theory is essentially lifting-line theory for the rotary wing.' (Section 9.2)."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "complex-qualification",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "## Final Comments"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ordered-switch",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "The above, seemingly dissociated, set of concepts are nothing new and have been covered at length in fundamental fluid dynamics courses. These would all be used in order to build the framework for unsteady aerodynamics analyses. \n",
    "\n",
    "One thing worth emphasising again before moving on is that the Kutta condition is an engineering condition that is imposed on the flow in order to obtain deterministic properties of lift and moment around an airfoil. It is not derived mathematically from first principles. While it has successfully been used for steady airfoil analysis, unsteady airfoil analysis presents a challenge because experimental observations suggest that indeed the flow can briefly turn around the airfoil trailing-edge when rapid motion is involved. So the Kutta condition isn't strictly applicable for unsteady flows {cite}`A29_R51`. Neither is it helpful for cases when a sharp edge doesn't exist, for e.g. the airship-like body from the tests of Fuhrmann and Prandtl above, or when more than one sharp edges exist for some random profile (i.e. not necessarily an airfoil-like profile). However, early unsteady airfoil theory was developed accepting the validity of the Kutta condition even for unsteady airfoil motion and has been shown to provide acceptable prediction of unsteady airfoil lift and moments.    \n",
    "\n",
    "Potential flow theory-based results remain meaningful only in when the flow remains attached to the airfoil. In flow conditions where there is large scale separation, the linearity imposed on the solutions by the potential flow theory introduces inaccuracies. If you are still wondering, how on Earth does a flow model based on zero viscosity give results that are representative of viscous real-world flows well you are not alone since some of the brightest minds of the 20th century struggled with the same ideas. This struggle is chronicled in {cite}`BA13` as a clash between the British and the German school of aerodynamicists. The theories proposed by the latter school were eventually universally accepted once supporting experimental evidence was obtained and the entire saga has been artfully concluded in the final passages of {cite}`BA13`."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "thrown-index",
   "metadata": {
    "slideshow": {
     "slide_type": "skip"
    }
   },
   "source": [
    "```{admonition} David Bloor in 'Enigma of the Aerofoil'\n",
    "The German advances in the understanding of lift and the properties of wings depended on the use of abstract and unreal concepts that were sometimes employed with questionable logic. Progress in aerodynamics thus depended on the triumph of falsity over truth. Everyone knows that false premises can sometimes lead to true conclusions and that evidence can sometimes support false theories, but the story of the aerofoil involved more than this. The successful strategy involved the deliberate use of known falsehoods poised in artful balance with accepted truths. The supporters of the theory of circulation showed how simple falsehoods could yield dependable conclusions when dealing with a complex and otherwise intractable reality. This is the real enigma of the aerofoil... The enigma of the aerofoil is the enigma of all knowledge.\n",
    "```"
   ]
  }
 ],
 "metadata": {
  "celltoolbar": "Slideshow",
  "hide_input": false,
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.5"
  },
  "latex_envs": {
   "LaTeX_envs_menu_present": true,
   "autoclose": true,
   "autocomplete": false,
   "bibliofile": "biblio.bib",
   "cite_by": "apalike",
   "current_citInitial": 1,
   "eqLabelWithNumbers": false,
   "eqNumInitial": 1,
   "hotkeys": {
    "equation": "Ctrl-E",
    "itemize": "Ctrl-I"
   },
   "labels_anchors": false,
   "latex_user_defs": false,
   "report_style_numbering": false,
   "user_envs_cfg": false
  },
  "rise": {
   "autolaunch": false,
   "backimage": "../assets/tum_bg_top_right.png",
   "enable_chalkboard": true,
   "footer": "Sumeet Kumar, PhD candidate",
   "scroll": true
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}