Vortex sheet simulation of a plane �canonical� mixing layer

Discrete vortex simulations of the mixing layer carried out in the past have usually involved large induced velocity fluctuations, and thus demanded rather long time-averaging to obtain satisfactory values of Reynolds stresses and third-order moments. This difficulty has been traced here, in part, to the use of discrete vortices to model what in actuality are continuous vortex sheets. We propose here a novel two-dimensional vortex sheet technique for computing mixing layer flow in the limit of infinite Reynolds number. The method divides the vortex sheet into constant-strength linear elements, whose motions are computed using the Biot-Savart law. The downstream far-field is modelled by a steady vorticity distribution derived by application of conical similarity from the solution obtained in a finite computational domain. The boundary condition on the splitter plate is satisfied rigorously using a doublet sheet. The computed large-scale roll-up of the vortex sheet is qualitatively similar to experimentally obtained shadow-graphs of the plane turbulent mixing layer. The mean streamwise velocity profile and the growth rate agree well with experimental data. The presently computed Reynolds stresses and third-order moments are comparable with experimental and previous vortex-dynamical results, without using any external parameter (such as the vortex core-size) of the kind often used in the latter. The computed autocorrelations are qualitatively similar to experimental results along the top and bottom edges of the mixing layer, and show a well-defined periodicity along the centreline. The accuracy of the present computation is independently established by demonstrating negligibly small changes in the five invariants (including the Hamiltonian) in vortex dynamics.

Vaporization dynamics of functional droplets in a hot laminar air jet

The vaporization characteristics of pendant droplets of various chemical compositions (like conventional fuels, alternative fuels and nanosuspensions) subjected to convective heating in a laminar air jet have been analyzed. Different heating conditions were achieved by controlling the air temperature and velocity fields around the droplet. A hybrid timescale has been proposed which incorporates the effects of latent heat of vaporization, saturation vapor pressure and thermal diffusivity. This timescale in essence encapsulates the different parameters that influence the droplet vaporization rate. The analysis further permits the evaluation of the effect of various parameters such as surrounding temperature, Reynolds number, far-field vapor presence, impurity content and agglomeration dynamics (nanosuspensions) in the droplet. Flow visualization has been carried out to understand the role of internal recirculation on the vaporization rate. The visualization indicates the presence of a single vortex cell within the droplet on account of the rotation and oscillation of the droplet due to aerodynamic load. External heating induced agglomeration in nanofluids leads to morphological changes during the vaporization process. These morphological changes and alteration in vaporization behavior have been assessed using high speed imaging of the diameter regression and Scanning Electron Microscopy images of the resultant precipitate. (C) 2012 Elsevier Ltd. All rights reserved.

Parabolized stability equation models for predicting large-scale mixing noise of turbulent round jets

Parabolized stability equation (PSE) models are being deve loped to predict the evolu-tion of low-frequency, large-scale wavepacket structures and their radiated sound in high-speed turbulent round jets. Linear PSE wavepacket models were previously shown to be in reasonably good agreement with the amplitude envelope and phase measured using a microphone array placed just outside the jet shear layer. 1,2 Here we show they also in very good agreement with hot-wire measurements at the jet center line in the potential core,for a different set of experiments. 3 When used as a model source for acoustic analogy, the predicted far field noise radiation is in reasonably good agreement with microphone measurements for aft angles where contributions from large -scale structures dominate the acoustic field. Nonlinear PSE is then employed in order to determine the relative impor-tance of the mode interactions on the wavepackets. A series of nonlinear computations with randomized initial conditions are use in order to obtain bounds for the evolution of the modes in the natural turbulent jet flow. It was found that n onlinearity has a very limited impact on the evolution of the wavepackets for St≥0. 3. Finally, the nonlinear mechanism for the generation of a low-frequency mode as the difference-frequency mode 4,5 of two forced frequencies is investigated in the scope of the high Reynolds number jets considered in this paper.

An Outward-Wave-Favouring Finite Element-Based Strategy for Exterior Acoustical Problems

This work presents a finite element-based strategy for exterior acoustical problems based on an assumed pressure form that favours outgoing waves. The resulting governing equation, weak formulation, and finite element formulation are developed both for coupled and uncoupled problems. The developed elements are very similar to conventional elements in that they are based on the standard Galerkin variational formulation and use standard Lagrange interpolation functions and standard Gaussian quadrature. In addition and in contrast to wave envelope formulations and their extensions, the developed elements can be used in the immediate vicinity of the radiator/scatterer. The method is similar to the perfectly matched layer (PML) method in the sense that each layer of elements added around the radiator absorbs acoustical waves so that no boundary condition needs to be applied at the outermost boundary where the domain is truncated. By comparing against strategies such as the PML and wave-envelope methods, we show that the relative accuracy, both in the near and far-field results, is considerably higher.

