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.

A Theory of the Lifted Temperature Minimum on Calm Clear Nights

Numerous reports from several parts of the world have confirmed that on calm clear nights a minimum in air temperature can occur just above ground, at heights of the order of $\frac{1}{2}$ m or less. This phenomenon, first observed by Ramdas & Atmanathan (1932), carries the associated paradox of an apparently unstable layer that sustains itself for several hours, and has not so far been satisfactorily explained. We formulate here a theory that considers energy balance between radiation, conduction and free or forced convection in humid air, with surface temperature, humidity and wind incorporated into an appropriate mathematical model as parameters. A complete numerical solution of the coupled air-soil problem is used to validate an approach that specifies the surface temperature boundary condition through a cooling rate parameter. Utilizing a flux-emissivity scheme for computing radiative transfer, the model is numerically solved for various values of turbulent friction velocity. It is shown that a lifted minimum is predicted by the model for values of ground emissivity not too close to unity, and for sufficiently low surface cooling rates and eddy transport. Agreement with observation for reasonable values of the parameters is demonstrated. A heuristic argument is offered to show that radiation substantially increases the critical Rayleigh number for convection, thus circumventing or weakening Rayleigh-Benard instability. The model highlights the key role played by two parameters generally ignored in explanations of the phenomenon, namely surface emissivity and soil thermal conductivity, and shows that it is unnecessary to invoke the presence of such particulate constituents as haze to produce a lifted minimum.

Excitation and relaxation energies of trans-stilbene: Confined singlet, triplet, and charged bipolarons

The π-electronic excitations and excited-state geometries of trans-stilbene (tS) are found by combining exact solutions of the Pariser-Parr-Pople (PPP) model and semiempirical Parametric Method 3 (PM3) calculations. Comprehensive comparisons with tS spectra are obtained and related to the fluorescence and topological alternation of poly(paraphenylenevinylene) (PPV). The one-photon absorption and triplet of tS correspond, respectively, to singlet and triplet bipolarons confined to two phenyls, while the tS2- ground state is a confined charged bipolaron. Independent estimates of the relaxation energy between vertical and adiabatic excitation show the bipolaron binding energy to depend on both charge and spin, as expected for interacting π electrons in correlated or molecular states. Complete configuration interaction within the PPP model of tS accounts for the singlet-triplet gap, for the fine-structure constants and triplet-triplet spectra, for two-photon transitions and intensities, and for one-photon spectra and the radiative lifetime, although the relative position of nearly degenerate covalent and ionic singlets is not resolved. The planar PM3 geometry and low rotational barrier of tS agree with resolved rotational and vibrational spectra in molecular beams. PM3 excitation and relaxation energies for tS bipolarons are consistent with experiment and with PPP results. Instead of the exciton model, we interpret tS excitations in terms of states that are localized on each ring or extended over an alternating chain, as found exactly in Hückel theory, and find nearly degenerate transitions between extended and localized states in the singlet, triplet, and dianion manifolds. The large topological alternation of the extended system increases the ionicity and interchanges the order of the lowest one- and two-photon absorption of PPV relative to polyenes.

Solvation and axial ligation properties of (2,3,7,8,12,13,17,18-octabromo-5,10,15,20-tetraphenylporphyrinato)zinc(II)

he solvation of (2,3,7,8,12,13,17,18-octabromo-5,10,15,20-tetraphenylporphyrinato)zinc(II)[Zn(obtpp)], in twelve different solvents results in large red shifts of the B and Q bands of the porphyrin accompanied by enhanced absorbance ratios of the Q bands. These observations are ascribed to the destabilisation of the highest occupied molecular orbital a2u of the porphyrin arising from a flow of charge from the axial ligand to the porphyrin ring through the zinc(II) ion. The binding constants of adducts of [Zn(obtpp)] with neutral bases have been found to be an order of magnitude greater than those observed for the corresponding adducts of (5,10,15,20-tetraphenylporphyrinato)-zinc and vary in the order piperidine > imidazole > pyridine > 3-methylpyridine > pyridine-3-carbaldehyde. The enhanced binding constants and large spectral shifts are interpreted in terms of the electrophilicity of [Zn(obtpp)] induced by the electron-withdrawing bromine substituents in the porphyrin core. The structure of [Zn(obtpp)(PrCN)2] has been determined; it reveals six-co-ordinated zinc(II) with two long Zn–N distance [2.51(4), 2.59(3)Å]. The porphyrin is non-planar and displays a saddle-shaped conformation.

