An official website of the United States government
The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($$Re$$)$$2000 \le Re \le 5000$$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$$Re^\gamma$$with the critical exponent$$\gamma \approx -1.25$$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$$\gamma \approx -1.25$$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.
more »
« less
This content will become publicly available on March 1, 2026
Convergence speed and approximation accuracy of numerical MCMC
Abstract When implementing Markov Chain Monte Carlo (MCMC) algorithms, perturbation caused by numerical errors is sometimes inevitable. This paper studies how the perturbation of MCMC affects the convergence speed and approximation accuracy. Our results show that when the original Markov chain converges to stationarity fast enough and the perturbed transition kernel is a good approximation to the original transition kernel, the corresponding perturbed sampler has fast convergence speed and high approximation accuracy as well. Our convergence analysis is conducted under either the Wasserstein metric or the$$\chi^2$$metric, both are widely used in the literature. The results can be extended to obtain non-asymptotic error bounds for MCMC estimators. We demonstrate how to apply our convergence and approximation results to the analysis of specific sampling algorithms, including Random walk Metropolis, Metropolis adjusted Langevin algorithm with perturbed target densities, and parallel tempering Monte Carlo with perturbed densities. Finally, we present some simple numerical examples to verify our theoretical claims.
more »
« less
- Award ID(s):
- 1720433
- PAR ID:
- 10574552
- Publisher / Repository:
- Cambridge University Press on behalf of Applied Probability Trust
- Date Published:
- Journal Name:
- Advances in Applied Probability
- Volume:
- 57
- Issue:
- 1
- ISSN:
- 0001-8678
- Page Range / eLocation ID:
- 101 to 133
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Stochastic embeddings of finite metric spaces into graph-theoretic trees have proven to be a vital tool for constructing approximation algorithms in theoretical computer science. In the present work, we build out some of the basic theory of stochastic embeddings in the infinite setting with an aim toward applications to Lipschitz free space theory. We prove that proper metric spaces stochastically embedding into$$\mathbb {R}$$-trees have Lipschitz free spaces isomorphic to$$L^1$$-spaces. We then undergo a systematic study of stochastic embeddability of Gromov hyperbolic metric spaces into$$\mathbb {R}$$-trees by way of stochastic embeddability of their boundaries into ultrametric spaces. The following are obtained as our main results: (1) every snowflake of a compact, finite Nagata-dimensional metric space stochastically embeds into an ultrametric space and has Lipschitz free space isomorphic to$$\ell ^1$$, (2) the Lipschitz free space over hyperbolicn-space is isomorphic to the Lipschitz free space over Euclideann-space and (3) every infinite, finitely generated hyperbolic group stochastically embeds into an$$\mathbb {R}$$-tree, has Lipschitz free space isomorphic to$$\ell ^1$$, and admits a proper, uniformly Lipschitz affine action on$$\ell ^1$$.more » « less
-
We examine the circular, self-similar expansion of frictional rupture due to fluid injected at a constant rate. Fluid migrates within a thin permeable layer parallel to and containing the fault plane. When the Poisson ratio$$\nu =0$$, self-similarity of the fluid pressure implies fault slip also evolves in an axisymmetric, self-similar manner, reducing the three-dimensional problem for the evolution of fault slip to a single self-similar dimension. The rupture radius grows as$$\lambda \sqrt {4\alpha _{hy} t}$$, where$$t$$is time since the start of injection and$$\alpha _{hy}$$is the hydraulic diffusivity of the pore fluid pressure. The prefactor$$\lambda$$is determined by a single parameter,$$T$$, which depends on the pre-injection stress state and injection conditions. The prefactor has the range$$0\lt \lambda \lt \infty$$, the lower and upper limits of which correspond to marginal pressurisation of the fault and critically stressed conditions, in which the fault-resolved shear stress is close to the pre-injection fault strength. In both limits, we derive solutions for slip by perturbation expansion, to arbitrary order. In the marginally pressurised limit ($$\lambda \rightarrow 0$$), the perturbation is regular and the series expansion is convergent. For the critically stressed limit ($$\lambda \rightarrow \infty$$), the perturbation is singular, contains a boundary layer and an outer solution, and the series is divergent. In this case, we provide a composite solution with uniform convergence over the entire rupture using a matched asymptotic expansion. We provide error estimates of the asymptotic expansions in both limits and demonstrate optimal truncation of the singular perturbation in the critically stressed limit.more » « less
-
Abstract In this paper, we consider random dynamical systems formed by concatenating maps acting on the unit interval$$[0,1]$$in an independent and identically distributed (i.i.d.) fashion. Considered as a stationary Markov process, the random dynamical system possesses a unique stationary measure$$\nu $$. We consider a class of non-square-integrable observables$$\phi $$, mostly of form$$\phi (x)=d(x,x_0)^{-{1}/{\alpha }}$$, where$$x_0$$is a non-recurrent point (in particular a non-periodic point) satisfying some other genericity conditions and, more generally, regularly varying observables with index$$\alpha \in (0,2)$$. The two types of maps we concatenate are a class of piecewise$$C^2$$expanding maps and a class of intermittent maps possessing an indifferent fixed point at the origin. Under conditions on the dynamics and$$\alpha $$, we establish Poisson limit laws, convergence of scaled Birkhoff sums to a stable limit law, and functional stable limit laws in both the annealed and quenched case. The scaling constants for the limit laws for almost every quenched realization are the same as those of the annealed case and determined by$$\nu $$. This is in contrast to the scalings in quenched central limit theorems where the centering constants depend in a critical way upon the realization and are not the same for almost every realization.more » « less
-
We extend the Matsuno–Gill model, originally developed on the equatorial$$\beta$$-plane, to the surface of the sphere. While on the$$\beta$$-plane the non-dimensional model contains a single parameter, the damping rate$$\gamma$$, on a sphere the model contains a second parameter, the rotation rate$$\epsilon ^{1/2}$$(Lamb number). By considering the different combinations of damping and rotation, we are able to characterize the solutions over the$$(\gamma, \epsilon ^{1/2})$$plane. We find that the$$\beta$$-plane approximation is accurate only for fast rotation rates, where gravity waves traverse a fraction of the sphere's diameter in one rotation period. The particular solutions studied by Matsuno and Gill are accurate only for fast rotation and moderate damping rates, where the relaxation time is comparable to the time on which gravity waves traverse the sphere's diameter. Other regions of the parameter space can be described by different approximations, including radiative relaxation, geostrophic, weak temperature gradient and non-rotating approximations. The effect of the additional parameter introduced by the sphere is to alter the eigenmodes of the free system. Thus, unlike the solutions obtained by Matsuno and Gill, where the long-term response to a symmetric forcing consists solely of Kelvin and Rossby waves, the response on the sphere includes other waves as well, depending on the combination of$$\gamma$$and$$\epsilon ^{1/2}$$. The particular solutions studied by Matsuno and Gill apply to Earth's oceans, while the more general$$\beta$$-plane solutions are only somewhat relevant to Earth's troposphere. In Earth's stratosphere, Venus and Titan, only the spherical solutions apply.more » « less
