Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

New integrable boundary conditions for the q-deformed supersymmetric U model and Bethe ansatz equations

New classes of integrable boundary conditions for the q-deformed (or two-parameter) supersymmetric U model are presented. The boundary systems are solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations are derived. (C) 1998 Elsevier Science B.V.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Critical values for unit root and cointegration test statistics - the use of response surface equations

This note considers the value of surface response equations which can be used to calculate critical values for a range of unit root and cointegration tests popular in applied economic research.

An evaluation of the time dependence of the anisotropy of the water diffusion tensor in acute human ischemia

We have performed MRI examinations to determine the water diffusion tensor in the brain of six patients who were admitted to the hospital within 12 h after the onset of cerebral ischemic symptoms. The examinations have been carried out immediately after admission, and thereafter at varying intervals up to 90 days post admission. Maps of the trace of the diffusion tensor, the fractional anisotropy and the lattice index, as well as maps of cerebral blood perfusion parameters, were generated to quantitatively assess the character of the water diffusion tensor in the infarcted area. In patients with significant perfusion deficits and substantial lesion volume changes, four of six cases, our measurements show a monotonic and significant decrease in the diffusion anisotropy within the ischemic lesion as a function of time. We propose that retrospective analysis of this quantity, in combination with brain tissue segmentation and cerebral perfusion maps, may be used in future studies to assess the severity of the ischemic event. (C) 1999 Elsevier Science Inc.

Water diffusion in hydroxyethyl methacrylate (HEMA)-based hydrogels formed by gamma-radiolysis

Polymer hydrogels based upon methacrylates are used extensively in the pharmaceutical industry, particularly as controlled release drug delivery systems. These materials are generally prepared by chemically initiated polymerization, but this can lead to the presence of unwanted initiator fragments in the polymer matrix. In the present work, initiation of polymerization by gamma-irradiation of hydroxyethyl methacrylate, with and without added crosslinkers, has been investigated, and the diffusion coefficients for water in the resulting polymers have been measured through mass uptake by the polymers. The diffusion of water in poly(hydroxyethyl methacrylate) at 310 K was found to be Fickian, with a diffusion coefficient of 1.96 +/- 0.1 x 10(11) m(2) s(-1) and an equilibrium water content of 58%, NMR imaging analyses confirmed the adherance to a Fickian model of the diffusion of water into polymer cylinders. The incorporation of small amounts (0.2-0.5 wt%) of added ethyleneglycol-dimethacrylate-based crosslinkers was found to have only a small effect on the diffusion coefficient and the equilibrium water content for the copolymers. (C) 1999 Society of Chemical Industry.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Loss of connectivity in Alzheimer's disease: an evaluation of white matter tract integrity with colour coded MR diffusion tensor imaging

A novel MRI method-diffusion tensor imaging-was used to compare the integrity of several white matter fibre tracts in patients with probable Alzheimer's disease. Relative to normal controls, patients with probable Alzheimer's disease showed a highly significant reduction in the integrity of the association white matter fibre tracts, such as the splenium of the corpus callosum, superior longitudinal fasciculus, and cingulum. By contrast, pyramidal tract integrity seemed unchanged. This novel finding is consistent with the clinical presentation of probable Alzheimer's disease, in which global cognitive decline is a more prominent feature than motor disturbance.

Mathematical models in percutaneous absorption

