Unsteady Incompressible Two-Dimensional and Axisymmetric Turbulent Boundary Layer Flows

The unsteady turbulent incompressible boundary-layer flow over two-dimensional and axisymmetric bodies with pressure gradient has been studied. An eddy-viscosity model has been used to model the Reynolds shear stress. The unsteadiness is due to variations in the free stream velocity with time. The nonlinear partial differential equation with three independent variables governing the flow has been solved using Keller's Box method. The results indicate that the free stram velocity distribution exerts strong influence on the boundary-layer characteristics. The point of zero skin friction is found to move upstream as time increases.

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.

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.

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.

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.

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.

Large Eddy Simulation

Large eddy simulation (LES) is an emerging technique for obtaining an approximation to turbulent flow fields It is an improvement over the widely prevalent practice of obtaining means of turbulent flows when the flow has large scale, low frequency, unsteadiness An introduction to the method, its general formulation, and the more common modelling for flows without reaction, is discussed Some attempts at extension to flows with combustion have been made Examples from present work for flows with and without combustion are given The final example of the LES of the combustor of a helicopter engine illustrates the state-of-the-art in application of the technique

Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

Relationship between the controllability grammian and closed-loop eigenvalues: the single input case

The controllability grammian is important in many control applications. Given a set of closed-loop eigenvalues the corresponding controllability grammian can be obtained by computing the controller which assigns the eigenvalues and then by solving the Lyapunov equation that defines the grammian. The relationship between the controllability grammian, resulting from state feedback, and the closed-loop eigenvalues of a single input linear time invariant (LTI) system is obtained. The proposed methodology does not require the computation of the controller that assigns the specified eigenvalues. The closed-loop system matrix is obtained from the knowledge of the open-loop system matrix, control influence matrix and the specified closed-loop eigenvalues. Knowing the closed-loop system matrix, the grammian is then obtained from the solution of the Lyapunov equation that defines it. Finally the proposed idea is extended to find the state covariance matrix for a specified set of closed-loop eigenvalues (without computing the controller), due to impulsive input in the disturbance channel and to solve the eigenvalue assignment problem for the single input case.

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588]

Numerical Electromagnetic Analysis of Current and Voltage Distribution along Down Conductors and Towers Hit by Direct Lightning

In order to answer the practically important question of whether the down conductors of lightning protection systems to tall towers and buildings can be electrically isolated from the structure itself, this work is conducted. As a first step in this regard, it is presumed that the down conductor placed on metallic tower will be a pessimistic representation of the actual problem. This opinion was based on the fact that the proximity of heavy metallic structure will have a large damping effect. The post-stroke current distributions along the down conductors and towers, which can be quite different from that in the lightning channel, govern the post-stroke near field and the resulting gradient in the soil. Also, for a reliable estimation of the actual stroke current from the measured down conductor currents, it is essential to know the current distribution characteristics along the down conductors. In view of these, the present work attempts to deduce the post-stroke current and voltage distribution along typical down conductors and towers. A solution of the governing field equations on an electromagnetic model of the system is sought for the investigation. Simulation results providing the spatio-temporal distribution of the post-stroke current and voltage has provided very interesting results. It is concluded that it is almost impossible to achieve electrical isolation between the structure and the down conductor. Furthermore, there will be significant induction into the steel matrix of the supporting structure.

Weighted subspace methods and spatial smoothing: analysis and comparison

The effect of using a spatially smoothed forward-backward covariance matrix on the performance of weighted eigen-based state space methods/ESPRIT, and weighted MUSIC for direction-of-arrival (DOA) estimation is analyzed. Expressions for the mean-squared error in the estimates of the signal zeros and the DOA estimates, along with some general properties of the estimates and optimal weighting matrices, are derived. A key result is that optimally weighted MUSIC and weighted state-space methods/ESPRIT have identical asymptotic performance. Moreover, by properly choosing the number of subarrays, the performance of unweighted state space methods can be significantly improved. It is also shown that the mean-squared error in the DOA estimates is independent of the exact distribution of the source amplitudes. This results in a unified framework for dealing with DOA estimation using a uniformly spaced linear sensor array and the time series frequency estimation problems.

Turbulent heat flux variation over the Monsoon-Trough region during MONTBLEX-90

Tower data collected during the Monsoon-Trough Boundary Layer Experiment (MONTBLEX-90) have been analysed to understand the observed structure of the surface layer over an arid region (Jodhpur) and a moist region (Kharagpur) during active and weak phases of the 1990 southwest monsoon. Turbulent heat and momentum fluxes are estimated by the eddy correlation method using sonic data. The turbulent momentum flux at both Jodhpur and Kharagpur was larger when the winds were stronger, reaching a maximum of the order of 0.5 N m(-2) on 5 and 6 August when a low pressure system was located over the region. The heat flux at Jodhpur is high during weak monsoon days, the maximum being 450 W m(-2), whereas during active days the flux never exceeds 200 W m(-2). At Kharagpur, the flux does not vary significantly between active and weak monsoon days, the maximum in either phase being 160 W m(-2) At Jodhpur, there is significant contrast in the near-surface air temperature, being higher during weak monsoon days as compared to active days. Cloud cover did not vary significantly in both the regions. The turbulent heat flux variation at both the sites appears to be correlated mainly with soil mixture, and less sensitive to cloud cover.

Large eddy simulation of a premixed flame with approximate deconvolution modeling

Approximate deconvolution modeling is a very recent approach to large eddy simulation of turbulent flows. It has been applied to compressible flows with success. Here, a premixed flame which forms in the wake of a flameholder has been selected to examine the subgrid-scale modeling of reaction rate by this new method because a previous plane two-dimensional simulation of this wake flame, using a wrinkling function and artificial flame thickening, had revealed discrepancies when compared with experiment. The present simulation is of the temporal evolution of a round wakelike flow at two Reynolds numbers, Re = 2000 and 10,000, based on wake defect velocity and wake diameter. A Fourier-spectral code has been used. The reaction is single-step and irreversible, and the rate follows an Arrhenius law. The reference simulation at the lower Reynolds number is fully resolved. At Re = 10,000, subgrid-scale contributions are significant. It was found that subgrid-scale modeling in the present simulation agrees more closely with unresolved subgrid-scale effects observed in experiment. Specifically, the highest contributions appeared in thin folded regions created by vortex convection. The wrinkling function approach had not selected subgrid-scale effects in these regions.