Stability analysis of core-annular flow and neutral stability wave number

A general transition criterion is proposed in order to locate the core-annular flow pattern in horizontal and vertical oil-water flows. It is based on a rigorous one-dimensional two-fluid model of liquid-liquid two-phase flow and considers the existence of critical interfacial wave numbers related to a non-negligible interfacial tension term to which the linear stability theory still applies. The viscous laminar-laminar flow problem is fully resolved and turbulence effects on the stability are analyzed through experimentally obtained shape factors. The proposed general transition criterion includes in its formulation the inviscid Kelvin-Helmholtz`s discriminator. If a theoretical maximum wavelength is considered as a necessary condition for stability, a stability criterion in terms of the Eotvos number is achieved. Effects of interfacial tension, viscosity ratio, density difference, and shape factors on the stability of core-annular flow are analyzed in detail. The more complete modeling allowed for the analysis of the neutral-stability wave number and the results strongly suggest that the interfacial tension term plays an indispensable role in the correct prediction of the stable region of core-annular flow pattern. The incorporation of a theoretical minimum wavelength into the transition model produced significantly better results. The criterion predictions were compared with recent data from the literature and the agreement is encouraging. (C) 2007 American Institute of Chemical Engineers.

Bidimensional codes recorded on an oxide glass surface using a continuous wave CO(2) laser

Surface heat treatment in glasses and ceramics, using CO(2) lasers, has attracted the attention of several researchers around the world due to its impact in technological applications, such as lab-on-a-chip devices, diffraction gratings and microlenses. Microlens fabrication on a glass surface has been studied mainly due to its importance in optical devices (fiber coupling, CCD signal enhancement, etc). The goal of this work is to present a systematic study of the conditions for microlens fabrications, along with the viability of using microlens arrays, recorded on the glass surface, as bidimensional codes for product identification. This would allow the production of codes without any residues (like the fine powder generated by laser ablation) and resistance to an aggressive environment, such as sterilization processes. The microlens arrays were fabricated using a continuous wave CO(2) laser, focused on the surface of flat commercial soda-lime silicate glass substrates. The fabrication conditions were studied based on laser power, heating time and microlens profiles. A He-Ne laser was used as a light source in a qualitative experiment to test the viability of using the microlenses as bidimensional codes.

CONCENTRATION FIELDS NEAR AIR-WATER INTERFACES DURING INTERFACIAL MASS TRANSPORT: OXYGEN TRANSPORT AND RANDOM SQUARE WAVE ANALYSIS

Mass transfer across a gas-liquid interface was studied theoretically and experimentally, using transfer of oxygen into water as the gas-liquid system. The experimental results support the conclusions of a theoretical description of the concentration field that uses random square waves approximations. The effect of diffusion over the concentration records was quantified. It is shown that the peak of the normalized rills concentration fluctuation profiles must be lower than 0.5, and that the position of the peak of the rms value is an adequate measure of the thickness of the diffusive layer. The position of the peak is the boundary between the regions more subject to molecular diffusion or to turbulent transport of dissolved mass.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Experimental study of a Tollmien-Schlichting wave interacting with a shallow 3D roughness element

The paper is devoted to an experimental study of the effect of a shallow 3D roughness element on the evolution of a 2D Tollmien-Schlichting wave in a Blasius boundary layer. The experiments were carried out under controlled disturbance conditions on an airfoil section which could provide a long run with zero pressure gradient flow. A pneumatically driven slit source was used to introduce the Tollmien-Schilichting wave upstream of the lower branch of the neutral stability curve. A few wavelengths downstream, the T-S wave interacts with a cylindrical roughness element. The height of the roughness was slowly oscillating in time, which allows a continuous measurement of the T-S wave response downstream the roughness. The oscillation frequency was approximately 1500 times lower than the frequency of the studied Tollmien-Schlichting wave and therefore, behaved as a steady roughness with respect to the T-S wave. Hot wire anemometry was used to measure wall normal profiles and spanwise scans close to the maximum of the eigenfunction of the T-S wave. The oscillation of the roughness and the synchronization of all-equipments permitted the use of ensemble average techniques. Two different amplitudes of T-S waves with a non-dimensional frequency of F120E-06 were studied. They show a strong amplification of the disturbances in a small spanwise wave number range. The analysis of the wall normal T-S profiles suggests the growth of oblique modes.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

Ultrasonic Viscosity Measurement Using the Shear-Wave Reflection Coefficient with a Novel Signal Processing Technique

Real-time viscosity measurement remains a necessity for highly automated industry. To resolve this problem, many studies have been carried out using an ultrasonic shear wave reflectance method. This method is based on the determination of the complex reflection coefficient`s magnitude and phase at the solid-liquid interface. Although magnitude is a stable quantity and its measurement is relatively simple and precise, phase measurement is a difficult task because of strong temperature dependence. A simplified method that uses only the magnitude of the reflection coefficient and that is valid under the Newtonian regimen has been proposed by some authors, but the obtained viscosity values do not match conventional viscometry measurements. In this work, a mode conversion measurement cell was used to measure glycerin viscosity as a function of temperature (15 to 25 degrees C) and corn syrup-water mixtures as a function of concentration (70 to 100 wt% of corn syrup). Tests were carried out at 1 MHz. A novel signal processing technique that calculates the reflection coefficient magnitude in a frequency band, instead of a single frequency, was studied. The effects of the bandwidth on magnitude and viscosity were analyzed and the results were compared with the values predicted by the Newtonian liquid model. The frequency band technique improved the magnitude results. The obtained viscosity values came close to those measured by the rotational viscometer with percentage errors up to 14%, whereas errors up to 96% were found for the single frequency method.

