A k-epsilon model for turbulent weld pool convection in gas metal arc welding process

In this paper, we present a modified k - epsilon model capable of addressing turbulent weld-pool convection in a GMAW process, taking into account the morphology of the phase change interface during a Gas Metal Arc Welding (GMAW) process. A three-dimensional turbulence mathematical model has been developed to study the heat transfer and fluid flow within the weld pool by considering the combined effect of three driving forces, viz., buoyancy, Lorentz force and surface tension (Marangoni convection). Mass and energy transports by the droplets are considered through the thermal analysis of the electrode. The falling droplet's heat addition to the molten pool is considered to be a volumetric heat source distributed in an imaginary cylindrical cavity ("cavity model") within the weld pool. This nature of heat source distribution takes into account the momentum and the thermal, energy of the falling droplets. The numerically predicted weld pool dimensions both from turbulence and laminar models are then compared with the experimental post-weld results sectioned across the weld axis. The above comparison enables us to analyze the overall effects of turbulent convection on the nature of heat and fluid flow and hence on the weld pool shape/size during the arc welding processes.

An EDC-based turbulent premixed combustion model

In the present study, a new turbulent premixed combustion model is proposed by integrating the Coherent Flame Model with the modified eddy dissipation concept, and relating the fine structure mass fraction to the flame surface density. First, experimental results of turbulent flame speed available from literature are compared with the predicted results at different turbulence intensities to validate the flame surface density model. It is observed that the model is able to predict the turbulent burning speeds accurately. Then, a comprehensive validation is carried out utilizing data on a turbulent lifted methane flame issuing into a vitiated co-flow. Detailed comparison of temperature and species concentrations between experiment and simulation is performed at different heights of the flame. Overall, the model is found to predict both the spatial variation and peak values of the scalars at various heights satisfactorily.

Generalization of the collision cone approach for motion safety in 3-D environments

Avoidance of collision between moving objects in a 3-D environment is fundamental to the problem of planning safe trajectories in dynamic environments. This problem appears in several diverse fields including robotics, air vehicles, underwater vehicles and computer animation. Most of the existing literature on collision prediction assumes objects to be modelled as spheres. While the conservative spherical bounding box is valid in many cases, in many other cases, where objects operate in close proximity, a less conservative approach, that allows objects to be modelled using analytic surfaces that closely mimic the shape of the object, is more desirable. In this paper, a collision cone approach (previously developed only for objects moving on a plane) is used to determine collision between objects, moving in 3-D space, whose shapes can be modelled by general quadric surfaces. Exact collision conditions for such quadric surfaces are obtained and used to derive dynamic inversion based avoidance strategies.

Flow visualization in turbulent free convection over horizontal smooth and grooved surfaces

The present paper discusses the flow visualization for turbulent free convection in a tank of water with the bottom surface being a smooth or a grooved surface and the top of the water surface exposed to ambient. The grooved surface is of parallel 90 degrees V-grooves with groove height of 10 mm and groove width of 20 mm. The experiment is carried out with aspect ratio (AR) of 2.9 and Rayleigh number (Ra) in the range, 1.3 x 10(7) - 4 x 10(7). Here AR is the aspect ratio (= width of fluid layer/height of fluid layer). Heat flux at the bottom surface is from electrical heating. From the pH-dye visualization, interesting flow structures are observed and these structures are analyzed with the help of plumes dynamics and temperature variations with time. (C) 2011 Elsevier Ltd. All rights reserved.

Cumulus clouds and convective boundary layers: a tropical perspective on two turbulent shear flows

This paper is a review prepared for the second Marseille Colloquium on the mechanics of turbulence, held in 2011, 50 years after the first. The review covers recent developments in our understanding of the large-scale dynamics of cumulus cloud flows and of the atmospheric boundary layer in the low-wind convective regime that is often encountered in the tropics. It has recently been shown that a variety of cumulus cloud forms and life cycles can be experimentally realized in the laboratory, with the transient diabatic plume taken as the flow model for a cumulus cloud. The plume is subjected to diabatic heating scaled to be dynamically similar to heat release from phase changes in clouds. The experiments are complemented by exact numerical solutions of the Navier-Stokes-Boussinesq equations for plumes with scaled off-source heating. The results show that the Taylor entrainment coefficient first increases with heating, reaches a positive maximum and then drops rapidly to zero or even negative values. This reduction in entrainment is a consequence of structural changes in the flow, smoothing out the convoluted boundaries in the non-diabatic plume, including the tongues engulfing the ambient flow. This is accompanied by a greater degree of mixedness in the core flow because of lower dilution by the ambient fluid. The cloud forms generated depend strongly on the history of the diabatic heating profile in the vertical direction. The striking effects of heating on the flow are attributable to the operation of the baroclinic torque due to the temperature field. The mean baroclinic torque is shown to peak around a quasi-cylindrical sheet situated midway between the axis of the flow and the edges. This torque is shear-enhancing and folds down the engulfment tongues. The increase in mixedness can be traced to an explosive growth in the enstrophy, triggered by a strong fluctuating baroclinic torque that acts as a source, especially at the higher wave numbers, thus enhancing the mixedness. In convective boundary layers field measurements show that, under conditions prevailing in the tropics, the eddy fluxes of momentum and energy do not follow the Monin-Obukhov similarity. Instead, the eddy momentum flux is found to be linear in the wind speed at low winds; and the eddy heat flux is, to a first approximation, governed by free convection laws, with wind acting as a small perturbation on a regime of free convection. A new boundary layer code, based on heat flux scaling rather than wall-stress scaling, shows promising improvements in predictive skills of a general circulation model.

