Characteristics of flow fields induced by interfacial waves in two-layer fluid

When designing deep ocean structures, it is necessary to estimate the effects of internal waves on the platform and auxiliary parts such as tension leg, riser and mooring lines. Up to now, only a few studies are concerned with the internal wave velocity fields. By using the most representative two-layer model, we have analyzed the behavior of velocity field induced by interfacial wave in the present paper. We find that there may exist velocity shear of fluid particles in the upper and lower layers so that any structures in the ocean are subjected to shear force nearby the interface. In the meantime, the magnitude of velocity for long internal wave appears spatially uniform in the respective layer although they still decay exponentially. Finally, the temporal variation for Stokes and solitary waves are shown to be of periodical and pulse type.

Acoustic Calculation for Supersonic Turbulent Boundary Layer Flow

An approach which combines direct numerical simulation (DNS) with the Lighthill acoustic analogy theory is used to study the potential noise sources during the transition process of a Mach 2.25 flat plate boundary layer. The quadrupole sound sources due to the flow fluctuations and the dipole sound sources due to the fluctuating surface stress are obtained. Numerical results suggest that formation of the high shear layers leads to a dramatic amplification of amplitude of the fluctuating quadrupole sound sources. Compared with the quadrupole sound source, the energy of dipole sound source is concentrated in the relatively low frequency range.

Theoretical investigation of turbulent boundary layer over a wavy surface

The important features of the two-dimensional incompressible turbulent flow over a wavy surface of wavelength comparable with the boundary layer thickness are analyzed.

A turbulent field method using model equation for turbulent shear stress similar to the scheme of Bradshaw, Ferriss and Atwell (1967) is employed with suitable modification to cover the viscous sublayer. The governing differential equations are linearized based on the small but finite amplitude to wavelength ratio. An orthogonal wavy coordinate system, accurate to the second order in the amplitude ratio, is adopted to avoid the severe restriction to the validity of linearization due to the large mean velocity gradient near the wall. Analytic solution up to the second order is obtained by using the method of matched-asymptotic-expansion based on the large Reynolds number and hence the small skin friction coefficient.

In the outer part of the layer, the perturbed flow is practically "inviscid." Solutions for the velocity, Reynolds stress and also the wall pressure distributions agree well with the experimental measurement. In the wall region where the perturbed Reynolds stress plays an important role in the process of momentum transport, only a qualitative agreement is obtained. The results also show that the nonlinear second-order effect is negligible for amplitude ratio of 0.03. The discrepancies in the detailed structure of the velocity, shear stress, and skin friction distributions near the wall suggest modifications to the model are required to describe the present problem.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Visualization of turbulent wedges under favourable pressure gradients using both shear-sensitive and temperature-sensitive liquid crystals

Turbulent wedges induced by a 3D surface roughness placed in a laminar boundary layer over a flat plate were visualised for the first time using both shear-sensitive and temperature-sensitive liquid crystals. The experiments were carried out at three different levels of favourable pressure gradients. The purpose of this investigation was to examine the spreading angles of the turbulent wedges indicated by their associated surface shear stresses and heat transfer characteristics and hence obtain further insight about the difference in the behaviour of transitional momentum and thermal boundary layers when a streamwise pressure gradient exists. It was shown that under a zero pressure gradient the spreading angles indicated by the two types of liquid crystals are the same, but the difference increases as the level of favourable pressure gradient increases. The result from the present study could have an important implication to the transition modelling of thermal boundary layers over gas turbine blades.

The shear mode of multilayer graphene

The quest for materials capable of realizing the next generation of electronic and photonic devices continues to fuel research on the electronic, optical and vibrational properties of graphene. Few-layer graphene (FLG) flakes with less than ten layers each show a distinctive band structure. Thus, there is an increasing interest in the physics and applications of FLGs. Raman spectroscopy is one of the most useful and versatile tools to probe graphene samples. Here, we uncover the interlayer shear mode of FLGs, ranging from bilayer graphene (BLG) to bulk graphite, and suggest that the corresponding Raman peak measures the interlayer coupling. This peak scales from ∼43cm -1 in bulk graphite to ∼31cm -1 in BLG. Its low energy makes it sensitive to near-Dirac point quasiparticles. Similar shear modes are expected in all layered materials, providing a direct probe of interlayer interactions. © 2012 Macmillan Publishers Limited. All rights reserved.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

Reliability study of the electroless Ni–P layer against solder alloy

