Shock velocity measurements using microwave Doppler shift

A commercially available 10.587 GHz microwave Doppler module is used for the measurement of shock velocity in a conventional shock tube. With proper electronic circuits the Doppler frequency obtained is found to be quite noisefree and consistent for shock velocities in the range of 1.8 mm/sec to 2.0 mm/sec.     

A multiplicative ergodic theorem for rotation numbers

Given a vector fieldX on a Riemannian manifoldM of dimension at least 2 whose flow leaves a probability measure invariant, the multiplicative ergodic theorem tells us that-a.s. every tangent vector possesses a Lyapunov exponent (exponential growth rate) that is equal to one of finitely many basic exponents corresponding toX and. We prove that, in the case of a simple Lyapunov spectrum, every tangent plane-a.s. possesses a rotation number that is equal to one of finitely many basic rotation numbers corresponding toX and. Rotation in a plane is measured as the time average of the infinitesimal changes of the angle between a frame moved by the linearized flow and the same frame parallel-transported by a (canonical) connection.     

Feasability study of wall shear stress imaging using microstructured surfaces with flexible micropillars

A new optical sensor technique based on a sensor film with arrays of hair-like flexible micropillars on the surface is presented to measure the temporal and spatial wall shear stress field in boundary layer flows. The sensor principle uses the pillar tip deflection in the viscous sublayer as a direct measure of the wall shear stress. The pillar images are recorded simultaneously as a grid of small bright spots by high-speed imaging of the illuminated sensor film. Two different ways of illumination were tested, one of which uses the fact that the transparent pillars act as optical microfibres, which guide the light to the pillar tips. The other method uses pillar tips which were reflective coated. The tip displacement field of the pillars is measured by image processing with subpixel accuracy. With a typical displacement resolution on the order of 0.2 m, the minimum resolvable wall friction value is _w20 mPa. With smaller pillar structures than those used in this study, one can expect even smaller resolution limits.     

Tensile testing of polysilicon

Tensile specimens of polysilicon are deposited on a silicon wafer; one end remains affixed to the wafer and the other end has a relatively large paddle that can be gripped by an electrostatic probe. The overall length of the specimen is less than 2 mm, but the smooth tensile portion can be as small as 1.5×2m in cross section and 50m long. The specimen is pulled by a computer-controlled translation stage. Force is recorded with a 100-g load cell, whereas displacement is recorded with a capacitance-based transducer. Strain can be measured directly on wider specimens with laser-based interferometry from two small gold markers deposited on the smooth portion of the specimen. The strength of this linear and brittle material is measured with relative ease. Young's modulus measurement is more difficult; it can be determined from either the stress-strain curve, the record of force versus displacement or the comparison of the records of two specimens of different sizes. Specimens of different sizes—thicknesses of 1.5 or 3.5 m, widths from 2 to 50 m and lengths from 50 to 500 m—were tested. The average tensile strength of this polysilicon is 1.45±0.19 GPa (210 ±28 ksi) for the 27 specimens that could be broken with electrostatic gripping. The average Young's modulus from force displacement records of 43 specimens is 162±14 GPa (23.5 ±2.0×10³ ksi). This single value is misleading because the modulus values tend to increase with decreasing specimen width; that is not the case for the strength. The three methods for determining the modulus agree in general, although the scatter can be large.     

The effect of background particulates on the delayed transition of a heated 9:1 ellipsoid

