Front shock behavior of stable curved detonation waves in rectangular-cross-section curved channels

The propagation of curved detonation waves of gaseous explosives stabilized in rectangular-cross-section curved channels is investigated. Three types of stoichiometric test gases, C₂H₄ + 3O₂, 2H₂ + O₂, and 2C₂H₂ + 5O₂ + 7Ar, are evaluated. The ratio of the inner radius of the curved channel (r_i) to the normal detonation cell width (λ) is an important factor in stabilizing curved detonation waves. The lower boundary of stabilization is around r_i/λ = 23, regardless of the test gas. The stabilized curved detonation waves eventually attain a specific curved shape as they propagate through the curved channels. The specific curved shapes of stabilized curved detonation waves are approximately formulated, and the normal detonation velocity (D_n)?curvature (κ) relations are evaluated. The D_n nondimensionalized by the Chapman–Jouguet (CJ) detonation velocity (D_CJ) is a function of the κ nondimensionalized by λ. The D_n/D_CJ?λκ relation does not depend on the type of test gas. The propagation behavior of the stabilized curved detonation waves is controlled by the D_n/D_CJ?λκ relation. Due to this propagation characteristic, the fully-developed, stabilized curved detonation waves propagate through the curved channels while maintaining a specific curved shape with a constant angular velocity. Self-similarity is seen in the front shock shapes of the stabilized curved detonation waves with the same r_i/λ, regardless of the curved channel and test gas.     

Critical conditions of hydrogen-air detonation in partially confined geometry

This work presents the results of the large scale experiments with detonation propagating in hydrogen–air mixtures in partially confined geometries. The main aim of the work was to find the critical conditions for detonation propagation in semi-confined geometries with uniform and non-uniform hydrogen–air mixtures. The experimental facility consisted of rectangular 9 × 3 × 0.6 m channel open from the bottom, acceleration section and test section, safety vessel, gas injection and data acquisition system. Sooted plates technique was used as a witness of the detonation. The rectangular channel was placed in a 100 m³ safety vessel. For uniform hydrogen–air mixtures experiments with four different channel heights h were performed: 8, 5, 3 and 2 cm. The critical hydrogen–air mixture height h* for which the detonation may propagate in a layer is close to the 3 cm which corresponds to approximately three detonation cell sizes. For non-uniform hydrogen–air mixture with hydrogen concentration slope equal approximately ?1.1%H₂/cm the critical hydrogen concentration at the top of the layer is approximately equal 26% and the mean detonation layer height is close to the 8.5 cm corresponding to the hydrogen concentration at the bottom of the layer approximately equal 16–17%.     

Fundamental burning velocities of meso-scale propagating spherical flames with H2, CH4 and C3H8 mixtures

This study is performed to experimentally examine the fundamental burning velocity characteristics of meso-scale outwardly propagating spherical laminar flames in the range of flame radius r_f approximately from 1 to 5 mm for hydrogen, methane and propane mixtures, in order to make clear a method for improving combustion of micro–meso scale flames. Macro-scale laminar flames with r_f > 7 mm are also examined for comparison. The mixtures have nearly the same laminar burning velocity (S_L0 = 25 cm/s) for unstretched flames and different equivalence ratios ?. The radius r_f and the burning velocity S_Ll of meso-scale flames are estimated by using sequential schlieren images recorded under appropriate ignition conditions. It is found that S_Ll of hydrogen and methane premixed meso-scale flames at the same r_f or the Karlovitz number Ka shows a tendency to increase with decreasing ?, whereas S_Ll of propane flames increases with ?. However, S_Ll tends to decrease with the Lewis number Le and the Markstein number Ma, irrespective of the type of fuel and ?. It also becomes clear that the optimum flame size and Ka to improve the burning velocity exist for some mixtures depending on Le and fuel types.     

Experimental investigation of flame spreading over pure methane hydrate in a laminar boundary layer

