Computational Aerodynamics Questions & Answers  
Question by Student 201527110
Professor, I don't understand how $\Gamma,(\frac{\partial \tau}{\partial t})$ can be chosen. You set $\Gamma=1$ because $\Gamma$ is constant when we trying to get $\Phi_1=0$ but, actually I don't understand how it really can be like this.
03.16.17
This is a good question. First note that $\Gamma$ must be chosen such that it respects the dependencies we specified during the derivation. By definition, $$ \Gamma\equiv\frac{\partial \tau}{\partial t} $$ Thus, because $\tau=\tau(t)$, $\Gamma$ can not depend on $x$ or $y$. $\Gamma$ could depend on $t$ thus. We here choose it to be simply $\Gamma=1$ because that is the simplest expression that satisfies the above dependency.. I'll give you 2 points bonus boost.
03.17.17
Question by Student 201127151
Professor, I have a question about $\Omega$. You said that ${\Omega}={x_{\xi}}{y_{\eta}} + y_{\xi}x_{\eta}$ and it means the area of a cell. I understood that $\Omega$ represents area because its dimension is $[m^2]$. But I don't understand that it is the area of a cell. Is it just the definition, or can we demonstrate it?
Good question. I have pondered this as well when I first encountered the generalized coordinates.. The definition of $\Omega$ is simply ${x_{\xi}}{y_{\eta}} + y_{\xi}x_{\eta}$, not the cell area. But it is very close to the cell area thus. In the next assignment, I will ask you to compare the cell area with $\Omega$ and assess the differences. I'll give you 2 points bonus.
Question by Student 201427142
Professor, I still cannot understand about $\Gamma$ in quastion by student #201527110 $$ \Phi_1 = U \frac{\partial}{\partial\tau}\left(\Omega\Gamma\right) = U \frac{\partial}{\partial\tau}\left(\frac{\partial x}{\partial\xi}\frac{\partial y}{\partial\eta}\frac{\partial\tau}{\partial t}-\frac{\partial y}{\partial\xi}\frac{\partial x}{\partial\eta}\frac{\partial\tau}{\partial t} \right) $$ at here $$\frac{\partial}{\partial\tau}\left(\frac{\partial x}{\partial\xi}\frac{\partial y}{\partial\eta}\frac{\partial\tau}{\partial t}\right) = \left(\frac{\partial}{\partial\tau}\left(\frac{\partial x}{\partial\xi}\right)\right)\frac{\partial y}{\partial\eta}\frac{\partial\tau}{\partial t}+\frac{\partial x}{\partial\xi}\left(\frac{\partial}{\partial\tau}\left(\frac{\partial y}{\partial\eta}\right)\right)\frac{\partial\tau}{\partial t}+\frac{\partial x}{\partial\xi}\frac{\partial y}{\partial\eta}\left(\frac{\partial}{\partial\tau}\left(\frac{\partial\tau}{\partial t}\right)\right)$$ and I understand 1st and 2nd in right term. But I think that the 3rd term($\left(\frac{\partial}{\partial\tau}\left(\frac{\partial\tau}{\partial t}\right)\right)$) cannot be omitted. Because $\frac{\partial\tau}{\partial t}$ is still depending on $\tau$. How can we assume $\Gamma = 1$??
We can set $\Gamma$ to any value as long as it respects the dependency we originally specified. That is, $\Gamma=\partial\tau/\partial t$ should be such that $\tau=\tau(t)$. Thus, we can choose $\Gamma$ to be $1$ or $t$ or $t^2$ but not $x$ or $tx$, etc. Here, it is perfectly fine to choose $\Gamma=1$. One advantage of choosing $\Gamma=1$ is that $$ \frac{\partial}{\partial\tau}\left(\frac{\partial\tau}{\partial t}\right) = \frac{\partial}{\partial\tau}\left(\Gamma\right) = \frac{\partial}{\partial\tau}\left(1\right) = 0 $$ If $\Gamma$ would have been chosen as $t$ or $t^2$, then the latter wouldn't be 0, and this would lead to $\Phi_1\ne 0$, and this would prevent us to write the Euler equations in generalized coordinates in conservative form. I'll give you 2 points bonus boost.
Question by Student 201427564
.PNG
Professor, When I try to open warp file through 'Pluma', this error occured. How should I do?
03.26.17
Warp is a binary file (executable), so it can not be opened in pluma. With pluma you can open text files like .wrp or .txt or .c etc. Please don't do assignment #4 now: the tgz package is not complete yet. I'll explain what to do in assignment 4 tomorrow.
Question by Student 201427564
Professor, I wonder about difference between 'EE' and 'NO'. 'E' makes mesh equally space and 'NO' also. Can we adjust distance with 'E' like below? Block(is,js,ie,je, EE,0.5,2,2, NO); Also, is it right that 'NO' only makes mesh equally 1? (like default) I know this is very basic question but I want to know clearly. Thank you.
03.27.17
Hmm, I am not sure what your first question is.. As for the second, 'NO' does not make the mesh equally 1, it simply needs to be used when the block is 1D (when there are no segments along a certain dimension). I'll give you 1 point bonus boost.
Question by Student 201427102
1.jpg  ./download/file.php?id=3370&sid=698d0a70f64fddfdc488392c22b14538  ./download/file.php?id=3370&t=1&sid=698d0a70f64fddfdc488392c22b14538
for middle block, you use "BLOCK(    ,GG,   );"
2.jpg  ./download/file.php?id=3368&sid=698d0a70f64fddfdc488392c22b14538  ./download/file.php?id=3368&t=1&sid=698d0a70f64fddfdc488392c22b14538
for middle block, you use "JOINT(    ,GG,   );"
3.jpg  ./download/file.php?id=3369&sid=698d0a70f64fddfdc488392c22b14538  ./download/file.php?id=3369&t=1&sid=698d0a70f64fddfdc488392c22b14538
However, this case, for middle block, you use "JOINT(   ,EE,   );" Why I shoud use 'EE' rather than 'GG'??
First, the command is Join(), not Joint(), and is case sensitive. So, if you use JOIN or join it will not work: it has to be spelled Join. For the third mesh, you can not use Join(is,js,ie,je,j,GG..) because the GG segment grabs the spacing from the mesh along the j coordinate. So in this case, he will set the first grid spacing (between j=js and j=js+1) to the distance between the point (i,js-1) and the point (i,js). However, there is no point (i,js-1) because it is outside of the domain. Good question: 2 points bonus boost.
1,  2 ,  3  ...  14    Next  •  PDF 1✕1 2✕1 2✕2  •  New Question
$\pi$