Structure-function relationships in the nasal cavity of Arctic and subtropical seals.

The heating and moistening of inhaled air, and the cooling and moisture removal from exhaled air, are crucial for the survival of animals under severe environmental conditions. Arctic mammals have evolved specific adaptive mechanisms to retain warmth and water and restrict heat loss during breathing. Here, the role of the porous turbinates of the nasal cavities of Arctic and subtropical seals is studied with this in mind. Mass and energy balance equations are used to compute the time-dependent temperature and water vapor profiles along the nasal passage. A quasi-1D model based on computed tomography images of seal nasal cavities is used in numerical simulations. Measured cross-sectional areas of the air channel and the perimeters of the computed tomography slices along the nasal cavities of the two seal species are used. The model includes coupled heat and vapor transfer at the air-mucus interface and heat transfer at the interfaces between the tissues and blood vessels. The model, which assumes constant blood flow to the nose, can be used to predict the temperature of the exhaled air as a function of ambient temperature. The energy dissipation (entropy production) in the nasal passages was used to measure the relative importance of structural parameters for heat and water recovery. We found that an increase in perimeter led to significant decreases in the total energy dissipation. This is explained by improved conditions for heat and water transfer with a larger complexity of turbinates. Owing to differences in their nasal cavity morphology, the Arctic seal is expected to be advantaged in these respects relative to the subtropical seal.

Capillary washboarding during slow drainage of a frictional fluid.

Numerous natural and industrial processes involve the mixed displacement of liquids, gases and granular materials through confining structures. However, understanding such three-phase flows remains a formidable challenge, despite their tremendous economic and environmental impact. To unveil the complex interplay of capillary and granular stresses in such flows, we consider here a model configuration where a frictional fluid (an immersed sedimented granular layer) is slowly drained out of a horizontal capillary. Analyzing how liquid/air menisci displace particles from such granular beds, we reveal various drainage patterns, notably the periodic formation of dunes, analogous to road washboard instability. Considering the competitive role of friction and capillarity, a 2D theoretical approach supported by numerical simulations of a meniscus bulldozing a front of particles provides quantitative criteria for the emergence of those dunes. A key element is the strong increase of the frictional forces, as the bulldozed particles accumulate and bend the meniscus horizontally. Interestingly, this frictional enhancement with the attack angle is also crucial in small-legged animals' locomotion over granular media.

Frictional fluid instabilities shaped by viscous forces.

Multiphase flows involving granular materials are complex and prone to pattern formation caused by competing mechanical and hydrodynamic interactions. Here we study the interplay between granular bulldozing and the stabilising effect of viscous pressure gradients in the invading fluid. Injection of aqueous solutions into layers of dry, hydrophobic grains represent a viscously stable scenario where we observe a transition from growth of a single frictional finger to simultaneous growth of multiple fingers as viscous forces are increased. The pattern is made more compact by the internal viscous pressure gradient, ultimately resulting in a fully stabilised front of frictional fingers advancing as a radial spoke pattern.

Thermal modeling of the respiratory turbinates in arctic and subtropical seals.

Mammals possess complex structures in their nasal cavities known as respiratory turbinate bones, which help the animal to conserve body heat and water during respiratory gas exchange. We considered the function of the maxilloturbinates of two species of seals, one arctic (Erignathus barbatus), one subtropical (Monachus monachus). By means of a thermo-hydrodynamic model that describes the heat and water exchange in the turbinate region we are able to reproduce the measured values of expired air temperatures in grey seals (Halichoerus grypus), a species for which experimental data are available. At the lowest environmental temperatures, however, this is only possible in the arctic seal, and only if we allow for the possibility of ice forming on the outermost turbinate region. At the same time the model predicts that for the arctic seals, the inhaled air is brought to deep body temperature and humidity conditions in passing the maxilloturbinates. The modeling shows that heat and water conservation go together in the sense that one effect implies the other, and that the conservation is most efficient and most flexible in the typical environment of both species. By controlling the blood flow through the turbinates the arctic seal is able to vary the heat and water conservation substantially at its average habitat temperatures, but not at temperatures around -40 °C. The subtropical species has simpler maxilloturbinates, and our model predicts that it is unable to bring inhaled air to deep body conditions, even in its natural environment, without some congestion of the vascular mucosa covering the maxilloturbinates. Physiological control of both blood flow rate and mucosal congestion is expected to have profound effects on the heat exchange function of the maxilloturbinates in seals.

