Thermodynamic consistency of the ternary functions

A ternary thermodynamic function has been developed based on statistico-thermodynamic considerations, with a particular emphasis on the higher-order terms indicating the effects of truncation at the various stages of the treatment. Although the truncation of a series involved in the equation introduces inconsistency, the latter may be removed by imposing various thermodynamic boundary conditions. These conditions are discussed in the paper. The present equation with higher-order terms shows that the α function of a component reduces to a quadratic function of composition at constant compositional paths involving the other two components in the system. The form of the function has been found to be representative of various experimental observations.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Elastic property estimation using Ultrasound Assisted Optical Elastography through Remote Palpation - A simulation study - art. no. 61640Q

We propose an effective elastography technique in which an acoustic radiation force is used for remote palpation to generate localized tissue displacements, which are directly correlated to localized variations of tissue stiffness and are measured using a light probe in the same direction of ultrasound propagation. The experimental geometry has provision to input light beam along the ultrasound propagation direction, and hence it can be prealigned to ensure proper interception of the focal region by the light beam. Tissue-mimicking phantoms with homogeneous and isotropic mechanical properties of normal and malignant breast tissue are considered for the study. Each phantom is insonified by a focusing ultrasound transducer (1 MHz). The focal volume of the transducer and the ultrasound radiation force in the region are estimated through solving acoustic wave propagation through medium assuming average acoustic properties. The forward elastography problem is solved for the region of insonification assuming the Lame's parameters and Poisson's ratio, under Dirichlet boundary conditions which gives a distribution of displacement vectors. The direction of displacement, though presented spatial variation, is predominantly towards the ultrasound propagation direction. Using Monte Carlo (MC) simulation we have traced the photons through the phantom and collected the photons arriving at the detector on the boundary of the object in the direction of ultrasound. The intensity correlations are then computed from detected photons. The intensity correlation function computed through MC simulation showed a modulation whose strength is found to be proportional to the amplitude of displacement and inversely related to the storage (elastic) modulus. It is observed that when the storage modulus in the focal region is increased the computed displacement magnitude, as indicated by the depth of modulation in the intensity autocorrelation, decreased and the trend is approximately exponential.

Flow Generated By Pitched Blade Turbines Ii: Simulation Using Κ-Ε Model

Experimental data on average velocity and turbulence intensity generated by pitched blade downflow turbines (PTD) were presented in Part I of this paper. Part II presents the results of the simulation of flow generated by PTD The standard κ-ε model along with the boundary conditions developed in the Part 1 have been employed to predict the flow generated by PTD in cylindrical baffled vessel. This part describes the new software FIAT (Flow In Agitated Tanks) for the prediction of three dimensional flow in stirred tanks. The basis of this software has been described adequately. The influence of grid size, impeller boundary conditions and values of model parameters on the predicted flow have been analysed. The model predictions successfully reproduce the three dimensionality and the other essential characteristics of the flow. The model can be used to improve the overall understanding about the relative distribution of turbulence by PTD in the agitated tank

Circular elastic inclusions in cylindrical shells

The problem of a circular elastic inclusion in a cylindrical shell subjected to internal pressure or thermal loading is studied. The two shallow-shell equations governing the behaviour of a cylindrical shell are transformed into a single differential equation involving a curvature parameter and a complex potential function in a non-dimensional form. In the shell region, the solution is represented by Hankel functions of first kind, whereas in the inclusion region it is represented by Bessel functions of first kind. Boundary conditions at the shell-inclusion junction are expressed in a simple form involving in-plane strains and change in curvature. The effect of such inclusion parameters as extensional rigidity, bending rigidity, and thermal expansion coefficients on the stress concentrations has been determined. The results are presented in non-dimensional form for ready use.

Theoretical Estimation of Length Dependent In-Plane Stiffness of Single Walled Carbon Nanotubes Using the Nonlocal Elasticity Theory

In the present paper, Eringen's nonlocal elasticity theory is employed to evaluate the length dependent in-plane stiffness of single-walled carbon nanotubes (SWCNTs). The SWCNT is modeled as an Euler-Bernoulli beam and is analyzed for various boundary conditions to evaluate the length dependent in-plane stiffness. It has been found that the nonlocal scaling parameter has a significant effect on the length dependent in-plane stiffness of SWCNTs. It has been observed that as the nonlocal scale parameter increases the stiffness ratio of SWCNT decreases. In nonlocality, the cantilever SWCNT has high in-plane stiffness as compared to the simply-supported and the clamped cases.