Flame spreading over pure methane hydrate in a laminar boundary layer is investigated experimentally. The free stream velocity (U_∞) was set constant at 0.4 m/s and the surface temperature of the hydrate at the ignition (T_s) was varied between ?10 and ?80 °C. Hydrate particle sizes were smaller than 0.5 mm. Two types of flame spreading were observed; “low speed flame spreading” and “high speed flame spreading”. The low speed flame spreading was observed at low temperature conditions (T_s = ?80 to ?60 °C) and temperatures in which anomalous self-preservation took place (T_s = ?30 to ?10 °C). In this case, the heat transfer from the leading flame edge to the hydrate surface plays a key role for flame spreading. The high speed flame spreading was observed when T_s = ?50 and ?40 °C. At these temperatures, the dissociation of hydrate took place and the methane gas was released from the hydrate to form a thin mixed layer of methane and air with a high concentration gradient over the hydrate. The leading flame edge spread in this premixed gas at a spread speed much higher than laminar burning velocity, mainly due to the effect of burnt gas expansion.     

Near limit behavior of the detonation velocity

The behavior of the detonation velocity near the limits is investigated. Circular tubes of diameters 65, 44 and 13 mm are used. To simulate a quasi two-dimensional rectangular geometry thin annular channels are also used. The annular channels are formed by a 1.5 m long insert of a smaller diameter tube into the larger outer diameter detonation tube. Premixed mixtures of C₂H₂ + 2.5O₂ + 70%Ar, CH₄ + 2O₂ and C₂H₂ + 5N₂O + 50%Ar are used in the present study. The high argon dilution stoichiometric C₂H₂ + 2.5O₂ mixture has a regular cell size and piecewise laminar reaction zone and thus referred to as “stable”. The other two mixtures give highly irregular cell pattern and a turbulent reaction zone and are hence, referred to as “unstable” mixtures. Pressure transducers and optical fibers spaced 10 cm apart along the tube are used for pressure and velocity measurements. Cell size of the three mixtures studied is also determined using smoked foils in both the circular tubes and annular channels. The ratio d/λ (representing the number of cells across the tube diameter) is found to be an appropriate sensitivity parameter to characterize the mixture. The present results indicate that well within the limit, the detonation velocity is generally a few percent below the theoretical Chapman–Jouguet (CJ) value. As the limit is approached, the velocity decreases rapidly to a minimum value before the detonation fails. The narrow range of values of d/λ of the mixture where the velocity drops rapidly is found to correspond to the range of values for the onset of single headed spinning detonations. Thus we may conclude that the onset of single headed spin can be used as a criterion for defining the limits. Spinning detonations are also observed near the limits in annular channels.     

Reaction propagation modes in millimeter-scale tubes for ethylene/oxygen mixtures

Effects of tube diameter and equivalence ratio on reaction front propagations of ethylene/oxygen mixtures in capillary tubes were experimentally analyzed using high speed cinematography. The inner diameters of the tubes investigated were 0.5, 1, 2 and 3 mm. The flame was ignited at the center of the 1.5 m long smooth tube under ambient pressure and temperature before propagated towards the exits in the opposite directions. A total of five reaction propagation scenarios, including deflagration-to-detonation transition followed by steady detonation wave transmission (DDT/C–J detonation), oscillating flame, steady deflagration, galloping detonation and quenching flame, were identified. DDT/C–J detonation mode was observed for all tubes for equivalence ratios in the vicinity of stoichiometry. The velocity for the steady detonation wave propagation was approximately Chapman–Jouguet velocity for 1, 2, and 3 mm I.D. tubes; however, a velocity deficit of 5% was found for the case in 0.5 mm I.D. tube. For leaner mixtures, an oscillating flame mode was found for tubes with diameters of 1 to 3 mm, and the reaction front travelled in a steady deflagrative flame mode with velocities around 2–3 m/s when the mixture equivalence ratio becomes even leaner. Galloping detonation wave propagation was the dominant mode for the fuel lean regime in the 0.5 mm I.D. tube. For rich mixtures beyond the detonation limits, a fast flame followed by flame quenching was observed.     

Investigation of a graded channel InGaAs/GaAs heterostructure transistor

An InGaAS/GaAs heterostructure transistor utilizing a gradedIn_xGa_1 − xAs channel grown by low-pressure metal-olorganic chemical vapor deposition has been demonstrated. A negative differential resistance (NDR) phenomenon is observed. Electron mobilities are significantly improved by using the graded InGaAs channel. For the In composition varying fromx =  0.25 (at the buffer–channel interface) to x =  0.1 (at the spacer–channel interface) structure, a peak extrinsic transconductance of 24.6 S mm ^− 1(atV_DS =  6.5 V,V_GSstep =  − 0.5 mV) and a saturation current density as high as 555 mA mm ^− 1for a gate length of 1.5 μ m are obtained.     