Dynamics of Dendritic Ice Freezing in Confinement.

We use high-speed photography to observe the dendritic freezing of ice between two closely spaced parallel plates. Measuring the propagation speeds of dendrites, we investigate whether there is a confinement-induced thermal influence upon the speed beyond that provided by a single surface. Plates of thermally insulating plastic and moderately thermally conductive glass are used alone and in combination, at temperatures between -10.6 and -4.8 °C, with separations between 17 and 135 µm wide. No effect of confinement was detected for propagation on glass surfaces, but a possible slowing of propagation speed was seen between insulating plates. The pattern of dendritic growth was also studied, with a change from curving to straight dendrites being strongly associated with a switch from a glass to a plastic substrate.

Thermal dissipation as both the strength and weakness of matter. A material failure prediction by monitoring creep.

In any domain involving some stressed solids, that is, from seismology to general engineering, the strength of matter is a paramount feature to understand. We here discuss the ability of a simple thermally activated sub-critical model, which includes the auto-induced thermal evolution of cracks tips, to predict the catastrophic failure of a vast range of materials. It is in particular shown that the intrinsic surface energy barrier, for breaking the atomic bonds of many solids, can be easily deduced from the slow creeping dynamics of a crack. This intrinsic barrier is however higher than the macroscopic load threshold at which brittle matter brutally fails, possibly as a result of thermal activation and of a thermal weakening mechanism. We propose a novel method to compute this macroscopic energy release rate of rupture, Ga, solely from monitoring slow creep, and we show that this reproduces the experimental values within 50% accuracy over twenty different materials, and over more than four decades of fracture energy.

Active Brownian particles moving through disordered landscapes.

Disordered media are ubiquitous in systems where self-propelled particles are present, ranging from biological settings to synthetic systems, like in active microfluidic devices. Here we investigate the behavior of active Brownian particles that have an internal energy depot and move through a landscape with a quenched frictional disorder. We consider the cases of very fast internal relaxation processes and the limit of strong disorder. Analytical calculations of the mean-square displacement in the fast-relaxation approximation is shown to agree well with numerically integrated energy depot dynamics and predict normal dispersion for a bounded drag coefficient and anomalous dispersion for power-law dependence of the drag on spatial coordinates. Furthermore, we show that in the strongly disordered limit the self-propulsion speed can, for practical purposes, be considered a fluctuating quantity. Distributions of self-propulsion speeds are investigated numerically for different parameter choices.

How heat controls fracture: the thermodynamics of creeping and avalanching cracks.

While of paramount importance in material science, the dynamics of cracks still lacks a complete physical explanation. The transition from their slow creep behavior to a fast propagation regime is a notable key, as it leads to full material failure if the size of a fast avalanche reaches that of the system. We here show that a simple thermodynamics approach can actually account for such complex crack dynamics, and in particular for the non-monotonic force-velocity curves commonly observed in mechanical tests on various materials. We consider a thermally activated failure process that is coupled with the production and the diffusion of heat at the fracture tip. In this framework, the rise in temperature only affects the sub-critical crack dynamics and not the mechanical properties of the material. We show that this description can quantitatively reproduce the rupture of two different polymeric materials (namely, the mode I opening of polymethylmethacrylate (PMMA) plates, and the peeling of pressure sensitive adhesive (PSA) tapes), from the very slow to the very fast fracturing regimes, over seven to nine decades of crack propagation velocities. In particular, the fastest regime is obtained with an increase of temperature of thousands of Kelvins, on the molecular scale around the crack tip. Although surprising, such an extreme temperature is actually consistent with different experimental observations that accompany the fast propagation of cracks, namely, fractoluminescence (i.e., the emission of visible light during rupture) and a complex morphology of post-mortem fracture surfaces, which could be due to the sublimation of bubbles.

