Exact free-surface flows for shallow-water equations .1. - the incompressible case

A comprehensive exact treatment of free surface flows governed by shallow water equations (in sigma variables) is given. Several new families of exact solutions of the governing PDEs are found and are shown to embed the well-known self-similar or traveling wave solutions which themselves are governed by reduced ODEs. The classes of solutions found here are explicit in contrast to those found earlier in an implicit form. The height of the free surface for each family of solutions is found explicitly. For the traveling or simple wave, the free surface is governed by a nonlinear wave equation, but is arbitrary otherwise. For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed; in another case, the free surface is a horizontal plane while the flow underneath is a sine wave. The existence of simple waves on shear flows is analytically proved. The interaction of large amplitude progressive waves with shear flow is also studied.

Exact free-surface flows for shallow-water equations .2. the compressible case

Exact free surface flows with shear in a compressible barotropic medium are found, extending the authors' earlier work for the incompressible medium. The barotropic medium is of finite extent in the vertical direction, while it is infinite in the horizontal direction. The ''shallow water'' equations for a compressible barotropic medium, subject to boundary conditions at the free surface and at the bottom, are solved in terms of double psi-series, Simple wave and time-dependent solutions are found; for the former the free surface is of arbitrary shape while for the latter it is a damping traveling wave in the horizontal direction, For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed. In the case of an isothermal medium, when gamma = 1, we again find simple wave and time-dependent solutions.

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.

The relationship between classification of multi-domain proteins using an alignment-free approach and their functions: a case study with immunoglobulins

Establishing functional relationships between multi-domain protein sequences is a non-trivial task. Traditionally, delineating functional assignment and relationships of proteins requires domain assignments as a prerequisite. This process is sensitive to alignment quality and domain definitions. In multi-domain proteins due to multiple reasons, the quality of alignments is poor. We report the correspondence between the classification of proteins represented as full-length gene products and their functions. Our approach differs fundamentally from traditional methods in not performing the classification at the level of domains. Our method is based on an alignment free local matching scores (LMS) computation at the amino-acid sequence level followed by hierarchical clustering. As there are no gold standards for full-length protein sequence classification, we resorted to Gene Ontology and domain-architecture based similarity measures to assess our classification. The final clusters obtained using LMS show high functional and domain architectural similarities. Comparison of the current method with alignment based approaches at both domain and full-length protein showed superiority of the LMS scores. Using this method we have recreated objective relationships among different protein kinase sub-families and also classified immunoglobulin containing proteins where sub-family definitions do not exist currently. This method can be applied to any set of protein sequences and hence will be instrumental in analysis of large numbers of full-length protein sequences.

Microstructural characterization of ultrafine-grain interstitial-free steel by X-ray diffraction line profile analysis

This paper highlights the microstructural features of commercially available interstitial free (IF) steel specimens deformed by equal channel angular pressing (ECAP) up to four passes following the route A. The microstructure of the samples was studied by different techniques of X-ray diffraction peak profile analysis as a function of strain (epsilon). It was found that the crystallite size is reduced substantially already at epsilon=2.3 and it does not change significantly during further deformation. At the same time, the dislocation density increases gradually up to epsilon=4.6. The dislocation densities estimated from X-ray diffraction study are found to correlate very well with the experimentally obtained yield strength of the samples.

Singularity-free path planning for the Stewart platform manipulator

This paper addresses the problem of singularity-free path planning for the six-degree-of-freedom parallel manipulator known as the Stewart platform manipulator. Unlike serial manipulators, the Stewart platform possesses singular configurations within the workspace where the manipulator is uncontrollable. An algorithm has been developed to construct continuous paths within the workspace of the manipulator by avoiding singularities and ill-conditioning. Given two end-poses of the manipulator, the algorithm finds out safe (well-conditioned) via points and plans a continuous path from the initial pose to the final one. When the two end-poses belong to different branches and no singularity-free path is possible, the algorithm indicates the impossibility of a valid path. A numerical example has also been presented as illustration of the path planning strategy.

Unsteady free convection flow in the stagnation-point region of a rotating sphere

The unsteady free convection boundary-layer flow in the forward stagnation-point region of a sphere, which is rotating with time-dependent angular velocity in an ambient fluid, has been studied. Both constant wall temperature and constant hear flux conditions have been considered. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The skin friction and the heat transfer are enhanced by the buoyancy force. The effect of the buoyancy force is found to be more pronounced for smaller Prandtl numbers than for larger Prandtl numbers. For a given buoyancy force, the heat transfer increases with an increase in Prandtl number, but the skin friction decreases.

Unsteady free convection flow in the stagnation-point region of a three-dimensional body

The unsteady free convection flow in the stagnation-point region of a heated three-dimensional body placed in an ambient fluid is studied under boundary layer approximations. We have considered the case where there is an initial steady state that is perturbed by a step-change in the wall temperature. The non-linear coupled partial differential equations governing the free convection flow are solved numerically using a finite difference scheme. The presented results show the temporal development of the momentum and thermal boundary layer characteristics.

Free Energy Gap Dependence of the Electron-Transfer Rate from the Inverted to the Normal Region

