scholarly journals Lyapunov dimension of elastic turbulence

2017 ◽  
Vol 822 ◽  
Author(s):  
Emmanuel Lance Christopher VI M. Plan ◽  
Anupam Gupta ◽  
Dario Vincenzi ◽  
John D. Gibbon
Keyword(s):  
Reynolds Number ◽  
Constitutive Model ◽  
Numerical Simulations ◽  
Mathematical Analysis ◽  
Critical Value ◽  
Weissenberg Number ◽  
Direct Numerical Simulations ◽  
Chaotic Behaviour ◽  
Two Dimensional ◽  
Lyapunov Dimension

Low-Reynolds-number polymer solutions exhibit a chaotic behaviour known as ‘elastic turbulence’ when the Weissenberg number exceeds a critical value. The two-dimensional Oldroyd-B model is the simplest constitutive model that reproduces this phenomenon. To make a practical estimate of the resolution scale of the dynamics, one requires the assumption that an attractor of the Oldroyd-B model exists; numerical simulations show that the quantities on which this assumption is based are bounded. We estimate the Lyapunov dimension of this assumed attractor as a function of the Weissenberg number by combining a mathematical analysis of the model with direct numerical simulations.

The dynamics of a laminar flow in a symmetric channel with a sudden expansion

2001 ◽  
Vol 436 ◽  
pp. 283-320 ◽  
Author(s):  
T. HAWA ◽  
Z. RUSAK
Keyword(s):  
Reynolds Number ◽  
Laminar Flow ◽  
Numerical Simulations ◽  
Direct Numerical Simulations ◽  
Sudden Expansion ◽  
Asymptotically Stable ◽  
Two Dimensional ◽  
Stable Mode ◽  
Symmetric States ◽  
Asymmetric Equilibrium

Bifurcation analysis, linear stability study, and direct numerical simulations of the dynamics of a two-dimensional, incompressible, and laminar flow in a symmetric long channel with a sudden expansion with right angles and with an expansion ratio D/d (d is the width of the channel inlet section and D is the width of the outlet section) are presented. The bifurcation analysis of the steady flow equations concentrates on the flow states around a critical Reynolds number Rec(D/d) where asymmetric states appear in addition to the basic symmetric states when Re [ges ] Rec(D/d). The bifurcation of asymmetric states at Rec has a pitchfork nature and the asymmetric perturbation grows like √Re − Rec(D/d). The stability analysis is based on the linearized equations of motion for the evolution of infinitesimal two-dimensional disturbances imposed on the steady symmetric as well as asymmetric states. A neutrally stable asymmetric mode of disturbance exists at Rec(D/d) for both the symmetric and the asymmetric equilibrium states. Using asymptotic methods, it is demonstrated that when Re < Rec(D/d) the symmetric states have an asymptotically stable mode of disturbance. However, when Re > Rec(D/d), the symmetric states are unstable to this mode of asymmetric disturbance. It is also shown that when Re > Rec(D/d) the asymmetric states have an asymptotically stable mode of disturbance. The direct numerical simulations are guided by the theoretical approach. In order to improve the numerical simulations, a matching with the asymptotic solution of Moffatt (1964) in the regions around the expansion corners is also included. The dynamics of both small- and large-amplitude disturbances in the flow is described and the transition from symmetric to asymmetric states is demonstrated. The simulations clarify the relationship between the linear stability results and the time-asymptotic behaviour of the flow. The current analyses provide a theoretical foundation for previous experimental and numerical results and shed more light on the transition from symmetric to asymmetric states of a viscous flow in an expanding channel. It is an evolution from a symmetric state, which loses its stability when the Reynolds number of the incoming flow is above Rec(D/d), to a stable asymmetric equilibrium state. The loss of stability is a result of the interaction between the effects of viscous dissipation, the downstream convection of perturbations by the base symmetric flow, and the upstream convection induced by two-dimensional asymmetric disturbances.