Seeking minimum entropy production for a tree-like flow-field in a fuel cell.

Common for tree-shaped, space-filling flow-field plates in polymer electrolyte fuel cells is their ability to distribute reactants uniformly across the membrane area, thereby avoiding excess concentration polarization or entropy production at the electrodes. Such a flow field, as predicted by Murray's law for circular tubes, was recently shown experimentally to give a better polarization curve than serpentine or parallel flow fields. In this theoretical work, we document that a tree-shaped flow-field, composed of rectangular channels with T-shaped junctions, has a smaller entropy production than the one based on Murray's law. The width w0 of the inlet channel and the width scaling parameter, a, of the tree-shaped flow-field channels were varied, and the resulting Peclet number at the channel outlets was computed. We show, using 3D hydrodynamic calculations as a reference, that pressure drops and channel flows can be accounted for within a few percents by a quasi-1D model, for most of the investigated geometries. Overall, the model gives lower energy dissipation than Murray's law. The results provide new tools and open up new possibilities for flow-field designs characterized by uniform fuel delivery in fuel cells and other catalytic systems.

Pressure evolution and deformation of confined granular media during pneumatic fracturing.

By means of digital image correlation, we experimentally characterize the deformation of a dry granular medium confined inside a Hele-Shaw cell due to air injection at a constant overpressure high enough to deform it (from 50 to 250 kPa). Air injection at these overpressures leads to the formation of so-called pneumatic fractures, i.e., channels empty of beads, and we discuss the typical deformations of the medium surrounding these structures. In addition we simulate the diffusion of the fluid overpressure into the medium, comparing it with the Laplacian solution over time and relating pressure gradients with corresponding granular displacements. In the compacting medium we show that the diffusing pressure field becomes similar to the Laplace solution on the order of a characteristic time given by the properties of the pore fluid, the granular medium, and the system size. However, before the diffusing pressure approaches the Laplace solution on the system scale, we find that it resembles the Laplacian field near the channels, with the highest pressure gradients on the most advanced channel tips and a screened pressure gradient behind them. We show that the granular displacements more or less always move in the direction against the local pressure gradients, and when comparing granular velocities with pressure gradients in the zone ahead of channels, we observe a Bingham type of rheology for the granular paste (the mix of air and beads), with an effective viscosity µ_{B} and displacement thresholds ∇[over ⃗]P_{c} evolving during mobilization and compaction of the medium. Such a rheology, with disorder in the displacement thresholds, could be responsible for placing the pattern growth at moderate injection pressures in a universality class like the dielectric breakdown model with η=2, where fractal dimensions are found between 1.5 and 1.6 for the patterns.

Pneumatic fractures in confined granular media.

We perform experiments where air is injected at a constant overpressure P_{in}, ranging from 5 to 250 kPa, into a dry granular medium confined within a horizontal linear Hele-Shaw cell. The setup allows us to explore compacted configurations by preventing decompaction at the outer boundary, i.e., the cell outlet has a semipermeable filter such that beads are stopped while air can pass. We study the emerging patterns and dynamic growth of channels in the granular media due to fluid flow, by analyzing images captured with a high speed camera (1000 images/s). We identify four qualitatively different flow regimes, depending on the imposed overpressure, ranging from no channel formation for P_{in} below 10 kPa, to large thick channels formed by erosion and fingers merging for high P_{in} around 200 kPa. The flow regimes where channels form are characterized by typical finger thickness, final depth into the medium, and growth dynamics. The shape of the finger tips during growth is studied by looking at the finger width w as function of distance d from the tip. The tip profile is found to follow w(d)∝d^{ß}, where ß=0.68 is a typical value for all experiments, also over time. This indicates a singularity in the curvature d^{2}d/dw^{2}∼κ∼d^{1-2ß}, but not of the slope dw/dd∼d^{ß-1}, i.e., more rounded tips rather than pointy cusps, as they would be for the case ß>1. For increasing P_{in}, the channels generally grow faster and deeper into the medium. We show that the channel length along the flow direction has a linear growth with time initially, followed by a power-law decay of growth velocity with time as the channel approaches its final length. A closer look reveals that the initial growth velocity v_{0} is found to scale with injection pressure as v_{0}∝P_{in}^{3/2}, while at a critical time t_{c} there is a cross-over to the behavior v(t)∝t^{-α}, where α is close to 2.5 for all experiments. Finally, we explore the fractal dimension of the fully developed patterns. For example, for patterns resulting from intermediate P_{in} around 100-150 kPa, we find that the box-counting dimensions lie within the range D_{B}∈[1.53,1.62], similar to viscous fingering fractals in porous media.