Extinguishment of propane/air co-flowing diffusion flames by fine water droplets

In the present study, extinguishment of propane/air co-flowing diffusion flame by fine water droplets was investigated experimentally. Water droplets are generated by piezoelectric atomizers with the maximum droplets flow rate of 1500 ml/h. When the fuel injection velocity U_f is low, an attached laminar diffusion flame with a premixed flame at the base is stabilized. At some distance from the burner rim, a transition from laminar to turbulent diffusion flame occurs, and a turbulent diffusion flame is formed in the downstream region. When the fuel injector rim is thin (δ = 0.5 mm), the flame stability deteriorates with increase of the co-flowing air stream velocity U_a and the water droplets flow rate Q_m. The stability mechanism can be explained by the balance of the gas velocity and the burning velocity of premixed flame formed at the base. However, when the injector rim is thick (δ = 5 mm), a recirculation zone is produced downstream of the injector rim. The dependence of the quenching distance H_q on U_f and Q_m is relatively weak, and the stability diagram shows curious features. It was shown that U_a is crucially important since it determines flow residence time; if U_a < 0.4 m/s, water droplets can evaporate when they go by the recirculation zone, and the water vapor can diffuse into the recirculation zone. However, if U_a > 0.4 m/s, the water droplets should pass by the recirculation zone without sufficiently evaporated and are not so effective to extinguish the flame. The supply velocity of droplet-laden air should be low enough so that water droplets can evaporate and water vapor can diffuse into the premixed region at the base to obtain sufficient effectiveness of water droplets for fire suppression.     

Derivation and solution of multifrequency radiation diffusion equations for homogeneous refractive lossy media

Starting from the radiation transport equation for homogeneous, refractive lossy media, we derive the corresponding time-dependent multifrequency diffusion equations. Zeroth and first moments of the transport equation couple the energy density, flux and pressure tensor. The system is closed by neglecting the temporal derivative of the flux and replacing the pressure tensor by its diagonal analogue. The radiation equations are coupled to a diffusion equation for the matter temperature. We are interested in modeling heating and cooling of silica (SiO₂), at possibly rapid rates. Hence, in contrast to related work, we retain the temporal derivative of the radiation field. We derive boundary conditions at a planar air–silica interface taking account of reflectivities obtained from the Fresnel relations that include absorption. The spectral dimension is discretized into a finite number of intervals leading to a system of multigroup diffusion equations. Three simulations are presented. One models cooling of a silica slab, initially at 2500 K, for 10 s. The other two are 1D and 2D simulations of irradiating silica with a CO₂ laser, λ = 10.59 μm. In 2D, a laser beam (Gaussian profile, r₀ = 0.5 mm for 1/e decay) shines on a disk (radius = 0.4, thickness = 0.4 cm).     

Microfluidic acoustic trapping force and stiffness measurement using viscous drag effect

It has recently been demonstrated that it was possible to individually trap 70 μm droplets flowing within a 500 μm wide microfluidic channel by a 24 MHz single element piezo-composite focused transducer. In order to further develop this non-invasive approach as a microfluidic particle manipulation tool of high precision, the trapping force needs to be calibrated to a known force, i.e., viscous drag force arising from the fluid flow in the channel. However, few calibration studies based on fluid viscosity have been carried out with focused acoustic beams for moving objects in microfluidic environments.In this paper, the acoustic trapping force (F_trapping) and the trap stiffness (or compliance k) are experimentally determined for a streaming droplet in a microfluidic channel. F_trapping is calibrated to viscous drag force produced from syringe pumps. Chebyshev-windowed chirp coded excitation sequences sweeping the frequency range from 18 MHz to 30 MHz is utilized to drive the transducer, enabling the beam transmission through the channel/fluid interface for interrogating the droplets inside the channel. The minimum force (F_min,trapping) required for initially immobilizing drifting droplets is determined as a function of pulse repetition frequency (PRF), duty factor (DTF), and input voltage amplitude (V_in) to the transducer. At PRF = 0.1 kHz and DTF = 30%, F_min,trapping is increased from 2.2 nN for V_in = 22 V_pp to 3.8 nN for V_in = 54 V_pp. With a fixed V_in = 54 V_pp and DTF = 30%, F_min,trapping can be varied from 3.8 nN at PRF = 0.1 kHz to 6.7 nN at PRF = 0.5 kHz. These findings indicate that both higher driving voltage and more frequent beam transmission yield stronger traps for holding droplets in motion.The stiffness k can be estimated through linear regression by measuring the trapping force (F_trapping) corresponding to the displacement (x) of a droplet from the trap center. By plotting F_trapping – x curves for certain values of V_in (22/38/54 V_pp) at DTF = 10% and PRF = 0.1 kHz, k is measured to be 0.09, 0.14, and 0.20 nN/μm, respectively. With variable PRF from 0.1 to 0.5 kHz at V_in = 54 V_pp, k is increased from 0.20 to 0.42 nN/μm. It is shown that a higher PRF leads to a more compliant trap formation (or a stronger F_trapping) for a given displacement x. Hence the results suggest that this acoustic trapping method has the potential as a noninvasive manipulation tool for individual moving targets in microfluidics by adjusting the transducer’s excitation parameters.     

