Stiffly accurate Runge-Kutta methods for stiff stochastic differential equations

In this paper we discuss implicit methods based on stiffly accurate Runge-Kutta methods and splitting techniques for solving Stratonovich stochastic differential equations (SDEs). Two splitting techniques: the balanced splitting technique and the deterministic splitting technique, are used in this paper. We construct a two-stage implicit Runge-Kutta method with strong order 1.0 which is corrected twice and no update is needed. The stability properties and numerical results show that this approach is suitable for solving stiff SDEs. (C) 2001 Elsevier Science B.V. All rights reserved.

A new wavelet-based method for the solution of the population balance equation

A new wavelet-based method for solving population balance equations with simultaneous nucleation, growth and agglomeration is proposed, which uses wavelets to express the functions. The technique is very general, powerful and overcomes the crucial problems of numerical diffusion and stability that often characterize previous techniques in this area. It is also applicable to an arbitrary grid to control resolution and computational efficiency. The proposed technique has been tested for pure agglomeration, simultaneous nucleation and growth, and simultaneous growth and agglomeration. In all cases, the predicted and analytical particle size distributions are in excellent agreement. The presence of moving sharp fronts can be addressed without the prior investigation of the characteristics of the processes. (C) 2001 Published by Elsevier Science Ltd.

Simultaneous administration of a cocktail of markers to measure renal drug elimination pathways: absence of a pharmacokinetic interaction between fluconazole and sinistrin, p-aminohippuric acid and pindolol

Aims Previous studies suggest that estimated creatinine clearance, the conventional measure of renal function, does not adequately reflect charges in renal drug handling in some patients, including the immunosuppressed. The aim of this study was to develop and validate a cocktail of markers. to be given in a single administration, capable of detecting alterations in the renal elimination pathways of glomerular filtration, tubular secretion and tubular reabsorption. Methods Healthy male subjects (n = 12) received intravenously infused 2500 mg sinistrin (glomerular filtration) and 440 mg p-aminohippuric acid (PAH; anion secretion), and orally administered 100 mg fluconazole (reabsorption) and 15 mg rac-pindolol (cation secretion). The potential interaction between these markers was investigated in a pharmacokinetic study where markers (M) or fluconazole (F) were administered alone or together (M + F). Validated analytical methods were used to measure plasma and urine concentrations in order to quantify the renal handling of each marker. Plasma protein binding of fluconazole was measured by ultrafiltration. All subjects had an estimated creatinine clearance within the normal range. The renal clearance of each marker (Mean +/- s.d.) was calculated as the ratio of the amount excreted in urine and thearea-under-the-concentration-time curve. Statistical comparisons were made using a paired t-test and 95% confidence intervals were reported. Results The renal clearances of sinistrin (M: 119 +/- 31 ml min(-1); M + F: 130 +/- 40 ml min(-1); P = 0.32), PAH (M: 469 +/- 145 ml min(-1); M + F: 467 +/- 146 ml min(-1); P = 0.95), R-pindolol (M: 204 +/- 41 ml min(-1); M + F: 190 +/- 41 ml min(-1); P = 0.39; n = 11), S-pindolol (M: 225 +/- 55 ml min(-1); M + F: 209 +/- 60 ml min(-1); P = 0.27; n = 11) and fluconazole (F: 14.9 +/-3.8 ml min(-1); M + F: 13.6 +/- 3.4 ml min(-1); P = 0.16) were similar when the markers or fluconazole were administered alone (M or F) or as a cocktail (M + F). Conclusions This study found no interaction between markers and fluconazole in healthy male subjects, suggesting that a single administration of this cocktail of markers of different renal processes call be used to simultaneously investigate pathways of renal drug elimination.

Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

Theory and applications of fuzzy Volterra integral equations

Formulations of fuzzy integral equations in terms of the Aumann integral do not reflect the behavior of corresponding crisp models. Consequently, they are ill-adapted to describe physical phenomena, even when vagueness and uncertainty are present. A similar situation for fuzzy ODEs has been obviated by interpretation in terms of families of differential inclusions. The paper extends this formalism to fuzzy integral equations and shows that the resulting solution sets and attainability sets are fuzzy and far better descriptions of uncertain models involving integral equations. The investigation is restricted to Volterra type equations with mildly restrictive conditions, but the methods are capable of extensive generalization to other types and more general assumptions. The results are illustrated by integral equations relating to control models with fuzzy uncertainties.

Controlling the complex Lorenz equations by modulation

We demonstrate that a system obeying the complex Lorenz equations in the deep chaotic regime can be controlled to periodic behavior by applying a modulation to the pump parameter. For arbitrary modulation frequency and amplitude there is no obvious simplification of the dynamics. However, we find that there are numerous windows where the chaotic system has been controlled to different periodic behaviors. The widths of these windows in parameter space are narrow, and the positions are related to the ratio of the modulation frequency of the pump to the average pulsation frequency of the output variable. These results are in good agreement with observations previously made in a far-infrared laser system.

On positive solutions of nonlinear elliptic equations involving concave and critical nonlinearities

We study the existence of nonnegative solutions of elliptic equations involving concave and critical Sobolev nonlinearities. Applying various variational principles we obtain the existence of at least two nonnegative solutions.

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.

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.

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.