Location and scale mixtures of Gaussians with flexible tail behaviour: Properties, inference and application to multivariate clustering

The family of location and scale mixtures of Gaussians has the ability to generate a number of flexible distributional forms. The family nests as particular cases several important asymmetric distributions like the Generalized Hyperbolic distribution. The Generalized Hyperbolic distribution in turn nests many other well known distributions such as the Normal Inverse Gaussian. In a multivariate setting, an extension of the standard location and scale mixture concept is proposed into a so called multiple scaled framework which has the advantage of allowing different tail and skewness behaviours in each dimension with arbitrary correlation between dimensions. Estimation of the parameters is provided via an EM algorithm and extended to cover the case of mixtures of such multiple scaled distributions for application to clustering. Assessments on simulated and real data confirm the gain in degrees of freedom and flexibility in modelling data of varying tail behaviour and directional shape.

The sensitivity of the northern hemisphere summer mean meridional circulation to changes in the large scale eddy forcing

A zonally averaged version of the Goddard Laboratory for Atmospheric Sciences (GLAS) climate model is used to study the sensitivity of the northern hemisphere (NH) summer mean meridional circulation to changes in the large scale eddy forcing. A standard solution is obtained by prescribing the latent heating field and climatological horizontal transports of heat and momentum by the eddies. The radiative heating and surface fluxes are calculated by model parameterizations. This standard solution is compared with the results of several sensitivity studies. When the eddy forcing is reduced to 0.5 times or increased to 1.5 times the climatological values, the strength of the Ferrel cells decrease or increase proportionally. It is also seen that such changes in the eddy forcing can influence the strength of theNH Hadley cell significantly. Possible impact of such changes in the large scale eddy forcing on the monsoon circulation via changes in the Hadley circulation is discussed. Sensitivity experiments including only one component of eddy forcing at a time show that the eddy momentum fluxes seem to be more important in maintaining the Ferrel cells than the eddy heat fluxes. In the absence of the eddy heat fluxes, the observed eddy momentum fluxes alone produce subtropical westerly jets which are weaker than those in the standard solution. On the other hand, the observed eddy heat fluxes alone produce subtropical westerly jets which are stronger than those in the standard solution.

MACHAL [Mitnadvei Hutz LaAretz] and Aliyah Bet Records, undated, 1930-2011

The MACHAL, the acronym for “Mitnadvei Hutz LaAretz” ("Volunteers from Abroad"), consisted of about 3500 men and women from over 40 countries from a variety of social and religious backgrounds who volunteered to fight for the establishment of Israel. This collection is unique in that it deals specifically with the experience of MACHAL and Aliyah Bet volunteers from Canada and the United States and others living in the United States. The collections consists of files on 500 volunteers, over 2000 original and reproduction photographs, numerous audio-visual material, books, manuscripts, and memoirs.

Laminar and turbulent natural convection in rectangular cavities

A computer program has been developed for the prediction of buoyancy-driven laminar and turbulent flow in rectangular air-filled two-dimensional cavities with differentially heated side walls. Laminar flow predictions for a square cavity and Rayleigh numbers from Ra = 10^3 up to the onset of unsteady flow have been obtained. Accurate solutions for Ra = 5 x 10^6, 10^7, 5 x 10^7 and 10^8 are presented and an estimate for the critical Rayleigh number at which the steady laminar flow becomes unsteady is given for this geometry. Numerical predictions of turbulent flow have been obtained for RaH~0(10^9 -10^11 ) and compared with existing experimental data. A previously developed second moment closure model (Behnia et al. 1987) has been used to model the turbulence. Results indicate that a second moment closure model is capable of predicting the observed flow features.

Laminar and turbulent natural convection in rectangular cavities

"This chapter discusses laminar and turbulent natural convection in rectangular cavities. Natural convection in rectangular two-dimensional cavities has become a standard problem in numerical heat transfer because of its relevance in understanding a number of problems in engineering. Current research identified a number of difficulties with regard to the numerical methods and the turbulence modeling for this class of flows. Obtaining numerical predictions at high Rayleigh numbers proved computationally expensive such that results beyond Ra ∼ 1014 are rarely reported. The chapter discusses a study in which it was found that turbulent computations in square cavities can't be extended beyond Ra ∼ O (1012) despite having developed a code that proved very efficient for the high Ra laminar regime. As the Rayleigh number increased, thin boundary layers began to form next to the vertical walls, and the central region became progressively more stagnant and highly stratified. Results obtained for the high Ra laminar regime were in good agreement with existing studies. Turbulence computations, although of a preliminary nature, indicated that a second moment closure model was capable of predicting the experimentally observed flow features."--Publisher Summary

Computation of confined swirling flow using RNG and k-epsilon turbulence models