Dispersive transport and symmetry of the dispersion tensor in porous media.

The macroscopic laws controlling the advection and diffusion of solute at the scale of the porous continuum are derived in a general manner that does not place limitations on the geometry and time evolution of the pore space. Special focus is given to the definition and symmetry of the dispersion tensor that is controlling how a solute plume spreads out. We show that the dispersion tensor is not symmetric and that the asymmetry derives from the advective derivative in the pore-scale advection-diffusion equation. When flow is spatially variable across a voxel, such as in the presence of a permeability gradient, the amount of asymmetry can be large. As first shown by Auriault [J.-L. Auriault et al. Transp. Porous Med. 85, 771 (2010)TPMEEI0169-391310.1007/s11242-010-9591-y] in the limit of low Péclet number, we show that at any Péclet number, the dispersion tensor D_{ij} satisfies the flow-reversal symmetry D_{ij}(+q)=D_{ji}(-q) where q is the mean flow in the voxel under analysis; however, Reynold's number must be sufficiently small that the flow is reversible when the force driving the flow changes sign. We also demonstrate these symmetries using lattice-Boltzmann simulations and discuss some subtle aspects of how to measure the dispersion tensor numerically. In particular, the numerical experiments demonstrate that the off-diagonal components of the dispersion tensor are antisymmetric which is consistent with the analytical dependence on the average flow gradients that we propose for these off-diagonal components.

Onsager symmetry from mesoscopic time reversibility and the hydrodynamic dispersion tensor for coarse-grained systems.

Onsager reciprocity relations derive from the fundamental time reversibility of the underlying microscopic equations of motion. This gives rise to a large set of symmetric cross-coupling phenomena. We here demonstrate that different reciprocity relations may arise from the notion of mesoscopic time reversibility, i.e., reversibility of intrinsically coarse-grained equations of motion. We use Brownian dynamics as an example of such a dynamical description and show how it gives rise to reciprocity in the hydrodynamic dispersion tensor as long as the background flow velocity is reversed as well.

Minimal model for anomalous diffusion.

A random walk model with a local probability of removal is solved exactly and shown to exhibit subdiffusive behavior with a mean square displacement the evolves as 〈x^{2}(t)〉∼t^{1/2} at late times. This model is shown to be well described by a diffusion equation with a sink term, which also describes the evolution of a pressure or temperature field in a leaky environment. For this reason a number of physical processes are shown to exhibit anomalous diffusion. The presence of the sink term is shown to change the late time behavior of the field from 1/t^{1/2} to 1/t^{3/2}.

Note: Localization based on estimated source energy homogeneity.

Acoustic signal localization is a complex problem with a wide range of industrial and academic applications. Herein, we propose a localization method based on energy attenuation and inverted source amplitude comparison (termed estimated source energy homogeneity, or ESEH). This inversion is tested on both synthetic (numerical) data using a Lamb wave propagation model and experimental 2D plate data (recorded with 4 accelerometers sensitive up to 26 kHz). We compare the performance of this technique with classic source localization algorithms: arrival time localization, time reversal localization, and localization based on energy amplitude. Our technique is highly versatile and out-performs the conventional techniques in terms of error minimization and cost (both computational and financial).