Comprehensive visualization of detonation-diffraction structures and sizes in unstable and stable mixtures

The geometry and characteristic length of diffraction and re-initiation during a two-dimensional detonation propagation were revealed by visualization. C₂H₄ + 3O₂ (unstable), 2C₂H₂ + 5O₂ + 7Ar (stable) and 2C₂H₂ + 5O₂ + 21Ar (stable) were used as the test mixtures. Experiments were performed over the deviation angle range from 30° to 150° and the initial pressure range from 15.8 to 102.3 kPa. By self-emitting photography, we confirmed that the geometry and the characteristic length of diffraction are not different among test gases, with the exception of the fan-like structure of re-initiation that occurred regardless of whether the mixture was unstable or stable. We conducted a compensative experiment by changing the deviation angle and initial pressure, and summarized the detonation diffraction by shadowgraph. At deviation angles larger than 60°, we measured the distances from the vertex of the channel corner to the point where the transverse detonation wave reflected on the under wall (= wall reflection distance) and confirmed that wall reflection distances are approximately in the range of 10–15 times the cell width, whether the mixture is unstable or stable.     

Dissociation of dimethyl ether at high temperatures

The decomposition of dimethyl ether (CH₃OCH₃) has been investigated behind incident shock waves in a diaphragmless shock tube using laser schlieren densitometry, LS (T = 1500–2450 K, P = 57 ± 4, 125 ± 5 and 253 ± 12 Torr). The LS density gradient profiles were simulated and excellent agreement was found between the simulations and experimental profiles. Rate coefficients for CH₃OCH₃ → CH₃O + CH₃ were obtained. They showed strong fall-off, and at the lower end of the experimental temperature range are close to the low pressure limit. First order rate coefficient expressions were determined over 1500 < T < 2450 K. k_57Torr = (3.10 ± 1.0) × 10⁷⁹T^?19.03 exp(?54417/T) s^?1, k_125Torr = (1.12 ± 0.3) × 10⁸³T^?19.94 exp(?55554/T) s^?1and k_253Torr = (1.02 ± 0.3) × 10⁷³T^?17.09 exp(?51500/T) s^?1. The effect of a roaming channel for decomposition of dimethyl ether was assessed and the best agreement was obtained with 1% dissociation of DME via the roaming path.     

Comparison of characterization methods in high frequency sonochemical reactors of differing configurations

The aim of this study was to compare different characterization methods in order to evaluate the sonochemical efficiency of a cavitational reactor. The selected characterization methods were calorimetry and dosimetry based on potassium iodide oxidation or nitrite and nitrate ion formation. The effects of experimental parameters on physical and chemical effects of ultrasound were quantified with two transducers at a frequency of 366 kHz. The studied factors comprised temperature (16–28 °C), acoustic power (6–38 W), power density (4–61 W L^?1) and reactor configuration (D_{reactor 1} = 65 mm, D_{reactor 2} = 102 mm). Spectrophotometry was compared to ionic chromatography as a method to quantify nitrite and nitrate ions. Spectrometry was shown to be as representative as ionic chromatography. The reaction system based on the formation of both nitrite and nitrate ions was demonstrated to be as reliable as a potassium iodide dosimeter. The representativity of calorimetry was limited since part of acoustic energy was assumed to be used in the chemical reactions observed by dosimetry. Similar sonochemical efficiencies resulted from an increase of sonified surface (D_{reactor 1} = 65 mm vs. D_{reactor 2} = 102 mm) coupled to a 2-time decrease in power density at a constant emitting surface. The effect of emitting-to-sonified surface area ratio on the acoustic field was apparently limited by the height of the liquid.     

