Boundary value problems for systems of difference equations associated with systems of second-order ordinary differential equations

We study difference equations which arise as discrete approximations to two-point boundary value problems for systems of second-order ordinary differential equations. We formulate conditions which guarantee a priori bounds on first differences of solutions to the discretized problem. We establish existence results for solutions to the discretized boundary value problems subject to nonlinear boundary conditions. We apply our results to show that solutions to the discrete problem converge to solutions of the continuous problem in an aggregate sense. (C) 2002 Elsevier Science Ltd. All rights reserved.

Difference equations associated with fully nonlinear boundary value problems for second order ordinary differential equations

We study the continuous problem y"=f(x,y,y'), xc[0,1], 0=G((y(0),y(1)),(y'(0), y'(1))), and its discrete approximation (y(k+1)-2y(k)+y(k-1))/h(2) =f(t(k), y(k), v(k)), k = 1,..., n-1, 0 = G((y(0), y(n)), (v(1), v(n))), where f and G = (g(0), g(1)) are continuous and fully nonlinear, h = 1/n, v(k) = (y(k) - y(k-1))/h, for k =1,..., n, and t(k) = kh, for k = 0,...,n. We assume there exist strict lower and strict upper solutions and impose additional conditions on f and G which are known to yield a priori bounds on, and to guarantee the existence of solutions of the continuous problem. We show that the discrete approximation also has solutions which approximate solutions of the continuous problem and converge to the solution of the continuous problem when it is unique, as the grid size goes to 0. Homotopy methods can be used to compute the solution of the discrete approximation. Our results were motivated by those of Gaines.

Prediction of outcomes in hypertensive patients with suspected coronary disease

Stress echocardiography has been shown to improve the diagnosis of coronary artery disease in the presence of hypertension, but its value in prognostic evaluation is unclear. We sought to determine whether stress echocardiography could be used to predict mortality in 2363 patients with hypertension, who were followed for up to 10 years (mean 4.0+/-1.8) for death and revascularization. Stress echocardiograms were normal in 1483 patients (63%), 16% had resting left ventricular (LV) dysfunction alone, and 21% had ischemia. Abnormalities were confined to one territory in 489 patients (21%) and to multiple territories in 365 patients (15%). Cardiac death was less frequent among the patients able to exercise than among those undergoing dobutamine echocardiography (4% versus 7%, P<0.001). The risk of death in patients with a negative stress echocardiogram was <1% per year. Ischemia identified by stress echocardiography was an independent predictor of mortality in those able to exercise (hazard ratio 2.21, 95% confidence intervals 1.10 to 4.43, P=0.0001) as well as those undergoing dobutamine echo (hazard ratio 2.39, 95% confidence intervals 1.53 to 3.75, P=0.0001); other predictors were age, heart failure, resting LV dysfunction, and the Duke treadmill score. In stepwise models replicating the sequence of clinical evaluation, the results of stress echocardiography added prognostic power to models based on clinical and stress-testing variables. Thus, the results of stress echocardiography are an independent predictor of cardiac death in hypertensive patients with known or suspected coronary artery disease, incremental to clinical risks and exercise results.

Supersymmetric t-J Gaudin models and KZ equations

Supersymmetric t-J Gaudin models with open boundary conditions are investigated by means of the algebraic Bethe ansatz method. Off-shell Bethe ansatz equations of the boundary Gaudin systems are derived, and used to construct and solve the KZ equations associated with sl (2\1)((1)) superalgebra.

Prediction of absolute rate coefficients and product branching ratios for the C(P-3) plus allene reaction system

Complex chemical reactions in the gas phase can be decomposed into a network of elementary (e.g., unimolecular and bimolecular) steps which may involve multiple reactant channels, multiple intermediates, and multiple products. The modeling of such reactions involves describing the molecular species and their transformation by reaction at a detailed level. Here we focus on a detailed modeling of the C(P-3)+allene (C3H4) reaction, for which molecular beam experiments and theoretical calculations have previously been performed. In our previous calculations, product branching ratios for a nonrotating isomerizing unimolecular system were predicted. We extend the previous calculations to predict absolute unimolecular rate coefficients and branching ratios using microcanonical variational transition state theory (mu-VTST) with full energy and angular momentum resolution. Our calculation of the initial capture rate is facilitated by systematic ab initio potential energy surface calculations that describe the interaction potential between carbon and allene as a function of the angle of attack. Furthermore, the chemical kinetic scheme is enhanced to explicitly treat the entrance channels in terms of a predicted overall input flux and also to allow for the possibility of redissociation via the entrance channels. Thus, the computation of total bimolecular reaction rates and partial capture rates is now possible. (C) 2002 American Institute of Physics.

