A numerical study of division of flow in open channels

A two-dimensional numerical model which employs the depth-averaged forms of continuity and momentum equations along with k-e turbulence closure scheme is used to simulate the flow at the open channel divisions. The model is generalised to flows of arbitrary geometries and MacCormack finite volume method is used for solving governing equations. Application of cartesian version of the model to analyse the flow at right-angled junction is presented. The numerical predictions are compared with experimental data of earlier investigators and measurements made as part of the present study. Performance of the model in predicting discharge distribution, surface profiles, separation zone parameters and energy losses is evaluated and discussed in detail. To illustrate the application of the numerical model to analyse the flow in acute angled offtakes and streamlined branch entries, a few computational results are presented.

Transient analysis of mammalian cell growth in hollow fibre bioreactor

A mathematical model describing the dynamics of mammalian cell growth in hollow fibre bioreactor operated in closed shell mode is developed. Mammalian cells are assumed to grow as an expanding biofilm in the extra-capillary space surrounding the fibre. Diffusion is assumed to be the dominant process in the radial direction while axial convection dominates in the lumen of the bioreactor. The transient simulation results show that steep gradients in the cell number are possible under the condition of substrate limitation. The precise conditions which result in nonuniform growth of cells along the length of the bioreactor are delineated. The effect of various operating conditions, such as substrate feed rate, length of the bioreactor and diffusivity of substrate in different regions of the bioreactor, on the bioreactor performance are evaluated in terms of time required to attain the steady-state. The rime of growth is introduced as a measure of effectiveness factor for the bioreactor and is found to be dependent on two parameters, a modified Peclet number and a Thiele modulus. Diffusion, reaction and/or convection control regimes are identified based on these two parameters. The model is further extended to include dual substrate growth limitations, and the relative growth limiting characteristics of two substrates are evaluated. (C) 1997 Elsevier Science Ltd.

A reinforcement learning approach to automatic generation control

This paper formulates the automatic generation control (AGC) problem as a stochastic multistage decision problem. A strategy for solving this new AGC problem formulation is presented by using a reinforcement learning (RL) approach This method of obtaining an AGC controller does not depend on any knowledge of the system model and more importantly it admits considerable flexibility in defining the control objective. Two specific RL based AGC algorithms are presented. The first algorithm uses the traditional control objective of limiting area control error (ACE) excursions, where as, in the second algorithm, the controller can restore the load-generation balance by only monitoring deviation in tie line flows and system frequency and it does not need to know or estimate the composite ACE signal as is done by all current approaches. The effectiveness and versatility of the approaches has been demonstrated using a two area AGC model. (C) 2002 Elsevier Science B.V. All rights reserved.

Study of gas-fines in decreasing gas velocity in a packed bed

It has been well recognized now that the blast furnace aerodynamics can be represented more accurately under the decreasing gas velocity condition. Therefore, gas-fines study has been carried out in a packed bed under the decreasing gas velocity condition. Gas and fines flow equations have been developed and solved, for two-dimensional case using finite volume method. To take into account the turbulence, k-e turbulent flow model has also been developed in two-dimension. The model's predictions have been validated against the published experimental data for the increasing gas velocity case, as no experimental data are available in open literature for the decreasing gas velocity. This study shows the difference in the results for increasing and decreasing gas velocity cases under various conditions which have been reported here. Implication of the results to the blast furnace condition has also been discussed.

Numerical Study Of Heat Transfer From Pin Fin Heat Sink Using Steady And Pulsated Impinging Jets

