Influence of \text{ CF }_{3}\text{ H } and \text{ CCl }_{4} additives on acetylene detonation

The influence of $\text{ CF }_{3}\text{ H }$ and $\text{ CCl }_{4}$ admixtures (known as detonation suppressors for combustible mixtures) on the development of acetylene detonation was experimentally investigated in a shock tube. The time-resolved images of detonation wave development and propagation were registered using a high-speed streak camera. Shock wave velocity and pressure profiles were measured by five calibrated piezoelectric gauges and the formation of condensed particles was detected by laser light extinction. The induction time of detonation development was determined as the moment of a pressure rise at the end plate of the shock tube. It was shown that $\text{ CF }_{3}\text{ H }$ additive had no influence on the induction time. For $\text{ CCl }_{4}$ , a significant promoting effect was observed. A simplified kinetic model was suggested and characteristic rates of diacetylene $\text{ C }_{4}\text{ H }_{2}$ formation were estimated as the limiting stage of acetylene polymerisation. An analysis of the obtained data indicated that the promoting species is atomic chlorine formed by $\text{ CCl }_{4}$ pyrolysis, which interacts with acetylene and produces $\text{ C }_{2}\text{ H }$ radical, initiating a chain mechanism of acetylene decomposition. The results of kinetic modelling agree well with the experimental data.     

Comparative characteristics of strong shock and detonation waves in bubble media by an electrical wire explosion

Strong shock and detonation waves in inert and chemically active bubble media, which are generated by a wire explosion initiated by a capacitor with a stored energy $W_0 =12.3$ –1,600 J, is experimentally studied. The measurements are performed near the wire and far from the wire in a vertical shock tube 4.5 m long with a volume fraction of the gas in the medium $\beta _0 =1$ –4 %. It is shown that in inert bubble medium, a short intensely decaying shock wave (SW) with intense pressure oscillations is formed in the vicinity of wire explosion point; near the explosion point at $\beta _0 \le 2$  % the SW propagates with the velocity of sound in a liquid. In chemically active bubble medium, an unsteady detonation wave generated by a wire explosion is formed. The pressure amplitude and the velocity of this wave are greater and the length is smaller than those of SW in an inert bubble medium in the same range of explosion energy. It is found that in the interval of low energy explosion from ${\sim }12$ to 64 J, the formation of the bubble detonation wave occurs faster than that at high energies ( $3\times 10^{2}$ – $10^{3}$  J).     

Expansivity and Cone-fields in Metric Spaces

Due to the results of Lewowicz and Tolosa expansivity can be characterized with the aid of Lyapunov function. In this paper we study a similar problem for uniform expansivity and show that it can be described using generalized cone-fields on metric spaces. We say that a function \(f:X\rightarrow X\) is uniformly expansive on a set \(\varLambda \subset X\) if there exist \(\varepsilon >0\) and \(\alpha \in (0,1)\) such that for any two orbits \(\hbox {x}:\{-N,\ldots ,N\} \rightarrow \varLambda \) , \(\hbox {v}:\{-N,\ldots ,N\} \rightarrow X\) of \(f\) we have $$\begin{aligned} \sup _{-N\le n\le N}d(\hbox {x}_n,\hbox {v}_n) \le \varepsilon \implies d(\hbox {x}_0,\hbox {v}_0) \le \alpha \sup _{-N\le n\le N}d(\hbox {x}_n,\hbox {v}_n). \end{aligned}$$ It occurs that a function is uniformly expansive iff there exists a generalized cone-field on \(X\) such that \(f\) is cone-hyperbolic.     

Experimental study on a heavy-gas cylinder accelerated by cylindrical converging shock waves

