Mixtures of foam and paste: suspensions of bubbles in yield stress fluids

We study the rheological behavior of mixtures of foams and pastes, which can be described as suspensions of bubbles in yield stress fluids. Model systems are designed by mixing monodisperse aqueous foams and concentrated emulsions. The elastic modulus of the bubble suspensions is found to depend on the elastic capillary number $\textit{Ca}_{_G}$ , defined as the ratio of the paste elastic modulus to the bubble capillary pressure. For values of $\textit{Ca}_{_G}$ larger than $\simeq 0.5$ , the dimensionless elastic modulus of the aerated material decreases as the bubble volume fraction $\phi $ increases, suggesting that bubbles behave as soft elastic inclusions. Consistently, this decrease is all the sharper as $\textit{Ca}_{_G}$ is high, which accounts for the softening of the bubbles as compared to the paste. By contrast, we find that the yield stress of most studied materials is not modified by the presence of bubbles. This suggests that their plastic behavior is governed by the plastic capillary number $\textit{Ca}_{\tau_y}$ , defined as the ratio of the paste yield stress to the bubble capillary pressure. At low $\textit{Ca}_{\tau_y}$ values, bubbles behave as nondeformable inclusions, and we predict that the suspension dimensionless yield stress should remain close to unity, in agreement with our data up to $\textit{Ca}_{\tau_y}=0.2$ . When preparing systems with a larger target value of $\textit{Ca}_{\tau_y}$ , we observe bubble breakup during mixing, which means that they have been deformed by shear. It then seems that a critical value $\textit{Ca}_{\tau_y}\simeq 0.2$ is never exceeded in the final material. These observations might imply that, in bubble suspensions prepared by mixing a foam and a paste, the suspension yield stress is always close to that of the paste surrounding the bubbles. Finally, at the highest $\phi $ investigated, the yield stress is shown to increase abruptly with $\phi $ : this is interpreted as a “foamy yield stress fluid” regime, which takes place when the paste mesoscopic constitutive elements (here, the oil droplets) are strongly confined in the films between the bubbles.     

Scaling the normal stresses in concentrated non-colloidal suspensions of spheres

We have collected six sets of results for the ratio of the second normal stress difference to the shear stress $(N_{2}/\tau )$ in non-colloidal suspensions of spheres in Newtonian matrices. They all show a near-cubic dependence on the volume fraction $\varphi $ in the range $0.1 < \varphi < 0.5$ , in contrast to the square law predictions of Brady and Morris (J Fluid Mech 348:103–139, 1997) for dilute suspensions. We suggest that the difference can be resolved by using a dependence on the square of the effective volume fraction $\varphi _{\textrm e}$ , and good agreement is then found.     

Unsteady flow with separation behind a shock wave diffracting over curved walls

The unsteady separation of the compressible flow field behind a diffracting shock wave was investigated along convex curved walls, using shock tube experimentation at large length and time scales, complemented by numerical computation. Tests were conducted at incident shock Mach numbers of $M_{\hbox {s}} =$ 1.5 and 1.6 over a 100 mm radius wall over a dimensionless time range up to $\tau \le $ 6.45. The development of the near wall flow at $M_{\hbox {s}} =$ 1.5 has been described in detail and is very similar to that observed for slightly lower $\tau $ ’s at $M_{\hbox {s}} =$ 1.6. Computations were performed at wall radii of 100 and 200 mm and for incident shock Mach numbers from 1.5 up to and including Mach 2.0. Comparing dimensionless times for different size walls shows that for a given value of $\tau $ the flow field is very similar for the various wall radii published to date and tested in this study. Previously published results that were examined alongside the results from this study had typical values of $1.6 < \tau < 3.2$ . At the later times presented here, flow features were observed that previously had only been observed at higher Mach numbers. The larger length scales allowed for a degree of Reynolds number independence in the results published here. The effect of turbulence on the numerical and experimental results could not be adequately examined due to limitations of the flow imaging system used and a number of questions remain unanswered.     

A Model for Spontaneous Imbibition as a Mechanism for Oil Recovery in Fractured Reservoirs

