FISH: A 3D parallel MHD code for astrophysical applications

FISH is a fast and simple ideal magneto-hydrodynamics code that scales to ~10 000 processes for a Cartesian computational domain of ~1000^3 cells. The simplicity of FISH has been achieved by the rigorous application of the operator splitting technique, while second order accuracy is maintained by the symmetric ordering of the operators. Between directional sweeps, the three-dimensional data is rotated in memory so that the sweep is always performed in a cache-efficient way along the direction of contiguous memory. Hence, the code only requires a one-dimensional description of the conservation equations to be solved. This approach also enable an elegant novel parallelisation of the code that is based on persistent communications with MPI for cubic domain decomposition on machines with distributed memory. This scheme is then combined with an additional OpenMP parallelisation of different sweeps that can take advantage of clusters of shared memory. We document the detailed implementation of a second order TVD advection scheme based on flux reconstruction. The magnetic fields are evolved by a constrained transport scheme. We show that the subtraction of a simple estimate of the hydrostatic gradient from the total gradients can significantly reduce the dissipation of the advection scheme in simulations of gravitationally bound hydrostatic objects. Through its simplicity and efficiency, FISH is as well-suited for hydrodynamics classes as for large-scale astrophysical simulations on high-performance computer clusters. In preparation for the release of a public version, we demonstrate the performance of FISH in a suite of astrophysically orientated test cases.Comment: 27 pages, 11 figure

Breakdown of Kolmogorov scaling in models of cluster aggregation with deposition

The steady state of the model of cluster aggregation with deposition is characterized by a constant flux of mass directed from small masses towards large masses. It can therefore be studied using phenomenological theories of turbulence, such as Kolmogorov's 1941 theory. On the other hand, the large scale behavior of the aggregation model in dimensions lower than or equal to two is governed by a perturbative fixed point of the renormalization group flow, which enables an analytic study of the scaling properties of correlation functions in the steady state. In this paper, we show that the correlation functions have multifractal scaling, which violates linear Kolmogorov scaling. The analytical results are verified by Monte Carlo simulations.Comment: 5 pages 4 figure

Nonequilibrium phase transitions in models of adsorption and desorption

