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.     

Full-field representation of discretely sampled surface deformation for displacement and strain analysis

A detailed evaluation of the feasibility of determining displacements and displacement gradients from measured surface displacement fields is presented. An improved methodology for both the estimation and elimination of noise is proposed. The methodology is used to analyze the gradients for three tests: (1) uniform rotation, (2) uniform strain, and (3) crack-tip displacement fields. Results of the study indicate that the proposed methodology can be used to extract the underlying two-dimensional displacements and their corresponding gradients from the noisy data with reasonable accuracy. Specifically, it is shown that (a) the digital correlation method for acquiring displacement fields has an error in strain of approximately 150 strain at each point, (b) the average strain in a region of uniform strain has much less error, typically on the order of 20 strain, (c) the displacement nolse present in digital correlation is very small, approximately 0.01 pixels, (d) the proposed methodology for reducing noise in the data is essential to the accurate evaluation of displacement gradients and (e) the successful evaluation of displacement and displacement gradients for all three cases indicates that the proposed methodology can be used both to quantify the displacement fields and to reasonably estimate the overall gradient trends.     

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.     

Submicron deformation field measurements: Part 3. Demonstration of deformation determinations

This is the third and last paper is a sequence devoted to an experimental investigation of deformation mechanisms at the submicron scale through the use of a specially designed scanning tunneling microscope. Its application, when used jointly with digital image correlation, as a tool for strain and deformation determinations is explored by way of two demonstrations. First, deformations in a uniaxially stressed, unplasticized (poly)vinylchloride sample are analyzed to yield the three-dimensional surface displacement field over a 10 m×10m area. Homogeneous deformations occur at the micrometer and large size scales. However, at the 100-nm scale, inhomogeneous deformations embedded in a homogeneous deformation field appear. The second example addresses the deformation field in the vicinity of an interface between a carbon fiber and the surrounding matrix under shear stresses along the fiber. This loading leads to shearing a sheath from the carbon fiber that is about half a micron thick.     

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.     

COD measurements at various positions along a crack

Load versus crack-opening displacement (COD) was measured at various positions along the border of a fatigue crack as it grew from a small surface crack on the edge of an aluminum specimen into a through-the-thickness crack. Displacements were measured with a laser-based interferometric system with a gage length of 70 m and a resolution of 0.01 m. These load-COD curves can be used to determine opening loads and thereby investigate the effect of closure on the growth of cracks. In general, the opening loads decrease as the crack grows. The compliances after the crack is fully open can also be measured, and show generally good agreement with the predictions of linear-elastic fracture mechanics.Paper was presented at the 1988 SEM Spring Conference on Experimental Mechanics held in Portland, OR on June 5–10.     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

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.     

An extended numerical method of calculating twodimensional or axisymmetric heated flows allowing for dissipation

Summary The analysis of [1] is extended to cover viscous flow for both plane twodimensional and axisymmetric flow. The equations used are the Navier-Stokes equations in streamwise coordinates, except that for simplicity terms of order are neglected, where is a coefficient of viscosity and the flow direction. The supersonic combustion chamber of [1] is used as an example, but the technique of specifying the streamlines in advance has to be modified for a boundary layer with pressure gradient. Thermal conductivity is included but does not have a very important effect. Flow separation is indicated in some circumstances.

---

Übersicht Das Rechenverfahren von [1] wird erweitert auf Strömungen mit Reibung und zwar für zweidimensionale ebene und rotationssymmetrische Strömungen. Benutzt werden die Navier-Stokesschen Gleichungen in Stromlinienkoordinaten; zur Vereinfachung werden Glieder der Größenordnung vernachlässigt, wobei der Reibungskoeffizient und die Richtung der Strömung sind. Als Beispiel dient die Überschallbrennkammer von [1]. Das Verfahren, im voraus Stromlinien festzulegen, muß bei Grenzschichten mit Druckgradient abgeändert werden. Wärmeleitfähigkeit wird berücksichtigt, aber ihr Einfluß ist gering. In einigen Fällen wird Strömungsablösung vorausgesagt.

     

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.     