The flow of oil and water in naturally fractured reservoirs (NFR) can be highly complex and a simplified model is presented to illustrate some main features of this flow system. NFRs typically consist of low-permeable matrix rock containing a high-permeable fracture network. The effect of this network is that the advective flow bypasses the main portions of the reservoir where the oil is contained. Instead capillary forces and gravity forces are important for recovering the oil from these sections. We consider a linear fracture which is symmetrically surrounded by porous matrix. Advective flow occurs only along the fracture, while capillary driven flow occurs only along the axis of the matrix normal to the fracture. For a given set of relative permeability and capillary pressure curves, the behavior of the system is completely determined by the choice of two dimensionless parameters: (i) the ratio of time scales for advective flow in fracture to capillary flow in matrix $\alpha =\tau ^f/\tau ^m$ ; (ii) the ratio of pore volumes in matrix and fracture $\beta =V^m/V^f$ . A characteristic property of the flow in the coupled fracture–matrix medium is the linear recovery curve (before water breakthrough) which has been referred to as the “filling fracture” regime Rangel-German and Kovscek (J Pet Sci Eng 36:45–60, 2002), followed by a nonlinear period, referred to as the “instantly filled” regime, where the rate is approximately linear with the square root of time. We derive an analytical solution for the limiting case where the time scale $\tau ^{m}$ of the matrix imbibition becomes small relative to the time scale $\tau ^{f}$ of the fracture flow (i.e., $\alpha \rightarrow \infty $ ), and verify by numerical experiments that the model will converge to this limit as $\alpha $ becomes large. The model provides insight into the role played by parameters like saturation functions, injection rate, volume of fractures versus volume of matrix, different viscosity relations, and strength of capillary forces versus injection rate. Especially, a scaling number $\omega $ is suggested that seems to incorporate variations in these parameters. An interesting observation is that at $\omega =1$ there is little to gain in efficiency by reducing the injection rate. The model can be used as a tool for interpretation of laboratory experiments involving fracture–matrix flow as well as a tool for testing different transfer functions that have been suggested to use in reservoir simulators.     

Numerical and experimental investigations of pseudo-shock systems in a planar nozzle: impact of bypass mass flow due to narrow gaps

During previous investigations on pseudo-shock systems, we have observed reproducible differences between measurement and simulations for the pressure distribution as well as for size and shape of the pseudo-shock system. A systematic analysis of the deviations leads to the conclusion that small gaps of $\Delta z=O(10^{-4})$  m between quartz glass side walls and metal contour of the test section are responsible for this mismatch. This paper describes a targeted experimental and numerical study of the bypass mass flow within these gaps and its interaction with the main flow. In detail, we analyze how the pressure distribution within the channel as well as the size, shape and oscillation of the pseudo-shock system are affected by the gap size. Numerical simulations are performed to display the flow inside the gaps and to reproduce and explain the experimental results. Numerical and experimental schlieren images of the pseudo-shock system are in good agreement and show that especially the structure of the primary shock is significantly altered by the presence of small gaps. Extensive unsteady flow simulations of the geometry with gaps reveal that the shear layer between subsonic gap flow and supersonic core flow is subject to a Kelvin–Helmholtz instability resulting in small pressure fluctuations. This leads to a shock oscillation with a frequency of $f= O(10^5) \hbox {s}^{-1}$ . The corresponding time scale $\tau $  (s) is 16 times higher than the characteristic time scale $\tau _\delta =\delta /U_\infty $ of the boundary layer given by the ratio of the boundary layer thickness $\delta $ directly ahead of the shock and the undisturbed free stream velocity $U_\infty $ . To assess the reliability of our numerical investigations, the paper includes a grid study as well as an extensive comparison of several RANS turbulence models and their impact on the predicted shape of pseudo-shock systems.     

The role of angular momentum in the laminar motion of viscous fluids