An experimental study was done to quantify the effects of a variety of background particulates on the delayed laminar-turbulent transition of a thermally stabilized boundary layer in water. A Laser-Doppler Velocimeter system was used to measure the location of boundary layer transition on a 50 mm diameter, 9:1 fineness ratio ellipsoid. The ellipsoid had a 0.15 m RMS surface finish. Boundary layer transition locations were determined for length Reynolds numbers ranging from 3.0 × 10⁶ to 7.5 × 106. The ellipsoid was tested in three different heating conditions in water seeded with particles of four distinct size ranges. For each level of boundary layer heating, measurements of transition were made for clean water and subsequently, water seeded with 12.5 m, 38.9 m, 85.5 m and 123.2 m particles, alternately. The three surface heating conditions tested were no heating, T = 10°C and T = 15°C where T is the difference between the inlet model heating water temperature, T _i, and free stream water temperature, T . The effects of particle concentration were studied for 85.5 m and 123.2 m particulates.The results of the study can be summarized as follows. The 12.5 m and 38.9 m particles has no measurable effect on transition for any of the test conditions. However, transition was significantly affected by the 85.5 m and 123.2 m particles. Above a length Reynolds number of 4 × 10⁶ the boundary layer transition location moved forward on the body due to the effect of the 85.5 m particles for all heating conditions. The largest percentage changes in transition location from clean water, were observed for 85.5 m particles seeded water.Transition measurements made with varied concentrations of background particulates indicated that the effect of the 85.5 m particles on the transition of the model reached a plateau between 2.65 particulates/ml concentration and 4.2 particles/ml. Measurements made with 123.3 m particles at concentrations up to 0.3 part/ml indicated no similar plateau.     

An experimental method for the investigation of droplet oscillations in a gaseous medium

An experimental method for the investigation of droplet oscillations in a gaseous medium is presented. The droplets are produced using vibrating orifice droplet generators. Experiments are carried out with droplets in the diameter range from 91 m to 288 m using propanol-2, water, and n-hexadecane; the gaseous host medium is air. Oscillatory motions of the fundamental mode n = 2 and of the first higher order mode n = 3 occur during the disintegration of the liquid jet produced by the droplet generator. The periodical production of the droplets allows the observation and evaluation of each phase of the motion under quasi-steady conditions. Surface energies are determined from the droplet shapes on photos. The periods of the oscillations are found to be very close to the prediction of the linear theory.     

Some problems of nonisothermal steady flow of a viscous fluid

The paper presents solutions to the problems of plane Couette flow, axial flow in an annulus between two infinite cylinders, and flow between two rotating cylinders. Taking into account energy dissipation and the temperature dependence of viscosity, as given by Reynolds's relation =₀ exp (–T) (₀, =const). Two types of boundary conditions are considered: a) the two surfaces are held at constant (but in general not equal) temperatures; b) one surface is held at a constant temperature, the other surface is insulated.Nonisothermal steady flow in simple conduits with dissipation of energy and temperature-dependent viscosity has been studied by several authors [1–11]. In most of these papers [1–6] viscosity was assumed to be a hyperbolic function of temperature, viz. =_m 1/1+²(T–T_m.Under this assumption the energy equation is linear in temperature and can he easily integrated. Couette flow with an exponential viscosity-temperature relation. =₀ e ^–T (₀, =const), (0.1) was studied in [7, 8]. Couette flow with a general (T) relation was studied in (9).Forced flow in a plane conduit and in a circular tube with a general (T) relation was studied in [10]. In particular, it has been shown in [10] that in the case of sufficiently strong dependence of viscosity on temperature there can exist a critical value of the pressure gradient, such that a steady flow is possible only for pressure gradients below this critical value.In a previous work [11] the authors studied Polseuille flow in a circular tube with an exponential (T) relation. This thermohydrodynamic problem was reduced to the problem of a thermal explosion in a cylindrical domain, which led to the existence of a critical regime. The critical conditions for the hydrodynamic thermal explosion and the temperature and velocity profiles were calculated.In this paper we treat the problems of Couette flow, pressureless axial flow in an annulus, and flow between two rotating cylinders taking into account dissipation and the variation of viscosity with temperature according to Reynolds's law (0.1). The treatment of the Couette flow problem differs from that given in [8] in that the constants of integration are found by elementary methods, whereas in [8] this step involved considerable difficulties. The solution to the two other problems is then based on the Couette problem.     

Measurement of the Coefficient of Transverse Dispersion in Flow Through Packed Beds for a Wide Range of Values of the Schmidt Number

