\documentclass[a4paper,12pt,twoside]{article} \usepackage[latin1]{inputenc} \usepackage{ngerman} \usepackage{graphicx} \usepackage{subfigure} \usepackage{a4wide} \usepackage{esvect} \usepackage{pstricks} \usepackage{pst-plot} \usepackage{pstricks-add} \usepackage{multicol,epsf} \usepackage{ae,aecompl} \usepackage{helvet} \usepackage[latin1]{inputenc} \usepackage{amsmath} \usepackage{amssymb} \usepackage{bm} \usepackage{bbm} \usepackage[mathcal]{euscript} \usepackage[T1]{fontenc} \usepackage{hyperref} \definecolor{linkcolor}{rgb}{0,0,0} \hypersetup{colorlinks=true, breaklinks=true, linkcolor=linkcolor, menucolor=linkcolor, urlcolor=linkcolor, citecolor=linkcolor, bookmarksopen=true, bookmarksnumbered=false} \numberwithin{equation}{section} \numberwithin{figure}{section} \pagestyle{headings} \begin{document} \newcommand{\diff}[1]{ \, \mathrm{d} #1 \, } \newcommand{\re}{ R \! \hspace{1mm} \! e } \newcommand{\st}{ S \! \hspace{1mm} \! t } %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Titelseite % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \normalsize \rule{0mm}{0,4cm} \thispagestyle{empty} \begin{figure}[htbp] \centering \includegraphics[scale=0.4]{DLR.png} \label{DLR} \end{figure} \normalsize \begin{center} \Huge Versuch 9 \\[5mm] \LARGE {\bf Bestimmung des Auftriebs, der Zirkulation und des Widerstandes für das Tragflächenprofil Gö 818} \end{center} \normalsize \vspace{5mm} \centerline{\epsfysize=2cm \hfill {\bf Strömungsmechanisches Praktikum} \hfill \epsfxsize=2cm} \centerline{\epsfysize=2cm \hfill des Deutschen Zentrums für Luft- und Raumfahrt \hfill \epsfxsize=2cm} \centerline{\epsfysize=2cm \hfill Georg-August-Universität Göttingen \hfill \epsfxsize=2cm} \vspace{10mm} \begin{center} 15. August 2008 \end{center} \vspace{7mm} \begin{center} \begin{tabular}{rl} \rule{5cm}{0mm} & \rule{5,5cm}{0mm} \\ {\it Praktikanten } \ \ & {\bf Johannes Dörr} \\ & \href{mailto:mail@johannesdoerr.de} {mail@johannesdoerr.de} \\ & \href{http://physik.johannesdoerr.de} {physik.johannesdoerr.de} \vspace{2mm} \\ & {\bf Jan Schumann-Bischoff} \\ & \href{mailto:jansb.stud@googlemail.com} {jansb.stud@googlemail.com} \\ \rule{5cm}{0mm} & \rule{5,5cm}{0mm} \\ {\it Durchführung } \ \ & am 02.07.2008 \\ \rule{5cm}{0mm} & \rule{5,5cm}{0mm} \\ {\it Betreuer } \ \ & Christoph Wolf \\ \rule{5cm}{0mm} & \rule{5,5cm}{0mm} \\ \end{tabular} \end{center} \newpage %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Inhaltsverzeichnis % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \tableofcontents \newpage %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % (Dokument) % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%% \section{Einleitung} %%%%%%%%%%%%%%%%%%%% Dieser Versuch befasst sich mit dem Phänomen des Auftriebs von Tragflächen. Dazu wird die Druckverteilung an einem Modellflügel gemessen und ausgewertet. %%%%%%%%%%%%%%%%% \section{Theorie} %%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%% \subsection{Druckmessung} %%%%%%%%%%%%%%%%%%%%%%%%% Der von uns untersuchte Tragflügel wird von Luft mit der Anströmgeschwindigkeit $u_\infty$ umströmt. Der Flügel enthält 36 Bohrungen an der Oberfläche (siehe Abb. \ref{fig:fluegel_profil}). \begin{figure}[htbp] \centering \includegraphics[width=0.80\textwidth]{Bilder/fluegel_profil.PNG} \caption{Tragflügelprofil mit nummerierten Bohrungen an der Oberfläche zur Druckmessung} \label{fig:fluegel_profil} \end{figure} Von jeder Bohrung geht ein Schlauch auf eine Seite eines mit Wasser gefüllten U-Rohrmanometers. Für $u_\infty=0$ gibt es keine Druckdifferenz zwischen beiden Seiten des U-Rohres, also ist das Manometer nicht aus der Ruhelage ausgeschlagen ($\Delta h=0$). An der Spitze des Flügels gibt es einen Punkt, den Staupunkt, an dem die Strömung zur Ruhe kommt, $u=0$. Nehmen wir die Luftströmung als stationär und inkompressibel an, so kann die Umströmung anhand der Bernoulligleichung \begin{equation} p+\frac{1}{2}\rho u^2=\mathrm{const} \label{eq:bernoulli} \end{equation} beschrieben werden. Dabei hat die Konstante für alle Punkte auf einer Stromlinie den gleichen Wert. In Abb. \ref{fig:fluegel_druck} ist ein Flügel dargestellt. Es ist zu erkennen, dass eine Stromlinie praktisch auf der Oberfläche (bzw. kurz darüber) entlang läuft. \begin{figure}[htbp] \centering \includegraphics[width=0.60\textwidth]{Bilder/fluegel_druck.PNG} \caption{Tragflügel mit eingezeichneten Stromlinien} \label{fig:fluegel_druck} \end{figure} Es gilt also für den Druck $p$ und der Geschwindigkeit $u$ auf einer Stromlinie, welche durch den Staupunkt verläuft, und dem Druck $p_{\mathrm{ges}}$ im Staupunkt (Geschwindigkeit im Staupunkt ist Null) \begin{equation} p_{\mathrm{ges}}=\underbrace{p}_{p_{\mathrm{stat}}}+\underbrace{\frac{1}{2}\rho_{\mathrm{Wasser}} u^2}_{p_{\mathrm{dyn}}} \ . \label{eq:bernoulli_einfach} \end{equation} Diese Drücke lassen sich wie im Folgenden beschrieben mit dem U-Rohrmanometer bestimmen: \\ \\ \textbf{Statischer Druck:}\\ Dieser Druck kann direkt aus dem Ausschlag der Manometer berechnet werden. Die eine Seite des U-Rohres ist mit einer Wandbohrung des Flügels verbunden (siehe Abb. \ref{fig:statischer_druck}). \begin{figure}[htbp] \centering \includegraphics[width=0.40\textwidth]{Bilder/statischer_druck.PNG} \caption{Bestimmung des statischen Drucks mittels U-Rohtmanometer} \label{fig:statischer_druck} \end{figure} Auf der entsprechenden Wasseroberfläche wirkt also der statische Druck $p_{\mathrm{stat}}$ des Punktes auf der Flügeloberfläche und der vernachlässigbare Druck der Luftsäule. Auf der anderen Seite wirkt der Referenzdruck $p_{\mathrm{ref}}$ (Druck für $u_\infty =0$) und der Druck der entstandenen Wassersäule, also \begin{equation} p_{\mathrm{stat}}=p_{\mathrm{ref}}+\rho_{\mathrm{Wasser}}g\Delta h \label{eq:stat_druck} \end{equation} \textbf{Gesamtdruck:} \\ Messen wir den Statischen Druck am Staupunkt, so verschwindet dort die Geschwindigkeit und somit ist $p_{\mathrm{dyn}}=0$. Nach \eqref{eq:bernoulli_einfach} enstpricht hier der statische Druck dem Gesamtdruck. \\ \\ \textbf{Dynamischer Druck:} \\ Nach Bernoulli entspricht der Druck am Staupunkt ($p_{\mathrm{ges}}$) genau der Konstanten der Stromlinie. Ist dieser bekannt, so kann durch die Messung des statischen Drucks $p_{\mathrm{stat}}$ nach \eqref{eq:bernoulli_einfach} der dynamische Druck $p_{\mathrm{dyn}}$ berechnet werden. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \subsection{Auftriebs- und Widerstandsbeiwert} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Man unterscheidet zwei Arten des Auftriebs: Der {\it statische Auftrieb} kommt durch das "`Archimedische Prinzip"' zu Stande, das besagt, dass die Auftriebskraft, die ein Körper in einem Medium erfährt, genauso groß ist wie die Gewichtskraft des von ihm verdrängten Mediums. Die Kraft, die ein Flugzeug zum Fliegen bringt, ist hingegen der {\it dynamische Auftrieb}. Er entsteht dadurch, dass der Körper, in unserem Fall eine Tragfläche, von dem Medium umströmt wird. Sie errechnet sich aus: \begin{align} F_\mathrm{A} = \frac{1}{2} \cdot \rho \cdot c_\mathrm{A} \cdot A \cdot u_\infty^2 \ . \label{eq:Auftriebskraft} \end{align} Dabei ist $\rho$ die Dichte des Mediums, $A$ die Querschnittsfläche des Körpers quer zur Strömungsrichtung, $u_\infty$ die Strömungsgeschwindigkeit und $c_\mathrm{A}$ der (einheitenlose) {\it Auftriebsbeiwert} der Tragfläche, gegeben durch: \begin{align} c_\mathrm{A} = \oint \frac{p_\mathrm{stat} - p_\mathrm{ref}}{p_\mathrm{dyn}} \diff{\!\left(\frac{x}{l}\right)} \ . \label{sec:Auftriebsbeiwert} \end{align} Dabei ist $l$ die Länge der Tragfläche und $x$ der Ort in Richtung der Profilsehne. Die Auftriebskraft wirkt senkrecht zur Strömungsrichtung. Den Integrand nennt man {\it Druckbeiwert}. Eine weitere Kraft, die auf die Tragfläche wirkt, diesmal jedoch der Strömungsrichtung entgegengesetzt, ist die {\it Widerstandskraft}, die sich ergibt aus: \begin{align*} F_\mathrm{W} = \frac{1}{2} \cdot \rho \cdot c_\mathrm{W} \cdot A \cdot u_\infty^2 \end{align*} mit dem (einheitenlosen) {\it Widerstandsbeiwert}: \begin{align*} c_\mathrm{W} = 2 \int \sqrt{\frac{p_\mathrm{ges}^1 - p_\mathrm{stat}^1}{p_\mathrm{dyn}^\infty}} \left(1-\sqrt{\frac{p_\mathrm{ges}^1 - p_\mathrm{stat}^\infty}{p_\mathrm{dyn}^\infty}} \right) \diff{\!\left(\frac{y}{l}\right)} \ . \end{align*} Dabei markiert $\infty$ die Größen (Drücke), die vor der Tragfläche gemessen werden, und $1$ die Werte dahinter. %%%%%%%%%%%%%%%%%%%%%%%% \subsection{Zirkulation} %%%%%%%%%%%%%%%%%%%%%%%% Die Zirkulation enthält Informationen über die Rotation in einem System. Sie ist definiert als: \begin{align*} \Gamma = \oint_C \mathbf{v} \cdot \diff{\mathbf{s}} = \int_A \mathrm{rot} (\mathbf{v}) \cdot \diff{\mathbf{A}} \ . \end{align*} Nach dem {\it Satz von Kutta-Joukowski} errechnet sich die Auftriebskraft aus der Zikulation über: \begin{align*} F_\mathrm{A} = \rho \cdot u_\infty \cdot b \cdot \Gamma \ , \end{align*} wobei $b$ die Breite der Tragfläche ist. Mit \eqref{eq:Auftriebskraft} ergibt sich daraus: \begin{align} \Gamma = \frac{1}{2} \cdot c_\mathrm{A} \cdot b \cdot u_\infty \ . \label{eq:Zirkulation} \end{align} Nach dem {\it Satz von Thomson} ist die Zirkulation in einem System örtlich und zeitlich erhalten. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \subsection{Strömung an der Tragfläche} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Bewegt sich die Tragfläche langsam in dem Medium, so fließt letzteres auf der Unterseite am hinteren Ende um die Kante auf die Oberseite, wo es das oben herum anfließende Medium im sogenannten {\it hinteren Staupunkt} trifft (siehe Abb. \ref{fig:tragfl_zirkulation}). \begin{figure}[htbp] \centering \includegraphics[width=0.60\textwidth]{Bilder/tragflzirkulation.jpg} \caption{Wirbel bei der langsam umströmten Tragfläche} \label{fig:tragfl_zirkulation} \end{figure} Wird die Geschwindigkeit größer, so lösen sich die auf diese Weise entstehenden Wirbel (Anfahrwirbel) ab, da das Medium zu träge ist, um um die untere Kante zu fließen. Auf Grund der Erhaltung der Zirkulation, die vor dem Flügel ja Null ist, muss sich nun ein gegenläufiger Wirbel bilden. Dieser entsteht um die Tragfläche (Abb. \ref{fig:tragfl_zirkulation2}). Er ist jedoch nicht so stark, dass er die Anströmung kompensieren könnte. Er sorgt allerdings dafür, dass die Geschwindigkeit oberhalb der Tragfläche höher ist als unterhalb, was nach dem Gesetz von Bernulli oben zu einem niedrigeren Druck führt. Es resultiert daraus die Auftriebskraft nach oben. \begin{figure}[htbp] \centering \includegraphics[width=0.4\textwidth]{Bilder/tragflzirkulation2.jpg} \caption{Strömung an der Tragfläche} \label{fig:tragfl_zirkulation2} \end{figure} %%%%%%%%%%%%%%%%%%%%%% \section{Durchführung} %%%%%%%%%%%%%%%%%%%%%% Das Tragflächenmodell (Abb. \ref{fig:fluegel_profil}) ist in einem Windkanal positioniert. Für verschiedene Anstellwinkel wird bei gleicher Anströmungsgeschwindigkeit die Druckverteilung mit Hilfe der U-Rohr-Manometer gemessen. \numberwithin{figure}{subsection} %%%%%%%%%%%%%%%%%%%% \section{Auswertung} %%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \subsection{Druckverteilung auf der Flügeloberfläche und Druckbeiwert} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Aus den gemessenen Höhen $\Delta h$ der Wassersäulen lassen sich nach \eqref{eq:stat_druck} die statischen Drücke für jede entsprechende Bohrung auf der Flügeloberfläche berechnen. I.A. sind die ermittelten Werte für $\Delta h$ unterhalb des Flügels positiv und oberhalb negativ. Dies würde einen größeren Druck oberhalb als unterhalb des Flügels bedeuten, was physikalisch aber nicht zu erwarten ist. Aus Konsistenzgründen haben wir deshalb die Vorzeichen von $\Delta h$ für alle Werte umgedreht (vermutlich sind die Schläuche an den U-Rohrmanometern genau anders hermum angeschlossen als von uns beschrieben). Mit den berechneten statischen Drücken lässt sich für jeden Wert mit \begin{displaymath} c_{\mathrm{p}}=\frac{p_{\mathrm{stat}}-p_{\mathrm{ref}}}{p_{\mathrm{dyn}}} \end{displaymath} der Druckbeiwert $c_{\mathrm{p}}$ berechnen. Für alle Anstellwinkel sind in Abb. \ref{fig:cp_xl} $c_{\mathrm{p}}$ über der dimensionslosen Länge $x/l$ des Flügels aufgetragen. \begin{figure}[htbp] \subfigure[-5° Anstellwinkel]{ \includegraphics[width=0.50\textwidth]{Bilder/Graph_-5.pdf}} \subfigure[0° Anstellwinkel]{ \includegraphics[width=0.50\textwidth]{Bilder/Graph_0.pdf}} \subfigure[5° Anstellwinkel]{ \includegraphics[width=0.50\textwidth]{Bilder/Graph_5.pdf}} \subfigure[10° Anstellwinkel]{ \includegraphics[width=0.50\textwidth]{Bilder/Graph_10.pdf}} \subfigure[12.5° Anstellwinkel]{ \includegraphics[width=0.50\textwidth]{Bilder/Graph_125.pdf}} \subfigure[15° Anstellwinkel]{ \includegraphics[width=0.50\textwidth]{Bilder/Graph_15.pdf}} \caption{$c_p$ in Abhängigkeit von $x/l$ für verschiedene Anstellwinkel} \label{fig:cp_xl} \end{figure} $l=40\mathrm{cm}$ ist die tatsächliche Länge des Flügels und $x$ der tatsächliche Abstand einer Bohrung vom linken Ende des Flügels. $x$ kann für jede Bohrung aus Tab. 1 aus dem Versuchsskript (Versuch 9) abgelesen werden. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \subsection{Auftriebsbeiwert und Zirkulation} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Der Auftriebsbeiwert berechnet sich mit \eqref{sec:Auftriebsbeiwert}. Die Integration realisieren wir über eine Aufaddierung der im vorigen Abschnitt berechneten Druckbeiwerte multipliziert mit $\diff{\!\left(\frac{x}{l}\right)} \equiv \Delta \!\left(\frac{x}{l}\right)$. Dabei berechnen wir \begin{align*} \Delta_n \!\left(\frac{x}{l}\right) = \frac{x_{n+1} - x_{n}}{2l} + \frac{x_{n} - x_{n-1}}{2l} = \frac{x_{n+1} - x_{n-1}}{2l} \ , \end{align*} wobei wir die Löcher jetzt entgegen des Uhrzeigersinns um die gesamte Tragfläche herum indiziert haben. Die Zirkulation ergibt sich dann mit $u_\infty=31.1\frac{\mathrm{m}}{\mathrm{s}}$ und $b=280\mathrm{mm}$ aus \eqref{eq:Zirkulation}. Die Werte für die verschiedenen Anstellwinkel sind in der folgenden Tabelle dargestellt: \begin{center} \begin{tabular}{lll} {\bf Anstellwinkel} & {\bf Auftriebsbeiwert} & {\bf Zirkulation} \\ \hline -5° & -0.21 & -0.89 \\ 0° & 0.04 & 0.19 \\ 5° & 0.44 & 1.90 \\ 10° & 0.61 & 2.64 \\ 12.5° & 0.66 & 2.85 \\ 15° & 1.72 & 7.62 \\ \end{tabular} \end{center} Dass der Auftriebsbeiwert bei -5° negativ ist, deckt sich mit den Graphen im vorigen Abschnitt, bei denen sichtbar wird, dass die Druckbeiwerte bei diesem Winkel auf der Oberseite größer sind als unten. Bei den übrigen Winkeln ist es anders herum. Dies zeigt, dass bei -5° ein Abtrieb vorhanden ist, bei den anderen hingegen ein Auftrieb. \begin{figure}[htbp] \centering \includegraphics[width=0.8\textwidth]{Bilder/graph_winkel_cp.pdf} \caption{Auftriebsbeiwert in Abhängigkeit vom Anstellwinkel} \label{fig:auftriebsbeiwert} \end{figure} %%%%%%%%%%%%%%%%%%%% \section{Diskussion} %%%%%%%%%%%%%%%%%%%% Die Versuchsdurchführung verlief problemlos und zügig. Bei der Auswertung haben wir eigentlich erwartet, bei 12.5° den maximalen Auftriebsbeiwert zu finden. Dies bestätigt sich jedoch nicht. \end{document}