The Richtmyer–Meshkov instability behavior of a heavy-gas $(\text{ SF }_6)$ cylinder accelerated by a cylindrical converging shock wave is studied experimentally. A curved wall profile is well-designed based on the shock dynamics theory [Phys. Fluids, 22: 041701 (2010)] with an incident planar shock Mach number of 1.2 and a converging angle of $15^\circ $ in a $95\,\text{ mm }\times 95$ mm square cross-section shock tube. The $\text{ SF }_6$ cylinder mixed with the glycol droplets flows vertically through the test section and is illuminated horizontally by a laser sheet. The images obtained only one per run by an ICCD (intensified charge coupled device) combined with a pulsed Nd:YAG laser are first presented and the complete evolution process of the $\text{ SF }_6$ cylinder is then captured in a single test shot by a high-speed video camera combined with a high-power continuous laser. In this way, both the developments of the first counter-rotating vortex pair and the second counter-rotating vortex pair with an opposite rotating direction from the first one are observed. The experimental results indicate that the phenomena induced by the converging shock wave and the reflected shock formed from the center of convergence are distinct from those found in the planar shock case.     

Effect of an impinging shock wave on a partially opened door

The behavior of an aluminum door hanging at the exit of an open shock tube at different angles, from 5 $^\circ $ to 85 $^\circ $ , and thereby providing partially open space for the exiting flow, was investigated experimentally. Experiments were conducted with an incident shock wave Mach number of $M_\mathrm{is}=1.1$ impinging on the partially opened door. Both pressure measurements in the vicinity of the door, on its center and inside the shock tube, and schlieren visualization were undertaken for studying the door movement and its maximum opening angle relative to its initial position. It was found that for an initial opening angle smaller than 25 $^\circ $ the door opened completely while for larger angles its motion is marginal. In addition, for an initial door opening angle of about 10 $^\circ $ the lowest pressures were recorded inside the shock tube behind the evolving waves after exiting of the incident shock wave. The present experimental results may be useful to numerical studies of fluid–structure interactions, e.g., in designing safety valves in jet engines. Such a device is needed for preventing rupture in the case when a sudden overpressure pulse is generated inside the aircraft engine compartment.     

Background-oriented schlieren with natural background for quantitative visualization of open-air explosions

This study describes an attempt of quantitative visualization of open-air explosions via the background-oriented schlieren method (BOS). The shock wave propagation curve and overpressure distribution were extracted from the obtained images and compared with the results of the numerical analysis. The potential of extracting the density distribution behind the shock front is also demonstrated. Two open-air explosions were conducted; one with a $36$ -kg emulsion explosive and the other with a $7.89$ -kg composition C4 explosive. A high-speed digital video camera was used with a frame rate of $10{,}000\,\mathrm{Hz}$ and a pixel size of $800 \times 600$ . A natural background, including trees and grass, was used for BOS measurements instead of the random dots used in a laboratory. The overpressure distribution given by the passing shock was estimated from the visualized images. The estimated overpressures agreed with the values recorded by pressure transducers in the test field. The background displacement caused by light diffraction inside the spherical shock waves was in good agreement, except at the shock front. The results shown here suggest that the BOS method for open-air experiments could provide increasingly better quantitative and conventional visualization results with increasing spatial resolution of high-speed cameras.     

Attenuation of shock waves propagating in polyurethane foams

Shock wave attenuation in polyurethane foams is investigated experimentally and numerically. This study is a part of research project regarding shock propagation in polyurethane foams with high-porosities = 0.951 ~ 0.977 and low densities of ρc = 27.6 ~55.8 kg/m³. Sixty Millimeter long cylindrical foams with various cell numbers and foam insertion condition were installed in a horizontal shock tube of 50 mm i.d. and 5.4 mm in length. Results of pressure measurements in air/foam combination are compared with CFD simulation solving the one-dimensional Euler equations. In the case of a foam B fixed on shock tube wall, pressures at the shock tube end wall increases relatively slowly comparing to non-fixed foam, free to move and a foam A fixed on shock tube wall. This implies that elastic inertia hardly contributes to pressure build up. Pressures behind a foam C fixed on shock tube wall decrease indicating that shock wave is degenerated into compression wave. Dimensionless impulse and attenuation factor decrease as the initial cell number increases. The momentum loss varies depending on cell structure and cell number.     

Microsize and initial pressure effects on shock wave propagation in a tube