TURBULENT PUMPING OF MAGNETIC FLUX REDUCES SOLAR CYCLE MEMORY AND THUS IMPACTS PREDICTABILITY OF THE SUN'S ACTIVITY

Prediction of the Sun's magnetic activity is important because of its effect on space environment and climate. However, recent efforts to predict the amplitude of the solar cycle have resulted in diverging forecasts with no consensus. Yeates et al. have shown that the dynamical memory of the solar dynamo mechanism governs predictability, and this memory is different for advection- and diffusion-dominated solar convection zones. By utilizing stochastically forced, kinematic dynamo simulations, we demonstrate that the inclusion of downward turbulent pumping of magnetic flux reduces the memory of both advection- and diffusion-dominated solar dynamos to only one cycle; stronger pumping degrades this memory further. Thus, our results reconcile the diverging dynamo-model-based forecasts for the amplitude of solar cycle 24. We conclude that reliable predictions for the maximum of solar activity can be made only at the preceding minimum-allowing about five years of advance planning for space weather. For more accurate predictions, sequential data assimilation would be necessary in forecasting models to account for the Sun's short memory.

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.

Visualization of coherent structures in turbulent subsonic jet using planar laser induced fluorescence of acetone

In this paper, we present the molecular density distribution measurement in turbulent nitrogen jet (Re approximate to 3 x 10(3)), using acetone as molecular tracer. The tracer was seeded in the nitrogen jet by purging through the liquid acetone at ambient temperature. Planar laser sheet of 266 nm wavelength from frequency quadrupled, Q-switched, Nd:YAG laser was used as an excitation source. Emitted fluorescence images of jet flow field were recorded on CMOS camera. The dependence of planar laser induced fluorescence (PLIF) intensity on acetone vapor density was used to convert PLIF image of nitrogen jet into the density image on pixel-by-pixel basis. Instantaneous quantitative density image of nitrogen jet, seeded with acetone, was obtained. The arrowhead-shaped coherent turbulent structures were observed in the present work. It was found that coherent structures were non-overlapping with separate boundaries. Breaking of coherent structures into turbulence was clearly observed above four times jet width.

Scaling of turbulent flame speed for expanding flames with Markstein diffusion considerations

In this paper we clarify the role of Markstein diffusivity, which is the product of the planar laminar flame speed and the Markstein length, on the turbulent flame speed and its scaling, based on experimental measurements on constant-pressure expanding turbulent flames. Turbulent flame propagation data are presented for premixed flames of mixtures of hydrogen, methane, ethylene, n-butane, and dimethyl ether with air, in near-isotropic turbulence in a dual-chamber, fan-stirred vessel. For each individual fuel-air mixture presented in this work and the recently published iso-octane data from Leeds, normalized turbulent flame speed data of individual fuel-air mixtures approximately follow a Re-T,f(0.5) scaling, for which the average radius is the length scale and thermal diffusivity is the transport property of the turbulence Reynolds number. At a given Re-T,Re-f, it is experimentally observed that the normalized turbulent flame speed decreases with increasing Markstein number, which could be explained by considering Markstein diffusivity as the leading dissipation mechanism for the large wave number flame surface fluctuations. Consequently, by replacing thermal diffusivity with the Markstein diffusivity in the turbulence Reynolds number definition above, it is found that normalized turbulent flame speeds could be scaled by Re-T,M(0.5) irrespective of the fuel, equivalence ratio, pressure, and turbulence intensity for positive Markstein number flames.

Environments of extended radio sources in the Australia Telescope Low-Brightness Survey