In this work we numerically model isothermal turbulent swirling flow in a cylindrical burner. Three versions of the RNG k-epsilon model are assessed against performance of the standard k-epsilon model. Sensitivity of numerical predictions to grid refinement, differing convective differencing schemes and choice of (unknown) inlet dissipation rate, were closely scrutinised to ensure accuracy. Particular attention is paid to modelling the inlet conditions to within the range of uncertainty of the experimental data, as model predictions proved to be significantly sensitive to relatively small changes in upstream flow conditions. We also examine the characteristics of the swirl--induced recirculation zone predicted by the models over an extended range of inlet conditions. Our main findings are: - (i) the standard k-epsilon model performed best compared with experiment; - (ii) no one inlet specification can simultaneously optimize the performance of the models considered; - (iii) the RNG models predict both single-cell and double-cell IRZ characteristics, the latter both with and without additional internal stagnation points. The first finding indicates that the examined RNG modifications to the standard k-e model do not result in an improved eddy viscosity based model for the prediction of swirl flows. The second finding suggests that tuning established models for optimal performance in swirl flows a priori is not straightforward. The third finding indicates that the RNG based models exhibit a greater variety of structural behaviour, despite being of the same level of complexity as the standard k-e model. The plausibility of the predicted IRZ features are discussed in terms of known vortex breakdown phenomena.

A Kalman filter based strategy for linear structural system identification based on multiple static and dynamic test data

The problem of identification of stiffness, mass and damping properties of linear structural systems, based on multiple sets of measurement data originating from static and dynamic tests is considered. A strategy, within the framework of Kalman filter based dynamic state estimation, is proposed to tackle this problem. The static tests consists of measurement of response of the structure to slowly moving loads, and to static loads whose magnitude are varied incrementally; the dynamic tests involve measurement of a few elements of the frequency response function (FRF) matrix. These measurements are taken to be contaminated by additive Gaussian noise. An artificial independent variable τ, that simultaneously parameterizes the point of application of the moving load, the magnitude of the incrementally varied static load and the driving frequency in the FRFs, is introduced. The state vector is taken to consist of system parameters to be identified. The fact that these parameters are independent of the variable τ is taken to constitute the set of ‘process’ equations. The measurement equations are derived based on the mechanics of the problem and, quantities, such as displacements and/or strains, are taken to be measured. A recursive algorithm that employs a linearization strategy based on Neumann’s expansion of structural static and dynamic stiffness matrices, and, which provides posterior estimates of the mean and covariance of the unknown system parameters, is developed. The satisfactory performance of the proposed approach is illustrated by considering the problem of the identification of the dynamic properties of an inhomogeneous beam and the axial rigidities of members of a truss structure.

Stationary processes for object detection and non-rigid structure-from-motion

Stationary processes are random variables whose value is a signal and whose distribution is invariant to translation in the domain of the signal. They are intimately connected to convolution, and therefore to the Fourier transform, since the covariance matrix of a stationary process is a Toeplitz matrix, and Toeplitz matrices are the expression of convolution as a linear operator. This thesis utilises this connection in the study of i) efficient training algorithms for object detection and ii) trajectory-based non-rigid structure-from-motion.

Investigations on voltages and currents in the lightning protection system of the Indian satellite launch pad-I during a stroke interception

A reliable protection against direct lightning hit is very essential for satellite launch pads. In view of this, suitable protection systems are generally employed. The evaluation of efficacy of the lightning protection schemes among others requires an accurate knowledge of the consequential potential rise at the struck point and the current injected into soil at the earth termination. The present work has made a detailed effort to deduce these quantities for the lightning protection scheme of the Indian satellite launch pad-I. A reduced scale model of the system with a frequency domain approach is employed for the experimental study. For further validation of the experimental approach, numerical simulations using numerical electromagnetic code-2 are also carried out on schemes involving single tower. The study results on the protection system show that the present design is quite safe with regard to top potential rise. It is shown that by connecting ground wires to the tower, its base current and, hence, the soil potential rise can be reduced. An evaluation of an alternate design philosophy involving insulated mast scheme is also made. The potential rise in that design is quantified and the possibility of a flashover to supporting tower is briefly looked into. The supporting tower is shown to have significant induced currents.

Eddy diffusivity: A single dispersion analysis of high resolution drifters in a tidal shallow estuary

