Calculation of percolation thresholds in high dimensions for fcc, bcc, and diamond lattices

In a recent article, Galam and Mauger proposed an invariant for site and bond percolation thresholds, based on known values for twenty lattices (Eur. Phys. J. B 1 (1998) 255-258). Here we give a larger list of values for more than forty lattices in two to six dimensions. In this list are new results for fcc, bcc, and diamond lattices in 4, 5, and 6 dimensions. The list contains examples of lattices with equal site percolation thresholds, but different bond percolation thresholds. These and other examples show that there are deviations from the proposed invariant of up to 12% in two dimensions, increasing to 69% in higher dimensions.Comment: 12 pages, 3 figures (EPS), LaTe

HP-multigrid as smoother algorithm for higher order discontinuous Galerkin discretizations of advection dominated flows. Part I. Multilevel Analysis

The hp-Multigrid as Smoother algorithm (hp-MGS) for the solution of higher order accurate space-(time) discontinuous Galerkin discretizations of advection dominated flows is presented. This algorithm combines p-multigrid with h-multigrid at all p-levels, where the h-multigrid acts as smoother in the p-multigrid. The performance of the hp-MGS algorithm is further improved using semi-coarsening in combination with a new semi-implicit Runge-Kutta method as smoother. A detailed multilevel analysis of the hp-MGS algorithm is presented to obtain more insight into the theoretical performance of the algorithm. As model problem a fourth order accurate space-time discontinuous Galerkin discretization of the advection-diffusion equation is considered. The multilevel analysis shows that the hp-MGS algorithm has excellent convergence rates, both for low and high cell Reynolds numbers and on highly stretched meshes

HP-multigrid as smoother algorithm for higher order discontinuous Galerkin discretizations of advection dominated flows. Part II. Optimization of the Runge-Kutta smoother

Using a detailed multilevel analysis of the complete hp-Multigrid as Smoother algorithm accurate predictions are obtained of the spectral radius and operator norms of the multigrid error transformation operator. This multilevel analysis is used to optimize the coefficients in the semi-implicit Runge-Kutta smoother, such that the spectral radius of the multigrid error transformation operator is minimal under properly chosen constraints. The Runge-Kutta coefficients for a wide range of cell Reynolds numbers and a detailed analysis of the performance of the hp-MGS algorithm are presented. In addition, the computational complexity of the hp-MGS algorithm is investigated. The hp-MGS algorithm is tested on a fourth order accurate space-time discontinuous Galerkin finite element discretization of the advection-diffusion equation for a number of model problems, which include thin boundary layers and highly stretched meshes, and a non-constant advection velocity. For all test cases excellent multigrid convergence is obtained

Discrete Fourier analysis of multigrid algorithms

The main topic of this report is a detailed discussion of the discrete Fourier multilevel analysis of multigrid algorithms. First, a brief overview of multigrid methods is given for discretizations of both linear and nonlinear partial differential equations. Special attention is given to the hp-Multigrid as Smoother algorithm, which is a new algorithm suitable for higher order accurate discontinuous Galerkin discretizations of advection dominated flows. In order to analyze the performance of the multigrid algorithms the error transformation operator for several linear multigrid algorithms are derived. The operator norm and spectral radius of the multigrid error transformation are then computed using discrete Fourier analysis. First, the main operations in the discrete Fourier analysis are defined, including the aliasing of modes. Next, the Fourier symbol of the multigrid operators is computed and used to obtain the Fourier symbol of the multigrid error transformation operator. In the multilevel analysis, two and three level h-multigrid, both for uniformly and semi-coarsened meshes, are considered, and also the analysis of the hp-Multigrid as Smoother algorithm for three polynomial levels and three uniformly and semi-coarsened meshes. The report concludes with a discussion of the multigrid operator norm and spectral radius. In the appendix some useful auxiliary results are summarized

State-dependence of climate sensitivity: attractor constraints and palaeoclimate regimes