Enriched behavioral prediction equation and its impact on structured learning and the dynamic calculus

This theoretical note describes an expansion of the behavioral prediction equation, in line with the greater complexity encountered in models of structured learning theory (R. B. Cattell, 1996a). This presents learning theory with a vector substitute for the simpler scalar quantities by which traditional Pavlovian-Skinnerian models have hitherto been represented. Structured learning can be demonstrated by vector changes across a range of intrapersonal psychological variables (ability, personality, motivation, and state constructs). Its use with motivational dynamic trait measures (R. B. Cattell, 1985) should reveal new theoretical possibilities for scientifically monitoring change processes (dynamic calculus model; R. B. Cattell, 1996b), such as encountered within psycho therapeutic settings (R. B. Cattell, 1987). The enhanced behavioral prediction equation suggests that static conceptualizations of personality structure such as the Big Five model are less than optimal.

Prediction of coal ash fusion temperatures with the F*A*C*T thermodynamic computer package

This paper outlines research on the processes taking place within the coal mineral matter at high temperatures and development of the relationship between ash fusion temperatures (AFT) and phase equilibria of the coal ash slags. A new thermodynamic database for the Al-Ca-Fe-O-Si system developed by the author was used in conjunction with the thermodynamic computer package F*A*C*T for these purposes. In addition, high temperature experimental studies were undertaken that involved heat treatment and quenching of the ash cones followed by the analyses using different techniques. The study provided new information on the processes taking place during AFT test and demonstrated the validity of the AFTs predictions with F*A*C*T. Examples of practical applications of the AFT prediction method are given in the paper. The results of this study are important not only for the AFT predictions, but also in general for the application of phase equilibrium science to the characterisation of the coal mineral matter interactions at high temperature. (C) 2002 Elsevier Science Ltd. All rights reserved.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

A variable stepsize implementation for stochastic differential equations

Stochastic differential equations (SDEs) arise from physical systems where the parameters describing the system can only be estimated or are subject to noise. Much work has been done recently on developing higher order Runge-Kutta methods for solving SDEs numerically. Fixed stepsize implementations of numerical methods have limitations when, for example, the SDE being solved is stiff as this forces the stepsize to be very small. This paper presents a completely general variable stepsize implementation of an embedded Runge Kutta pair for solving SDEs numerically; in this implementation, there is no restriction on the value used for the stepsize, and it is demonstrated that the integration remains on the correct Brownian path.

Rumors and stable-cause attribution in prediction and behavior

Two stock-market simulation experiments investigated the notion that rumors that invoke stable-cause attributions spawn illusory associations and less regressive predictions and behavior. In Study 1, illusory perceptions of association and stable causation (rumors caused price changes on the day after they appeared) existed despite rigorous conditions of nonassociation (price changes were unrelated to rumors). Predictions (recent price trends will continue) and trading behavior (departures from a strong buy-low-sell-high strategy) were both anti-regressive. In Study 2, stability of attribution was manipulated via a computerized tutorial. Participants taught to view price-changes as caused by stable forces predicted less regressively and departed more from buy-low-sell-high trading patterns than those taught to perceive changes as caused by unstable forces. Results inform a social cognitive and decision theoretic understanding of rumor by integrating it with causal attribution, covariation detection, and prediction theory. (C) 2002 Elsevier Science (USA). All rights reserved.

The prediction of sharkskin instability observed during film blowing

The characteristics of sharkskin surface instability for linear low density polyethylene are studied as a function of film blowing processing conditions. By means of scanning electron microscopy and surface profilometry, is it found that for the standard industrial die geometry studied, sharkskin only occurs on the inside of the film bubble. Previous work suggests that this instability may be due to critical extensional stress levels at the exit of the die. Isothermal integral viscoelastic simulations of the annular extrusion process are reported, and confirm that the extensional stress at the die exit is large enough to cause local melt rupture. However the extensional stress level at the outer die wall predicts melt rupture of the outside bubble surface also, which contradicts the experimental findings. A significant temperature gradient is expected to exist across the die gap at the exit of the die, due to the external heating of the die and the low conductivity, of the polymer melt. It is shown that a gradient of 20 degreesC is required to cause sharkskin to only appear on the inner bubble surface.