Experimental values of the coefficient of transverse dispersion (D _T) were measured with the system 2-naphthol/water, over a range of temperatures between 293K and 373K, which corresponds to a range of values of viscosity () between 2.83×10^–4 Ns/m² and 1.01×10^–3 Ns/m² and of molecular diffusion coefficient (D _m) between 1.03×10^–9 m²/s and 5.49×10^–9 m²/s. Since the density () of water is close to 10³ kg/m³, the corresponding variation of the Schmidt number (Sc=/D _m) was in the range 1000 – 50. More than 200 experimental values of the transverse dispersion coefficient were obtained using beds of silica sand with average particle sizes (d) of 0.297 and 0.496mm, operated over a range of interstitial liquid velocities (u) between 0.1mm/s and 14mm/s and this gave a variation of the Reynolds number (Re=du/) between 0.01 and 3.5.Plots of the dimensionless coefficient of transverse dispersion (D _T/D _m) vs. the Peclet number (Pe_m=ud/D _m) based on molecular diffusion bring into evidence the influence of Sc on transverse dispersion. As the temperature is increased, the value of Sc decreases and the values of D _T/D _m gradually approach the line corresponding to gas behaviour (i.e. Sc 1), which is known to be well approximated by the equation D _T/D _m=1/+ud/12D _m, where is the tortuosity with regard to diffusion.     

Digital image plane holography (DIPH) for two-phase flow diagnostics in multiple planes

A technique for measuring the size and displacement of the disperse phase in two planes of a two-phase flow is presented. Digital image plane holography (DIPH) is used for the simultaneous recording and independent reconstruction of both planes. Each fluid plane is illuminated with two laser sheets propagating in opposite directions. The defocused image fields are holographically recorded at 90°, and can be reconstructed either in a defocused or in the best-focused plane. The analysis of the images in a defocused plane provides the sizes, while the cross-correlation of the focused images provides the velocity field, as in a regular particle image velocimetry (PIV) experiment. For air bubbles freely drifting in glycerine, diameters from 50 m to 400 m and displacements of up to 300 m have been measured.     

Asymmetric snap-buckling of a column restrained by a stiff wire

Summary A relatively simple example of asymmetric snap-through buckling in a continuous structure is the nonlinear problem of a cantilevered column restrained at its tip by a stiff wire, which is inclined at an acute angle to the column centerline, and loaded at its tip by a force perpendicular to the centerline. A parameter called , which is the nondimensional ratio of the flexural rigidity of the column to the combined extensional stiffness of the wire and the column centerline, determines the essential features of the buckling. If is zero, or is small compared to unity, the bending of the column is small enough to justify the use of linear bending theory for the column. Hence, even though the constraint is nonlinear, the solution to this problem is obtained in closed form. The critical point for the structure is found to be an asymmetric branching point for =0, while for positive, the critical point is a snap-through type. The effect of is similar to that induced by initial imperfections in more complex structures. For very small , the critical load is markedly decreased from the value for =0. Moreover, the graph of the load vs. tip deflection has the appearance of having an acute discontinuity in slope at the critical point for very small, although it is actually found that the graph has a horizontal tangent there.

---

Sommario Un esempio relativamente semplice di collasso asimmetrico per carico di punta in una struttura continua è il problema non lineare di una colonna incastrata ad una estremità e vincolata all'altra estremità da un filo rigido che forma un angolo acuto con l'asse della colonna e caricata a quella estremità con una forza perpendicolare allo stesso asse. Un parametro chiamato , che è il rapporto adimensionale fra la rigidità flessionale della colonna e la rigidezza longitudinale del filo e dell'asse della colonna, determina le caratteristiche essenziali del cedimento. Se è zero o è piccolo rispetto all'unità, l'inflessione della colonna è sufficientemente piccola per giustificare l'uso della teoria lineare di inflessione per la colonna. Di conseguenza anche se il vincolo non è lineare la soluzione del problema è ottenuta in forma compatta. Il punto critico della struttura si trova nel punto di biforcazione asimmetrica per =0, mentre per positivo, il punto critico rappresenta un punto di collasso. L'effetto di è simile a quello prodotto da imperfezioni iniziali in strutture più complesse. Per molto piccolo il carico critico è notevolmente ridotto rispetto al valore per =0. Inoltre il grafico del carico in funzione della curvatura all'estremità sembra avere una netta discontinuità nella pendenza dal punto di biforcazione per molto piccolo benchè, in realtà si trovi che il grafico ha lì una tangente orizzontale.

