From: Qiao Lei Date: Fri, 13 Feb 2015 20:20:57 +0000 (-0600) Subject: Fix bug in energy source term according to Praveen C's report. X-Git-Tag: v8.3.0-rc1~464^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F546%2Fhead;p=dealii.git Fix bug in energy source term according to Praveen C's report. The energy source term due to gravity should be rho*g*v=g*momentum. --- diff --git a/examples/step-33/step-33.cc b/examples/step-33/step-33.cc index c3162bed51..6ae98203de 100644 --- a/examples/step-33/step-33.cc +++ b/examples/step-33/step-33.cc @@ -314,9 +314,7 @@ namespace Step33 forcing[c] = gravity * W[density_component]; break; case energy_component: - forcing[c] = gravity * - W[density_component] * - W[first_momentum_component+dim-1]; + forcing[c] = gravity * W[first_momentum_component+dim-1]; break; default: forcing[c] = 0;