Analysis of the free vibration of rectangular plates with central cut-outs using the discrete Ritz method

We present an analysis of the free vibration of plates with internal discontinuities due to central cut-outs. A numerical formulation for a basic L-shaped element which is divided into appropriate sub-domains that are dependent upon the location of the cut-out is used as the basic building element. Trial functions formed to satisfy certain boundary conditions are employed to define the transverse deflection of each sub-domain. Mathematical treatments in terms of the continuities in displacement, slope, moment, and higher derivatives between the adjacent sub-domains are enforced at the interconnecting edges. The energy functional results, from the proper assembly of the coupled strain and kinetic energy contributions of each sub-domain, are minimized via the Ritz procedure to extract the vibration frequencies and. mode shapes of the plates. The procedures are demonstrated by considering plates with central cut-outs that are subjected to two types of boundary conditions. (C) 2003 Elsevier Ltd. All rights reserved.

Uniqueness for boundary value problems for second order finite difference equations

We establish maximum principles for second order difference equations and apply them to obtain uniqueness for solutions of some boundary value problems.

Solution of the dual reflection equation for A(n-1)((1)) solid-on-solid model

We obtain a diagonal solution of the dual reflection equation for the elliptic A(n-1)((1)) solid-on-solid model. The isomorphism between the solutions of the reflection equation and its dual is studied. (C) 2004 American Institute of Physics.

Evaluation of the BEASY program using linear and piecewise linear approaches for the boundary conditions

The boundary element method (BEM) was used to study galvanic corrosion using linear and logarithmic boundary conditions. The linear boundary condition was implemented by using the linear approach and the piecewise linear approach. The logarithmic boundary condition was implemented by the piecewise linear approach. The calculated potential and current density distribution were compared with the prior analytical results. For the linear boundary condition, the BEASY program using the linear approach and the piecewise linear approach gave accurate predictions of the potential and the galvanic current density distributions for varied electrolyte conditions, various film thicknesses, various electrolyte conductivities and various area ratio of anode/cathode. The 50-point piecewise linear method could be used with both linear and logarithmic polarization curves.

Recent developments in the quantum dynamical characterization of unimolecular resonances

We give a selective review of quantum mechanical methods for calculating and characterizing resonances in small molecular systems, with an emphasis on recent progress in Chebyshev and Lanczos iterative methods. Two archetypal molecular systems are discussed: isolated resonances in HCO, which exhibit regular mode and state specificity, and overlapping resonances in strongly bound HO2, which exhibit irregular and chaotic behavior. Recent progresses for non-zero total angular momentum J calculations of resonances including parallel computing models are also included and future directions in this field are discussed.

Gauge Poisson representations for birth/death master equations

Poisson representation techniques provide a powerful method for mapping master equations for birth/death processes -- found in many fields of physics, chemistry and biology -- into more tractable stochastic differential equations. However, the usual expansion is not exact in the presence of boundary terms, which commonly occur when the differential equations are nonlinear. In this paper, a gauge Poisson technique is introduced that eliminates boundary terms, to give an exact representation as a weighted rate equation with stochastic terms. These methods provide novel techniques for calculating and understanding the effects of number correlations in systems that have a master equation description. As examples, correlations induced by strong mutations in genetics, and the astrophysical problem of molecule formation on microscopic grain surfaces are analyzed. Exact analytic results are obtained that can be compared with numerical simulations, demonstrating that stochastic gauge techniques can give exact results where standard Poisson expansions are not able to.

GCMC-surface area of carbonaceous materials with N-2 and Ar adsorption as an alternative to the classical BET method

In this paper we apply a new method for the determination of surface area of carbonaceous materials, using the local surface excess isotherms obtained from the Grand Canonical Monte Carlo simulation and a concept of area distribution in terms of energy well-depth of solid–fluid interaction. The range of this well-depth considered in our GCMC simulation is from 10 to 100 K, which is wide enough to cover all carbon surfaces that we dealt with (for comparison, the well-depth for perfect graphite surface is about 58 K). Having the set of local surface excess isotherms and the differential area distribution, the overall adsorption isotherm can be obtained in an integral form. Thus, given the experimental data of nitrogen or argon adsorption on a carbon material, the differential area distribution can be obtained from the inversion process, using the regularization method. The total surface area is then obtained as the area of this distribution. We test this approach with a number of data in the literature, and compare our GCMC-surface area with that obtained from the classical BET method. In general, we find that the difference between these two surface areas is about 10%, indicating the need to reliably determine the surface area with a very consistent method. We, therefore, suggest the approach of this paper as an alternative to the BET method because of the long-recognized unrealistic assumptions used in the BET theory. Beside the surface area obtained by this method, it also provides information about the differential area distribution versus the well-depth. This information could be used as a microscopic finger-print of the carbon surface. It is expected that samples prepared from different precursors and different activation conditions will have distinct finger-prints. We illustrate this with Cabot BP120, 280 and 460 samples, and the differential area distributions obtained from the adsorption of argon at 77 K and nitrogen also at 77 K have exactly the same patterns, suggesting the characteristics of this carbon.

Adsorption of argon on homogeneous graphitized thermal carbon black and heterogeneous carbon surface