Direct numerical simulations of bubbly flows Part 2. Moderate Reynolds number arrays

1999 ◽  
Vol 385 ◽  
pp. 325-358 ◽  
Author(s):  
ASGHAR ESMAEELI ◽  
GRÉTAR TRYGGVASON
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Volume Fraction ◽  
Three Dimensional ◽  
Direct Numerical Simulations ◽  
Dimensional System ◽  
Two Dimensional ◽  
Reynolds Stresses ◽  
Bubble Distribution ◽  
Moderate Reynolds Number

Direct numerical simulations of the motion of two- and three-dimensional finite Reynolds number buoyant bubbles in a periodic domain are presented. The full Navier–Stokes equations are solved by a finite difference/front tracking method that allows a fully deformable interface between the bubbles and the ambient fluid and the inclusion of surface tension. The rise Reynolds numbers are around 20–30 for the lowest volume fraction, but decrease as the volume fraction is increased. The rise of a regular array of bubbles, where the relative positions of the bubbles are fixed, is compared with the evolution of a freely evolving array. Generally, the freely evolving array rises slower than the regular one, in contrast to what has been found earlier for low Reynolds number arrays. The structure of the bubble distribution is examined and it is found that while the three-dimensional bubbles show a tendency to line up horizontally, the two-dimensional bubbles are nearly randomly distributed. The effect of the number of bubbles in each period is examined for the two-dimensional system and it is found that although the rise Reynolds number is nearly independent of the number of bubbles, the velocity fluctuations in the liquid (the Reynolds stresses) increase with the size of the system. While some aspects of the fully three-dimensional flows, such as the reduction in the rise velocity, are predicted by results for two-dimensional bubbles, the structure of the bubble distribution is not. The magnitude of the Reynolds stresses is also greatly over-predicted by the two-dimensional results.

Direct numerical simulations of bubbly flows. Part 1. Low Reynolds number arrays

1998 ◽  
Vol 377 ◽  
pp. 313-345 ◽  
Author(s):  
ASGHAR ESMAEELI ◽  
GRÉTAR TRYGGVASON
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Stokes Flow ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Direct Numerical Simulations ◽  
Two Dimensional ◽  
Rise Velocity ◽  
Regular Array ◽  
Average Rise

Direct numerical simulations of the motion of two- and three-dimensional buoyant bubbles in periodic domains are presented. The full Navier–Stokes equations are solved by a finite difference/front tracking method that allows a fully deformable interface between the bubbles and the ambient fluid and the inclusion of surface tension. The governing parameters are selected such that the average rise Reynolds number is O(1) and deformations of the bubbles are small. The rise velocity of a regular array of three-dimensional bubbles at different volume fractions agrees relatively well with the prediction of Sangani (1988) for Stokes flow. A regular array of two- and three-dimensional bubbles, however, is an unstable configuration and the breakup, and the subsequent bubble–bubble interactions take place by ‘drafting, kissing, and tumbling’. A comparison between a finite Reynolds number two-dimensional simulation with sixteen bubbles and a Stokes flow simulation shows that the finite Reynolds number array breaks up much faster. It is found that a freely evolving array of two-dimensional bubbles rises faster than a regular array and simulations with different numbers of two-dimensional bubbles (1–49) show that the rise velocity increases slowly with the size of the system. Computations of four and eight three-dimensional bubbles per period also show a slight increase in the average rise velocity compared to a regular array. The difference between two- and three-dimensional bubbles is discussed.

Numerical Simulations of Electric Field Driven Self-Assembly of Monolayers of Mixtures of Nanoparticles

2017 ◽  
Author(s):  
E. Amah ◽  
N. Musunuri ◽  
Ian S. Fischer ◽  
Pushpendra Singh
Keyword(s):  
Electric Field ◽  
Physical Properties ◽  
Numerical Simulations ◽  
Self Assembly ◽  
Composite Particle ◽  
Critical Value ◽  
Direct Numerical Simulations ◽  
Composite Particles ◽  
Liquid Interfaces ◽  
Induced Dipole