The Reproducing Kernel DMS-FEM: 3D Shape Functions and Applications to Linear Solid Mechanics

We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and ID NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of polynomial reproduction whilst deriving the shape functions. Nevertheless, given the higher complexity in forming the knotclouds for tetrahedral elements especially when higher demand is placed on the order of continuity of the shape functions across inter-element boundaries, we presently emphasize an exploration of the triangular prism based formulation in the context of several benchmark problems of interest in linear solid mechanics. In the absence of a more rigorous study on the convergence analyses, the numerical exercise, reported herein, helps establish the method as one of remarkable accuracy and robust performance against numerical ill-conditioning (such as locking of different kinds) vis-a-vis the conventional FEM.

Top-spin analysis of new scalar and tensor interactions in e(+)e(-) collisions with beam polarization

We utilize top polarization in the process e(+)e(-) -> t (t) over bar at the International Linear Collider ( ILC) with transverse beam polarization to probe interactions of the scalar and tensor type beyond the standard model and to disentangle their individual contributions. Ninety percent confidence level limits on the interactions with realistic integrated luminosity are presented and are found to improve by an order of magnitude compared to the case when the spin of the top quark is not measured. Sensitivities of the order of a few times 10(-3) TeV-2 for real and imaginary parts of both scalar and tensor couplings at root s = 500 and 800 GeV with an integrated luminosity of 500 fb(-1) and completely polarized beams are shown to be possible. A powerful model-independent framework for inclusive measurements is employed to describe the spin-momentum correlations, and their C, P, and T properties are presented in a technical appendix.

Transformation formula for exponential sums involving fourier coefficients of modular forms

In 1984 Jutila [5] obtained a transformation formula for certain exponential sums involving the Fourier coefficients of a holomorphic cusp form for the full modular group SL(2, Z). With the help of the transformation formula he obtained good estimates for the distance between consecutive zeros on the critical line of the Dirichlet series associated with the cusp form and for the order of the Dirichlet series on the critical line, [7]. In this paper we follow Jutila to obtain a transformation formula for exponential sums involving the Fourier coefficients of either holomorphic cusp forms or certain Maass forms for congruence subgroups of SL(2, Z) and prove similar estimates for the corresponding Dirichlet series.

Dye sensitized visible light degradation of phenolic compounds

The present research work reports the eosin Y (EY) and fluorescein (FL) sensitized visible light degradation of phenol, 4-chlorophenol (CP), 2,4-dichlorophenol (DCP) and 2,4,6-trichlorophenol (TCP) using combustion synthesized nano-TiO2 (CS TiO2). The rate of degradation of the phenolic compounds was higher in the presence of EY/CS TiO2 compared to FL/CS TiO2 system. A detailed mechanism of sensitized degradation was proposed and a mechanistic model for the rate of degradation of the phenolic compound was derived using the pyramidal network reduction technique. It was found that at low initial dye concentrations, the rate of degradation of the phenolic compound is first order in the concentration of the dye, while at high initial dye concentrations, the rate is first order in the concentration of the phenolic compound. The order of degradation of the different phenolic compounds follows: CP > TCP > DCP > phenol. The different phenolic and dye intermediates that were formed during the degradation were identified by liquid chromatography-mass spectrometry (LC-MS) and the most probable pathway of degradation is proposed. (C) 2010 Elsevier B.V. All rights reserved.

Mechanical properties of B-modified Ti-6Al-4V alloys

