A computer study of the miscibility limits of Sb-Te-Bi

The quaternary system Sb1bTe1bBi1bSe with small amounts of suitable dopants is of interest for the manufacture of thermoelectric modules which exhibit the Peltier and Seebeck effects. This property could be useful in the production of energy from the thermoelectric effect. Other substances are bismuth telluride (Bi2Te3) and Sb1bTe1bBi and compounds such as ZnIn2Se4. In the present paper the application of computer programs such as MIGAP of Kaufman is used to indicate the stability of the ternary limits of Sb1bTe1bBi within the temperature ranges of interest, namely 273 K to 300 K.

Free vibration of circular plates of arbitrary thickness

Free vibration of circular plates of arbitrary thickness is investigated using the method of initial functions. State-space approach is used to derive the governing equations of the above method. The formulation is such that theories of any desired order can be obtained by deleting higher terms in the infinite-order differential equations. Numerical results are obtained for flexural and extensional vibration of circular plates. Results are also computed using Mindlin's theory and they are in agreement with the present analysis.

Integration-free algorithms for fixed end-point regulator problems

A new formulation is suggested for the fixed end-point regulator problem, which, in conjunction with the recently developed integration-free algorithms, provides an efficient means of obtaining numerical solutions to such problems.

Circulation and salinity distribution in coastal inlets

The mathematical model developed by Hansen and Rattray based on Pritchard's equations for a coastal-plain estuary has been analysed to study the circulation and salinity distributions in coastal inlets with constant width and depth. Numerical solutions of the basic equations have been obtained without placing any restriction on Rayleigh numbers. A noteworthy contribution of the present analysis is that solutions of equations have been obtained for higher Rayleigh numbers, which was not possible in the earlier model. It is found that the effect of higher Rayleigh numbers is to increase the vertical advection, making the salinities in the upper and lower layers more uniform with a distinct halocline near the mid-depths. Solutions are discussed for some special cases of practical interest.

Unsteady mixed convection over two-dimensional and axisymmetric bodies

The unsteady laminar incompressible mixed convection flow over a two-dimensional body (cylinder) and an axisymmetric body (sphere) has been studied when the buboyancy forces arise from both thermal and mass diffusion and the unsteadiness in the flow field is introduced by the time dependent free stream velocity. The nonlinear partial differential equations with three independent variables governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The results indicate that for the thermally assisting flow the local skin friction, heat transfer and mass diffusion are enhanced when the buoyancy force from mass diffusion assists the thermal buoyancy force. But this trend is opposite for the thermally opposing flow. The point of zero skin friction moves upstream due to unsteadiness. No singularity is observed at the point of zero skin friction for unsteady flow unlike steady flow. The flow reversal is observed after a certain instant of time. The velocity overshoot occurs for assisting flows.

Unsteady boundary layer flow on a cylinder in a channel

The unsteady laminar incompressible nonsimilar boundary layer flow over a circular cylinder placed symmetrically inside a channel has been studied when the unsteadiness and nonsimilarity are due to the free stream velocity. The nonlinear partial differential equations with three independent variables have been solved numerically using an implicit finite-difference in combination with the quasilinearization technique. It is found that the channel blockage parameter controls the transfer of heat from the cylinder and delays separation. The skin friction and heat transfer are significantly affected by the free stream velocity distributions.

An accurate numerical integration scheme for finite rotations using rotation vector parametrization