In this paper, the shock wave propagation in a channel with a micrometric hydraulic diameter is numerically simulated for an initial Mach number \(M_{s}=2.61\) . The obtained values of the Mach number along the tube are compared to experimental and numerical data given in the literature. The microscale effects on the flow behavior, such as shock wave attenuation and pressure increase behind the shock wave, are amplified by further reducing the scaling ratio (or Reynolds number) of the flow. This reduction is obtained by either decreasing the hydraulic diameter \(D_\mathrm{H}\) or the initial driven gas pressure \(P_1\) . Under these conditions, the flow behavior changes drastically.     

Experiments on the Richtmyer–Meshkov instability with an imposed, random initial perturbation

A vertical shock tube is used to perform experiments on the Richtmyer–Meshkov instability with a three-dimensional random initial perturbation. A membraneless flat interface is formed by opposed gas flows in which the light and heavy gases enter the shock tube from the top and from the bottom of the shock tube driven section. An air/SF $_{6}$ gas combination is used and a Mach number $ M = 1.2$ incident shock wave impulsively accelerates the interface. Initial perturbations on the interface are created by vertically oscillating the gas column within the shock tube to produce Faraday waves on the interface resulting in a short wavelength, three-dimensional perturbation. Planar Mie scattering is used to visualize the flow in which light from a laser sheet is scattered by smoke seeded in the air, and image sequences are captured using three high-speed video cameras. Measurements of the integral penetration depth prior to reshock show two growth behaviors, both having power law growth with growth exponents in the range found in previous experiments and simulations. Following reshock, all experiments show very consistent linear growth with a growth rate in good agreement with those found in previous studies.     

Several new types of bounded wave solutions for the generalized two-component Camassa–Holm equation

