A space-time adaptive method for simulating complex cardiac dynamics

For plane-wave and many-spiral states of the experimentally based Luo-Rudy 1 model of heart tissue in large (8 cm square) domains, we show that a space-time-adaptive time-integration algorithm can achieve a factor of 5 reduction in computational effort and memory-but without a reduction in accuracy-when compared to an algorithm using a uniform space-time mesh at the finest resolution. Our results indicate that such an algorithm can be extended straightforwardly to simulate quantitatively three-dimensional electrical dynamics over the whole human heart.     

A weight-function approach for determining watershed initial conditions for combustion waves

In many generic combustion models, one finds that a combustionwave will develop with a specific wave speed. However, thereare possible initial temperature profiles which do not evolveinto such waves, but rather die out to the ambient temperature.There can exist, in some models, a clear distinction betweenthose initial conditions that do evolve into combustion wavesand those that do not; this is sometimes referred to as thewatershed initial condition. When fuel consumption is consideredto be negligible, analytical methods can be used to obtain theexact watershed. In this paper, we consider the problem of determiningpseudo-watersheds and ascertaining the relationship betweenthese pseudo-watersheds and the exact watersheds. In the processa novel weight-function approach for infinite spatial domainsis developed.     

Efficient simulation of three-dimensional anisotropic cardiac tissue using an adaptive mesh refinement method

A recently developed space-time adaptive mesh refinement algorithm (AMRA) for simulating isotropic one- and two-dimensional excitable media is generalized to simulate three-dimensional anisotropic media. The accuracy and efficiency of the algorithm is investigated for anisotropic and inhomogeneous 2D and 3D domains using the Luo-Rudy 1 (LR1) and FitzHugh-Nagumo models. For a propagating wave in a 3D slab of tissue with LR1 membrane kinetics and rotational anisotropy comparable to that found in the human heart, factors of 50 and 30 are found, respectively, for the speedup and for the savings in memory compared to an algorithm using a uniform space-time mesh at the finest resolution of the AMRA method. For anisotropic 2D and 3D media, we find no reduction in accuracy compared to a uniform space-time mesh. These results suggest that the AMRA will be able to simulate the 3D electrical dynamics of canine ventricles quantitatively for 1 s using 32 1-GHz Alpha processors in approximately 9 h.     

Meson production in p + d reactions

Total and differential cross sections for the reactions p + d → ³He + m ⁰ with m=π, η and p + d → ³H+π⁺ were measured with the GEM detector at COSY for beam momenta between threshold and the maximum of the corresponding baryon resonance. For both reactions a strong forward-backward asymmetry was found. The data were compared with model calculations. The aspect of isospin symmetry breaking is studied. Representing the GEM Collaboration     

Power-law behavior of power spectra in low Prandtl number Rayleigh-Bénard convection

The origin of the power-law decay measured in the power spectra of low Prandtl number Rayleigh-Bénard convection near the onset of chaos is addressed using long time numerical simulations of the three-dimensional Boussinesq equations in cylindrical domains. The power law is found to arise from quasidiscontinuous changes in the slope of the time series of the heat transport associated with the nucleation of dislocation pairs and roll pinch-off events. For larger frequencies, the power spectra decay exponentially as expected for time continuous deterministic dynamics.