A numerical integration procedure for rotational motion using a rotation vector parametrization is explored from an engineering perspective by using rudimentary vector analysis. The incremental rotation vector, angular velocity and acceleration correspond to different tangent spaces of the rotation manifold at different times and have a non-vectorial character. We rewrite the equation of motion in terms of vectors lying in the same tangent space, facilitating vector space operations consistent with the underlying geometric structure. While any integration algorithm (that works within a vector space setting) may be used, we presently employ a family of explicit Runge-Kutta algorithms to solve this equation. While this work is primarily motivated out of a need for highly accurate numerical solutions of dissipative rotational systems of engineering interest, we also compare the numerical performance of the present scheme with some of the invariant preserving schemes, namely ALGO-C1, STW, LIEMIDEA] and SUBCYC-M. Numerical results show better local accuracy via the present approach vis-a-vis the preserving algorithms. It is also noted that the preserving algorithms do not simultaneously preserve all constants of motion. We incorporate adaptive time-stepping within the present scheme and this in turn enables still higher accuracy and a `near preservation' of constants of motion over significantly longer intervals. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Spectral Invariance of SG Pseudo-Differential Operators on L-P(R-n)

We prove the spectral invariance of SG pseudo-differential operators on L-P(R-n), 1 < p < infinity, by using the equivalence of ellipticity and Fredholmness of SG pseudo-differential operators on L-p(R-n), 1 < p < infinity. A key ingredient in the proof is the spectral invariance of SC pseudo-differential operators on L-2(R-n).

Mixed Convection Induced in a Gas Flowing Past a Hot Horizontal Flat Plate

A boundary layer analysis of mixed convective motion over a hot horizontal flat plate is performed under the conditions of steady flow and low speed. Use of the Howarth-Dorodnytsyn transformation makes it possible to dispense with the usual Boussinesq approximation, and variable gas properties are accounted for via the assumption that dynamic viscosity and thermal conductivity are proportional to the absolute temperature. The formulation presented enables the entire mixed convection regime to be described by a single set of equations. Finite difference solutions when the Prandtl number is 0.72 are obtained over the entire range of the mixed convection parameter ξ from 0 (free convection) to 1 (forced convection) and heating parameter ▵ values from 2 to 12. The effects of both ξ and ▵on the velocity profiles, the temperature profiles, and the variation of skin friction and heat transfer functions are clearly illustrated in tables and graphs.

Analysis of an Aquifer-water table aquitard system

The problem of pumping an aquifer in an aquifer-water table aquitard system is considered, accounting for the elastic properties of both the aquifer and the aquitard, the gravity drainage in the aquitard and treating the water table as an unknown boundary. The coupled partial differential equations are nondimensionalised, yielding three principal parameters governing the problem. The numerical solution of these equations is obtained for a wide range of parameter values. Type curves are generated and their use is illustrated through a field application.

An efficient marching algorithm for waterhammer analysis by the method of characteristics

A new fast and efficient marching algorithm is introduced to solve the basic quasilinear, hyperbolic partial differential equations describing unsteady, flow in conduits by the method of characteristics. The details of the marching method are presented with an illustration of the waterhammer problem in a simple piping system both for friction and frictionless cases. It is shown that for the same accuracy the new marching method requires fewer computational steps, less computer memory and time.

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.

Transient Flow of a Compressible Viscous Fluid Near an Infinite Disk

We have consider ed the transient motion of art electrically conducting viscous compressible fluid which is in contact with an insulated infinite disk. The initial motion is considered to be due to the uniform rotation of the disk in an otherwise stationary fluid or due to the uniform rigid rotation of the fluid over a stationary disk. Different cases of transient motion due to finite impulse imparted either to the disk or to the distant fluid have been investigated. Effects of the imposed axial magnetic field and the disk temperature on the transient flow are included. The nonlinear partial differential equations governing the motion are solved numerically using an implicit finite-difference scheme along with the Newton's linearisation technique.

Risk Minimizing Option Pricing for a Class of Exotic Options in a Markov-Modulated Market

We address risk minimizing option pricing in a regime switching market where the floating interest rate depends on a finite state Markov process. The growth rate and the volatility of the stock also depend on the Markov process. Using the minimal martingale measure, we show that the locally risk minimizing prices for certain exotic options satisfy a system of Black-Scholes partial differential equations with appropriate boundary conditions. We find the corresponding hedging strategies and the residual risk. We develop suitable numerical methods to compute option prices.

Accretion of Chaplygin gas upon black holes: formation of faster outflowing winds

We study the accretion of modified Chaplygin gas upon different types of black holes. Modified Chaplygin gas is one of the best candidates for a combined model of dark matter and dark energy. In addition, from a field theoretical point of view the modified Chaplygin gas model is equivalent to that of a scalar field having a self-interacting potential. We formulate the equations related to both spherical accretion and disc accretion, and respective winds. The corresponding numerical solutions of the flow, particularly of velocity, are presented and analysed. We show that the accretion-wind system of modified Chaplygin gas dramatically alters the wind solutions, producing faster winds, upon changes in physical parameters, while accretion solutions qualitatively remain unaffected. This implies that modified Chaplygin gas is more prone to produce outflow which is the natural consequence of the dark energy into the system.