Representation of thermodynamic properties of ternary systems and its application to the system silver-gold-copper at 1350 K

The article presents a generalized analytical expression for description of the integral excess Gibbs free energy of mixing of a ternary system. Twelve constants of the equation are assessed by the least mean squares regressional analysis of the experimental integral excess data of the constituent binaries; three ternary parameters are evaluated by a regressional analysis based on the partial experimental data of a component of the ternary system. The assessed values of the ternary parameters describe the nature of the ternary interaction in the system. Activities and isoactivities of the components in the Ag-Au-Cu system at 1350 K are calculated and found to be in good agreement with the experimental data. This analytical treatment is particularly useful to ternary systems where the thermodynamic data are available from different sources.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Transformation of a laterally diverging boundary layer flow to a two-dimensional boundary layer flow

Laterally diverging boundary layer flow over a plate is shown to be reducible to a two-dimensional flow by modelling the diverging streamlines by a source flow.

Unsteady compressible boundary layer flow in the stagnation region of a sphere with a magnetic field

Abstract: An analysis is performed to study the unsteady compressible laminar boundary layer flow in the forward stagnation-point region of a sphere with a magnetic field applied normal, to the surface. We have considered the case where there is an initial steady state that is perturbed by the step change in the total enthalpy at the wall. The nonlinear coupled parabolic partial differential equations governing the flow and heat transfer have been solved numerically using a finite-difference scheme. The numerical results are presented, which show the temporal development of the boundary layer. The magnetic field in the presence of variable electrical conductivity causes an overshoot in the velocity profile. Also, when the total enthalpy at the wall is suddenly increased, there is a change in the direction of transfer of heat in a small interval of time.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

Docking and free energy simulations to predict conformational domains involved in hCG-LH receptor interactions using recombinant antibodies

Single chain fragment variables (ScFvs) have been extensively employed in studying the protein-protein interactions. ScFvs derived from phage display libraries have an additional advantage of being generated against a native antigen, circumventing loss of information on conformational epitopes. In the present study, an attempt has been made to elucidate human chorionic gonadotropin (hCG)-luteinizing hormone (LH) receptor interactions by using a neutral and two inhibitory ScFvs against hCG. The objective was to dock a computationally derived model of these ScFvs onto the crystal structure of hCG and understand the differential roles of the mapped epitopes in hCG-LH receptor interactions. An anti-hCG ScFv, whose epitope was mapped previously using biochemical tools, served as the positive control for assessing the quality of docking analysis. To evaluate the role of specific side chains at the hCG-ScFv interface, binding free energy as well as residue interaction energies of complexes in solution were calculated using molecular mechanics Poisson-Boltzmann/surface area method after performing the molecular dynamic simulations on the selected hCG-ScFv models and validated using biochemical and SPR analysis. The robustness of these calculations was demonstrated by comparing the theoretically determined binding energies with the experimentally obtained kinetic parameters for hCG-ScFv complexes. Superimposition of hCG-ScFv model onto a model of hCG complexed with the 51-266 residues of LH receptor revealed importance of the residues previously thought to be unimportant for hormone binding and response. This analysis provides an alternate tool for understanding the structure-function analysis of ligand-receptor interactions. Proteins 2011;79:3108-3122. (C) 2011 Wiley-Liss, Inc.

Variation of the gas and radiation content in the sub-Keplerian accretion disk around black holes and its impact to the solutions

We investigate the variation of the gas and the radiation pressure in accretion disks during the infall of matter to the black hole and its effect to the flow. While the flow far away from the black hole might be non-relativistic, in the vicinity of the black hole it is expected to be relativistic behaving more like radiation. Therefore, the ratio of gas pressure to total pressure (beta) and the underlying polytropic index (gamma) should not be constant throughout the flow. We obtain that accretion flows exhibit significant variation of beta and then gamma, which affects solutions described in the standard literature based on constant beta. Certain solutions for a particular set of initial parameters with a constant beta do not exist when the variation of beta is incorporated appropriately. We model the viscous sub-Keplerian accretion disk with a nonzero component of advection and pressure gradient around black holes by preserving the conservations of mass, momentum, energy, supplemented by the evolution of beta. By solving the set of five coupled differential equations, we obtain the thermo-hydrodynamical properties of the flow. We show that during infall, beta of the flow could vary up to similar to 300%, while gamma up to similar to 20%. This might have a significant impact to the disk solutions in explaining observed data, e.g. super-luminal jets from disks, luminosity, and then extracting fundamental properties from them. Hence any conclusion based on constant gamma and beta should be taken with caution and corrected. (C) 2011 Elsevier B.V. All rights reserved.

Relevance of local parallel theory to the linear stability of laminar separation bubbles

Laminar separation bubbles are thought to be highly non-parallel, and hence global stability studies start from this premise. However, experimentalists have always realized that the flow is more parallel than is commonly believed, for pressure-gradient-induced bubbles, and this is why linear parallel stability theory has been successful in describing their early stages of transition. The present experimental/numerical study re-examines this important issue and finds that the base flow in such a separation bubble becomes nearly parallel due to a strong-interaction process between the separated boundary layer and the outer potential flow. The so-called dead-air region or the region of constant pressure is a simple consequence of this strong interaction. We use triple-deck theory to qualitatively explain these features. Next, the implications of global analysis for the linear stability of separation bubbles are considered. In particular we show that in the initial portion of the bubble, where the flow is nearly parallel, local stability analysis is sufficient to capture the essential physics. It appears that the real utility of the global analysis is perhaps in the rear portion of the bubble, where the flow is highly non-parallel, and where the secondary/nonlinear instability stages are likely to dominate the dynamics.

