Numerical integrators for stiff and highly oscillatory differential equations /

Bibliography: p. 27-28.

Evaluation of an integration technique for the solution of nonlinear equations.

Thesis (M.S.)--University of Illinois at Urbana-Champaign.

Future developments in stiff integration techniques : stability of methods that do not use an exact Jacobian /

"COO-2383-0035."

Automatic integration of the heat equation /

Thesis (M. Sc.)--University of Illinois at Urbana-Champaign.

The automatic integration package for ordinary differential equations,

"COO-1469-0103."

Geometric integration theory

Mode of access: Internet.

Mechanism and model of the oscillatory electro-oxidation of methanol

A mechanism for the kinetic instabilities observed in the galvanostatic electro-oxidation of methanol is suggested and a model developed. The model is investigated using stoichiometric network analysis as well as concepts from algebraic geometry (polynomial rings and ideal theory) revealing the occurrence of a Hopf and a saddle-node bifurcation. These analytical solutions are confirmed by numerical integration of the system of differential equations. (C) 2010 American Institute of Physics

Gravitational quasinormal modes of AdS black branes in d spacetime dimensions

The AdS/CFT duality has established a mapping between quantities in the bulk AdS black-hole physics and observables in a boundary finite-temperature field theory. Such a relationship appears to be valid for an arbitrary number of spacetime dimensions, extrapolating the original formulations of Maldacena`s correspondence. In the same sense properties like the hydrodynamic behavior of AdS black-hole fluctuations have been proved to be universal. We investigate in this work the complete quasinormal spectra of gravitational perturbations of d-dimensional plane-symmetric AdS black holes (black branes). Holographically the frequencies of the quasinormal modes correspond to the poles of two-point correlation functions of the field-theory stress-energy tensor. The important issue of the correct boundary condition to be imposed on the gauge-invariant perturbation fields at the AdS boundary is studied and elucidated in a fully d-dimensional context. We obtain the dispersion relations of the first few modes in the low-, intermediate- and high-wavenumber regimes. The sound-wave (shear-mode) behavior of scalar (vector)-type low- frequency quasinormal mode is analytically and numerically confirmed. These results are found employing both a power series method and a direct numerical integration scheme.

Buckling and post-buckling of extensible rods revisited: A multiple-scale solution

An exact non-linear formulation of the equilibrium of elastic prismatic rods subjected to compression and planar bending is presented, electing as primary displacement variable the cross-section rotations and taking into account the axis extensibility. Such a formulation proves to be sufficiently general to encompass any boundary condition. The evaluation of critical loads for the five classical Euler buckling cases is pursued, allowing for the assessment of the axis extensibility effect. From the quantitative viewpoint, it is seen that such an influence is negligible for very slender bars, but it dramatically increases as the slenderness ratio decreases. From the qualitative viewpoint, its effect is that there are not infinite critical loads, as foreseen by the classical inextensible theory. The method of multiple (spatial) scales is used to survey the post-buckling regime for the five classical Euler buckling cases, with remarkable success, since very small deviations were observed with respect to results obtained via numerical integration of the exact equation of equilibrium, even when loads much higher than the critical ones were considered. Although known beforehand that such classical Euler buckling cases are imperfection insensitive, the effect of load offsets were also looked at, thus showing that the formulation is sufficiently general to accommodate this sort of analysis. (c) 2008 Elsevier Ltd. All rights reserved.

Likelihood-based estimation of the regression model with scrambled responses

A significant problem in the collection of responses to potentially sensitive questions, such as relating to illegal, immoral or embarrassing activities, is non-sampling error due to refusal to respond or false responses. Eichhorn & Hayre (1983) suggested the use of scrambled responses to reduce this form of bias. This paper considers a linear regression model in which the dependent variable is unobserved but for which the sum or product with a scrambling random variable of known distribution, is known. The performance of two likelihood-based estimators is investigated, namely of a Bayesian estimator achieved through a Markov chain Monte Carlo (MCMC) sampling scheme, and a classical maximum-likelihood estimator. These two estimators and an estimator suggested by Singh, Joarder & King (1996) are compared. Monte Carlo results show that the Bayesian estimator outperforms the classical estimators in almost all cases, and the relative performance of the Bayesian estimator improves as the responses become more scrambled.

Periodic electric field enhanced transport through membranes

We examine the mean flux across a homogeneous membrane of a charged tracer subject to an alternating, symmetric voltage waveform. The analysis is based on the Nernst-Planck flux equation, with electric field subject to time dependence only. For low frequency electric fields the quasi steady-state flux can be approximated using the Goldman model, which has exact analytical solutions for tracer concentration and flux. No such closed form solutions can be found for arbitrary frequencies, however we find approximations for high frequency. An approximation formula for the average flux at all frequencies is also obtained from the two limiting approximations. Numerical integration of the governing equation is accomplished by use of the numerical method of lines and is performed for four different voltage waveforms. For the different voltage profiles, comparisons are made with the approximate analytical solutions which demonstrates their applicability. (c) 2005 Elsevier B.V. All rights reserved.

Relative dispersion of intravascular transit times during isolated human limb perfusions for recurrent melanoma

Aims We have characterized the relative dispersion of vascular and extravascular markers in the limbs of three patients undergoing isolated limb perfusions with the cytotoxic melphalan for recurrent malignant melanoma both before and after melphalan dosing. Methods A bolus of injectate containing [Cr-51] labelled red blood cells, [C-14]-sucrose and [H-3]-water was injected into an iliac or femoral artery and outflow samples collected at 1 s intervals by a fraction collector. The radioactivity due to each isotype was analysed by either gamma [Cr-51] or beta [C-14 and H-3] counting. The moments of the outflow fraction-time profiles were estimated by a nonparametric (numerical integration) method and a parametric model (sum of two inverse Gaussian functions). Results The availability, mean transit time and normalised variance (CV2) obtained for labelled red blood cells, sucrose and water were similar before and after melphalan dosing and with the two methods of calculation but varied between the patients. Conclusions The vascular space is not well-stirred but characterized by a CV2 similar that reported previously for in situ rat hind limb and rat liver perfusions. A flow-limited blood-tissue exchange was observed for the permeating indicators. Administration of melphalan did not influence the distribution characteristics of the indicators.

Model-Predicted Performance of Second-Trimester Down Syndrome Screening With Sonographic Prenasal Thickness

Objective. The purpose of this study was to estimate the Down syndrome detection and false-positive rates for second-trimester sonographic prenasal thickness (PT) measurement alone and in combination with other markers. Methods. Multivariate log Gaussian modeling was performed using numerical integration. Parameters for the PT distribution, in multiples of the normal gestation-specific median (MoM), were derived from 105 Down syndrome and 1385 unaffected pregnancies scanned at 14 to 27 weeks. The data included a new series of 25 cases and 535 controls combined with 4 previously published series. The means were estimated by the median and the SDs by the 10th to 90th range divided by 2.563. Parameters for other markers were obtained from the literature. Results. A log Gaussian model fitted the distribution of PT values well in Down syndrome and unaffected pregnancies. The distribution parameters were as follows: Down syndrome, mean, 1.334 MoM; log(10) SD, 0.0772; unaffected pregnancies, 0.995 and 0.0752, respectively. The model-predicted detection rates for 1%, 3%, and 5% false-positive rates for PT alone were 35%, 51%, and 60%, respectively. The addition of PT to a 4 serum marker protocol increased detection by 14% to 18% compared with serum alone. The simultaneous sonographic measurement of PT and nasal bone length increased detection by 19% to 26%, and with a third sonographic marker, nuchal skin fold, performance was comparable with first-trimester protocols. Conclusions. Second-trimester screening with sonographic PT and serum markers is predicted to have a high detection rate, and further sonographic markers could perform comparably with first-trimester screening protocols.

Studies of solute self-association by sedimentation equilibrium: allowance for effects of thermodynamic non-ideality beyond the consequences of nearest-neighbor interactions

A sedimentation equilibrium study of a-chymotrypsin self-association in acetate-chloride buffer, pH 4.1 I 0.05, has been used to illustrate determination of a dimerization constant under conditions where thermodynamic non-ideality is manifested beyond the consequences of nearest-neighbor interactions. Because the expressions for the experimentally determinable interaction parameters comprise a mixture of equilibrium constant and excluded volume terms, the assignment of reasonable magnitudes to the relevant virial coefficients describing non-associative cluster formation is essential for the evaluation of a reliable estimate of the dimerization constant. Determination of these excluded volume parameters by numerical integration over the potential-of-mean-force is shown to be preferable to their calculation by approximate analytical solutions of the integral for this relatively small enzyme monomer with high net charge (+ 10) under conditions of low ionic strength (0.05 M). (C) 2001 Elsevier Science B.V. All rights reserved.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.