Minor addition of B to the Ti-6Al-4V alloy reduces the prior beta grain size by more than an order of magnitude. TiB formed in-situ in the process has been noted to decorate the grain boundaries. This microstructural modification influences the mechanical behavior of the Ti-6Al-4V alloy significantly. In this paper, an overview of our current research on tensile properties, fracture toughness as well as notched and un-notched fatigue properties of Ti-6Al-4V-xB with x varying between 0.0 to 0.55 wt.% is presented. A quantitative relationship between the microstructural length scales and the various mechanical properties have been developed. Moreover, the effect of the presence of hard and brittle TiB has also been studied.

Copper( 11) Complexes of Novel Tripodal Ligands containing Phenolate and Benzimidazole/Pyridine Pendants: Synthesis, Structure, Spectra and Electrochemical Behaviour

Mononuclear copper(II) complexes of tri- and tetra-dentate tripodal ligands containing phenolic hydroxyl and benzimidazole or pyridine groups have been isolated. They are of the type (CuL(X)].nH2O, [CuL(H2O)]X.nH2O or [CuL].nH2O where X = Cl-, ClO4-, N3- or NCS- and n = 0-4. The electronic spectra of all the complexes exhibit a broad absorption band around 14000 cm-1 and the polycrystalline as well as the frozen-solution EPR spectra are axial, indicating square-based geometries. The crystal structure of [CuL(Cl)] [HL = (2-hydroxy-5-nitrobenzyl)bis(2-pyridyl-methyl)amine] revealed a square-pyramidal geometry around Cu(II). The mononuclear complex crystallises in the triclinic space group P1BAR with a = 6.938(1), b = 11.782(6), c = 12.678(3) angstrom and alpha = 114.56(3), beta = 92.70(2), gamma = 95.36(2)-degrees. The co-ordination plane is comprised of one tertiary amine and two pyridine nitrogens and a chloride ion. The phenolate ion unusually occupies the axial site, possibly due to the electron-withdrawing p-nitro group. The enhanced pi delocalisation involving the p-nitrophenolate donor elevates the E1/2 values. The spectral and electrochemical results suggest the order of donor strength as nitrophenolate < pyridine < benzimidazole in the tridentate and nitrophenolate < benzimidazole < pyridine in the tetradentate ligand complexes.

Data compression by linear prediction for storage and transmission of EEG signals

The EEG time series has been subjected to various formalisms of analysis to extract meaningful information regarding the underlying neural events. In this paper the linear prediction (LP) method has been used for analysis and presentation of spectral array data for the better visualisation of background EEG activity. It has also been used for signal generation, efficient data storage and transmission of EEG. The LP method is compared with the standard Fourier method of compressed spectral array (CSA) of the multichannel EEG data. The autocorrelation autoregressive (AR) technique is used for obtaining the LP coefficients with a model order of 15. While the Fourier method reduces the data only by half, the LP method just requires the storage of signal variance and LP coefficients. The signal generated using white Gaussian noise as the input to the LP filter has a high correlation coefficient of 0.97 with that of original signal, thus making LP as a useful tool for storage and transmission of EEG. The biological significance of Fourier method and the LP method in respect to the microstructure of neuronal events in the generation of EEG is discussed.

Effects of Ultrafast Solvation on the Rate of Adiabatic Outer-Sphere Electron Transfer Reactions

