From b5a9e599de704456db71e80d77949838b9b09cbc Mon Sep 17 00:00:00 2001 From: Billy Cao Date: Sat, 5 Sep 2026 14:09:05 +0000 Subject: [PATCH] Fix LF volume-fraction flux and interface velocity --- src/simulation/m_riemann_solver_lf.fpp | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/src/simulation/m_riemann_solver_lf.fpp b/src/simulation/m_riemann_solver_lf.fpp index 1ea7ffa12..eaee6d6ab 100644 --- a/src/simulation/m_riemann_solver_lf.fpp +++ b/src/simulation/m_riemann_solver_lf.fpp @@ -306,11 +306,12 @@ contains ! Advection flux and source: interface velocity for volume fraction transport $:GPU_LOOP(parallelism='[seq]') do i = eqn_idx%adv%beg, eqn_idx%adv%end - flux_rsx_vf(${SF('')}$, i) = (qL_prim_rsx_vf(${SF('')}$, i) - qR_prim_rsx_vf(${SF(' + 1')}$, & - & i))*s_M*s_P/(s_M - s_P) - flux_src_rsx_vf(${SF('')}$, i) = (s_M*qR_prim_rsx_vf(${SF(' + 1')}$, & - & i) - s_P*qL_prim_rsx_vf(${SF('')}$, i))/(s_M - s_P) + flux_rsx_vf(${SF('')}$, i) = (s_M*qR_prim_rsx_vf(${SF(' + 1')}$, & + & i)*vel_R(norm_dir) - s_P*qL_prim_rsx_vf(${SF('')}$, & + & i)*vel_L(norm_dir) + s_M*s_P*(qL_prim_rsx_vf(${SF('')}$, & + & i) - qR_prim_rsx_vf(${SF(' + 1')}$, i)))/(s_M - s_P) end do + flux_src_rsx_vf(${SF('')}$, eqn_idx%adv%beg) = (s_M*vel_R(norm_dir) - s_P*vel_L(norm_dir))/(s_M - s_P) if (bubbles_euler) then ! From HLLC: Kills mass transport @ bubble gas density @@ -345,7 +346,7 @@ contains ! Geometrical source of the void fraction(s) is zero $:GPU_LOOP(parallelism='[seq]') do i = eqn_idx%adv%beg, eqn_idx%adv%end - flux_gsrc_rsx_vf(${SF('')}$, i) = flux_rsx_vf(${SF('')}$, i) + flux_gsrc_rsx_vf(${SF('')}$, i) = 0._wp end do end if #:endif