Source author record

Ian A. Frigaard

Ian A. Frigaard appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn cond-mat.soft math.NA math.OC Numerical Analysis

Catalog footprint

What is connected

4works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Fully turbulent flows of viscoplastic fluids in a rectangular duct

Turbulent flows of viscoplastic fluids at high Reynolds numbers have been investigated recently with direct numerical simulations (DNS) but experimental results have been limited. For this reason, we carry out an experimental study of fully turbulent flows of a yield stress fluid in a rectangular aspect ratio channel with a high-resolution laser doppler velocimetry (LDA) setup. We employ aqueous Carbopol solutions, often considered to be a simple yield stress fluid. We formulate different concentrations to address the effect of the rheology of the fluid on the turbulence statistics at an approximately constant Reynolds number. Additionally, we also perform experiments with a single Carbopol formulation at different Reynolds numbers to study its effect. The flow analysis is performed via rheology measurements, turbulence statistics and power spectral densities of velocity fluctuations. The addition of Carbopol to the flow increases turbulence anisotropy, with an enhancement of streamwise velocity fluctuations and a decrease in wall normal velocity fluctuations in comparison to water at the same mean velocity. This change is reflected on the power spectral densities of streamwise velocity fluctuations, where we observe a large increase in energy of large scale turbulent structures. Conversely, the energy of smaller scales is decreased in comparison to water, where the energy drops with a steeper scale than the Newtonian power law of $k_x^{-5/3}$. As we increase the Reynolds number with a Carbopol solution, the streamwise Reynolds stresses approach Newtonian values in the core, which suggests diminishing effects of shear-thinning. The power spectral densities reveal that the energy content at larger scales decreases slightly with the Reynolds number. However, the shear thinning effects do not disappear even as the Reynolds number approaches 50000.

preprint2021arXiv

Flow onset for a single bubble in yield-stress fluids

We use computational methods to determine the minimal yield-stress required in order to hold static a buoyant bubble in a yield-stress liquid. The static limit is governed by the bubble shape, the dimensionless surface tension ($γ$) and the ratio of the yield-stress to the buoyancy stress ($Y$). For a given geometry, bubbles are static for $Y > Y_c$, which we determine for a range of shapes. Given that surface tension is negligible, long prolate bubbles require larger yield-stress to hold static compared to oblate bubbles. Non-zero $γ$ increases $Y_c$ and for large $γ$ the yield-capillary number ($Y/γ$) determines the static boundary. In this limit, although bubble shape is important, bubble orientation is not. 2D planar and axisymmetric bubbles are studied.

preprint2020arXiv

Computing the Yield Limit in Three-dimensional Flows of a Yield Stress Fluid About a Settling Particle

Calculating the yield limit $Y_c$ (the critical ratio of the yield stress to the driving stress), of a viscoplastic fluid flow is a challenging problem, often needing iteration in the rheological parameters to approach this limit, as well as accurate computations that account properly for the yield stress and potentially adaptive meshing. For particle settling flows, in recent years calculating $Y_c$ has been accomplished analytically for many antiplane shear flow configurations and also computationally for many geometries, under either two dimensional (2D) or axisymmetric flow restrictions. Here we approach the problem of 3D particle settling and how to compute the yield limit directly, i.e. without iteratively changing the rheology to approach the yield limit. The presented approach develops tools from optimization theory, taking advantage of the fact that $Y_c$ is defined via a minimization problem. We recast this minimization in terms of primal and dual variational problems, develop the necessary theory and finally implement a basic but workable algorithm. We benchmark results against accurate axisymmetric flow computations for cylinders and ellipsoids, computed using adaptive meshing. We also make comparisons of accuracy in calculating $Y_c$ on comparable fixed meshes. This demonstrates the feasibility and benefits of directly computing $Y_c$ in multiple dimensions. Lastly, we present some sample computations for complex 3D particle shapes.

preprint2010arXiv

On the stability of plane Couette-Poiseuille flow with uniform cross-flow

We present a detailed study of the linear stability of plane Couette-Poiseuille flow in the presence of a cross-flow. The base flow is characterised by the cross flow Reynolds number, $R_{inj}$ and the dimensionless wall velocity, $k$. Squire's transformation may be applied to the linear stability equations and we therefore consider 2D (spanwise-independent) perturbations. Corresponding to each dimensionless wall velocity, $k\in[0,1]$, two ranges of $R_{inj}$ exist where unconditional stability is observed. In the lower range of $R_{inj}$, for modest $k$ we have a stabilisation of long wavelengths leading to a cut-off $R_{inj}$. This lower cut-off results from skewing of the velocity profile away from a Poiseuille profile, shifting of the critical layers and the gradual decrease of energy production. Cross-flow stabilisation and Couette stabilisation appear to act via very similar mechanisms in this range, leading to the potential for robust compensatory design of flow stabilisation using either mechanism. As $R_{inj}$ is increased, we see first destabilisation and then stabilisation at very large $R_{inj}$. The instability is again a long wavelength mechanism. Analysis of the eigenspectrum suggests the cause of instability is due to resonant interactions of Tollmien-Schlichting waves. A linear energy analysis reveals that in this range the Reynolds stress becomes amplified, the critical layer is irrelevant and viscous dissipation is completely dominated by the energy production/negation, which approximately balances at criticality. The stabilisation at very large $R_{inj}$ appears to be due to decay in energy production, which diminishes like $R_{inj}^{-1}$. Our study is limited to two dimensional, spanwise independent perturbations.