Trim verbose paper references and remove duplicate comments

Author: sbryngelson

src/common/m_helper.fpp

@@ -174,15 +174,15 @@ contains
         integer                 :: ir
         real(wp), dimension(nb) :: chi_vw0, cp_m0, k_m0, rho_m0, x_vw, omegaN, rhol0
         real(wp), parameter     :: k_poly = 1._wp !< polytropic index used to compute isothermal natural frequency
-        ! Chapman-Enskog transport coefficients for vapor-gas mixture, Ando JAS (2010) Eq. 2.22
+        ! Chapman-Enskog transport coefficients for vapor-gas mixture, Ando JAS (2010)
 
         phi_vg = (1._wp + sqrt(mu_v/mu_g)*(M_g/M_v)**(0.25_wp))**2/(sqrt(8._wp)*sqrt(1._wp + M_v/M_g))
         phi_gv = (1._wp + sqrt(mu_g/mu_v)*(M_v/M_g)**(0.25_wp))**2/(sqrt(8._wp)*sqrt(1._wp + M_g/M_v))
 
         ! Initial internal bubble pressure (Euler number + Laplace pressure)
         pb0 = Eu + 2._wp/Web/R0
 
-        ! Vapor mass fraction at bubble wall, Ando JAS (2010) Eq. 2.19
+        ! Vapor mass fraction at bubble wall, Ando JAS (2010)
         chi_vw0 = 1._wp/(1._wp + R_v/R_g*(pb0/pv - 1._wp))
 
         ! Mixture specific heat from mass-weighted vapor/gas contributions
@@ -206,7 +206,7 @@ contains
         ! Peclet numbers
         Pe_T(:) = rho_m0*cp_m0(:)/k_m0(:)
 
-        ! Bubble natural frequency, Ando JAS (2010) Eq. B.1
+        ! Bubble natural frequency, Ando JAS (2010)
         omegaN(:) = sqrt(3._wp*k_poly*Ca + 2._wp*(3._wp*k_poly - 1._wp)/(Web*R0))/R0/sqrt(rho0ref)
         do ir = 1, nb
             call s_transcoeff(omegaN(ir)*R0(ir), Pe_T(ir)*R0(ir), Re_trans_T(ir), Im_trans_T(ir))

src/common/m_phase_change.fpp

@@ -54,9 +54,7 @@ contains
     !! selecting the phase change module that will be used (pT- or pTg-equilibrium)
     impure subroutine s_initialize_phasechange_module
 
-        ! Saturation curve coefficients derived from Clausius-Clapeyron relation via stiffened gas EOS. Saurel et al. JCP (2008), Le
-        ! Metayer et al. JFE (2004) gs_min = gamma-1, cvs = specific heat at constant volume, qvps = reference entropy, qvs =
-        ! reference energy
+        ! Saturation curve coefficients via stiffened gas EOS. Saurel et al. JCP (2008), Le Metayer et al. JFE (2004)
         A = (gs_min(lp)*cvs(lp) - gs_min(vp)*cvs(vp) + qvps(vp) - qvps(lp))/((gs_min(vp) - 1.0_wp)*cvs(vp))
 
         B = (qvs(lp) - qvs(vp))/((gs_min(vp) - 1.0_wp)*cvs(vp))

src/common/m_variables_conversion.fpp

@@ -572,8 +572,7 @@ contains
 #endif
                     end if
 
-                    ! Recover primitive variables from relativistic MHD conserved variables via Newton-Raphson iteration on total
-                    ! enthalpy W. Mignone & Bodo A&A (2006) W = total enthalpy, Ga = Lorentz factor, B2 = |B|^2, m2 = |m|^2, S = m.B
+                    ! Relativistic MHD primitive variable recovery, Mignone & Bodo A&A (2006)
                     if (relativity) then
                         if (n == 0) then
                             B(1) = Bx0

src/simulation/m_acoustic_src.fpp

@@ -4,8 +4,7 @@
 
 #:include 'macros.fpp'
 
-!> @brief One-way acoustic source injection following Maeda and Colonius, JCP 2017. Supports planar, focused transducer, transducer
-!! array, and broadband waveforms.
+!> @brief One-way acoustic source injection, Maeda and Colonius JCP (2017)
 module m_acoustic_src
 
     use m_derived_types        !< Definitions of the derived types

src/simulation/m_bubbles.fpp

@@ -75,8 +75,7 @@ contains
 
     end function f_rddot
 
-    !> Bubble wall pressure: stiffened gas with Laplace pressure and viscous stress Ca = cavitation number, Web = Weber number,
-    !! Re_inv = inverse Reynolds number Rayleigh PRSLA (1917), Plesset JAM (1949), Keller-Miksis JASA (1980)
+    !> Bubble wall pressure: stiffened gas with Laplace pressure and viscous stress
     !! @param fR0 Equilibrium bubble radius
     !! @param fR Current bubble radius
     !! @param fV Current bubble velocity