Effect of contact stresses on shape recovery of NiTiCu thin films

Permanent plastic deformation induced by mechanical contacts affects the shape recovery of shape memory alloys. To understand the shape recovery of NiTiCu thin films subjected to local contact stresses, systematic investigations are carried out by inducing varying levels of contact stresses using nanoindentation. The resulting indents are located precisely for imaging using a predetermined array consisting of different sized indents. Morphology and topography of these indents before and after shape recovery are characterized using Scanning Electron Microscope and Atomic Force Microscope quantitatively. Shape recovery is found to be dependent on the contact stresses at the low loads while the recovery ratio remains constant at 0.13 for higher loads. Shape recovery is found to occur mainly in depth direction of the indent, while far field residual stresses play very little role in the recovery. (C) 2014 Elsevier B.V. All rights reserved.

A general methodology for calculating mixed mode stress intensity factors and fracture toughness of solder joints with interfacial cracks

Solder joints in electronic packages undergo thermo-mechanical cycling, resulting in nucleation of micro-cracks, especially at the solder/bond-pad interface, which may lead to fracture of the joints. The fracture toughness of a solder joint depends on material properties, process conditions and service history, as well as strain rate and mode-mixity. This paper reports on a methodology for determining the mixed-mode fracture toughness of solder joints with an interfacial starter-crack, using a modified compact mixed mode (CMM) specimen containing an adhesive joint. Expressions for stress intensity factor (K) and strain energy release rate (G) are developed, using a combination of experiments and finite element (FE) analysis. In this methodology, crack length dependent geometry factors to convert for the modified CMM sample are first obtained via the crack-tip opening displacement (CTOD)-based linear extrapolation method to calculate the under far-field mode I and II conditions (f(1a) and f(2a)), (ii) generation of a master-plot to determine a(c), and (iii) computation of K and G to analyze the fracture behavior of joints. The developed methodology was verified using J-integral calculations, and was also used to calculate experimental fracture toughness values of a few lead-free solder-Cu joints. (C) 2014 Elsevier Ltd. All rights reserved.

Unsteady-Flow Of A Viscous-Fluid Between 2 Parallel Disks With A Time-Varying Gap Width And A Magnetic-Field

The unsteady incompressible viscous fluid flow between two parallel infinite disks which are located at a distance h(t*) at time t* has been studied. The upper disk moves towards the lower disk with velocity h'(t*). The lower disk is porous and rotates with angular velocity Omega(t*). A magnetic field B(t*) is applied perpendicular to the two disks. It has been found that the governing Navier-Stokes equations reduce to a set of ordinary differential equations if h(t*), a(t*) and B(t*) vary with time t* in a particular manner, i.e. h(t*) = H(1 - alpha t*)(1/2), Omega(t*) = Omega(0)(1 - alpha t*)(-1), B(t*) = B-0(1 - alpha t*)(-1/2). These ordinary differential equations have been solved numerically using a shooting method. For small Reynolds numbers, analytical solutions have been obtained using a regular perturbation technique. The effects of squeeze Reynolds numbers, Hartmann number and rotation of the disk on the flow pattern, normal force or load and torque have been studied in detail

Time dependent rotational flow of a viscous fluid over an infinite porous disk with a magnetic field

Both the semi-similar and self-similar flows due to a viscous fluid rotating with time dependent angular velocity over a porous disk of large radius at rest with or without a magnetic field are investigated. For the self-similar case the resulting equations for the suction and no mass transfer cases are solved numerically by quasilinearization method whereas for the semi-similar case and injection in the self-similar case an implicit finite difference method with Newton's linearization is employed. For rapid deceleration of fluid and for moderate suction in the case of self-similar flow there exists a layer of fluid, close to the disk surface where the sense of rotation is opposite to that of the fluid rotating far away. The velocity profiles in the absence of magnetic field are found to be oscillatory except for suction. For the accelerating freestream, (semi-similar flow) the effect of time is to reduce the amplitude of the oscillations of the velocity components. On the other hand the effect of time for the oscillating case is just the opposite.

Rural energy consumption patterns - A field study