In this paper we investigate the effects of surface mediation on the adsorption behavior of argon at different temperatures on homogeneous graphitized thermal carbon black and on heterogeneous nongraphitized carbon black surface. The grand canonical Monte Carlo (GCMC) simulation is used to study the adsorption, and its performance is tested against a number of experimental data on graphitized thermal carbon black (which is known to be highly homogeneous) that are available in the literature. The surface-mediation effect is shown to be essential in the correct description of the adsorption isotherm because without accounting for that effect the GCMC simulation results are always greater than the experimental data in the region where the monolayer is being completed. This is due to the overestimation of the fluid–fluid interaction between particles in the first layer close to the solid surface. It is the surface mediation that reduces this fluid–fluid interaction in the adsorbed layers, and therefore the GCMC simulation results accounting for this surface mediation that are presented in this paper result in a better description of the data. This surface mediation having been determined, the surface excess of argon on heterogeneous carbon surfaces having solid–fluid interaction energies different from the graphite can be readily obtained. Since the real heterogeneous carbon surface is not the same as the homogeneous graphite surface, it can be described by an area distribution in terms of the well depth of the solid–fluid energy. Assuming a patchwise topology of the surface with patches of uniform well depth of solid–fluid interaction, the adsorption on a real carbon surface can be determined as an integral of the local surface excess of each patch with respect to the differential area. When this is matched against the experimental data of a carbon surface, we can derive the area distribution versus energy and hence the geometrical surface area. This new approach will be illustrated with the adsorption of argon on a nongraphitized carbon at 87.3 and 77 K, and it is found that the GCMC surface area is different from the BET surface area by about 7%. Furthermore, the description of the isotherm in the region of BET validity of 0.06 to 0.2 is much better with our method than with the BET equation.

Solvability of singular second order m-point boundary value problems of Dirichlet type

Let f : [0, 1] x R2 -> R be a function satisfying the Caxatheodory conditions and t(1 - t)e(t) epsilon L-1 (0, 1). Let a(i) epsilon R and xi(i) (0, 1) for i = 1,..., m - 2 where 0 < xi(1) < xi(2) < (...) < xi(m-2) < 1 - In this paper we study the existence of C[0, 1] solutions for the m-point boundary value problem [GRAPHICS] The proof of our main result is based on the Leray-Schauder continuation theorem.

Multiplicity results for second-order two-point boundary value problems with nonlinearities across several eigenvalues

We consider the boundary value problems for nonlinear second-order differential equations of the form u '' + a(t)f (u) = 0, 0 < t < 1, u(0) = u (1) = 0. We give conditions on the ratio f (s)/s at infinity and zero that guarantee the existence of solutions with prescribed nodal properties. Then we establish existence and multiplicity results for nodal solutions to the problem. The proofs of our main results are based upon bifurcation techniques. (c) 2004 Elsevier Ltd. All rights reserved.

Multiplicity results for second-order two-point boundary value problems with superlinear or sublinear nonlinearities

We consider boundary value problems for nonlinear second order differential equations of the form u + a(t) f(u) = 0, t epsilon (0, 1), u(0) = u(1) = 0, where a epsilon C([0, 1], (0, infinity)) and f : R --> R is continuous and satisfies f (s)s > 0 for s not equal 0. We establish existence and multiplicity results for nodal solutions to the problems if either f(0) = 0, f(infinity) = infinity or f(0) = infinity, f(0) = 0, where f (s)/s approaches f(0) and f(infinity) as s approaches 0 and infinity, respectively. We use bifurcation techniques to prove our main results. (C) 2004 Elsevier Inc. All rights reserved.

On the Schrodinger equation involving a critical Sobolev exponent and magnetic field

We consider the semilinear Schrodinger equation -Delta(A)u + V(x)u = Q(x)vertical bar u vertical bar(2* -2) u. Assuming that V changes sign, we establish the existence of a solution u not equal 0 in the Sobolev space H-A,V(1) + (R-N). The solution is obtained by a min-max type argument based on a topological linking. We also establish certain regularity properties of solutions for a rather general class of equations involving the operator -Delta(A).

Projected Gross-Pitaevskii equation for harmonically confined Bose gases at finite temperature

We extend the projected Gross-Pitaevskii equation formalism of Davis [Phys. Rev. Lett. 87, 160402 (2001)] to the experimentally relevant case of thermal Bose gases in harmonic potentials and outline a robust and accurate numerical scheme that can efficiently simulate this system. We apply this method to investigate the equilibrium properties of the harmonically trapped three-dimensional projected Gross-Pitaevskii equation at finite temperature and consider the dependence of condensate fraction, position, and momentum distributions and density fluctuations on temperature. We apply the scheme to simulate an evaporative cooling process in which the preferential removal of high-energy particles leads to the growth of a Bose-Einstein condensate. We show that a condensate fraction can be inferred during the dynamics even in this nonequilibrium situation.

Properties of the stochastic Gross-Pitaevskii equation: finite temperature Ehrenfest relations and the optimal plane wave representation

We present Ehrenfest relations for the high temperature stochastic Gross-Pitaevskii equation description of a trapped Bose gas, including the effect of growth noise and the energy cutoff. A condition for neglecting the cutoff terms in the Ehrenfest relations is found which is more stringent than the usual validity condition of the truncated Wigner or classical field method-that all modes are highly occupied. The condition requires a small overlap of the nonlinear interaction term with the lowest energy single particle state of the noncondensate band, and gives a means to constrain dynamical artefacts arising from the energy cutoff in numerical simulations. We apply the formalism to two simple test problems: (i) simulation of the Kohn mode oscillation for a trapped Bose gas at zero temperature, and (ii) computing the equilibrium properties of a finite temperature Bose gas within the classical field method. The examples indicate ways to control the effects of the cutoff, and that there is an optimal choice of plane wave basis for a given cutoff energy. This basis gives the best reproduction of the single particle spectrum, the condensate fraction and the position and momentum densities.