---

This research was supported by the Advanced Research Projects Agency (Project DEFENDER) and was monitored by the U.S. Army Research Office, Durham, under Contract DA-31-124-ARO-D-257.     

Turbulent dispersion of droplets for air flow in a pipe

An optical technique was used to study the dispersion of 50 m, 90 m and 150 m droplets downstream of a source located in the center of a vertical pipe through which turbulent air is flowing. A turbulent dispersion coefficient, _P, and a mean-square of the fluctuations in the turbulent velocity, v _p ² , are determined from the change of the measured mean-square displacement of the droplets over the pipe cross section with time. The interesting aspect of the experiments is that they explored conditions where the inertia of the particles is believed to be a much more important effect than that of the crossing of trajectories associated with the inequality of the average velocities of the particles and the fluid.     

Temperature measurement near the triple line during phase change using thermochromic liquid crystal thermography

The temperature generated by the evaporation of a volatile liquid in a confined space (tube =1,630 m) was mapped onto the tube surface with the use of unsealed thermochromic liquid crystals (TLCs). The strong evaporative cooling effect located near the meniscus triple line generates a temperature dip. Despite the thermal diffusion through the tubes thickness and its geometry, the TLC thickness and the inherent difficulties of working with unsealed TLCs, the present technique has revealed to be a suitable tool for accurate temperature measurement at the microscale size. The evaporation flux is deduced from the profile of temperature and comparison with the experimental measurement shows a very good agreement. The role of the nature and thickness of the tube wall material on the diffusion of the temperature profile from the inside to the outside is also investigated.     

Die Temperaturabhängigkeit der Oberflächenspannung von reinen Kältemitteln vom Tripelpunkt bis zum kritischen Punkt

Zusammenfassung Die Oberflächenspannung von sechs reinen Substanzen — SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃ und CHClF₂ — wurde mit Hilfe einer modifizierten Kapillarmethode gemessen. Die zur Berechnung der Oberflächenspannung erforderlichen Sättigungsdichten und wurden teils aus vorhandenen Zustandsgleichungen, teils aus ebenfalls gemessenen Brechungsindizes bestimmt. Die Temperaturabhängigkeit der Oberflächenspannung läßt sich durch einen erweiterten Ansatz nach van der Waals =_O (T_c-T)(1+...) darstellen, wobei bei einfachen Stoffen ein eingliedriger, bei polaren und assoziierenden Stoffen ein zweigliedriger Ansatz notwendig und ausreichend ist. Für den kritischen Exponenten der Oberflächenspannung wurde ein von der molekularen Substanz weitgehend unabhängiger Wert von =1.284±0.005 gefunden.

---

Temperature dependence of surface tension of pure refrigerants from triple point up to the critical point

The surface tension of six fluids (SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃, CHClF₂) have been measured by means of a modified capillary rise method. The liquid vapor densities, which are needed to calculate the surface tension, have partly been determined by means of refractive indices simultaneously measured in the same apparatus. The temperature dependence of the surface tension is described by an extended van der Waals power law =_O(T_c-T)(1+...). For simple fluids one term and for polar and associating fluids two terms are necessary and sufficient. The critical exponent is found to be 1.284 ± 0.005 and nearly independent of the molecular structure.

---

Formelzeichen a² Laplace-Koeffizient - a Parameter - B_O, B_on Koeffizient der Koexistenzkurve - g Erdbeschleunigung - H Höhe, kapillare Steighöhe - LL Lorentz-Lorenz-Funktion oder Refraktionskonstante - M molare Masse - M Zahl der Meßwerte - N Zahl der unbekannten Parameter - n Brechungsindex - p Druck - R,r Radius - s Entropie - SD Standardabweichung - T, t Temperatur - u innere Energie Griechische Formelzeichen Exponent des Laplace-Koeffizienten - Exponent der Koexistenzkurve - 2. Exponent der Oberflächenspannung - Wellenlänge des Lichts - Exponent der Oberflächenspannung - _D Dipolmoment - , Dichte der Flüssigkeit bzw. des Dampfes - Oberflächenspannung - reduzierte Temperatur (1-T/T_c) - ² gewichtete Varianz Indizes c kritischer Zustand - D Differenz - m Mittelwert - T Isotherme - t Zustand am Tripelpunkt - S Zustand am Schmelzpunkt - bezogen auf Oberfläche     