@@ -234,7 +233,7 @@ contains
 
     end function f_rddot_G
 
-    !> Function that computes the bubble wall pressure for Keller--Miksis bubbles
+    !> Keller-Miksis bubble wall pressure
     !! @param fR0 Equilibrium bubble radius
     !! @param fR Current bubble radius
     !! @param fV Current bubble velocity

src/simulation/m_ibm.fpp

@@ -111,9 +111,7 @@ contains
         @:ALLOCATE(ghost_points(1:max_num_gps))
 
         $:GPU_ENTER_DATA(copyin='[ghost_points]')
-        ! Ghost-cell immersed boundary method procedure: 1. Identify ghost points adjacent to IB surface 2. Apply levelset to
-        ! determine inside/outside 3. Compute image points (reflection across IB surface) 4. Interpolate flow variables at image
-        ! points Tseng & Ferziger JCP (2003), Mittal & Iaccarino ARFM (2005)
+        ! Ghost-cell IBM, Tseng & Ferziger JCP (2003), Mittal & Iaccarino ARFM (2005)
         call s_find_ghost_points(ghost_points)
         call s_apply_levelset(ghost_points, num_gps)
 

src/simulation/m_riemann_solvers.fpp

@@ -526,7 +526,7 @@ contains
                                     s_L = min(vel_L(dir_idx(1)) - c_fast%L, vel_R(dir_idx(1)) - c_fast%R)
                                     s_R = max(vel_R(dir_idx(1)) + c_fast%R, vel_L(dir_idx(1)) + c_fast%L)
                                 else if (hypoelasticity) then
-                                    ! Elastic longitudinal wave speed: sqrt(c^2 + (4G/3 + tau_e)/rho), Rodriguez et al. JCP (2019)
+                                    ! Elastic wave speed, Rodriguez et al. JCP (2019)
                                     s_L = min(vel_L(dir_idx(1)) - sqrt(c_L*c_L + (((4._wp*G_L)/3._wp) + tau_e_L(dir_idx_tau(1))) &
                                               & /rho_L), &
                                               & vel_R(dir_idx(1)) - sqrt(c_R*c_R + (((4._wp*G_R)/3._wp) + tau_e_R(dir_idx_tau(1))) &
@@ -2026,8 +2026,7 @@ contains
                                 ! COMPUTING THE DIRECT WAVE SPEEDS
                                 if (wave_speeds == 1) then
                                     if (elasticity) then
-                                        ! Elastic longitudinal wave speed: sqrt(c^2 + (4G/3 + tau_e)/rho), Rodriguez et al. JCP
-                                        ! (2019)
+                                        ! Elastic wave speed, Rodriguez et al. JCP (2019)
                                         s_L = min(vel_L(dir_idx(1)) - sqrt(c_L*c_L + (((4._wp*G_L)/3._wp) + tau_e_L(dir_idx_tau(1) &
                                                   & ))/rho_L), &
                                                   & vel_R(dir_idx(1)) - sqrt(c_R*c_R + (((4._wp*G_R)/3._wp) &
@@ -3101,8 +3100,7 @@ contains
 
                                 if (wave_speeds == 1) then
                                     if (elasticity) then
-                                        ! Elastic longitudinal wave speed: sqrt(c^2 + (4G/3 + tau_e)/rho), Rodriguez et al. JCP
-                                        ! (2019)
+                                        ! Elastic wave speed, Rodriguez et al. JCP (2019)
                                         s_L = min(vel_L(dir_idx(1)) - sqrt(c_L*c_L + (((4._wp*G_L)/3._wp) + tau_e_L(dir_idx_tau(1) &
                                                   & ))/rho_L), &
                                                   & vel_R(dir_idx(1)) - sqrt(c_R*c_R + (((4._wp*G_R)/3._wp) &

src/simulation/m_surface_tension.fpp

@@ -160,7 +160,6 @@ contains
                                                 & j, l, i)
                                 end do
 
-                                ! Continuum surface force capillary stress, Schmidmayer et al. JCP (2017)
                                 flux_src_vf(E_idx)%sf(j, k, l) = flux_src_vf(E_idx)%sf(j, k, &
                                             & l) + sigma*c_divs(num_dims + 1)%sf(j, k, l)*vSrc_rsy_vf(k, j, l, 2)
                             end if
@@ -203,7 +202,6 @@ contains
                                                 & k, j, i)
                                 end do
 
-                                ! Continuum surface force capillary stress, Schmidmayer et al. JCP (2017)
                                 flux_src_vf(E_idx)%sf(j, k, l) = flux_src_vf(E_idx)%sf(j, k, &
                                             & l) + sigma*c_divs(num_dims + 1)%sf(j, k, l)*vSrc_rsz_vf(l, k, j, 3)
                             end if