2
$\begingroup$

The interested equation is advection-diffusion equation. One of the canonical example is Navier-Stokes equations. However, I would like to let the coefficient of diffusion constant goes to zero, $\epsilon \rightarrow 0$, while mesh sizes, $ \Delta x, \Delta dy \rightarrow 0 $, and time steps,$\Delta t \rightarrow 0$, decreasing towards zero. Then, I can obtain inviscid compressible equation which can be an example for the Euler equation of gas dynamics. The following equation is scalar advection-diffusion equation in 2D cartesian coordinates.
$$ \frac{df}{dt} + U\frac{df}{dx} + V\frac{df}{dy} = \epsilon \left( \frac{d^2f}{dx^2} + \frac{d^2f}{dy^2} \right)$$

The main interest is to find exact formulation for the Euler equations of gas dynamics in cylindrical coordinates $(r, \theta)$ with artificial diffusion part.

Edit : The main equation that I am interested in the Euler equations of gas dynamics in 2D. The cartesian vector form can be stated as $$\begin{pmatrix} \rho \\ \rho u \\ \rho v \\ E \end{pmatrix}_t + \begin{pmatrix} \rho u \\ \rho u^2 + p \\ \rho u v \\ u(E + p) \end{pmatrix}_x + \begin{pmatrix} \rho v\\ \rho u v \\ \rho v^2 + p \\ v(E + p) \end{pmatrix}_y = \textbf{0}$$

Here, conservation of mass, momentums in -x and -y directions and energy equations can be seen in a matrix form. However, I would like to add artificial viscosity with small constant factor as in Navier-Stokes equations to smoothen the possible oscillations that will occur during numerical simulation. For this reason, I need to know the vector Laplacian form in cylindrical coordinates for Navier-Stokes equations. Then I can apply the above theory that was stated.

$\endgroup$
7
  • 2
    $\begingroup$ Welcome to Scicomp! I think your question is on-topic, but a little hard to understand. Are you really interested in exactly the equation you give, but transformed into polar coordinates? If so, what are $U$ and $V$? Are they given? If not, can you edit your question to show the exact equation (in Cartesian coordinates) you are interested in? $\endgroup$ Oct 14, 2014 at 14:39
  • $\begingroup$ @ChristianClason Thank you for your interest. I tried to edit my question to be more specific. I hope it makes the question more clear. In my example U and V are velocities which depend on spatial, time and solution variables. $\endgroup$
    – Loading...
    Oct 15, 2014 at 14:16
  • $\begingroup$ Yes, it's clearer now. The vector Laplacian in polar coordinates (I assume that's what you mean; cylindrical coordinates usually refer to $(r,\theta,z)$ in 3D) is actually easy to find online, e.g., mathworld.wolfram.com/VectorLaplacian.html or en.wikipedia.org/wiki/… $\endgroup$ Oct 15, 2014 at 14:40
  • $\begingroup$ You should uncheck my answer. Since it is probably not complete anymore. $\endgroup$
    – Jan
    Oct 15, 2014 at 15:28
  • $\begingroup$ And, it is a particularity in cartesian coordinates and for divergence free vector fields (+ further assumptions) that the diffusive term in the Navier Stokes equations can be written as a Laplacian. You probably will have to use the more general formulation given in the references in my answer. $\endgroup$
    – Jan
    Oct 15, 2014 at 15:38

1 Answer 1

1
$\begingroup$

Have a look at Chapter 3 of the bachelor thesis "Modellierung und Simulation von Dispersionen in turbulenten Strömungen" by Manuel Baumann.

It is in German but the equations and pictures can be understood.

In Chapter 3.1.1, the coordinate transformation that in particular affects the convection is explained.

In Chapter 3.2.1, the transformed RANS equations are formulated. Set the right terms to zero and you get the Euler equations (equation (19) on page 29).

For the Euler equation, the main work will be the convective term that in polar coordinates writes as

enter image description here

Note that in polar coordinates, the coordinate unit vectors $e_\theta$, $e_r$ depend on the angular coordinate $\theta$. This means, e.g., $$ \partial _\theta (U_\theta \cdot e_\theta) = \partial _\theta (U_\theta) \cdot e_\theta - U_\theta \cdot e_r $$ what makes the additional terms appear in the last equality of the equation above.

As a reference, I can recommend "Mechanics of Fluid Flow" by Longwell.

Disclaimer: I was supervising the thesis mentioned above.

$\endgroup$
3
  • $\begingroup$ It will help me to obtain detailed formulation for the Euler equations of gas dynamics in cylindrical coordinate. Thank you. $\endgroup$
    – Loading...
    Oct 13, 2014 at 18:32
  • $\begingroup$ @Jan I've edited your post to add a link to the online version of the thesis and the exact equation therein so that people don't have to hunt for it -- I hope that's OK, otherwise feel free to revert. Still, your answer would be much more useful if you included the actual equation in it (since such PDFs can vanish from the net without warning). $\endgroup$ Oct 14, 2014 at 14:33
  • $\begingroup$ Thanks for the edit, @Christian. I don't have the sources for copy and paste, but I will include a screenshot tomorrow. $\endgroup$
    – Jan
    Oct 14, 2014 at 20:44

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.