This paper is a condensed version of the final report of a detailed field study of rural energy consumption patterns in six villages located west of Bangalore in the dry belt of Karnataka State in India. The study was carried out in two phases; first, a pilot study of four villages and second, the detailed study of six villages, the populations of which varied from around 350 to about 950. The pilot survey ended in late 1976, and most of the data was collected for the main project in 1977. Processing of the collected data was completed in 1980. The aim was to carry out a census survey, rather than a sample study. Hence, considerable effort was expended in production of both a suitable questionnaire, ensuring that all respondents were contacted, and devising methods which would accurately reflect the actual energy use in various energy-utilising activities. In the end, 560 households out of 578 (97%) were surveyed. The following ranking was found for the various energy sources in order of average percentage contribution to the annual total energy requirement: firewood, 81·6%; human energy, 7·7%; animal energy, 2·7%; kerosene, 2·1%; electricity, 0·6% and all other sources (rice husks, agro-wastes, coal and diesel fuel), 5·3%. In other words commercial fuels made only a small contribution to the overall energy use. It should be noted that dung cakes are not burned in this region. The average energy use pattern, sector by sector, again on a percentage basis, was as follows: domestic, 88·3%; industry, 4·7%; agriculture, 4·3%; lighting, 2·2% and transport, 0·5%. The total annual per capita energy consumption was 12·6 ± 1·2 GJ, giving an average annual household consumption of around 78·6 GJ.

Rural energy consumption patterns - A field study

This paper is a condensed version of the final report of a detailed field study of rural energy consumption patterns in six villages located west of Bangalore in the dry belt of Karnataka State in India. The study was carried out in two phases; first, a pilot study of four villages and second, the detailed study of six villages, the populations of which varied from around 350 to about 950. The pilot survey ended in late 1976, and most of the data was collected for the main project in 1977. Processing of the collected data was completed in 1980. The aim was to carry out a census survey, rather than a sample study. Hence, considerable effort was expended in production of both a suitable questionnaire, ensuring that all respondents were contacted, and devising methods which would accurately reflect the actual energy use in various energy-utilising activities. In the end, 560 households out of 578 (97%) were surveyed. The following ranking was found for the various energy sources in order of average percentage contribution to the annual total energy requirement: firewood, 81A·6%; human energy, 7A·7%; animal energy, 2A·7%; kerosene, 2A·1%; electricity, 0A·6% and all other sources (rice husks, agro-wastes, coal and diesel fuel), 5A·3%. In other words commercial fuels made only a small contribution to the overall energy use. It should be noted that dung cakes are not burned in this region. The average energy use pattern, sector by sector, again on a percentage basis, was as follows: domestic, 88A·3%; industry, 4A·7%; agriculture, 4A·3%; lighting, 2A·2% and transport, 0A·5%. The total annual per capita energy consumption was 12A·6 A± 1A·2 GJ, giving an average annual household consumption of around 78A·6 GJ.

Phase-Field Simulation of Fusion Interface Events during Solidification of Dissimilar Welds: Effect of Composition Inhomogeneity

We investigate the events near the fusion interfaces of dissimilar welds using a phase-field model developed for single-phase solidification of binary alloys. The parameters used here correspond to the dissimilar welding of a Ni/Cu couple. The events at the Ni and the Cu interface are very different, which illustrate the importance of the phase diagram through the slope of the liquidus curves. In the Ni side, where the liquidus temperature decreases with increasing alloying, solutal melting of the base metal takes place; the resolidification, with continuously increasing solid composition, is very sluggish until the interface encounters a homogeneous melt composition. The growth difficulty of the base metal increases with increasing initial melt composition, which is equivalent to a steeper slope of the liquidus curve. In the Cu side, the initial conditions result in a deeply undercooled melt and contributions from both constrained and unconstrained modes of growth are observed. The simulations bring out the possibility of nucleation of a concentrated solid phase from the melt, and a secondary melting of the substrate due to the associated recalescence event. The results for the Ni and Cu interfaces can be used to understand more complex dissimilar weld interfaces involving multiphase solidification.

Tendency toward crossover of the effective susceptibility exponent from its doubled Ising value to its doubled mean-field value near a double critical point