We numerically study the process of self-assembly of particle mixtures on fluid-liquid interfaces when an electric field is applied in the direction normal to the interface. The force law for the dependence of the electric field induced dipole-dipole and capillary forces on the distance between the particles and their physical properties obtained in an earlier study by performing direct numerical simulations is used for conducting simulations. The inter-particle forces cause mixtures of nanoparticles to self-assemble into molecular-like hierarchical arrangements consisting of composite particles which are organized in a pattern. However, there is a critical electric intensity value below which particles move under the influence of Brownian forces and do not self-assemble. Above the critical value, when the particles sizes differed by a factor of two or more, the composite particle has a larger particle at its core and several smaller particles forming a ring around it. Approximately same sized particles, when their concentrations are approximately equal, form binary particles or chains (analogous to polymeric molecules) in which positively and negatively polarized particles alternate, but when their concentrations differ the particles whose concentration is larger form rings around the particles with smaller concentration.

scholarly journals Waves and vortices in the inverse cascade regime of stratified turbulence with or without rotation

2016 ◽  
Vol 806 ◽  
pp. 165-204 ◽  
Author(s):  
Corentin Herbert ◽  
Raffaele Marino ◽  
Duane Rosenberg ◽  
Annick Pouquet
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Direct Numerical Simulations ◽  
Vertical Shear ◽  
Energy Spectra ◽  
Stratified Turbulence ◽  
Viscous Effects ◽  
Inverse Cascade ◽  
Main Effects ◽  
The Waves

We study the partition of energy between waves and vortices in stratified turbulence, with or without rotation, for a variety of parameters, focusing on the behaviour of the waves and vortices in the inverse cascade of energy towards the large scales. To this end, we use direct numerical simulations in a cubic box at a Reynolds number $Re\approx 1000$, with the ratio between the Brunt–Väisälä frequency $N$ and the inertial frequency $f$ varying from $1/4$ to 20, together with a purely stratified run. The Froude number, measuring the strength of the stratification, varies within the range $0.02\leqslant Fr\leqslant 0.32$. We find that the inverse cascade is dominated by the slow quasi-geostrophic modes. Their energy spectra and fluxes exhibit characteristics of an inverse cascade, even though their energy is not conserved. Surprisingly, the slow vortices still dominate when the ratio $N/f$ increases, also in the stratified case, although less and less so. However, when $N/f$ increases, the inverse cascade of the slow modes becomes weaker and weaker, and it vanishes in the purely stratified case. We discuss how the disappearance of the inverse cascade of energy with increasing $N/f$ can be interpreted in terms of the waves and vortices, and identify the main effects that can explain this transition based on both inviscid invariants arguments and viscous effects due to vertical shear.

scholarly journals Peristaltic transport of Johnson–Segalman fluid in a curved channel with compliant walls

2012 ◽  
Vol 17 (3) ◽  
pp. 297-311 ◽  
Author(s):  
Sadia Hina ◽  
Tasawar Hayat ◽  
Saleem Asghar
Keyword(s):  
Reynolds Number ◽  
Mathematical Analysis ◽  
Channel Wall ◽  
Peristaltic Flow ◽  
Weissenberg Number ◽  
Peristaltic Transport ◽  
Curved Channel ◽  
Long Wavelength ◽  
Channel Effects ◽  
Compliant Walls

The present investigation deals with the peristaltic flow of an incompressible Johnson–Segalman fluid in a curved channel. Effects of the channel wall properties are taken into account. The associated equations for peristaltic flow in a curved channel are modeled. Mathematical analysis is simplified under long wavelength and low Reynolds number assumptions. The solution expressions are established for small Weissenberg number. Effects of several embedded parameters on the flow quantities are discussed.

Sign in / Sign up

Export Citation Format

Share Document