Simplified Mathematical Model for Effects of Freezing on the Low-Temperature Performance of the Lead-Acid Battery

Discharge periods of lead-acid batteries are significantly reduced at subzero centigrade temperatures. The reduction is more than what can he expected due to decreased rates of various processes caused by a lowering of temperature and occurs despite the fact that active materials are available for discharge. It is proposed that the major cause for this is the freezing of the electrolyte. The concentration of acid decreases during battery discharge with a consequent increase in the freezing temperature. A battery freezes when the discharge temperature falls below the freezing temperature. A mathematical model is developed for conditions where charge-transfer reaction is the rate-limiting step. and Tafel kinetics are applicable. It is argued that freezing begins from the midplanes of electrodes and proceeds toward the reservoir in-between. Ionic conduction stops when one of the electrodes freezes fully and the time taken to reach that point, namely the discharge period, is calculated. The predictions of the model compare well to observations made at low current density (C/5) and at -20 and -40 degrees C. At higher current densities, however, diffusional resistances become important and a more complicated moving boundary problem needs to be solved to predict the discharge periods. (C) 2009 The Electrochemical Society.

Gauge theory of a group of diffeomorphisms. I. General principles

Any (N+M)-parameter Lie group G with an N-parameter subgroup H can be realized as a global group of diffeomorphisms on an M-dimensional base space B, with representations in terms of transformation laws of fields on B belonging to linear representations of H. The gauged generalization of the global diffeomorphisms consists of general diffeomorphisms (or coordinate transformations) on a base space together with a local action of H on the fields. The particular applications of the scheme to space-time symmetries is discussed in terms of Lagrangians, field equations, currents, and source identities. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Representations and properties of para-Bose oscillator operators. I. Energy position and momentum eigenstates

Para-Bose commutation relations are related to the SL(2,R) Lie algebra. The irreducible representation [script D]alpha of the para-Bose system is obtained as the direct sum Dbeta[direct-sum]Dbeta+1/2 of the representations of the SL(2,R) Lie algebra. The position and momentum eigenstates are then obtained in this representation [script D]alpha, using the matrix mechanical method. The orthogonality, completeness, and the overlap of these eigenstates are derived. The momentum eigenstates are also derived using the wave mechanical method by specifying the domain of the definition of the momentum operator in addition to giving it a formal differential expression. By a careful consideration in this manner we find that the two apparently different solutions obtained by Ohnuki and Kamefuchi in this context are actually unitarily equivalent. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Mathematical model of arterial stenosis

A mathematical model for pulsatile flow in a partially occluded tube is presented. The problem has applications in studying the effects of blood flow characteristics on atherosclerotic development. The model brings out the importance of the pulsatility of blood flow on separation and the stress distribution. The results obtained show fairly good agreement with the available experimental results.

Generalized Burgers equations and Euler-Painleve transcendents. II