Copper (Cu) has been widely used in the under bump metallurgy of chip and substrate metallization for chip packaging. However, due to the rapid formation of Cu–Sn intermetallic compound (IMC) at the tin-based solder/Cu interface during solder reaction, the reliability of this type of solder joint is a serious concern. In this work, electroless nickel–phosphorous (Ni–P) layer was deposited on the Cu pad of the flexible substrate as a diffusion barrier between Cu and the solder materials. The deposition was carried out in a commercial acidic sodium hypophosphite bath at 85 °C for different pH values. It was found that for the same deposition time period, higher pH bath composition (mild acidic) yields thicker Ni–P layer with lower phosphorous content. Solder balls having composition 62%Sn–36%Pb–2%Ag were reflowed at 240 °C for 1 to 180 min on three types of electroless Ni–P layers deposited at the pH value of 4, 4.8 and 6, respectively. Thermal stability of the electroless Ni–P barrier layer against the Sn–36%Pb–2%Ag solder reflowed for different time periods was examined by scanning electron microscopy equipped with energy dispersed X-ray. Solder ball shear test was performed in order to find out the relationship between the mechanical strength of solder joints and the characteristics of the electroless Ni–P layer deposited. The layer deposited in the pH 4 acidic bath showed the weak barrier against reflow soldering whereas layer deposited in pH 6 acidic bath showed better barrier against reflow soldering. Mechanical strength of the joints were deteriorated quickly in the layer deposited at pH 4 acidic bath, which was found to be thin and has a high phosphorous content. From the cross-sectional studies and fracture surface analyses, it was found that the appearance of the dark crystalline phosphorous-rich Ni layer weakened the interface and hence lower solder ball shear strength. Ni–Sn IMC formed at the interfaces was found to be more stable at the low phosphorous content (∼14 at.%) layer. Electroless Ni–P deposited at mild acidic bath resulting phosphorous content of around 14 at.% is suggested as the best barrier layer for Sn–36%Pb–2%Ag solder.

A multi-fibre multi-layer representative volume element (M2RVE) for prediction of matrix and interfacial damage in composite laminates

Multiscale micro-mechanics theory is extensively used for the prediction of the material response and damage analysis of unidirectional lamina using a representative volume element (RVE). Th is paper presents a RVE-based approach to characterize the materi al response of a multi-fibre cross-ply laminate considering the effect of matrix damage and fibre-matrix interfacial strength. The framework of the homogenization theory for periodic media has been used for the analysis of a 'multi-fibre multi-layer representative volume element' (M2 RVE) representing cross-ply laminate. The non-homogeneous stress-strain fields within the M2RVE are related to the average stresses and strains by using Gauss theorem and the Hill-Mandal strain energy equivalence principle. The interfacial bonding strength affects the in-plane shear stress-strain response significantl y. The material response predicted by M2 RVE is in good agreement with the experimental results available in the literature. The maximum difference between the shear stress predicted using M2 RVE and the experimental results is ~15% for the bonding strength of 30MPa at the strain value of 1.1%

Effect of the surface water layer on the optical signal in apertureless scanning near field optical microscopy

Optical signals measured in apertureless scanning near field optical microscopy (ASNOM) under ambient conditions are found to be affected significantly by the thin water layer absorbed on the surface under investigation, the presence of which is detected through measurements of the shear force experienced by the tip. This water layer also results in a large hysteresis between optical signals measured during approach and withdrawal of the tip to the sample surface. The role of this effect in ASNOM is anticipated to be significant, with the possibility of resultant topographically induced artefacts for ASNOM involving intermittent contact of tip and sample, but also providing a potential mechanism for nanoscale optical resolution.

High Strength Thermoplastic Bonding for Multi-channel, Multi-layer Lab-on-Chip Devices for Ocean and Environmental applications

A solvent-vapour thermoplastic bonding process is reported which provides high strength bonding of PMMA over a large area for multi-channel and multi-layer microfluidic devices with shallow high resolution channel features. The bond process utilises a low temperature vacuum thermal fusion step with prior exposure of the substrate to chloroform (CHCl3) vapour to reduce bond temperature to below the PMMA glass transition temperature. Peak tensile and shear bond strengths greater than 3 MPa were achieved for a typical channel depth reduction of 25 µm. The device-equivalent bond performance was evaluated for multiple layers and high resolution channel features using double-side and single-side exposure of the bonding pieces. A single-sided exposure process was achieved which is suited to multi-layer bonding with channel alignment at the expense of greater depth loss and a reduction in peak bond strength. However, leak and burst tests demonstrate bond integrity up to at least 10 bar channel pressure over the full substrate area of 100 mm x 100 mm. The inclusion of metal tracks within the bond resulted in no loss of performance. The vertical wall integrity between channels was found to be compromised by solvent permeation for wall thicknesses of 100 µm which has implications for high resolution serpentine structures. Bond strength is reduced considerably for multi-layer patterned substrates where features on each layer are not aligned, despite the presence of an intermediate blank substrate. Overall a high performance bond process has been developed that has the potential to meet the stringent specifications for lab-on-chip deployment in harsh environmental conditions for applications such as deep ocean profiling.