Assessment of the link between quantitative biexponential diffusion-weighted imaging and contrast-enhanced MRI in the liver

PurposeTo investigate if intravoxel incoherent motion (IVIM) modeled diffusion-weighted imaging (DWI) can be linked to contrast-enhanced (CE-)MRI in liver parenchyma and liver lesions.MethodsTwenty-five patients underwent IVIM-DWI followed by multiphase CE-MRI using Gd-EOB-DTPA (n = 20) or Gd-DOTA (n = 5) concluded with IVIM-DWI. Diffusion (D_slow), microperfusion (D_fast), its fraction (f_fast), wash-in-rate (R_early) and late-enhancement-rate (R_late) of Gd-EOB-DTPA were calculated voxel-wise for the liver. Parenchyma and lesions were segmented. Pre-contrast IVIM was compared 1) between low, medium and high R_early for parenchyma 2) to post-contrast IVIM substantiated with simulations 3) between low and high R_late per lesion type.ResultsD_fast and f_fast increased (P < 0.001) with 25.6% and 33.8% between low and high R_early of Gd-EOB-DTPA. D_slow decreased (− 15.0%; P < 0.001) with increasing R_early. Gd-DOTA demonstrated similar observations. f_fast (+ 10%; P < 0.001) and D_fast (+ 6.6%; P < 0.001) increased after Gd-EOB-DTPA, while decreasing after Gd-DOTA (− 4.2% and − 5.7%, P < 0.001) and were confirmed by simulations. For focal nodular hyperplasia lesions (n = 5) D_fast and f_fast increased (P < 0.001) with increasing R_late, whereas for hepatocellular carcinoma (n = 4) and adenoma (n = 7) no differences were found.ConclusionMicroperfusion measured by IVIM reflects perfusion in a way resembling CE-MRI. Also IVIM separated intra- and extracellular MR contrast media. This underlines the potential of IVIM in quantitative liver imaging.     

Superconductivity of the platinum doped 122 iron arsenide SrFe2−xPtxAs2

Pt doped 122 iron arsenide SrFe_1?xPt_xAs₂ (0 ? x ? 0.4) was successfully synthesized. The tetragonal unit-cell volume and the lattice constant a increase with increasing the Pt content, while c decreases, suggesting that the Fe ions are indeed replaced by Pt ions. By the Pt doping, the magnetic order of the parent phase is suppressed, and superconductivity emerges at approximately x = 0.15. T_c reaches the maximum of 16 K at x = 0.2. The compounds series can be a suitable subject to investigate role of the doped 5d state in the superconducting 3d Fe–As layer.     

Electrical conductivity of (Ba0.3Sr0.2La0.5)(In1−xFex)O3−δ

We have synthesized the perovskite oxides of the (Ba_0.3Sr_0.2La_0.5)(In_1−xFe_x)O_3−δ system and measured the total electrical conductivity as a function of temperature and oxygen partial pressure. It was found that the single-phase composition region extended from x = 0.0 to x = 1.0, and that the Fe valence increased from 3.06 to 3.50 in that region. The electrical conductivity was semiconducting from x = 0.0 to x = 0.40 and metallic from x = 0.50 to x = 1.0. The total electrical conductivity at 800 °C also increased with the Fe content and achieved a maximum value of 140 (S/cm) at x = 1.0. From the dependence of the electrical conductivity on the oxygen partial pressure, we conclude that above x = 0.50, the majority carriers are holes. The estimated hole conductivity increased exponentially with the amount of Fe⁴⁺ cation present. The oxide ion conductivity was dependent on the oxygen vacancy content.     

Compact WDM demultiplexer for seven channels in photonic crystal