It was proposed earlier [P. L. Sachdev, K. R. C. Nair, and V. G. Tikekar, J. Math. Phys. 27, 1506 (1986)] that the Euler Painlevé equation yy[script `]+ay[script ']2+ f(x)yy[script ']+g(x) y2+by[script ']+c=0 represents the generalized Burgers equations (GBE's) in the same manner as Painlevé equations do the KdV type. The GBE was treated with a damping term in some detail. In this paper another GBE ut+uaux+Ju/2t =(gd/2)uxx (the nonplanar Burgers equation) is considered. It is found that its self-similar form is again governed by the Euler Painlevé equation. The ranges of the parameter alpha for which solutions of the connection problem to the self-similar equation exist are obtained numerically and confirmed via some integral relations derived from the ODE's. Special exact analytic solutions for the nonplanar Burgers equation are also obtained. These generalize the well-known single hump solutions for the Burgers equation to other geometries J=1,2; the nonlinear convection term, however, is not quadratic in these cases. This study fortifies the conjecture regarding the importance of the Euler Painlevé equation with respect to GBE's. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Resonances, scattering theory, and rigged Hilbert spaces

The problem of decaying states and resonances is examined within the framework of scattering theory in a rigged Hilbert space formalism. The stationary free,''in,'' and ''out'' eigenvectors of formal scattering theory, which have a rigorous setting in rigged Hilbert space, are considered to be analytic functions of the energy eigenvalue. The value of these analytic functions at any point of regularity, real or complex, is an eigenvector with eigenvalue equal to the position of the point. The poles of the eigenvector families give origin to other eigenvectors of the Hamiltonian: the singularities of the ''out'' eigenvector family are the same as those of the continued S matrix, so that resonances are seen as eigenvectors of the Hamiltonian with eigenvalue equal to their location in the complex energy plane. Cauchy theorem then provides for expansions in terms of ''complete'' sets of eigenvectors with complex eigenvalues of the Hamiltonian. Applying such expansions to the survival amplitude of a decaying state, one finds that resonances give discrete contributions with purely exponential time behavior; the background is of course present, but explicitly separated. The resolvent of the Hamiltonian, restricted to the nuclear space appearing in the rigged Hilbert space, can be continued across the absolutely continuous spectrum; the singularities of the continuation are the same as those of the ''out'' eigenvectors. The free, ''in'' and ''out'' eigenvectors with complex eigenvalues and those corresponding to resonances can be approximated by physical vectors in the Hilbert space, as plane waves can. The need for having some further physical information in addition to the specification of the total Hamiltonian is apparent in the proposed framework. The formalism is applied to the Lee–Friedrichs model and to the scattering of a spinless particle by a local central potential. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Doping Dependence of Semiconductor-Metal Transition in InP at High Pressures

The resistivity of selenium-doped n-InP single crystal layers grown by liquid-phase epitaxy with electron concentrations varying from 6.7 x 10$^18$ to 1.8 x 10$^20$ cm$^{-3}$ has been measured as a function of hydrostatic pressure up to 10 GPa. Semiconductor-metal transitions were observed in each case with a change in resistivity by two to three orders of magnitude. The transition pressure p$_c$ decreased monotonically from 7.24 to 5.90 GPa with increasing doping concentration n according to the relation $p_c = p_o [1 - k(n/n_m)^a]$, where n$_m$ is the concentration (per cubic centimetre) of phosphorus donor sites in InP atoms, p$_o$ is the transition pressure at low doping concentrations, k is a constant and $\alpha$ is an exponent found experimentally to be 0.637. The decrease in p$_c$ is considered to be due to increasing internal stress developed at high concentrations of ionized donors. The high-pressure metallic phase had a resistivity (2.02-6.47) x 10$^{-7}$ $\Omega$ cm, with a positive temperature coefficient dependent on doping.

Generalized Burgers equations and Euler-Painlev transcendents. I

Initial-value problems for the generalized Burgers equation (GBE) ut+u betaux+lambdaualpha =(delta/2)uxx are discussed for the single hump type of initial data both continuous and discontinuous. The numerical solution is carried to the self-similar ``intermediate asymptotic'' regime when the solution is given analytically by the self-similar form. The nonlinear (transformed) ordinary differential equations (ODE's) describing the self-similar form are generalizations of a class discussed by Euler and Painlevé and quoted by Kamke. These ODE's are new, and it is postulated that they characterize GBE's in the same manner as the Painlev equations categorize the Kortweg-de Vries (KdV) type. A connection problem for some related ODE's satisfying proper asymptotic conditions at x=±[infinity], is solved. The range of amplitude parameter is found for which the solution of the connection problem exists. The other solutions of the above GBE, which display several interesting features such as peaking, breaking, and a long shelf on the left for negative values of the damping coefficient lambda, are also discussed. The results are compared with those holding for the modified KdV equation with damping. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Singularity-free cosmology: A simple model

The nonminimal coupling of a self-interacting complex scalar field with gravity is studied. For a Robertson-Walker open universe the stable solutions of the scalar-field equations are time dependent. As a result of this, a novel spontaneous symmetry breaking occurs which leads to a varying effective gravitational coupling coefficient. It is found that the coupling coefficient changes sign below a critical ‘‘radius’’ of the Universe implying the appearance of repulsive gravity. The occurrence of the repulsive interaction at an early epoch facilitates singularity avoidance. The model also provides a solution to the horizon problem.

Application of Mathematical Programming to the Plywood Design and Manufacturing Problem

Plywood manufacture includes two fundamental stages. The first is to peel or separate logs into veneer sheets of different thicknesses. The second is to assemble veneer sheets into finished plywood products. At the first stage a decision must be made as to the number of different veneer thicknesses to be peeled and what these thicknesses should be. At the second stage, choices must be made as to how these veneers will be assembled into final products to meet certain constraints while minimizing wood loss. These decisions present a fundamental management dilemma. Costs of peeling, drying, storage, handling, etc. can be reduced by decreasing the number of veneer thicknesses peeled. However, a reduced set of thickness options may make it infeasible to produce the variety of products demanded by the market or increase wood loss by requiring less efficient selection of thicknesses for assembly. In this paper the joint problem of veneer choice and plywood construction is formulated as a nonlinear integer programming problem. A relatively simple optimal solution procedure is developed that exploits special problem structure. This procedure is examined on data from a British Columbia plywood mill. Restricted to the existing set of veneer thicknesses and plywood designs used by that mill, the procedure generated a solution that reduced wood loss by 79 percent, thereby increasing net revenue by 6.86 percent. Additional experiments were performed that examined the consequences of changing the number of veneer thicknesses used. Extensions are discussed that permit the consideration of more than one wood species.

Representation and properties of para-Bose oscillator operators. II. Coherent states and the minimum uncertainty states

The energy, position, and momentum eigenstates of a para-Bose oscillator system were considered in paper I. Here we consider the Bargmann or the analytic function description of the para-Bose system. This brings in, in a natural way, the coherent states ||z;alpha> defined as the eigenstates of the annihilation operator ?. The transformation functions relating this description to the energy, position, and momentum eigenstates are explicitly obtained. Possible resolution of the identity operator using coherent states is examined. A particular resolution contains two integrals, one containing the diagonal basis ||z;alpha><−z;alpha||. We briefly consider the normal and antinormal ordering of the operators and their diagonal and discrete diagonal coherent state approximations. The problem of constructing states with a minimum value of the product of the position and momentum uncertainties and the possible alpha dependence of this minimum value is considered. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Continuum percolation with discs having a distribution of radii

It is shown that for continuum percolation with overlapping discs having a distribution of radii, the net areal density of discs at percolation threshold depends non-trivially on the distribution, and is not bounded by any finite constant. Results of a Monte Carlo simulation supporting the argument are presented.

Self-Induced Transparency in Excitons

The interaction of intense coherent light with Frenkel excitons has been studied for investigating the self-induced transparency. Some nonlinear effects neglected before have been included. It is found that the frequency spectrum consistent with the pulse propagation is wider by two orders of magnitude compared with the previous result.

A phase space approach to generalized Hamilton–Jacobi theory

Generalizations of H–J theory have been discussed before in the literature. The present approach differs from others in that it employs geometrical ideas on phase space and classical transformation theory to derive the basic equations. The relation between constants of motion and symmetries of the generalized H–J equations is then clarified. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Cosmological Constant and Scalar Gravitons

If a cosmological term is included in the equations of general relativity, the linearized equations can be interpreted as a tensor-scalar theory of finite-range gravitation. The scalar field cannot be transformed away be a gauge transformation (general co-ordinate transformation) and so must be interpreted as a physically significant degree of freedom. The hypothesis that a massive spin-two meson (mass m2) satisfied equations identical in form to the equations of general relativity leads to the prediction of a massive spin-zero meson (mass m0), the ratio of masses being m0 / m2 = 3*3.