Search

Most studies of triple flames in counterflowing streams of fuel and oxidizer have been focused on the symmetric problem in which the stoichiometric mixture fraction is $1/2$. There then exist lean and rich premixed flames of roughly equal strengths, with a diffusion flame trailing behind from the stoichiometric point at which they meet. In the majority of realistic situations, however, the stoichiometric mixture fraction departs appreciably from unity, typically being quite small. With the objective of clarifying the influences of stoichiometry, attention is focused on one of the simplest possible models, addressed here mainly by numerical integration. When the stoichiometric mixture fraction departs appreciably from $1/2$, one of the premixed wings is found to be dominant to such an extent that the diffusion flame and the other premixed flame are very weak by comparison. These curved, partially premixed flames are expected to be relevant in realistic configurations. In addition, a simple kinematic balance is shown to predict the shape of the front and the propagation velocity reasonably well in the limit of low stretch and low curvature.

Numerical computations employing detailed chemistry are used to characterize the different combustion modes emerging in mixing layers separating nitrogen-diluted counterflowing planar streams of hydrogen and oxygen. Attention is focused on high degrees of dilution, resulting in near-limit flames, with peak temperatures close to the crossover temperature. A bifurcation diagram is presented in a plane, having the stoichiometric mixture fraction and normalized strain rate as coordinates, that identifies six different combustion regimes involving four different flame types, namely, diffusion-flame sheets, advancing and retreating edge flames, multiple flame tubes, and single isolated flame tubes. Multiple-tube flame configurations vary from small, round, widely separated flame strings at high strain rates to wide, flat, densely packed flame strips, with narrow flame-free gaps between them, at lower strain rates, and they are steady and stable in various arrays over a continuum of tube-separation distances. The observed flame behavior exhibits hysteresis in a certain range of parameters, with the structure that is established depending on the ignition mechanism, as it also does at high strain rates, and a continuum of different stable steady-state flame configurations exists, each accessed from a different initial condition.

Directional fluid flow in perivascular spaces surrounding cerebral arteries is hypothesized to play a key role in brain solute transport and clearance. While various drivers for pulsatile flow, such as cardiac or respiratory pulsations, are well quantified, the question remains as to which mechanisms could induce directional flow within physiological regimes. To address this question, we develop theoretical and numerical reduced-order models to quantify the directional (net) flow induceable by peristaltic pumping in periarterial networks. Each periarterial element is modeled as a slender annular space bounded internally by a circular tube supporting a periodic traveling (peristaltic) wave. Under the reasonable assumptions of small Reynolds number flow, small radii, and small-amplitude peristaltic waves, we use lubrication theory and regular perturbation methods to derive theoretical expressions for the directional net flow and pressure distribution in the perivascular network. The reduced model is used to derive closed-form analytical expressions for the net flow for simple network configurations of interest, including single elements, two elements in tandem, and a three element bifurcation, with results compared with numerical predictions. In particular, we provide a computable theoretical estimate of the net flow induced by peristaltic motion in perivascular networks as a function of physiological parameters, notably wave length, frequency, amplitude and perivascular dimensions. Quantifying the maximal net flow for specific physiological regimes, we find that vasomotion may induce net pial periarterial flow velocities on the order of a few to tens of mum/s and that sleep-related changes in vasomotion pulsatility may drive a threefold flow increase.

The buoyancy-driven boundary-layer flow that develops over a semi-infinite inclined hot plate is known to become unstable at a finite distance from the leading edge, characterized by a critical value of the Grashof number Gr based on the local boundary-layer thickness. The nature of the resulting instability depends on the inclination angle $\phi$, measured from the vertical direction. For values of $\phi$ below a critical value $\phi_c$ the instability is characterized by the appearance of spanwise traveling waves, whereas for $\phi>\phi_c$ the bifurcated flow displays G\"ortler-like streamwise vortices. The Boussinesq approximation, employed in previous linear stability analyses, ceases to be valid for gaseous flow when the wall-to-ambient temperature ratio $\Theta_w$ is not close to unity. The corresponding non-Boussinesq analysis is presented here, accounting also for the variation with temperature of the different transport properties. A temporal stability analysis including nonparallel effects of the base flow is used to determine curves of neutral stability, which are then employed to delineate the dependences of the critical Grashof number and of its associated wave length on the inclination angle $\phi$ and on the temperature ratio $\Theta_w$ for the two instability modes, giving quantitative information of interest for configurations with $\Theta_w-1\sim 1$. The analysis provides in particular the predicted dependence of the crossover inclination angle $\phi_c$ on $\Theta_w$, indicating that for gaseous flow with $\Theta_w-1\sim 1$ spanwise traveling waves are predominant over a range of inclination angles $0 \le \phi \le \phi_c$ that is significantly wider than that predicted in the Boussinesq approximation.