The moiré method and the evaluation of principal-moment and stress directions

This paper discusses, in relation to the moiré method as used for the solution of plate bending and two-dimensional stress problems, two graphical techniques for the determination of the directions of principal moments and stresses.The so-called isoclinic method and the point method are described.The application of these new techniques on three different models—a circular disk under diametrically opposite loads and two different circular plates subjected to a lateral load—are fully discussed.The graphically determined principal-stress and moment directions show excellent agreement with analytically determined comparable values.Paper was presented at 1963 SESA Annual Meeting held in Boston, Mass., on November 6–8.     

A Wollaston prism interferometer implemented with a digitizer

Objective of this work was to implement a Wollaston prism interferometer with a digitizer so as to perform the evaluation procedure of interferograms automatically by means of a microcomputer. The adopted opto-electronic sensing device is a linear array of 1,024 bar-shaped (15 m × 26 m) silicon photodiodes with center-to-center spacing 25 m. A stepping motor can move the array board in order to investigate the whole interferometric image. Scanning time is about 10 ms/line. Application of the system to measurements of concentration profiles in some mass diffusion processes of a binary gas mixture proves the ability of the digitizer to accurately measure gas densities. The performance and the potentiality of the sensor are analyzed and discussed.     

Measurement of liquid film in microchannels using a laser focus displacement meter

This paper presents a new method for measuring the interfacial displacement of a liquid film in microchannels using a laser focus displacement meter (LFD). The purpose of the study is to clarify the effectiveness of the new method for obtaining detailed information concerning interfacial displacement, especially in the case of a thin liquid film, in microchannels and minichannels. To prevent the tube wall signal from disturbing that of the gas–liquid interface, a fluorocarbon tube with a water box was used; the refraction index of this device is the same as that for water. With this method, accurate instantaneous measurements of the interfacial displacement of the liquid film were achieved. The error caused by refraction of the laser beam passing through the acrylic water box and fluorocarbon tube was estimated analytically and experimentally. The formulated analytical equation can estimate the real interface displacement by using the measured displacement in a fluorocarbon tube of 25 m to 2.0 mm I.D. A preliminary test using fluorocarbon tubes of 1 mm and 2 mm I.D. showed that the corrected interface displacement calculated by the equation agreed with the real displacement to within a 1% margin of error. It was also confirmed that the LFD in the system could measure a liquid film of 0.25 m at the thinnest. We made simultaneous measurements of the interface in fluorocarbon tubes of 0.5 mm and 1 mm I.D. using the LFD and a high-speed video camera with a microscope. These showed that the LFD could measure the interface of a liquid film with high spatial and temporal resolution during annular, slug, and piston flow regimes. The data also clarified the existence of a thin liquid film of less than 1 m in thickness in the slug and annular flow regimes.     

Errors in particle sizing by LDA due to turbidity in the incident laser beams