In an estuary, mixing and dispersion resulting from turbulence and small scale fluctuation has strong spatio-temporal variability which cannot be resolved in conventional hydrodynamic models while some models employs parameterizations large water bodies. This paper presents small scale diffusivity estimates from high resolution drifters sampled at 10 Hz for periods of about 4 hours to resolve turbulence and shear diffusivity within a tidal shallow estuary (depth < 3 m). Taylor's diffusion theorem forms the basis of a first order estimate for the diffusivity scale. Diffusivity varied between 0.001 – 0.02 m2/s during the flood tide experiment. The diffusivity showed strong dependence (R2 > 0.9) on the horizontal mean velocity within the channel. Enhanced diffusivity caused by shear dispersion resulting from the interaction of large scale flow with the boundary geometries was observed. Turbulence within the shallow channel showed some similarities with the boundary layer flow which include consistency with slope of 5/3 predicted by Kolmogorov's similarity hypothesis within the inertial subrange. The diffusivities scale locally by 4/3 power law following Okubo's scaling and the length scale scales as 3/2 power law of the time scale. The diffusivity scaling herein suggests that the modelling of small scale mixing within tidal shallow estuaries can be approached from classical turbulence scaling upon identifying pertinent parameters.

Dynamical processes in the stratosphere and upper troposphere and their influence on the distribution of trace gases in the polar atmosphere

Transport plays an important role in the distribution of long-lived gases such as ozone and water vapour in the atmosphere. Understanding of observed variability in these gases as well as prediction of the future changes depends therefore on our knowledge of the relevant atmospheric dynamics. This dissertation studies certain dynamical processes in the stratosphere and upper troposphere which influence the distribution of ozone and water vapour in the atmosphere. The planetary waves that originate in the troposphere drive the stratospheric circulation. They influence both the meridional transport of substances as well as parameters of the polar vortices. In turn, temperatures inside the polar vortices influence abundance of the Polar Stratospheric Clouds (PSC) and therefore the chemical ozone destruction. Wave forcing of the stratospheric circulation is not uniform during winter. The November-December averaged stratospheric eddy heat flux shows a significant anticorrelation with the January-February averaged eddy heat flux in the midlatitude stratosphere and troposphere. These intraseasonal variations are attributable to the internal stratospheric vacillations. In the period 1979-2002, the wave forcing exhibited a negative trend which was confined to the second half of winter only. In the period 1958-2002, area, strength and longevity of the Arctic polar vortices do not exhibit significant long-term changes while the area with temperatures lower than the threshold temperature for PSC formation shows statistically significant increase. However, the Arctic vortex parameters show significant decadal changes which are mirrored in the ozone variability. Monthly ozone tendencies in the Northern Hemisphere show significant correlations (|r|=0.7) with proxies of the stratospheric circulation. In the Antarctic, the springtime vortex in the lower stratosphere shows statistically significant trends in temperature, longevity and strength (but not in area) in the period 1979-2001. Analysis of the ozone and water vapour vertical distributions in the Arctic UTLS shows that layering below and above the tropopause is often associated with poleward Rossby wave-breaking. These observations together with calculations of cross-tropopause fluxes emphasize the importance of poleward Rossby wave breaking for the stratosphere-troposphere exchange in the Arctic.

Aerosol characteristics and radiative impacts over the Arabian Sea during the intermonsoon season: Results from ARMEX field campaign

During the second phase of the Arabian Sea Monsoon Experiment (ARMEX-II), extensive measurements of spectral aerosol optical depth, mass concentration, and mass size distribution of ambient aerosols as well as mass concentration of aerosol black carbon (BC) were made onboard a research vessel during the intermonsoon period (i.e., when the monsoon winds are in transition from northeasterlies to westerlies/ southwesterlies) over the Arabian Sea (AS) adjoining the Indian Peninsula. Simultaneous measurements of spectral aerosol optical depths (AODs) were made at different regions over the adjoining Indian landmass. Mean AODs (at 500-nm wavelength) over the ocean (similar to0.44) were comparable to those over the coastal land (similar to0.47), but were lower than the values observed over the plateau regions of central Indian Peninsula (similar to0.61). The aerosol properties were found to respond distinctly with respect to change in the trajectories, with higher optical depths and flatter AOD spectra associated with trajectories indicating advection from west Asia, and northwest and west-coastal India. On average, BC constituted only similar to2.2% to total aerosol mass compared to the climatological values of similar to6% over the coastal land during the same season. These data are used to characterize the physical properties of aerosols and to assess the resulting short-wave direct aerosol forcing. The mean values were similar to27 W m(-2) at the surface and -12 W m(-2) at the top of the atmosphere (TOA), resulting in a net atmospheric forcing of +15 W m(-2). The forcing also depended on the region from where the advection predominates. The surface and atmospheric forcing were in the range -40 to -57 W m(-2) and +27 to +39 W m(-2), respectively, corresponding to advection from the west Asian and western coastal India where they were as low as -19 and +10 W m(-2), respectively, when the advection was mainly from the Bay of Bengal and from central/peninsular India. In all these cases, the net atmospheric forcing (heating) efficiency was lower than the values reported for northern Indian Ocean during northern winter, which is attributed to the reduced BC mass fraction.