Comparison between DC(+) and Square Wave AC SAW Current Outputs to Weld AISI 304 for Low-Temperature Applications

Welded equipment for cryogenic applications is utilized in chemical, petrochemical, and metallurgical industries. One material suitable for cryogenic application is austenitic stainless steel, which usually doesn`t present ductile/brittle transition temperature, except in the weld metal, where the presence of ferrite and micro inclusions can promote a brittle failure, either by ferrite cleavage or dimple nucleation and growth, respectively. A 25-mm- (1-in.-) thick AISI 304 stainless steel base metal was welded with the SAW process using a 308L solid wire and two kinds of fluxes and constant voltage power sources with two types of electrical outputs: direct current electrode positive and balanced square wave alternating current. The welded joints were analyzed by chemical composition, microstructure characterization, room temperature mechanical properties, and CVN impact test at -100 degrees C (-73 degrees F). Results showed that an increase of chromium and nickel content was observed in all weld beads compared to base metal. The chromium and nickel equivalents ratio for the weld beads were always higher for welding with square wave AC for the two types of fluxes than for direct current. The modification in the Cr(eq)/Ni(eq) ratio changes the delta ferrite morphology and, consequently, modifies the weld bead toughness at lower temperatures. The oxygen content can also affect the toughness in the weld bead. The highest absorbed energy in a CVN impact test was obtained for the welding condition with square wave AC electrical output and neutral flux, followed by DC(+) electrical output and neutral flux, and square wave AC electrical output and alloyed flux.

Estimating wave spectra from the motions of moored vessels: Experimental validation

The practicability of estimating directional wave spectra based on a vessel`s 1st order response has been recently addressed by several researchers. Different alternatives regarding statistical inference methods and possible drawbacks that could arise from their application have been extensively discussed, with an apparent preference for estimations based on Bayesian inference algorithms. Most of the results on this matter, however, rely exclusively on numerical simulations or at best on few and sparse full-scale measurements, comprising a questionable basis for validation purposes. This paper discusses several issues that have recently been debated regarding the advantages of Bayesian inference and different alternatives for its implementation. Among those are the definition of the best set of input motions, the number of parameters required for guaranteeing smoothness of the spectrum in frequency and direction and how to determine their optimum values. These subjects are addressed in the light of an extensive experimental campaign performed with a small-scale model of an FPSO platform (VLCC hull), which was conducted in an ocean basin in Brazil. Tests involved long and short crested seas with variable levels of directional spreading and also bimodal conditions. The calibration spectra measured in the tank by means of an array of wave probes configured the paradigm for estimations. Results showed that a wide range of sea conditions could be estimated with good precision, even those with somewhat low peak periods. Some possible drawbacks that have been pointed out in previous works concerning the viability of employing large vessels for such a task are then refuted. Also, it is shown that a second parameter for smoothing the spectrum in frequency may indeed increase the accuracy in some situations, although the criterion usually proposed for estimating the optimum values (ABIC) demands large computational effort and does not seem adequate for practical on-board systems, which require expeditious estimations. (C) 2009 Elsevier Ltd. All rights reserved.

