skip to main content

Title: Stability of a Subcritical Fluid Model for Fair Bandwidth Sharing with General File Size Distributions
This work concerns the asymptotic behavior of solutions to a (strictly) subcritical fluid model for a data communication network, where file sizes are generally distributed and the network operates under a fair bandwidth-sharing policy. Here we consider fair bandwidth-sharing policies that are a slight generalization of the [Formula: see text]-fair policies introduced by Mo and Walrand [Mo J, Walrand J (2000) Fair end-to-end window-based congestion control. IEEE/ACM Trans. Networks 8(5):556–567.]. Since the year 2000, it has been a standing problem to prove stability of the data communications network model of Massoulié and Roberts [Massoulié L, Roberts J (2000) Bandwidth sharing and admission control for elastic traffic. Telecommunication Systems 15(1):185–201.], with general file sizes and operating under fair bandwidth sharing policies, when the offered load is less than capacity (subcritical conditions). A crucial step in an approach to this problem is to prove stability of subcritical fluid model solutions. In 2012, Paganini et al. [Paganini F, Tang A, Ferragut A, Andrew LLH (2012) Network stability under alpha fair bandwidth allocation with general file size distribution. IEEE Trans. Automatic Control 57(3):579–591.] introduced a Lyapunov function for this purpose and gave an argument, assuming that fluid model solutions are sufficiently smooth in time and space that they are strong solutions of a partial differential equation and assuming that no fluid level on any route touches zero before all route levels reach zero. The aim of the current paper is to prove stability of the subcritical fluid model without these strong assumptions. Starting with a slight generalization of the Lyapunov function proposed by Paganini et al., assuming that each component of the initial state of a measure-valued fluid model solution, as well as the file size distributions, have no atoms and have finite first moments, we prove absolute continuity in time of the composition of the Lyapunov function with any subcritical fluid model solution and describe the associated density. We use this to prove that the Lyapunov function composed with such a subcritical fluid model solution converges to zero as time goes to infinity. This implies that each component of the measure-valued fluid model solution converges vaguely on [Formula: see text] to the zero measure as time goes to infinity. Under the further assumption that the file size distributions have finite pth moments for some p > 1 and that each component of the initial state of the fluid model solution has finite pth moment, it is proved that the fluid model solution reaches the measure with all components equal to the zero measure in finite time and that the time to reach this zero state has a uniform bound for all fluid model solutions having a uniform bound on the initial total mass and the pth moment of each component of the initial state. In contrast to the analysis of Paganini et al., we do not need their strong smoothness assumptions on fluid model solutions and we rigorously treat the realistic, but singular situation, where the fluid level on some routes becomes zero, whereas other route levels remain positive.  more » « less
Award ID(s):
1712974 1206772
Author(s) / Creator(s):
Date Published:
Journal Name:
Stochastic Systems
Page Range / eLocation ID:
251 to 273
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We consider a connection-level model proposed by Massoulié and Roberts for bandwidth sharing among file transfer flows in a communication network. We study weighted proportionally fair sharing policies and establish explicit-form bounds on the weighted sum of the expected numbers of flows on different routes in heavy traffic. The bounds are linear in the number of critically loaded links in the network, and they hold for a class of phase-type file-size distributions; that is, the bounds are heavy-traffic insensitive to the distributions in this class. Our approach is Lyapunov drift based, which is different from the widely used diffusion approximation approach. A key technique we develop is to construct a novel inner product in the state space, which then allows us to obtain a multiplicative type of state-space collapse in steady state. Furthermore, this state-space collapse result implies the interchange of limits as a byproduct for the diffusion approximation of the unweighted proportionally fair sharing policy under phase-type file-size distributions, demonstrating the heavy-traffic insensitivity of the stationary distribution. 
    more » « less
  2. We consider the well-known Lieb-Liniger (LL) model for \begin{document}$ N $\end{document} bosons interacting pairwise on the line via the \begin{document}$ \delta $\end{document} potential in the mean-field scaling regime. Assuming suitable asymptotic factorization of the initial wave functions and convergence of the microscopic energy per particle, we show that the time-dependent reduced density matrices of the system converge in trace norm to the pure states given by the solution to the one-dimensional cubic nonlinear Schrödinger equation (NLS) with an explict rate of convergence. In contrast to previous work [3] relying on the formalism of second quantization and coherent states and without an explicit rate, our proof is based on the counting method of Pickl [65,66,67] and Knowles and Pickl [44]. To overcome difficulties stemming from the singularity of the \begin{document}$ \delta $\end{document} potential, we introduce a new short-range approximation argument that exploits the Hölder continuity of the \begin{document}$ N $\end{document}-body wave function in a single particle variable. By further exploiting the \begin{document}$ L^2 $\end{document}-subcritical well-posedness theory for the 1D cubic NLS, we can prove mean-field convergence when the limiting solution to the NLS has finite mass, but only for a very special class of \begin{document}$ N $\end{document}-body initial states.

    more » « less
  3. We study the focusing NLS equation in $R\mathbb{R}^N$ in the mass-supercritical and energy-subcritical (or intercritical ) regime, with $H^1$ data at the mass-energy threshold $\mathcal{ME}(u_0)=\mathcal{ME}(Q)$, where Q is the ground state. Previously, Duyckaerts–Merle studied the behavior of threshold solutions in the $H^1$-critical case, in dimensions $N = 3, 4, 5$, later generalized by Li–Zhang for higher dimensions. In the intercritical case, Duyckaerts–Roudenko studied the threshold problem for the 3d cubic NLS equation. In this paper, we generalize the results of Duyckaerts–Roudenko for any dimension and any power of the nonlinearity for the entire intercritical range. We show the existence of special solutions, $Q^\pm$, besides the standing wave $e^{it}Q$, which exponentially approach the standing wave in the positive time direction, but differ in its behavior for negative time. We classify solutions at the threshold level, showing either blow-up occurs in finite (positive and negative) time, or scattering in both time directions, or the solution is equal to one of the three special solutions above, up to symmetries. Our proof extends to the $H^1$-critical case, thus, giving an alternative proof of the Li–Zhang result and unifying the critical and intercritical cases. These results are obtained by studying the linearized equation around the standing wave and some tailored approximate solutions to the NLS equation. We establish important decay properties of functions associated to the spectrum of the linearized Schrödinger operator, which, in combination with modulational stability and coercivity for the linearized operator on special subspaces, allows us to use a fixed-point argument to show the existence of special solutions. Finally, we prove the uniqueness by studying exponentially decaying solutions to a sequence of linearized equations. 
    more » « less
  4. Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena. 
    more » « less
  5. This work is motivated by a longstanding interest in the long time behavior of flow‐structure interaction (FSI) PDE dynamics. We consider a linearized compressible flow structure interaction (FSI) PDE model with a view of analyzing the stability properties of both the compressible flow and plate solution components. In our earlier work, we gave an answer in the affirmative to question of uniform stability for finite energy solutions of said compressible flow‐structure system, by means of a “frequency domain” approach. However, the frequency domain method of proof in that work is not “robust” (insofar as we can see), when one wishes to study longtime behavior of solutions of compressible flow‐structure PDE models, which track the appearance of the ambient state onto the boundary interface. Nor is a frequency domain approach in this earlier work availing when one wishes to consider the dynamics, in long time, of solutions to physically relevant nonlinear versions of the compressible flow‐structure PDE system under present consideration (e.g., the Navier–Stokes nonlinearity in the PDE flow component or a nonlinearity of Berger/Von Karman type in the plate equation). Accordingly, in the present work, we operate in the time domain by way of obtaining the necessary energy estimates, which culminate in an alternative proof for the uniform stability of finite energy compressible flow‐structure solutions. Since there is a need to avoid steady states in our stability analysis, as a prerequisite result, we also show here that zero is an eigenvalue for the generators of flow‐structure systems, whether the material derivative term be absent or present. Moreover, we provide a clean characterization of the (one dimensional) zero eigenspace, with or without material derivative, under an appropriate assumption on the underlying ambient vector field.

    more » « less