We demonstrate a new device concept for wavelength division demultiplexing based on planar photonic crystal waveguides. The filtering of wavelength channels is realized by shifting the cutoff frequency of the fundamental photonic bandgap mode in consecutive sections of the waveguide. The shift is realized by modifying the size of the border holes.The proposed demultiplexer has an area equal to (16.5 μm × 6.5 μm) and thus it is verified that this structure is very small and can be integrated easily into optical integrated circuits with nanophotonic technologies. The output wavelengths of designed structure can be tuned for communication applications, around 1550 nm. The wavelengths of demultiplexer channels are λ₁ = 1.590 μm, λ₂ = 1.566 μm, λ₃ = 1.525 μm, λ₄ = 1.510 μm, λ₅ = 1.484 μm, λ₆ = 1.450 μm, λ₇ = 1.400 μm respectively. Designs offering improvement of number of the separate wavelengths (seven), miniaturization of the structure (107.25 μm²) is our aim in this work.In our structure, we consider that the 2D triangular lattice photonic crystal is composed of air holes surrounded by dielectric. Its parameters are: radius of holes (r = 0.130 μm), lattice constant (a = 0.380 μm), and index of membrane (n = 3.181:InP). The numerical model used to simulate the structure of the demultiplexer is based on the finite difference time domain (FDTD).     

Crystal growth and structural properties of the spinel-type Li1 + xMn2 − xO4 (x = 0.10, 0.14)

Single crystals of the lithium-rich lithium manganese oxide spinels Li_1 + xMn_2 − xO₄ with x = 0.10 and 0.14 have been successfully synthesized in high-temperature molten chlorides at 1023 K. The single-crystal X-ray diffraction study confirmed the cubic Fd3¯m space group and the lattice parameters of a = 8.2401(9) Å for x = 0.10 and a = 8.2273(10) Å for x = 0.14 at 300 K, respectively. The crystal structures have been refined to the conventional values R = 3.7% for x = 0.10 and R = 3.1% for x = 0.14, respectively. Low-temperature single-crystal X-ray diffraction experiments revealed that these single crystal samples showed no phase transition between 100 and 300 K. The electron-density distribution images in these compounds by the single-crystal MEM analysis clearly showed strong covalent bonding features between the Mn and O atoms due to the Mn–3d and O–2p interaction.     

Excess currents larger than the point contact limit in normal-metal/superconducting junctions

In a point contact NS junction, perfect Andreev reflection occurs over a range of voltages equal to the superconducting energy gap, producing an excess current of I_exc =  (4 / 3)(2 eΔ / h). If the superconductor has a finite width, rather than the infinite width of the point contact, one cannot neglect superfluid flow inside the superconducting contact. The energy range available for perfect Andreev reflections then becomes larger than the superconducting gap, since superfluid flow alters the dispersion relation inside the finite width superconductor. We find a maximum excess current of approximately (7 / 3)(2 eΔ / h) when the width of the superconductor is approximately 7 / 3 times the width of the normal metal.     

Diffusion-promoted-desorption mechanism for D2 desorption from Si(100) surfaces

Temperature-programmed-desorption (TPD) spectra and isothermal desorption rates of D₂ molecules from a Si(100) surface have been calculated to reproduce experimental β_1, A-TPD spectra and isothermal desorption rate curves. In the diffusion-promoted-desorption (DPD) mechanism, hydrogen desorption from the Si(100) (2 × 1) surfaces takes place via D atom diffusion from doubly-occupied Si dimers (DODs) to their adjacent unoccupied Si dimers (UODs). Taking a clustering interaction among DODs into consideration, coverages θ_DU of desorption sites consisting of a pair of a DOD and UOD are evaluated by a Monte Carlo (MC) method. The TPD spectra for the β_1, A peak are obtained by numerically integrating the desorption rate equation R = ν_A exp(? E_d, A / k_BT)θ_DU, where ν_A is the pre-exponential factor and E_d, A is the desorption barrier. The TPD spectra calculated for E_d, A = 1. 6 eV and ν_A = 2.7 × 10⁹ /s are found to be in good agreement with the experimental TPD data for a wide coverage range from 0.01 to 0.74 ML. Namely, the deviation from first-order kinetics observed in the coverage dependent TPD spectra as well as in the isothermal desorption rate curves can be reproduced by the model simulations. This success in reproducing both the experimental TPD data and the very low desorption barrier validates the proposed DPD mechanism.