Residual noise covariance for Planck low-resolution data analysis

Aims: Develop and validate tools to estimate residual noise covariance in Planck frequency maps. Quantify signal error effects and compare different techniques to produce low-resolution maps. Methods: We derive analytical estimates of covariance of the residual noise contained in low-resolution maps produced using a number of map-making approaches. We test these analytical predictions using Monte Carlo simulations and their impact on angular power spectrum estimation. We use simulations to quantify the level of signal errors incurred in different resolution downgrading schemes considered in this work. Results: We find an excellent agreement between the optimal residual noise covariance matrices and Monte Carlo noise maps. For destriping map-makers, the extent of agreement is dictated by the knee frequency of the correlated noise component and the chosen baseline offset length. The significance of signal striping is shown to be insignificant when properly dealt with. In map resolution downgrading, we find that a carefully selected window function is required to reduce aliasing to the sub-percent level at multipoles, ell > 2Nside, where Nside is the HEALPix resolution parameter. We show that sufficient characterization of the residual noise is unavoidable if one is to draw reliable contraints on large scale anisotropy. Conclusions: We have described how to compute the low-resolution maps, with a controlled sky signal level, and a reliable estimate of covariance of the residual noise. We have also presented a method to smooth the residual noise covariance matrices to describe the noise correlations in smoothed, bandwidth limited maps.

Large-scale structures in turbulent and reverse-transitional sink flow boundary layers

Aspects of large-scale organized structures in sink flow turbulent and reverse-transitional boundary layers are studied experimentally using hot-wire anemometry. Each of the present sink flow boundary layers is in a state of 'perfect equilibrium' or 'exact self-preservation' in the sense of Townsend (The Structure of Turbulent Shear Flow, 1st and 2nd edns, 1956, 1976, Cambridge University Press) and Rotta (Progr. Aeronaut. Sci., vol. 2, 1962, pp. 1-220) and conforms to the notion of 'pure wall-flow' (Coles, J. Aerosp. Sci., vol. 24, 1957, pp. 495-506), at least for the turbulent cases. It is found that the characteristic inclination angle of the structure undergoes a systematic decrease with the increase in strength of the streamwise favourable pressure gradient. Detectable wall-normal extent of the structure is found to be typically half of the boundary layer thickness. Streamwise extent of the structure shows marked increase as the favourable pressure gradient is made progressively severe. Proposals for the typical eddy forms in sink flow turbulent and reverse-transitional flows are presented, and the possibility of structural self-organization (i.e. individual hairpin vortices forming streamwise coherent hairpin packets) in these flows is also discussed. It is further indicated that these structural ideas may be used to explain, from a structural viewpoint, the phenomenon of soft relaminarization or reverse transition of turbulent boundary layers when subjected to strong streamwise favourable pressure gradients. Taylor's 'frozen turbulence' hypothesis is experimentally shown to be valid for flows in the present study even though large streamwise accelerations are involved, the flow being even reverse transitional in some cases. Possible conditions, which are required to be satisfied for the safe use of Taylor's hypothesis in pressure-gradient-driven flows, are also outlined. Measured convection velocities are found to be fairly close to the local mean velocities (typically 90% or more) suggesting that the structure gets convected downstream almost along with the mean flow.

Interannual Variations of Sea Surface Temperature over the Arabian Sea and the Indian Monsoon: A New Perspective

The interannual variation of surface fields over the Arabian Sea and Bay of Bengal are studied using data between 1900 and 1979. It is emphasized that the monthly mean sea surface temperature (SST) over the north Indian Ocean and monsoon rainfall are significantly affected by synoptic systems and other intraseasonal variations. To highlight the interannual signals it is important to remove the large-amplitude high-frequency noise and very low frequency long-term trends, if any. By suitable spatial and temporal averaging of the SST and the rainfall data and by removing the long-term trend from the SST data, we have been able to show that there exists a homogeneous region in the southeastern Arabian Sea over which the March�April (MA) SST anomalies are significantly correlated with the seasonal (June�September) rainfall over India. A potential of this premonsoon signal for predicting the seasonal rainfall over India is indicated. It is shown that the correlation between the SST and the seasonal monsoon rainfall goes through a change of sign from significantly positive with premonsoon SST to very small values with SST during the monsoon season and to significantly negative with SST during the post-monsoon months. For the first time, we have demonstrated that heavy or deficient rainfall years are associated with large-scale coherent changes in the SST (although perhaps of small amplitude) over the north Indian 0cean. We also indicate possible reasons for the apparent lack of persistence of the premonsoon SST anomalies.