Three-dimensional aspects of fluid flows in channels. II. Effects of meniscus and thin film regimes on viscous fingers

We perform a three-dimensional study of steady state viscous fingers that develop in linear channels. By means of a three-dimensional lattice-Boltzmann scheme that mimics the full macroscopic equations of motion of the fluid momentum and order parameter, we study the effect of the thickness of the channel in two cases. First, for total displacement of the fluids in the channel thickness direction, we find that the steady state finger is effectively two-dimensional and that previous two-dimensional results can be recovered by taking into account the effect of a curved meniscus across the channel thickness as a contribution to surface stresses. Second, when a thin film develops in the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement with experimental results. The effect of the thin film renders the problem three-dimensional and results deviate from the two-dimensional prediction.

Three-dimensional aspects of fluid flows in channels. I. Meniscus and thin film regimes

We study the forced displacement of a fluid-fluid interface in a three-dimensional channel formed by two parallel solid plates. Using a lattice-Boltzmann method, we study situations in which a slip velocity arises from diffusion effects near the contact line. The difference between the slip and channel velocities determines whether the interface advances as a meniscus or a thin film of fluid is left adhered to the plates. We find that this effect is controlled by the capillary and Péclet numbers. We estimate the crossover from a meniscus to a thin film and find good agreement with numerical results. The penetration regime is examined in the steady state. We find that the occupation fraction of the advancing finger relative to the channel thickness is controlled by the capillary number and the viscosity contrast between the fluids. For high viscosity contrast, lattice-Boltzmann results agree with previous results. For zero viscosity contrast, we observe remarkably narrow fingers. The shape of the finger is found to be universal.

Synovial fluid examination for the diagnosis of synovial amyloidosis in patients with chronic renal failure undergoing haemodialysis

The diagnosis of synovial amyloidosis is based upon synovial biopsy. Synovial fluid (SF) in seven patients with amyloid arthropathy associated with chronic renal failure undergoing haemodialysis were studied. The SF and synovial samples of 10 consecutive patients with seronegative mono- or oligoarthritis served as controls. Six of the seven patients with amyloid positive synovial biopsy specimens showed amyloid in their SF. No amyloid was found in the synovial tissue or fluid of the 10 patients in the control group, the sensitivity being 87.7%. The finding of amyloid in SF was highly reproducible, showing its presence in the same joint on several occasions. The deposits were Congophilia resistant to potassium permanganate pretreatment, and the immunohistochemical analysis proved that they contained beta 2 microglobulin. The high sensitivity and good reproducibility of the method shows that the finding of amyloid in SF is sufficient for the diagnosis of synovial amyloidosis. It is possible to perform immunohis

Quasi-two-dimensional electrodeposition under forced fluid flow

Experimental quasi-two-dimensional Zn electrodeposits are grown under forced convection conditions. Large-scale effects, with preferential growth towards the impinging flow, together with small-scale roughness suppression effects are evidenced and separately analyzed by using two different radial cell configurations. Interpretations are given in terms of primary concepts concerning current and concentration distributions.

Inertial effects on reactive particles advected by turbulence

We study the problem of the advection of passive particles with inertia in a two-dimensional, synthetic, and stationary turbulent flow. The asymptotic analytical result and numerical simulations show the importance of inertial bias in collecting the particles preferentially in certain regions of the flow, depending on their density relative to that of the flow. We also study how these aggregates are affected when a simple chemical reaction mechanism is introduced through a Eulerian scheme. We find that inertia can be responsible for maintaining a stationary concentration pattern even under nonfavorable reactive conditions or destroying it under favorable ones.

Chemoconvection patterns in the Methylene-blue-glucose system: Weakly nonlinear analysis

The oxidation of solutions of glucose with methylene-blue as a catalyst in basic media can induce hydrodynamic overturning instabilities, termed chemoconvection in recognition of their similarity to convective instabilities. The phenomenon is due to gluconic acid, the marginally dense product of the reaction, which gradually builds an unstable density profile. Experiments indicate that dominant pattern wavenumbers initially increase before gradually decreasing or can even oscillate for long times. Here, we perform a weakly nonlinear analysis for an established model of the system with simple kinetics, and show that the resulting amplitude equation is analogous to that obtained in convection with insulating walls. We show that the amplitude description predicts that dominant pattern wavenumbers should decrease in the long term, but does not reproduce the aforementioned increasing wavenumber behavior in the initial stages of pattern development. We hypothesize that this is due to horizontally homogeneous steady states not being attained before pattern onset. We show that the behavior can be explained using a combination of pseudo-steady-state linear and steady-state weakly nonlinear theories. The results obtained are in qualitative agreement with the analysis of experiments.