The critical behavior of osmotic susceptibility in an aqueous electrolyte mixture 1-propanol (1P)+water (W)+potassium chloride is reported. This mixture exhibits re-entrant phase transitions and has a nearly parabolic critical line with its apex representing a double critical point (DCP). The behavior of the susceptibility exponent is deduced from static light-scattering measurements, on approaching the lower critical solution temperatures (TL’s) along different experimental paths (by varying t) in the one-phase region. The light-scattering data analysis substantiates the existence of a nonmonotonic crossover behavior of the susceptibility exponent in this mixture. For the TL far away from the DCP, the effective susceptibility exponent γeff as a function of t displays a nonmonotonic crossover from its single limit three-dimensional (3D)-Ising value ( ∼ 1.24) toward its mean-field value with increase in t. While for that closest to the DCP, γeff displays a sharp, nonmonotonic crossover from its nearly doubled 3D-Ising value toward its nearly doubled mean-field value with increase in t. The renormalized Ising regime extends over a relatively larger t range for the TL closest to the DCP, and a trend toward shrinkage in the renormalized Ising regime is observed as TL shifts away from the DCP. Nevertheless, the crossover to the mean-field limit extends well beyond t>10−2 for the TL’s studied. The observed crossover behavior is attributed to the presence of strong ion-induced clustering in this mixture, as revealed by various structure probing techniques. As far as the critical behavior in complex or associating mixtures with special critical points (like the DCP) is concerned, our results indicate that the influence of the DCP on the critical behavior must be taken into account not only on the renormalization of the critical exponent but also on the range of the Ising regime, which can shrink with decrease in the influence of the DCP and with the extent of structuring in the system. The utility of the field variable tUL in analyzing re-entrant phase transitions is demonstrated. The effective susceptibility exponent as a function of tUL displays a nonmonotonic crossover from its asymptotic 3D-Ising value toward a value slightly lower than its nonasymptotic mean-field value of 1. This behavior in the nonasymptotic, high tUL region is interpreted in terms of the possibility of a nonmonotonic crossover to the mean-field value from lower values, as foreseen earlier in micellar systems.

Capacitive-Resistive Transients in terms of field quantities

Capacitive-resistive transients in extended media are discussed in tenns of electric field quantities. Obviously, in rhese problems, the contribution of the magnetlc field to the electric field is deemed negligible. For a simple lllusfratlve example, the field solution is compared with the circuit-theoretical resuit for the voltage and current. An algorithm for solving such transients in space and time doman with the help of a Laplace solver is presented. Any other Laplace solver can also be used far this purpose. Its applicability is demonstrated with three examples, one of which is chosen to have a circuit-theoretical solution.

Crystal structure of quenched and in-field electroluminescent phosphors

Electroluminescent zinc sulfide doped with copper and chloride (ZnS:Cu, Cl) powder was heated to 400°C and rapidly quenched to room temperature. Comparison between the quenched and non-quenched phosphors using synchrotron radiation X-ray powder diffraction (XRPD) (λ = 0.828692 Å) and X-ray absorption spectroscopy (XAS) was made. XRPD shows that the expected highly faulted structure is observed with excellent resolution out to 150° 2θ (or to (12 2 2) of the sphalerite phase). The quenched sample compared to the unheated sample shows a large change in peak ratios between 46.7° and 46.9°, which is thought to correspond to the wurtzite (0 0 6), (0 3 2) and sphalerite (3 3 3)/(5 1 1) peaks. Hence, a large proportion of this sphalerite diffraction is lost from the material upon rapid quenching, but not when the material is allowed to cool slowly. The Zn K-edge XAS data indicate that the crystalline structures are indistinguishable using this technique, but do give an indication that the electronic structure has altered due to changing intensity of the white line. It is noted that the blue electroluminescence (EL) emission bands are lost upon quenching: however, a large amount of total EL emission intensity is also removed, which is consistent with our findings. We report the XRPD of a working alternating-current electroluminescence device in the synchrotron X-ray beam, which exhibits a new diffraction pattern when the device is powered in an AC field even though the phosphor is fixed in the binder. Significantly, only a few crystals are required to yield the diffraction data because of the high flux X-ray source. These in panel data show multiple sharp diffraction lines spread out under the region, where capillary data show broad diffraction intensity indicating that the phosphor powder is comprised of unique crystals, each having different structures.

Network Rhythms Influence the Relationship between Spike-Triggered Local Field Potential and Functional Connectivity

Characterizing the functional connectivity between neurons is key for understanding brain function. We recorded spikes and local field potentials (LFPs) from multielectrode arrays implanted in monkey visual cortex to test the hypotheses that spikes generated outward-traveling LFP waves and the strength of functional connectivity depended on stimulus contrast, as described recently. These hypotheses were proposed based on the observation that the latency of the peak negativity of the spike-triggered LFP average (STA) increased with distance between the spike and LFP electrodes, and the magnitude of the STA negativity and the distance over which it was observed decreased with increasing stimulus contrast. Detailed analysis of the shape of the STA, however, revealed contributions from two distinct sources-a transient negativity in the LFP locked to the spike (similar to 0 ms) that attenuated rapidly with distance, and a low-frequency rhythm with peak negativity similar to 25 ms after the spike that attenuated slowly with distance. The overall negative peak of the LFP, which combined both these components, shifted from similar to 0 to similar to 25 ms going from electrodes near the spike to electrodes far from the spike, giving an impression of a traveling wave, although the shift was fully explained by changing contributions from the two fixed components. The low-frequency rhythm was attenuated during stimulus presentations, decreasing the overall magnitude of the STA. These results highlight the importance of accounting for the network activity while using STAs to determine functional connectivity.