Bifurcations and complex instability in a 4-dimensional symplectic mapping

Summary We study the main periodic solutions of a 4-dimensional symplectic mapping composed of two coupled 2-dimensional mappings. Their bifurcations were calculated numerically and also theoretically for small values of the coupling parameter . Most bifurcating families of period 2ⁿ (n0) have complex unstable regions that extend from =0 to the maximum allowed value of for each family. These complex unstable regions do not allow the transmisssion of the stability of the corresponding families to families of higher order. We found only one family with a complex unstable region not extending to the maximum , but in this case also the transmission of the stability is stopped at an inverse bifurcation. Thus although there are infinite sequences of bifurcations (of the Feigenbaum type) in the limiting 2-dimensional case =0, all such sequences are interrupted when the system is 4-dimensional (i.e. for 0). The appearance of complex instability for =0 can be predicted by studying the cases =0 and applying the Krein-Moser theorem.

---

Sommario Si svolge uuno studio dettagliato delle soluzioni periodiche principali di due mappe simplettiche bidimensionali accoppiate, calcolandone sia analiticamente che numericamente le biforcazioni per piccoli valori del parametro di accoppiamento . Quasi tutte le famiglie di periodo 2ⁿ (n0) prodotte dalle biforcazioni presentano regioni di instabilità complessa che si estendono da =0 fino al massimo valore di considerato. Queste regioni di instabilità complessa impediscono il trasferimento della stabilità di una famiglia a famiglie di ordine più elevato. In un solo caso si osserva una famiglia la cui regione di instabilità complessa non arriva ad estendersi fino al valore massimo di ; in questo caso però il trasferimento della stabilità viene interrotto da una biforcazione inversa. Se ne conclude che, nonostante I'esistenza di una famiglia di infinite biforcazioni di tipo Feigenbaum nel caso limite bidimensionale (=0), tutte le sequenze si interrompono se il sistema è a quattro dimensioni. Il formarsi di regioni di instabilità complessa per 0 può essere previsto studiando il caso =0 ed appplicando il teorema di Krein-Moser.

     

On the evolutions of the discontinuities of orders ⩾ m in solutions for quasi-linear systems of PDEs of order m,or for the linearizations of these systems

Summary One considers a system L[u]=0 of PDEs, quasi-linear (according to [1]) and of order m, which possesses a bicharacteristic line , as it happens in the hyperbolic case. For v=0, , –m (>0) let u(^v) be a discontinuity wave of order m+v that solves the system above and whose discontinuity hypersurface includes . The corresponding transport equations along are considered. Furthermore some interesting cases are pointed out, in which these equations turn out to be mutually equivalent in a suitable sense. Some theorems are stated to compare the transport equations for the discontinuities of the above kinds, that are connected with the systems d^hL[u]/dt^h=0 (h=0, , –m) and/or the linearization of the system L[u]=0 around any regular solution of it.

---

Sommario Si considera un sistema L[u]=0 di equazioni alle derivate parziali, quasi lineare (secondo [1]) e di ordine m, il quale sia dotato di qualche bicaratteristica , come accade nel caso iperbolico. Per v=0, , –m(>0) sia u(^v) un'onda di discontinuità di ordine m+v risolvente il detto sistema e avente ipersuperficie di discontinuità contenente Si considerano le relative equazioni di trasporto lungo e si determinano casi interessanti in cui queste equazioni sono mutuamente equivalenti in senso opportuno. Si stabiliscono teoremi di confronto per il trasporto delle discontinuità del tipo suddetto, relative ai sistemi d^hL[u]/dt^h=0 (h=0, , –m) e/o alla linearizazione del sistema L[u]=0 attorno a qualche sua soluzione regolare.

     