This paper investigates the initiation of a deflagration in a premixed boundary-layer stream by continuous heat deposition from a line energy source placed perpendicular to the flow on the wall surface, a planar flow configuration relevant for small-scale combustion applications, including portable rotary engines. Ignition is investigated in the constant density approximation with a one-step irreversible reaction with large activation energy adopted for the chemistry description. The ratio of the characteristic strain time, given by the inverse of the wall velocity gradient, to the characteristic deflagration residence time defines the relevant controlling Damk\"ohler number $D$. The time-dependent evolution following the activation of the heat source is obtained by numerical integration of the energy and fuel conservation equations. For sufficiently small values of $D$, the solution is seen to evolve towards a steady flow in which the chemical reaction remains confined to a finite near-source reactive kernel. This becomes increasingly slender for increasing values of $D$, corresponding to smaller near-wall velocities, until a critical value $D_{c1}$ is reached at which the confined kernel is replaced by a steady anchored deflagration, assisted by the source heating rate, which develops indefinitely downstream. As the boundary-layer velocity gradient is further decreased, a second critical Damk\"ohler number $D_{c2}>D_{c1}$ is reached at which the energy deposition results in a flashback deflagration propagating upstream against the incoming flow along the base of the boundary layer. The computations are extended to investigate the dependence of $D_{c1}$ and $D_{c2}$ on the fuel diffusivity as well as the additional dependence of $D_{c1}$ on the source heating rate, thereby delineating the boundaries that define the relevant regime diagram for these combustion systems.

A method for computing minimum ignition energies for gaseous fuel mixtures with detailed and reduced chemistry, by numerical integration of time-dependent conservation equations in a spherically symmetrical configuration, is presented and discussed, testing its general characteristics and accuracy. The method is applied to methanol-air mixtures described by a 38-step Arrhenius chemistry description and by an 8-step chemistry description based on steady-state approximations for reaction intermediaries. Comparisons of predictions with results of available experimental measurements produced reasonable agreements and supported both the robustness of the computational method and the usefulness of the 8-step reduction in achieving accurate predictions.

This paper investigates the steady axisymmetric structure of the cold boundary-layer flow surrounding fire whirls developing over localized fuel sources lying on a horizontal surface. The inviscid swirling motion found outside the boundary layer, driven by the entrainment of the buoyant turbulent plume of hot combustion products that develops above the fire, is described by an irrotational solution, obtained by combining Taylor's self-similar solution for the motion in the axial plane with the azimuthal motion induced by a line vortex of circulation $2 \pi \Gamma$. The development of the boundary layer from a prescribed radial location is determined by numerical integration for different swirl levels, measured by the value of the radial-to-azimuthal velocity ratio $\sigma$ at the initial radial location. As in the case $\sigma=0$, treated in the seminal boundary-layer analysis of Burggraf et al. (Phys. Fluids, 1971), the pressure gradient associated with the centripetal acceleration of the inviscid flow is seen to generate a pronounced radial inflow. Specific attention is given to the terminal shape of the boundary-layer velocity near the axis, which displays a three-layered structure that is described by matched asymptotic expansions. The resulting composite expansion, dependent on the level of ambient swirl through the parameter $\sigma$, is employed as boundary condition to describe the deflection of the boundary-layer flow near the axis to form a vertical swirl jet. Numerical solutions of the resulting non-slender collision region for different values of $\sigma$ are presented both for inviscid flow and for viscous flow with moderately large values of the controlling Reynolds number $\Gamma/\nu$. The velocity description provided is useful in mathematical formulations of localized fire-whirl flows, providing consistent boundary conditions accounting for the ambient swirl level.