Recent studies have demonstrated that solvation dynamics in many common dipolar liquids contain an initial, ultrafast Gaussian component which may contribute even more than 60% to the total solvation energy. It is also known that adiabatic electron transfer reactions often probe the high-frequency components of the relevant solvent friction (Hynes, J. T. J. Phys. Chem. 1986, 90, 3701). In this paper, we present a theoretical study of the effects of the ultrafast solvent polar modes on the adiabatic electron transfer reactions by using the formalism of Hynes. Calculations have been carried out for a model system and also for water and acetonitrile. It is found that, in general, the ultrafast modes can greatly enhance the rate of electron transfer, even by more than an order of magnitude, over the rate obtained by using only the slow overdamped modes usually considered. For water, this acceleration of the rate can be attributed to the high-frequency intermolecular vibrational and librational modes. For a weakly adiabatic reaction, the rate is virtually indistinguishable from the rate predicted by the Marcus transition state theory. Another important result is that even in this case of ultrafast underdamped solvation, energy diffusion appears to be efficient so that electron transfer reaction in water is controlled essentially by the barrier crossing dynamics. This is because the reactant well frequency is-directly proportional to the rate of the initial Gaussian decay of the solvation time correlation function. As a result, the value of the friction at the reactant well frequency rarely falls below the value required for the Kramers turnover except when the polarizability of the water molecules may be neglected. On the other hand, in acetonitrile, the rate of electron transfer reaction is found to be controlled by the energy diffusion dynamics, although a significant contribution to the rate comes also from the barrier crossing rate. Therefore, the present study calls for a need to understand the relaxation of the high-frequency modes in dipolar liquids.

Characterization of the DNA-binding domain of beta protein, a component of phage lambda Red-pathway by UV catalyzed cross-linking

beta protein, a key component of Red-pathway of phage lambda is necessary for its growth and general genetic recombination in recombination-deficient mutants of Escherichia coli. To facilitate studies on structure-function relationships, we overexpressed beta protein and purified it to homogeneity. A chemical cross-linking reagent, glutaraldehyde, was used to stabilize the physical association of beta protein in solution. A 67-kDa band, corresponding to homodimer, was identified after separation by SDS-polyacrylamide gel electrophoresis. Stoichiometric measurements indicated a site-size of 1 monomer of beta protein/5 nucleotide residues. Electrophoretic gel mobility shift assays suggested that beta protein formed stable nucleoprotein complexes with 36-mer, but not with 27- or 17-mer DNA. Interestingly, the interaction of beta protein with DNA and the stability of nucleoprotein complexes was dependent on the presence of MgCl2, and the binding was abolished by 250 mM NaCl. The K-d of beta protein binding to 36-mer DNA was on the order of 1.8 x 10(-6) M. Photochemical cross-linking of native beta protein or its fragments, generated by chymotrypsin, to 36-mer DNA was performed to identify its DNA-binding domain. Characterization of the cross-linked peptide disclosed that amino acids required for DNA-binding specificity resided within a 20-kDa peptide at the N-terminal end. These findings provide a basis for further understanding oi the structure and function of beta protein.

Experimental validation of the 1-D acoustical model for conical concentric tube resonators with moving medium

Variable cross-sectional area ducts are often used for attenuation at lower frequencies (of the order of firing frequency), whereas concentric tube resonators provide attenuation at relatively higher frequencies. In this paper, analysis of one dimensional control volume approach of conical concentric tube resonators is validated experimentally. The effects of mean flow and taper are investigated. The experimental setup is specially designed to measure the pressure transfer function in the form of Level Difference or Noise Reduction across the test muffler. It is shown that there is a reasonably good agreement between the predicted values of the Noise Reduction and the measured ones for incompressible mean flow as well as stationary medium. (C) 2011 Institute of Noise Control Engineering.

Field evaluation of unsaturated hydraulic conductivity models and parameter estimation from retention data

Predictions of two popular closed-form models for unsaturated hydraulic conductivity (K) are compared with in situ measurements made in a sandy loam field soil. Whereas the Van Genuchten model estimates were very close to field measured values, the Brooks-Corey model predictions were higher by about one order of magnitude in the wetter range. Estimation of parameters of the Van Genuchten soil moisture characteristic (SMC) equation, however, involves the use of non-linear regression techniques. The Brooks-Corey SMC equation has the advantage of being amenable to application of linear regression techniques for estimation of its parameters from retention data. A conversion technique, whereby known Brooks-Corey model parameters may be converted into Van Genuchten model parameters, is formulated. The proposed conversion algorithm may be used to obtain the parameters of the preferred Van Genuchten model from in situ retention data, without the use of non-linear regression techniques.