From 53bd3d4aa449e898a37bb02998b0b7f498e53995 Mon Sep 17 00:00:00 2001 From: Miikka Vaisala Date: Fri, 16 Aug 2019 11:13:21 +0800 Subject: [PATCH] Correction factor for the equation. --- acc/mhd_solver/stencil_process.sps | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps index 3f17d92..2103454 100644 --- a/acc/mhd_solver/stencil_process.sps +++ b/acc/mhd_solver/stencil_process.sps @@ -162,9 +162,12 @@ momentum(int3 globalVertexIdx, in VectorField uu, in ScalarField lnrho, in Scala ) + zeta * gradient_of_divergence(uu) #if LSINK - + sink_gravity(globalVertexIdx) - - sink_accretion(globalVertexIdx, lnrho, dt) / exp(value(lnrho)) - * value(uu); //TODO: Confirm + //Gravity term + + sink_gravity(globalVertexIdx) + //Corresponding loss of momentum + - (Scalar(1.0) / Scalar( (dsx*dsy*dsz) * exp(value(lnrho)))) // Correction factor by unit mass + * (sink_accretion(globalVertexIdx, lnrho, dt) * value(uu)) // As in Lee et al.(2014) + ; #else ; #endif