Experimental investigation of the interaction between weak shock waves and granular layers

 The paper describes new experimental results regarding the pressure fields in front of and inside granular layers of different materials during their collision with weak shock waves. A variety of waves result from the shock wave-granular layer interaction. The pressure behind the reflected wave from the material interface approaches the equilibrium value, P ₅, which would have been reached had the shock wave reflected from a solid end-wall. The wave succession inside the layer depends solely on two processes: the complex interaction of the compaction wave with the granular material and the gas filtration, which affects the particles by the drag forces between the two phases. Inside a material with a permeability coefficient f>0.001 mm² the transmitted wave moves with a constant velocity which is largely governed by the gas filtration. For low permeability materials ( f<0.0003 mm²) the transmitted wave trajectory strongly depends on the compaction wave propagation. In such cases the compaction wave was found to be unsteady and its acceleration was higher in material having low material densities. The maximum compressive stress values, P _c , reached at the shock tube end-wall, covered by the materials under investigation, manifested as an unsteady pressure peak twice as large as the gas pressure P ₅, measured ahead of the layer. Comparing the present data with those available in the literature showed that the amplitude of the unsteady pressure peak was higher in materials having low effective densities, γ, and small permeability coefficients f. Contrary to flexible foams where the available experimental data indicated that the compressive stress in the post peak period converges to P ₅=P _g , the results obtained in the present study indicated that during the test time the compressive stress, P _s , was well preserved in the material and for most of the sample length its value was within the range P _s >P ₅>P _g . Received: 4 March 1996 / Accepted:26 September 1996     

Bifurcations of Heteroclinic Orbits

The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ ₁ as t → −∞ and u(t) → γ ₂ as t → ∞ where γ ₁, γ ₂ are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.     

An investigation of electromagnetic wave propagation in plasma by shock tube

