Interplay among unstable modes in films over permeable walls

The stability of open-channel flows (or film flows) has been extensively investigated for the case of impermeable smooth walls. In contrast, despite its relevance in many geophysical and industrial flows, the case that considers a permeable rather than an impermeable wall is almost unexplored. In the present work, a linear stability analysis of a film falling over a permeable and inclined wall is developed and discussed. The focus is on the mutual interaction between three modes of instability, namely, the well-known free-surface and hydrodynamic (i.e. shear) modes, which are commonly observed in open-channel flows over impermeable walls, plus a new one associated with the flow within the permeable wall (i.e. the porous mode). The flow in this porous region is modelled by the volume-averaged Navier-Stokes equations and, at the wall interface, the surface and subsurface flow are coupled through a stress-jump condition, which allows one to obtain a continuous velocity profile throughout the whole flow domain. The generalized eigenvalue problem is then solved via a novel spectral Galerkin method, and the whole spectrum of eigenvalues is presented and physically interpreted. The results show that, in order to perform an analysis with a full coupling between surface and subsurface flow, the convective terms in the volume-averaged equations have to be retained. In previous studies, this aspect has never been considered. For each kind of instability, the critical Reynolds number (${\mathit{Re}}_{c}$) is reported for a wide range of bed slopes ($\theta$) and permeabilities ($\sigma$). The results show that the free-surface mode follows the behaviour that was theoretically predicted by Benjamin and Yih for impermeable walls and is independent of wall permeability. In contrast, the shear mode shows a high dependence on $\sigma$: at $\sigma = 0$ the behaviour of ${\mathit{Re}}_{c} (\theta )$ recovers the well-known non-monotonic behaviour of the impermeable-wall case, with a minimum at \theta \sim 0. 05\textdegree . However, with an increase in wall permeability, ${\mathit{Re}}_{c}$ gradually decreases and eventually recovers a monotonic decreasing behaviour. At high values of $\sigma$, the porous mode of instability also occurs. A physical interpretation of the results is presented on the basis of the interplay between the free-surface-induced perturbation of pressure, the increment of straining due to shear with the increase in slope, and the shear stress condition at the free surface. Finally, the paper investigates the extent to which Squire's theorem is applicable to the problem presented herei

Effect of Formation Hydrodynamics on Mechanical Properties of Container Materials

The objectives of this study were to compare the mechanical and physical properties of the sheets made using the Vortigen technology (a non-conventional technique that creates very high number vortices in a fluid flow mixture of water, fibers, and chemical additives) with those produced from a conventional method of papermaking and to provide insight into the impact of formation hydrodynamics on sheet properties. The results of formation, ultrasonic stiffness, and creep/accelerated creep measurements of the Vortigen sheets as compared with the standard sheets are presented. Samples of Vortigen (V) and standard (S) sheets (4 samples from each group) were obtained from papers produced on a pilot machine. Formation measurements (that provides a measure of density distribution in a sheet) were performed using a formation tester which is based on beta particle absorption. Measurements of creep and accelerated creep were made at a constant relative humidity (RH) of 80% and a cyclic RH between 30% and 80% for strips cut along the machine direction (MD) and cross machine direction (CD, which is perpendicular to MD) directions. There was a significant difference between the distributions of basis weights for the two types of papers. The mean coefficient of variation in grammage for the V samples was 8.97 while that for the S samples was 12.60. The mean MD/CD stiffness ratios for the V and S samples were 1.1 and 1.6, respectively. The mean Z-direction longitudinal specific stiffness corresponding to the V samples were 18% greater than the corresponding value for the S samples. The MD strips from the S samples exhibited the smallest creep while the CD strips from the S samples exhibited the largest creep. Creep values corresponding to the Vortigen sheets were between the extreme values of the standard samples. The results of this study indicated that because of the influence of formation hydrodynamics on fiber orientation and formation, in general, the stiffness properties (and specifically the CD stiffness) of the Vortigen samples were greater than those of the standard samples

Coaxial Jets with Disparate Viscosity: Mixing and Laminarization Characteristics

Mixing of fluids in a coaxial jet is studied under four distinct viscosity ratios, m = 1, 10, 20 and 40, using highly resolved large-eddy simulations (LES), particle image velocimetry and planar laser-induced fluorescence. The accuracy of predictions is tested against data obtained by the simultaneous experimental measurements of velocity and concentration fields. For the highest and lowest viscosity ratios, standard RANS models with unclosed terms pertaining to viscosity variations are employed. We show that the standard Reynolds-averaged Navier-Stokes (RANS) approach with no explicit modelling for variable-viscosity terms is not applicable whereas dynamic LES models provide high-quality agreement with the measurements. To identify the underlying mixing physics and sources of discrepancy in RANS predictions, two distinct mixing modes are defined based on the viscosity ratio. Then, for each mode, the evolution of mixing structures, momentum budget analysis with emphasis on variable-viscosity terms, analysis of the turbulent activity and decay of turbulence are investigated using highly resolved LES data. The mixing dynamics is found to be quite distinct in each mixing mode. Variable viscosity manifests multiple effects that are working against each other. Viscosity gradients induce additional instabilities while increasing overall viscosity decreases the effective Reynolds number leading to laminarization of the turbulent jet, explaining the lack of dispersion and turbulent diffusion. Momentum budget analysis reveals that variable-viscosity terms are significant to be neglected. The scaling of the energy spectrum cascade suggests that in the TLL mode the unsteady laminar shedding is responsible for the eddies observed