numerical study of the free energy gap (FEG) dependence of the electron-transfer rate in polar solvents is presented. This study is based on the generalized multidimensional hybrid model, which not only includes the solvent polarization and the molecular vibration modes, but also the biphasic polar response of the solvent. The free energy gap dependence is found to be sensitive to several factors, including the solvent relaxation rate, the electronic coupling between the surfaces, the frequency of the high-frequency quantum vibrational mode, and the magnitude of the solvent reorganization energy. It is shown that in some cases solvent relaxation can play an important role even in the Marcus normal regime. The minimal hybrid model involves a large number of parameters, giving rise to a diverse non-Marcus FEG behavior which is often determined collectively by these parameters. The model gives the linear free energy gap dependence of the logarithmic rate over a substantial range of FEG, spanning from the normal to the inverted regime. However, even for favorable values of the relevant parameters, a linear free energy gap dependence of the rate could be obtained only over a range of 5000-6000 cm(-1) (compared to the experimentally observed range of 10000 cm(-1) reported by Benniston et al.). The present work suggests several extensions/generalizations of the hybrid model which might be necessary to fully understand the observed free energy gap dependence.

Free energy landscape of a dense hard-sphere system

The topography of the free energy landscape in phase space of a dense hard-sphere system characterized by a discretized free energy functional of the Ramakishnan-Yussouff form is investigated numerically using a specially devised Monte Carlo procedure. We locate a considerable number of glassy local minima of the free energy and analyze the distributions of the free energy at a minimum and an appropriately defined phase-space "distance" between different minima. We find evidence for the existence of pairs of closely related glassy minima("two-level systems"). We also investigate the way the system makes transitions as it moves from the basin of attraction of a minimum to that of another one after a start under nonequilibrium conditions. This allows us to determine the effective height of free energy barriers that separate a glassy minimum from the others. The dependence of the height of free energy barriers on the density is investigated in detail. The general appearance of the free energy landscape resembles that of a putting green: relatively deep minima separated by a fairly flat structure. We discuss the connection of our results with the Vogel-Fulcher law and relate our observations to other work on the glass transition.

Free vibration of a kinked cantilever with attached masses

In this paper free vibration characteristics of a centrally kinked cantilever beam of unit mass carrying masses at the kink (m(k)) and at the tip (m(t)) are analyzed. Frequency factors are presented for the first two modes for different combinations of m(k),m(t) and the kink angle delta. A relationship of the form f(m(k),m(t), delta) = m(k) + m(t)(4 + 10/3 cos delta+ 2/3 cos(2) delta)=const appears to give the same fundamental frequency for a given delta and different combinations of [m(k), m(t)]. Mode shapes as well as bending moments at the support and at the kink are also discussed. The utility of a discrete beam model in understanding the free vibration characteristics is also highlighted.

Extended Kalman filters using explicit and derivative-free local linearizations

We propose three variants of the extended Kalman filter (EKF) especially suited for parameter estimations in mechanical oscillators under Gaussian white noises. These filters are based on three versions of explicit and derivative-free local linearizations (DLL) of the non-linear drift terms in the governing stochastic differential equations (SDE-s). Besides a basic linearization of the non-linear drift functions via one-term replacements, linearizations using replacements through explicit Euler and Newmark expansions are also attempted in order to ensure higher closeness of true solutions with the linearized ones. Thus, unlike the conventional EKF, the proposed filters do not need computing derivatives (tangent matrices) at any stage. The measurements are synthetically generated by corrupting with noise the numerical solutions of the SDE-s through implicit versions of these linearizations. In order to demonstrate the effectiveness and accuracy of the proposed methods vis-à-vis the conventional EKF, numerical illustrations are provided for a few single degree-of-freedom (DOF) oscillators and a three-DOF shear frame with constant parameters.

Controversy on the Free Energy of Formation of CaO Additional Evidence in Support of Thermochemical Data

The standard free energies of formation of CaO derived from a variety of high-temperature equilibrium measurements made by seven groups of experimentalists are significantly different from those given in the standard compilations of thermodynamic data. Indirect support for the validity of the compiled data comes from new solid-state electrochemical measurements using single-crystal CaF2 and SrF2 as electrolytes. The change in free energy for the following reactions are obtained: CaO + MgF2 --> MgO + CaF2 Delta G degrees = -68,050 -2.47 T(+/-100) J mol(-1) SrO + CaF2 --> SrF2 + CaO Delta G degrees = -35,010 + 6.39 T (+/-80) J mol(-1) The standard free energy changes associated with cell reactions agree with data in standard compilations within +/- 4 kJ mol(-1). The results of this study do not support recent suggestions for a major revision in thermodynamic data for CaO.

Unsteady Free-Convection Flow on a Rotating flat-Plate Under Nonuniform Gravity

The free convection problem with nonuniform gravity finds applications in several fields. For example, centrifugal gravity fieldsarisein many rotating machinery applications. A gravity field is also created artificially in an orbital space station by rotation. The effect of nonuniform gravity due to the rotation of isothermal or nonisothermal plates has been studied by several authors [l-5] using various mathematical techniques.