Thermal stresses in a spherical shell with a conical nozzle

A thermal stress problem of a spherical shell with a conical nozzle is solved using a continuum approach. The thermal loading consists of a steady temperature which is uniform on the inner and outer surfaces of the shell and the conical nozzle but may vary linearly across the thickness. The thermal stress problem is converted to an equivalent boundary value problem and boundary conditions are specified at the junction of the spherical shell and conical nozzle. The stresses are obtained for a uniform increase in temperature and for a linear variation of temperature across the thickness of the shell, and are presented in graphical form for ready use.

Svekofennialaisen orogenian ylityönnön terminen mallinnus Itä-Suomen arkeeisella kratonilla

This study reviews the thermophysical properties of rocks and the lithosphere and describes a one dimensional thermal numeric model of hypothetical 20 km thick overthrust plate obtruded on to the Archeaen craton in the Svecofennian orogeny (1.92-1.77 Ga). The objective is to find out if the overthrust plate and its radiogenic heat sources were able to produce the thermal effects observed on the current erosion level of the Archaean craton. Heat transfer in lithosphere is assumed conductive, and advective heat transfer due to melting and melt transfer is supposed negligible. The study area is located in the Eastern Finland, approximately current Kainuu and Northern Karelia regions, east from the most active orogenic belt (Raahe-Ladoga zone), so that orogenic magmatism can be neglected. Physical parameters and boundary conditions for the model are from different earlier published sources: deep seismic profiles (rock variation in depth), laboratory measurements (heat production and conductivity of rocks), field measurements (heat flow densities), and pT(t) estimations from the Finnish precambrian to estimate the size and thickness of the sheet. Comparison of the modelling results to previous K-Ar datings and other pTt estimations show, that the effect of the overthrust sheet has been adequate (max. T 450°C at 4 kbar) to produce the K-Ar resetting ages measured from the Archaean bedrock at current erosion level. No other kind of thermal activation in lithosphere is required. Results show possibly very minor partial melting in upper middle crust underneath the overthrust sheet.

Modelling of symmetric laminates under extension

A higher-order theory of laminated composites under in-plane loads is developed. The displacement field is expanded in terms of the thickness co-ordinate, satisfying the zero shear stress condition at the surfaces of the laminate. Using the principle of virtual displacement, the governing equations and boundary conditions are established. Numerical results for interlaminar stresses arising in the case of symmetric laminates under uniform extension have been obtained and are compared with similar results available in the literature.

Analysis of performance of a hot gas injection thrust vector control system

The complex three-dimensional flowfield produced by secondary injection of hot gases in a rocket nozzle for thrust vector control is analyzed by solving unsteady three-dimensional Euler equations with appropriate boundary conditions. Various system performance parameters like secondary jet amplification factor and axial thrust augmentation are deduced by integrating the nozzle wall pressure distributions obtained as part of the flowfield solution and compared with measurements taken in actual static tests. The agreement is good within the practical range of secondary injectant flow rates for thrust vector control applications.

Analysis of transients in a canal network

A generalised formulation of the mathematical model developed for the analysis of transients in a canal network, under subcritical flow, with any realistic combination of control structures and their multiple operations, has been presented. The model accounts for a large variety of control structures such as weirs, gates, notches etc. discharging under different conditions, namely submerged and unsubmerged. A numerical scheme to compute and approximate steady state flow condition as the initial condition has also been presented. The model can handle complex situations that may arise from multiple gate operations. This has been demonstrated with a problem wherein the boundary conditions change from a gate discharge equation to an energy equation and back to a gate discharge equation. In such a situation the wave strikes a fixed gate and leads to large and rapid fluctuations in both discharge and depth.

Overwintering ecology of northern field layer plants : snow and photosynthesis in Vaccinium vitis-idaea L.