Measurements have been made of the effect of flow turbidity on the visibility and pedestal amplitude of an anemometer signal when incident laser beams are interrupted by particulate flow. The purpose is to assess the likely accuracy of particle sizing and the reliability of discrimination between continuous and particulate phase velocities. Optical depths of field were varied between 2.5 × 10^–2 and 14 × 10^–2 mm the diameter of the interrupting particles ranged between 14 and 800 m in six discrete ranges and the corresponding void fractions lay between 0.003% and 0.378%. The incident beam diameter was approximately 400 m.The measured size is subject to both systematic and random errors when inferred from measurements of pedestal amplitude: the random error increases as the ratio of the incident beam diameter to that of the particulate phase decreases. Systematic errors corresponding to a 10% underestimation of diameter occur at void fractions of 0.003%, 0.01% and 0.018% for particles below 40 m 75 m and 105 m respectively over a 5 cm depth of field. The r.m.s. error is smaller than 7% for particles below 40 m for all conditions studied but increases with increasing diameter and exceeds 10% at void fractions greater than 0.1% for particles above about 100 m. The random error in measured diameter derived from measurements of visibility is influenced mostly by the flow turbidity over the 5 cm of the incident beams closest to the measuring volume. For interrupting particles smaller than about 100 m the r.m.s. error is similar to that for measurements based on the pedestal amplitude.Discrimination of the velocity signal between the particulate and dispersed phase, based on the separation of pedestal amplitudes, is likely to be unreliable if the particle diameter is comparable to the diameter of the incident beams and if the probability of two particles simultaneously present in each beam is not negligible. A method for estimating the level of turbidity at which discrimination is no longer possible is described.List of symbols b beam diameter - d particle diameter - D amplitude of high passed Doppler signal (equation 7) - D ₁ depth of field of water tank - E particle extinction coefficient (= 2) - I _1,2 instantaneous light intensity at LDA measuring volume of each incident beam - n particle concentration (number per unit volume) - N number of particles in the incident beam - P maximum amplitude of pedestal component of Doppler signal - T Transmittance of beam {(attenuated/unattenuated) beam intensity} - V signal visibility (equation 8) - relative error in measured particle diameter {(truemeasured)/true)} - wavelength of laser beam (632.8 nm)     

Singular solutions of some quasilinear elliptic equations

We study isolated singularities of the quasilinear equation in an open set of ^N , where 1 < p N, p -1 q < N(p — 1)/ (N -p). We prove that, for any positive solution, if a singularity at the origin is not removable then either or u(x)/(x) any positive constant as x 0 where is the fundamental solution of the p-harmonic equation: . Global positive solutions are also classified.     

Spin-up in a saturated porous medium

It is shown that, on the Brinkman model, spin-up is confined to boundary layers whose thickness is of order k ^1/2, and the spin-up is established in a time of order k/, where k, , and denote permeability, density, porosity and dynamic viscosity, respectively.     

Solitary waves at the interface of a two-layer fluid

In this paper, we discuss the solitary waves at the interface of a two-layer incompressible inviscid fluid confined by two horizontal rigid walls, taking the effect of surface tension into account. First of all, we establish the basic equations suitable for the model considered, and hence derive the Korteweg-de Vries (KdV) equation satisfied by the first-order elevation of the interface with the aid of the reductive perturbation method under the approximation of weak dispersion. It is found that the KdV solitary waves may be convex upward or downward. It depends on whether the signs of the coefficients and of the KdV equation are the same or not. Then we examine in detail two critical cases, in which the nonlinear effect and the dispersion effect cannot balance under the original approximation. Applying other appropriate approximations, we obtain the modified KdV equation for the critical case of first kind (=0), and conclude that solitary waves cannot exist in the case considered as >0, but may still occur as <0, being in the form other than that of the KdV solitary wave.As for the critical case of second kind (=0), we deduce the generalized KdV equation, for which a kind of oscillatory solitary waves may occur. In addition, we discuss briefly the near-critical cases. The conclusions in this paper are in good agreement with some classical results which are extended considerably.     

Measurements in microscopic flow with a solid-state LDA

The application of a solid-state laser doppler anemometer (LDA) to laminar flow occurring in a 175 m channel is reported. Details of the optical system are described.     

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.

     

Codimension two bifurcations of symmetric cycles in Hamiltonian systems with an antisymplectic involution

Consider a Hamiltonian system with parameters, such that there exists an involution which reverses this Hamiltonian system. Let us assume the linear part L at =0 has only nonzero purely imaginary eigen-values ±ib₁,..., ±ib_n. In this paper, we classify the typical bifurcations of families of symmetric periodic solutions of this system at resonance if b_i/b_j=±1, ±2, or ±1/2 and the number of parameters needed is one or two. First, one puts the Hamiltonians into a convenient normal form. Next, applying a Lyapunov-Schmidt reduction and making further manipulations, one can geta reduced bifurcation equation which can possess certain symmetry. Finally, by using elementary methods from singularity theory or isotopy methods, one obtains the desired bifurcation diagrams.