The nonequilibrium phase transition in a system of diffusing, coagulating particles in the presence of a steady input and evaporation of particles is studied. The system undergoes a transition from a phase in which the average number of particles is finite to one in which it grows linearly in time. The exponents characterizing the mass distribution near the critical point are calculated in all dimensions.Comment: 10 pages, 2 figures (To appear in Phys. Rev. E

Prolonged decrease in heart rate variability after elective hip arthroplasty

The pattern of postoperative heart rate variability may provide insight into the response of the autonomic nervous system to anaesthesia and surgery. We have obtained spectral (fast Fourier transform) and non-spectral indices of heart rate variability from electrocardiographic recordings, sampled during continuous perioperative Holter monitoring in 15 otherwise healthy patients with an uncomplicated postoperative course, undergoing elective hip arthroplasty with either spinal or general anaesthesia. In both groups, total spectral energy (0.01-1 Hz), low-frequency spectral energy (0.01-0.15 Hz) and high-frequency spectral energy (0.15-0.40 Hz) decreased after surgery to 32% (95% confidence interval (Cl) 10.5; P <0.01), 29% (95% Cl 12.5; P <0.07; and 33% (95% Cl 12.5; P <0.01) of their preoperative values, respectively, and these indices remained suppressed for up to 5 days. Non-spectral indices decreased to a similar extent. These findings indicate a substantial and prolonged postoperative decrease in both parasympathetic and sympathetic influence on the sinus nod

Gravitational waves from supernova matter

We have performed a set of 11 three-dimensional magnetohydrodynamical core collapse supernova simulations in order to investigate the dependencies of the gravitational wave signal on the progenitor's initial conditions. We study the effects of the initial central angular velocity and different variants of neutrino transport. Our models are started up from a 15 solar mass progenitor and incorporate an effective general relativistic gravitational potential and a finite temperature nuclear equation of state. Furthermore, the electron flavour neutrino transport is tracked by efficient algorithms for the radiative transfer of massless fermions. We find that non- and slowly rotating models show gravitational wave emission due to prompt- and lepton driven convection that reveals details about the hydrodynamical state of the fluid inside the protoneutron stars. Furthermore we show that protoneutron stars can become dynamically unstable to rotational instabilities at T/|W| values as low as ~2 % at core bounce. We point out that the inclusion of deleptonization during the postbounce phase is very important for the quantitative GW prediction, as it enhances the absolute values of the gravitational wave trains up to a factor of ten with respect to a lepton-conserving treatment.Comment: 10 pages, 6 figures, accepted, to be published in a Classical and Quantum Gravity special issue for MICRA200

Macroscopic Equations of Motion for Two Phase Flow in Porous Media

The established macroscopic equations of motion for two phase immiscible displacement in porous media are known to be physically incomplete because they do not contain the surface tension and surface areas governing capillary phenomena. Therefore a more general system of macroscopic equations is derived here which incorporates the spatiotemporal variation of interfacial energies. These equations are based on the theory of mixtures in macroscopic continuum mechanics. They include wetting phenomena through surface tensions instead of the traditional use of capillary pressure functions. Relative permeabilities can be identified in this approach which exhibit a complex dependence on the state variables. A capillary pressure function can be identified in equilibrium which shows the qualitative saturation dependence known from experiment. In addition the new equations allow to describe the spatiotemporal changes of residual saturations during immiscible displacement.Comment: 15 pages, Phys. Rev. E (1998), in prin

Phase Transition in the Takayasu Model with Desorption

We study a lattice model where particles carrying different masses diffuse, coalesce upon contact, and also unit masses adsorb to a site with rate $q$ or desorb from a site with nonzero mass with rate $p$. In the limit $p=0$ (without desorption), our model reduces to the well studied Takayasu model where the steady-state single site mass distribution has a power law tail $P(m)\sim m^{-\tau}$ for large mass. We show that varying the desorption rate $p$ induces a nonequilibrium phase transition in all dimensions. For fixed $q$, there is a critical $p_c(q)$ such that if $p<p_c(q)$, the steady state mass distribution, $P(m)\sim m^{-\tau}$ for large $m$ as in the Takayasu case. For $p=p_c(q)$, we find $P(m)\sim m^{-\tau_c}$ where $\tau_c$ is a new exponent, while for $p>p_c(q)$, $P(m)\sim \exp(-m/m^*)$ for large $m$. The model is studied analytically within a mean field theory and numerically in one dimension.Comment: RevTex, 11 pages including 5 figures, submitted to Phys. Rev.

Using state variables to model the response of tumour cells to radiation and heat: a novel multi-hit-repair approach

In order to overcome the limitations of the linear-quadratic model and include synergistic effects of heat and radiation, a novel radiobiological model is proposed. The model is based on a chain of cell populations which are characterized by the number of radiation induced damages (hits). Cells can shift downward along the chain by collecting hits and upward by a repair process. The repair process is governed by a repair probability which depends upon state variables used for a simplistic description of the impact of heat and radiation upon repair proteins. Based on the parameters used, populations up to 4-5 hits are relevant for the calculation of the survival. The model describes intuitively the mathematical behaviour of apoptotic and nonapoptotic cell death. Linear-quadratic-linear behaviour of the logarithmic cell survival, fractionation, and (with one exception) the dose rate dependencies are described correctly. The model covers the time gap dependence of the synergistic cell killing due to combined application of heat and radiation, but further validation of the proposed approach based on experimental data is needed. However, the model offers a work bench for testing different biological concepts of damage induction, repair, and statistical approaches for calculating the variables of state

Ionic liquids at electrified interfaces

Until recently, “room-temperature” (<100–150 °C) liquid-state electrochemistry was mostly electrochemistry of diluted electrolytes(1)–(4) where dissolved salt ions were surrounded by a considerable amount of solvent molecules. Highly concentrated liquid electrolytes were mostly considered in the narrow (albeit important) niche of high-temperature electrochemistry of molten inorganic salts(5-9) and in the even narrower niche of “first-generation” room temperature ionic liquids, RTILs (such as chloro-aluminates and alkylammonium nitrates).(10-14) The situation has changed dramatically in the 2000s after the discovery of new moisture- and temperature-stable RTILs.(15, 16) These days, the “later generation” RTILs attracted wide attention within the electrochemical community.(17-31) Indeed, RTILs, as a class of compounds, possess a unique combination of properties (high charge density, electrochemical stability, low/negligible volatility, tunable polarity, etc.) that make them very attractive substances from fundamental and application points of view.(32-38) Most importantly, they can mix with each other in “cocktails” of one’s choice to acquire the desired properties (e.g., wider temperature range of the liquid phase(39, 40)) and can serve as almost “universal” solvents.(37, 41, 42) It is worth noting here one of the advantages of RTILs as compared to their high-temperature molten salt (HTMS)(43) “sister-systems”.(44) In RTILs the dissolved molecules are not imbedded in a harsh high temperature environment which could be destructive for many classes of fragile (organic) molecules