Lattice Boltzmann for Binary Fluids with Suspended Colloids

A new description of the binary fluid problem via the lattice Boltzmann method is presented which highlights the use of the moments in constructing two equilibrium distribution functions. This offers a number of benefits, including better isotropy, and a more natural route to the inclusion of multiple relaxation times for the binary fluid problem. In addition, the implementation of solid colloidal particles suspended in the binary mixture is addressed, which extends the solid-fluid boundary conditions for mass and momentum to include a single conserved compositional order parameter. A number of simple benchmark problems involving a single particle at or near a fluid-fluid interface are undertaken and show good agreement with available theoretical or numerical results.Comment: 10 pages, 4 figures, ICMMES 200

Lattice-Boltzmann Method for Geophysical Plastic Flows

We explore possible applications of the Lattice-Boltzmann Method for the simulation of geophysical flows. This fluid solver, while successful in other fields, is still rarely used for geotechnical applications. We show how the standard method can be modified to represent free-surface realization of mudflows, debris flows, and in general any plastic flow, through the implementation of a Bingham constitutive model. The chapter is completed by an example of a full-scale simulation of a plastic fluid flowing down an inclined channel and depositing on a flat surface. An application is given, where the fluid interacts with a vertical obstacle in the channel.Comment: in W. Wu, R.I. Borja (Edts.) Recent advances in modelling landslides and debris flow, Springer Series in Geomechanics and Geoengineering (2014), ISBN 978-3-319-11052-3, pp. 131-14

Real-time flow simulation of indoor environments using lattice Boltzmann method

A novel lattice Boltzmann method (LBM) based 3D computational fluid dynamics (CFD) technique has been implemented on the graphics processing unit (GPU) for the purpose of simulating the indoor environment in real-time. We study the time evolution of the turbulent airflow and temperature inside a test chamber and in a simple model of a four-bed hospital room. The predicted results from LBM are compared with traditional CFD based large eddy simulations (LES). Reasonable agreement between LBM results and LES method is observed with significantly faster computational times

Lattice Boltzmann simulations of soft matter systems

This article concerns numerical simulations of the dynamics of particles immersed in a continuum solvent. As prototypical systems, we consider colloidal dispersions of spherical particles and solutions of uncharged polymers. After a brief explanation of the concept of hydrodynamic interactions, we give a general overview over the various simulation methods that have been developed to cope with the resulting computational problems. We then focus on the approach we have developed, which couples a system of particles to a lattice Boltzmann model representing the solvent degrees of freedom. The standard D3Q19 lattice Boltzmann model is derived and explained in depth, followed by a detailed discussion of complementary methods for the coupling of solvent and solute. Colloidal dispersions are best described in terms of extended particles with appropriate boundary conditions at the surfaces, while particles with internal degrees of freedom are easier to simulate as an arrangement of mass points with frictional coupling to the solvent. In both cases, particular care has been taken to simulate thermal fluctuations in a consistent way. The usefulness of this methodology is illustrated by studies from our own research, where the dynamics of colloidal and polymeric systems has been investigated in both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures, 76 page

Numerical simulations of complex fluid-fluid interface dynamics

Interfaces between two fluids are ubiquitous and of special importance for industrial applications, e.g., stabilisation of emulsions. The dynamics of fluid-fluid interfaces is difficult to study because these interfaces are usually deformable and their shapes are not known a priori. Since experiments do not provide access to all observables of interest, computer simulations pose attractive alternatives to gain insight into the physics of interfaces. In the present article, we restrict ourselves to systems with dimensions comparable to the lateral interface extensions. We provide a critical discussion of three numerical schemes coupled to the lattice Boltzmann method as a solver for the hydrodynamics of the problem: (a) the immersed boundary method for the simulation of vesicles and capsules, the Shan-Chen pseudopotential approach for multi-component fluids in combination with (b) an additional advection-diffusion component for surfactant modelling and (c) a molecular dynamics algorithm for the simulation of nanoparticles acting as emulsifiers.Comment: 24 pages, 12 figure

The USDA Barley Core Collection:Genetic Diversity, Population Structure, and Potential for Genome-Wide Association Studies

New sources of genetic diversity must be incorporated into plant breeding programs if they are to continue increasing grain yield and quality, and tolerance to abiotic and biotic stresses. Germplasm collections provide a source of genetic and phenotypic diversity, but characterization of these resources is required to increase their utility for breeding programs. We used a barley SNP iSelect platform with 7,842 SNPs to genotype 2,417 barley accessions sampled from the USDA National Small Grains Collection of 33,176 accessions. Most of the accessions in this core collection are categorized as landraces or cultivars/breeding lines and were obtained from more than 100 countries. Both STRUCTURE and principal component analysis identified five major subpopulations within the core collection, mainly differentiated by geographical origin and spike row number (an inflorescence architecture trait). Different patterns of linkage disequilibrium (LD) were found across the barley genome and many regions of high LD contained traits involved in domestication and breeding selection. The genotype data were used to define 'mini-core' sets of accessions capturing the majority of the allelic diversity present in the core collection. These 'mini-core' sets can be used for evaluating traits that are difficult or expensive to score. Genome-wide association studies (GWAS) of 'hull cover', 'spike row number', and 'heading date' demonstrate the utility of the core collection for locating genetic factors determining important phenotypes. The GWAS results were referenced to a new barley consensus map containing 5,665 SNPs. Our results demonstrate that GWAS and high-density SNP genotyping are effective tools for plant breeders interested in accessing genetic diversity in large germplasm collections