The work reported in this thesis is an attempt to enhance heat transfer in electronic devices with the use of impinging air jets on pin-finned heat sinks. The cooling per-formance of electronic devices has attracted increased attention owing to the demand of compact size, higher power densities and demands on system performance and re-liability. Although the technology of cooling has greatly advanced, the main cause of malfunction of the electronic devices remains overheating. The problem arises due to restriction of space and also due to high heat dissipation rates, which have increased from a fraction of a W/cm2to 100s of W /cm2. Although several researchers have at-tempted to address this at the design stage, unfortunately the speed of invention of cooling mechanism has not kept pace with the ever-increasing requirement of heat re- moval from electronic chips. As a result, efficient cooling of electronic chip remains a challenge in thermal engineering. Heat transfer can be enhanced by several ways like air cooling, liquid cooling, phase change cooling etc. However, in certain applications due to limitations on cost and weight, eg. air borne application, air cooling is imperative. The heat transfer can be increased by two ways. First, increasing the heat transfer coefficient (forced convec- tion), and second, increasing the surface area of heat transfer (finned heat sinks). From previous literature it was established that for a given volumetric air flow rate, jet im-pingement is the best option for enhancing heat transfer coefficient and for a given volume of heat sink material pin-finned heat sinks are the best option because of their high surface area to volume ratio. There are certain applications where very high jet velocities cannot be used because of limitations of noise and presence of delicate components. This process can further be improved by pulsating the jet. A steady jet often stabilizes the boundary layer on the surface to be cooled. Enhancement in the convective heat transfer can be achieved if the boundary layer is broken. Disruptions in the boundary layer can be caused by pulsating the impinging jet, i.e., making the jet unsteady. Besides, the pulsations lead to chaotic mixing, i.e., the fluid particles no more follow well defined streamlines but move unpredictably through the stagnation region. Thus the flow mimics turbulence at low Reynolds number. The pulsation should be done in such a way that the boundary layer can be disturbed periodically and yet adequate coolant is made available. So, that there is not much variation in temperature during one pulse cycle. From previous literature it was found that square waveform is most effective in enhancing heat transfer. In the present study the combined effect of pin-finned heat sink and impinging slot jet, both steady and unsteady, has been investigated for both laminar and turbulent flows. The effect of fin height and height of impingement has been studied. The jets have been pulsated in square waveform to study the effect of frequency and duty cycle. This thesis attempts to increase our understanding of the slot jet impingement on pin-finned heat sinks through numerical investigations. A systematic study is carried out using the finite-volume code FLUENT (Version 6.2) to solve the thermal and flow fields. The standard k-ε model for turbulence equations and two layer zonal model in wall function are used in the problem Pressure-velocity coupling is handled using the SIMPLE algorithm with a staggered grid. The parameters that affect the heat transfer coefficient are: height of the fins, total height of impingement, jet exit Reynolds number, frequency of the jet and duty cycle (percentage time the jet is flowing during one complete cycle of the pulse). From the studies carried out it was found that: a) beyond a certain height of the fin the rate of enhancement of heat transfer becomes very low with further increase in height, b) the heat transfer enhancement is much more sensitive to any changes at low Reynolds number than compared to high Reynolds number, c) for a given total height of impingement the use of fins and pulsated jet, increases the effective heat transfer coefficient by almost 200% for the same average Reynolds number, d) for all the cases it was observed that the optimum frequency of impingement is around 50 − 100 Hz and optimum duty cycle around 25-33.33%, e) in the case of turbulent jets the enhancement in heat transfer due to pulsations is very less compared to the enhancement in case of laminar jets.

Effect of Electromagnetic Stirring on Macrosegregation during Solidification of a Binary Alloy

The effect of electromagnetic stirring of melt on the final macrosegregation in the continuous casting of an aluminium alloy billet is studied numerically. A continuum mixture model for solidification in presence of electromagnetic stirring is presented. As a case study, simulations are performed for direct chill (DC) casting of an Al-Cu alloy and the effect of electromagnetic stirring on macrosegregation is analysed. The model predicts the temperature, velocity, and species distribution in the mold. As a special case, we have also studied the case in which dendritic particles are fragmented at the interface due to vigorous electromagnetic stirring. For this case, an additional conservation equation for the transport of solid fraction is solved. For modeling the resistance offered by moving solid crystals, a switching function in the momentum equations is used for variation of viscosity. The fragmentation and transport of dendritic particles has a profound effect on the final macrosegregation and microstructure of the solidified billet. It is found that the application of electromagnetic stirring in continuous casting of billets results in better temperature uniformity and macrosegregation pattern.