Equilibrium climate sensitivity (ECS) is a key predictor of climate change. However, it is not very well constrained, either by climate models or by observational data. The reasons for this include strong internal variability and forcing on many time scales. In practise this means that the 'equilibrium' will only be relative to fixing the slow feedback processes before comparing palaeoclimate sensitivity estimates with estimates from model simulations. In addition, information from the late Pleistocene ice age cycles indicates that the climate cycles between cold and warm regimes, and the climate sensitivity varies considerably between regime because of fast feedback processes changing relative strength and time scales over one cycle. In this paper we consider climate sensitivity for quite general climate dynamics. Using a conceptual Earth system model of Gildor and Tziperman (2001) (with Milankovich forcing and dynamical ocean biogeochemistry) we explore various ways of quantifying the state-dependence of climate sensitivity from unperturbed and perturbed model time series. Even without considering any perturbations, we suggest that climate sensitivity can be usefully thought of as a distribution that quantifies variability within the 'climate attractor' and where there is a strong dependence on climate state and more specificially on the 'climate regime' where fast processes are approximately in equilibrium. We also consider perturbations by instantaneous doubling of CO$_2$ and similarly find a strong dependence on the climate state using our approach.Comment: 32 pages, 10 figure

On the Kinetics of Body versus End Evaporation and Addition of Supramolecular Polymers

Although pathway-specific kinetic theories are fundamentally important to describe and understand reversible polymerisation kinetics, they come in principle at a cost of having a large number of system-specific parameters. Here, we construct a dynamical Landau theory to describe the kinetics of activated linear supramolecular self-assembly, which drastically reduces the number of parameters and still describes most of the interesting and generic behavior of the system in hand. This phenomenological approach hinges on the fact that if nucleated, the polymerisation transition resembles a phase transition. We are able to describe hysteresis, overshooting, undershooting and the existence of a lag time before polymerisation takes off, and pinpoint the conditions required for observing these types of phenomenon in the assembly and disassembly kinetics. We argue that the phenomenological kinetic parameter in our theory is a pathway controller, i.e., it controls the relative weights of the molecular pathways through which self-assembly takes place

SGR J1550–5418 Bursts Detected with the Fermi Gamma-Ray Burst Monitor during its Most Prolific Activity

We have performed detailed temporal and time-integrated spectral analysis of 286 bursts from SGR J1550–5418 detected with the Fermi Gamma-ray Burst Monitor (GBM) in 2009 January, resulting in the largest uniform sample of temporal and spectral properties of SGR J1550–5418 bursts. We have used the combination of broadband and high time-resolution data provided with GBM to perform statistical studies for the source properties. We determine the durations, emission times, duty cycles, and rise times for all bursts, and find that they are typical of SGR bursts. We explore various models in our spectral analysis, and conclude that the spectra of SGR J1550–5418 bursts in the 8-200 keV band are equally well described by optically thin thermal bremsstrahlung (OTTB), a power law (PL) with an exponential cutoff (Comptonized model), and two blackbody (BB) functions (BB+BB). In the spectral fits with the Comptonized model, we find a mean PL index of –0.92, close to the OTTB index of –1. We show that there is an anti-correlation between the Comptonized E_(peak) and the burst fluence and average flux. For the BB+BB fits, we find that the fluences and emission areas of the two BB functions are correlated. The low-temperature BB has an emission area comparable to the neutron star surface area, independent of the temperature, while the high-temperature BB has a much smaller area and shows an anti-correlation between emission area and temperature. We compare the properties of these bursts with bursts observed from other SGR sources during extreme activations, and discuss the implications of our results in the context of magnetar burst models

A micromechanical fracture analysis to investigate the effect of healing particles on the overall mechanical response of a self-healing particulate composite

A computational fracture analysis is conducted on a self‐healing particulate composite employing a finite element model of an actual microstructure. The key objective is to quantify the effects of the actual morphology and the fracture properties of the healing particles on the overall mechanical behaviour of the (MoSi2) particle‐dispersed Yttria Stabilised Zirconia (YSZ) composite. To simulate fracture, a cohesive zone approach is utilised whereby cohesive elements are embedded throughout the finite element mesh allowing for arbitrary crack initiation and propagation in the microstructure. The fracture behaviour in terms of the composite strength and the percentage of fractured particles is reported as a function of the mismatch in fracture properties between the healing particles and the matrix as well as a function of particle/matrix interface strength and fracture energy. The study can be used as a guiding tool for designing an extrinsic self‐healing material and understanding the effect of the healing particles on the overall mechanical properties of the material