{ "cells": [ { "cell_type": "markdown", "id": "e9f9b7c8", "metadata": { "hide_input": false, "slideshow": { "slide_type": "slide" } }, "source": [ " \n", "# Momentum Theory" ] }, { "cell_type": "markdown", "id": "6c33d237", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "```{epigraph}\n", "\"Fortunately, engineers and technologists do not wait until everything is completely understood before building and trying new devices...\"\n", " \n", "```\n", "

W. J. McCROSKEY

" ] }, { "cell_type": "markdown", "id": "d009ed8f", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "In case you were wondering that helicopters have been around for almost a century, surely engineers know everything there is to know about them. What contribution could you make to the state-of-the art in this field. The above quote is particularly relevant." ] }, { "cell_type": "markdown", "id": "5a36fbfd", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Objective" ] }, { "cell_type": "markdown", "id": "e7399701", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "First order of business is to develop a framework to quantify a relation between the power necessary to develop a unit amount of thrust. One can imagine that this would be a highly relevant excercise in the early days of the invention of a rotorcraft - *is rotor-based VTOL flight even allowed within the laws of physics?!* The simplest framework that can be used to model the action of a rotor, which is also the framework that was used/developed by early practitioners in this field, and quantify the power consumption is called momentum theory. The idea is to model the dominant flow physics, which in the case of a rotor is a *push* to the flow opposite to the direction of the generated thrust, and evaluate its effect. As the course progresses, more and more from the rich feature-set exhibited by rotorcraft aerodynamics will be introduced along with ways to account for/model them. " ] }, { "cell_type": "markdown", "id": "d7bd8568", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "The objective is to detail rotor analysis model called **momentum theory** and derive a relation between the following quantities\n", "- rotor radius R\n", "- thrust T generated\n", "- power P required \n", "- induced velocity $v$" ] }, { "cell_type": "markdown", "id": "dcff4a39", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Why is a new model even necessary?" ] }, { "cell_type": "markdown", "id": "b7d2ba11", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "Alternately, why can we not just solve the Navier-Stokes equations (NSE)?" ] }, { "cell_type": "markdown", "id": "669709aa", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "An aerodynamicist in the early 20th century was very much familiar with the fundamental conservation laws governing fluid dynamics - the Navier-Stokes equations (NSE), a set of nonlinear partial differential field equations in space and time that should, in principle, give you the flow behavior for rotors (in hover, or forward flight state) much as they would for wings in forward flight. However, analytical solutions to the NSE are available only for a handful of very simple flows (or rather simple boundary conditions). Instead of waiting around for a powerful enough computational device to be invented to solve the most complicated form of the problem, aerodynamicists at the time solved the problem by introducing assumptions/idealisations to the rotor flow problem. The result is the momentum theory whose equations are much easier to solve. " ] }, { "cell_type": "markdown", "id": "832fc75d", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "Momentum theory-based equation are\n", "- simple and intuitive\n", "- solved using a hand calculator\n", "- match true results *pretty* closely" ] }, { "cell_type": "markdown", "id": "20667f81", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "Turns out the dominant flow characterisitcs of a rotor are indeed dictated by these simpler equations. If one is seeking the relationship between T and P for a rotor of given R, momentum theory is a very simple alternative. Momentum theory provides us a simple framework of the flow phenomenology of lifting rotors. The flow problem is idealized so that it still models the dominant flow physics but simplified enough so that the number of problem variables are minimized. Notwithstanding how computationally intensive it is solve the full NSE for a rotor flowfield, staring at the end result or the equations that led to it is not going to be very helpful in understanding relationships between, for e.g., thrust and power." ] }, { "cell_type": "markdown", "id": "287f7bfa", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "```{note}\n", "Lower level models such as momentum theory are great for all the reasons mentioned above. But with great model simplicity comes great responsibility (on part of the user) to always keep in mind the limitations/assumptions used, and whether they even apply for the problem at hand. One major pitfall of having a simple model is you will *always* get an answer. Whether it is appropriate or outright wrong depends on the application context. More on the application context below.\n", "```" ] }, { "cell_type": "markdown", "id": "bfb0e0b5", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Fundamentals of momentum theory" ] }, { "cell_type": "markdown", "id": "5dbcc40a", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "How would you go about quantifying the amount of work done in order to accelerate said mass of air and keeping the helicopter aloft? Thats where the momentum theory comes in. Sure, there are issues of whether the aircraft can be controlled, structural stresses are not exceeded etc. but if there isn't sufficient power to even operate, everything else is pointless. To evaluate required power, we adopt this model which was in fact used first for analysing ship propellers (ships were around way before aircraft were!). \n", "\n", "For this, the rotor is modelled following the Rankine & Froude actuator disk model. As per this physical model, the propeller or rotor plane is assumed as a disk of zero thickness that supports a pressure difference/discontinuity but not in velocity. More on this below." ] }, { "cell_type": "markdown", "id": "1ad8df13", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "
\n", " \n", "

Heirarchy of physical theories [taken from \"Understanding Aerodynamics\" by Doug McLean]\n", "

\n", "
\n", "\n", "- ideal representation of the rotor \n", " - no losses\n", "- mass, momentum and energy conservation \n", " - ...but so do the NSE!\n", "- based on Bernoulli’s fluid model \n", "- results in closed-form expressions\n", " - not possible with NSE\n", "- simplest rotor analysis theory!!\n", " - and quickest\n", "- dominant rotor flow physics modelled\n", " - some effects left out\n", " - tip vortices, viscosity, $v\\neq f(r,\\psi)$\n", " \n", " \n", " " ] }, { "cell_type": "markdown", "id": "a091349f", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "### Model details" ] }, { "cell_type": "markdown", "id": "a6bae87b", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "
\n", " \n", "

Hovering rotor flowfield model in momentum theory

\n", "
\n", "\n", "- models only thrust, no torque\n", "- rotor as an actuator disk\n", "- pressure jump across rotor disk assumed uniform\n", "- velocity uniform and continuous across the rotor disk\n", " - continuity of (incompressible) fluid\n", "- rotor details \n", " - only radius accounted for\n", " - number of blade, chord, airfoil etc. ignored\n", "- control volume analysis for average flow velocity \n", "- flow is axisymmetric, incompressible, steady, inviscid" ] }, { "cell_type": "markdown", "id": "fefb14c1", "metadata": { "rise": { "scroll": true }, "slideshow": { "slide_type": "slide" } }, "source": [ "### Mathematical formulation: Hover" ] }, { "cell_type": "markdown", "id": "bb79091f", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "The details of momentum theory are discussed here as they apply in the context of a helicopter rotor. In hover, the rotor generates thrust by propelling the surrounding fluid in the opposite direction. This requires energy, since the kinetic energy imparted to the surrounding fluid has to come from somewhere, and chemical energy of the fuel burnt onboard the helicopter provides this. In addition to the energy needed to 'push' the air (to hover), energy gets expended separately to overcome friction encountered by the rotor blades as well as imparting swirl to the flow (this takes additional kinetic energy). The latter is just a by-product of the rotational motion of the helicopter rotor. There are different variations to momentum theory, we'll just be concerned with the simplest version. So, we ignore the viscous and swirl losses and consider the rotor operates by just pushing the air down." ] }, { "cell_type": "markdown", "id": "f7190b12", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "A control volume as shown in the figure below is considered and relation between T, P, R and v is obtained via the application of mass, momentum and energy conservation laws.\n", "The derivation below is for a rotor in hover i.e. " ] }, { "cell_type": "markdown", "id": "c980b92e", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "$$ \\overline{V}_{freestream}=0 $$" ] }, { "cell_type": "markdown", "id": "869ee73a", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "This includes vertical as well as horizontal velocities!" ] }, { "cell_type": "markdown", "id": "c8843ff3", "metadata": { "cell_style": "center", "slideshow": { "slide_type": "subslide" } }, "source": [ "\n", " \n", " \n", "\n", "
\n", " \n", "

Hovering rotor flowfield model in momentum theory \n", "

\n", "
\n", "\n", "$$\\text{Mass Conservation:}\\quad \\quad \\quad\\dot{m} = \\rho A v = \\rho A_{\\infty} w $$ \n", "$$\\text{Momentum Conservation:}\\quad T = \\dot{m}w = (\\rho A v) w $$ \n", "$$\\text{Energy Conservation:} T v = \\frac{1}{2} \\dot{m} w^2$$ \n", "\n", "\n", "$$\\Rightarrow T / A = \\frac{1}{2} \\rho w^2$$ \n", "$$ 2v = w$$\n", "\n", "$$ v= \\sqrt{\\frac{T}{2 \\rho A}}=\\sqrt{\\frac{DL}{2\\rho}}$$\n", "\n", "$$\\quad C_{T}=\\frac{T}{\\rho A(\\Omega R)^{2}}$$\n", "$$P=Tv=\\frac{T^{3 / 2}}{\\sqrt{2 \\rho A}}$$\n", "$$\\quad C_{P}=\\frac{P}{\\rho A(\\Omega R)^{3}}$$\n", "\n", "\n", "$$\\frac{P}{T}=v=\\sqrt{\\frac{(T / A)}{2\\rho}}=\\sqrt{\\frac{DL}{2\\rho}}$$ \n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "b96efff0", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "The power estimate obtained above has been obtained without assuming any viscosity in the flow. So the power consumption cannot be attributed to any viscous losses and indeed is just one of the components of total power consumed in a real rotor. This power component is referred to as induced power. The quantity $v$ is referred to as induced velocity, i.e. flow velocity at the rotor disk. This quantity is not just a coincidental by product of the formulation of equations shown above but a very fundamental physical entity to quantity accurately if one has any hope of obtaining not just power consumption, but also aerodynamic load distribution of the rotor disk. Of course momentum theory, in its simplest form presented above, deals with the entire rotor" ] }, { "cell_type": "markdown", "id": "11ffbfb6", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "The following observations can be made from the results obtained-\n", "- flow attains half the speed from rest at station 0 to the rotor disk at station 1\n", "- the other half is attained travelling from station 2 to $\\infty$\n", "- a larger rotor disk is more efficient\n", "- operation at higher altitudes requires more power" ] }, { "cell_type": "markdown", "id": "ae67e85b", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "Here $\\frac{T}{P}$ is called power loading PL and DL refers to disk loading. Naturally, a high power loading is desirable and one of the measures of how efficient a rotor is. Disk loading is a fundamental attribute of a rotor and in-general is to be kept low - from the expression above it is obvious that this leads to high power loading. In reality, having a extremely large rotor (in order to keep DL low) is not always a viable option. " ] }, { "cell_type": "markdown", "id": "8f96c386", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "Applying the Bernoulli's principle, a relation between velocity at a point in the slipstream and the corresponding pressure can be obtained. " ] }, { "cell_type": "markdown", "id": "704664d6", "metadata": { "slideshow": { "slide_type": "subslide" } }, "source": [ "
\n", " \n", "

Hover slipstream pressure distribution\n", "

\n", "
\n", "\n", "\n", "$$p_1 +\\frac{1}{2} \\rho v^2 = p_0 = p_{\\infty}$$\n", "$$p_2 +\\frac{1}{2} \\rho v^2 = p_{\\infty} +\\frac{1}{2} \\rho w^2$$\n", "$$p_2-p_1= \\frac{T}{A}$$\n", "$$p_0+\\frac{1}{2} \\rho w^2-p_0 = \\frac{T}{A}= \\frac{1}{2} \\rho w^2$$\n", "$$\\Rightarrow p_1 = p_{\\infty}-\\frac{1}{2} \\rho v^2 = p_{\\infty}-\\frac{1}{4}\\left ( \\frac{T}{A} \\right )$$\n", "$$\\Rightarrow p_2 = p_{\\infty}+\\frac{1}{2}\\rho w^2 -\\frac{1}{2}\\rho v^2 = p_{\\infty}+\\frac{3}{4}\\left(\\frac{T}{A}\\right)$$\n", "\n" ] }, { "cell_type": "markdown", "id": "e5408a4d", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "```{note} \n", "Momentum theory, with its uniform inflow velocity assumption gives the lowest estimate of induced power. Its not due to the idealised flow conditions (i.e. neglecting certain phenomena that contribute to disspative power loss) but due to the uniform $v$ assumption itself.\n", "```" ] }, { "cell_type": "markdown", "id": "cab47732", "metadata": { "slideshow": { "slide_type": "subslide" } }, "source": [ "### Aside: T/P, DL, MTOW" ] }, { "cell_type": "markdown", "id": "733ac4b5", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "
\n", " \n", " \n", " \n", "\n", " \n", "
\n", "\n", "
\n", "

DL vs MTOW [taken from Helicopter Aerodynamic II by RW Prouty]\n", "

\n", "

DL vs Speed [taken from Helicopter Aerodynamic II by RW Prouty]\n", "

\n", "

T/P vs DL [taken from Helicopter Aerodynamic II by RW Prouty]\n", "

\n", "\n", "
\n" ] }, { "cell_type": "markdown", "id": "b15c7150", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "As is evident from the figures above, large helicopters tend to have smaller rotors (relative to their size) and thereby higher DL when compared to lighter helicopters. The rules of physics apply to all helicopters equally, so based on the above momentum theory derivation they must consume relatively higher power in hover than lighter helicopters (and they do!). This is because overall helicopter design, and the rotor design in particular, is not dictated by the sole criterion of having minimum power consumption during hover. A helicopter should exhibit efficient forward flight performance as well and having a smaller main rotor facilitates a relatively smaller fuselage [^SmallFuselage] which leads to lower drag in forward flight. However, the DL vs speed chart indicates that the higher speed performance for helicopters has grown in tandem with the increase in MTOW. This is because of the presence of relatively larger engines on these designs that are necessary to deliver power during hover when the rotors have a high DL. \n", " \n", "[^SmallFuselage]: In the case of a conventional design, for e.g., the tail rotor cannot in the wake of the main rotor since ingesting the distrubed wake of the main rotor can lead to fluctuating performance of the tail rotor. So the tail boom length will grow in size together with the main rotor. In order to keep the helicopter CG close to the rotor chaft, the front of the fuselage will need to be expanded. This not only makes the overall design much heavier but also increases its parasite drag. " ] }, { "cell_type": "markdown", "id": "84656880", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Summary" ] }, { "cell_type": "markdown", "id": "31496f36", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "The rotor power obtained via this analysis is the best case scenario pediction. In reality it will always be higher." ] }, { "cell_type": "markdown", "id": "0f56ee2f", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "- no flow losses assumed " ] }, { "cell_type": "markdown", "id": "e7b84efb", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "- induced velocity is uniform over the disk \n", " - minimum induced power " ] }, { "cell_type": "markdown", "id": "2b300b59", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "- it makes the rotor analysis problem computationally simple\n", " - in contrast, NSE implementation is anything but simple" ] }, { "cell_type": "markdown", "id": "432361a2", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "- provides basic intuitive understanding of thrust-power relation " ] }, { "cell_type": "markdown", "id": "77fe558f", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Relevance of induced velocity?" ] }, { "cell_type": "markdown", "id": "40c9052e", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "The flow velocity at the rotor disk, $v$, is called the induced velocity and is needed to accurately estimate blade aerodynamic loads.\n", "The flow velocity at far downstream, $w$, is relevant for accounting for different rotor wake interaction phenomena - rotor-wing or rotor-fuselage interference effects, brownout, whiteout etc." ] }, { "cell_type": "markdown", "id": "42b5f5bd", "metadata": { "slideshow": { "slide_type": "subslide" } }, "source": [ "### Detour: interference effects" ] }, { "cell_type": "markdown", "id": "bbea3c8d", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "
\n", " \n", "

Mi-6 [source]\n", "

\n", "
\n", "\n", " \n", " \n", " \n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\\begin{gather*}\n", "T = W + k_{dl}T \\\\\n", "T - k_{dl}T = W \\\\\n", "(1-k_{dl})T = W \\\\\n", "T = \\frac{W}{(1-k_{dl})}\n", "\\end{gather*}\n" ] }, { "cell_type": "markdown", "id": "0372153d", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ " \n", " \n", " \n", " \n", "where $k_{dl}$ is called the download factor.\n", "The induced velocity $v$ is flow velocity at the rotor disk itself. Flow velocity impinging on the fuselage is not going to be very different from this value. If a lift-augmented design of the Mi-6 is considered, knowledge of $v$ would directly give us a quantitative understanding of the *download* on the wing. Similarly, the download force experienced by the helicopter because of fuselage can be roughly ascertained. " ] }, { "cell_type": "markdown", "id": "37c8c802", "metadata": { "hide_input": true, "slideshow": { "slide_type": "subslide" } }, "source": [ "### Detour: brownout/whiteout" ] }, { "cell_type": "markdown", "id": "138acb99", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "A fully loaded CH-53K (40,000 kg) would generate *hurricane force winds*" ] }, { "cell_type": "markdown", "id": "c17cc60d", "metadata": { "cell_style": "split", "slideshow": { "slide_type": "-" } }, "source": [ "$$w = 2\\sqrt{\\frac{T}{2 \\rho A}} = 37.6 m/s = 135.4 km/hr$$\n", "\n", "\n", "
\n", " \n", "

Induced flow velocity in the far wake (taken from Helicopter Aerodynamics Vol I by Ray Prouty)\n", "

\n", "
\n", "\n" ] }, { "cell_type": "code", "execution_count": 2, "id": "13284896", "metadata": { "cell_style": "split", "hide_input": false, "slideshow": { "slide_type": "-" }, "tags": [ "hide-input" ] }, "outputs": [ { "data": { "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEAAUDBBANDQ0NDRANDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDRANDQ0ODQ0NDRUNDhERExMTDQ0WGBYSGBASExIBBQUFCAcIDwkJDxUPDw8VFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFRUVFf/AABEIAWgB4AMBIgACEQEDEQH/xAAcAAACAwEBAQEAAAAAAAAAAAAABAIDBQYBBwj/xABEEAACAQMBBAYGCAUCBQUBAAAAAQIDBBEhBRIxQQcTUXGS0xdSU2GBkQYUIjKhscHRI0Lh8PFisxWCk6PDQ2Ny0tQI/8QAGQEBAQEBAQEAAAAAAAAAAAAAAAECAwQF/8QAJxEBAQACAgIDAAEDBQAAAAAAAAECERIhE1EDMUEUImFxBJGhseH/2gAMAwEAAhEDEQA/APxkAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAB0lD6G1ZcHT+Ll/9BqPR/W9aj4p+WdPFl6Z5RyIHYPo8retR8U/LBdHlf1qPin5ZPHl6OUceB2Po8revR8U/LPPR5W9aj4p+WXxZejlHHgdj6O63r0PFPyyUejmv69DxVPLHiy9HKOMA7RdG9f16HiqeUe+jWv69DxVPKHiy9HOOKA7X0a1/WoeKflgujWv69DxVPKHiy9HOOKA7b0aV/XoeKp5R7Hoyr+vQ8VTyh4s/RzjiAO39GVx69DxVPKD0Z1/XoeKp5Q8WXo5xxAHavo1r+tQ8VTyw9Gtf16HiqeUPFn6OccUB2/ozr+vQ8VTyjx9Glf16HiqeUPFn6OccSB3VLouuHwnQ8VTyiU+iu5X81DxVPKHiy9HOODA72PRXcevb+Kp5Rf6Ibn17fx1fJJ48vRzj52B9Fo9EFzLhUtfHV8ktXQxde0tvHV8k14M/SX5MZ+vmoH0v0L3XtLbx1fJD0L3XtLbx1fJL/H+T0nlw9vmgH0v0L3XtLbx1fJPV0LXXtLXx1fJJ4Pk9Hmw9vmYH0x9Ct17S18dXyQXQrde0tfHV8kvg+T0ebD2+ZgfTH0K3XtLXx1fJPPQtde0tfHV8kfx/k9J5sPb5oB9L9C117S18dXyT30K3XtLXx1fJH8f5PR5sPb5mB9L9C917S18dXyQ9C917S18dXyR/H+T0vmw9vmgH0v0L3XtLbx1fJD0L3XtLXx1fJHg+T0ebD2+aAfS/Qvde0tfHV8kPQtde0tfHV8kng+T0eXD2+aAfS/Qtde0tfHV8k99C117S18dXyR4M/R5cPb5mB9L9C917S28dXyQ9C917S28dXyR4M/R5cPb5oB9L9C917S18dXyT30LXXtLXx1fJHgz9Hlw9vmYH0z0K3XtLXx1fJD0K3XtLXx1fJHgz9Hmw9vmYH0z0K3XtLXx1fJB9Ct17S18dXyR4M/R5sPb5mB9L9C917S18dXyQ9C117S18dXyR4M/R5sPb5oB9L9C917S18dXyQ9C117S18dXyR4c/RPlxv6+aAfTI9C117S18dXySxdCF37S18dXyDPDL01zj5eB9RXQdd+0tfHW8guh0DXj/APUtP+pW8gnGnONPZ1HBpUxe1GVE+ha4pI9ZEj1hBMD1Mn1b7BsQROMiCRNl2mllOoi+KEycKjQDE2eRZGDyScSmk0icSpIkmTZpaRWpHdJQiJRZuFUlgvyVuGRsUOROJPqyWMDezS6m8GlZ1FweuTLgmzb2PZ8G+BjK6g06GyYvGmBi92diOEjY2dSykakLbKOUzNPm8rGWeBqbNz92XHt7TqK9giuOzsm/LWeOy1PZueQPZx0NhS01GXbJ8ifyKX445VbKRF7Ja951cbRBKzTL/IqeJxk7Rk3ZnWuyS5ITvbLs4G5/qE8Tlq1DCE5aG3dWzRze2bV51ba7OR3w+aMeLZiE0+DT7iQhG3aWmVkrk5pPVm580/Wb8NaZ4znZ7RnF8fg0SjtyXNJ/M3zlZ8db+QTMmG2M/wAv4jNttBPR6P8AAnI42HwRFM9DKTYNkcg2Qe5PAAXse5LEVFkWFSI7pJHqZCaQweSJkWXaVFgokkj3A2mkMEZE5Fcmc8nXFOghuDFqIxE8uf27Q1RjkcjDArY8TQZyyumo+GRiiQddIkpP1c9x05VtCVMI0T3rFzTRdbpSeMl5I8p0h63fIZhYLGhD6s4k5bUtcW2eBUrfHE2I2r3coQrLUvIKVkU5GKk0eQgjcyTSEZF1OZ5uE6VDPAuzSyEz3KLI2T7CqpTJuLpJUy1IV3iXWA0Yyet+4pjJ9o3TeSbFUdT3qi/d7iVGGqJyQ3sjZrby1odbRs2loR2HTW6jd6nQ4ZZ7XSmzhjQ1KDM2OjGqM+w52taX3ESjDHoRzxJzt0ORpC0HqaM2Lw8D1FnPK9tSL9w8cT1sGYuSyK2U1KeS5RJQpFxzNM6tZKSMS/2c1y07WdhGgjytQTO0+XTNxfO7rZkuKEoWzTPok7JGFtjZzWqOuPybYuLj9o7I39Voc/tKxcDv3onk5D6R1c6dh6PjyrFjCjIsjVKsHp3YP2e0HH3rs/Y3La4UllfLsOULqNxg1tm47dXkDl5Xz7X82Ww2tPtXyRZlGODo0BmWW1k9JaPt5fujSyNs3HSSPUewgDQNJZPUQJKQSPWeM9QNGbWtPIkyMUTYtJFUiVKnklGBZBGcq1IrnDBfQiTihyjQPJlXaR7axGynBPePNll26SPzZS2lJc8jVLbLXIurfR9ReIuT+H7MVnsCWcRl801+h7OeNXjTlDa6fHP5nTbGoKS3sHK230emuP4HbfR+k1FJo5fLlLejS2EMDsIZ4j9C0TwsDFzY8H2E3EqiVHEe85y8p6/0NzbFLTOWc8od5qakQhc0H7vwI7jX+GJbRvZQfFtZ9z/Me2dcznww+zRfpg1ZY1FsKW8PW6cUQoXG6/tRXv4/1HpVotZx+JN0EZtoQryxyOh2TBNc/ihi/s4cXgzyHHSeSVNe46uns+D4YGqGxYrUtzXTlqdu+wcoWb7Dr7a1T5Ibo2qzwRnmacxZ7Kb4mjb7Ca/U66x2ePwscnPmcWBs+ywakKeNDTo2OAVucrk1Iy5WmRy3ssGhSpF8IGebWisLXJZKhoOxiSaM80052tS1HaccIsvLcWo6EyrUM7xJorchmmY21pWixFypkZxERE9wRcQcTUqaeSYlc00x2RVURqZFjk9s0sJ8ji9sWuctH0y+t85OO2/b7rfYz3fFm45RwdWlgrcTYr0M8iipbI9fNz4s6SKUsmmrYhKjjkOSaJOGD2lHLL5Uy6wpajkaLdUzY2ZeNYUtV281+5JW+uWX7nYieRLi2LZZRGvAx4Xrg/d2DlK73uBqZfrFxWNnkqqIzo5I/VzVzjPBdCqTVQrpUy3qjnlnG5gkponlEoWrZGtb4M84vF62XW8MltjZ51Y/Stjnlm1MSMqYxbzwh5W6Kalvh6HHKtwvF5Lksk3TwSp0jz2tx84VrTf3WviWQ2dH3MzLa0cXob+zLXmdWrUKezUaVlZlkaY7COBWU6NBEqlElGReXaac5tvGMHM1prgv2O5v6CZl1NjxlrgtrUjj6+zIyWH+Zo/QzYqc3jgjep7CR0GwtlKHDmcss7rRoldfRmMtfmc/f7GxLTRLkfSow5GTtDZmXk3M7+ppz2yrbCHLvZrmuBsWmzsGrQtVgzPku2tOW2TsNc+JtUtldhrQtcDNKmLmaZFtsjA9bbOSZp06ZcoHO5dtaLU6WC6KJTInO5tSPWzxo8J5McmuLxFkUVZJxZja6XokimMiyCLtzseVo6Gd1OpqsVnAu1xUwoZHKcCrgTjUEq3a4rkexmV1ZltZkqMmeqRS6hHrRK3pZVRBs8dUWr1cGpGbEbjmcV9JaD1OwlMwfpHRymer4qxlHC7uAaG6lLAuenk56UuJ5u95c0yqomTkaVSgW2NDUgkaOzqZeSaPUrHKJ1bbCHbRaFVxU7TnyNMW6t8kre1wSq1XklSqm+SaTlPARuSNSGRSrEvJGnQeTSt6KMWylg1bGvrhmcsjTYtLYunZk7CRpqBy5taY31fdPaKNKpTKqVDUnNdF3MshDIy7QsVHBm5EjPuKRCmh6rTM+4jg5Wtxwdnba6o2aNI8pUy6J32ynTiWbpXFnqkZ5K9kTjMg5EXInL9a0umynB5KRCUyczRqjI17R5RgUpm1YT0M3I004spuKb5EI1B2gZuWmpNlqNA0begewQzSM814pRok4W5Kmy9MTNLEHSBxJ5ISmZtgqkRR6e4OdydJNK5Iiy2SKpoza1Az3eKlI9gyba0Zgy2MimBGU8DbnZumshkojUIbxOTMwMTiVTCVQqnM1yaxlDkVzqldSQpKqXG7b0uqVRedyL16ghKep1xx2xWv9YFbq6EqlwVb+Ttjg52tChdZFtpUd5EKNMub5G8eqlcvdU2hWMXzSwdDf23cLUrN8PeddsFaFJPkK3dDB0lDZ+EZW1rZ6iVGJOkhmxpYM+vPBbaXaNtOjt8Yxkpr0hKNYvp1zGkK3lJ4FaLHK1xqQlUwXZo5QjoIXD1H7aomheulkmzRakPWUsi8YdhKyqYlqLdpp0exqrzh8DpKa0M3Z9DRPBq0zjksjxQCdMuRCozG2tIJBJHjYZM72mlM6YjcRNKchScRaunIRJIqcifWHTZpPJ5vHjkVtmbWl28RcymUytzM7XS6cyqUyvrCuUjPJridt5GtQnhGNZGlGZmZdnFq2epoUjI2fI1abM5ZNSGoTGKczPcySrnO5N6acKpdGqZEaxYqxmZrcdtXfIsShWL6NQvNi4r90MHkpgpl2zqvZi9aWC/JXNZM3JrG6Vw1LaMD2jEZUROzLJTKJTcLQdcRS4hoWxML2opSL8i8KZeoHLHbpkhNkJEqrFrmrobkWI1quBWvLKFa9chVraG8UyVTqEJMUq1inrz1Y47ccqZlVCm8sSlV1G7SodbGK0KUBhU+ZXSjkfoWzJyTReNEs6pF8qWCtl5GkFIzr+kPSiRWHxLKmmLX2SpLhqIvYOOSf4HWQlgYism9o5OpsxcMClS2Ue07SrSWDDvrDLeBKOZlQ14lV3F8sMe2jabvNoUjbt8/wLsTs8kq1ElSotdgy4sxaI7Ih9rDXE37LZaznGX7xfZtHOHg6iwt9DFoqoQwuBM0FRK6lAypKNc9U8kuoPY25jVVFMmkT6o9jAsQvOBQ4YH9wrlAcTb531iKalQyXcEuvO/FNtSFU96wzo1iaqkuCynJSK3IqhVCUjzZ42O2N2m5EWyDZKDObR6ykOJmdRmXfWkjWMv4lsbezJGpGRzFjtHDNilfJ4M54WLjYarVCPWFFa7i9MlfWHnzljrNU+qx67gyLy83e8U/4v3Gvj+HLKbhllJ9uot6pbWuN2OfkcVU2m+0tW2vs4eTtP8ATWd1jyxtXe2Zc3jHJGpsbbUZ6c/zPnVzd5HNkXyi02ds/ixsc5n2+pRmSwY+ztpRktHkf648Vuuq6cPRtMsUxFVi2NQnPTNwpuKI1aeSFKqWqZ0lljnqxVTpBKBa5IXr1hJI1N2lruRnXVXQZuJmVfrKM/bvOiM6upRXrkdRa4lg9Hx4sZ1VVqFXWEZsqbPVHnq2dQYsa2GZ7kWwehuyI6e0vFojoLauj53Gsaezr9po5ZYLt2ly1gzpV9eBXG73lkVqV+JnDFbpoKaKZYM53TK5V2dpixaf63UsldmbCuWb+TbMOUrsujUyZikW06uCUS2ps/eWhifU2n2fgb9vtOMc51M2/wBq6+4zZutbLVafYTtllCdxf6ha7R3c8ycTbq9j7vB4TNe1r6tdh88t72Tlk63Yc21r+Jz+SahHSwPZRKaMyamYmS6DpEHEsnMq3hyNJRieSiCkeSZdmlbRVItkyEiWpp8QyBVGRLeO/La6WxkTVQX3z3eJyTRunULXISiyxTFnKLLpepE1MW3j1zM+KRbnaY64hKqLuRB1Tcsn0mjkao1TvWY/Wlkaw5ReLWV0dBsS+jutSxpw7fccSqwxTrGMpjn1VluLW27eZfwMeNYruq4q6wtmM1F1vtodYVVqpQq5GpVyYuTUxe9eW0qpn4wW0qhJmtxb1htBx4M6XZW12+LOFp1By0u8MuWMynaTKx9KV1wGKdc5Oy2kpRS4ND0b/HFnzfkwuF7evGzKN6V3gKO0TmJ7S7C6jKT4GdtcXURuslU64nYyNShS7Udcd1yvSp094rr2+6h9pISvK56McY52se+Ri3jzk2No1smFWkejCONpSrIWdUuqsXqM7Ss1ZvEt4T3jxVsPjgB6LL7eWpk1b7d48O0Y2NfRqfd+Kaw1/nkXjftHR0K+hNyKqESbZm5TE09weNnrkUzkJ8iaSciKrMrnIqq1C8jRiV1gou9oaGfO7M+5r5NodnfClS5FZTKnVM3JqQ59YJU6opGWSaZJd9o2LG8UdeZu7H2297X5HGqoP2l0k0zOUlV9To18rJYqxxtpt7RJGls/aeXg82WNjcsb1SqRjVFZVDyNQxtdHlUJdaIxqEnVJtbDMpkXUFJVDO2xWeiRrHtivk0pBGYvKsDmd6aNQqFsZmepl0KprHIsPKQbwqqh65mrkzoypHrmJ9aSdQnLca4rZVCipIk2eqJi3bUilTLoyPZRTKnHBitLXMtdUro0c8S5Wy7TWO2Kzru4E51zTvLFPVMyLuGBl26YrI3gfW8CUabfBMdhsiUlnh7jnpuaiX1onC5M2tQlB4kmvyChPI7XTZo3GRilVMyCQxSkalYsdBYXGDrNkWXWrP4HB21U6vYG09zGGy/LhzmmcMuFbn/Btzii+lR3dSFHbe+8SZpW9LTU8WfwXG9PTj8m0bTimbdHgZE66XBYK47Rxk6YYXFyyu2vcSRjX1bAO6M2/rHWM60TvK/Ey61Q8uq2os5nfGOVezF2WzZSbQNnNbcvk5dXuqSa4qWqku3Dysdh0VSGVg4vaewJwk92LqOSbTjo1zy0tcrhpxOvw6t7Zy+ulX0anOc+rc1JrOFLg2v5c5Wvz4Hd7KrbmkqW6/dJZb/593efc2cv0YbGkqsnUh93XMspp8nhr8T6TcQWMNJrsepPmznLX4Y40q9ppcVKPfF/msovtq6lqnkThav+VtL1XrH5Ph8BmEMcUk+1Hmys/G5DMplMpldWoKVaox7Smp1hG6rFNSsK1qp2mmdPJzFKlQJ1ROtMlyWQxOoVweWJTrchy1kc5e9N66Np4RDrj1zTISplyy9Mye1kZ5JxmIzljQlSZjyNXFqUa+Dd2LcrKbfA5amxyyud1nS/1Rj6fRo3SxkjG8TeDk5bQbSWdBjZ9xlrt7zjPi9tXJ1Ea5465TTloimvWxwJwntOVMzucGdd3ZVWuDIvbwSd9Ja+db571wjGWT2C7TTrGnCoWb4nTkT3zeLNMqqSjVyZ06pbRnnBMq1MWjk86wrlU0KWy43pmzs3CedB6nRM6zqYH6dck1+mScqBX1GpYqob5bYTbyRW6hKrVEq8zGVakWTvMMjGG+Zlw8jWyq+MoztrXTao00hlVEjO+sHsYtm+emOK29aa5GDWtdW1objjyaL6dFPlkb23OnM7rGLSpnibv1COcnleyj2YLMUucKUGaVtWMuNHAxSZv+zOm5a19e7U6bZ+1HjVnI2THVWJlOTMum7d3mSiFyZTrnnXnO4VvbbV1xKa1wmjKlc6FDrlmJtbWKJM8dQg5HaOe3uTxIjkspCrE4QH7OkU00M2pmocprBTcMvRTXMZLihGWERbyQUiurUObpoTYndIvlMoq1jWNrNZ9aYtXqF83lid1SfI6/TKmrVFa1bQXu6mBWpcnO10kTq1iyhdmVUrldOvhmb7b06ahcDlOqc/b1DRt6gl7YsP1Z4KnXKriroJQrGcqs+mtGoWQqiEZltKRuZsXFoRucDFnd4ZmOR5QkaxzYuLtNmbW0w/geVbrLyc5b1OGByVwYz1skNXN2ZVe4K7q4M25uMHO5adJiwaPAhUmZ7ucHsavM9GtDZtY6EakhOjdaFc7osnTOuzLmW0KxlyrDVhNM55x0acZ8zydYpnMVq1DWtTplrwqaDVOrgW2bcRcUufM8r1NcIxLsNSq6k6lbArTp9oUlhkahyLZXdw095a6xXKfaSqzpvBfa9vaV3kUJTr40Eit2lVQ/bzycvZxm2tHum1SqY0G04tFxZ5TlyKqdwV1bjGpYH3MhKRj1Nodh7SuWzcykZ4NacUyuGgtRqZPHCTZ1xu4xY1oVi6NQzqNPBfSNRk3vBkoyS3i8RKZVOWCbmRkSweQmTZXFJEnUIPSyEinfJRYs6D9EeoLCMqhVL3cM51ZGnOuKXNcUncC1xXONrpIajXJykZdOpqNdaZa0uqyM+7mX1Khn3LOmNrNjynLJ7NlNFajEqZ205sfa9BtZMGpQb5M7Ca5MXcUng55OmNci7fiuYnRTydZeWq1a0ZhXVPHIzpuVKjVHrWpr3mB17TNnZFdN5Y0VsQoZWvAz7uvjgjWTyiipYZ1XEljEpGhcZRqWUk0IvZkuzBXs+o1LdZFa6jgFTzqXIlwNY9OeT2EsHk6p7FHjgbuMrMui1bLM2tRbkk00nzNSpkFV7Tlfijc+SuFtdmqKw8t9rk/wBi2Vsl7u9m7Oy7yt2ef8Hbe1c7UwueRSpV14nSVNkp9grW2Cn/AJNy4jFnXCjfuLHan0exwb/Bmdd7GmuGo1KrTjtDJbGehzNtSnGSytM+86KDM3HV0G9nT4jtCX2jIoXO7oatGtpkxobVCQzCKMKF2X07/wCRmkP3bwKKMpcEFOspPXgaKuYpaYM7bVq3UY/axkQhSU5r+9Cd5d5fuHrCkorPMb9JTsKuEK16u8L3Ndp+48+saEqyFK7nHVYwVqjKeNfkaVOnvlit3HVlXeidez3Y9p5aRbNbKawyhU1nQtSVCjReUP4wVw04kamvA1MkW9YXUpCaly4nrlg6Ss2Q91hCVQWt8ltSJrluMcdLesRGdTIjGWHhljqGZlVuKUqpD6wL3csCyqFtXjGpGsX06hk0pjtKRqXbFh+EyUqop1hJMtm03pPrTyq8lUg3jhcHTkplPBZCvkhOORerFo55Y6dJdnZVhavVKFVKq1Q1jUyN2s8jXWGZQfYeusdOVkc+O60aiyI15cicK3IK1OXwOdu25NEak8iVyzU6v3GLtaEk1plGWoy762zwPNmUnnCL4s9sbvcfuZva10dppFJ6scoywV20N6KfasjNChniSzbltbF5Eru3Wd5cR+UUhOtImU0kKRuNR2FbQRr086riUWcJ72Hokc5a6aa6kRVQsyit0Me83tz1Fm+ii5joUXDcT2FXOhuZb6Ys0qaZGL/vBrStCp27NNbZ3wIj8qJDqfcXS7IygLV7NP8Ayac6feVSp/20T6WVgXOxc8GxSpZTjyTOqVD++BXWtM82i8j7c7ZbPctMYf5Gk9mOI9bUXDP8xf8AW1/kky9pdsapHCKYVljHM1LuCaMC8p4zga21DtunJ4XzH608afic5s29lCerST7TVrbTTkkc7j23KYjW4d50FC2eOJhqopNZxhf2jVp1Own0lG0rV4M6nVwtTXrvKMiu1qnx4F+1lPbPvMIYdbeOUqVGnu54mnYzxxYvQ2Wmjz60k8CVS70EaVxvS7hroblSvkthLBi16r49hfG70AZqXSQfWjNtJKUmmO11jga2hindtjG97zBpXeGWSuWyypYburhZDrk1oyNOimtVkupW6XBD7q9JUK2dBqMOxCio41WC+lXNzpzr2pSxqi2nDtIyqFVWsa5M6MVJkqZnqZbSngvJePRwGilTIqpqS5JItnoLVpl9RFNOhn3GLNtS6LOJ7Stt566IbqU0iFS4RnjprdMxikYt9W+3hdw/KoZO0KUt7eWPcS5dLjG1aW+OOow2Z1jWeE2+I1G55HPa1fJZMi/hy+RrJi9xjgzVqRzUz23slvLHLU1K1iuQuo7rJKu2xSaSLoVdDIpVy+FdvgLWLi0W8lNWJXCaRLrhctmtIqOCqN0s64GN3IrUtFnhk57aNqafA8SYjUilqsolTvs95qXaa6SvXozOta+o1f3S0QrGGTf70y6hp8sv35IPPb+Ry+x9qYl99yclnD5d2NMD9TaUnLdaSi8rO9qaxwtNxqPPaQ337jMtrnc+xrJLVN64XZniJU7+Sbkszi+3l7kb8dNxvuR4+4z1tDKzjDxwz/QUsNouTe8sLt48CautjXbXYR30J3F3GPF92Suhdbyy8x78YFmjbQVRdjPHKIoqifd7gUtRoMSUfcUTs6cuw8mREUle/R6nLhLD+ZmVPorJawms8s6G+4B1OS7q7c1RsatLLmt7Pq5fzLobQk1z7joOr7GyudHPHD+CFkq7I2W03zEds1ZZ3uRsStE+WD2pZqSw8Y7MkxkhtgWs01l8RuF6s6jH/A1wTlFfNfuZu0fo9N6Ka+TX45ZZjKcmhUus6cEW27ws9pj2mxasX9qUZLvefyNijTkljGvzJcV2lc1dMGdQuGm1xLbqjPjgzZzcJJ/MnEldRTtXFKXN6kK83jiZS+lGdMd43QvN5aIWa+yUmqrT1GaNXLQpdx3pYPOsUcLma0u3R20+Q+6uDB2RPjJjVW8WUkY3pDs6xXF88ile6SRnO+1NRnTeldLmxWpdGPK4b1Iu60L9U4teFyWyuUjD6/J5KuT/AA1p0lK6PZXJmbPqDF0s8DNpJGja3GRmTMmxlhZGHcFl0zYun7yiPEVu7l72ESp1sGWtJ3csalM5b+OJHadXMewTp32ES0jQo5WgfWnFlf1pNIWqfaeuUZqtantJCdO8zJ5EYRSZZti4Shnmmv7Yo1esE5Pelgy7LaO9ojfsrDTOdSYpS0KHFDNCnhYK3Pdk0yxVcloYgtC2nDHEgqpTWukiVk06+CmtcpCUauSFSnvcTLWkL28S19wna9vEp+kdvuwbTM3YW0NN19pZD8dfSw1wKK1Pd1FKFzgfVZNHSON6c9siu45hu8OHY+7P7j/1jMknF9uez4lFnTnCW402vWxoOXC3VnktT1a/qqK6knvaYx/M+f8AUhcSS4vC4Y7TKvr2bTcfu55JtsetHvx1xn3mvqxClOruyaXCPDKy38R+3vN9e/n/AEF/qUmsSnh54o9js/g3JtrmmkvjxLlFJ7aTn93P2f5cZz/gr2HWnPMZL7HB54/uicJ4k/tReW8JZWvvHLehJNNYjn7y7TOffX7EhhzVNP1VwQrRuJb+9FNxfN8EW3C33iMsbr104jNaoor+hNbWLbio3HRtP3FlrdprRqWOLM+0hlye9vJ8n+RG4WIt01qs8P7xlC47XbVubndWcZEbPau9LDi1px1/DKE7G/zHEmnhZf7YGrGSeJe54S4Y967ScYbrRjW1/wBOOzXJfGS9wokCJcVlWXV7GGN58eBapZ1/ISq0m+z5FlNPn+wuEibMM8a9/wC4s6f2s7z4cORP7WuMPszjUcWtrlFniWCO80llavszxPXN4zh9xNG3suHIjOimvtLl3lNa9Slu89NC+TWnHUmtVSM9l0uO7juIR2fH+VyivmNzrRyo5znsw+HcXS0WnzNW0lZ9lstRbec+9nlzZfaTWHjl2/M0Os/l+L4FlPtM7XbHvqs1ooS/vuE7ycsJ4kmvczp6b5kZa8y6htxlPa288PPx0Krq81SidbWtYvik+9GdW2FTf8uO54NyQ5M+1qLd4lNesOy2Bj7spLv1K3sJ9uScTkrjPTOSv6xrjOhbV2VPGNBWpsua5Z7mOPS8mtQusDsrltHPwtqi0cWl28TRop41OdxNtGjc4LI1zHVRdpRWv914TM3FWvUq/mNwWUZE7z7JRbbS95BtXMOCxo+3gIu4UW4hV2m93XgZlnU6yW9w159g0batCg8qT4PgM3Mmk2eQq7ySQhd3uH1fzMUaVCW8l2hXtMxafMhbVFjsKLi8fBPPcQZcbV05Z+PuN3Z230+LS9wtWpJx1OUu3uzyvwNa2O52vcppNceCZLZVOTRnW0sw144z3DuztoYws5MfQ1Hby+BGtbZGqN0mi1ajcZIQt13FlN8h/qURqW6xoQZW0bFTi4s4u5snSevzO1uJOLw/mZG2qsdx72NeHeaibZVG60NSyudDDox10G7aLTNfSZPpkdhtRWuXzb/E5+72dipJZy8ZSb0f4cz6FsCvCf8ACbbqxWZPHD4J8Dl7ilvXDjuupiXNOG6u1LOWdMc9Y5b6s+5/lNbrBtqOZbjhuyxlr/GhVebKVNubb7ny/U+m3XVxkqeYb+OC4rv5mHtXYsY1OulrHd1y9FjsXvNYZ7/23o04HZtk5SbTbXZjHyeXktp2iUnBceOO872rslVFGUHurjouJkbR2PGnNVXhxekm0k0+WObNc970mnI1NmYnlyxnklh/Mbp08rCee7idNtPZ0XHM0kuTXHBi2Ozln+FLKzqnp34WNRctzaaYUaGJtQT3ueXp+I9OjLHHd7Xj8jqVs6K1wk3xZkSkqu9GC3sZy3w+Y5bu1kYN1Te8lDdzjk1qu4btaTf3ljOm7z7y2l9H5vEsqL92qx8smxbWazjjJLDZbkSOZX0dw3uaKXHOcjdOhKDjFYaxrj8+46TqSmVst/K+8o8MfqLlv7XTI2lab0Gs7meZn05OKUYrfksrONX79XwOj2jbSdN/Zi5dnL4Z5mXsqTlJRW7mOkm/vY7CZXraVKtS0TxuySzrnHvLa9FuOFo+TNOGysuTlrnguOEIULeVPMm809eOcr5k2pS3tcSy8trT3FVr9jelNvDfB8tfyNue7LEZYTl/LnURvNjtzX2lu4w488e7RovKWGk7W4jL7rTLpU0xfZ+xVTbcW9e3D/Y0OpJKE52kW1LmuD/vieV7bOMYS59oxOSTSfF8Ct729jTdx8R+qx19H3FycG8N53c4/ERrXMYTUJKTcf5m3lfudDHaCTxrnu+BK8tI1M70Yt4441X6mtpf7MNXalP7EeHGWNfmXW9bek1h9mNGu/kU3dLq4qMnlyejWWl8GzR2VsxQ+1nefb/bYvSROcZZxurGPvEZ1MN5Sfdy7y20uXh7+Mp/y66cuBC6ymnCKeePJ/mNLUJXEcJ9vJahTnF55d+hC/jJxzHSXDt4vkUvMY7stdPvLi32aZwJE2aVNPg0e9W0Y1rXTabUoKOnDR683g1pSzhxcWuef00Gl2kof3gg6aCtJ4W5rr2kbmuo43njLLpLR1faeyporncRljVNdmmvfzLqVwm93muRNGyv1JPOYx+HZ8kKz2HTfGK/H9zYqt8ln8izqidryc9V2FHGMyX/ADfuLL6PY4SfxWfyaOn6j4nroCX2cnKPYtThvrHZhkqeyZQ4a9z/AHOq6kPq4sTk5GVWrFNdXLv4/kIpyy3JST7WmvzO9dIHbmeMa5OSt9tJYi8fgM9apfdwmuZ0E9mRfGMX3pEFsiHqr5IxcY1ycdtTabj9l8TBvLw+iX30Wpz4rXvMe7+gUX92TRrGRN7V/Rm9i6Oc5lzR7a3CjPIpQ+iNSnnGvZhim0aNWKw4T054b/Izlhv6WV1yrvRr7vabljepxOP2btrMFDGuMdwtabVdOWGcbhY19u5ubzDRZSvE+DMVX0XHjlsTo1GnhP5EkStfa0k0snD/AElu2pqL+7xWmvYdZc1WoNy7Px5HBbdquUlnkdZN9MNGwjnDRu0LJvng5TZlw1pyOktL8qV9w2BZN7jxFyjxklhy71k1fpEoQUZSUoyTWMJLV6a5WvE/L2yOna8opKMLWWOcoVc/hXS/Afqf/wBFXr40bH/pV/8A9B0k73f8M7fdqmz6ShKVdY+02pSlHea7YtYf/Kz11t+lmnncw0mlmSxp93mfmXafSzc1cucKDecr7NX7Puj/ABsJG5adP95CCgqVnosbzp197v0uMZ+BjhlMp++//GuT7tsCb3PtZbzhKSafwUsZH7q0UsN5010f5n53qdPt60luWunPq6uX3vruZGv083kk06VprzVOsn/viYWXr9XlH23b8qu8owUYw0bqPDXdjkZK2jQjPCWKj03l91vu979x8b9M91uuLp20ovk4VvzVdP8AERvOlGtP/wBG0j7406ifz65stwtnr/pNx+go7KqSlmbhKD5JvPw0PKVvGGYqLil2rR/F8T4ta9Ol3CKiqVphLCzTreeU3fTZdTTUqVph/wDt1v1rjHG2f1f8G4+w3MnFyUsRjj7GMfEXtqjbc1jq8cVx07dD43X6XrmWkoWzWn8lXTHZisi2l0x3KWOrtWvfCr+lZGrDk+xSl10cwljD44/A8s7V/ey03ph8F3HxqfS7cYaVK1jn1YVU/wDeKodLFyo7u7Q792pn/dJcam4+5UZLLSbb55/QXo1t1tuGOx4WfwyfD5dKNzlPFFNf6Z+YTqdKly8ZVHT/AEz8wuqvKPu9Svom3u8BTa7nJYhutPnz/Y+KXPSpcSWHC3+EanmnkOlS5XCNDhj7tTC7v4vEkxsW5R9iVtFtZeKuOPFae409xJpvjjHE+DrpPuN7e3KGcYzu1Ne/+KTr9KdxLGYUNOyNRf8AlHGpyj7ncZzhLTHHTiKVauFjGZ8lpr8j42+lm59Wh4J+aRfSrcZ3urtsrnuVPNNSLyfY4TcYJzy2+SWpOE1nMkk+WcZwfHPS3c+pb+Cp5pT6UrjO84W7a4Nxqaf90nFOT7dWgt3ewtNVx/Qx7C8zPda3VnnnP48T5j6Xbn1LfwVPNKLjpSuJNScLfK4Pdq+aa1uapyfaa1gpYys4Ea9zuTUHonwen79p8qpdLdylhU7fw1fOEa/SNWlJylTt23/pq4Xd/FJq60cn2arsrEt6DxnWSfPuGp2Lx+f9D4xQ6VrmP8tB98J/pUReul659nbeCr5xd39TcfQ7q3kniGYxzhqWVnubHKlnu0+O6scefefKL3pTr1Fuyp2+OOkaqf8AvC66SK+MONFx47rjUx/uZKdPptpWTbi3lY4vm+4XhTam01LHDLTS+fA+bS6Qa2MblBa50jU80vqdJ9y47uKOPdGef9wI+pbuYqPDPDHBNGTeWD4VG32Na/NYOCh0lV0kurt3u8Hu1M/7pZPpQrvjTt33wqeaWUdpb7GlFqW8sf17DqKdBccLKWM41+Z8eq9JNZ4/h26xwxGosf8Ad/MLXpKrxy92i89saj/8pKj69XS0y32+4YprPA+O+k+vnO5Q8NTzT2n0oV1wp2+uv3avmjY+xzi8aY+P9/oeWtXVp8Vz5d3YfIfStcepb+Cp5oLpXuPUt/BU80g+zwjlvK4cGezopNN6JrD10+B8ZXSxcepb+Cp5p5U6VrhrDp27T/01fOGh9no04yb3Wn2Y1enHPxJKhrg+KbO6VLilndp2+vbGq/8AzFy6XLjj1dtn/wCFbzyWUfZZwSzqtOJ7GlwPiy6V7jOXStX3wqecSqdLVw8fwrZY4JRrJfLrxqm32uVNJa6InTpJ6r4Hw/0tXOc7lv3btXd+XXFq6YLn2dt4KvnGeNXb7d1Hceq1PiL6Yrn2dt4KvnHvpiuvZ23gq+cLjV2+z1tlQecxi33Izr76LUan3o/FNp/gz5T6Y7r2dt4KvnAumO69nbeCr5xdU2+kVvoTBZUJzj8c/pkhS+j7jwkm128T51LpiuvUtvBV84rn0vXT/kt13QqfrVGqcndX+zpSwpNrDeEtcilx9G8cU3pxeP2OM9LFxx3LbwVPNPKvStcPjC3/AOnPzC6TbrYbA1wml7sfkxa7sVB4UpN54bv6nLPpOr+zt+/dq/pWFrnpBrS4woeGenc3UbFm0ceAAaAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAf/9k=\n", "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from IPython.display import YouTubeVideo\n", "YouTubeVideo('RODMMX3SOBs', width=700, height=500)" ] }, { "cell_type": "markdown", "id": "3e03f723", "metadata": { "hide_input": true, "slideshow": { "slide_type": "skip" } }, "source": [ "High induced velocities from the rotor can kick up dust/snow on unprepared landing grounds. This leads to degradation in visual cues that are necessary for pilots to maintain aircraft orientation. " ] }, { "cell_type": "markdown", "id": "4c00e863", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "```{admonition} Aside\n", "Dust is even otherwise a major concern due to moving mechanical components involved. Heavy maintainance is warranted when operating in a dusty environment. \n", "\n", "\"... (during the Vietnam war) dust put more helicopters out of commission than enemy action\"\n", "\n", "--- [Military archives](https://www.usace.army.mil/About/History/Historical-Vignettes/Military-Construction-Combat/054-Helicopter-Dust/) \n", "\n", "\"... America’s elite rescue force had lost eight men, seven helicopters, and a C-130, and had not even made contact with the enemy. It was a debacle. It defined the word “debacle.”\"\n", "\n", "--- [article](https://www.theatlantic.com/magazine/archive/2006/05/the-desert-one-debacle/304803/) detailing the failed Iran hostage rescue attempt by American Delta forces (an unexpected dust cloud encountered en-route was one of the leading causes of the failure)\n", "\n", "Its not all *bad* though. Some really pretty pictures can be obtained when sand particles strike rotor blades.\n", "--- [Kopp–Etchells effect](https://en.wikipedia.org/wiki/Kopp%E2%80%93Etchells_effect)\n", "```" ] }, { "cell_type": "markdown", "id": "1d84ad16", "metadata": { "hide_input": true, "slideshow": { "slide_type": "slide" } }, "source": [ "## Mathematical formulation: Climb" ] }, { "cell_type": "markdown", "id": "0643227f", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "
\n", " \n", "

Climbing rotor flowfield model in momentum theory

\n", "
\n", "\n", "\n", "$$\\text{Mass Conservation:} \\quad \\dot{m}=\\rho A(V+v) $$\n", "\n", "$$\\text{Momentum Conservation:} T=\\dot{m} w $$\n", "\n", "$$\\text{Energy balance:} T(V+v)=\\frac{1}{2} \\dot{m}(V+w)^{2}-\\frac{1}{2} \\dot{m} V^{2}$$\n", "\n", "$$=\\frac{1}{2} \\dot{m} w(w+2 V)$$\n", "\n", "\n", "$$\\Rightarrow w=2 v $$\n", "\n", "$$p_{0}+\\frac{1}{2} \\rho(V+w)^{2}-p_{0}+\\frac{1}{2} \\rho V^{2}=T / A$$ \n", "\n", "$$ \\Rightarrow \\frac{1}{2} \\rho w^{2}+\\rho Vw = T / A$$\n", "$$ \\Rightarrow \\rho w(\\frac{w}{2}+V)=T / A $$\n", "\n", "substitute $2v$ and solve quadratic equation in $v$\n", "\n", "\\begin{align} v &=-\\frac{V}{2}+\\sqrt{\\left(\\frac{V}{2}\\right)^{2}+\\frac{T}{2 \\rho A}} \\\\ &=-\\frac{V}{2}+\\sqrt{\\left(\\frac{V}{2}\\right)^{2}+v_{h}^{2}} \\text{ } \\rho A v_h$ so momentum imparted to the air is the same but the velocity is smaller." ] }, { "cell_type": "markdown", "id": "370b381f", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "Total power required in climb = Climb power ($TV$) + induced power ($Tv$)" ] }, { "cell_type": "markdown", "id": "23f34fb6", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "$$P = T (V+v) > P_h = Tv_h$$" ] }, { "cell_type": "markdown", "id": "e1a1699a", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "$$\\frac{P}{T} = \\frac{V}{2}+ \\sqrt{ \\left(\\frac{V}{2} \\right)^2 + v_h^2}$$" ] }, { "cell_type": "markdown", "id": "5454816d", "metadata": { "slideshow": { "slide_type": "subslide" } }, "source": [ "- Other root indicates that flow at the rotor disc is in opposite direction (bottom to top) which is not physically viable\n", "\\begin{align}\n", "v + V &=-\\frac{V}{2}-\\sqrt{\\left(\\frac{V}{2}\\right)^{2}+\\frac{T}{2 \\rho A}} + V = \\frac{V}{2}-\\sqrt{\\left(\\frac{V}{2}\\right)^{2}+\\frac{T}{2 \\rho A}} \\quad < 0\n", "\\end{align}" ] }, { "cell_type": "markdown", "id": "675acdf2", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Further applications " ] }, { "cell_type": "markdown", "id": "748099c8", "metadata": { "cell_style": "split", "slideshow": { "slide_type": "-" } }, "source": [ "- propeller driven aircraft\n", "- VAWT/HAWT wind turbines\n", "- tidal turbines" ] }, { "cell_type": "markdown", "id": "737e8940", "metadata": { "cell_style": "split", "slideshow": { "slide_type": "-" } }, "source": [ "
\n", " \n", "

Turbines of different kinds [source]

\n", "\n", "
" ] }, { "cell_type": "markdown", "id": "94fa3e7e", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Limitations" ] }, { "cell_type": "markdown", "id": "f9d011d7", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "\n", "\n", "Momentum theory only takes into account the rotor disc area. It does not take into account other rotor geometry parameters - blade twist, taper, chord, sweep angle (forward and back), anhedral/dihedral angle, or aerodynamic parameters - airfoil type." ] }, { "cell_type": "markdown", "id": "68d973a3", "metadata": { "cell_style": "center", "slideshow": { "slide_type": "-" } }, "source": [ " - rotor geometry not modelled\n", " - rotor torque not modelled\n", " - i.e. blade drag due to rotational velocity \n", " - no losses modelled\n", " - swirl flow, viscous drag\n", " - Does not include aerodynamic limitations of an airfoil\n", " - only maps thrust to induced velocity" ] }, { "cell_type": "markdown", "id": "20bd996b", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Modelling losses in momentum theory" ] }, { "cell_type": "markdown", "id": "c14915c2", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "The fact that the rotor is not an actuator disk but rather a geomteric entity with finite number of lifting surfaces (blades) leads to physical phenomena that have not been mentioned in the above discussion or the derivation of the results. Consequently, the derived momentum theory results do not include their effects. Turns out these effects are non-negligible and can result in an error of atleast 10-15% in the power calculation of most rotors if not accounted for. So, the way they are accounted for is using an empirical (i.e. based on experimental results) parameter called induced power factor $\\kappa$. It is an adjustment adopted because momentum theory is not perfect. " ] }, { "cell_type": "markdown", "id": "cbb7dad0", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "- tip losses\n", "- blade root cutout\n", "- non-uniform induced velocity\n", "- swirl in rotor wake\n", "\n", "$$P_i=\\kappa Tv_i=\\kappa\\frac{T^{3 / 2}}{\\sqrt{2 \\rho A}}$$\n", "$$\\kappa\\approx1.15$$\n", "\n", "
\n", " \n", " \n", "
\n", "\n", "
\n", "

Wing tip vortex \n", "

\n", "

Root cutout\n", "

\n", "
\n" ] }, { "cell_type": "markdown", "id": "0c01090d", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "While $\\kappa$ helps improve the power prediction using momentum theory, all the shortcomings of momentum theory discussed earlier still persist. For e.g. nothing in the equation above tells us whether a rotor designed for (say) a 5-ton helicopter would also be able to generate 100-tons of thrust! Obviously, the answer is no. The fact that airfoil blades would stall after a certain point and fail to generate any further lift is simply not included in any form within the assumptions or the derivation (go ahead, scroll up and convince yourself that that is the case!)\n", "There are no free lunches. Since the momentum theory is based on approximations introduced to the NSE, we cannot obtain all information about the rotor operation or as quantitatively accurate results as can potentially be obtained using NSE. In other words the momentum theory approximations are appropriate in certain scenarios and as a consequence the theory will lead to pretty reliable results. In a lot of other scenarios where the approximations don't apply or are violated, the results will be *bad*. " ] }, { "cell_type": "markdown", "id": "8a65c966", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "### Right application context" ] }, { "cell_type": "markdown", "id": "b692bed1", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "* Provides *good* estimate in steady conditions - climb, hover, (limited range) descent\n", " - models only the *dominant* physics exhibited by rotors in steady state operation\n", " - used extensively in preliminary design" ] }, { "cell_type": "markdown", "id": "1ec2c6c7", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "* relation between T, P, $v$, DL is desired" ] }, { "cell_type": "markdown", "id": "e9c62f33", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "* easy to see direct cause and effect relations\n", " - \"...momentum approach may provide a clarity of the overall picture that could be lacking when more complicated theories are applied\" {cite}`B28`" ] }, { "cell_type": "markdown", "id": "dbe348bb", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "### Not so right application context" ] }, { "cell_type": "markdown", "id": "e27ad0d8", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "- heavily loaded rotor\n", " - for e.g. tiltrotors" ] }, { "cell_type": "markdown", "id": "b02d8767", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "- sudden changes in blade geometry - twist, taper, sweep etc." ] }, { "cell_type": "markdown", "id": "32be95ef", "metadata": { "slideshow": { "slide_type": "fragment" } }, "source": [ "- Transient conditions - for e.g. maneuvers" ] }, { "cell_type": "markdown", "id": "4a63e2c8", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "A higher fidelity (i.e. better) model is required to overcome the limitations of momentum theory. This is where the blade element momentum theory (BEMT) comes in!" ] } ], "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": true, "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 }