Directional Wave Spectrum Estimation Based on Vessel`s 1st-order Motions: Field Results

In this paper, 2 different approaches for estimating the directional wave spectrum based on a vessel`s 1st-order motions are discussed, and their predictions are compared to those provided by a wave buoy. The real-scale data were obtained in an extensive monitoring campaign based on an FPSO unit operating at Campos Basin, Brazil. Data included vessel motions, heading and tank loadings. Wave field information was obtained by means of a heave-pitch-roll buoy installed in the vicinity of the unit. `two of the methods most widely used for this kind of analysis are considered, one based on Bayesian statistical inference, the other consisting of a parametrical representation of the wave spectrum. The performance of both methods is compared, and their sensitivity to input parameters is discussed. This analysis complements a set of previous validations based on numerical and towing-tank results and allows for a preliminary evaluation of reliability when applying the methodology at full scale.

Generalized Coupled Algebraic Riccati Equations for Discrete-time Markov Jump with Multiplicative Noise Systems

In this paper we consider the existence of the maximal and mean square stabilizing solutions for a set of generalized coupled algebraic Riccati equations (GCARE for short) associated to the infinite-horizon stochastic optimal control problem of discrete-time Markov jump with multiplicative noise linear systems. The weighting matrices of the state and control for the quadratic part are allowed to be indefinite. We present a sufficient condition, based only on some positive semi-definite and kernel restrictions on some matrices, under which there exists the maximal solution and a necessary and sufficient condition under which there exists the mean square stabilizing solution fir the GCARE. We also present a solution for the discounted and long run average cost problems when the performance criterion is assumed be composed by a linear combination of an indefinite quadratic part and a linear part in the state and control variables. The paper is concluded with a numerical example for pension fund with regime switching.

New Analytic Solution Related to the Richards, Philip, and Green-Ampt Equations for Infiltration

Based on physical laws of similarity, an analytic solution of the soil water potential form of the Richards equation was derived for water infiltration into a homogeneous sand. The derivation assumes a similarity between the soil water retention function and that of the soil water content profiles taken at fixed times. The new solution successfully described soil water content profiles experimentally measured for water infiltrating downward, upward, and horizontally into a homogeneous sand and agrees with that presented by Philip in 1957. The utility of this analysis is still to be verified, but it is expected to hold for soils that have a narrow pore-size distribution before wetting and that manifest a sharp increase of water content at the wetting front during infiltration. The effect of van Genuchten`s parameters alpha and n on the application of the solution to other porous media was investigated. The solution also improves and provides a more realistic description of the infiltration process than that pioneered by Green and Ampt in 1911.

Effect of follicular wave synchronization on in vitro embryo production in heifers

Aiming to achieve the ideal time of ovum pick-up (OPU) for in vitro embryo production (IVP) in crossbred heifers, two Latin square design studies investigated the effect of ovarian follicular wave synchronization with estradiol benzoate (EB) and progestins. For each experiment, crossbred heifers stage of estrous cycle was synchronized either with a norgestomet ear implant (Experiment 1) or a progesterone intravaginal device (Experiment 2) for 7d, followed by the administration of 150 mu g D-cloprostenol. On Day 7, all follicles >3 mm in diameter were aspirated and implants/devices were replaced by new ones. Afterwards, implant/device replacement was conducted every 14 d. Each experiment had three treatment groups. In Experiment I (n = 12), heifers in Group 2X had their follicles aspirated twice a week and those in Groups 1X and 1X-EB were submitted to OPU once a week for a period of 28 d. Heifers from Group 1X-EB also received 2 mg EB i.m. immediately after each OPU session. In Experiment 2 (n = 11), animals from Group 0EB did not receive EB while heifers in Groups 2EB and 5EB received 2 and 5 mg of EB respectively, immediately after OPU. The OPU sessions were performed once weekly for 28 d. Therefore, in both experiments, four OPU sessions were performed in heifers aspirated once a week and in Experiment 1, eight OPU sessions were done in heifers aspirated twice a week. Additionally, during the 7-d period following follicular aspiration, ovarian ultrasonography examinations were conducted to measure diameter of the largest follicle and blood samples were collected for FSH quantification by RIA. In Experiment 1, all viable oocytes recovered were in vitro matured and fertilized. Results indicated that while progestin and EB altered follicular wave patterns, this treatment did not prevent establishment of follicular dominance on the ovaries of heifers during OPU at 7-d intervals. Furthermore, the proposed stage of follicular wave synchronization strategies did not improve the number and quality of the recovered oocytes, or the number of in vitro produced embryos. (C) 2009 Elsevier B.V. All rights reserved.

Standing wave bed forms at Plage de la Grande Roche, le Verdelet, Pléneuf-Val-André, Cote d’Armor (France)

Coastal Photograph by Hubert Chanson This photograph of standing wave bed forms was taken at very low tide. The tidal range was 10 m. The bed forms were located on the island of Le Verdelet, in a channel between Le Grande Jaune and Le Verdelet. It is likely that these standing wave bed forms were formed during transcritical shallow water flows at the end of ebb tide. The author’s watch is in the foreground for scale. (Coastal Photograph by Hubert Chanson, Division of Civil Engineering, the University of Queensland, Brisbane, Queensland 4072, Australia.)

Quantum noise in optical fibers. I. Stochastic equations

We analyze the quantum dynamics of radiation propagating in a single-mode optical fiber with dispersion, nonlinearity, and Raman coupling to thermal phonons. We start from a fundamental Hamiltonian that includes the principal known nonlinear effects and quantum-noise sources, including linear gain and loss. Both Markovian and frequency-dependent, non-Markovian reservoirs are treated. This treatment allows quantum Langevin equations, which have a classical form except for additional quantum-noise terms, to be calculated. In practical calculations, it is more useful to transform to Wigner or 1P quasi-probability operator representations. These transformations result in stochastic equations that can be analyzed by use of perturbation theory or exact numerical techniques. The results have applications to fiber-optics communications, networking, and sensor technology.