Fracture analysis of concrete using singular fractal functions with lattice beam network and confirmation with acoustic emission study

In the present study singular fractal functions (SFF) were used to generate stress-strain plots for quasibrittle material like concrete and cement mortar and subsequently stress-strain plot of cement mortar obtained using SFF was used for modeling fracture process in concrete. The fracture surface of concrete is rough and irregular. The fracture surface of concrete is affected by the concrete's microstructure that is influenced by water cement ratio, grade of cement and type of aggregate 11-41. Also the macrostructural properties such as the size and shape of the specimen, the initial notch length and the rate of loading contribute to the shape of the fracture surface of concrete. It is known that concrete is a heterogeneous and quasi-brittle material containing micro-defects and its mechanical properties strongly relate to the presence of micro-pores and micro-cracks in concrete 11-41. The damage in concrete is believed to be mainly due to initiation and development of micro-defects with irregularity and fractal characteristics. However, repeated observations at various magnifications also reveal a variety of additional structures that fall between the `micro' and the `macro' and have not yet been described satisfactorily in a systematic manner [1-11,15-17]. The concept of singular fractal functions by Mosolov was used to generate stress-strain plot of cement concrete, cement mortar and subsequently the stress-strain plot of cement mortar was used in two-dimensional lattice model [28]. A two-dimensional lattice model was used to study concrete fracture by considering softening of matrix (cement mortar). The results obtained from simulations with lattice model show softening behavior of concrete and fairly agrees with the experimental results. The number of fractured elements are compared with the acoustic emission (AE) hits. The trend in the cumulative fractured beam elements in the lattice fracture simulation reasonably reflected the trend in the recorded AE measurements. In other words, the pattern in which AE hits were distributed around the notch has the same trend as that of the fractured elements around the notch which is in support of lattice model. (C) 2011 Elsevier Ltd. All rights reserved.

Top polarisation studies in H(-)t and Wt production

The polarisation of top quarks produced in high energy processes can be a very sensitive probe of physics beyond the Standard Model. The kinematical distributions of the decay products of the top quark can provide clean information on the polarisation of the produced top and thus can probe new physics effects in the top quark sector. We study some of the recently proposed polarisation observables involving the decay products of the top quark in the context of H(-)t and Wt production. We show that the effect of the top polarisation on the decay lepton azimuthal angle distribution, studied recently for these processes at leading order in QCD, is robust with respect to the inclusion of next-to-leading order and parton shower corrections. We also consider the leptonic polar angle, as well as recently proposed energy-related distributions of the top decay products. We construct asymmetry parameters from these observables, which can be used to distinguish the new physics signal from the Wt background and discriminate between different values of tan beta and m(H)- in a general type II two-Higgs doublet model. Finally, we show that similar observables may be useful in separating a Standard Model Wt signal from the much larger QCD induced top pair production background.

Automatic energy-momentum conserving time integrators for hyperelastic waves

An energy-momentum conserving time integrator coupled with an automatic finite element algorithm is developed to study longitudinal wave propagation in hyperelastic layers. The Murnaghan strain energy function is used to model material nonlinearity and full geometric nonlinearity is considered. An automatic assembly algorithm using algorithmic differentiation is developed within a discrete Hamiltonian framework to directly formulate the finite element matrices without recourse to an explicit derivation of their algebraic form or the governing equations. The algorithm is illustrated with applications to longitudinal wave propagation in a thin hyperelastic layer modeled with a two-mode kinematic model. Solution obtained using a standard nonlinear finite element model with Newmark time stepping is provided for comparison. (C) 2012 Elsevier B.V. All rights reserved.

NMR Analysis of Cross Strand Aromatic Interactions in an 8 Residue Hairpin and a 14 Residue Three Stranded beta-Sheet Peptide

Cross strand aromatic interactions between a facing pair of phenylalanine residues in antiparallel beta-sheet structures have been probed using two structurally defined model peptides. The octapeptide Boc-(LFVPPLFV)-P-D-P-L-OMe (peptide 1) favors the beta-hairpin conformation nucleated by the type II' beta-turn formed by the (D)Pro-(L)Pro segment, placing Phe2 and Phe7 side chains in proximity. Two centrally positioned (D)Pro-(L)Pro segments facilitate the three stranded beta-sheet formation in the 14 residue peptide Boc-LFV(D)P(L)PLFVA(D)P(L)PLFV-OMe (peptide 2) in which the Phe2/Phe7 orientations are similar to that in the octapeptide. The anticipated folded conformations of peptides 1 and 2 are established by the delineation of intramolecularly hydrogen bonded NH groups and by the observation of specific cross strand NOEs. The observation of ring current shifted aromatic protons is a diagnostic of close approach of the Phe2 and Phe7 side chains. Specific assignment of aromatic proton resonances using HSQC and HSQC-TOCSY methods allow an analysis of interproton NOEs between the spatially proximate aromatic rings. This approach facilitates specific assignments in systems containing multiple aromatic rings in spectra at natural abundance. Evidence is presented for a dynamic process which invokes a correlated conformational change about the C-alpha-C-beta(chi(1)) bond for the pair of interacting Phe residues. NMR results suggest that aromatic ring orientations observed in crystals are maintained in solution. Anomalous temperature dependence of ring current induced proton chemical shifts suggests that solvophobic effects may facilitate aromatic ring clustering in apolar solvents.

Simulation of Influence of Bilayer Melting on Dynamics and Thermodynamics of Interfacial Water

We investigate the effect of bilayer melting transition on thermodynamics and dynamics of interfacial water using molecular dynamics simulation with the two-phase thermodynamic model. We show that the diffusivity of interface water depicts a dynamic crossover at the chain melting transition following an Arrhenius behavior until the transition temperature. The corresponding change in the diffusion coefficient from the bulk to the interface water is comparable with experimental observations found recently for water near 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles Phys. Chem. Chem. Phys. 13, 7732 (2011)]. The entropy and potential energy of interfacial water show distinct changes at the bilayer melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane. DOI: 10.1103/PhysRevLett.110.018303

Wave propagation in a porous composite beam: Porosity determination, location and quantification

A wave-based method is developed to quantify the defect due to porosity and also to locate the porous regions, in a composite beam-type structure. Wave propagation problem for a porous laminated composite beam is modeled using spectral finite element method (SFEM), based on the modified rule of mixture approach, which is used to include the effect of porosity on the stiffness and density of the composite beam structure. The material properties are obtained from the modified rule of mixture model, which are used in a conventional SFEM to develop a new model for solving wave propagation problems in porous laminated composite beam. The influence of the porosity content on the group speed and also the effect of variation in theses parameters on the time responses are studied first, in the forward problem. The change in the time responses with the change in the porosity of the structure is used as a parameter to find the porosity content in a composite beam. The actual measured response from a structure and the numerically obtained time responses are used for the estimation of porosity, by solving a nonlinear optimization problem. The effect of the length of the porous region (in the propagation direction), on the time responses, is studied. The damage force indicator technique is used to locate the porous region in a beam and also to find its length, using the measured wave propagation responses. (C) 2012 Elsevier Ltd. All rights reserved.

Short- and long-term temporal variations in salinity and the oxygen, carbon and hydrogen isotopic compositions of the Hooghly Estuary water, India

The Hooghly River estuary provides a unique experimental site to understand the effect of monsoonal river discharge on freshwater and seawater mixing. Water samples collected bi-weekly for a duration of 17 months were analyzed for salinity, delta O-18,delta C-13(DIC), as well as delta D to investigate the differential mixing of freshwater and seawater. The differences in salinity and delta O-18 of samples collected during low and high tides on the same day are strongly correlated suggesting a well mixed water column at our sampling site. Low salinity and depleted delta O-18 during monsoon is consistent with increased river discharge as well as high rainfall. We identified different slopes in a delta O-18 versus salinity plot for the estuary water samples collected during monsoon and non-monsoon seasons. This is driven by composition of the freshwater source which is dominated by rainwater during monsoon and rivers during non-monsoon months. Selected delta D analyses of samples indicate that groundwater contributes significantly to the Hooghly Estuary during low rainfall times of the year. delta C-13(DIC) measured in the water recorded low values towards the end of monsoon indicating low productivity (i.e. increased organic respiration) while progressively increasing delta C-13(DIC) values from October till January as well as during some of the pre-monsoon months can be explained by increasing productivity. Very low delta C-13(DIC) (similar to-20%0) suggests involvement of carbon derived from anaerobic oxidation of organics and/or methane with potential contribution from increased anthropogenic water supply. An estimate of seawater incursion into the Hooghly Estuary at different times of the year is obtained by using salinity data in a two-component mixing model. Presence of seawater was found maximum (31-37%) during February till July and lowest (less than or equal to 6%) from September till November. We notice a temporal offset between Ganges River discharge farther upstream at Farakka and salinity variation at the Hooghly Estuary. We believe that this time lag is a result of the physical distance between Farakka and Kakdweep (our sampling location) and put constraints on the travel time of river water during early monsoon. (c) 2012 Published by Elsevier B.V.

Probing the indefinite CP nature of the Higgs boson through decay distributions in the process e(+)e(-) -> t(t)over-bar Phi

The recently discovered scalar resonance at the Large Hadron Collider is now almost confirmed to be a Higgs boson, whose CP properties are yet to be established. At the International Linear Collider with and without polarized beams, it may be possible to probe these properties at high precision. In this work, we study the possibility of probing departures from the pure CP-even case, by using the decay distributions in the process e(+)e(-) -> t (t) over bar Phi, with Phi mainly decaying into a b (b) over bar pair. We have compared the case of a minimal extension of the Standard Model case (model I) with an additional pseudoscalar degree of freedom, with a more realistic case namely the CP-violating two-Higgs doublet model (model II) that permits a more general description of the couplings. We have considered the International Linear Collider with root s = 800 GeV and integrated luminosity of 300 fb(-1). Our main findings are that even in the case of small departures from the CP-even case, the decay distributions are sensitive to the presence of a CP-odd component in model II, while it is difficult to probe these departures in model I unless the pseudoscalar component is very large. Noting that the proposed degrees of beam polarization increase the statistics, the process demonstrates the effective role of beam polarization in studies beyond the Standard Model. Further, our study shows that an indefinite CP Higgs would be a sensitive laboratory to physics beyond the Standard Model.

Remoulded shear strength at plastic and semi-solid states

The undrained shear strength of remoulded soils is of great concern in geotechnical engineering applications. This study aims to develop a reliable approach for determining the undrained shear strength of remoulded fine-grained soils, through the use of index test results, at both the plastic and semi-solid states of consistency. Experimental investigation and subsequent analysis involving a number of fine-grained soils of widely varying plasticity and geological origin have led to a two-parameter linear model of the relationship between logarithm of remoulded undrained shear strength and liquidity index. The numerical values of the parameters are found to be dependent to a lesser extent on the soil group and to a greater extent on the soil state. Based on the values of regression coefficient, ranking index and ranking distance, it seems that the relationship represents the experimental results well. It may be pointed out that the possibility of such a relationship in the semi-solid state of a soil has not been explored in the past. It is also shown that the shear strength at the plastic limit is about 32-34 times that at the liquid limit.