Dynamic aerofracture of dense granular packings.

A transition in hydraulically induced granular displacement patterns is studied by means of discrete numerical molecular dynamics simulations. During this transition the patterns change from fractures and fingers to finely dispersed bubbles. The dynamics of the displacement patterns are studied in a rectangular Hele-Shaw cell filled with a dense but permeable two-dimensional granular layer. At one side of the cell the pressure of the compressible interstitial gas is increased. At the opposite side from the inlet of the cell a semipermeable boundary is located. This boundary is only permeable towards the gas phase while preventing grains from leaving the cell. The imposed pressure gradient compacts the grains. In the process we can identify and describe a mechanism that controls the transition of the emerging displacement patterns from fractures and fingers to finely dispersed bubbles as a function of the interstitial gas's properties and the characteristics of the granular phase.

Mixing of a granular layer falling through a fluid.

We analyze the granular Rayleigh-Taylor instability of densely packed grains immersed in a compressible or an incompressible fluid using numerical simulations and two types of experiments. The simulations are based on a two-dimensional (2D) molecular dynamics model and the experiments have been carried out in systems of grains immersed in water/glycerol (incompressible fluid) and in air (compressible fluid). The variation of the interstitial fluid is shown to generate different dynamical patterns and mixing properties of the granular systems. The results have been quantified using 2D autocorrelation functions, the power spectrum of the velocity field and velocity field histograms. Excellent agreement is found between the numerical simulations and the experiments.

Oscillation-induced displacement patterns in a two-dimensional porous medium: a lattice Boltzmann study.

We present a numerical study of the statistical behavior of a two-phase flow in a two-dimensional porous medium subjected to an oscillatory acceleration transverse to the overall direction of flow. A viscous nonwetting fluid is injected into a porous medium filled with a more viscous wetting fluid. During the whole process sinusoidal oscillations of constant amplitude and frequency accelerates the porous medium sideways, perpendicular to the overall direction of flow. The invasion process displays a transient behavior where the saturation of the defending fluid decreases, before it enters a state of irreducible wetting fluid saturation, where there is no net transport of defending fluid toward the outlet of the system. In this state the distribution of sizes of the remaining clusters are observed to obey a power law with an exponential cutoff. The cutoff cluster size is found to be determined by the flow and oscillatory stimulation parameters. This cutoff size is also shown to be directly related to the extracted amount of defending fluid. Specifically, the results show that the oscillatory acceleration of the system leads to potentially a large increase in extracted wetting fluid.

Size invariance of the granular Rayleigh-Taylor instability.

The size scaling behavior of the granular Rayleigh-Taylor instability [J. L. Vinningland, Phys. Rev. Lett. 99, 048001 (2007)] is investigated experimentally, numerically, and theoretically. An upper layer of grains displaces a lower gap of air by organizing into dense fingers of falling grains separated by rising bubbles of air. The dependence of these structures on the system and grain sizes is investigated. A spatial measurement of the finger structures is obtained by the Fourier power spectrum of the wave number k. As the size of the grains increases the wave number decreases accordingly which leaves the dimensionless product of wave number and grain diameter, dk, invariant. A theoretical interpretation of the invariance, based on the scaling properties of the model equations, suggests a gradual breakdown of the invariance for grains smaller than approximately 70 microm or greater than approximately 570 microm in diameter.

Sedimentation instabilities: impact of the fluid compressibility and viscosity.

The effect of an interstitial fluid on the mixing of sedimenting grains is studied numerically in a closed rectangular Hele-Shaw cell. We investigate the impact of the fluid compressibility and fluid viscosity on the dynamics and structures of the granular Rayleigh-Taylor instability. First we discuss the effect of the fluid compressibility on the initial fluid pressure evolution and on the dynamics of the particles. Here, the emerging patterns do not seem highly affected by the compressibility change studied. To characterize the patterns and motion the combined length of the particle trajectories in relation to the movement of the center of mass is analyzed, and the separation of particle pairs is measured as a function of the fluid viscosity.