Ocean Boundary Layer Dynamics and Air-Sea Interaction

Thesis (Ph.D.)--University of Washington, 2005-12

Studies on the Vertical Structure of Horizontal Wind Variability in the Surface Boundary Layer Over Sriharikota

The present study on the vertical structure of horizontal wind variability in the surface boundary layer over Sriharikota. Based on clock wind speed and direction measuring meteorological tower facility from seven levels in the 100 m layer. The study on wind variability and elliptical approximation of wind hodographs investigated for this tropical coastal station established that Sriharikota is of meso-scale weather entity. Wind variability ratio increases from lower levels to upper levels. In South West monsoon months the station is of high ratio values and it gets affected with meso-scale weather features like thunderstorms. Average total shears are observed greater values than scalar shears. Scalar shears are high in the lowest shear levels compared to upper levels. Semi diurnal types of oscillation in average total shears are found in south west monsoon months. During cyclonic storm passage it is observed that there can be significant difference in mean wind speed from 10 m to 100 m level, but it is not so for peak wind speeds. The variations in wind variability ratio in different months is clearly depicted its strong link to define meso-scale or synoptic –scale forcing domination for this station. Meso-scale forcing is characterized by diurnal wind variability and synoptic- scale forcing by interdiurnal wind variability.

Diurnal Structure of Monsoon Boundary Layer

The characteristics of the monsoon boundary layer are imperative to understand in the perception of the tropical regions. The southwest monsoon is associated with a strong wind in the lower troposphere near 1.5 km and is referred to as Low Level Jet stream (LLJ). The boundary layer structure associated with the LLJ during monsoon can be studied using L-band Ultra High Frequency (UHF) radar. This L-band wind profiler-commonly referred as lower atmospheric wind profiler (LAWP), was installed at NARL, Gadanki. Zonal, meridional and vertical wind components are used to understand the diurnal variation of the wind in the Atmospheric Boundary Layer (ABL) and associated features. From the analysis during non rainy days of the southwest monsoon, it is found that the LLJ has maximum strength during the early morning hours at lower level and the height increases as day progresses. The vertical wind shows the transfer of momentum from the LLJ towards the surface, indicating the sinking motion during the daytime. Vertical gradient of the wind shear shows the intensity of clear air turbulence is moderate and no severe clear air turbulence is noticed during the monsoon period

Langmuir turbulence and deeply penetrating jets in an unstratified mixed layer

The influence of surface waves and an applied wind stress is studied in an ensemble of large eddy simulations to investigate the nature of deeply penetrating jets into an unstratified mixed layer. The influence of a steady monochromatic surface wave propagating parallel to the wind direction is parameterized using the wave-filtered Craik-Leibovich equations. Tracer trajectories and instantaneous downwelling velocities reveal classic counterrotating Langmuir rolls. The associated downwelling jets penetrate to depths in excess of the wave's Stokes depth scale, δs. Qualitative evidence suggests the depth of the jets is controlled by the Ekman depth scale. Analysis of turbulent kinetic energy (tke) budgets reveals a dynamical distinction between Langmuir turbulence and shear-driven turbulence. In the former, tke production is dominated by Stokes shear and a vertical flux term transports tke to a depth where it is dissipated. In the latter, tke production is from the mean shear and is locally balanced by dissipation. We define the turbulent Langmuir number Lat = (v*/Us)0.5 (v* is the ocean's friction velocity and Us is the surface Stokes drift velocity) and a turbulent anisotropy coefficient Rt = /( + ). The transition between shear-driven and Langmuir turbulence is investigated by varying external wave parameters δs and Lat and by diagnosing Rt and the Eulerian mean and Stokes shears. When either Lat or δs are sufficiently small the Stokes shear dominates the mean shear and the flow is preconditioned to Langmuir turbulence and the associated deeply penetrating jets.