From a1c07502ecb7030454c116d13ee63aac913c4bb2 Mon Sep 17 00:00:00 2001 From: Qiao Lei Date: Fri, 13 Feb 2015 14:20:57 -0600 Subject: [PATCH] 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. --- examples/step-33/step-33.cc | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) 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; -- 2.39.5