The effect of lattice defects on Hele-Shaw flow over an etched lattice

We examine the patterns formed by injecting nitrogen gas into the center of a horizontal, radial Hele-Shaw cell filled with paraffin oil. We use smooth plates and etched plates with lattices having different amounts of defects (010 %). In all cases, a quantitative measure of the pattern ramification shows a regular trend with injection rate and cell gap, such that the dimensionless perimeter scales with the dimensionless time. By adding defects to the lattice, we observe increased branching in the pattern morphologies. However, even in this case, the scaling behavior persists. Only the prefactor of the scaling function shows a dependence on the defect density. For different lattice defect densities, we examine the nature of the different morphology phases.

Nonlinear chemoconvection in the Methylene-blue-glucose system

Interfacial hydrodynamic instabilities arise in a range of chemical systems. One mechanism for instability is the occurrence of unstable density gradients due to the accumulation of reaction products. In this paper we conduct two-dimensional nonlinear numerical simulations for a member of this class of system: the methylene-blue¿glucose reaction. The result of these reactions is the oxidation of glucose to a relatively, but marginally, dense product, gluconic acid, that accumulates at oxygen permeable interfaces, such as the surface open to the atmosphere. The reaction is catalyzed by methylene-blue. We show that simulations help to disassemble the mechanisms responsible for the onset of instability and evolution of patterns, and we demonstrate that some of the results are remarkably consistent with experiments. We probe the impact of the upper oxygen boundary condition, for fixed flux, fixed concentration, or mixed boundary conditions, and find significant qualitative differences in solution behavior; structures either attract or repel one another depending on the boundary condition imposed. We suggest that measurement of the form of the boundary condition is possible via observation of oxygen penetration, and improved product yields may be obtained via proper control of boundary conditions in an engineering setting. We also investigate the dependence on parameters such as the Rayleigh number and depth. Finally, we find that pseudo-steady linear and weakly nonlinear techniques described elsewhere are useful tools for predicting the behavior of instabilities beyond their formal range of validity, as good agreement is obtained with the simulations.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.

Electronic and magnetic structure of LaMnO3 from hybrid periodic density-functional theory

The electronic and magnetic structures of the LaMnO3 compound have been studied by means of periodic calculations within the framework of spin polarized hybrid density-functional theory. In order to quantify the role of approximations to electronic exchange and correlation three different hybrid functionals have been used which mix nonlocal Fock and local Dirac-Slater exchange. Periodic Hartree-Fock results are also reported for comparative purposes. The A-antiferromagnetic ground state is properly predicted by all methods including Hartree-Fock exchange. In general, the different hybrid methods provide a rather accurate description of the band gap and of the two magnetic coupling constants, strongly suggesting that the corresponding description of the electronic structure is also accurate. An important conclusion emerging from this study is that the nature of the occupied states near the Fermi level is intermediate between the Hartree-Fock and local density approximation descriptions with a comparable participation of both Mn and O states.

Surface structure and stability of PdZn and PtZn alloys: Density-functional slab model studies

Geometric parameters of binary (1:1) PdZn and PtZn alloys with CuAu-L10 structure were calculated with a density functional method. Based on the total energies, the alloys are predicted to feature equal formation energies. Calculated surface energies of PdZn and PtZn alloys show that (111) and (100) surfaces exposing stoichiometric layers are more stable than (001) and (110) surfaces comprising alternating Pd (Pt) and Zn layers. The surface energy values of alloys lie between the surface energies of the individual components, but they differ from their composition weighted averages. Compared with the pure metals, the valence d-band widths and the Pd or Pt partial densities of states at the Fermi level are dramatically reduced in PdZn and PtZn alloys. The local valence d-band density of states of Pd and Pt in the alloys resemble that of metallic Cu, suggesting that a similar catalytic performance of these systems can be related to this similarity in the local electronic structures.

Antiferromagnetic exchange interactions from hybrid density functional theory

A hybrid theory which combines the full nonlocal ¿exact¿ exchange interaction with the local spin-density approximation of density-functional theory is shown to lead to marked improvement in the description of antiferromagnetically coupled systems. Semiquantitative agreement with experiment is found for the magnitude of the coupling constant in La2CuO4, KNiF3, and K2NiF4. The magnitude of the unpaired spin population on the metal site is in excellent agreement with experiment for La2CuO4.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.