Winter is a significant period for the seasonality of northern plants, but is often overlooked when studying the interactions of plants and their environment. This study focuses on the effects of overwintering conditions, including warm winter periods, snow, and snowmelt on boreal and sub-Arctic field layer plants. Wintertime photosynthesis and related physiological factors of evergreen dwarf shrubs, particularly of Vaccinium vitis-idaea, are emphasised. The work combines experiments both in the field and in growth chambers with measurements in natural field conditions. Evergreen dwarf shrubs are predominantly covered by snow in the winter. The protective snow cover provides favourable conditions for photosynthesis, especially during the spring before snowmelt. The results of this study indicate that photosynthesis occurs under the snow in V. vitis-idaea. The light response of photosynthesis determined in field conditions during the period of snow cover shows that positive net CO2 exchange is possible under the snow in the prevailing light and temperature. Photosynthetic capacity increases readily during warm periods in winter and the plants are thus able to replenish carbohydrate reserves lost through respiration. Exposure to low temperatures in combination with high light following early snowmelt can set back photosynthesis as sustained photoprotective measures are activated and photodamage begins to build up. Freezing may further decrease the photosynthetic capacity. The small-scale distribution of many field layer plants, including V. vitis-idaea and other dwarf shrubs, correlates with the snow distribution in a forest. The results of this study indicate that there are species-specific differences in the snow depth affinity of the field and ground layer species. Events and processes taking place in winter can have a profound effect on the overall performance of plants and on the interactions between plants and their environment. Understanding the processes involved in the overwintering of plants is increasingly important as the wintertime climate in the north is predicted to change in the future.

The 26 July 2005 heavy rainfall event over Mumbai: numerical modeling aspects

The performance of the Advanced Regional Prediction System (ARPS) in simulating an extreme rainfall event is evaluated, and subsequently the physical mechanisms leading to its initiation and sustenance are explored. As a case study, the heavy precipitation event that led to 65 cm of rainfall accumulation in a span of around 6 h (1430 LT-2030 LT) over Santacruz (Mumbai, India), on 26 July, 2005, is selected. Three sets of numerical experiments have been conducted. The first set of experiments (EXP1) consisted of a four-member ensemble, and was carried out in an idealized mode with a model grid spacing of 1 km. In spite of the idealized framework, signatures of heavy rainfall were seen in two of the ensemble members. The second set (EXP2) consisted of a five-member ensemble, with a four-level one-way nested integration and grid spacing of 54, 18, 6 and 1 km. The model was able to simulate a realistic spatial structure with the 54, 18, and 6 km grids; however, with the 1 km grid, the simulations were dominated by the prescribed boundary conditions. The third and final set of experiments (EXP3) consisted of a five-member ensemble, with a four-level one-way nesting and grid spacing of 54, 18, 6, and 2 km. The Scaled Lagged Average Forecasting (SLAF) methodology was employed to construct the ensemble members. The model simulations in this case were closer to observations, as compared to EXP2. Specifically, among all experiments, the timing of maximum rainfall, the abrupt increase in rainfall intensities, which was a major feature of this event, and the rainfall intensities simulated in EXP3 (at 6 km resolution) were closest to observations. Analysis of the physical mechanisms causing the initiation and sustenance of the event reveals some interesting aspects. Deep convection was found to be initiated by mid-tropospheric convergence that extended to lower levels during the later stage. In addition, there was a high negative vertical gradient of equivalent potential temperature suggesting strong atmospheric instability prior to and during the occurrence of the event. Finally, the presence of a conducive vertical wind shear in the lower and mid-troposphere is thought to be one of the major factors influencing the longevity of the event.

Diffraction of a cylindrical pulse by a metallic half plane

A direct transform technique is applied to the initial and boundary value problem involving diffraction of a cylindrical pulse by a half plane, on which impedance type of boundary conditions must be met by the total field. The solution to the time harmonic incident plane wave is deduced as a particular case of the general time-dependent problem considered here and we avoid the Wiener–Hopf technique which leads to very complicated factorization and which masks the role of the impedance factor Z′ (a small quantity) in the expression for the scattered field.

On the possibility of an alpha-sq omega-type dynamo in a thin layer inside the sun

If the solar dynamo operates in a thin layer of 10,000-km thickness at the interface between the convection zone and the radiative core, using the facts that the dynamo should have a period of 22 years and a half-wavelength of 40 deg in the theta-direction, it is possible to impose restrictions on the values which various dynamo parameters are allowed to have. It is pointed out that the dynamo should be of alpha-sq omega nature, and kinematical calculations are presented for free dynamo waves and for dynamos in thin rectangular slabs with appropriate boundary conditions. An alpha-sq omega dynamo is expected to produce a significant poloidal field which does not leak to the solar surface. It is found that the turbulent diffusity eta and alpha-coefficient are restricted to values within about a factor of 10, the median values being eta of about 10 to the 10th sq cm/sec and alpha of about 10 cm/sec. On the basis of mixing length theory, it is pointed out that such values imply a reasonable turbulent velocity of the order 30 m/s, but rather small turbulent length scales like 300 km.