A number of mathematical models have been used to describe percutaneous absorption kinetics. In general, most of these models have used either diffusion-based or compartmental equations. The object of any mathematical model is to a) be able to represent the processes associated with absorption accurately, b) be able to describe/summarize experimental data with parametric equations or moments, and c) predict kinetics under varying conditions. However, in describing the processes involved, some developed models often suffer from being of too complex a form to be practically useful. In this chapter, we attempt to approach the issue of mathematical modeling in percutaneous absorption from four perspectives. These are to a) describe simple practical models, b) provide an overview of the more complex models, c) summarize some of the more important/useful models used to date, and d) examine sonic practical applications of the models. The range of processes involved in percutaneous absorption and considered in developing the mathematical models in this chapter is shown in Fig. 1. We initially address in vitro skin diffusion models and consider a) constant donor concentration and receptor conditions, b) the corresponding flux, donor, skin, and receptor amount-time profiles for solutions, and c) amount- and flux-time profiles when the donor phase is removed. More complex issues, such as finite-volume donor phase, finite-volume receptor phase, the presence of an efflux. rate constant at the membrane-receptor interphase, and two-layer diffusion, are then considered. We then look at specific models and issues concerned with a) release from topical products, b) use of compartmental models as alternatives to diffusion models, c) concentration-dependent absorption, d) modeling of skin metabolism, e) role of solute-skin-vehicle interactions, f) effects of vehicle loss, a) shunt transport, and h) in vivo diffusion, compartmental, physiological, and deconvolution models. We conclude by examining topics such as a) deep tissue penetration, b) pharmacodynamics, c) iontophoresis, d) sonophoresis, and e) pitfalls in modeling.

Diffusion- and perfusion-weighted MRI response to thrombolysis in stroke

Diffusion- and perfusion-weighted magnetic resonance imaging provides important pathophysiological information in acute bra-in ischemia. We performed a prospective study in 19 sub-6-hour stroke patients using serial diffusion- and perfusion-weighted imaging before intravenous thrombolysis, with repeat studies, both subacutely and at outcome. For comparison of ischemic lesion evolution and clinical outcome, we used a historical control group of 21 sub-6-hour ischemic stroke patients studied serially with diffusion- and perfusion-weighted imaging. The two groups were well matched for the baseline National Institutes of Health Stroke Scale and magnetic resonance parameters. Perfusion-weighted imaging-diffusion-weighted imaging mismatch was present in 16 of 19 patients treated with tissue plasminogen activator, and 16 of 21 controls. Perfusion-weighted imaging-diffusion-weighted imaging mismatch patients treated with tissue plaminogen activator had higher recanalization rates and enhanced reperfusion at day 3 (81% vs 47% in controls), and a greater proportion of severely hypoperfused acute mismatch tissue not progressing to infarction (82% vs -25% in controls). Despite similar baseline diffusion-weighted imaging lesions, infarct expansion was less in the recombinant tissue plaminogen activator group (14cm(3) vs 56cm(3) in controls). The positive effect of thrombolysis on lesion growth in mismatch patients translated into a greater improvement in baseline to outcome National Institutes of Health Stroke Scale in the group treated with recombinant tissue plaminogen activator, and a significantly larger proportion of patients treated with recombinant tissue plaminogen activator having a clinically meaningful improvement in National Institutes of Health Stroke Scale of;2:7 points. The natural evolution of acute perfusion-weighted imaging-diffusion-weighted imaging mismatch tissue may be altered by thrombolysis, with improved stroke outcome. This has implications for the use of diffusion- and perfusion-weighted imaging in selecting and monitoring patients for thrombolytic therapy.

Quantum error correction for continuously detected errors

We show that quantum feedback control can be used as a quantum-error-correction process for errors induced by a weak continuous measurement. In particular, when the error model is restricted to one, perfectly measured, error channel per physical qubit, quantum feedback can act to perfectly protect a stabilizer codespace. Using the stabilizer formalism we derive an explicit scheme, involving feedback and an additional constant Hamiltonian, to protect an (n-1)-qubit logical state encoded in n physical qubits. This works for both Poisson (jump) and white-noise (diffusion) measurement processes. Universal quantum computation is also possible in this scheme. As an example, we show that detected-spontaneous emission error correction with a driving Hamiltonian can greatly reduce the amount of redundancy required to protect a state from that which has been previously postulated [e.g., Alber , Phys. Rev. Lett. 86, 4402 (2001)].

Asymptotically almost periodic solutions for abstract neutral integro-differential equations

We study the existence of asymptotically almost periodic classical solutions for a class of abstract neutral integro-differential equation with unbounded delay. A concrete application to partial neutral integro-differential equations which arise in the study of heat conduction in fading memory material is considered. (C) 2011 Elsevier Inc. All rights reserved.