diff --git a/src/incflo_apply_predictor.cpp b/src/incflo_apply_predictor.cpp index f353cadf..34f3b71e 100644 --- a/src/incflo_apply_predictor.cpp +++ b/src/incflo_apply_predictor.cpp @@ -11,7 +11,10 @@ using namespace amrex; // else // conv_u = - del dot (rho u u) // conv_r = - div( u rho ) -// conv_t = - div( u trac ) +// if (m_iconserv_tracer) then +// conv_t = - div( u trac ) +// else +// conv_t = - u dot grad trac // eta_old = visosity at m_cur_time // if (m_diff_type == DiffusionType::Explicit) // divtau _old = div( eta ( (grad u) + (grad u)^T ) ) / rho^n @@ -168,7 +171,6 @@ void incflo::ApplyPredictor (bool incremental_projection) // Compute the explicit advective terms R_u^(n+1/2), R_s^(n+1/2) and R_t^(n+1/2) // if (advection_type == "MOL" ) // Compute the explicit advective terms R_u^n , R_s^n and R_t^n - // Note that "get_conv_tracer_old" returns div(rho u tracer) // ************************************************************************************* compute_convective_term(get_conv_velocity_old(), get_conv_density_old(), get_conv_tracer_old(), get_velocity_old_const(), get_density_old_const(), get_tracer_old_const(),