Regarding the numerical stabilization of two-dimensional advection equation, \begin{equation} \dfrac{\partial f}{\partial t} + \Big(\dfrac{d\varepsilon_1(k)}{dk}\Big)\dfrac{\partial f}{\partial z} - \Big( \dfrac{d\varepsilon_2(z)}{dz} \Big)\dfrac{\partial f}{\partial k} = 0 \end{equation} I have encountered some texts that propose using entropy functions. For example, it is mentioned that the choice of \begin{equation} h(f,z,k) = fe^{\varepsilon_1+\varepsilon_2} \end{equation} will do the job for the above equation.

How is the given function (or any other entropy function) increasing the numerical stability? What is the main idea behind using such functions?

  • 1
    $\begingroup$ What do the papers that "propose using entropy functions" suggest? $\endgroup$ May 4 '16 at 19:49
  • 1
    $\begingroup$ Would you add a reference to the texts that you have read ? $\endgroup$
    – BrunoLevy
    May 5 '16 at 7:12
  • 1
    $\begingroup$ math.la.asu.edu/~chris/ttsp02.pdf This paper, for example. $\endgroup$ May 5 '16 at 10:06
  • $\begingroup$ Well, that's a reference. But what do they propose to use the entropy function for? $\endgroup$ May 5 '16 at 14:55
  • $\begingroup$ For stabilization. When I try to directly discretize the two-dimensional advection equation, I get crazy oscillations and negative values for $f(z,k)$. I don't know how such entropy functions help with stabilizing the numerical system. $\endgroup$ May 6 '16 at 8:25

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.