Rapid proton diffusion in microfluidic devices by means of micro-LIF technique

The micro-LIF (laser-induced fluorescence) technique was applied to the measurement of pH distributions in a chemically reacting flow in a microfluidic device. Two liquid streams were combined at the junction of a Y-shaped microchannel (100-m width and 33-m depth), and allowed to diffuse into each other and react. The results for non-reacting fluids (hydrochloric acid and water) show good agreement with theoretical values calculated using conventional diffusion. When a reaction occurred (hydrochloric acid and sodium hydroxide), a large difference between the measurement results and the theoretical values was observed, indicating rapid proton diffusion compared with the theoretically calculated values.     

A technique for measuring Lagrangian and Eulerian particle statistics in a turbulent flow

An experimental technique is described which has been developed to study particle dispersion in a round turbulent jet. Droplets are injected on the jet axis, and a laser sheet and position sensitive photomultiplier tube are used to track their radial displacement. Data processing is greatly simplified compared to video or photo imaging techniques which provide similar measurements. Statistically large samples are used to calculate dispersion and axial velocity as a function of axial downstream distance or particle time-of-flight. Dispersion and velocity statistics can be computed which are Lagrangian or Eulerian in nature. The technique has been demonstrated with 69 m droplets of hexadecane in a jet of air with a Reynolds number of 15,000; in principle it could be used to study the motion of very small, quasi-fluid particles.     

Microbubble lensing-induced photobleaching (μ-BLIP) with application to microflow visualization

The curvature of gas–liquid interfaces and the step change in properties across these interfaces in microchannels are shown here to create a powerful lens/mirror effect. In a hydrophilic system, light incident on the bubble is focused into the surrounding liquid, resulting in a locally increased total light exposure. The optical phenomena leading to this are discussed, and the effect is demonstrated experimentally by imaging the increased photobleaching rate of fluorophores in the near-bubble region. Numerical simulations of the system are performed to investigate the electrical potential and flow fields resulting from the application of an axial electric field. Microbubble lensing-induced photobleaching (-BLIP) is then applied as a method to inject a negative scalar flow marker for flow visualization in microchannels. Once formed, the electrokinetic transport of this marker is analyzed to determine the cross-channel velocity profile of the liquid phase and the liquid velocity in the film. Experimental data is verified by comparison with numerical predictions and previous experimental studies. This contribution represents both a new application of microscale gas–liquid interfacial phenomena, and a new technique for microfluidic flow visualization, particularly applicable (though not limited) to the study of multiphase microchannel flows.     

Shear Instability at the “Explosion Product–Metal” Interface for Sliding Detonation of an Explosive Charge

Periodic perturbations at the explosion product–metal interface were studied experimentally. Experiments were performed for both spherical and plane geometry. Critical conditions of wave formation (detonation velocity of an explosive charge D 6.9 mm/sec) are determined, and an explanation of this effect is given. It is found experimentally that a dynamic pulse causes intense plastic strains at the explosion products–metal interface, leading to thermal softening of the steel boundary layer. In this layer, Kelvin–Helmholtz instability occurs. Calculationanalytical estimates of the critical boundary unstable wavelength agree satisfactorily with experimental results.     

Submicron deformation field measurements: Part 2. Improved digital image correlation

This is the second paper in a series of three devoted to the application of scanning tunneling microscopy (STM) to mechanics problems. In this paper, improvements to the digital image correlation method are outlined, a technique that compares digital images of a specimen surface before and after deformation to deduce its two-dimensional surface displacement field and strains. The necessity of using the framework of large deformation theory for accurately addressing rigid body rotations to reduce associated errors in the strain components is pointed out. In addition, the algorithm is extended to compute the three-dimensional surface displacement field from STM data; also, significant improvements are achieved in the rate as well as the robustness of the convergence. For (STM) topographs, the resolution yields 4.8 nm for the in-plane and 1.5 nm for the out-of-plane displacement components spanning an area of 10 m×10m.