This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters he, v, w, L, wb. The φ800 mm high temperature shock tube has been used to produce a uniform plasma. In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave, the microwave transmission has been used to measure the relative change of the wave power. The working frequency is f = (2-35)GHz (ω=2πf, wave length A =15cm-8mm). The electron density in the plasma is ne = (3&#215;10^10-1&#215;10^14) cm^-3. The collision frequency v = (1&#215;10^8-6&#215;10^10) Hz. The thickness of the plasma layer L = (2-80)cm. The electron circular frequency ωb=eBo/me, magnetic flux density B0 = (0-0.84)T. The experimental results show that when the plasma layer is thick (such as L/λ≥10), the correlation between the attenuation coefficients of the electromagnetic waves and the parameters ne,v,ω, L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma. When the plasma layer is thin (such as when both L and A are of the same order), the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters ne,v,ω, L. In fact, if ω&lt;ωp, v^2&lt;&lt;ω^2, the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions. On the other hand, if ω&gt;ωp, v^2&lt;&lt;ω^2 (just v≈f), the measurements are much larger than the theoretical results. Also, we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field. The result indicates that the value measured under the magnetic field shows a distinct improvement.     

Irradiation-Induced Solute Clustering in a Low Nickel FeMnNi Ferritic Alloy

Understanding the radiation embrittlement of reactor pressure vessel (RPV) steels is required to be able to operate safely a nuclear power plant or to extend its lifetime. The mechanical properties degradation is partly due to the clustering of solute under irradiation. To gain knowledge about the clustering process, a Fe−1.1 Mn−0.7 Ni (at.%) alloy was irradiated in a test reactor at two fluxes of 0.15 and 9 ×10¹⁷ n _E > 1MeV .m^− 2.s^− 1 and at increasing doses from 0.18 to 1.3 ×10²⁴ n _E > 1MeV .m^− 2 at 300°C. Atom probe tomography (APT) experiments revealed that the irradiation promotes the formation in the α iron matrix of Mn/Mn and/or Ni/Ni pair correlations at low dose and Mn–Ni enriched clusters at high dose. These clusters dissolve partially after a thermal treatment at 400°C. Based on a comparison with thermodynamic calculations, we show that the solute clustering under irradiation can just result from an induced mechanism.     

Regularity of the Inverse of a Planar Sobolev Homeomorphism

Let be a domain. Suppose that f ∈ W^1,1_loc(Ω,R²) is a homeomorphism such that Df(x) vanishes almost everywhere in the zero set of J_f. We show that f^-1 ∈ W^1,1_loc(f(Ω),R²) and that Df⁻¹(y) vanishes almost everywhere in the zero set of Sharp conditions to quarantee that f⁻¹ ∈ W^1,q(f(Ω),R²) for some 1<q≤2 are also given.     

Reynolds' analogy for the thermal convection driven by nonisothermal stretching surfaces

In the present paper the steady boundary-layer flows induced by permeable stretching surfaces with variable temperature distribution are investigated under the aspect of Reynolds' analogy r = St _x /C _f(x). It is shown that for certain stretching velocities and wall temperature distributions, “Reynolds' function”r, i.e. the ratio of the local Stanton number St _x and the skin friction coefficient C _f(x) equals −1/2 for any value of the Prandtl number Pr and of the dimensionless suction/injection velocity f _w. In all of these cases, the dimensionless temperature field ϑ is connected to the dimensionless downstream velocity f ^′ by the simple relationship ϑ=(f ^′)^Pr. It is also shown that in the general case, Reynolds' function r may possess several singularities in f _w. The largest of them represents a critical value, so that for f _w<f _w,crit the solutions of the energy equation (although they still satisfy all the boundary conditions) become nonphysical.     

Precursor ionization ahead of strong shock waves in argon

The mechanism of precursor ionization ahead of strong shock waves has been studied in a low density shock tube. The experimental results are illustrated with Arrhenius plots with kink points dividing them into two parts with apparent activation energy ratio 1:2, namely with the values 7.7 eV and 15.3 eV, and varying with first and third power of the density respectively. A model is proposed to interpret the facts where the process taking place in the precursor region, is a two step photo ionization accompanied with the drift flow effect of the gas relative to the shock wave or the ionization recombination effect according to whether the shock speed and initial density are low enough. The product of the A-A collision excitation cross section coefficientS ^* multiplied by the radiation cross sectionQ ^* of ArgonS ^*×Q ^*=1×10⁻³⁶ (cm⁴eV⁻¹) and the three body recombination coefficient of Argon at room temperaturek _ra =1×10⁻²⁴ (cm⁻⁶s⁻¹). The project supported by the National Natural Science Foundation of China     

A pulse-expansion wave tube for nucleation studies at high pressures

 The design and performance of a new pulse-expansion wave tube for nucleation studies at high pressures are described. The pulse-expansion wave tube is a special shock tube in which a nucleation pulse is formed at the endwall of the high pressure section. The nucleation pulse is due to reflections of the initial shock wave at a local widening situated in the low pressure section at a short distance from the diaphragm. The nucleation pulse has a duration of the order of 200 μs, while nucleation pressures that can be achieved range from 1 to 50 bar total pressure. Droplet size and droplet number density can accurately be determined by a 90°-Mie light scattering method and a light extinction method. The range of nucleation rates that can be measured is 10⁸ cm^-3 s^-1<J<10¹¹ cm^-3 s^-1. We will illustrate the functioning and possibilities of the new pulse-expansion wave tube by nucleation rate measurements in the gas-vapour mixture nitrogen/water in the temperature range 200–260 K, and in the mixture methane/n-nonane as a function of supersaturation S at various total pressures up to 40 bar and temperatures around 240 K. Received: 5 June 1996/Accepted: 9 December 1996     

Shear layers of a circular cylinder with rotary oscillation

The behavior of the separated shear layers and the near wake of a circular cylinder with small-amplitude rotary oscillations (Ω₁ = 0.05−0.15 for f _f/f _o ≤ 1.25) were investigated experimentally at Re = 3,700. Measurements of an unforced cylinder were also made for 2,000 ≤ Re ≤ 10,000 to better understand the effects of rotary oscillations. The results show that the shear-layer vortices formed closer to the cylinder and the distance separating them was found to decrease with cylinder oscillations. The shear-layer frequency, however, increased with increasing forcing frequency f _f. The formation-region length l _f decreased significantly with increasing f _f while decreased to a lesser extent with increasing normalized oscillation amplitude Ω₁. The shear layer also diffused to a length L _d larger than that of an unforced cylinder, while the l _f-L _d-Strouhal frequency offsetting mechanism was generally maintained. The near wake was of lower momentum compared to an unforced cylinder, and the transverse velocity fluctuations associated with the unforced vortex-shedding frequency f _o always presented a local peak at f _f/f _o = 0.5, regardless of Ω₁ tested.     

Miniaturized Compression Test at Very High Strain Rates by Direct Impact

A modified miniaturized version of the Direct Impact Compression Test (DICT) technique is described in this paper. The method permits determination of the rate-sensitive plastic properties of materials up to strain rate ∼10⁵ s⁻¹. Miniaturization of the experimental setup with specimen dimensions: diameter d _S = 2.0 mm and thickness l _S = 1.0 mm, Hopkinson bar diameter 5.2 mm, with application of a novel optical arrangement in measurement of specimen strain, makes possible compression tests at strain rates from ∼10³ s⁻¹ to ∼10⁵ s⁻¹. In order to estimate the rate sensitivity of a low-alloy construction steel, quasi-static, Split Hopkinson Pressure Bar (SHPB) and DICT tests have been performed at room temperature within the rate spectrum ranging from 5*10⁻⁴ s⁻¹ to 5*10⁴ s⁻¹. Adiabatic heating and friction effects are analyzed and the final true stress versus true strain curves at different strain rates are corrected to a constant temperature and zero friction. The results have been analyzed in the form of true stress versus the logarithm of strain rate and they show two regions of a constant rate sensitivity : relatively low up to the strain rate threshold ∼50 s⁻¹, and relatively high above the threshold, up to strain rate ∼4.5*10⁴ s⁻¹.     

Investigation of liquid–liquid drop coalescence using tomographic PIV

High-speed tomographic PIV was used to investigate the coalescence of drops placed on a liquid/liquid interface; the coalescence of a single drop and of a drop in the presence of an adjacent drop (side-by-side drops) was investigated. The viscosity ratio between the drop and surrounding fluids was 0.14, the Ohnesorge number (Oh = μ_d/(ρ_dσD)^1/2) was 0.011, and Bond numbers (Bo = (ρ _d  − ρ _s )gD ²/σ) were 3.1–7.5. Evolving volumetric velocity fields of the full coalescence process allowed for quantification of the velocity scales occurring over different time scales. For both single and side-by-side drops, the coalescence initiates with an off-axis film rupture and film retraction speeds an order of magnitude larger than the collapse speed of the drop fluid. This is followed by the formation and propagation of an outward surface wave along the coalescing interface with wavelength of approximately 2D. For side-by-side drops, the collapse of the first drop is asymmetric due to the presence of the second drop and associated interface deformation. Overall, tomographic PIV provides insight into the flow physics and inherent three-dimensionalities in the coalescence process that would not be achievable with flow visualization or planar PIV only.     

Estimating the viscoelastic moduli of complex fluids using the generalized Stokes–Einstein equation

We obtain the linear viscoelastic shear moduli of complex fluids from the time-dependent mean square displacement, <Δr ²(t)>, of thermally-driven colloidal spheres suspended in the fluid using a generalized Stokes–Einstein (GSE) equation. Different representations of the GSE equation can be used to obtain the viscoelastic spectrum, G˜(s), in the Laplace frequency domain, the complex shear modulus, G ^*(ω), in the Fourier frequency domain, and the stress relaxation modulus, G _r (t), in the time domain. Because trapezoid integration (s domain) or the Fast Fourier Transform (ω domain) of <Δr ²(t)> known only over a finite temporal interval can lead to errors which result in unphysical behavior of the moduli near the frequency extremes, we estimate the transforms algebraically by describing <Δr ²(t)> as a local power law. If the logarithmic slope of <Δr ²(t)> can be accurately determined, these estimates generally perform well at the frequency extremes. Received: 8 September 2000/Accepted: 9 March 2000     

Radial filtration during constant-force squeeze flow of soft solids

Various soft solid suspensions were squeezed at constant force between polished and roughened circular glass plates and the time-dependence of the interplate separation was measured. The filterability of suspensions was quantified by their desorptivity S obtained from measurements of capillary suction time. The squeeze flow (SF) of suspensions for which S < 2 μm s^−1/2 was largely consistent with rheological theory, which neglected radial filtration: the relative motion between the liquid and solid phases of the suspension in the radial pressure gradient. Suspensions having S > 2 μm s^−1/2 showed SF behavior that was consistent with the presence of radial filtration.     

On the completeness of oscillation spaces

The oscillation spaces introduced by Jaffard are a variation on the definition of Besov spaces for either s ≥ 0 or s ≤ −d/p. On the contrary, the spaces for −d/p < s < 0 cannot be sharply imbedded between Besov spaces with almost the same exponents, and, thus, they are new spaces of really different nature. Their norms take into account correlations between the positions of large wavelet coefficients through the scales. Several numerical studies uncovered such correlations in several settings including turbulence, image processing, traffic, finance, etc. These spaces allow one to capture oscillatory behaviors that are left undetected by Sobolev or Besov spaces. Unlike Sobolev spaces (respectively, Besov spaces B _p ^s,q (ℝ^d)), which are expressed by simple conditions on wavelet coefficients as ℓ^p norms (respectively, mixed ℓ^p − ℓ^q norms), oscillation spaces are written as ℓ^p averages of local C ^s ′ norms. In this paper, we prove the completeness of oscillation spaces in spite of such a mixture of two norms of different kinds. Published in Neliniini Kolyvannya, Vol. 8, No. 4, pp. 435–443, October–December, 2005.     

ON THE CRITICAL POINTS OF THE MAP Fp:X→‖AXB-C‖pp

Suppose A,B and C are the bounded linear operators on a Hilbert space H, when A has a generalized inverse A^- such that (AA^-)^*=AA^- and B has a generalized inverse B^- such that (B^-B)^*=B^-B,the general characteristic forms for the critical points of the map F_p:X^→‖AXB-C‖_p^p(1p=2. Similarly, the same question has been discussed for several operators.     

Experiments on the lift of a spinning sphere in a range of intermediate Reynolds numbers

 The lift force experienced by a spinning sphere moving in a viscous fluid, with constant linear and angular velocities, is measured by means of a trajectographic technique. Measurements are performed in the range of dimensionless angular velocities γ=aω/V lying between 1 and 6, and in the range of Reynolds numbers Re=2aV/ν lying between 10 and 140 (a sphere radius, ω angular velocity, V relative velocity of the sphere centre, ν fluid kinematic viscosity). A notable departure from the theoretical relationship at low Reynolds number, C _L =2γ, is obtained, the ratio C _L /γ being found to significantly decrease with increasing γ and increasing Re. The following correlation is finally proposed to estimate the lift coefficient in the range 10<Re<140: C _L ≅0.45+(2γ−0.45) exp (−0.075γ^0.4 Re ^0.7) Received: 20 May 1996/Accepted: 9 November 1997     

Stability of Degenerate Stationary Waves for Viscous Gases

This paper is concerned with the asymptotic stability of degenerate stationary waves for viscous gases in the half space. We discuss the following two cases: (1) viscous conservation laws and (2) damped wave equations with nonlinear convection. In each case, we prove that the solution converges to the corresponding degenerate stationary wave at the rate t ^−α/4 as t → ∞, provided that the initial perturbation is in the weighted space L²_a=L²(\mathbb R₊; (1+x)^a dx){L^2_\alpha=L^2({\mathbb R}_+;\,(1+x)^\alpha dx)} . This convergence rate t ^−α/4 is weaker than the one for the non-degenerate case and requires the restriction α < α_*(q), where α_*(q) is the critical value depending only on the degeneracy exponent q. Such a restriction is reasonable because the corresponding linearized operator for viscous conservation laws cannot be dissipative in L²_a{L^2_\alpha} for α > α^*(q) with another critical value α^*(q). Our stability analysis is based on the space–time weighted energy method in which the spatial weight is chosen as a function of the degenerate stationary wave.     

On the generalized Cauchy function and new Conjecture on its exterior singularities

This article studies on Cauchy’s function f (z) and its integral, (2πi)J[ f (z)] ≡■C f (t)dt/(t z) taken along a closed simple contour C, in regard to their comprehensive properties over the entire z = x + iy plane consisted of the simply connected open domain D + bounded by C and the open domain D outside C. (1) With f (z) assumed to be C n (n < ∞-times continuously differentiable) z ∈ D + and in a neighborhood of C, f (z) and its derivatives f (n) (z) are proved uniformly continuous in the closed domain D + = [D + + C]. (2) Cauchy’s integral formulas and their derivatives z ∈ D + (or z ∈ D ) are proved to converge uniformly in D + (or in D = [D +C]), respectively, thereby rendering the integral formulas valid over the entire z-plane. (3) The same claims (as for f (z) and J[ f (z)]) are shown extended to hold for the complement function F(z), defined to be C n z ∈ D and about C. (4) The uniform convergence theorems for f (z) and F(z) shown for arbitrary contour C are adapted to find special domains in the upper or lower half z-planes and those inside and outside the unit circle |z| = 1 such that the four general- ized Hilbert-type integral transforms are proved. (5) Further, the singularity distribution of f (z) in D is elucidated by considering the direct problem exemplified with several typ- ical singularities prescribed in D . (6) A comparative study is made between generalized integral formulas and Plemelj’s formulas on their differing basic properties. (7) Physical sig- nificances of these formulas are illustrated with applicationsto nonlinear airfoil theory. (8) Finally, an unsolved inverse problem to determine all the singularities of Cauchy function f (z) in domain D , based on the continuous numerical value of f (z) z ∈ D + = [D + + C], is presented for resolution as a conjecture.     

Characterization of shock accelerometers using Davies bar and laser interferometer

In this paper, we propose a novel method for evaluating the frequency response of shock accelerometers using Davies bar and interferometry. The method adopts elastic wave pulses propagating in a thin circular bar for the generation of high accelerations. The accelerometer to be examined is attached to one end of the bar and experiences high accelerations of the order of 10³∼10⁵ m/s². A laser interferometer system is newly designed for the absolute measurement of the bar end motion. It can measure the motion of a diffuse surface specimen at a speed of 10⁻³ ∼10⁰ m/s. Uncertainty of the velocity measurement is estimated to be±6×10⁻⁴ m/s, proving a high potential for use in the primary calibration of shock accelerometers. Frequency characteristics of the accelerometer are determined by comparing the accelerometer's output with velocity data of the interferometry in the frequency domain. Two piezoelectric-type accelerometers are tested in the experiment, and their frequency characteristics are obtained over a wide frequency range up to several ten kilohertz. It is also shown that the results obtained using strain gages are consistent with those by this new method. Paper was presented at the 1992 SEM Spring Conference on Experimental Mechanics held in Las Vegas, NV on June 8–11.     

Study of airborne particles produced by normal impact of millimetric droplets onto a liquid film

The aim of this work is to carry out an experimental investigation into the generation of airborne microparticles when millimetric droplets of aqueous solutions impact onto a liquid film. Impact experiments using 3.9 mm diameter droplets were carried out for Weber numbers between 159 and 808, with a fixed Ohnesorge number of 2 × 10⁻³ and film parameters S _f (the ratio between the thickness of the liquid film h _film and the diameter of the impacting droplet d _i) between 0.3 and 1. Observed results show that the deposition/splashing threshold is independent of the parameter S _f in agreement with the data in the literature. The aerosol measurement results demonstrate the production of solid particles from the evaporation of secondary microdroplets with diameters less than 30 μm formed when splash occurs. The median diameter of these microdroplets is around 20 μm, corresponding to a value of d ₅₀/d _i = 5 × 10⁻³. Taken together, the results show that the mass and the number of particles emitted increase as the Weber number increases. Moreover, at a Weber number of 808, the results show that the mass and number of particles emitted increases as the parameter S _f decreases. In this case, the mean number of microdroplets emitted per impact is equal to 14 for S _f = 1 and equal to 76 for S _f = 0.3.