In laminar flow, viscous fluids must exert appropriate elastic shear stresses normal to the flow direction. This is a direct consequence of the balance of angular momentum. There is a limit, however, to the maximum elastic shear stress that a fluid can exert. This is the ultimate shear stress, \(\tau _\mathrm{y}\), of the fluid. If this limit is exceeded, laminar flow becomes dynamically incompatible. The ultimate shear stress of a fluid can be determined from experiments on plane Couette flow. For water at \(20\,^{\circ }\hbox {C}\), the data available in the literature indicate a value of \(\tau _\mathrm{y}\) of about \(14.4\times 10^{-3}\, \hbox {Pa}\). This study applies this value to determine the Reynolds numbers at which flowing water reaches its ultimate shear stress in the case of Taylor–Couette flow and circular pipe flow. The Reynolds numbers thus obtained turn out to be reasonably close to those corresponding to the onset of turbulence in the considered flows. This suggests a connection between the limit to laminar flow, on the one hand, and the occurrence of turbulence, on the other.     

A hereditary partial differential equation with applications in the theory of simple fluids

The linearized boundary-initial history value problem for simple fluids obeying the Coleman-Noll constitutive equation $$S + p\delta = 2\int\limits_0^\infty {m(s)(E(t - s} ) - E(t))ds$$ is considered. Here S is the stress tensor, δ the Kronecker delta, p the constitutively indeterminate mean normal stress, E the infinitesimal strain tensor, and m(s) a material function. The shear relaxation modulus G is defined as (i) $$G(s) = \int\limits_\infty ^s {m(\xi )d\xi .}$$ In this paper it is shown that if G satisfies the assumptions (i) $$G \in C^2 [0,\infty ),{\text{ }}G(s) \to 0{\text{ as }}s \to \infty,$$ (ii) $$( - 1)^k \frac{{d^k G(s)}}{{ds^k }} > 0,{\text{ }}k = 0,1,$$ (iii) $$G''(s) \geqq 0,$$ then the rest state of the fluid is stable in an appropriate “fading memory” norm. The additional assumption (iv) $$ - \int\limits_0^\infty {G'} (s)s^2 ds < \infty$$ yields asymptotic stability.     

Unsteady Conjugate Natural Convection in a Vertical Cylinder Containing a Horizontal Porous Layer: Darcy Model and Brinkman-Extended Darcy Model

Transient natural convection in a vertical cylinder partially filled with a porous media with heat-conducting solid walls of finite thickness in conditions of convective heat exchange with an environment has been studied numerically. The Darcy and Brinkman-extended Darcy models with Boussinesq approximation have been used to solve the flow and heat transfer in the porous region. The Oberbeck–Boussinesq equations have been used to describe the flow and heat transfer in the pure fluid region. The Beavers–Joseph empirical boundary condition is considered at the fluid–porous layer interface with the Darcy model. In the case of the Brinkman-extended Darcy model, the two regions are coupled by equating the velocity and stress components at the interface. The governing equations formulated in terms of the dimensionless stream function, vorticity, and temperature have been solved using the finite difference method. The main objective was to investigate the influence of the Darcy number $10^{-5}\le \hbox {Da}\le 10^{-3}$ , porous layer height ratio $0\le d/L\le 1$ , thermal conductivity ratio $1\le k_{1,3}\le 20$ , and dimensionless time $0\le \tau \le 1000$ on the fluid flow and heat transfer on the basis of the Darcy and non-Darcy models. Comprehensive analysis of an effect of these key parameters on the Nusselt number at the bottom wall, average temperature in the cylindrical cavity, and maximum absolute value of the stream function has been conducted.     

Asymptotic analysis of linearly elastic shells. III. Justification of Koiter's shell equations

We consider as in Parts I and II a family of linearly elastic shells of thickness 2?, all having the same middle surfaceS=?(?)?R ³, whereω?R ² is a bounded and connected open set with a Lipschitz-continuous boundary, and? ∈ ?³ (?;R ³). The shells are clamped on a portion of their lateral face, whose middle line is?(γ ₀), whereγ ₀ is a portion of?ω withlength γ ₀>0. For all?>0, let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $u_i^\varepsilon g^{i,\varepsilon }$ of the points of the shell, obtained by solving the three-dimensional problem; let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $\zeta _i^\varepsilon$ a ⁱ of the points of the middle surfaceS, obtained by solving the two-dimensional model ofW.T. Koiter, which consists in finding $$\zeta ^\varepsilon = \left( {\zeta _i^\varepsilon } \right) \in V_K (\omega ) = \left\{ {\eta = (\eta _\iota ) \in {\rm H}^1 (\omega ) \times H^1 (\omega ) \times H^2 (\omega ); \eta _i = \partial _v \eta _3 = 0 on \gamma _0 } \right\}$$ such that $$\begin{gathered} \varepsilon \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \gamma _{\sigma \tau } (\zeta ^\varepsilon )\gamma _{\alpha \beta } (\eta )\sqrt a dy + \frac{{\varepsilon ^3 }}{3} \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \rho _{\sigma \tau } (\zeta ^\varepsilon )\rho _{\alpha \beta } (\eta )\sqrt a dy \hfill \\ = \mathop \smallint \limits_\omega p^{i,\varepsilon } \eta _i \sqrt a dy for all \eta = (\eta _i ) \in V_K (\omega ), \hfill \\ \end{gathered}$$ where $a^{\alpha \beta \sigma \tau }$ are the components of the two-dimensional elasticity tensor ofS, $\gamma _{\alpha \beta }$ (η) and $\rho _{\alpha \beta }$ (η) are the components of the linearized change of metric and change of curvature tensors ofS, and $p^{i,\varepsilon }$ are the components of the resultant of the applied forces. Under the same assumptions as in Part I, we show that the fields $\frac{1}{{2_\varepsilon }}\smallint _{ - \varepsilon }^\varepsilon u_i^\varepsilon g^{i,\varepsilon } dx_3^\varepsilon$ and $\zeta _i^\varepsilon$ a ⁱ , both defined on the surfaceS, have the same principal part as? → 0, inH ¹ (ω) for the tangential components, and inL ²(ω) for the normal component; under the same assumptions as in Part II, we show that the same fields again have the same principal part as? → 0, inH ¹ (ω) for all their components. For “membrane” and “flexural” shells, the two-dimensional model ofW.T. Koiter is therefore justified.     

Remnant functions

Remnant functions are defined, with \(\kappa = \sigma + \tau + \tfrac{1}{2}\) , by $$R_{\sigma \tau } (z) = [{{\Gamma (\sigma - [\kappa ])} \mathord{\left/ {\vphantom {{\Gamma (\sigma - [\kappa ])} {\Gamma (\sigma )}}} \right. \kern-\nulldelimiterspace} {\Gamma (\sigma )}}]\sum\limits_{r = 1}^\infty {r^{2\tau } \left[\kern-0.15em\left[ {(r^2 + z)^{\sigma - 1} } \right]\kern-0.15em\right]_\kappa }$$ where \(\left[\kern-0.15em\left[ \right]\kern-0.15em\right]_\kappa\) denotes subtraction of sufficiently many terms of the Taylor series in powers of z to yield a convergent sum; for integral σ a factor \([1 + ({z \mathord{\left/ {\vphantom {z {r^2 }}} \right. \kern-0em} {r^2 }})]\) may also enter. These functions arise in various contexts, in particular, in the calculation of uniform remainder terms for the approximation by integrals of sums with singular summands. Differential recurrence relations, Taylor expansions, and various integral representations are obtained. The full asymptotic expansions for ¦z¦→∞ with ¦arg z¦ <π are derived, and it is shown that for integral τ these converge exponentially fast.     

A Comparative Study of Turbulence in Ramp-Up and Ramp-Down Unsteady Flows

DNS of a turbulent channel flow subjected to a step change in pressure gradient are performed to facilitate a direct comparison between ramp-up and ramp-down flows. Strong differences are found between behaviours of turbulence in the two flows. The wall shear stress in the ramp-up flow first overshoots, and then strongly undershoots the quasi-steady value in the initial stage of the excursion, before approaching the quasi-steady value. In a strongly decelerating flow, the wall shear stress tends to first undershoot but then overshoot the quasi-steady value. ??Slow?? response of turbulence as well as flow inertia is responsible for these behaviours. In the ramp-up flow, the response of turbulence is similar to that observed in uniformly accelerating flows from previous studies, exhibiting a three-stage development. However, the transition between the various stages is more gradual and the responding stage is much longer and slower in the flows considered here. It has been shown that the delay in the near wall region is longer than that in the buffer layer confirming that turbulence response first occurs at the location of peak turbulence production. In a strongly decelerating flow, the response of turbulence exhibits a two-stage development. In both ramp-up and ramp down flows, the energy distribution in the three components of turbulent kinetic energy deviates from that of the steady flow. In a ramp-up flow, more energy is in $u_1^\prime $ and less in $u_2^\prime $ and $u_3^\prime $ , whereas the trend is reversed in a ramp-down flow. This is a reflection of the redistribution of turbulence from $u_1^\prime $ to $u_2^\prime $ and $u_3^\prime $ .     

Drainage of Capillary-Trapped Oil by an Immiscible Gas: Impact of Transient and Steady-State Water Displacement on Three-Phase Oil Permeability

In a previous paper (Dehghanpour et al., Phys Rev E 83:065302, 2011a), we showed that relative permeability of mobilized oil, $k_\mathrm{ro}$ , measured during tertiary gravity drainage, is significantly higher than that of the same oil saturation in other tests where oil is initially a continuous phase. We also showed that tertiary $k_\mathrm{ro}$ strongly correlates to both water saturation, $S_\mathrm{w}$ , water flux (water relative permeability), $k_\mathrm{rw}$ , and the change in water saturation with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . To develop a model and understanding of the enhanced oil transport, identifying which of these parameters ( $S_\mathrm{w},\,k_{\mathrm{rw}}$ , or $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ ) plays the controlling role is necessary, but in the previous experiments these could not be deconvolved. To answer the remaining question, we conduct specific three-phase displacement experiments in which $k_{\mathrm{rw}}$ is controlled by applying a fixed water influx, and $S_\mathrm{w}$ develops naturally. We obtain $k_{\mathrm{ro}}$ by using the saturation data measured in time and space. The results suggest that steady-state water influx, in contrast to transient water displacement, does not enhance $k_{\mathrm{ro}}$ . Instead, reducing water influx rate results in excess oil flow. Furthermore, according to our pore scale hydraulic conductivity calculations, viscous coupling and fluid positioning do not sufficiently explain the observed correlation between $k_{\mathrm{ro}}$ and $S_{\mathrm{w}}$ . We conclude that tertiary $k_{\mathrm{ro}}$ is controlled by the oil mobilization rate, which in turn is linked to the rate of water saturation decrease with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . Finally, we develop a simple model which relates tertiary $k_{\mathrm{ro}}$ to transient two-phase gas/water relative permeability.     

Double-Diffusive Natural Convection in Anisotropic Porous Medium Bounded by Finite Thickness Walls: Validity of Local Thermal Equilibrium Assumption

Double-diffusive natural convection in fluid-saturated porous medium inside a vertical enclosure bounded by finite thickness walls with opposing temperature, concentration gradients on vertical walls as well as adiabatic and impermeable horizontal ones has been performed numerically. The Darcy model was used to predict fluid flow inside the porous material, while thermal fields are simulated based on two-energy equations for fluid and solid phases on the basis of a local thermal non-equilibrium model. Computations have been performed for different controlling parameters such as the buoyancy ratio $N$ , the Lewis number Le, the anisotropic permeability ratio $R_\mathrm{p}$ , the fluid-to-solid thermal conductivity ratio $R_\mathrm{c}$ , the interphase heat transfer coefficient $\mathcal{H}$ , the ratio of the wall thickness to its height $D$ , the wall-to-porous medium thermal diffusivity ratio $R_\mathrm{w}$ , and the solid-to-fluid heat capacity ratio $\gamma $ . Thus, the effects of the controlling parameters on heat and mass transfer characteristics are discussed in detail. Moreover, the validity domain of the local thermal equilibrium (LTE) assumption has been delimited for different set of the governing parameters. It has been shown that Le has a noticeable significant effect on fluid temperature profiles and that higher $N$ values lead to a significant enhancement in heat and mass transfer rates. Moreover, for higher $\mathcal{H}, R_\mathrm{c}$ , $R_\mathrm{p}, R_\mathrm{w}$ , or $D$ values and/or lower $\gamma $ values, the solid and fluid phases tend toward LTE.     

Besov Space Regularity Conditions for Weak Solutions of the Navier–Stokes Equations

Consider a bounded domain ${{\Omega \subseteq \mathbb{R}^3}}$ with smooth boundary, some initial value ${{u_0 \in L^2_{\sigma}(\Omega )}}$ , and a weak solution u of the Navier–Stokes system in ${{[0,T) \times\Omega,\,0 < T \le \infty}}$ . Our aim is to develop regularity and uniqueness conditions for u which are based on the Besov space $$B^{q,s}(\Omega ):=\left\{v\in L^2_{\sigma}(\Omega ); \|v\|_{B^{q,s}(\Omega )} := \left(\int\limits^{\infty}_0 \left\|e^{-\tau A}v\right\|^s_q {\rm d} \tau\right)^{1/s}<\infty \right\}$$ with ${{2 < s < \infty,\,3 < q <\infty,\,\frac2{s}+\frac{3}{q} = 1}}$ ; here A denotes the Stokes operator. This space, introduced by Farwig et al. (Ann. Univ. Ferrara 55:89–110, 2009 and J. Math. Fluid Mech. 14: 529–540, 2012), is a subspace of the well known Besov space ${{{\mathbb{B}}^{-2/s}_{q,s}(\Omega )}}$ , see Amann (Nonhomogeneous Navier–Stokes Equations with Integrable Low-Regularity Data. Int. Math. Ser. pp. 1–28. Kluwer/Plenum, New York, 2002). Our main results on the regularity of u exploits a variant of the space ${{B^{q,s}(\Omega )}}$ in which the integral in time has to be considered only on finite intervals (0, δ ) with ${{\delta \to 0}}$ . Further we discuss several criteria for uniqueness and local right-hand regularity, in particular, if u satisfies Serrin’s limit condition ${{u\in L^{\infty}_{\text{loc}}([0,T);L^3_{\sigma}(\Omega ))}}$ . Finally, we obtain a large class of regular weak solutions u defined by a smallness condition ${{\|u_0\|_{B^{q,s}(\Omega )} \le K}}$ with some constant ${{K=K(\Omega, q)>0}}$ .     

Water Vapor Sorption in Cementitious Materials—Measurement,Modeling and Interpretation

The rate and extent of uptake and release of moisture are critical in controlling the behavior of cementitious materials ranging from fluid transport to hygral deformations. While classically determined using an equilibrium (static) salt solution method (Baroghel-Bouny in Cem Concr Res 37:414–437, 2007), advanced capabilities offered by gravimetric dynamic vapor sorption (DVS) analyzers, are now permitting acquisition of sorption spectra at microgram ( $\upmu \hbox {g}$ ) resolution on the order of a few weeks. This work highlights new multicycle determinations of adsorption/desorption isotherms, acquired using a custom-built DVS analyzer for well-hydrated alite and ordinary portland cement pastes over a range of water-to-solid ratios ( $w/s$ , mass basis). Special focus is paid to describe measurement aspects relevant to acquiring reliable spectra, and their interpretation. Sorption isotherms are used to assess transport properties, and sorption hysteresis and its irreversibility following first drying. Based on an optimization-based criterion, the Young-Nelson model is selected to simulate sorption evolutions, including the effects of hysteresis. Sensitivity analyses carried out using this model are used to understand the role of parameters, including porosity and $w/s$ , on the hysteresis that develops from the first to subsequent sorption cycles.     

Numerical stability analysis of numerical methods for volterra integral equations with delay argument

The present paper deals with the stability properties of numerical methods for Volterra integral equations with delay argument. We assess the numerical stability of numerical methods with respect to the following test equations (0.1a) $$y\left( t \right) = \psi \left( 0 \right) + \int_0^t {\left( {py\left( s \right) + q\left( {s - \tau } \right)} \right)ds (0 \leqslant t \leqslant X)}$$ (0.1b) $$y\left( t \right) = \psi \left( t \right) \left( {t \in [ - \tau ,0)} \right)$$ where τ is a positive constant, and P and q are complex valued. We investigate the stability properties of reducible quadrature methods and θ-methods in the case of the above test equations     

Stability and Neimark-Sacker bifurcation analysis for a discrete single genetic negative feedback autoregulatory system with delay

In this paper, we deal with a discrete single genetic negative feedback autoregulatory system with delay by using Euler method. Choosing the delay $\tau $ as the bifurcation parameter and analyzing the associated characteristic equation corresponding to the unique positive fixed point, it is found that the stability of the positive equilibrium and Neimark-Sacker bifurcation may occur when $\tau $ crosses some critical values. Then the explicit formula which determines the stability, direction, and other properties of bifurcating periodic solution is derived by using the center manifold theorem and normal form theory. Finally, in order to illustrate our theoretical analysis, numerical simulations are also included in the end.     

A simple mixture theory for ν Newtonian and generalized Newtonian constituents

This work presents the development of mathematical models based on conservation laws for a saturated mixture of ν homogeneous, isotropic, and incompressible constituents for isothermal flows. The constituents and the mixture are assumed to be Newtonian or generalized Newtonian fluids. Power law and Carreau–Yasuda models are considered for generalized Newtonian shear thinning fluids. The mathematical model is derived for a ν constituent mixture with volume fractions ${\phi_\alpha}$ using principles of continuum mechanics: conservation of mass, balance of momenta, first and second laws of thermodynamics, and principles of mixture theory yielding continuity equations, momentum equations, energy equation, and constitutive theories for mechanical pressures and deviatoric Cauchy stress tensors in terms of the dependent variables related to the constituents. It is shown that for Newtonian fluids with constant transport properties, the mathematical models for constituents are decoupled. In this case, one could use individual constituent models to obtain constituent deformation fields, and then use mixture theory to obtain the deformation field for the mixture. In the case of generalized Newtonian fluids, the dependence of viscosities on deformation field does not permit decoupling. Numerical studies are also presented to demonstrate this aspect. Using fully developed flow of Newtonian and generalized Newtonian fluids between parallel plates as a model problem, it is shown that partial pressures p _α of the constituents must be expressed in terms of the mixture pressure p. In this work, we propose ${p_\alpha=\phi_\alpha p}$ and ${\sum_\alpha^\nu p_\alpha = p}$ which implies ${\sum_\alpha^\nu \phi_\alpha = 1}$ which obviously holds. This rule for partial pressure is shown to be valid for a mixture of Newtonian and generalized Newtonian constituents yielding Newtonian and generalized Newtonian mixture. Modifications of the currently used constitutive theories for deviatoric Cauchy stress tensor are proposed. These modifications are demonstrated to be essential in order for the mixture theory for ν constituents to yield a valid mathematical model when the constituents are the same. Dimensionless form of the mathematical models is derived and used to present numerical studies for boundary value problems using finite element processes based on a residual functional, that is, least squares finite element processes in which local approximations are considered in ${H^{k,p}\left(\bar{\Omega}^e\right)}$ scalar product spaces. Fully developed flow between parallel plates and 1:2 asymmetric backward facing step is used as model problems for a mixture of two constituents.     

Shear-layer instability in the Mach reflection of shock waves

The Kelvin–Helmholtz instability (KHI) is an instability that takes the form of repeating wave-like structures which forms on a shear layer where two adjacent fluids are moving at a relative velocity to one another. Such a shear layer forms in the Mach reflection of shock waves. This work focuses on experimentally visualising the presence of the KHI in Mach reflection as well as its evolution. Experimentation was performed at shock Mach numbers of 1.34, 1.46 and 1.61. Plane test pieces and parabolic profiled pieces followed by a plane section having wedge angles of 30 \(^\circ \) and 38 \(^\circ \) were tested. Flow field visualisation was performed with a schlieren optical system. The KHI was best visualised with the camera-side knife edge perpendicular to the shear layer (i.e. the axis of sensitivity along the length of the shear layer). The structure and growth of the instability were readily identified. The KHI forms more readily with increasing Mach number and wedge angle. Second-order Euler, and Navier–Stokes numerical simulations of the flow field were also conducted. It was found that the Euler and laminar Navier–Stokes solvers achieved very similar results, both producing the KHI, but at a much less developed state than the experimental cases. The k \(-\epsilon \) solver, however, did not produce the instability.