We present a study of the environments of extended radio sources in the Australia Telescope Low-Brightness Survey (ATLBS). The radio sources were selected from the ATLBS Extended Source Sample, which is a well defined sample containing the most extended of radio sources in the ATLBS sky survey regions. The environments were analysed using 4-m Cerro-Tololo Inter-American Observatory Blanco telescope observations carried out for ATLBS fields in the Sloan Digital Sky Survey r(') band. We have estimated the properties of the environments using smoothed density maps derived from galaxy catalogues constructed using these optical imaging data. The angular distribution of galaxy density relative to the axes of the radio sources has been quantified by defining anisotropy parameters that are estimated using a new method presented here. Examining the anisotropy parameters for a subsample of extended double radio sources that includes all sources with pronounced asymmetry in lobe extents, we find good evidence for environmental anisotropy being the dominant cause for lobe asymmetry in that higher galaxy density occurs almost always on the side of the shorter lobe, and this validates the usefulness of the method proposed and adopted here. The environmental anisotropy parameters have been used to examine and compare the environments of Fanaroff-Riley Class I (FRI) and Fanaroff-Riley Class II (FRII) radio sources in two redshift regimes (z < 0.5 and z > 0.5). Wide-angle tail sources and head-tail sources lie in the most overdense environments. The head-tail source environments (for the HT sources in our sample) display dipolar anisotropy in that higher galaxy density appears to lie in the direction of the tails. Excluding the head-tail and wide-angle tail sources, subsamples of FRI and FRII sources from the ATLBS appear to lie in similar moderately overdense environments, with no evidence for redshift evolution in the regimes studied herein.

On the k(1)(-1) scaling in sink-flow turbulent boundary layers

Scaling of the streamwise velocity spectrum phi(11)(k(1)) in the so-called sink-flow turbulent boundary layer is investigated in this work. The present experiments show strong evidence for the k(1)(-1) scaling i.e. phi(11)(k(1)) = Lambda(1)U(tau)(2)k(1)(-1), where k(1)(-1) is the streamwise wavenumber and U-tau is the friction velocity. Interestingly, this k(1)(-1) scaling is observed much farther from the wall and at much lower flow Reynolds number (both differing by almost an order of magnitude) than what the expectations from experiments on a zero-pressure-gradient turbulent boundary layer flow would suggest. Furthermore, the coefficient A(1) in the present sink-flow data is seen to be non-universal, i.e. A(1) varies with height from the wall; the scaling exponent -1 remains universal. Logarithmic variation of the so-called longitudinal structure function, which is the physical-space counterpart of spectral k(1)(-1) scaling, is also seen to be non-universal, consistent with the non-universality of A(1). These observations are to be contrasted with the universal value of A(1) (along with the universal scaling exponent of 1) reported in the literature on zero-pressure-gradient turbulent boundary layers. Theoretical arguments based on dimensional analysis indicate that the presence of a streamwise pressure gradient in sink-flow turbulent boundary layers makes the coefficient A(1) non-universal while leaving the scaling exponent -1 unaffected. This effect of the pressure gradient on the streamwise spectra, as discussed in the present study (experiments as well as theory), is consistent with other recent studies in the literature that are focused on the structural aspects of turbulent boundary layer flows in pressure gradients (Harun etal., J. Flui(d) Mech., vol. 715, 2013, pp. 477-498); the present paper establishes the link between these two. The variability of A(1) accommodated in the present framework serves to clarify the ideas of universality of the k(1)(-1) scaling.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Flame surface statistics of constant-pressure turbulent expanding premixed flames

In this paper we investigate the local flame surface statistics of constant-pressure turbulent expanding flames. First the statistics of local length ratio is experimentally determined from high-speed planar Mie scattering images of spherically expanding flames, with the length ratio on the measurement plane, at predefined equiangular sectors, defined as the ratio of the actual flame length to the length of a circular-arc of radius equal to the average radius of the flame. Assuming isotropic distribution of such flame segments we then convolute suitable forms of the length-ratio probability distribution functions (pdfs) to arrive at the corresponding area-ratio pdfs. It is found that both the length ratio and area ratio pdfs are near log-normally distributed and shows self-similar behavior with increasing radius. Near log-normality and rather intermittent behavior of the flame-length ratio suggests similarity with dissipation rate quantities which stimulates multifractal analysis. (C) 2014 AIP Publishing LLC.

Scaling of pressure spectrum in turbulent boundary layers

Scaling of pressure spectrum in zero-pressure-gradient turbulent boundary layers is discussed. Spatial DNS data of boundary layer at one time instant (Re-theta = 4500) are used for the analysis. It is observed that in the outer regions the pressure spectra tends towards the -7/3 law predicted by Kolmogorov's theory of small-scale turbulence. The slope in the pressure spectra varies from -1 close to the wall to a value close to -7/3 in the outer region. The streamwise velocity spectra also show a -5/3 trend in the outer region of the flow. The exercise carried out to study the amplitude modulation effect of the large scales on the smaller ones in the near-wall region reveals a strong modulation effect for the streamwise velocity, but not for the pressure fluctuations. The skewness of the pressure follows the same trend as the amplitude modulation coefficient, as is the case for the velocity. In the inner region, pressure spectra were seen to collapse better when normalized with the local Reynolds stress (-(u'v') over bar) than when scaled with the local turbulent kinetic energy (q(2) = (u'(2)) over bar + (v'(2)) over bar + (w'(2)) over bar)