Li and Qiao studied the bifurcations and exact traveling wave solutions for the generalized two-component Camassa–Holm equation $$\begin{aligned} \left\{ \begin{array}{l} m_{t}+\sigma um_{x}-Au_{x}+2m \sigma u_{x}+3(1-\sigma )uu_{x}\\ \quad +\rho \rho _{x}=0, \\ \rho _{t} +(\rho u)_{x}=0, \end{array} \right. \end{aligned}$$ \(m=u-u_{xx}, A>0\) . They showed that there exist solitary wave solutions, cusp wave solutions, and periodic wave solutions for the equation, and their analysis focused on the bifurcations when \(\sigma >0\) . In this paper, we first complement the bifurcations when \(\sigma <0\) by following the same procedure as that of Li, and then show the existence and implicit expressions of several new types of bounded wave solutions, including solitary waves, periodic waves, compacton-like waves, and kink-like waves. In addition, the numerical simulations of the bounded wave solutions are given to show the correctness of our results.     

Semianalytical Solution for \text{ CO}_{2} Plume Shape and Pressure Evolution During \text{ CO}_{2} Injection in Deep Saline Formations

The injection of supercritical carbon dioxide ( $\text{ CO}_{2})$ in deep saline aquifers leads to the formation of a $\text{ CO}_{2}$ rich phase plume that tends to float over the resident brine. As pressure builds up, $\text{ CO}_{2}$ density will increase because of its high compressibility. Current analytical solutions do not account for $\text{ CO}_{2}$ compressibility and consider a volumetric injection rate that is uniformly distributed along the whole thickness of the aquifer, which is unrealistic. Furthermore, the slope of the $\text{ CO}_{2}$ pressure with respect to the logarithm of distance obtained from these solutions differs from that of numerical solutions. We develop a semianalytical solution for the $\text{ CO}_{2}$ plume geometry and fluid pressure evolution, accounting for $\text{ CO}_{2}$ compressibility and buoyancy effects in the injection well, so $\text{ CO}_{2}$ is not uniformly injected along the aquifer thickness. We formulate the problem in terms of a $\text{ CO}_{2}$ potential that facilitates solution in horizontal layers, with which we discretize the aquifer. Capillary pressure is considered at the interface between the $\text{ CO}_{2}$ rich phase and the aqueous phase. When a prescribed $\text{ CO}_{2}$ mass flow rate is injected, $\text{ CO}_{2}$ advances initially through the top portion of the aquifer. As $\text{ CO}_{2}$ is being injected, the $\text{ CO}_{2}$ plume advances not only laterally, but also vertically downwards. However, the $\text{ CO}_{2}$ plume does not necessarily occupy the whole thickness of the aquifer. We found that even in the cases in which the $\text{ CO}_{2}$ plume reaches the bottom of the aquifer, most of the injected $\text{ CO}_{2}$ enters the aquifer through the layers at the top. Both $\text{ CO}_{2}$ plume position and fluid pressure compare well with numerical simulations. This solution permits quick evaluations of the $\text{ CO}_{2}$ plume position and fluid pressure distribution when injecting supercritical $\text{ CO}_{2}$ in a deep saline aquifer.     

Pressure drop in alternating curved tubes

Pressure drop measurements in the laminar and turbulent regions for water flowing through an alternating curved circular tube (x=h sin 2πz/λ) are presented. Using the minimum radius of curvature of this curved tube in place of that of the toroidally curved one in calculating the Dean number (ND=Re(D/2R _c )², it is found that the resulting Dean number can help in characterizing this flow. Also, the ratio between the height and length of the tube waves which represents the degree of waveness affects significantly the pressure drop and the transition Dean number. The following correlations have been found:

1. For laminar flow: $$F_w \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} = F_s \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} + 0.03,\operatorname{Re}< 2000.$$
2. For turbulent flow: $$F_w \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} = F_s \left( {\frac{{2R_c }}{D}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} + 0.005,2000< \operatorname{Re}< 15000.$$
3. The transition Dean number: $$ND_{crit} = 5.012 \times 10^3 \left( {\frac{D}{{2R}}} \right)^{2.1} ,0.0111< {D \mathord{\left/ {\vphantom {D {2R_c }}} \right. \kern-\nulldelimiterspace} {2R_c }}< 0.71.$$

     

Dense particle cloud dispersion by a shock wave

A dense particle flow is generated by the interaction of a shock wave with an initially stationary packed granular bed. High-speed particle dispersion research is motivated by the energy release enhancement of explosives containing solid particles. The initial packed granular bed is produced by compressing loose powder into a wafer with a particle volume fraction of $\phi _\mathrm{p} = 0.48$ . The wafer is positioned inside the shock tube, uniformly filling the entire cross-section. This results in a clean experiment where no flow obstructing support structures are present. Through high-speed shadowgraph imaging and pressure measurements along the length of the channel, detailed information about the particle shock interaction was obtained. Due to the limited strength of the incident shock wave, no transmitted shock wave is produced. The initial solid-like response of the particle wafer acceleration forms a series of compression waves that eventually coalesce to form a shock wave. Breakup is initiated along the periphery of the wafer as the result of shear that forms due to the fixed boundary condition. Particle breakup is initiated by local failure sites that result in the formation of particle jets that extend ahead of the accelerating, largely intact, wafer core. In a circular tube, the failure sites are uniformly distributed along the wafer circumference. In a square channel, the failure sites, and the subsequent particle jets, initially form at the corners due to the enhanced shear. The wafer breakup subsequently spreads to the edges forming a highly non-uniform particle cloud.     

Verlauf der Fluidtemperaturen im Querstromrohrbündel

A new analytical method is presented for the determination of temperature distribution and effectiveness of heat transfer in different cross-flow arrangements. As an improvement over the known analytical solutions simplified energy balance equations are used. This simplification consists only in their notation but not in representation of the exact energy balance. The demanded accuracy of the thermal analysis is achieved by proper selection of a mean value of fluid temperature outside the tube. For each tubei in the bundle a mean value ? _i of the dimensionless fluid temperature outside the tube is introduced according to $$\vartheta _i = \omega \vartheta _{i + {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} + (1 - \omega )\vartheta _{i - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} $$ where ? _i?1/2 and ? _i+1/2 are the local temperatures in front and behind thei-th tube, respectively. With symbol ω a weight coefficient is denoted $$\omega = {1 \mathord{\left/ {\vphantom {1 {(1 - e^{ - NTU_{{ \bot \mathord{\left/ {\vphantom { \bot n}} \right. \kern-\nulldelimiterspace} n}} } )}}} \right. \kern-\nulldelimiterspace} {(1 - e^{ - NTU_{{ \bot \mathord{\left/ {\vphantom { \bot n}} \right. \kern-\nulldelimiterspace} n}} } )}} - {n \mathord{\left/ {\vphantom {n {NTU_ \bot }}} \right. \kern-\nulldelimiterspace} {NTU_ \bot }}$$ The number of tube rows in the bundle isn and the number of transfer units of the outside stream is NTU_⊥. Through the introduction of the weight coefficient ω, the mathematical operations related to calculation of the temperature field are radically simplified. This enabled the development of the procedure, valid for three codirected cross-flow arrangements (Fig. 1) and any numbern of the rows.     

Parabolic Systems with p, q-Growth: A Variational Approach

We consider the evolution problem associated with a convex integrand ${f : \mathbb{R}^{Nn}\to [0,\infty)}$ satisfying a non-standard p, q-growth assumption. To establish the existence of solutions we introduce the concept of variational solutions. In contrast to weak solutions, that is, mappings ${u\colon \Omega_T \to \mathbb{R}^n}$ which solve $$ \partial_tu-{\rm div} Df(Du)=0 $$ weakly in ${\Omega_T}$ , variational solutions exist under a much weaker assumption on the gap q ? p. Here, we prove the existence of variational solutions provided the integrand f is strictly convex and $$\frac{2n}{n+2} < p \le q < p+1.$$ These variational solutions turn out to be unique under certain mild additional assumptions on the data. Moreover, if the gap satisfies the natural stronger assumption $$ 2\le p \le q < p+ {\rm min}\big \{1,\frac{4}{n} \big \},$$ we show that variational solutions are actually weak solutions. This means that solutions u admit the necessary higher integrability of the spatial derivative Du to satisfy the parabolic system in the weak sense, that is, we prove that $$u\in L^q_{\rm loc}\big(0,T; W^{1,q}_{\rm loc}(\Omega,\mathbb{R}^N)\big).$$      

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.     

Experimental exploration of underexpanded supersonic jets

Two underexpanded free jets at fully expanded Mach numbers $M_\mathrm{j}$ = 1.15 and 1.50 are studied. Schlieren visualizations as well as measurements of static pressure, Pitot pressure and velocity are performed. All these experimental techniques are associated to obtain an accurate picture of the jet flow development. In particular, expansion, compression and neutral zones have been identified in each shock cell. Particle lag is considered by integrating the equation of motion for particles in a fluid flow and it is found that the laser Doppler velocimetry is suitable for investigating shock-containing jets. Even downstream of the normal shock arising in the $M_\mathrm{j}$ = 1.50 jet, the measured gradual velocity decrease is shown to be relevant.     

Singular Perturbation of Reduced Wave Equation and Scattering from an Embedded Obstacle

We consider time-harmonic wave scattering from an inhomogeneous isotropic medium supported in a bounded domain ${\Omega \subset \mathbb{R}^N}$ (N ≥?2). In a subregion ${D \Subset \Omega}$ , the medium is supposed to be lossy and have a large mass density. We study the asymptotic development of the wave field as the mass density ρ → +?∞ and show that the wave field inside D will decay exponentially while the wave filed outside the medium will converge to the one corresponding to a sound-hard obstacle ${D \Subset \Omega}$ buried in the medium supported in ${\Omega \backslash \overline{D}}$ . Moreover, the normal velocity of the wave field on ? D from outside D is shown to be vanishing as ρ → +?∞. We derive very accurate estimates for the wave field inside and outside D and on ? D in terms of ρ, and show that the asymptotic estimates are sharp. The implication of the obtained results is given for an inverse scattering problem of reconstructing a complex scatterer.