Computational Study of the Migration of Rhenium from One Enantioface of an Olefin to the Other Facilitated by (C-H)center dot center dot center dot Re Interactions

The migration of a metal atom in a metal olefin complex from one pi face of the olefin to the opposite pi face has been rarely documented. Gladysz and co-workers showed that such a movement is indeed possible in monosubstituted chiral Re olefin complexes, resulting in diastereomerization. Interestingly, this isomerization occurred without dissociation, and on the basis of kinetic isotope effects, the involvement of a trans C-H bond was indicated. Either oxidative addition or an agostic interaction of the vinylic C-H(D) bond with the metal could account for the experimentally observed kinetic isotope effect. In this study we compute the free energy of activation for the migration of Re from one enantioface of the olefin to the other through various pathways. On the basis of DFT calculations at the B3LYP level we show that a trans (C-H)center dot center dot center dot Re interaction and trans C-H oxidative addition provide a nondissociative path for the diastereomerization. The trans (C-H)center dot center dot center dot Re interaction path is computed to be more favorable by 2.3 kcal mol(-1) than the oxidative addition path. While direct experimental evidence was not able to discount the migration of the metal through the formation of a eta(2)-arene complex (conducted tour mechanism), computational results at the B3LYP level show that it is energetically more expensive. Surprisingly, a similar analysis carried out at the M06 level computes a lower energy path for the conducted tour mechanism and is not consistent with the experimental isotope effects observed. Metal-(C-H) interactions and oxidative additions of the metal into C-H bonds are closely separated in energy and might contribute to unusual fluxional processes such as this diastereomerization.

Two-pion low-energy contribution to the muon g-2 with improved precision from analyticity and unitarity

The two-pion contribution from low energies to the muon magnetic moment anomaly, although small, has a large relative uncertainty since in this region the experimental data on the cross sections are neither sufficient nor precise enough. It is therefore of interest to see whether the precision can be improved by means of additional theoretical information on the pion electromagnetic form factor, which controls the leading-order contribution. In the present paper, we address this problem by exploiting analyticity and unitarity of the form factor in a parametrization-free approach that uses the phase in the elastic region, known with high precision from the Fermi-Watson theorem and Roy equations for pi pi elastic scattering as input. The formalism also includes experimental measurements on the modulus in the region 0.65-0.70 GeV, taken from the most recent e(+)e(-) ->pi(+)pi(-) experiments, and recent measurements of the form factor on the spacelike axis. By combining the results obtained with inputs from CMD2, SND, BABAR, and KLOE, we make the predictions a(mu)(pi pi,LO)2m(pi), 0.30 GeV] = (0.553 +/- 0.004) x 10(-10) and a(mu)(pi pi,LO)0.30 GeV; 0.63 GeV] = (133.083 +/- 0.837) x 10(-10). These are consistent with the other recent determinations and have slightly smaller errors.

Experimental studies on the flow through soft tubes and channels

Experiments conducted in channels/tubes with height/diameter less than 1 mm with soft walls made of polymer gels show that the transition Reynolds number could be significantly lower than the corresponding value of 1200 for a rigid channel or 2100 for a rigid tube. Experiments conducted with very viscous fluids show that there could be an instability even at zero Reynolds number provided the surface is sufficiently soft. Linear stability studies show that the transition Reynolds number is linearly proportional to the wall shear modulus in the low Reynolds number limit, and it increases as the 1/2 and 3/4 power of the shear modulus for the `inviscid' and `wall mode' instabilities at high Reynolds number. While the inviscid instability is similar to that in the flow in a rigid channel, the mechanisms of the viscous and wall mode instabilities are qualitatively different. These involve the transfer of energy from the mean flow to the fluctuations due to the shear work done at the interface. The experimental results for the viscous instability mechanism are in quantitative agreement with theoretical predictions. At high Reynolds number, the instability mechanism has characteristics similar to the wall mode instability. The experimental transition Reynolds number is smaller, by a factor of about 10, than the theoretical prediction for the parabolic flow through rigid tubes and channels. However, if the modification in the tube shape due to the pressure gradient, and the consequent modification in the velocity profile and pressure gradient, are incorporated, there is quantitative agreement between theoretical predictions and experimental results. The transition has important practical consequences, since there is a significant enhancement of mixing after transition.

Application of Similarity Theory to the Characterization of Non-Transferred Laminar Plasma Jet Generation

Laminar-flow non-transferred DC plasma jets were generated by a torch with an inter-electrode insert by which the arc column was limited to a length of about 20 mm. Current–voltage characteristics, thermal efficiency and jet length, a parameter which changes greatly with the generating parameters in contrast with the almost unchangeable jet length of the turbulent plasma, were investigated systematically, by using the similarity theory combined with the corresponding experimental examination. Formulae in non-dimensional forms were derived for predicting the characteristics of the laminar plasma jet generation, within the parameter ranges where no transfer to turbulent flow occurs. Mean arc temperature in the torch channel and mean jet-flow temperature at the torch exit were obtained, and the results indicate that the thermal conductivity feature of the working gas seems to be an important factor affecting thermal efficiency of laminar plasma generation.