--- title: Analytical Branch-Wise Drag Model url: https://www.emergentmind.com/topics/analytical-branch-wise-drag-model type: topic --- # Analytical Branch-Wise Drag Model Searching arXiv for recent and relevant papers on analytical drag models and branch-wise formulations. I’ll look up arXiv entries related to analytical drag models, ensemble drag-based models, and branch-wise decompositions. “Analytical Branch-Wise Drag Model” (*Editor’s term*) denotes an analytical drag formulation in which the governing force law, reduced dynamics, or closure is decomposed into explicit branches associated with distinct sign conditions, thresholds, geometric classes, or scale-separated regimes. Across the literature considered here, branch structure appears in drag-based coronal mass ejection propagation, unsteady wave drag, sparse-data drag-reduction response modeling, landslide drag, rough-surface drag decomposition, canopy drag partition, filtered multiphase drag, solute drag, dusty-disc instability, atmospheric geostrophic drag, and flexible-blade reconfiguration [1801.07473], [2005.00857], [1909.08310], [2207.09337], [2204.12660], [2404.07893], [1811.00714], [1412.0297], [2212.04986], [2309.06650], [2501.03215]. This suggests that the term does not denote a single canonical equation; rather, it identifies a family of analytical constructions in which drag or drag-governed dynamics is written as a finite set of interpretable branches instead of a single undifferentiated empirical law. ## 1. Conceptual structure of branch-wise drag formulations The surveyed literature supports several recurrent branch types. A first class is **sign-explicit branching**, where the drag law changes its analytical expression according to the sign of a relative velocity, as in coronal mass ejection propagation or landslide motion. A second class is **history-dependent branching**, where the force is split into transient and asymptotic steady parts, or into local and nonlocal contributions, as in unsteady wave drag. A third class is **hierarchy- or manifold-based branching**, where a high-dimensional drag response is reduced to variation along a one-dimensional ridge and variation across it, or where unresolved morphology is mapped to an equivalent branch order. A fourth class is **component-wise partitioning**, where total drag is decomposed into viscous and pressure parts, ground and obstacle parts, or additive sub-grid corrections. A fifth class is **stability or regime branching**, where drag couples branches such as one-fluid versus two-fluid behavior, static versus fluttering motion, or conventionally neutral versus stable boundary layers. This suggests that “branch-wise” need not mean only piecewise formulas in one scalar variable. In the cited works, branches may be defined by the sign of $v-w$, by the existence of memory kernels, by parameter-space ridgelines, by hierarchy levels such as $n_{\mathrm{eff}}$, by drag partitions such as $\tau=\tau_R+\tau_S$, or by mode transitions in a dispersion relation. The common analytical feature is explicit regime separation with interpretable transition conditions. ## 2. Sign-explicit kinematic drag laws In heliospheric propagation, the drag-based model (DBM) assumes that beyond about $R \gtrsim 15\,R_{\odot}$ a CME is governed solely by magnetohydrodynamical drag with constant solar-wind speed $w$ and constant drag parameter $\gamma$. The governing law is \[ \frac{dv}{dt}=-\gamma(v-w)|v-w|. \] If $v_0>w$, the fast-CME branch is decelerating, \[ v(t)=w+\frac{v_0-w}{1+\gamma (v_0-w)(t-t_0)}, \] \[ r(t)=r_0 + w(t-t_0) + \frac{1}{\gamma}\ln\!\left[1+\gamma (v_0-w)(t-t_0)\right]. \] If $v_0\Lambda$, it is on a downwind branch. A second branch structure enters through the earth-pressure coefficient \[ K = 2\sec^2\phi\left(1 \mp \sqrt{1-\cos^2\phi\sec^2\delta}\right)-1, \] with the minus sign for expanding motion $(\partial_x u\ge 0)$ and the plus sign for contracting motion $(\partial_x u<0)$. Because \[ \lambda(K)=g^z\left\{\big[(1-\gamma)K+\gamma\big]\alpha_s+(1-\alpha_s)\right\}, \qquad \beta=\frac{\alpha}{\lambda(K)\,h}, \] the drag coefficient is branch-dependent: \[ \beta = \begin{cases} \dfrac{\alpha}{\lambda(K_{\text{exp}})\,h}, & \partial_x u \ge 0,\\[1.0em] \dfrac{\alpha}{\lambda(K_{\text{con}})\,h}, & \partial_x u < 0. \end{cases} \] The paper’s central claim is that drag is fundamentally different for expanding and contracting motions [2207.09337]. ## 3. Memory, asymptotic, and threshold branches For a disturbance moving on an inviscid, irrotational, incompressible fluid of infinite depth, the unsteady surface elevation obeys \[ \partial_t^2 \hat\zeta(\mathbf k,t)+\omega(k)^2\hat\zeta(\mathbf k,t) = -\frac{1}{\rho}k\,\hat P_e(\mathbf k)e^{-i\mathbf k\cdot \mathbf r_0(t)}, \] with dispersion relation \[ \omega(k)^2=gk+\frac{\gamma}{\rho}k^3. \] The exact wake is the history integral \[ \hat\zeta(\mathbf k,t) = -\int_0^t \sin\!\big(\omega(k)(t-\tau)\big)\, \frac{k\hat P_e(\mathbf k)}{\rho\,\omega(k)} e^{-i\mathbf k\cdot \mathbf r_0(\tau)} \,d\tau, \] and the exact unsteady drag is \[ \mathbf R_w(t) = \frac{\mathbf m}{4\pi^2} \iint d^2\mathbf k\; \frac{i\,k\,|\hat P_e(\mathbf k)|^2(\mathbf k\cdot \mathbf m)} {\rho\,\omega(k)} \int_0^t d\tau\, \sin\!\big(\omega(k)(t-\tau)\big) e^{i\mathbf k\cdot(\mathbf r_0(t)-\mathbf r_0(\tau))}. \] For straight-line motion with axisymmetric pressure, this reduces to \[ \mathbf R_w(t) = -\frac{\mathbf m}{2\pi} \int_0^\infty dk\, \frac{k^3|\hat P_e(k)|^2}{\rho\,\omega(k)} \int_0^t d\tau\, \sin\!\big(\omega(k)(t-\tau)\big)\, J_1\!\big(k(r_0(t)-r_0(\tau))\big). \] The resulting branch structure is explicit. If motion reaches constant velocity in finite time, the unsteady wake and drag converge to the steady Havelock solution; this defines a **transient versus asymptotic steady** split. The same formulas also distinguish **local/quasi-static versus nonlocal/history-dependent** drag, because $R_w(t)\approx R_H(v(t))$ is only an approximation to the exact functional of trajectory history. Additional branches arise from the dispersion relation: **pure gravity** versus **capillary-gravity** behavior; **below-radiation-threshold** versus **above-threshold** asymptotics in capillary-gravity flow, with \[ c_{\min}=\left(\frac{4g\gamma}{\rho}\right)^{1/4}; \] and **subcritical-transient** versus **supercritical-transient** decay, separated by \[ c_{\mathrm{crit}}\approx 0.77\,c_{\min}. \] The paper reports that if $v_\inftyc_{\mathrm{crit}}$. For pure gravity waves, quasi-static replacement by steady Havelock drag becomes reasonable when $a\ll g$, but in capillary-gravity cases significant oscillatory memory may persist even for slow ramps [2005.00857]. ## 4. Hierarchy- and manifold-based branch reductions In turbulent boundary-layer drag reduction by spanwise traveling transversal surface waves, the response surface \[ (\lambda^+,T^+,A^+) \mapsto \Delta c_d \] was reduced to a one-dimensional ridgeline \[ (T_r^+(\lambda^+),A_r^+(\lambda^+)) = \arg\max_{T^+,A^+} J(\lambda^+,T^+,A^+), \] with ridge response \[ J_r(\lambda^+) = J(\lambda^+,T_r^+(\lambda^+),A_r^+(\lambda^+)). \] For $\lambda^+\gtrsim 1000$, the fitted ridge coordinates are \[ T_r^+ = 44 - 46721 \exp(-0.0128 \lambda^+), \qquad A_r^+ = 100 - 113 \exp(-0.0009 \lambda^+), \] and the along-ridge drag-reduction law is \[ \widehat{\Delta c}_{d,r} = 0.95\,A_r^+\sqrt{\frac{2\pi}{T_r^+}} - 5.16. \] With normalized cross-ridge coordinates \[ T^*=\frac{T^+}{T_r^+},\qquad A^*=\frac{A^+}{A_r^+}, \] the full self-similar model becomes \[ \widehat{\Delta c}_d(\lambda^+,T^+,A^+)=\widehat{\Delta c}_{d,r}(\lambda^+)\, F(T^*,A^*). \] The model was trained on 71 LES cases over $\lambda^+\in[500,3000]$, $T^+\in[20,120]$, $A^+\in[10,78]$, achieved $R^2=0.93$ for the SVR interpolant and $R^2=0.92$ for the reduced self-similar model, and extrapolated successfully to $\lambda^+=5000$ with relative errors of $2.1\%$ at the on-ridge point $B=(5000,44,99)$ and $0.7\%$ at the off-ridge point $B'=(5000,44,92)$ [1909.08310]. In porous-media tree modeling, the analogous branch variable is not a trajectory branch but an equivalent hierarchy level. The local porous sink is \[ S_{u,i}=-\frac{1}{2}C_D\,\mathrm{PAD}\,|\mathbf{u}|u_i, \] with a cell-wise closure \[ C_D(\mathbf{x}) = C_D\!\left(n_{\mathrm{eff}}(\mathbf{x}),Re_{\mathrm{eff}}(\mathbf{x})\right). \] Here \[ \Omega(\mathbf{x})=\mathrm{SVR}\times R_g, \qquad n_{\mathrm{eff}}(\mathbf{x}) = \mathcal{F}_\Omega\!\left(\Omega(\mathbf{x})\right), \] and \[ Re_{\mathrm{eff}}(\mathbf{x})= \frac{|\mathbf{u}(\mathbf{x})|\,\lambda_R R_g(\mathbf{x})}{\nu}, \qquad \lambda_R=2.77. \] Performance was evaluated with aerodynamic porosity \[ \mathrm{AP} = \frac{\langle u_x \rangle_{\text{leeward plane}}}{\langle u_x \rangle_{\text{windward plane}}}, \] using planes at $x/H=-0.6$ and $x/H=+0.6$. Over grid resolutions $\Delta/H = 1,1/2,1/4,1/6,1/8,1/10$ and $Re_H \approx 2{,}600,6{,}600,66{,}000$, the proposed model reduced the standard deviation of $(1-\mathrm{AP})$ across grids to about $8\%$, compared with about $20\%$ for a general conventional constant-drag model and about $12\%$ for an advanced constant-drag model [2605.25096]. A direct branch-wise force model appears in fractal trees modeled as assemblies of cylindrical segments. For each branch, \[ F_d = C_d \left(\frac{1}{2} \rho U_\infty^2 \right) (L D) \sin^3 \alpha, \qquad F_f = C_f \left(\frac{1}{2} \rho U_\infty^2 \right) (L D) \sin^3 \alpha. \] Whole-tree coefficients are formed from branchwise sums, \[ C_{D\text{-est}} = \frac{F_{D\text{-est}}}{\frac{1}{2}\rho U_\infty^2 A_x}, \qquad C_{F\text{-est}} = \frac{F_{F\text{-est}}}{\frac{1}{2}\rho U_\infty^2 A_x}, \qquad C_{P\text{-est}} = C_{D\text{-est}} - C_{F\text{-est}}. \] Because branch Reynolds number scales with diameter, larger branches enter drag crisis first while smaller branches remain subcritical. For the tested L-system trees, the generation scaling laws are \[ L_k = 0.8\,L_{k-1}, \qquad D_k = 0.5^{0.5} D_{k-1}. \] The analysis indicates a whole-tree drag-crisis transition near $Re_H \approx 3\times10^6$ under uniform inflow and near $Re_H \approx 1.5\times10^5$ when inflow turbulence with $I_u\approx 8\%$ is imposed. Increasing structural complexity smooths the transition because smaller branches remain subcritical [2603.27954]. ## 5. Partition and sub-grid branch decompositions For rough surfaces, the Transpiration–Resistance framework augments slip-length modeling with two constitutive drag-partition parameters: the **shear correction factor** $\mathcal S_c$ and the **pressure correction factor** $\mathcal P_c$. In the flat-interface case, \[ F_v = \mathcal{S}_c\int_{\Gamma}{\tau_{xy} \, dS}, \qquad F_p = \mathcal{P}_c\int_{\Gamma}{\tau_{xy} \, dS}, \] with \[ \mathcal{P}_c+\mathcal{S}_c=1. \] For curved interfaces such as a rough cylinder, \[ F_v = \mathcal{S}_c F^{\mathrm{uncorr}}_v, \qquad F_p = F^{\mathrm{uncorr}}_p+\mathcal{P}_c F^{\mathrm{uncorr}}_v. \] The same microscale cell problem used to compute slip parameters supplies $\mathcal L$, $\mathcal M$, $\mathcal P_c$, and $\mathcal S_c$. In Couette flow, corrected viscous- and pressure-drag errors dropped from order-$100\%$ uncorrected levels to values such as approximately $0.275\%$ and $0.596\%$ for ellipse roughness and $2.04\%$ and $4.63\%$ for square roughness at $h=0.1l$, $\eta=0.1$, over $Re=0.1,3,9,15$ [2204.12660]. For canopy-like roughness, the Raupach partition is \[ \tau=\tau_R+\tau_S=\rho u_*^2. \] An isolated obstacle is characterized by shelter area and shelter volume, \[ A = c_A b h \frac{U_h}{u_*}, \qquad V = c_V b h^2 \frac{U_h}{u_*}, \] and the array density is expressed through \[ \lambda = \frac{n b h}{S}. \] In the infinite-area limit, \[ \tau_S = \tau_{S0}\exp\left(-\frac{\lambda A}{bh}\right), \qquad \tau_R = \lambda \frac{\Phi_0}{bh}\exp\left(-\frac{\lambda V}{bh^2}\right). \] If $V=Ah$ and $\beta=C_R/C_S$, the partition simplifies to \[ \frac{\tau_S}{\tau} = \frac{1}{1+\beta\lambda}, \qquad \frac{\tau_R}{\tau} = \frac{\beta\lambda}{1+\beta\lambda}. \] The revisiting analysis interprets the classical multiplicative sheltering law as a stochastic renormalization rather than literal geometric subtraction of wake area. Mean fitted values differ by roughness type: plants gave approximately $C_R\approx 0.24$ and $c_A\approx 0.19$, whereas cubes gave approximately $C_R\approx 0.53$ and $c_A\approx 0.63$ [2404.07893]. In filtered Eulerian gas–solid drag, the branch structure is additive and sub-grid. A second-order Taylor expansion of the microscopic drag coefficient yields \[ \frac{\bar{\beta}_i}{\beta^*}\approx 1+(i)+(ii)+(vi)+(vii)+\cdots \] and, in the low-Re reduction, \[ \frac{\bar{\beta}_i}{\beta^*}\approx 1 +(a)\frac{\tilde{v}_{d,g,i}}{\tilde{v}_i} +(b)\frac{\tilde{v}_{d,s,i}}{\tilde{v}_i} +(f)\frac{\overline{(\phi'_s)^2}}{(1-\bar{\phi}_s)\bar{\phi}_s} +(f)\frac{\overline{(\phi'_s)^2 v''_i}}{(1-\bar{\phi}_s)\bar{\phi}_s\,\tilde{v}_i}. \] Here the gas drift branch is primarily responsible for drag reduction, the solid drift branch tends to attenuate that reduction, the scalar variance term increases filtered drag over the full solid-volume-fraction range, and the third-order term has the same coefficient as the variance term but mostly decreases drag because it is typically negative. This is one of the clearest examples of a formally additive analytical branch decomposition [1811.00714]. ## 6. Thermodynamic and gravitational branch structures In rapid solidification, the key branching variable is \[ b(V)=\omega+3u\psi_s^2-\frac{\gamma V^2}{M}. \] Using a phase-field-crystal-based amplitude model and a tanh approximation for the order parameter profile, the liquid concentration correction is \[ \psi_1 \approx e^{ - \frac{V}{Mb} \left\{ z + \frac{\delta W}{b} \Phi_o(z) \right\} } \left\{1+ \mathcal{O}\left(\delta/b\right)+\cdots \right\}, \] with \[ \delta = 6 B_2^{\ell}\phi_s^2. \] The segregation coefficient is \[ K(V)=\frac{ \psi_s+1}{ \left(\psi_s+\epsilon \psi_1(W/2)\right) + 1 }. \] If $\gamma=0$, then $b$ never vanishes and complete trapping occurs only asymptotically as $V\to\infty$; this is the diffusive Aziz-like branch. If $\gamma\neq 0$, then $b(V)$ decreases and vanishes at the finite trapping speed \[ V^*=\left\{\frac{M\left(\omega +3u \psi_s^2\right)}{\gamma}\right\}^{1/2}, \] which defines the inertial or Sobolev-like branch. Solute drag is then represented thermodynamically through \[ \Delta G_c = \Delta G_s - f \Delta G_d, \qquad \Delta G_d = (\psi_{\ell} - \psi_s)( \mu_{\ell} - \mu_s), \] with the paper treating $f$ as velocity-independent for a given material setting and assigning the velocity dependence primarily to $\Delta G_d(V)$ [1412.0297]. In dusty self-gravitating discs, drag acts as a branch-coupling parameter between gas and dust. The control parameters are \[ \epsilon = \frac{\Sigma_d}{\Sigma_g}, \qquad \xi = \left(\frac{c_d}{c_g}\right)^2, \qquad \mathrm{St}=t_s\kappa. \] The drag-free two-fluid system has a quartic dispersion relation, \[ \omega^4 -\omega^2(\alpha_g+\alpha_d) + (\alpha_g\alpha_d -\beta_g\beta_d) = 0, \] where \[ \alpha_i = \kappa^2 + c_i^2 k^2 - 2\pi G \Sigma_i |k|, \qquad \beta_i = 2\pi G \Sigma_i |k|. \] With drag but without backreaction, the problem becomes fifth order, \[ D_{\mathrm{nbr}}(\omega,k)=0, \] and the branch interpretation is one-fluid versus two-fluid behavior with an additional drag-relaxation mode. In the weak-drag limit $\mathrm{St}\gg 1$, the model recovers the two-fluid threshold; in the strong-drag limit $\mathrm{St}\ll 1$, it approaches one-fluid behavior. The gas-to-dust transition is fitted by \[ \epsilon_{\mathrm{tr}} = f(\mathrm{St})\sqrt{\xi}, \qquad f(\mathrm{St}) = 1 + 0.72\,\mathrm{St}^{-1.36}, \] and, using \[ \xi = \frac{\alpha_{SS}}{1+\mathrm{St}}, \] the critical Stokes number is approximated as \[ \mathrm{St}_{\mathrm{crit}} \simeq \epsilon^{-2}\alpha_{SS}. \] The main physical consequence is a dust-driven short-wavelength branch with Jeans masses much smaller than in standard gas-only GI; the paper states that the predicted mass is about $10$ Earth masses in a representative regime [2212.04986]. ## 7. Boundary-layer and compliant-structure branches In atmospheric boundary-layer modeling, the drag law is branch-wise in both vertical structure and stability class. The total stress magnitude is taken as \[ T=\left(1-\frac{\hat\xi}{\hat h}\right)^{3/2}, \] and the spanwise stress is written piecewise: \[ \hat T_{yz}= \begin{cases} - g(\hat\xi)\left(1-\dfrac{\hat\xi}{\hat h}\right)^{3/2}, & \hat\xi>\hat\xi_m,\\[1.0em] \dfrac{\hat\xi}{\kappa}\ln\dfrac{\hat\xi}{\hat\xi_0} -(\hat\xi-\hat\xi_0)\Biggl[\dfrac{1}{\kappa}+\dfrac{U_g}{u_*} -\dfrac{1}{2}(5\mu+0.3\mu_N)(\hat\xi-\hat\xi_0)\Biggr], & \hat\xi<\hat\xi_m, \end{cases} \] with \[ g(\hat\xi)=M\left(1-e^{-\Gamma \hat\xi/\hat h}\right), \qquad M=1.43,\ \Gamma=1.2,\ c_m=0.20,\ \hat\xi_m=c_m\hat h. \] The inner branch is MOST-consistent, \[ \frac{U^{\mathrm{inner}}}{u_*} = \frac{1}{\kappa}\ln\frac{\hat\xi}{\hat\xi_0} +(5\mu+0.3\mu_N)(\hat\xi-\hat\xi_0), \] with $\mu=0$ in the conventionally neutral branch and $\mu>0$ in the stable branch. Matching the inner and outer solutions yields a self-consistent geostrophic drag law, \[ \frac{\kappa G\cos\alpha_0}{u_*}=\ln(Ro)-A, \qquad \frac{\kappa G\sin\alpha_0}{u_*}=-B, \qquad u_*=\frac{\kappa G}{\sqrt{[\ln(Ro)-A]^2+B^2}}. \] The paper validates the model against 41 LES cases and reports excellent agreement for CNBL and good overall agreement for SBL, with the main discrepancy near the SBL top where LES retains non-zero stress above $h$ [2309.06650]. For side-by-side flexible blades, branch structure separates a **static reconfiguration regime** from a **flutter regime**. The governing nondimensional parameters are \[ \upbeta = \frac{m_a}{\mu + m_a}, \qquad \mathrm{Ca} = \frac{1}{2}\frac{\rho C_D b U^2 l^3}{EI}, \qquad \mathrm{B} = \frac{(\rho_s - \rho)g b d l^3}{EI}, \qquad \uplambda = \frac{2}{\pi}\frac{C_D}{C_M}\frac{l}{b}, \] and \[ \mathrm{u}^2=\frac{\mathrm{Ca}}{\uplambda}=\frac{m_a U^2 l}{EI}. \] For a bunch of $N$ overlapping identical blades, the equivalent reduction is \[ d_{\mathrm{eq}}=Nd, \qquad (EI)_{\mathrm{eq}}=N EI. \] In the static branch, drag reduction is described by the reconfiguration number \[ \mathcal R=\frac{F}{F_{\mathrm{rigid}}}=\frac{l_e}{l}, \qquad \frac{l_e}{l}=\int_0^1 \cos^3\theta\,d\tilde s. \] The experiments show that $C_{D,\mathrm{bulk}}$ decreases in the static regime starting at \[ \frac{\mathrm{Ca}}{\mathrm{B}} > O(1), \] and settles to an almost constant value in the flutter regime at high $\mathrm{Ca}$. Increasing $\upbeta$, $\mathrm B$, or $\uplambda$ stabilizes the system and delays flutter. The paper’s further claim is that the reactive-force framework, which includes \[ \bm q_{am} = -m_a\left[ \frac{\partial (U_n \bm n)}{\partial t} -\frac{\partial (U_\tau U_n \bm n)}{\partial s} +\frac12\frac{\partial (U_n^2 \bm \tau)}{\partial s} \right], \] is superior to traditional Morison inertia because the convective reactive terms are necessary to capture both flutter onset and flutter-limited drag reduction [2501.03215]. Taken together, these works indicate that branch-wise analytical drag modeling is best understood as a general strategy of explicit regime separation. A common misconception is that such models are merely piecewise curve fits. The surveyed literature instead shows branch definitions derived from sign structure, dispersion relations, wake memory, hierarchy mappings, self-similar ridges, component partitions, and stability thresholds. Another common misconception is that branches are only kinematic; several of the cited models show branches generated by thermodynamic coupling, unresolved morphology, or sub-grid covariance structure. The main limitation, also repeated across domains, is that branchwise tractability is usually purchased by strong assumptions—constant background states, linear theory, independent-branch superposition, equivalent-element reductions, or precomputed constitutive tables—so the utility of any given branch model depends on whether those assumptions remain valid in the target regime. Source: https://www.emergentmind.com/topics/analytical-branch-wise-drag-model